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

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

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

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&