Normalized to: Hamel, S.
[1]
oai:arXiv.org:2006.16208 [pdf] - 2123772
Sensitivity of stellar physics to the equation of state
Submitted: 2020-06-29
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.
[2]
oai:arXiv.org:1709.04133 [pdf] - 1727642
Planetary Ices and the Linear Mixing Approximation
Bethkenhagen, Mandy;
Meyer, Edmund R.;
Hamel, Sebastien;
Nettelmann, Nadine;
French, Martin;
Scheibe, Ludwig;
Ticknor, Christopher;
Collins, Lee A.;
Kress, Joel D.;
Fortney, Jonathan J.;
Redmer, Ronald
Submitted: 2017-09-13
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).
[3]
oai:arXiv.org:1605.00171 [pdf] - 1399848
Uranus evolution models with simple thermal boundary layers
Submitted: 2016-04-30
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.
[4]
oai:arXiv.org:1001.4851 [pdf] - 451639
Mass-radius relationships for exoplanets
Submitted: 2010-01-26, last modified: 2011-09-15
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.