Interior Models of Saturn: Including the Uncertainties in Shape and Rotation

The accurate determination of Saturn's gravitational coefficients by Cassini could provide tighter constrains on Saturn's internal structure. Also, occultation measurements provide important information on the planetary shape which is often not considered in structure models. In this paper we explore how wind velocities and internal rotation affect the planetary shape and the constraints on Saturn's interior. We show that within the geodetic approach (Lindal et al., 1985, ApJ, 90, 1136) the derived physical shape is insensitive to the assumed deep rotation. Saturn's re-derived equatorial and polar radii at 100 mbar are found to be 54,445 $\pm$10 km and 60,365$\pm$10 km, respectively. To determine Saturn's interior we use {\it 1 D} three-layer hydrostatic structure models, and present two approaches to include the constraints on the shape. These approaches, however, result in only small differences in Saturn's derived composition. The uncertainty in Saturn's rotation period is more significant: with Voyager's 10h39mns period, the derived mass of heavy elements in the envelope is 0-7 M$_{\oplus}$. With a rotation period of 10h32mns, this value becomes $<4$ $M_{\oplus}$, below the minimum mass inferred from spectroscopic measurements. Saturn's core mass is found to depend strongly on the pressure at which helium phase separation occurs, and is estimated to be 5-20 M$_{\oplus}$. Lower core masses are possible if the separation occurs deeper than 4 Mbars. We suggest that the analysis of Cassini's radio occultation measurements is crucial to test shape models and could lead to constraints on Saturn's rotation profile and departures from hydrostatic equilibrium.Comment: Accepted for publication in Ap

The Effects of Metallicity, and Grain Growth and Settling on the Early Evolution of Gaseous Protoplanets

Giant protoplanets formed by gravitational instability in the outer regions of circumstellar disks go through an early phase of quasi-static contraction during which radii are large and internal temperatures are low. The main source of opacity in these objects is dust grains. We investigate two problems involving the effect of opacity on the evolution of planets of 3, 5, and 7 M_J. First, we pick three different overall metallicities for the planet and simply scale the opacity accordingly. We show that higher metallicity results in slower contraction as a result of higher opacity. It is found that the pre-collapse time scale is proportional to the metallicity. In this scenario, survival of giant planets formed by gravitational instability is predicted to be more likely around low-metallicity stars, since they evolve to the point of collapse to small size on shorter time scales. But metal-rich planets, as a result of longer contraction times, have the best opportunity to capture planetesimals and form heavy-element cores. Second, we investigate the effects of opacity reduction as a result of grain growth and settling, for the same three planetary masses and for three different values of overall metallicity. When these processes are included, the pre-collapse time scale is found to be of order 1000 years for the three masses, significantly shorter than the time scale calculated without these effects. In this case the time scale is found to be relatively insensitive to planetary mass and composition. However, the effects of planetary rotation and accretion of gas and dust, which could increase the timescale, are not included in the calculation. The short time scale we find would preclude metal enrichment by planetesimal capture, as well as heavy-element core formation, over a large range of planetary masses and metallicities.Comment: 22 pages, accepted to Icaru

Methane Planets and their Mass-Radius Relation

Knowledge of both the mass and radius of an exoplanet allows us to estimate its mean density, and therefore its composition. Exoplanets seem to fill a very large parameter space in terms of mass and composition, and unlike the solar-system's planets, exoplanets also have intermediate masses (~ 5 - 50 M_Earth) with various densities. In this letter, we investigate the behavior of the Mass-Radius relation for methane (CH_4) planets and show that when methane planets are massive enough (Mp >~ 15 M_Earth), the methane can dissociate and lead to a differentiated planet with a carbon core, a methane envelope, and a hydrogen atmosphere. The contribution of a rocky core to the behavior of CH_4 planet is considered as well. We also develop interior models for several detected intermediate-mass planets that could, in principle, be methane/methane-rich planets. The example of methane planets emphasizes the complexity of the Mass-Radius relation and the challenge involved in uniquely inferring the planetary composition.Comment: Published in ApJ

The formation of mini-Neptunes

Mini-Neptunes seem to be common planets. In this work we investigate the possible formation histories and predicted occurrence rates of mini-Neptunes assuming the planets form beyond the iceline. We consider pebble and planetesimal accretion accounting for envelope enrichment and two different opacity conditions. We find that the formation of mini-Neptunes is a relatively frequent output when envelope enrichment by volatiles is included, and that there is a "sweet spot" for mini-Neptune formation with a relatively low solid accretion rate of ~10^{-6} Earth masses per year. This rate is typical for low/intermediate-mass protoplanetary disks and/or disks with low metallicities. With pebble accretion, envelope enrichment and high opacity favor the formation of mini-Neptunes, with more efficient formation at large semi-major axes (~30 AU) and low disk viscosity. For planetesimal accretion, such planets can form also without enrichment, with the opacity being a key aspect in the growth history and favorable formation location. Finally, we show that the formation of Neptune-like planets remains a challenge for planet formation theories.Comment: Accepted for publication in Ap

The Fuzziness of Giant Planets’ Cores

Giant planets are thought to have cores in their deep interiors, and the division into a heavy-element core and hydrogen–helium envelope is applied in both formation and structure models. We show that the primordial internal structure depends on the planetary growth rate, in particular, the ratio of heavy elements accretion to gas accretion. For a wide range of likely conditions, this ratio is in one-to-one correspondence with the resulting post-accretion profile of heavy elements within the planet. This flux ratio depends sensitively on the assumed solid-surface density in the surrounding nebula. We suggest that giant planets' cores might not be distinct from the envelope and includes some hydrogen and helium, and the deep interior can have a gradual heavy-element structure. Accordingly, Jupiter's core may not be well defined. Accurate measurements of Jupiter's gravitational field by Juno could put constraints on Jupiter's core mass. However, as we suggest here, the definition of Jupiter's core is complex, and the core's physical properties (mass, density) depend on the actual definition of the core and on the planet's growth history

Two Empirical Regimes of the Planetary Mass-Radius Relation

Today, with the large number of detected exoplanets and improved measurements, we can reach the next step of planetary characterization. Classifying different populations of planets is not only important for our understanding of the demographics of various planetary types in the galaxy, but also for our understanding of planet formation. We explore the nature of two regimes in the planetary mass-radius (M-R) relation. We suggest that the transition between the two regimes of "small" and "large" planets, occurs at a mass of 124 \pm 7, M_Earth and a radius of 12.1 \pm 0.5, R_Earth. Furthermore, the M-R relation is R \propto M^{0.55\pm 0.02} and R \propto M^{0.01\pm0.02} for small and large planets, respectively. We suggest that the location of the breakpoint is linked to the onset of electron degeneracy in hydrogen, and therefore, to the planetary bulk composition. Specifically, it is the characteristic minimal mass of a planet which consists of mostly hydrogen and helium, and therefore its M-R relation is determined by the equation of state of these materials. We compare the M-R relation from observational data with the one derived by population synthesis calculations and show that there is a good qualitative agreement between the two samples.Comment: accepted for publication in A&