Calculation of site-specific carbon-isotope fractionation in pedogenic oxide minerals

Abstract

Ab initio molecular dynamics and quantum chemistry techniques are used to calculate the structure, vibrational frequencies, and carbon-isotope fractionation factors of the carbon dioxide component [CO2(m)] of soil (oxy)hydroxide minerals goethite, diaspore, and gibbsite. We have identified two possible pathways of incorporation of CO2(m) into (oxy)hydroxide crystal structures: one in which the C4+ substitutes for four H+ [CO2(m)A] and another in which C4+ substitutes for (Al3+,Fe3+) + H+ [CO2(m)B]. Calculations of isotope fractionation factors give large differences between the two structures, with the CO 2(m)A being isotopically lighter than CO 2(m)B by ≈10 per mil in the case of gibbsite and nearly 20 per mil in the case of goethite. The reduced partition function ratio of CO2(m)B structure in goethite differs from CO 2(g) by <1 per mil. The predicted fractionation for gibbsite is >10 per mil higher, close to those measured for calcite and aragonite. The surprisingly large difference in the carbon-isotope fractionation factor between the CO2(m)A and CO2(m)B structures within a given mineral suggests that the isotopic signatures of soil (oxy)hydroxide could be heterogeneous. © 2008 by The National Academy of Sciences of the USA.