-
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.
-
Atom-in-jellium predictions of the shear modulus at high pressure
Authors:
Damian C. Swift,
Thomas Lockard,
Sebastien Hamel,
Christine J. Wu,
Lorin X. Benedict,
Philip A. Sterne
Abstract:
Atom-in-jellium calculations of the Einstein frequency in condensed matter and of the equation of state were used to predict the variation of shear modulus from zero pressure to ~$10^7$ g/cm$^3$, for several elements relevant to white dwarf (WD) stars and other self-gravitating systems. This is by far the widest range reported electronic structure calculation of shear modulus, spanning from ambien…
▽ More
Atom-in-jellium calculations of the Einstein frequency in condensed matter and of the equation of state were used to predict the variation of shear modulus from zero pressure to ~$10^7$ g/cm$^3$, for several elements relevant to white dwarf (WD) stars and other self-gravitating systems. This is by far the widest range reported electronic structure calculation of shear modulus, spanning from ambient through the one-component plasma to extreme relativistic conditions. The predictions were based on a relationship between Debye temperature and shear modulus which we assess to be accurate at the o(10%) level, and is the first known use of atom-in-jellium theory to calculate a shear modulus. We assessed the overall accuracy of the method by comparing with experimental measurements and more detailed electronic structure calculations at lower pressures.
△ Less
Submitted 27 August, 2021; v1 submitted 25 May, 2021;
originally announced May 2021.
-
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.
-
Atom-in-jellium equations of state and melt curves in the white dwarf regime
Authors:
Damian C. Swift,
Thomas Lockard,
Sebastien Hamel,
Christine J. Wu,
Lorin X. Benedict,
Philip A. Sterne,
Heather D. Whitley
Abstract:
Atom-in-jellium calculations of the electron states, and perturbative calculations of the Einstein frequency, were used to construct equations of state (EOS) from around $10^{-5}$ to $10^7$g/cm$^3$ and $10^{-4}$ to $10^{6}$eV for elements relevant to white dwarf (WD) stars. This is the widest range reported for self-consistent electronic shell structure calculations. Elements of the same ratio of…
▽ More
Atom-in-jellium calculations of the electron states, and perturbative calculations of the Einstein frequency, were used to construct equations of state (EOS) from around $10^{-5}$ to $10^7$g/cm$^3$ and $10^{-4}$ to $10^{6}$eV for elements relevant to white dwarf (WD) stars. This is the widest range reported for self-consistent electronic shell structure calculations. Elements of the same ratio of atomic weight to atomic number were predicted to asymptote to the same $T=0$ isotherm, suggesting that, contrary to recent studies of the crystallization of WDs, the amount of gravitational energy that could be released by separation of oxygen and carbon is small. A generalized Lindemann criterion based on the amplitude of the ion-thermal oscillations calculated using atom-in-jellium theory, previously used to extrapolate melt curves for metals, was found to reproduce previous thermodynamic studies of the melt curve of the one component plasma with a choice of vibration amplitude consistent with low pressure results. For elements for which low pressure melting satisfies the same amplitude criterion, such as Al, this melt model thus gives a likely estimate of the melt curve over the full range of normal electronic matter; for the other elements, it provides a useful constraint on the melt locus.
△ Less
Submitted 4 March, 2021;
originally announced March 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.
-
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.
-
Mass-radius relationships for exoplanets
Authors:
Damian Swift,
Jon Eggert,
Damien Hicks,
Sebastien Hamel,
Kyle Caspersen,
Eric Schwegler,
Gilbert Collins,
Nadine Nettelmann,
Graeme Ackland
Abstract:
For planets other than Earth, interpretation of the composition and structure depends largely on comparing the mass and radius with the composition expected given their distance from the parent star. The composition implies a mass-radius relation which relies heavily on equations of state calculated from electronic structure theory and measured experimentally on Earth. We lay out a method for deri…
▽ More
For planets other than Earth, interpretation of the composition and structure depends largely on comparing the mass and radius with the composition expected given their distance from the parent star. The composition implies a mass-radius relation which relies heavily on equations of state calculated from electronic structure theory and measured experimentally on Earth. We lay out a method for deriving and testing equations of state, and deduce mass-radius and mass-pressure relations for key materials whose equation of state is reasonably well established, and for differentiated Fe/rock. We find that variations in the equation of state, such as may arise when extrapolating from low pressure data, can have significant effects on predicted mass- radius relations, and on planetary pressure profiles. The relations are compared with the observed masses and radii of planets and exoplanets. Kepler-10b is apparently 'Earth- like,' likely with a proportionately larger core than Earth's, nominally 2/3 of the mass of the planet. CoRoT-7b is consistent with a rocky mantle over an Fe-based core which is likely to be proportionately smaller than Earth's. GJ 1214b lies between the mass-radius curves for H2O and CH4, suggesting an 'icy' composition with a relatively large core or a relatively large proportion of H2O. CoRoT-2b is less dense than the hydrogen relation, which could be explained by an anomalously high degree of heating or by higher than assumed atmospheric opacity. HAT-P-2b is slightly denser than the mass-radius relation for hydrogen, suggesting the presence of a significant amount of matter of higher atomic number. CoRoT-3b lies close to the hydrogen relation. The pressure at the center of Kepler-10b is 1.5+1.2-1.0 TPa. The central pressure in CoRoT-7b is probably close to 0.8TPa, though may be up to 2TPa.
△ Less
Submitted 15 September, 2011; v1 submitted 26 January, 2010;
originally announced January 2010.