Gas-phase density functional theory calculations (B3LYP, double zeta plus polarization basis sets) are used to predict the solution-phase infrared spectra for a series of CO- and CN-containing iron complexes. It is shown that simple linear scaling of the computed C--O and C--N stretching frequencies yields accurate predictions of the the experimentally determined nu(CO) and nu(CN) values for a variety of complexes of different charges and in solvents of varying polarity. As examples of the technique, the resulting correlation is used to assign structures to spectroscopically observed but structurally ambiguous species in two different systems. For the (mu-SCH2CH2CH2S)[Fe(CO)3]2 complex in tetrahydrofuran solution, our calculations show that the initial electrochemical reduction process leads to a simple one-electron reduced product with a structure very similar to the (mu-SCH2CH2CH2S)[Fe(CO)3]2 parent complex. For the iron-iron hydrogenase enzyme active site, our computations show that the absence or presence of a water molecule near the distal iron center (the iron center further from the [4Fe4S] cluster and protein backbone) has very little effect on the predicted infrared spectra.