The magnitude of equilibrium iron isotope fractionation between Fe(H 2O) 6 3+ and Fe(H 2O) 6 2+ is calculated using density functional theory (DFT) and compared to prior theoretical and experimental results. DFT is a quantum chemical approach that permits a priori estimation of all vibrational modes and frequencies of these complexes and the effects of isotopic substitution. This information is used to calculate reduced partition function ratios of the complexes (10 3 · ln(β)), and hence, the equilibrium isotope fractionation factor (10 3 · ln(α)). Solvent effects are considered using the polarization continuum model (PCM). DFT calculations predict fractionations of several per mil in 56Fe/ 54Fe favoring partitioning of heavy isotopes in the ferric complex. Quantitatively, 10 3 · ln(α) predicted at 22°C, ∼ 3 ‰, agrees with experimental determinations but is roughly half the size predicted by prior theoretical results using the Modified Urey-Bradley Force Field (MUBFF) model. Similar comparisons are seen at other temperatures. MUBFF makes a number of simplifying assumptions about molecular geometry and requires as input IR spectroscopic data. The difference between DFT and MUBFF results is primarily due to the difference between the DFT-predicted frequency for the ν 4 mode (O-Fe-O deformation) of Fe(H 2O) 6 3+ and spectroscopic determinations of this frequency used as input for MUBFF models (185–190 cm −1 vs. 304 cm −1, respectively). Hence, DFT-PCM estimates of 10 3 · ln(β) for this complex are ∼ 20% smaller than MUBFF estimates. The DFT derived values can be used to refine predictions of equilibrium fractionation between ferric minerals and dissolved ferric iron, important for the interpretation of Fe isotope variations in ancient sediments. Our findings increase confidence in experimental determinations of the Fe(H 2O) 6 3+ − Fe(H 2O) 6 2+ fractionation factor and demonstrate the utility of DFT for applications in “heavy” stable isotope geochemistry.
Read full abstract