Immobilized metal ion affinity chromatography (IMAC) is an important chromatographic technique for biomolecules. In order to get a detailed understanding of the hydration of immobilized Fe(III), complexes of Fe(III) with methyl substituted iminodiacetate ([Fe(MSIDA)(H2O)(3)](+)) as well as with methyl substituted nitrilotriacetate ([Fe(MSNTA)(H2O)(2)]) were simulated in aqueous solutions with the quantum mechanical charge field molecular dynamics (QMCF MD) approach. The simulations were carried out at the Hartree-Fock (HF) level of theory, since cluster calculations at the HF, MP2, and B3LYP levels of theory showed that this method results in a good compromise between computational effort and accuracy. None of the coordinating water molecules were exchanged during the simulation period of 15 ps. The Fe-OH2O bond distances as well as the Fe-OH2O stretching motions differed among the coordinating water molecules, indicating different bond strengths. For the water molecules in the second hydration layer, mean residence times of 2.7 and 1.9 ps were obtained for [Fe(MSIDA)(H2O)(3)](+) and [Fe(MSNTA)(H2O)(2)], respectively. Furthermore, infrared measurements were carried out to characterize the most prominent bond features of aqueous Fe(III)-NTA and to discuss these results in conjunction with the computationally derived frequencies.