skip to main content

Title: Thermal equation of state of post-aragonite CaCO3-Pmmn
Abstract Calcium carbonate (CaCO3) is one of the most abundant carbonates on Earth's surface and transports carbon to Earth's interior via subduction. Although some petrological observations support the preservation of CaCO3 in cold slabs to lower mantle depths, the geophysical properties and stability of CaCO3 at these depths are not known, due in part to complicated polymorphic phase transitions and lack of constraints on thermodynamic properties. Here we measured thermal equation of state of CaCO3-Pmmn, the stable polymorph of CaCO3 through much of the lower mantle, using synchrotron X-ray diffraction in a laser-heated diamond-anvil cell up to 75 GPa and 2200 K. The room-temperature compression data for CaCO3-Pmmn are fit with third-order Birch-Murnaghan equation of state, yielding KT0 = 146.7 (±1.9) GPa and K′0 = 3.4(±0.1) with V0 fixed to the value determined by ab initio calculation, 97.76 Å3. High-temperature compression data are consistent with zero-pressure thermal expansion αT = a0 + a1T with a0 = 4.3(±0.3)×10-5 K-1, a1 = 0.8(±0.2)×10-8 K-2, temperature derivative of the bulk modulus (∂KT/∂T)P = –0.021(±0.001) GPa/K; the Grüneisen parameter γ0 = 1.94(±0.02), and the volume independent constant q = 1.9(±0.3) at a fixed Debye temperature θ0 = 631 K predicted via ab initio calculation. more » Using these newly determined thermodynamic parameters, the density and bulk sound velocity of CaCO3-Pmmn and (Ca,Mg)-carbonate-bearing eclogite are quantitatively modeled from 30 to 80 GPa along a cold slab geotherm. With the assumption that carbonates are homogeneously mixed into the slab, the results indicate the presence of carbonates in the subducted slab is unlikely to be detected by seismic observations, and the buoyancy provided by carbonates has a negligible effect on slab dynamics. « less
; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
American Mineralogist
Page Range or eLocation-ID:
1365 to 1374
Sponsoring Org:
National Science Foundation
More Like this
  1. Fe‐Al‐bearing bridgmanite may be the dominant host for ferric iron in Earth's lower mantle. Here we report the synthesis of (Mg0.5Fe3+0.5)(Al0.5Si0.5)O3 bridgmanite (FA50) with the highest Fe3+‐Al3+ coupled substitution known to date. X‐ray diffraction measurements showed that at ambient conditions the FA50 adopted the LiNbO3 structure. Upon compression at room temperature to 18 GPa, it transformed back into the bridgmanite structure, which remained stable up to 102 GPa and 2600 K. Fitting Birch‐Murnaghan equation of state of FA50 bridgmanite yields V 0 = 172.1(4) Å3, K 0 = 229(4) GPa with K 0′ = 4(fixed). The calculated bulk sound velocity of the FA50 bridgmanite is ~7.7% lower than MgSiO3 bridgmanite, mainly because the presence of ferric iron increases the unit‐cell mass by 15.5%. This difference likely represents the upper limit of sound velocity anomaly introduced by Fe3+‐Al3+ substitution. X‐ray emission and synchrotron Mössbauer spectroscopy measurements showed that after laser annealing ~6% of Fe3+ cations exchanged with Al3+ and underwent the high‐spin to low‐spin transition at 59 GPa. The low‐spin proportion of Fe3+ increased gradually with pressure and reached 17‐31% at 80 GPa. Since the cation exchange and spin transition in this Fe3+‐Al3+‐enriched bridgmanite do not cause resolvable unit‐cell volume reduction,more »and the increase of low‐spin Fe3+ fraction with pressure occurs gradually, the spin transition would not produce a distinct seismic signature in the lower mantle. However, it may influence iron partitioning and isotopic fractionation, thus introducing chemical heterogeneity in the lower mantle.« less
  2. Thermoelastic behavior of transition metal boride Os2B3 was studied under quasi-hydrostatic and isothermal conditions in a Paris-Edinburgh cell to 5.4 GPa and 1273 K. In-situ Energy Dispersive X-ray diffraction was used to determine interplanar spacings of the hexagonal crystal structure and thus the volume and axial compression. P-V-T data were fitted to a 3rd Order Birch-Murnaghan equation of state with a temperature modification to determine thermal elastic constants. The bulk modulus was shown to be K0 = 402 ± 21 GPa when the first pressure derivative was held to K0’ = 4.0 from the room temperature P-V curve. Under a quadratic fit α=α_0+α_1 T-α_2 T^(-2), the thermal expansion coefficients were determined to be α_0=1.862×10^(-5) K-1, α_1=0.841×10^(-9) K-2, and α_2=-0.525 K. Density functional theory (DFT) with the quasi-harmonic approximation (QHA) were further employed to study Os2B3, including its P-V-T curves, phonon spectra, bulk modulus, specific heat, thermal expansion, and the Grüneisen parameter. A good agreement between the first-principle theory and experimental observations was achieved, highlighting the success of the Armiento-Mattsson 2005 generalized gradient approximation functional employed in this study and QHA for describing thermodynamic properties of Os2B3.
  3. Abstract In this study, we have investigated the crystal structure and equation of state of tetragonal CaSiO3-perovskite up to 200 GPa using synchrotron X-ray diffraction in laser-heated diamond-anvil cells. X-ray diffraction patterns of the quenched CaSiO3-perovskite above 148 GPa clearly show that 200, 211, and 220 peaks of the cubic phase split into 004+220, 204+312, and 224+400 peak pairs, respectively, in the tetragonal structure, and their calculated full-width at half maximum (FWHM) exhibits a substantial increase with pressure. The distribution of diffraction peaks suggests that the tetragonal CaSiO3-perovskite most likely has an I4/mcm space group at 300 K between 148 and 199 GPa, although other possibilities might still exist. Using the Birch-Murnaghan equations, we have determined the equation of state of tetragonal CaSiO3-perovskite, yielding the bulk modulus K0T = 227(21) GPa with the pressure derivative of the bulk modulus, K0T′ = 4.0(3). Modeled sound velocities at 580 K and around 50 GPa using our results and literature values show the difference in the compressional (VP) and shear-wave velocity (VS) between the tetragonal and cubic phases to be 5.3 and 6.7%, respectively. At ~110 GPa and 1000 K, this phase transition leads to a 4.3 and 9.1% jump in VP andmore »VS, respectively. Since the addition of Ti can elevate the transition temperature, the transition from the tetragonal to cubic phase may have a seismic signature compatible with the observed mid-lower mantle discontinuity around the cold subduction slabs, which needs to be explored in future studies.« less
  4. Both seismic observations of dense low shear velocity regions and models of magma ocean crystallization and mantle dynamics support enrichment of iron in Earth’s lowermost mantle. Physical properties of iron-rich lower mantle heterogeneities in the modern Earth depend on distribution of iron between coexisting lower mantle phases (Mg,Fe)O magnesiowüstite, (Mg,Fe)SiO3 bridgmanite, and (Mg,Fe)SiO3 post-perovskite. The partitioning of iron between these phases was investigated in synthetic ferrous-iron-rich olivine compositions (Mg0.55Fe0.45)2SiO4 and (Mg0.28Fe0.72)2SiO4 at lower mantle conditions ranging from 33–128 GPa and 1900–3000 K in the laser-heated diamond anvil cell. The resulting phase assemblages were characterized by a combination of in situ X-ray diffraction and ex situ transmission electron microscopy. The exchange coefficient between bridgmanite and magnesiowüstite decreases with pressure and bulk Fe# and increases with temperature. Thermodynamic modeling determines that incorporation and partitioning of iron in bridgmanite are explained well by excess volume associated with Mg-Fe exchange. Partitioning results are used to model compositions and densities of mantle phase assemblages as a function of pressure, FeO-content and SiO2-content. Unlike average mantle compositions, iron-rich compositions in the mantle exhibit negative dependence of density on SiO2-content at all mantle depths, an important finding for interpretation of deep lower mantle structures.
  5. Deeply subducted carbonates likely cause low-degree melting of the upper mantle and thus play an important role in the deep carbon cycle. However, direct seismic detection of carbonate-induced partial melts in the Earth’s interior is hindered by our poor knowledge on the elastic properties of carbonate melts. Here we report the first experimentally determined sound velocity and density data on dolomite melt up to 5.9 GPa and 2046 K by in-situ ultrasonic and sink-float techniques, respectively, as well as first-principles molecular dynamics simulations of dolomite melt up to 16 GPa and 3000 K. Using our new elasticity data, the calculated V P /V S ratio of the deep upper mantle (∼180–330 km) with a small amount of carbonate-rich melt provides a natural explanation for the elevated V P /V S ratio of the upper mantle from global seismic observations, supporting the pervasive presence of a low-degree carbonate-rich partial melt (∼0.05%) that is consistent with the volatile-induced or redox-regulated initial melting in the upper mantle as argued by petrologic studies. This carbonate-rich partial melt region implies a global average carbon (C) concentration of 80–140 ppm. by weight in the deep upper mantle source region, consistent with the mantle carbon content determinedmore »from geochemical studies.« less