Structures of dolomite at ultrahigh pressure and their influence on the deep carbon cycle

Pressure-driven phase transformations on Mg3Ca(CO3)4 huntite carbonate

Magnesium and calcium carbonate minerals are significant reservoirs of Earth's carbon and understanding their behavior under different conditions is crucial for elucidating the mechanisms of deep carbon storage. Huntite, Mg3Ca(CO3)4, is one of the two stable calcium magnesium carbonate phases, together with dolomite. The distinctive cation coordination environment of Ca atoms compared to calcite-type and dolomite structures makes huntite a comparatively less dense phase. Here we examine the behavior of a polycrystalline natural huntite sample under room-temperature compression up to 38 GPa. Synchrotron X-ray diffraction and Raman spectroscopy experiments were carried out in a diamond-anvil cell using He as a highly hydrostatic pressure transmitting medium. XRD results suggest that the initial R32 huntite structure persists up to 21 GPa. The Raman experiment agrees with this result but also suggests the appearance of structural defects from 10 GPa on. Birch-Murnaghan equation of state parameters were fit to the pressure-volume huntite data resulting in zero-pressure volume V0 of 611.7(2) Å3, a bulk modulus B0 of 99.5(11) GPa and a pressure derivative of the bulk modulus of . At 21 GPa, huntite transforms to another trigonal phase (R3), designated here as huntite II. This phase persists up to at least 38 GPa, the maximum pressure reached in this study. The major structural differences between huntite and the huntite-II phase involve the tilting of the [CO3] units with respect to the basal plane and a rotation, which cause a progressive change in the coordination number of the Ca atoms, from 6 to 9. DFT calculations complement the experimental data, providing new insights into the structural response to high-pressure conditions of this magnesium-calcium double carbonate mineral.

Structural Behavior of Minrecordite Carbonate Mineral upon Compression: Effect of Mg → Zn Chemical Substitution in Dolomite-Type Compounds

We report the structural behavior and compressibility of minrecordite, a naturally occurring Zn-rich dolomite mineral, determined using diamond-anvil cell synchrotron X-ray diffraction. Our data show that this rhombohedral CaZn_0.52Mg_0.48(CO₃)₂ carbonate exhibits a highly anisotropic behavior, the c axis being 3.3 times more compressible than the a axis. The axial compressibilities and the equation of state are governed by the compression of the [CaO₆] and [ZnO₆] octahedra, which are the cations in larger proportion in each layer. We observe the existence of a dense polymorph above 13.4(3) GPa using Ne as a pressure-transmitting medium, but the onset pressure of the phase transition decreases with the appearance of deviatoric stresses in nonhydrostatic conditions. Our results suggest that the phase transition observed in minrecordite is strain-induced and that the high-pressure polymorph is intimately related to the CaCO₃-II-type structure. A comparison with other dolomite minerals indicates that the transition pressure decreases when the ratio Zn/Mg in the crystal lattice of pure dolomite is larger than 1. Density functional theory (DFT) calculations predict that a distorted CaCO₃-II-type structure is energetically more stable than dolomite-type CaZn(CO₃)₂ above 10 GPa. However, according to our calculations, the most stable structure above this pressure is a dolomite-V-type phase, a polymorph not observed experimentally.

Probing recycled carbonate in the lower mantle

Whether surficial carbonates can be carried into the Earth's lower mantle is key to global deep carbon cycles but remains poorly understood. New clues from magnesium and zinc isotopic systematics on rocks from deep-rooted mantle plumes are presented and discussed in this Perspective.

Phase stability and dense polymorph of the BaCa(CO3)2 barytocalcite carbonate

The double carbonate BaCa(CO₃)₂ holds potential as host compound for carbon in the Earth’s crust and mantle. Here, we report the crystal structure determination of a high-pressure BaCa(CO₃)₂ phase characterized by single-crystal X-ray diffraction. This phase, named post-barytocalcite, was obtained at 5.7 GPa and can be described by a monoclinic Pm space group. The barytocalcite to post-baritocalcite phase transition involves a significant discontinuous 1.4% decrease of the unit-cell volume, and the increase of the coordination number of 1/4 and 1/2 of the Ba and Ca atoms, respectively. High-pressure powder X-ray diffraction measurements at room- and high-temperatures using synchrotron radiation and DFT calculations yield the thermal expansion of barytocalcite and, together with single-crystal data, the compressibility and anisotropy of both the low- and high-pressure phases. The calculated enthalpy differences between different BaCa(CO₃)₂ polymorphs confirm that barytocalcite is the thermodynamically stable phase at ambient conditions and that it undergoes the phase transition to the experimentally observed post-barytocalcite phase. The double carbonate is significantly less stable than a mixture of the CaCO₃ and BaCO₃ end-members above 10 GPa. The experimental observation of the high-pressure phase up to 15 GPa and 300 ºC suggests that the decomposition into its single carbonate components is kinetically hindered.

Compressibility and Phase Stability of Iron-Rich Ankerite

The structure of the naturally occurring, iron-rich mineral Ca1.08(6)Mg0.24(2)Fe0.64(4)Mn0.04(1)(CO3)2 ankerite was studied in a joint experimental and computational study. Synchrotron X-ray powder diffraction measurements up to 20 GPa were complemented by density functional theory calculations. The rhombohedral ankerite structure is stable under compression up to 12 GPa. A third-order Birch–Murnaghan equation of state yields V0 = 328.2(3) Å3, bulk modulus B0 = 89(4) GPa, and its first-pressure derivative B’0 = 5.3(8)—values which are in good agreement with those obtained in our calculations for an ideal CaFe(CO3)2 ankerite composition. At 12 GPa, the iron-rich ankerite structure undergoes a reversible phase transition that could be a consequence of increasingly non-hydrostatic conditions above 10 GPa. The high-pressure phase could not be characterized. DFT calculations were used to explore the relative stability of several potential high-pressure phases (dolomite-II-, dolomite-III- and dolomite-V-type structures), and suggest that the dolomite-V phase is the thermodynamically stable phase above 5 GPa. A novel high-pressure polymorph more stable than the dolomite-III-type phase for ideal CaFe(CO3)2 ankerite was also proposed. This high-pressure phase consists of Fe and Ca atoms in sevenfold and ninefold coordination, respectively, while carbonate groups remain in a trigonal planar configuration. This phase could be a candidate structure for dense carbonates in other compositional systems.

Crystal Structure of BaCa(CO3)2 Alstonite Carbonate and Its Phase Stability upon Compression

New single-crystal X-ray diffraction experiments and density functional theory (DFT) calculations reveal that the crystal chemistry of the CaO-BaO-CO2 system is more complex than previously thought. We characterized the BaCa(CO3)2 alstonite structure at ambient conditions, which differs from the recently reported crystal structure of this mineral in the stacking of the carbonate groups. This structural change entails the existence of different cation coordination environments. The structural behavior of alstonite at high pressures was studied using synchrotron powder X-ray diffraction data and ab initio calculations up to 19 and 50 GPa, respectively. According to the experiments, above 9 GPa, the alstonite structure distorts into a monoclinic C2 phase derived from the initial trigonal structure. This is consistent with the appearance of imaginary frequencies and geometry relaxation in DFT calculations. Moreover, calculations predict a second phase transition at 24 GPa, which would cause the increase in the coordination number of Ba atoms from 10 to 11 and 12. We determined the equation of state of alstonite (V 0 = 1608(2) Å3, B 0 = 60(3) GPa, B'0 = 4.4(8) from experimental data) and analyzed the evolution of the polyhedral units under compression. The crystal chemistry of alstonite was compared to that of other carbonates and the relative stability of all known BaCa(CO3)2 polymorphs was investigated.

Insights on the deep carbon cycle from the electrical conductivity of carbon-bearing aqueous fluids

The dehydration and decarbonation in the subducting slab are intricately related and the knowledge of the physical properties of the resulting C-H-O fluid is crucial to interpret the petrological, geochemical, and geophysical processes associated with subduction zones. In this study, we investigate the C-H-O fluid released during the progressive devolatilization of carbonate-bearing serpentine-polymorph chrysotile, with in situ electrical conductivity measurements at high pressures and temperatures. The C-H-O fluid produced by carbonated chrysotile exhibits high electrical conductivity compared to carbon-free aqueous fluids and can be an excellent indicator of the migration of carbon in subduction zones. The crystallization of diamond and graphite indicates that the oxidized C-H-O fluids are responsible for the recycling of carbon in the wedge mantle. The carbonate and chrysotile bearing assemblages stabilize dolomite during the devolatilization process. This unique dolomite forming mechanism in chrysotile in subduction slabs may facilitate the transport of carbon into the deep mantle.

Phase Stability and Vibrational Properties of Iron-Bearing Carbonates at High Pressure

The spin transition of iron can greatly affect the stability and various physical properties of iron-bearing carbonates at high pressure. Here, we reported laser Raman measurements on iron-bearing dolomite and siderite at high pressure and room temperature. Raman modes of siderite FeCO3 were investigated up to 75 GPa in the helium (He) pressure medium and up to 82 GPa in the NaCl pressure medium, respectively. We found that the electronic spin-paring transition of iron in siderite occurred sharply at 42–44 GPa, consistent with that in the neon (Ne) pressure medium in our previous study. This indicated that the improved hydrostaticity from Ne to He had minimal effects on the spin transition pressure. Remarkably, the spin crossover of siderite was broadened to 38–48 GPa in the NaCl pressure medium, due to the large deviatoric stress in the sample chamber. In addition, Raman modes of iron-bearing dolomite Ca1.02Mg0.76Fe0.20Mn0.02(CO3)2 were explored up to 58 GPa by using argon as a pressure medium. The sample underwent phase transitions from dolomite-Ⅰ to -Ⅰb phase at ~8 GPa, and then to -Ⅱ at ~15 and -Ⅲb phase at 36 GPa, while no spin transition was observed in iron-bearing dolomite up to 58 GPa. The incorporation of FeCO3 by 20 mol% appeared to marginally decrease the onset pressures of the three phase transitions aforementioned for pure dolomite. At 55–58 GPa, the ν1 mode shifted to a lower frequency at ~1186 cm−1, which was likely associated with the 3 + 1 coordination in dolomite-Ⅲb. These results shed new insights into the nature of iron-bearing carbonates at high pressure.

A new BaCa(CO3)2 polymorph

A new polymorph of the double carbonate BaCa(CO₃)₂, `a C2 phase', has been synthesized. Its structure has been obtained by density-functional-theory-based (DFT-based) model calculations and has been refined by Rietveld analysis of X-ray powder diffraction data. The structure of the new polymorph differs significantly from those of the established polymorphs barytocalcite, paralstonite and alstonite. The unit-cell parameters of the new monoclinic (space group C2) compound are a = 6.6775 (5), b = 5.0982 (4), c = 4.1924 (3) Å, β = 109.259 (1)°. The new compound has been further characterized using Raman spectroscopy. This work shows that earlier studies have misidentified the products of an established synthesis route and that findings based on the incorrect identification of the synthesis product concerning the suitability of barytocalcite as a matrix for the retention of radioactive isotopes will need to be reconsidered.

Stability of iron-bearing carbonates in the deep Earth's interior

The presence of carbonates in inclusions in diamonds coming from depths exceeding 670 km are obvious evidence that carbonates exist in the Earth's lower mantle. However, their range of stability, crystal structures and the thermodynamic conditions of the decarbonation processes remain poorly constrained. Here we investigate the behaviour of pure iron carbonate at pressures over 100 GPa and temperatures over 2,500 K using single-crystal X-ray diffraction and Mössbauer spectroscopy in laser-heated diamond anvil cells. On heating to temperatures of the Earth's geotherm at pressures to ∼50 GPa FeCO₃ partially dissociates to form various iron oxides. At higher pressures FeCO₃ forms two new structures-tetrairon(III) orthocarbonate Fe₄³⁺C₃O₁₂, and diiron(II) diiron(III) tetracarbonate Fe₂²⁺Fe₂³⁺C₄O₁₃, both phases containing CO₄ tetrahedra. Fe₄C₄O₁₃ is stable at conditions along the entire geotherm to depths of at least 2,500 km, thus demonstrating that self-oxidation-reduction reactions can preserve carbonates in the Earth's lower mantle.

Time differentiated nuclear resonance spectroscopy coupled with pulsed laser heating in diamond anvil cells

Developments in pulsed laser heating applied to nuclear resonance techniques are presented together with their applications to studies of geophysically relevant materials. Continuous laser heating in diamond anvil cells is a widely used method to generate extreme temperatures at static high pressure conditions in order to study the structure and properties of materials found in deep planetary interiors. The pulsed laser heating technique has advantages over continuous heating, including prevention of the spreading of heated sample and/or the pressure medium and, thus, a better stability of the heating process. Time differentiated data acquisition coupled with pulsed laser heating in diamond anvil cells was successfully tested at the Nuclear Resonance beamline (ID18) of the European Synchrotron Radiation Facility. We show examples applying the method to investigation of an assemblage containing ε-Fe, FeO, and Fe3C using synchrotron Mössbauer source spectroscopy, FeCO3 using nuclear inelastic scattering, and Fe2O3 using nuclear forward scattering. These examples demonstrate the applicability of pulsed laser heating in diamond anvil cells to spectroscopic techniques with long data acquisition times, because it enables stable pulsed heating with data collection at specific time intervals that are synchronized with laser pulses.

High-pressure orthorhombic ferromagnesite as a potential deep-mantle carbon carrier

Knowledge of the physical and chemical properties of candidate deep-carbon carriers such as ferromagnesite [(Mg,Fe)CO₃] at high pressure and temperature of the deep mantle is necessary for our understanding of deep-carbon storage as well as the global carbon cycle of the planet. Previous studies have reported very different scenarios for the (Mg,Fe)CO₃ system at deep-mantle conditions including the chemical dissociation to (Mg,Fe)O+CO₂, the occurrence of the tetrahedrally-coordinated carbonates based on CO₄ structural units, and various high-pressure phase transitions. Here we have studied the phase stability and compressional behavior of (Mg,Fe)CO₃ carbonates up to relevant lower-mantle conditions of approximately 120 GPa and 2400 K. Our experimental results show that the rhombohedral siderite (Phase I) transforms to an orthorhombic phase (Phase II with Pmm2 space group) at approximately 50 GPa and 1400 K. The structural transition is likely driven by the spin transition of iron accompanied by a volume collapse in the Fe-rich (Mg,Fe)CO₃ phases; the spin transition stabilizes the high-pressure phase II at much lower pressure conditions than its Mg-rich counterpart. It is conceivable that the low-spin ferromagnesite phase II becomes a major deep-carbon carrier at the deeper parts of the lower mantle below 1900 km in depth.

Mantle-slab interaction and redox mechanism of diamond formation

Subduction tectonics imposes an important role in the evolution of the interior of the Earth and its global carbon cycle; however, the mechanism of the mantle-slab interaction remains unclear. Here, we demonstrate the results of high-pressure redox-gradient experiments on the interactions between Mg-Ca-carbonate and metallic iron, modeling the processes at the mantle-slab boundary; thereby, we present mechanisms of diamond formation both ahead of and behind the redox front. It is determined that, at oxidized conditions, a low-temperature Ca-rich carbonate melt is generated. This melt acts as both the carbon source and crystallization medium for diamond, whereas at reduced conditions, diamond crystallizes only from the Fe-C melt. The redox mechanism revealed in this study is used to explain the contrasting heterogeneity of natural diamonds, as seen in the composition of inclusions, carbon isotopic composition, and nitrogen impurity content.

Single-crystal diffraction at the Extreme Conditions beamline P02.2: procedure for collecting and analyzing high-pressure single-crystal data

Fast detectors employed at third-generation synchrotrons have reduced collection times significantly and require the optimization of commercial as well as customized software packages for data reduction and analysis. In this paper a procedure to collect, process and analyze single-crystal data sets collected at high pressure at the Extreme Conditions beamline (P02.2) at PETRA III, DESY, is presented. A new data image format called `Esperanto' is introduced that is supported by the commercial software package CrysAlis(Pro) (Agilent Technologies UK Ltd). The new format acts as a vehicle to transform the most common area-detector data formats via a translator software. Such a conversion tool has been developed and converts tiff data collected on a Perkin Elmer detector, as well as data collected on a MAR345/555, to be imported into the CrysAlis(Pro) software. In order to demonstrate the validity of the new approach, a complete structure refinement of boron-mullite (Al5BO9) collected at a pressure of 19.4 (2) GPa is presented. Details pertaining to the data collections and refinements of B-mullite are presented.