The Evolution and Internal Structure of Jupiter and Saturn with Compositional Gradients

The internal structure of gas giant planets may be more complex than the commonly assumed core-envelope structure with an adiabatic temperature profile. Different primordial internal structures as well as various physical processes can lead to non-homogenous compositional distributions. A non-homogenous internal structure has a significant impact on the thermal evolution and final structure of the planets. In this paper, we present alternative structure and evolution models for Jupiter and Saturn allowing for non-adiabatic primordial structures and the mixing of heavy elements by convection as these planets evolve. We present the evolution of the planets accounting for various initial composition gradients, and in the case of Saturn, include the formation of a helium-rich region as a result of helium rain. We investigate the stability of regions with composition gradients against convection, and find that the helium shell in Saturn remains stable and does not mix with the rest of the envelope. In other cases, convection mixes the planetary interior despite the existence of compositional gradients, leading to the enrichment of the envelope with heavy elements. We show that non-adiabatic structures (and cooling histories) for both Jupiter and Saturn are feasible. The interior temperatures in that case are much higher that for standard adiabatic models. We conclude that the internal structure is directly linked to the formation and evolution history of the planet. These alternative internal structures of Jupiter and Saturn should be considered when interpreting the upcoming Juno and Cassini data.Comment: accepted for publication in Ap

Metallicity of the Massive Protoplanets Around HR 8799 If Formed by Gravitational Instability

The final composition of giant planets formed as a result of gravitational instability in the disk gas depends on their ability to capture solid material (planetesimals) during their 'pre-collapse' stage, when they are extended and cold, and contracting quasi-statically. The duration of the pre-collapse stage is inversely proportional roughly to the square of the planetary mass, so massive protoplanets have shorter pre-collapse timescales and therefore limited opportunity for planetesimal capture. The available accretion time for protoplanets with masses of 3, 5, 7, and 10 Jupiter masses is found to be 7.82E4, 2.62E4, 1.17E4 and 5.67E3 years, respectively. The total mass that can be captured by the protoplanets depends on the planetary mass, planetesimal size, the radial distance of the protoplanet from the parent star, and the local solid surface density. We consider three radial distances, 24, 38, and 68 AU, similar to the radial distances of the planets in the system HR 8799, and estimate the mass of heavy elements that can be accreted. We find that for the planetary masses usually adopted for the HR 8799 system, the amount of heavy elements accreted by the planets is small, leaving them with nearly stellar compositions.Comment: accepted for publication in Icaru

Explaining the low luminosity of Uranus: A self-consistent thermal and structural evolution

The low luminosity of Uranus is a long-standing challenge in planetary science. Simple adiabatic models are inconsistent with the measured luminosity, which indicates that Uranus is non-adiabatic because it has thermal boundary layers and/or conductive regions. A gradual composition distribution acts as a thermal boundary to suppress convection and slow down the internal cooling. Here we investigate whether composition gradients in the deep interior of Uranus can explain its low luminosity, the required composition gradient, and whether it is stable for convective mixing on a timescale of some billion years. We varied the primordial composition distribution and the initial energy budget of the planet, and chose the models that fit the currently measured properties (radius, luminosity, and moment of inertia) of Uranus. We present several alternative non-adiabatic internal structures that fit the Uranus measurements. We found that convective mixing is limited to the interior of Uranus, and a composition gradient is stable and sufficient to explain its current luminosity. As a result, the interior of Uranus might still be very hot, in spite of its low luminosity. The stable composition gradient also indicates that the current internal structure of Uranus is similar to its primordial structure. Moreover, we suggest that the initial energy content of Uranus cannot be greater than 20% of its formation (accretion) energy. We also find that an interior with a mixture of ice and rock, rather than separated ice and rock shells, is consistent with measurements, suggesting that Uranus might not be "differentiated". Our models can explain the luminosity of Uranus, and they are also consistent with its metal-rich atmosphere and with the predictions for the location where its magnetic field is generated.Comment: 10 pages, 7 figures, accepted for publication in A&

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

The link between infall location, early disc size, and the fraction of self-gravitationally fragmenting discs

Context. Many protoplanetary discs are self-gravitating early in their lives. If they fragment under their own gravity, they form bound gaseous clumps that can evolve to become giant planets. Today, the fraction of discs that undergo fragmentation, and therefore also the frequency of conditions that may lead to giant planet formation via gravitational instability, is still unknown. Aims. We study the formation and evolution of a large number of star-disc systems, focusing on the early sizes of the discs and their likelihood to fragment. We investigate how the fraction of discs that fragments depends on the disc-size distribution at early times. Methods. We performed a population synthesis of discs from formation to dispersal. Whilst varying the infall radius, we study the relationship between early disc size and fragmentation. Furthermore, we investigate how stellar accretion heating affects the fragmentation fraction. Results. We find that discs fragment only if they become sufficiently large early in their lives. This size depends sensitively on where mass is added to the discs during the collapse of their parent molecular cloud core. Infall locations derived from pure hydrodynamic and non-ideal magnetised collapse simulations lead to large and small discs, respectively, and 22 and 0% fragmentation fractions, respectively, in populations representative of the initial mass function; however, the resulting synthetic disc size distribution is larger and smaller, respectively, than the observed Class 0 disc size distribution. By choosing intermediate infall locations, leading to a synthetic disc size distribution that is in agreement with the observed one, we find a fragmentation fraction of between 0.1 and 11%, depending on the efficiency of stellar accretion heating of the discs. Conclusions. We conclude that the frequency of fragmentation is strongly affected by the early formation process of the disc and its interaction with the star. The early disc size is mainly determined by the infall location during the collapse of the molecular cloud core and controls the population-wide frequency of fragmentation. Stellar accretion heating also plays an important role in fragmentation and must be studied further. Our work is an observationally informed step towards a prediction of the frequency of giant planet formation by gravitational instability. Upcoming observations and theoretical studies will further our understanding of the formation and early evolution of discs in the near future. This will eventually allow us to understand how infall, disc morphology, giant planet formation via gravitational instability, and the observed extrasolar planet population are linked

Revisited Mass-Radius relations for exoplanets below 120 Earth masses

The masses and radii of exoplanets are fundamental quantities needed for their characterisation. Studying the different populations of exoplanets is important for understanding the demographics of the different planetary types, which can then be linked to planetary formation and evolution. We present an updated exoplanet catalog based on reliable, robust and as much as possible accurate mass and radius measurements of transiting planets up to 120 $M_{\oplus}$. The resulting mass-radius (M-R) diagram shows two distinct populations, corresponding to rocky and volatile-rich exoplanets which overlap in both mass and radius. The rocky exoplanet population shows a relatively small density variability and ends at mass of $\sim25 M_{\oplus}$, possibly indicating the maximum core mass that can be formed. We use the composition line of pure-water to separate the two populations, and infer two new empirical M-R relations based on this data: $M = (0.9 \pm 0.06) \ R^{(3.45 \pm 0.12)}$ for the rocky population, and $M = (1.74 \pm 0.38) \ R^{(1.58 \pm 0.10)}$ for the volatile-rich population. While our results for the two regimes are in agreement with previous studies, the new M-R relations better match the population in the transition-region from rocky to volatile-rich exoplanets, which correspond to a mass range of 5-25 $M_{\oplus}$ and a radius range of 2-3 $R_{\oplus}$.Comment: 13 pages, 5 figures. Accepted for publication in A&

A generalized bayesian inference method for constraining the interiors of super Earths and sub-Neptunes

We aim to present a generalized Bayesian inference method for constraining interiors of super Earths and sub-Neptunes. Our methodology succeeds in quantifying the degeneracy and correlation of structural parameters for high dimensional parameter spaces. Specifically, we identify what constraints can be placed on composition and thickness of core, mantle, ice, ocean, and atmospheric layers given observations of mass, radius, and bulk refractory abundance constraints (Fe, Mg, Si) from observations of the host star's photospheric composition. We employed a full probabilistic Bayesian inference analysis that formally accounts for observational and model uncertainties. Using a Markov chain Monte Carlo technique, we computed joint and marginal posterior probability distributions for all structural parameters of interest. We included state-of-the-art structural models based on self-consistent thermodynamics of core, mantle, high-pressure ice, and liquid water. Furthermore, we tested and compared two different atmospheric models that are tailored for modeling thick and thin atmospheres, respectively. First, we validate our method against Neptune. Second, we apply it to synthetic exoplanets of fixed mass and determine the effect on interior structure and composition when (1) radius, (2) atmospheric model, (3) data uncertainties, (4) semi-major axes, (5) atmospheric composition (i.e., a priori assumption of enriched envelopes versus pure H/He envelopes), and (6) prior distributions are varied. Our main conclusions are: [...]Comment: Astronomy & Astrophysics, 597, A37, 17 pages, 11 figure

On the Evolution and Survival of Protoplanets Embedded in a Protoplanetary Disk

We model the evolution of a Jupiter-mass protoplanet formed by the disk instability mechanism at various radial distances accounting for the presence of the disk. Using three different disk models, it is found that a newly-formed Jupiter-mass protoplanet at radial distance of $\lesssim$ 5-10 AU cannot undergo a dynamical collapse and evolve further to become a gravitational bound planet. We therefore conclude that {\it giant planets, if formed by the gravitational instability mechanism, must form and remain at large radial distances during the first $\sim$ 10$^5-10^6$ years of their evolution}. The minimum radial distances in which protoplanets of 1 Saturn-mass, 3 and 5 Jupiter-mass protoplanets can evolve using a disk model with $\dot{M}=10^{-6} M_{Sun}/yr$ and $\alpha=10^{-2}$ are found to be 12, 9, and 7 AU, respectively. The effect of gas accretion on the planetary evolution of a Jupiter-mass protoplanet is also investigated. It is shown that gas accretion can shorten the pre-collapse timescale substantially. Our study suggests that the timescale of the pre-collapse stage does not only depend on the planetary mass, but is greatly affected by the presence of the disk and efficient gas accretion.Comment: 26 pages, 2 tables, 10 figures. Accepted for publication in Ap