Thermal evolution and structure models of the transiting super-Earth GJ 1214b

The planet GJ 1214b is the second known super-Earth with a measured mass and radius. Orbiting a quiet M-star, it receives considerably less mass-loss driving X-ray and UV radiation than CoRoT-7b, so that the interior may be quite dissimilar in composition, including the possibility of a large fraction of water. We model the interior of GJ 1214b assuming a two-layer (envelope+rock core) structure where the envelope material is either H/He, pure water, or a mixture of H/He and H2O. Within this framework we perform models of the thermal evolution and contraction of the planet. We discuss possible compositions that are consistent with Mp=6.55 ME, Rp=2.678 RE, an age tau=3-10 Gyr, and the irradiation level of the atmosphere. These conditions require that if water exists in the interior, it must remain in a fluid state, with important consequences for magnetic field generation. These conditions also require the atmosphere to have a deep isothermal region extending down to 80-800 bar, depending on composition. Our results bolster the suggestion of a metal-enriched H/He atmosphere for the planet, as we find water-world models that lack an H/He atmosphere to require an implausibly large water-to-rock ratio of more than 6:1. We instead favor a H/He/H2O envelope with high water mass fraction (~0.5-0.85), similar to recent models of the deep envelope of Uranus and Neptune. Even with these high water mass fractions in the H/He envelope, generally the bulk composition of the planet can have subsolar water:rock ratios. Dry, water-enriched, and pure water envelope models differ to an observationally significant level in their tidal Love numbers k2 of respectively ~0.018, 0.15, and 0.7.Comment: 11 pages, 6 figures, 1 table, accepted to Ap

New indication for a dichotomy in the interior structure of Uranus and Neptune from the application of modified shape and rotation data

Since the Voyager fly-bys of Uranus and Neptune, improved gravity field data have been derived from long-term observations of the planets' satellite motions, and modified shape and solid-body rotation periods were suggested. A faster rotation period (-40 min) for Uranus and a slower rotation period (+1h20) of Neptune compared to the Voyager data were found to minimize the dynamical heights and wind speeds. We apply the improved gravity data, the modified shape and rotation data, and the physical LM-R equation of state to compute adiabatic three-layer structure models, where rocks are confined to the core, and homogeneous thermal evolution models of Uranus and Neptune. We present the full range of structure models for both the Voyager and the modified shape and rotation data. In contrast to previous studies based solely on the Voyager data or on empirical EOS, we find that Uranus and Neptune may differ to an observationally significant level in their atmospheric heavy element mass fraction Z1 and nondimensional moment of inertia, nI. For Uranus, we find Z1 < 8% and nI=0.2224(1), while for Neptune Z1 < 65% and nI=0.2555(2) when applying the modified shape and rotation data, while for the unmodified data we compute Z1 < 17% and nI=0.230(1) for Uranus and Z1 < 54% and nI=0.2410(8) for Neptune. In each of these cases, solar metallicity models (Z1=0.015) are still possible. The cooling times obtained for each planet are similar to recent calculations with the Voyager rotation periods: Neptune's luminosity can be explained by assuming an adiabatic interior while Uranus cools far too slowly. More accurate determinations of these planets' gravity fields, shapes, rotation periods, atmospheric heavy element abundances, and intrinsic luminosities are essential for improving our understanding of the internal structure and evolution of icy planets.Comment: accepted to Planet. Space Sci., special editio

Uranus evolution models with simple thermal boundary layers

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.Comment: 13 pages, Accepted to Icaru

Correlations in Hot Dense Helium

Hot dense helium is studied with first-principles computer simulations. By combining path integral Monte Carlo and density functional molecular dynamics, a large temperature and density interval ranging from 1000 to 1000000 K and 0.4 to 5.4 g/cc becomes accessible to first-principles simulations and the changes in the structure of dense hot fluids can be investigated. The focus of this article are pair correlation functions between nuclei, between electrons, and between electrons and nuclei. The density and temperature dependence of these correlation functions is analyzed in order to describe the structure of the dense fluid helium at extreme conditions.Comment: accepted for publication in Journal of Physics

Probing the interiors of the ice giants: Shock compression of water to 700 GPa and 3.8 g/ccm

Recently there has been tremendous increase in the number of identified extra-solar planetary systems. Our understanding of their formation is tied to exoplanet internal structure models, which rely upon equations of state of light elements and compounds like water. Here we present shock compression data for water with unprecedented accuracy that shows water equations of state commonly used in planetary modeling significantly overestimate the compressibility at conditions relevant to planetary interiors. Furthermore, we show its behavior at these conditions, including reflectivity and isentropic response, is well described by a recent first-principles based equation of state. These findings advocate this water model be used as the standard for modeling Neptune, Uranus, and "hot Neptune" exoplanets, and should improve our understanding of these types of planets.Comment: Accepted to Phys. Rev. Lett.; supplementary material attached including 2 figures and 2 tables; to view attachments, please download and extract the gzipped tar source file listed under "Other formats

Interior structure models of GJ 436b

GJ 436b is the first extrasolar planet discovered that resembles Neptune in mass and radius. The particularly interesting property of Neptune-sized planets is that their mass Mp and radius Rp are close to theoretical M-R relations of water planets. Given Mp, Rp, and equilibrium temperature, however, various internal compositions are possible. A broad set of interior structure models is presented here that illustrates the dependence of internal composition and possible phases of water occurring in presumably water-rich planets, such as GJ 436b on the uncertainty in atmospheric temperature profile and mean density. We show how the set of solutions can be narrowed down if theoretical constraints from formation and model atmospheres are applied or potentially observational constraints for the atmospheric metallicity Z1 and the tidal Love number k2. We model the interior by assuming either three layers (hydrogen-helium envelope, water layer, rock core) or two layers (H/He/H2O envelope, rocky core). For water, we use the equation of state H2O-REOS based on FT-DFT-MD simulations. Some admixture of H/He appears mandatory for explaining the measured radius. For the warmest considered models, the H/He mass fraction can reduce to 10^-3, still extending over ~0.7 REarth. If water occurs, it will be essentially in the plasma phase or in the superionic phase, but not in an ice phase. Metal-free envelope models have 0.02<k2<0.2, and the core mass cannot be determined from a measurement of k2. In contrast, models with 0.3<k2<0.82 require high metallicities Z1<0.89 in the outer envelope. The uncertainty in core mass decreases to 0.4 Mp, if k2>0.3, and further to 0.2 Mp, if k2>0.5, and core mass and Z1 become sensitive functions of k2. To further narrow the set of solutions, a proper treatment of the atmosphere and the evolution is necessary.Comment: 9 pages, accepted to A&

Analysis of new high-precision transit light curves of WASP-10 b: starspot occultations, small planetary radius, and high metallicity

The WASP-10 planetary system is intriguing because different values of radius have been reported for its transiting exoplanet. The host star exhibits activity in terms of photometric variability, which is caused by the rotational modulation of the spots. Moreover, a periodic modulation has been discovered in transit timing of WASP-10 b, which could be a sign of an additional body perturbing the orbital motion of the transiting planet. We attempt to refine the physical parameters of the system, in particular the planetary radius, which is crucial for studying the internal structure of the transiting planet. We also determine new mid-transit times to confirm or refute observed anomalies in transit timing. We acquired high-precision light curves for four transits of WASP-10 b in 2010. Assuming various limb-darkening laws, we generated best-fit models and redetermined parameters of the system. The prayer-bead method and Monte Carlo simulations were used to derive error estimates. Three transit light curves exhibit signatures of the occultations of dark spots by the planet during its passage across the stellar disk. The influence of stellar activity on transit depth is taken into account while determining system parameters. The radius of WASP-10 b is found to be no greater than 1.03 Jupiter radii, a value significantly smaller than most previous studies indicate. We calculate interior structure models of the planet, assuming a two-layer structure with one homogeneous envelope atop a rock core. The high value of the WASP-10 b's mean density allows one to consider the planet's internal structure including 270 to 450 Earth masses of heavy elements. Our new mid-transit times confirm that transit timing cannot be explained by a constant period if all literature data points are considered. They are consistent with the ephemeris assuming a periodic variation of transit timing...Comment: Accepted for publication in A&

Ultrahigh compression of water using intense heavy ion beams: laboratory planetary physics

Intense heavy ion beams offer a unique tool for generating samples of high energy density matter with extreme conditions of density and pressure that are believed to exist in the interiors of giant planets. An international accelerator facility named FAIR (Facility for Antiprotons and Ion Research) is being constructed at Darmstadt, which will be completed around the year 2015. It is expected that this accelerator facility will deliver a bunched uranium beam with an intensity of 5x10(11) ions per spill with a bunch length of 50-100 ns. An experiment named LAPLAS (Laboratory Planetary Sciences) has been proposed to achieve a low-entropy compression of a sample material like hydrogen or water (which are believed to be abundant in giant planets) that is imploded in a multi-layered target by the ion beam. Detailed numerical simulations have shown that using parameters of the heavy ion beam that will be available at FAIR, one can generate physical conditions that have been predicted to exist in the interior of giant planets. In the present paper, we report simulations of compression of water that show that one can generate a plasma phase as well as a superionic phase of water in the LAPLAS experiments