Materials at high pressures and temperatures are of great interest for planetary science and astrophysics, warm dense matter physics, and inertial confinement fusion research. Planetary structure models rely on our understanding of the behaviour of elements (and their mixtures) at exotic conditions that do not exist on Earth, and at the same time planets serve as natural laboratories for studying materials at extreme conditions. In this review we discuss the connection between modelling planetary interiors and high-pressure physics of hydrogen and helium. First, we summarise key experiments for determining the equation of state and phase diagram of hydrogen and helium as well as state-of-the-art theoretical approaches. We next briefly review our current knowledge of the internal structures of the giant planets in the Solar System, Jupiter and Saturn, and the importance of high pressure physics to their characterisation.

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.

Atmospheres regulate the planetary heat loss and therefore influence planetary thermal evolution. Uncertainty in a giant planet's thermal state contributes to the uncertainty in the inferred abundance of heavy elements it contains. Within an analytic atmosphere model, we here investigate the influence that different cloud opacities and cloud depths can have on the metallicity of irradiated extrasolar gas giants, which is inferred from interior models. In this work, the link between inferred metallicity and assumed cloud properties is the thermal profile of atmosphere and interior. Therefore, we perform coupled atmosphere, interior, and evolution calculations. The atmosphere model includes clouds in a much simplified manner; it includes long-wave absorption but neglects shortwave scattering. Within that model, we show that optically thick, high clouds have negligible influence, whereas deep-seated, optically very thick clouds can lead to warmer deep tropospheres and therefore higher bulk heavy element mass estimates. For the young hot Jupiter WASP-10b, we find a possible enhancement in inferred metallicity of up to 10% due to possible silicate clouds at $\sim$0.3 bar. For WASP-39b, whose observationally derived metallicity is higher than predicted by cloudless models, we find an enhancement by at most 50%. However, further work on cloud properties and their self-consistent coupling to the atmospheric structure is needed in order to reduce uncertainties in the choice of model parameter values, in particular of cloud opacities.

The brightness of Neptune is often found to be in accordance with an adiabatic interior, while the low luminosity of Uranus challenges this assumption. Here we apply revised equation of state data of hydrogen, helium, and water and compute the thermal evolution of Uranus and Neptune assuming an adiabatic interior. For this purpose, we have developed a new planetary model and evolution code. We investigate the influence of albedo, solar energy influx, and equations of state of H and He, and water on the cooling time. Our cooling times of about $\tau_\text{U}=5.1\times 10^9\text{ years}$ for Uranus and $\tau_\text{N}=3.7\times 10^9\text{ years}$ for Neptune bracket the known age of the planets of $4.56\times 10^9\text{ years}$ implying that neither planet's present-day luminosity can be explained by adiabatic cooling. We also find that uncertainties on input parameters such as the level of irradiation matter generally more for Uranus than for Neptune. Our results suggest that in contrast to common assumptions, neither planet is fully adiabatic in the deeper interior.

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.

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).

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.

We present new equations of state (EOS) for hydrogen and helium covering a wide range of temperatures from 60 K to 10$^7$ K and densities from $10^{-10}$ g/cm$^3$ to $10^3$ g/cm$^3$. They include an extended set of ab initio EOS data for the strongly correlated quantum regime with an accurate connection to data derived from other approaches for the neighboring regions. We compare linear-mixing isotherms based on our EOS tables with available real-mixture data. A first important astrophysical application of this new EOS data is the calculation of interior models for Jupiter and the comparison with recent results. Secondly, mass-radius relations are calculated for Brown Dwarfs which we compare with predictions derived from the widely used EOS of Saumon, Chabrier and van Horn. Furthermore, we calculate interior models for typical Brown Dwarfs with different masses, namely Corot-3b, Gliese-229b and Corot-15b, and the Giant Planet KOI-889b. The predictions for the central pressures and densities differ by up to 10$\%$ dependent on the EOS used. Our EOS tables are made available in the supplemental material of this paper.

PLATO 2.0 has recently been selected for ESA's M3 launch opportunity (2022/24). Providing accurate key planet parameters (radius, mass, density and age) in statistical numbers, it addresses fundamental questions such as: How do planetary systems form and evolve? Are there other systems with planets like ours, including potentially habitable planets? The PLATO 2.0 instrument consists of 34 small aperture telescopes (32 with 25 sec readout cadence and 2 with 2.5 sec candence) providing a wide field-of-view (2232 deg2) and a large photometric magnitude range (4-16 mag). It focusses on bright (4-11 mag) stars in wide fields to detect and characterize planets down to Earth-size by photometric transits, whose masses can then be determined by ground-based radial-velocity follow-up measurements. Asteroseismology will be performed for these bright stars to obtain highly accurate stellar parameters, including masses and ages. The combination of bright targets and asteroseismology results in high accuracy for the bulk planet parameters: 2%, 4-10% and 10% for planet radii, masses and ages, respectively. The planned baseline observing strategy includes two long pointings (2-3 years) to detect and bulk characterize planets reaching into the habitable zone (HZ) of solar-like stars and an additional step-and-stare phase to cover in total about 50% of the sky. PLATO 2.0 will observe up to 1,000,000 stars and detect and characterize hundreds of small planets, and thousands of planets in the Neptune to gas giant regime out to the HZ. It will therefore provide the first large-scale catalogue of bulk characterized planets with accurate radii, masses, mean densities and ages. This catalogue will include terrestrial planets at intermediate orbital distances, where surface temperatures are moderate. Coverage of this parameter range with statistical numbers of bulk characterized planets is unique to PLATO 2.0.

The core mass of Saturn is commonly assumed to be 10-25 ME as predicted by interior models with various equations of state (EOSs) and the Voyager gravity data, and hence larger than that of Jupiter (0-10 ME). We here re-analyze Saturn's internal structure and evolution by using more recent gravity data from the Cassini mission and different physical equations of state: the ab initio LM-REOS which is rather soft in Saturn's outer regions but stiff at high pressures, the standard Sesame-EOS which shows the opposite behavior, and the commonly used SCvH-i EOS. For all three EOS we find similar core mass ranges, i.e. of 0-20 ME for SCvH-i and Sesame EOS and of 0-17 ME for LM-REOS. Assuming an atmospheric helium mass abundance of 18%, we find maximum atmospheric metallicities, Zatm of 7x solar for SCvH-i and Sesame-based models and a total mass of heavy elements, MZ of 25-30 ME. Some models are Jupiter-like. With LM-REOS, we find MZ=16-20 ME, less than for Jupiter, and Zatm less than 3x solar. For Saturn, we compute moment of inertia values lambda=0.2355(5). Furthermore, we confirm that homogeneous evolution leads to cooling times of only about 2.5 Gyr, independent on the applied EOS. Our results demonstrate the need for accurately measured atmospheric helium and oxygen abundances, and of the moment of inertia for a better understanding of Saturn's structure and evolution.

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.

The amount and distribution of heavy elements in Jupiter gives indications on the process of its formation and evolution. Core mass and metallicity predictions however depend on the equations of state used, and on model assumptions. We present an improved ab initio hydrogen equation of state, H-REOS.2 and compute the internal structure and thermal evolution of Jupiter within the standard three-layer approach. The advance over our previous Jupiter models with H-REOS.1 by Nettelmann et al.(2008) is that the new models are also consistent with the observed 2 or more times solar heavy element abundances in Jupiter's atmosphere. Such models have a rock core mass Mcore=0-8 ME, total mass of heavy elements MZ=28-32 ME, a deep internal layer boundary at 4 or more Mbar, and a cooling time of 4.4-5.0 Gyrs when assuming homogeneous evolution. We also calculate two-layer models in the manner of Militzer et al.(2008) and find a comparable large core of 16-21 ME, out of which ~11 ME is helium, but a significantly higher envelope metallicity of 4.5 times solar. According to our preferred three-layer models, neither the characteristic frequency (nu0 ~156 microHz) nor the normalized moment of inertia (~0.276) are sensitive to the core mass but accurate measurements could well help to rule out some classes of models.

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.

Transit and radial velocity observations continuously discover an increasing number of exoplanets. However, when it comes to the composition of the observed planets the data are compatible with several interior structure models. Thus, a planetary parameter sensitive to the planet's density distribution could help constrain this large number of possible models even further. We aim to investigate to what extent an exoplanet's interior can be constrained in terms of core mass and envelope metallicity by taking the tidal Love number k_2 into account as an additional possibly observable parameter. Because it is the only planet with an observationally determined k_2, we constructed interior models for the Hot Jupiter exoplanet HAT-P-13b by solving the equations of hydrostatic equilibrium and mass conservation for different boundary conditions. In particular, we varied the surface temperature and the outer temperature profile, as well as the envelope metallicity within the widest possible parameter range. We also considered atmospheric conditions that are consistent with nongray atmosphere models. For all these models we calculated the Love number k_2 and compared it to the allowed range of k_2 values that could be obtained from eccentricity measurements of HAT-P-13b. We use the example of HAT-P-13b to show the general relationships between the quantities temperature, envelope metallicity, core mass, and Love number of a planet. For any given k_2 value a maximum possible core mass can be determined. For HAT-P-13b we find Mcore < 27 ME, based on the latest eccentricity measurement. We are able to constrain both the envelope and bulk metallicity of HAT-P-13b to 1 -- 11 times stellar metallicity and the extension of the isothermal layer in the planet's atmosphere to 3 -- 44 bar. Assuming equilibrium tidal theory, we find lower limits on the tidal Q consistent with 10^3 - 10^5.

We present the Young Exoplanet Transit Initiative (YETI), in which we use several 0.2 to 2.6m telescopes around the world to monitor continuously young (< 100 Myr), nearby (< 1 kpc) stellar clusters mainly to detect young transiting planets (and to study other variability phenomena on time-scales from minutes to years). The telescope network enables us to observe the targets continuously for several days in order not to miss any transit. The runs are typically one to two weeks long, about three runs per year per cluster in two or three subsequent years for about ten clusters. There are thousands of stars detectable in each field with several hundred known cluster members, e.g. in the first cluster observed, Tr-37, a typical cluster for the YETI survey, there are at least 469 known young stars detected in YETI data down to R=16.5 mag with sufficient precision of 50 milli-mag rms (5 mmag rms down to R=14.5 mag) to detect transits, so that we can expect at least about one young transiting object in this cluster. If we observe 10 similar clusters, we can expect to detect approximately 10 young transiting planets with radius determinations. The precision given above is for a typical telescope of the YETI network, namely the 60/90-cm Jena telescope (similar brightness limit, namely within +/-1 mag, for the others) so that planetary transits can be detected. For planets with mass and radius determinations, we can calculate the mean density and probe the internal structure. We aim to constrain planet formation models and their time-scales by discovering planets younger than 100 Myr and determining not only their orbital parameters, but also measuring their true masses and radii, which is possible so far only by the transit method. Here, we present an overview and first results. (Abstract shortened)

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.

In order to accurately model giant planets, a whole set of observational constraints is needed. As the conventional constraints for extrasolar planets like mass, radius, and temperature allow for a large number of acceptable models, a new planetary parameter is desirable in order to further constrain planetary models. Such a parameter may be the tidal Love number k2. In this paper we aim to study the capability of k2 to reveal further information about the interior structure of a planet. With theoretical planetary models we investigate how the tidal Love number k2 responds to the internal density distribution of a planet. In particular, we demonstrate the effect of the degeneracy of k2 due to a density discontinuity in the envelope of a three-layer planetary model. The effect of a possible outer density discontinuity masks the effect of the core mass on the Love number k2. Hence, there is no unique relationship between the Love number k2 and the core mass of a planet. We show that the degeneracy of k2 with respect to a layer boundary in the envelope also occurs in existing planets, e.g. Saturn and the Hot Neptune GJ436b. As a result of the degeneracy, the planetary parameter k2 cannot be used to further constrain Saturnian models and for GJ436b only a maximum possible core mass can be derived from a given k2. To significantly narrow the uncertainty about the core mass of GJ436b the combined knowledge of k2 and atmospheric metallicity and temperature profile is necessary.

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.

We investigate the effects of Pauli blocking on the properties of hydrogen at high pressures, where recent experiments have shown a transition from insulating behavior to metal-like conductivity. Since the Pauli principle prevents multiple occupation of electron states (Pauli blocking), atomic states disintegrate subsequently at high densities (Mott effect). We calculate the energy shifts due to Pauli blocking and discuss the Mott effect solving an effective Schroedinger equation for strongly correlated systems. The ionization equilibrium is treated on the basis of a chemical approach. Results for the ionization equilibrium and the pressure in the region 4.000 K < T < 20.000 K are presented. We show that the transition to a highly conducting state is softer than found in earlier work. A first order phase transition is observed at T < 6.450 K, but a diffuse transition appears still up to 20.000 K.

The equation of state of hydrogen, helium, and water effects interior structure models of giant planets significantly. We present a new equation of state data table, LM-REOS, generated by large scale quantum molecular dynamics simulations for hydrogen, helium, and water in the warm dense matter regime, i.e.for megabar pressures and temperatures of several thousand Kelvin, and by advanced chemical methods in the complementary regions. The influence of LM-REOS on the structure of Jupiter is investigated and compared with state-of-the-art results within a standard three-layer model consistent with astrophysical observations of Jupiter. Our new Jupiter models predict an important impact of mixing effects of helium in hydrogen with respect to an altered compressibility and immiscibility.