There are still many open questions regarding the nature of Uranus and Neptune, the outermost planets in the Solar System. In this review we summarize the current-knowledge about Uranus and Neptune with a focus on their composition and internal structure, formation including potential subsequent giant impacts, and thermal evolution. We present key open questions and discuss the uncertainty in the internal structures of the planets due to the possibility of non-adiabatic and inhomogeneous interiors. We also provide the reasoning for improved observational constraints on their fundamental physical parameters such as their gravitational and magnetic fields, rotation rates, and deep atmospheric composition and temperature. Only this way will we be able to improve our understating of these planetary objects, and the many similar-sized objects orbiting other stars.

We explore the application of machine learning based on mixture density neural networks (MDNs) to the interior characterization of low-mass exoplanets up to 25 Earth masses constrained by mass, radius, and fluid Love number $k_2$. We create a dataset of 900$\:$000 synthetic planets, consisting of an iron-rich core, a silicate mantle, a high-pressure ice shell, and a gaseous H/He envelope, to train a MDN using planetary mass and radius as inputs to the network. For this layered structure, we show that the MDN is able to infer the distribution of possible thicknesses of each planetary layer from mass and radius of the planet. This approach obviates the time-consuming task of calculating such distributions with a dedicated set of forward models for each individual planet. While gas-rich planets may be characterized by compositional gradients rather than distinct layers, the method presented here can be easily extended to any interior structure model. The fluid Love number $k_2$ bears constraints on the mass distribution in the planets' interior and will be measured for an increasing number of exoplanets in the future. Adding $k_2$as an input to the MDN significantly decreases the degeneracy of the possible interior structures.

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.

The Juno Orbiter is measuring the three-dimensional gravity field perturbation of Jupiter induced by its rapid rotation, zonal flows, and tidal response to its major natural satellites. This paper aims to provide the contributions to the tesseral harmonics coefficients Cnm, Snm, and the Love numbers knm to be expected from static tidal response in the gravity field of rotating Jupiter. For that purpose, we apply the method of Concentric Maclaurin Ellipsoids (CMS). As we are interested in the variation of the tidal potential with the longitudes of the moons, we take into account the simultaneous presence of the satellites Io, Europa, and Ganymede. We assume co-planar, circular orbits with normals parallel to Jupiter's spin axis. The planet-centered longitude of Io in the three-moon case is arbitrarily assumed varphi = 0. Under these assumptions we find maximum amplitudes and fluctuations of 3.5 times 10^-8 +- 15% for C22. For the Love numbers, largest variation of 10% to 20% is seen in k42 and k62, whereas the values k2, k33, and k44 fall into narrow ranges of 0.1% uncertainty or less. In particular, we find k2=k2,Io(1 +- 0.02%) where k2,Io=0.5897 is the static tidal response to lone Io. Our obtained gravity field perturbation leads to a maximum equatorial shape deformation of up to 28m. We suggest that should Juno measurements of the knm deviate from those values, it may be due to dynamic or dissipative effects on Jupiter's tidal response. Finally, an analytic expression is provided to calculate the tesseral harmonics contribution from static tidal response for any configuration of the satellites.

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 Juno Orbiter has provided improved estimates of the even gravitational harmonics J2 to J8 of Jupiter. To compute higher-order moments, new methods such as the Concentric Maclaurin Spheroids (CMS) method have been developed which surpass the so far commonly used Theory of Figures (ToF) method in accuracy. This progress rises the question whether ToF can still provide a useful service for deriving the internal structure of giant planets in the Solar system. In this paper, I apply both the ToF and the CMS method to compare results for polytropic Jupiter and for the physical equation of state H/He-REOS.3 based models. An accuracy in the computed values of J2 and J4 of 0.1% is found to be sufficient in order to obtain the core mass safely within 0.5 Mearth numerical accuracy and the atmospheric metallicity within about 0.0004. ToF to 4th order provides that accuracy, while ToF to 3rd order does not for J4. Furthermore, I find that the assumption of rigid rotation yields J6 and J8 values in agreement with the current Juno estimates, and that higher order terms (J10 to J18) deviate by about 10% from predictions by polytropic models. This work suggests that ToF4 can still be applied to infer the deep internal structure, and that the zonal winds on Jupiter reach less deep than 0.9 RJup.

We review our current understanding of the interior structure and thermal evolution of Saturn, with a focus on recent results in the Cassini era. There has been important progress in understanding physical inputs, including equations of state of planetary materials and their mixtures, physical parameters like the gravity field and rotation rate, and constraints on Saturnian free oscillations. At the same time, new methods of calculation, including work on the gravity field of rotating fluid bodies, and the role of interior composition gradients, should help to better constrain the state of Saturn's interior, now and earlier in its history. However, a better appreciation of modeling uncertainties and degeneracies, along with a greater exploration of modeling phase space, still leave great uncertainties in our understanding of Saturn's interior. Further analysis of Cassini data sets, as well as precise gravity field measurements from the Cassini Grand Finale orbits, will further revolutionize our understanding of Saturn's interior over the next few years.

We investigate the physical characteristics of the Solar System's proposed Planet Nine using modeling tools with a heritage in studying Uranus and Neptune. For a range of plausible masses and interior structures, we find upper limits on the intrinsic Teff, from ~35-50 K for masses of 5-20 M_Earth, and we also explore lower Teff values. Possible planetary radii could readily span from 3 to 6 R_Earth depending on the mass fraction of any H/He envelope. Given its cold temperature, the planet encounters significant methane condensation, which dramatically alters the atmosphere away from simple Neptune-like expectations. We find the atmosphere is strongly depleted in molecular absorption at visible wavelengths, suggesting a Rayleigh scattering atmosphere with a high geometric albedo approaching 0.75. We highlight two diagnostics for the atmosphere's temperature structure, the first being the value of the methane mixing ratio above the methane cloud. The second is the wavelength at which cloud scattering can be seen, which yields the cloud-top pressure. Surface reflection may be seen if the atmosphere is thin. Due to collision-induced opacity of H2 in the infrared, the planet would be extremely blue (instead of red) in the shortest wavelength WISE colors if methane is depleted, and would, in some cases, exist on the verge of detectability by WISE. For a range of models, thermal fluxes from ~3-5 microns are ~20 orders of magnitude larger than blackbody expectations. We report a search of the AllWISE Source Catalog for Planet Nine, but find no detection.

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.

Jupiter's atmosphere has been observed to be depleted in helium (Yatm~0.24), suggesting active helium sedimentation in the interior. This is accounted for in standard Jupiter structure and evolution models through the assumption of an outer, He-depleted envelope that is separated from the He-enriched deep interior by a sharp boundary. Here we aim to develop a model for Jupiter's inhomogeneous thermal evolution that relies on a more self-consistent description of the internal profiles of He abundance, temperature, and heat flux. We make use of recent numerical simulations on H/He demixing, and on layered (LDD) and oscillatory (ODD) double diffusive convection, and assume an idealized planet model composed of a H/He envelope and a massive core. A general framework for the construction of interior models with He rain is described. Despite, or perhaps because of, our simplifications made we find that self-consistent models are rare. For instance, no model for ODD convection is found. We modify the H/He phase diagram of Lorenzen et al. to reproduce Jupiter's atmospheric helium abundance and examine evolution models as a function of the LDD layer height, from those that prolong Jupiter's cooling time to those that actually shorten it. Resulting models that meet the luminosity constraint have layer heights of about 0.1-1 km, corresponding to ~10,-20,000 layers in the rain zone between ~1 and 3-4.5 Mbars. Present limitations and directions for future work are discussed, such as the formation and sinking of He droplets.

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.

Remote sensing observations meet some limitations when used to study the bulk atmospheric composition of the giant planets of our solar system. A remarkable example of the superiority of in situ probe measurements is illustrated by the exploration of Jupiter, where key measurements such as the determination of the noble gases abundances and the precise measurement of the helium mixing ratio have only been made available through in situ measurements by the Galileo probe. This paper describes the main scientific goals to be addressed by the future in situ exploration of Saturn placing the Galileo probe exploration of Jupiter in a broader context and before the future probe exploration of the more remote ice giants. In situ exploration of Saturn's atmosphere addresses two broad themes that are discussed throughout this paper: first, the formation history of our solar system and second, the processes at play in planetary atmospheres. In this context, we detail the reasons why measurements of Saturn's bulk elemental and isotopic composition would place important constraints on the volatile reservoirs in the protosolar nebula. We also show that the in situ measurement of CO (or any other disequilibrium species that is depleted by reaction with water) in Saturn's upper troposphere would constrain its bulk O/H ratio. We highlight the key measurements required to distinguish competing theories to shed light on giant planet formation as a common process in planetary systems with potential applications to most extrasolar systems. In situ measurements of Saturn's stratospheric and tropospheric dynamics, chemistry and cloud-forming processes will provide access to phenomena unreachable to remote sensing studies. Different mission architectures are envisaged, which would benefit from strong international collaborations.

We perform modeling investigations to aid in understanding the atmospheres and composition of small planets of ~2-4 Earth radii, which are now known to be common in our galaxy. GJ 1214b is a well studied example whose atmospheric transmission spectrum has been observed by many investigators. Here we take a step back from GJ 1214b to investigate the role that planetary mass, composition, and temperature play in impacting the transmission spectra of these low-mass low-density (LMLD) planets. Under the assumption that these planets accrete modest hydrogen-dominated atmospheres and planetesimals, we use population synthesis models to show that predicted metal enrichments of the H/He envelope are high, with metal mass fraction Z_env values commonly 0.6 to 0.9, or ~100 to 400+ times solar. The high mean molecular weight of such atmospheres (\mu ~ 5-12) would naturally help to flatten the transmission spectrum of most LMLD planets. The high metal abundance would also provide significant condensible material for cloud formation. It is known that the H/He abundance in Uranus and Neptune decreases with depth, and we show that atmospheric evaporation of LMLD planets could expose atmospheric layers with gradually higher Z_env. However, values of Z_env close to solar composition can also arise, so diversity should be expected. Photochemically produced hazes, potentially due to methane photolysis, are another possibility for obscuring transmission spectra. Such hazes may not form above T_eq of ~800-1100 K, which is testable if such warm, otherwise low mean molecular weight atmospheres are stable against atmospheric evaporation. We find that available transmission data are consistent with relatively high mean molecular weight atmospheres for GJ 1214b and "warm Neptune" GJ 436b. We examine future prospects for characterizing GJ 1214b with Hubble and the James Webb Space Telescope.

Neptune-sized extrasolar planets that orbit relatively close to their host stars -- often called "hot Neptunes" -- are common within the known population of exoplanets and planetary candidates. Similar to our own Uranus and Neptune, inefficient accretion of nebular gas is expected produce hot Neptunes whose masses are dominated by elements heavier than hydrogen and helium. At high atmospheric metallicities of 10-10,000x solar, hot Neptunes will exhibit an interesting continuum of atmospheric compositions, ranging from more Neptune-like, H2-dominated atmospheres to more Venus-like, CO2-dominated atmospheres. We explore the predicted equilibrium and disequilibrium chemistry of generic hot Neptunes and find that the atmospheric composition varies strongly as a function of temperature and bulk atmospheric properties such as metallicity and the C/O ratio. Relatively exotic H2O, CO, CO2, and even O2-dominated atmospheres are possible for hot Neptunes. We apply our models to the case of GJ 436b, where we find that a CO-rich, CH4-poor atmosphere can be a natural consequence of a very high atmospheric metallicity. From comparisons of our results with Spitzer eclipse data for GJ 436b, we conclude that although the spectral fit from the high-metallicity forward models is not quite as good as the fit obtained from pure retrieval methods, the atmospheric composition predicted by these forward models is more physically and chemically plausible. High-metallicity atmospheres (orders of magnitude in excess of solar) should therefore be considered as a possibility for GJ 436b and other hot Neptunes.

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.

Jupiter is expected to pulsate in a spectrum of acoustic modes and recent re-analysis of a spectroscopic time series has identified a regular pattern in the spacing of the frequencies \citep{gaulme2011}. This exciting result can provide constraints on gross Jovian properties and warrants a more in-depth theoretical study of the seismic structure of Jupiter. With current instrumentation, such as the SYMPA instrument \citep{schmider2007} used for the \citet{gaulme2011} analysis, we assume that, at minimum, a set of global frequencies extending up to angular degree $\ell=25$ could be observed. In order to identify which modes would best constrain models of Jupiter's interior and thus help motivate the next generation of observations, we explore the sensitivity of derived parameters to this mode set. Three different models of the Jovian interior are computed and the theoretical pulsation spectrum from these models for $\ell\leq 25$ is obtained. We compute sensitivity kernels and perform linear inversions to infer details of the expected discontinuities in the profiles in the Jovian interior. We find that the amplitude of the sound-speed jump of a few percent in the inner/outer envelope boundary seen in two of the applied models should be reasonably inferred with these particular modes. Near the core boundary where models predict large density discontinuities, the location of such features can be accurately measured, while their amplitudes have more uncertainty. These results suggest that this mode set would be sufficient to infer the radial location and strength of expected discontinuities in Jupiter's interior, and place strong constraints on the core size and mass. We encourage new observations to detect these Jovian oscillations.

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.

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

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.

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.

More than 80 giant planets are known by mass and radius. Their interior structure in terms of core mass, number of layers, and composition however is still poorly known. An overview is presented about the core mass Mcore and envelope mass of metals MZ in Jupiter as predicted by various equations of state. It is argued that the uncertainty about the true H/He EOS in a pressure regime where the gravitational moments J2 and J4 are most sensitive, i.e. between 0.5 and 4 Mbar, is in part responsible for the broad range Mcore=0-18 Mearth, MZ=0-38 Mearth, and Mcore+MZ=14-38 Mearth currently offered for Jupiter. We then compare the Jupiter models obtained when we only match J2 with the range of solutions for the exoplanet GJ436b, when we match an assumed tidal Love number k2 value.

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.

We compute grids of radiative-convective model atmospheres for Jupiter, Saturn, Uranus, and Neptune over a range of intrinsic fluxes and surface gravities. The atmosphere grids serve as an upper boundary condition for models of the thermal evolution of the planets. Unlike previous work, we customize these grids for the specific properties of each planet, including the appropriate chemical abundances and incident fluxes as a function of solar system age. Using these grids, we compute new models of the thermal evolution of the major planets in an attempt to match their measured luminosities at their known ages. Compared to previous work, we find longer cooling times, predominantly due to higher atmospheric opacity at young ages. For all planets, we employ simple "standard" cooling models that feature adiabatic temperature gradients in the interior H/He and water layers, and an initially hot starting point for the calculation of subsequent cooling. For Jupiter we find a model cooling age 10% longer than previous work, a modest quantitative difference. This may indicate that the hydrogen equation of state used here overestimates the temperatures in the deep interior of the planet. For Saturn we find a model cooling age 20% longer than previous work. However, an additional energy source, such as that due to helium phase separation, is still clearly needed. For Neptune, unlike in work from the 1980s and 1990s, we match the measured Teff of the planet with a model that also matches the planet's current gravity field. This is predominantly due to advances in the equation of state of water. This may indicate that the planet possesses no barriers to efficient convection in its deep interior. However, for Uranus, our models exacerbate the well-known problem that Uranus is far cooler than calculations predict, which could imply strong barriers to interior convective cooling.

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.

The discovery of CoRoT-7b, a planet of radius 1.68 +/- 0.09 R_E, mass 4.8 +/- 0.8 M_E and orbital period of 0.854 days demonstrates that small planets can orbit extremely close to their star. We use knowledge of hot Jupiters, mass loss estimates and models for the interior structure and evolution of planets to understand its composition, structure and evolution. The inferred mass and radius of CoRoT-7b are consistent with a rocky planet that would be depleted in iron relative to Earth. However, a one sigma increase in mass (5.6 M_E) and decrease in size (1.59 R_E) are compatible with an Earth-like composition (33% iron, 67% silicates). Alternatively, it is possible that CoRoT-7b contains a significant amount of volatiles. An equally good fit to the data is found for a vapor envelope equal to 3% (and up to 10%) by mass above an Earth-like nucleus. Because of its intense irradiation and small size, the planet cannot possess an envelope of H and He of more than 1/10,000 of its total mass. We show that the mass loss is significant (~ 10^11 g/s) and independent of planetary composition. This is because the hydrodynamical escape rate is independent of the atmosphere's mean molecular mass, and owing to the intense irradiation, even a bare rocky planet would possess an equilibrium vapor atmosphere thick enough to capture stellar UV photons. This escape rate rules out the possibility of a H-He envelope as it would escape in only ~1 Ma. A water vapor atmosphere would escape in ~ 1 Ga, and thus it is a plausible scenario. The origin of CoRoT-7b cannot be inferred from present observations: It may have formed rocky; or be the remnant of a Uranus-like ice giant, or a gas giant with a small core that was stripped of its gaseous envelope.

We discuss our current understanding of the interior structure and thermal evolution of giant planets. This includes the gas giants, such as Jupiter and Saturn, that are primarily composed of hydrogen and helium, as well as the "ice giants," such as Uranus and Neptune, which are primarily composed of elements heavier than H/He. The effect of different hydrogen equations of state (including new first-principles computations) on Jupiter's core mass and heavy element distribution is detailed. This variety of the hydrogen equations of state translate into an uncertainty in Jupiter's core mass of 18 M_Earth. For Uranus and Neptune we find deep envelope metallicities up to 0.95, perhaps indicating the existence of an eroded core, as also supported by their low luminosity. We discuss the results of simple cooling models of our solar system's planets, and show that more complex thermal evolution models may be necessary to understand their cooling history. We review how measurements of the masses and radii of the ~50 transiting extrasolar giant planets are changing our understanding of giant planets. In particular a fraction of these planets appear to be larger than can be accommodated by standard models of planetary contraction. We review the proposed explanations for the radii of these planets. We also discuss very young giant planets, which are being directly imaged with ground- and space-based telescopes.

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.