-
The miscibility of hydrogen and water in planetary atmospheres and interiors
Authors:
Akash Gupta,
Lars Stixrude,
Hilke E. Schlichting
Abstract:
Many planets in the solar system and across the galaxy have hydrogen-rich atmospheres overlying more heavy element-rich interiors with which they interact for billions of years. Atmosphere-interior interactions are thus crucial to understanding the formation and evolution of these bodies. However, this understanding is still lacking in part because the relevant pressure-temperature conditions are…
▽ More
Many planets in the solar system and across the galaxy have hydrogen-rich atmospheres overlying more heavy element-rich interiors with which they interact for billions of years. Atmosphere-interior interactions are thus crucial to understanding the formation and evolution of these bodies. However, this understanding is still lacking in part because the relevant pressure-temperature conditions are extreme. We conduct molecular dynamics simulations based on Density Functional Theory to investigate how hydrogen and water interact over a wide range of pressure and temperature, encompassing the interiors of Neptune-sized and smaller planets. We determine the critical curve at which a single homogeneous phase exsolves into two separate, hydrogen-rich and water-rich phases, finding good agreement with existing experimental data. We find that the temperature along the critical curve increases with increasing pressure and shows the influence of a change in fluid structure from molecular to atomic near 30 GPa and 3000 K, which may impact magnetic field generation. The internal temperatures of many exoplanets, including TOI-270 d and K2-18 b may lie entirely above the critical curve: the envelope is expected to consist of a single homogeneous hydrogen-water fluid, that is much less susceptible to atmospheric loss as compared with a pure hydrogen envelope. As planets cool, they cross the critical curve, leading to rainout of water-rich fluid and an increase in internal luminosity. Compositions of the resulting outer, hydrogen-rich, and inner, water-rich envelopes depend on age and instellation and are governed by thermodynamics. Rainout of water may be occurring in Uranus and Neptune at present.
△ Less
Submitted 5 July, 2024;
originally announced July 2024.
-
Nonlinearity of the post-spinel transition and its expression in slabs and plumes worldwide
Authors:
Junjie Dong,
Rebecca A. Fischer,
Lars Stixrude,
Matthew C. Brennan,
Kierstin Daviau,
Terry-Ann Suer,
Katlyn M. Turner,
Yue Meng,
Vitali B. Prakapenka
Abstract:
At the interface of Earth's upper and lower mantle, the post-spinel transition boundary controls the dynamics and morphologies of downwelling slabs and upwelling plumes, and its Clapeyron slope is hence one of the most important constraints on mantle convection. In this study, we reported a new in situ experimental dataset on phase stability in Mg$_{2}$SiO$_{4}$ at mantle transition zone pressures…
▽ More
At the interface of Earth's upper and lower mantle, the post-spinel transition boundary controls the dynamics and morphologies of downwelling slabs and upwelling plumes, and its Clapeyron slope is hence one of the most important constraints on mantle convection. In this study, we reported a new in situ experimental dataset on phase stability in Mg$_{2}$SiO$_{4}$ at mantle transition zone pressures from laser-heated diamond anvil cell experiments, along with a compilation of corrected in situ experimental datasets from the literature. We presented a machine learning framework for high-pressure phase diagram determination and focused on its application to constrain the location and Clapeyron slope of the post-spinel transition: ringwoodite $\leftrightarrow$ bridgmanite + periclase. We found that the post-spinel boundary is nonlinear and its Clapeyron slope varies locally from $-2.3_{-1.4}^{+0.6}$ MPa/K at 1900 K, to $-1.0_{-1.7}^{+1.3}$ MPa/K at 1700 K, and to $0.0_{-2.0}^{+1.7}$ MPa/K at 1500 K. We applied the temperature-dependent post-spinel Clapeyron slope to estimate its lateral variation across the "660-km" seismic discontinuity in subducting slabs and hotspot-associated plumes worldwide, as well as the ambient mantle. We found that, in the present-day mantle, the average post spinel Clapeyron slope in the plumes is three times more negative than that in slabs, and we then discussed the effects of a nonlinear post-spinel transition on the dynamics of Earth's mantle.
△ Less
Submitted 18 August, 2022;
originally announced August 2022.
-
Heat and charge transport in H$_2$O at ice-giant conditions from ab initio molecular dynamics simulations
Authors:
Federico Grasselli,
Lars Stixrude,
Stefano Baroni
Abstract:
The impact of the inner structure and thermal history of planets on their observable features, such as luminosity or magnetic field, crucially depends on the poorly known heat and charge transport properties of their internal layers. The thermal and electric conductivities of different phases of water (liquid, solid, and super-ionic) occurring in the interior of ice giant planets, such as Uranus o…
▽ More
The impact of the inner structure and thermal history of planets on their observable features, such as luminosity or magnetic field, crucially depends on the poorly known heat and charge transport properties of their internal layers. The thermal and electric conductivities of different phases of water (liquid, solid, and super-ionic) occurring in the interior of ice giant planets, such as Uranus or Neptune, are evaluated from equilibrium ab initio molecular dynamics, leveraging recent progresses in the theory and data analysis of transport in extended systems. The implications of our findings on the evolution models of the ice giants are briefly discussed
△ Less
Submitted 14 July, 2020; v1 submitted 27 March, 2020;
originally announced March 2020.
-
Spin crossover in ferropericlase from first-principles molecular dynamics
Authors:
Eero Holmstrom,
Lars Stixrude
Abstract:
Ferropericlase, (Mg,Fe)O, is the second-most abundant mineral of the Earth's lower mantle. With increasing pressure, the Fe ions in the material begin to collapse from a magnetic to non-magnetic spin state. We present a finite-temperature first-principles phase diagram of this spin crossover, finding a broad pressure range with coexisting magnetic and non-magnetic ions due to favorable enthalpy of…
▽ More
Ferropericlase, (Mg,Fe)O, is the second-most abundant mineral of the Earth's lower mantle. With increasing pressure, the Fe ions in the material begin to collapse from a magnetic to non-magnetic spin state. We present a finite-temperature first-principles phase diagram of this spin crossover, finding a broad pressure range with coexisting magnetic and non-magnetic ions due to favorable enthalpy of mixing of the two. Furthermore, we find the electrical conductivity of the mineral to reach semi-metallic values inside the Earth.
△ Less
Submitted 13 January, 2015;
originally announced January 2015.
-
Evidence of local magnetic order in hcp iron from Raman mode splitting
Authors:
Gerd Steinle-Neumann,
Lars Stixrude,
R. E. Cohen,
Boris Kiefer
Abstract:
Experimental measurements of Raman spectra for hcp iron at high pressure show two modes over a considerable pressure range in contrast to the prediction of one doubly degenerate mode for the hcp lattice. We use density functional theory to investigate the influence of magnetic order on the Raman active modes of hcp iron. We find an antiferromagnetic state that lifts the degeneracy of the transve…
▽ More
Experimental measurements of Raman spectra for hcp iron at high pressure show two modes over a considerable pressure range in contrast to the prediction of one doubly degenerate mode for the hcp lattice. We use density functional theory to investigate the influence of magnetic order on the Raman active modes of hcp iron. We find an antiferromagnetic state that lifts the degeneracy of the transverse optical mode, and yields stable antiferromagnetic moments up to approximately 60 GPa (55 Bohr$^3$). The resulting frequencies of the two transverse optical modes are in good agreement with the experimental Raman shifts, lending support to the existence of local antiferromagnetic order in hcp iron.
△ Less
Submitted 26 November, 2001;
originally announced November 2001.
-
Importance of Magnetism in Phase Stability, Equations of State, and Elasticity
Authors:
R. E. Cohen,
S. Gramsch,
G. Steinle-Neumann,
L. Stixrude
Abstract:
The effects of magnetism on high pressure properties of transition metals and transition metal compounds can be quite important. In the case of Fe, magnetism is responsible for stability of the body-centered cubic (bcc) phase at ambient conditions, and the large thermal expansivity in face-centered cubic (fcc) iron, and also has large effects on the equation of state and elasticity of hexagonal…
▽ More
The effects of magnetism on high pressure properties of transition metals and transition metal compounds can be quite important. In the case of Fe, magnetism is responsible for stability of the body-centered cubic (bcc) phase at ambient conditions, and the large thermal expansivity in face-centered cubic (fcc) iron, and also has large effects on the equation of state and elasticity of hexagonal close-packed (hcp) iron. In transition metal oxides, local magnetic moments are responsible for their insulating behavior. LDA+U results are presented for CoO and FeO, and predictions are made for high pressure metallization. The inclusion of a local Coulomb repulsion, U, greatly inhibits the high-spin low-spin transitions found with conventional exchange-correlation functionals (i.e. generalized gradient corrections, GGA). We discuss theory and computations for the effects of magnetism on high pressure cohesive properties.
△ Less
Submitted 24 July, 2002; v1 submitted 1 October, 2001;
originally announced October 2001.
-
Absence of lattice strain anomalies at the electronic topological transition in zinc at high pressure
Authors:
Gerd Steinle-Neumann,
Lars Stixrude,
Ronald E. Cohen
Abstract:
High pressure structural distortions of the hexagonal close packed (hcp) element zinc have been a subject of controversy. Earlier experimental results and theory showed a large anomaly in lattice strain with compression in zinc at about 10 GPa which was explained theoretically by a change in Fermi surface topology. Later hydrostatic experiments showed no such anomaly, resulting in a discrepancy…
▽ More
High pressure structural distortions of the hexagonal close packed (hcp) element zinc have been a subject of controversy. Earlier experimental results and theory showed a large anomaly in lattice strain with compression in zinc at about 10 GPa which was explained theoretically by a change in Fermi surface topology. Later hydrostatic experiments showed no such anomaly, resulting in a discrepancy between theory and experiment. We have computed the compression and lattice strain of hcp zinc over a wide range of compressions using the linearized augmented plane wave (LAPW) method paying special attention to k-point convergence. We find that the behavior of the lattice strain is strongly dependent on k-point sampling, and with large k-point sets the previously computed anomaly in lattice parameters under compression disappears, in agreement with recent experiments.
△ Less
Submitted 16 October, 2000; v1 submitted 30 May, 2000;
originally announced May 2000.
-
First-Principles Elastic Constants for the hcp Transition Metals Fe, Co, and Re at High Pressure
Authors:
G. Steinle-Neumann,
L. Stixrude,
R. E. Cohen
Abstract:
The elastic constant tensors for the hcp phases of three transition metals (Co, Re, and Fe) are computed as functions of pressure using the Linearized Augmented Plane Wave method with both the local density and generalized gradient approximations. Spin-polarized states are found to be stable for Co (ferromagnetic) and Fe (antiferromagnetic at low pressure). The elastic constants of Co and Re are…
▽ More
The elastic constant tensors for the hcp phases of three transition metals (Co, Re, and Fe) are computed as functions of pressure using the Linearized Augmented Plane Wave method with both the local density and generalized gradient approximations. Spin-polarized states are found to be stable for Co (ferromagnetic) and Fe (antiferromagnetic at low pressure). The elastic constants of Co and Re are compared to experimental measurements near ambient conditions and excellent agreement is found. Recent measurements of the lattice strain in high pressure experiments when interpreted in terms of elastic constants for Re and Fe are inconsistent with the calculated moduli.
△ Less
Submitted 29 April, 1999;
originally announced April 1999.
-
Prediction of Phase Transition in CaSiO$_3$ Perovskite and Implications for Lower Mantle Structure
Authors:
Lars Stixrude,
Ronald E. Cohen,
Rici Yu,
Henry Krakauer
Abstract:
First principles linear response calculations are used to investigate the lattice dynamics of what is thought to be the third most abundant phase in the lower mantle, CaSiO_3 perovskite. The commonly assumed cubic structure (Pm3m) is found to be dynamically unstable at all pressures, exhibiting unstable modes along the Brillouin zone edge from the M-point to the R-point. Based on these results,…
▽ More
First principles linear response calculations are used to investigate the lattice dynamics of what is thought to be the third most abundant phase in the lower mantle, CaSiO_3 perovskite. The commonly assumed cubic structure (Pm3m) is found to be dynamically unstable at all pressures, exhibiting unstable modes along the Brillouin zone edge from the M-point to the R-point. Based on these results, we predict that the ground state structure of CaSiO_3 perovskite is a distorted phase with lower than cubic symmetry. We predict that a phase transition occurs in CaSiO_3 perovskite within the earth's lower mantle from the low temperature distorted phase to the cubic phase at high temperature. The predicted phase transition provides a possible explanation of some of the seismological observations of reflective features within the lower mantle.
△ Less
Submitted 31 October, 1996;
originally announced October 1996.