We present a new equation of state (EOS) for dense hydrogen/helium mixtures that covers a range of densities from 10−8 to , pressures from 10−9 to 1013 GPa, and temperatures from 102 to 108 K. The ...calculations combine the EOS of Saumon, Chabrier & van Horn in the low-density, low-temperature molecular/atomic domain, the EOS of Chabrier & Potekhin in the high-density, high-temperature fully ionized domain, the limits of which differ for H and He, and ab initio quantum molecular dynamics calculations in the regime of intermediate density and temperature, characteristic of pressure dissociation and ionization. The EOS for the H/He mixture is based on the so-called additive volume law and thus does not take into account the interactions between the two species. A major improvement of the present calculations over existing ones is that we calculate the entropy over the entire density-temperature domain, a necessary quantity for calculations of stellar or planetary evolution. The EOS results are compared with existing experimental data, namely Hugoniot shock experiments for pure H and He, and with first-principles numerical simulations for both the single elements and the mixture. This new EOS covers a wide range of physical and astrophysical conditions, from Jovian planets to solar-type stars, and recovers the existing relativistic EOS at very high densities, in the domains of white dwarfs and neutron stars. All the tables are made publicly available.
While conventional interior models for Jupiter and Saturn are based on the simplistic assumption of a solid core surrounded by a homogeneous gaseous envelope, we have derived new models with an ...inhomogeneous distribution of heavy elements within these planets. Such a compositional gradient hampers large-scale convection that turns into double-diffusive convection, yielding an inner thermal profile that departs from the traditionally assumed adiabatic interior and affecting these planets heat content and cooling history. To address this problem, we have developed an analytical approach to describe layered double-diffusive convection and apply this formalism to solar system gaseous giant planet interiors. These models satisfy all observational constraints and yield values for the metal enrichment of our gaseous giants that are up to 30% to 60% higher than previously thought. The models also constrain the size of the convective layers within the planets. Because the heavy elements tend to be redistributed within the gaseous envelope, the models predict smaller than usual central cores inside Saturn and Jupiter, with possibly no core for the latter. These models open a new window and raise new challenges to our understanding of the internal structure of giant (solar and extrasolar) planets, in particular on how to determine their heavy material content, a key diagnostic for planet formation theories.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
Aims. We study the relative importance of several recent updates of microphysics input to the neutron star cooling theory and the effects brought about by superstrong magnetic fields of magnetars, ...including the effects of the Landau quantization in their crusts. Methods. We use a finite-difference code for simulation of neutron-star thermal evolution on timescales from hours to megayears with an updated microphysics input. The consideration of short timescales (≲1 yr) is made possible by a treatment of the heat-blanketing envelope without the quasistationary approximation inherent to its treatment in traditional neutron-star cooling codes. For the strongly magnetized neutron stars, we take into account the effects of Landau quantization on thermodynamic functions and thermal conductivities. We simulate cooling of ordinary neutron stars and magnetars with non-accreted and accreted crusts and compare the results with observations. Results. Suppression of radiative and conductive opacities in strongly quantizing magnetic fields and formation of a condensed radiating surface substantially enhance the photon luminosity at early ages, making the life of magnetars brighter but shorter. These effects together with the effect of strong proton superfluidity, which slows down the cooling of kiloyear-aged neutron stars, can explain thermal luminosities of about a half of magnetars without invoking heating mechanisms. Observed thermal luminosities of other magnetars are still higher than theoretical predictions, which implies heating, but the effects of quantizing magnetic fields and baryon superfluidity help to reduce the discrepancy.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
Context. The modeling of planetary interiors requires accurate equations of state (EOSs) for the basic constituents with proven validity in the difficult pressure–temperature regime extending up to ...50 000 K and hundreds of megabars. While EOSs based on first-principles simulations are now available for the two most abundant elements, hydrogen and helium, the situation is less satisfactory for water where no wide-range EOS is available despite its requirement for interior modeling of planets ranging from super-Earths to planets several times the size of Jupiter. Aims. As a first step toward a multi-phase EOS for dense water, we develop a temperature-dependent EOS for dense water covering the liquid and plasma regimes and extending to the super-ionic and gas regimes. This equation of state covers the complete range of conditions encountered in planetary modeling. Methods. We use first-principles quantum molecular dynamics simulations and the Thomas-Fermi extension to reach the highest pressures encountered in giant planets several times the size of Jupiter. Using these results, as well as the data available at lower pressures, we obtain a parametrization of the Helmholtz free energy adjusted over this extended temperature and pressure domain. The parametrization ignores the entropy and density jumps at phase boundaries but we show that it is sufficiently accurate to model interior properties of most planets and exoplanets. Results. We produce an EOS given in analytical form that is readily usable in planetary modeling codes and dynamical simulations (a fortran implementation is provided). The EOS produced is valid for the entire density range relevant to planetary modeling, for densities where quantum effects for the ions can be neglected, and for temperatures below 50 000K. We use this EOS to calculate the mass-radius relationship of exoplanets up to 5000 MEarth, explore temperature effects in the wet Earth-like, ocean planets and pure water planets, and quantify the influence of the water EOS for the core on the gravitational moments of Jupiter.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
Aims. We study the long-term thermal evolution of neutron stars in soft X-ray transients (SXTs), taking the deep crustal heating into account consistently with the changes of the composition of the ...crust. We collect observational estimates of average accretion rates and thermal luminosities of such neutron stars and compare the theory with observations. Methods. We performed simulations of thermal evolution of accreting neutron stars, considering the gradual replacement of the original nonaccreted crust by the reprocessed accreted matter, the neutrino and photon energy losses, and the deep crustal heating due to nuclear reactions in the accreted crust. We also tested and compared results for different modern theoretical models. We updated a compilation of the observational estimates of the thermal luminosities in quiescence and average accretion rates in the SXTs and compared the observational estimates with the theoretical results. Results. The long-term thermal evolution of transiently accreting neutron stars is nonmonotonic. The quasi-equilibrium temperature in quiescence reaches a minimum and then increases toward the final steady state. The quasi-equilibrium thermal luminosity of a neutron star in an SXT can be substantially lower at the minimum than in the final state. This enlarges the range of possibilities for theoretical interpretation of observations of such neutron stars. The updates of the theory and observations leave the previous conclusions unchanged, namely that the direct Urca process operates in relatively cold neutron stars and that an accreted heat-blanketing envelope is likely present in relatively hot neutron stars in the SXTs in quiescence. The results of the comparison of theory with observations favor suppression of the triplet pairing type of nucleon superfluidity in the neutron-star matter.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
We present an implementation of the Hall term in the non-ideal magnetohydrodynamics equations into the adaptive-mesh-refinement code RAMSES to study its impact on star formation. Recent works show ...that the Hall effect heavily influences the regulation of the angular momentum in collapsing dense cores, strengthening or weakening the magnetic braking. Our method consists of a modification of the two-dimensional constrained transport scheme. Our scheme shows convergence of second order in space and the frequency of the propagation of whistler waves is accurate. We confirm previous results, namely that during the collapse, the Hall effect generates a rotation of the fluid with a direction in the mid-plane that depends on the sign of the Hall resistivity, while counter-rotating envelopes develop on each side of the mid-plane. However, we find that the predictability of our numerical results is severely limited. The angular momentum is not conserved in any of our dense core-collapse simulations with the Hall effect: a large amount of angular momentum is generated within the first Larson core, a few hundred years after its formation, without compensation by the surrounding gas. This issue is not mentioned in previous studies and may be correlated to the formation of the accretion shock on the Larson core. We expect that this numerical effect could be a serious issue in star formation simulations.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
Following up our recent analysis devoted to the impact of non steady accretion on the location of young low-mass stars or brown dwarfs in the Herzsprung-Russell diagram, we perform a detailed ...analysis devoted to the effect of burst accretion on the internal structure of low-mass and solar type stars. We find that episodic accretion can produce objects with significantly higher central temperatures than the ones of the non accreting counterparts of same mass and age. As a consequence, lithium depletion can be severely enhanced in these objects. This provides a natural explanation for the unexpected level of lithium depletion observed in young objects for the inferred age of their parent cluster. These results confirm the limited reliability of lithium abundance as a criterion for assessing or rejecting cluster membership. They also show that lithium is not a reliable age indicator, because its fate strongly depends on the past accretion history of the star. Under the assumption that giant planets primarily form in massive disks prone to gravitational instability and thus to accretion burst episodes, the same analysis also explains the higher Li depletion observed in planet hosting stars. At last, we show that, depending on the burst rate and intensity, accretion outbursts can produce solar mass stars with lower convective envelope masses, at ages less than a few tens of Myr, than predicted by standard (non or slowly accreting) pre-main sequence models. This result has interesting, although speculative, implications for the recently discovered depletion of refractory elements in the Sun.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
We develop a detailed chemical network relevant to calculate the conditions that are characteristic of prestellar core collapse. We solve the system of time-dependent differential equations to ...calculate the equilibrium abundances of molecules and dust grains, with a size distribution given by size-bins for these latter. These abundances are used to compute the different non-ideal magneto-hydrodynamics resistivities (ambipolar, Ohmic and Hall), needed to carry out simulations of protostellar collapse. For the first time in this context, we take into account the evaporation of the grains, the thermal ionisation of potassium, sodium, and hydrogen at high temperature, and the thermionic emission of grains in the chemical network, and we explore the impact of various cosmic ray ionisation rates. All these processes significantly affect the non-ideal magneto-hydrodynamics resistivities, which will modify the dynamics of the collapse. Ambipolar diffusion and Hall effect dominate at low densities, up to nH = 1012 cm-3, after which Ohmic diffusion takes over. We find that the time-scale needed to reach chemical equilibrium is always shorter than the typical dynamical (free fall) one. This allows us to build a large, multi-dimensional multi-species equilibrium abundance table over a large temperature, density and ionisation rate ranges. This table, which we make accessible to the community, is used during first and second prestellar core collapse calculations to compute the non-ideal magneto-hydrodynamics resistivities, yielding a consistent dynamical-chemical description of this process.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
The anomalously large radii of strongly irradiated exoplanets have remained a major puzzle in astronomy. Based on a two-dimensional steady-state atmospheric circulation model, the validity of which ...is assessed by comparison to three-dimensional calculations, we reveal a new mechanism, namely the advection of the potential temperature due to mass and longitudinal momentum conservation, a process occurring in the Earth's atmosphere or oceans. In the deep atmosphere, the vanishing heating flux forces the atmospheric structure to converge to a hotter adiabat than the one obtained with 1D calculations, implying a larger radius for the planet. Not only do the calculations reproduce the observed radius of HD 209458b, but also reproduce the observed correlation between radius inflation and irradiation for transiting planets. Vertical advection of potential temperature induced by non-uniform atmospheric heating thus provides a robust mechanism to explain the inflated radii of irradiated hot Jupiters.
Context.We examine the evolution of low-mass star and brown dwarf eclipsing binaries. These objects are rapid rotators and are believed to shelter large magnetic fields. Aims.We suggest that reduced ...convective efficiency, due to fast rotation and large field strengths, and/or to magnetic spot coverage of the radiating surface significantly affect their evolution, leading to a reduced heat flux and thus larger radii and cooler effective temperatures than for regular objects. Methods.We have considered such processes in our evolutionary calculations, using a phenomenological approach. Results.This yields mass–radius and effective temperature–radius relationships in agreement with the observations. We also reproduce the effective temperature ratio and the radii of the two components of the recently discovered puzzling eclipsing brown dwarf system. Conclusions.These calculations show that fast rotation and/or magnetic activity may significantly affect the evolution of eclipsing binaries and that the mechanical and thermal properties of these objects depart from the ones of non-active low-mass objects. We find that, for internal field strengths compatible with the observed surface value of a few kiloGauss, convection can be severely inhibited. The onset of a central radiative zone for rapidly rotating active low-mass stars might thus occur below the usual ~$0.35~M_\odot$ limit.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK