-
Diamond formation from hydrocarbon mixtures in planets
Authors:
Bingqing Cheng,
Sebastien Hamel,
Mandy Bethkenhagen
Abstract:
Hydrocarbon mixtures are extremely abundant in the Universe, and diamond formation from them can play a crucial role in shaping the interior structure and evolution of planets. With first-principles accuracy, we first estimate the diamond nucleation rate in pure liquid carbon, and then reveal the nature of chemical bonding in hydrocarbons at extreme conditions. We finally establish the pressure-te…
▽ More
Hydrocarbon mixtures are extremely abundant in the Universe, and diamond formation from them can play a crucial role in shaping the interior structure and evolution of planets. With first-principles accuracy, we first estimate the diamond nucleation rate in pure liquid carbon, and then reveal the nature of chemical bonding in hydrocarbons at extreme conditions. We finally establish the pressure-temperature phase boundary where diamond can form from hydrocarbon mixtures with different atomic fractions of carbon. Notably, we find a depletion zone at pressures above 200~GPa and temperatures below 3000~K-3500~K where diamond formation is thermodynamically favorable regardless of the carbon atomic fraction, due to a phase separation mechanism. The cooler condition of the interior of Neptune compared to Uranus means that the former is much more likely to contain the depletion zone. Our findings can help explain the dichotomy of the two ice giants manifested by the low luminosity of Uranus, and lead to a better understanding of (exo-)planetary formation and evolution.
△ Less
Submitted 2 August, 2022; v1 submitted 6 July, 2022;
originally announced July 2022.
-
Predicting the phase behaviors of superionic water at planetary conditions
Authors:
Bingqing Cheng,
Mandy Bethkenhagen,
Chris J. Pickard,
Sebastien Hamel
Abstract:
Most water in the universe may be superionic, and its thermodynamic and transport properties are crucial for planetary science but difficult to probe experimentally or theoretically. We use machine learning and free energy methods to overcome the limitations of quantum mechanical simulations, and characterize hydrogen diffusion, superionic transitions, and phase behaviors of water at extreme condi…
▽ More
Most water in the universe may be superionic, and its thermodynamic and transport properties are crucial for planetary science but difficult to probe experimentally or theoretically. We use machine learning and free energy methods to overcome the limitations of quantum mechanical simulations, and characterize hydrogen diffusion, superionic transitions, and phase behaviors of water at extreme conditions. We predict that a close-packed superionic phase with mixed stacking is stable over a wide temperature and pressure range, while a body-centered cubic phase is only thermodynamically stable in a small window but is kinetically favored. Our phase boundaries, which are consistent with the existing-albeit scarce-experimental observations, help resolve the fractions of insulating ice, different superionic phases, and liquid water inside of ice giants.
△ Less
Submitted 16 March, 2021;
originally announced March 2021.
-
Metallization of Shock-Compressed Liquid Ammonia
Authors:
A. Ravasio,
M. Bethkenhagen,
J. -A. Hernandez,
A. Benuzzi-Mounaix,
F. Datchi,
M. French,
M. Guarguaglini,
F. Lefevre,
S. Ninet,
R. Redmer,
T. Vinci
Abstract:
Ammonia is predicted to be one of the major components in the depths of the ice giant planets Uranus and Neptune. Their dynamics, evolution, and interior structure are insufficiently understood and models rely imperatively on data for equation of state and transport properties. Despite its great significance, the experimentally accessed region of the ammonia phase diagram today is still very limit…
▽ More
Ammonia is predicted to be one of the major components in the depths of the ice giant planets Uranus and Neptune. Their dynamics, evolution, and interior structure are insufficiently understood and models rely imperatively on data for equation of state and transport properties. Despite its great significance, the experimentally accessed region of the ammonia phase diagram today is still very limited in pressure and temperature. Here we push the probed regime to unprecedented conditions, up to $\sim$350 GPa and $\sim$40000 K. Along the Hugoniot, the temperature measured as a function of pressure shows a subtle change in slope at $\sim$7000 K and $\sim$90 GPa, in agreement with ab initio simulations we have performed. This feature coincides with the gradual transition from a molecular liquid to a plasma state. Additionally, we performed reflectivity measurements, providing the first experimental evidence of electronic conduction in high-pressure ammonia. Shock reflectance continuously rises with pressure above 50 GPa and reaches saturation values above 120 GPa. Corresponding electrical conductivity values are up to 1 order of magnitude higher than in water in the 100 GPa regime, with possible significant contributions of the predicted ammonia-rich layers to the generation of magnetic dynamos in ice giant interiors.
△ Less
Submitted 17 January, 2021;
originally announced January 2021.
-
Sensitivity of stellar physics to the equation of state
Authors:
D. C. Swift,
T. Lockard,
M. Bethkenhagen,
A. Kritcher,
S. Hamel,
D. Dearborn
Abstract:
The formation and evolution of stars depends on various physical aspects of stellar matter, including the equation of state (EOS) and transport properties. Although often dismissed as `ideal gas-like' and therefore simple, states occurring in stellar matter are dense plasmas, and the EOS has not been established precisely. EOS constructed using multi-physics approaches found necessary for laborato…
▽ More
The formation and evolution of stars depends on various physical aspects of stellar matter, including the equation of state (EOS) and transport properties. Although often dismissed as `ideal gas-like' and therefore simple, states occurring in stellar matter are dense plasmas, and the EOS has not been established precisely. EOS constructed using multi-physics approaches found necessary for laboratory studies of warm dense matter give significant variations in stellar regimes, and vary from the EOS commonly used in simulations of the formation and evolution of stars. We have investigated the sensitivity of such simulations to variations in the EOS, for sun-like and low-mass stars, and found a strong sensitivity of the lifetime of the Sun and of the lower luminosity limit for red dwarfs. We also found a significant sensitivity in the lower mass limit for red dwarfs. Simulations of this type are also used for other purposes in astrophysics, including the interpretation of absolute magnitude as mass, the conversion of inferred mass distribution to the initial mass function using predicted lifetimes, simulations of star formation from nebulae, simulations of galactic evolution, and the baryon census used to bound the exotic contribution to dark matter. Although many of the sensitivities of stellar physics to the EOS are large, some of the inferred astrophysical quantities are also constrained by independent measurements, although the constraints may be indirect and non-trivial. However, it may be possible to use such measurements to constrain the EOS more than presently possible by established plasma theory.
△ Less
Submitted 29 June, 2020;
originally announced June 2020.
-
Carbon ionization at Gbar pressures: an ab initio perspective on astrophysical high-density plasmas
Authors:
Mandy Bethkenhagen,
Bastian B. L. Witte,
Maximilian Schörner,
Gerd Röpke,
Tilo Döppner,
Dominik Kraus,
Siegfried H. Glenzer,
Philip A. Sterne,
Ronald Redmer
Abstract:
A realistic description of partially-ionized matter in extreme thermodynamic states is critical to model the interior and evolution of the multiplicity of high-density astrophysical objects. Current predictions of its essential property, the ionization degree, rely widely on analytical approximations that have been challenged recently by a series of experiments. Here, we propose a novel ab initio…
▽ More
A realistic description of partially-ionized matter in extreme thermodynamic states is critical to model the interior and evolution of the multiplicity of high-density astrophysical objects. Current predictions of its essential property, the ionization degree, rely widely on analytical approximations that have been challenged recently by a series of experiments. Here, we propose a novel ab initio approach to calculate the ionization degree directly from the dynamic electrical conductivity using the Thomas-Reiche-Kuhn sum rule. This Density Functional Theory framework captures genuinely the condensed matter nature and quantum effects typical for strongly-correlated plasmas. We demonstrate this new capability for carbon and hydrocarbon, which most notably serve as ablator materials in inertial confinement fusion experiments aiming at recreating stellar conditions. We find a significantly higher carbon ionization degree than predicted by commonly used models, yet validating the qualitative behavior of the average atom model Purgatorio. Additionally, we find the carbon ionization state to remain unchanged in the environment of fully-ionized hydrogen. Our results will not only serve as benchmark for traditional models, but more importantly provide an experimentally accessible quantity in the form of the electrical conductivity.
△ Less
Submitted 28 April, 2020;
originally announced April 2020.
-
Material Properties for the Interiors of Massive Giant Planets and Brown Dwarfs
Authors:
Andreas Becker,
Mandy Bethkenhagen,
Clemens Kellermann,
Johannes Wicht,
Ronald Redmer
Abstract:
We present thermodynamic material and transport properties for the extreme conditions prevalent in the interiors of massive giant planets and brown dwarfs. They are obtained from extensive \textit{ab initio} simulations of hydrogen-helium mixtures along the isentropes of three representative objects. In particular, we determine the heat capacities, the thermal expansion coefficient, the isothermal…
▽ More
We present thermodynamic material and transport properties for the extreme conditions prevalent in the interiors of massive giant planets and brown dwarfs. They are obtained from extensive \textit{ab initio} simulations of hydrogen-helium mixtures along the isentropes of three representative objects. In particular, we determine the heat capacities, the thermal expansion coefficient, the isothermal compressibility, and the sound velocity. Important transport properties such as the electrical and thermal conductivity, opacity, and shear viscosity are also calculated. Further results for associated quantities including magnetic and thermal diffusivity, kinematic shear viscosity, as well as the static Love number $k_2$ and the equidistance are presented. In comparison to Jupiter-mass planets, the behavior inside massive giant planets and brown dwarfs is stronger dominated by degenerate matter. We discuss the implications on possible dynamics and magnetic fields of those massive objects. The consistent data set compiled here may serve as starting point to obtain material and transport properties for other substellar H-He objects with masses above one Jovian mass and finally may be used as input for dynamo simulations.
△ Less
Submitted 13 August, 2018;
originally announced August 2018.
-
Planetary Ices and the Linear Mixing Approximation
Authors:
Mandy Bethkenhagen,
Edmund R. Meyer,
Sebastien Hamel,
Nadine Nettelmann,
Martin French,
Ludwig Scheibe,
Christopher Ticknor,
Lee A. Collins,
Joel D. Kress,
Jonathan J. Fortney,
Ronald Redmer
Abstract:
The validity of the widely used linear mixing approximation for the equations of state (EOS) of planetary ices is investigated at pressure-temperature conditions typical for the interior of Uranus and Neptune. The basis of this study are ab initio data ranging up to 1000 GPa and 20 000 K calculated via density functional theory molecular dynamics simulations. In particular, we calculate a new EOS…
▽ More
The validity of the widely used linear mixing approximation for the equations of state (EOS) of planetary ices is investigated at pressure-temperature conditions typical for the interior of Uranus and Neptune. The basis of this study are ab initio data ranging up to 1000 GPa and 20 000 K calculated via density functional theory molecular dynamics simulations. In particular, we calculate a new EOS for methane and EOS data for the 1:1 binary mixtures of methane, ammonia, and water, as well as their 2:1:4 ternary mixture. Additionally, the self-diffusion coefficients in the ternary mixture are calculated along three different Uranus interior profiles and compared to the values of the pure compounds. We find that deviations of the linear mixing approximation from the results of the real mixture are generally small; for the thermal EOS they amount to 4% or less. The diffusion coefficients in the mixture agree with those of the pure compounds within 20% or better. Finally, a new adiabatic model of Uranus with an inner layer of almost pure ices is developed. The model is consistent with the gravity field data and results in a rather cold interior ($\mathrm{T_{core}} \mathtt{\sim}$ 4000 K).
△ Less
Submitted 13 September, 2017;
originally announced September 2017.
-
Uranus evolution models with simple thermal boundary layers
Authors:
N. Nettelmann,
K. Wang,
J. J. Fortney,
S. Hamel,
S. Yellamilli,
M. Bethkenhagen,
R. Redmer
Abstract:
The strikingly low luminosity of Uranus (Teff ~ Teq) constitutes a long-standing challenge to our understanding of Ice Giant planets. Here we present the first Uranus structure and evolution models that are constructed to agree with both the observed low luminosity and the gravity field data. Our models make use of modern ab initio equations of state at high pressures for the icy components water,…
▽ More
The strikingly low luminosity of Uranus (Teff ~ Teq) constitutes a long-standing challenge to our understanding of Ice Giant planets. Here we present the first Uranus structure and evolution models that are constructed to agree with both the observed low luminosity and the gravity field data. Our models make use of modern ab initio equations of state at high pressures for the icy components water, methane, and ammonia. Proceeding step by step, we confirm that adiabatic models yield cooling times that are too long, even when uncertainties in the ice:rock ratio (I:R) are taken into account. We then argue that the transition between the ice/rock-rich interior and the H/He-rich outer envelope should be stably stratified. Therefore, we introduce a simple thermal boundary and adjust it to reproduce the low luminosity. Due to this thermal boundary, the deep interior of the Uranus models are up to 2--3 warmer than adiabatic models, necessitating the presence of rocks in the deep interior with a possible I:R of $1\times$ solar. Finally, we allow for an equilibrium evolution (Teff ~ Teq) that begun prior to the present day, which would therefore no longer require the current era to be a "special time" in Uranus' evolution. In this scenario, the thermal boundary leads to more rapid cooling of the outer envelope. When Teff ~ Teq is reached, a shallow, subadiabatic zone in the atmosphere begins to develop. Its depth is adjusted to meet the luminosity constraint. This work provides a simple foundation for future Ice Giant structure and evolution models, that can be improved by properly treating the heat and particle fluxes in the diffusive zones.
△ Less
Submitted 30 April, 2016;
originally announced May 2016.