Ab Initio Study of the Structure and Stability of High-Pressure Iron-Bearing Dolomite

Carbon is subducted into the mantle primarily in the form of metasomatically calcium-enriched basaltic rock, calcified serpentinites and carbonaceous ooze, all of which often contain dolomite. End-member CaMg(CO3)2 dolomite typically breaks down upon compression into two carbonates at 5-6 GPa in the temperature range of 800-1200 K [1]. However, high-pressure X-ray diffraction experiments have recently shown that the presence of iron may be sufficient to stabilize high-pressure dolomite over single-cation carbonates above 35 GPa [2,3]. The structure and equation of state of high-pressure dolomite phases have been debated, creating a need for theoretical calculations. Using density functional theory interfaced with a genetic algorithm that predicts crystal structures (USPEX), we have found a monoclinic phase with space group C2/c. The C2/c structure has a lower energy than previously reported dolomite structures at relevant pressures. It is possible that this phase is not achieved experimentally due to a large energy barrier and a correspondingly large required volume drop, resulting in the transformation to metastable dolomite II. We calculate the equation of state of trigonal dolomite, dolomite III and monoclinic C2/c dolomite to 80 GPa with 0 and 50 mol% CaFe(CO3)2 and compare their enthalpies to single-carbonate assemblages. Although end-member C2/c CaMg(CO3)2 dolomite is not stable relative to single-cation carbonates, C2/c CaMg0.5Fe0.5(CO3)2 is preferred over single-cation carbonates at high pressures. Thus, iron-bearing C2/c dolomite may be an important host phase for carbon in slabs subducted into the lower mantle. [1] Shirasaka, M., et al. (2002) American Mineralogist, 87, 922-930. [2] Mao, Z. et al. (2011) Geophysical Research Letters, 38. [3] Merlini, M. et al. (2012) Proceedings of the National Academy of Sciences, 109, 13509-13514.