Molecular dynamics simulations of single-ion Fe2+(aq) and Fe3+(aq) solutions have been performed with two rigid-water models (SPC and SPC/E) and a newly constructed SPC-based flexible-water model (SPC+CCL). The SPC+CCL water model in combination with effective Fe2+ and Fe3+ ion-water potentials manages to reproduce many experimental structural and dynamical properties of the solutions. Special attention is given to the large ion-induced frequency shifts of the OH stretching bands, which are also well reproduced by the SPC+CCL model.