  A discrepancy has been reported in the literature between experimental
measurements and density functional theory (DFT) calculations, on the one hand,
and simulations based on empirical potential energy functions on the other
hand, regarding the solvation shell of the aqueous $\mathrm{Fe}^{3+}$ ion; the
former studies reporting a coordination number of 6.0 corresponding to a well
defined octahedral solvation shell, while the latter studies report a higher
value, typically around 6.5. The simulations presented here were conducted for
long time intervals and a wide range in ion concentration, using various
potential energy functions as well as DFT calculations for several of the
generated configurations. The results indicate that the solvation shell can
undergo abrupt transitions between two well-separated states: an octahedral
state (OH) with 6-fold coordination and a high degree of sphericity, $\phi^6$,
and a capped trigonal prism (CTP) state with 7-fold coordination and
significantly smaller $\phi^6$. In dilute $\mathrm{FeCl_3}$ solutions, with ion
concentration on the order of 0.1 $\mathrm{mol/dm^3}$
($[\mathrm{Fe}^{3+}]/[\mathrm{H_2O}]$ ratio of 1:600), the lifetime of the two
states is similar, on the order of a ns, but the lifetime of the OH state
increases with ion concentration, while that of the CTP state remains roughly
constant. When a uniform negative background is used instead of explicit
counterions, as is often done in DFT calculations, the lifetime of the OH state
increases strongly, quadrupling as the ion concentration is increased to 0.9
$\mathrm{mol/dm^3}$. Experimental measurements of lower concentration
solutions, as well as corresponding simulations based on electronic structure
calculations spanning at least several nanoseconds, are needed to test the
existence of the predicted CTP state.
