Stability of calcium and magnesium carbonates at lower mantle thermodynamic conditions

The highest pressure form of the major Earth-forming mantle silicate is MgSiO3 post-perovskite (PPv). Understanding the fate of PPv at TPa pressures is the first step for understanding the mineralogy of super-Earths-type exoplanets, arguably the most interesting for their similarities with Earth. Modeling their internal structure requires knowledge of stable mineral phases, their properties under compression, and major element abundances. Several studies of PPv under extreme ...

Thermal conductivity ($\kappa$) of mantle minerals is a fundamental property in geodynamic modeling. It controls the style of mantle convection and the time scale of the mantle and core cooling. Cubic CaSiO$_3$ perovskite (CaPv) is the third most abundant mineral in the lower mantle (LM) (7 vol%). However, despite its importance, no theoretical or experimental estimate of CaPv's $\kappa$ exists. Theoretical investigations of its properties are challenging because of its stron...

We investigated the stability of polymeric CO2 over a wide range of pressures, temperatures, and chemical environments. We find that the I-42d polymeric structure, consisting of a three-dimensional network of corner sharing CO4 tetrahedra, forms at 40-140 GPa and from a CO-N2 mixture at 39 GPa. An exceptional stability field of 0 to 286 GPa and 100 to 2500 K is documented for this structure, making it a viable candidate for planetary interiors. The stability of the tetrahedra...

Using evolutionary crystal structure prediction algorithm USPEX, we showed that at pressures of the Earth's lower mantle CaAl2O4 is the only stable calcium aluminate. At pressures above 7.0 GPa it has the CaFe2O4-type structure and space group Pnma. This phase is one of prime candidate aluminous phases in the lower mantle of the Earth. We show that at low pressures 5CaO * 3Al2O3 (C5A3) with space group Cmc21, CaAl4O7 (C2/c) and CaAl2O4 (P21/m) structures are stable at pressur...

We report lattice thermal conductivities of MgO and MgSiO3 in the perovskite and post-perovskite structures at conditions of the Earth's lower mantle, obtained from equilibrium molecular dynamics simulations. Using an advanced ionic interaction potential, the full conductivity tensor was calculated by means of the Green-Kubo method, and the conductivity of MgSiO3 post-perovskite was found to be significantly anisotropic. The thermal conductivities of all three phases were par...

We have determined the lattice dynamics of MgSiO$_3$ perovskite (bridgmanite) by a combination of single-crystal inelastic x-ray scattering and ab initio calculations. We observe a remarkable agreement between experiment and theory, and provide accurate results for phonon dispersion relations, phonon density of states and the full elasticity tensor. The present work constitutes an important milestone to extend this kind of combined studies to extreme conditions of pressure an...

To address the effects of lattice anharmonicity across the perovskite to post-perovskite transition in MgSiO$_3$, we conduct calculations using the phonon quasiparticle (PHQ) approach. The PHQ is based on \textit{ab initio} molecular dynamics and, in principle, captures full anharmonicity. Free energies in the thermodynamic limit ($N \rightarrow \infty$) are computed using temperature-dependent quasiparticle dispersions within the phonon gas model. Systematic results on anhar...

MgO is an abundant mineral in the rocky mantle of terrestrial planets, but its melting behaviors remain enigmatic. Here we introduce a simple theoretical model to investigate the B1-liquid transition of MgO up to 370 GPa. Vibrational free energies of B1-MgO are fully computed by the moment recurrence technique in quantum statistical physics. On that basis, we associate the melting temperature with the isothermal bulk modulus via the work-heat equivalence principle. This strat...

Carbonate precipitation in oceans is essential for the carbonate-silicate cycle (inorganic carbon cycle) to maintain temperate climates. By considering the thermodynamics of carbonate chemistry, we demonstrate that the ocean pH decreases by approximately 0.5 for a factor of 10 increase in the atmospheric carbon dioxide content. The upper and lower limits of ocean pH are within 1-4 of each other, where the upper limit is buffered by carbonate precipitation and defines the ocea...

The family of binary compounds composed of magnesium and silicon is rather rich. In addition to the well-known magnesium silicide Mg$_2$Si, other compounds, including MgSi$_2$, Mg$_4$Si$_7$, Mg$_5$Si$_6$, MgSi, and Mg$_9$Si$_5$, have also been identified and/or proposed in precipitated Al-Mg-Si solid solutions. Nevertheless, computational studies show that only Mg$_2$Si is thermodynamically stable at ambient conditions while certain non-zero hydrostatic pressure can stabilize...