Oligarchic growth of giant planets

Runaway growth ends when the largest protoplanets dominate the dynamics of the planetesimal disk; the subsequent self-limiting accretion mode is referred to as “oligarchic growth.” Here, we begin by expanding on the existing analytic model of the oligarchic growth regime. From this, we derive global estimates of the planet formation rate throughout a protoplanetary disk. We find that a relatively high-mass protoplanetary disk (∼10 × minimum-mass) is required to produce giant planet core-sized bodies (∼10 M_⊕) within the lifetime of the nebular gas (?10 million years). However, an implausibly massive disk is needed to produce even an Earth mass at the orbit of Uranus by 10 Myrs. Subsequent accretion without the dissipational effect of gas is even slower and less efficient. In the limit of noninteracting planetesimals, a reasonable-mass disk is unable to produce bodies the size of the Solar System’s two outer giant planets at their current locations on any timescale; if collisional damping of planetesimal random velocities is sufficiently effective, though, it may be possible for a Uranus/Neptune to form in situ in less than the age of the Solar System. We perform numerical simulations of oligarchic growth with gas and find that protoplanet growth rates agree reasonably well with the analytic model as long as protoplanet masses are well below their estimated final masses. However, accretion stalls earlier than predicted, so that the largest final protoplanet masses are smaller than those given by the model. Thus the oligarchic growth model, in the form developed here, appears to provide an upper limit for the efficiency of giant planet formation.     

The internal activity and thermal evolution of Earth-like planets

We model the internal thermal evolution of planets with Earth-like composition and masses ranging from 0.1 to 10 Earth masses over a period of 10 billion years. We also characterize the internal activity of the planets by the velocity of putative tectonic plates, the rate at which mantle material is processed through melting zones, and the time taken to process one mantle mass. The more massive the planet the larger its processing rate (?), which scales approximately as ?∝M0.8-1.0. The processing times for all the planets increase with time as they cool and become less active. As would be expected, the surface heat flow scales with planet mass. All planets have similar declines in mantle temperature except for the largest, in which pressure effects cause a larger decline. The larger planets have higher mantle temperatures over all times. The less massive the planet, the larger the decrease in core temperature with time. The core heat flow is also found to decrease more rapidly for smaller planet masses. Finally, rough predictions are made for the time required to generate an atmosphere from estimates of the time to degas water and carbon dioxide in mantle melting zones. The degassing times depend strongly on the initial temperature of the planet, but for the temperatures used in our model all the planets degas within ∼32 Ma after their formation.     

Tidal interactions and disruptions of giant planets on highly eccentric orbits

We calculate the evolution of planets undergoing a strong tidal encounter using smoothed particle hydrodynamics (SPH), for a range of periastron separations. We find that outside the Roche limit, the evolution of the planet is well-described by the standard model of linear, non-radial, adiabatic oscillations. If the planet passes within the Roche limit at periastron, however, mass can be stripped from it, but in no case do we find enough energy transferred to the planet to lead to complete disruption. In light of the three new extrasolar planets discovered with periods shorter than two days, we argue that the shortest-period cases observed in the period-mass relation may be explained by a model whereby planets undergo strong tidal encounters with stars, after either being scattered by dynamical interactions into highly eccentric orbits, or tidally captured from nearly parabolic orbits. Although this scenario does provide a natural explanation for the edge found for planets at twice the Roche limit, it does not explain how such planets will survive the inevitable expansion that results from energy injection during tidal circularization.     

The full set of gas giant structures: I: On the origin of planetary masses and the planetary initial mass function

Context

Current planet search programs are detecting extrasolar planets at a rate of 60 planets per year. These planets show more diverse properties than was expected.

Aims

We try to get an overview of possible gas giant (proto-) planets for a full range of orbital periods and stellar masses. This allows the prediction of the full range of possible planetary properties which might be discovered in the near future.

Methods

We calculate the purely hydrostatic structure of the envelopes of proto-planets that are embedded in protoplanetary disks for all conceivable locations: combinations of different planetesimal accretion rates, host star masses, and orbital separations. At each location all hydrostatic equilibrium solutions to the planetary structure equations are determined by variation of core mass and pressure over many orders of magnitude. For each location we analyze the distribution of planetary masses.

Results

We get a wide spectrum of core-envelope structures. However, practically all calculated proto-planets are in the planetary mass range. Furthermore, the planet masses show a characteristic bimodal, sometimes trimodal, distribution. For the first time, we identify three physical processes that are responsible for the three characteristic planet masses: self-gravity in the Hill sphere, compact objects, and a region of very low adiabatic pressure gradient in the hydrogen equation of state. Using these processes, we can explain the dependence of the characteristic masses on the planet’s location: orbital period, host star mass, and planetesimal accretion rate (luminosity). The characteristic mass caused by the self-gravity effect at close proximity to the host star is typically one Neptune mass, thus producing the so-called hot Neptunes.

Conclusions

Our results suggest that hot Jupiters with orbital period less than 64 days (the exact location of the boundary depends on stellar type and accretion rate) have quite distinct properties which we expect to be reflected in a different mass distribution of these planets when compared to the “normal” planetary population. We use our theoretical survey to produce an upper mass limit for embedded planets: the maximum embedded equilibrium mass (MEEM). This naturally explains the lack of high mass planets between 3 and 64 days orbital period.     

A hybrid scenario for gas giant planet formation in rings

The core-accretion mechanism for gas giant formation may be too slow to create all observed gas giant planets during reasonable gas disk lifetimes, but it has yet to be firmly established that the disk instability model can produce permanent bound gaseous protoplanets under realistic conditions. Based on our recent simulations of gravitational instabilities in disks around young stars, we suggest that, even if instabilities due to disk self-gravity do not produce gaseous protoplanets directly, they may create persistent dense rings that are conducive to accelerated growth of gas giants through core accretion. The rings occur at and near the boundary between stable and unstable regions of the disk and appear to be produced by resonances with discrete spiral modes on the unstable side.     

The formation and habitability of terrestrial planets in the presence of close-in giant planets

‘Hot jupiters,’ giant planets with orbits very close to their parent stars, are thought to form farther away and migrate inward via interactions with a massive gas disk. If a giant planet forms and migrates quickly, the planetesimal population has time to re-generate in the lifetime of the disk and terrestrial planets may form [P.J. Armitage, A reduced efficiency of terrestrial planet formation following giant planet migration, Astrophys. J. 582 (2003) L47-L50]. We present results of simulations of terrestrial planet formation in the presence of hot/warm jupiters, broadly defined as having orbital radii ?0.5 AU. We show that terrestrial planets similar to those in the Solar System can form around stars with hot/warm jupiters, and can have water contents equal to or higher than the Earth's. For small orbital radii of hot jupiters (e.g., 0.15, 0.25 AU) potentially habitable planets can form, but for semi-major axes of 0.5 AU or greater their formation is suppressed. We show that the presence of an outer giant planet such as Jupiter does not enhance the water content of the terrestrial planets, but rather decreases their formation and water delivery timescales. We speculate that asteroid belts may exist interior to the terrestrial planets in systems with close-in giant planets.     

The hydrogen ortho-to-para ratio in the stratospheres of the giant planets

Observations of the H₂ S(0) and S(1) quadrupole lines in the four giant planets by the short-wavelength spectrometer of the Infrared Space Observatory are analyzed. These lines probe pressure levels located between 10 and 1 mbar and allow us to determine the stratospheric hydrogen para fraction for the first time. In Jupiter and Saturn, the stratospheric para fraction is close to its tropopause value. In the stratosphere of these planets as well as in Neptune’s, the para fraction presents a significant departure from thermodynamic equilibrium. This situation results from a lagged conversion between the ortho and the para states as molecular hydrogen is transported upward under the influence of turbulent eddy diffusion. In contrast, the uranian stratosphere lies close to thermodynamic equilibrium. The magnitude of the departure from thermodynamic equilibrium appears to be anti-correlated with the amount of stratospheric aerosols. To validate this assumption, we estimate the hydrogen equilibration time with a one-dimensional diffusion model for different conversion processes in the gas phase or on aerosols. The comparison between our results and the tropospheric estimates from Conrath et al. (1998, Icarus,135, 501-517) shows that paramagnetic conversion on aerosols matches the estimated tropospheric and stratospheric relaxation times in the four giant planets. In contrast, paramagnetic conversion in the gas phase can only explain the relaxation times measured in Jupiter and Saturn atmospheres. This situation provides quantitative evidence for an equilibration mechanism dominated by conversion on aerosols.     

Internal structure of massive terrestrial planets

Planetary formation models predict the existence of massive terrestrial planets and experiments are now being designed that should succeed in discovering them and measuring their masses and radii. We calculate internal structures of planets with one to ten times the mass of the Earth (Super-Earths) to obtain scaling laws for total radius, mantle thickness, core size and average density as a function of mass. We explore different compositions and obtain a scaling law of R∝M0.267-0.272 for Super-Earths. We also study a second family of planets, Super-Mercuries with masses ranging from one mercury-mass to ten mercury-masses with similar composition to the Earth's but with a larger core mass fraction. We explore the effect of surface temperature and core mass fraction on the scaling laws for these planets. The scaling law obtained for the Super-Mercuries is R∝M∼0.3.     

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 (∼1 AU) and internal temperatures are low (<2000 K). The main source of opacity in these objects is dust grains. We investigate two problems involving the effect of opacity on the evolution of isolated, non-accreting 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.     

The search for extrasolar giant planets using integral field spectroscopy: Simulations

We present here simulations of extrasolar planets detections obtained using a combination of extreme adaptive optics and integral field spectroscopy. The simulation code, written for IDL, provides images and, in particular, spectra, taking into account realistic Speckle Noise, AO correction effects and specific instrumental features. A detailed study has been done for ESO VLT telescopes (8.2 m), within the Phase A of the CHEOPS project, but the code is particularly flexible and can be updated for larger telescope diameters (ELTs) in order to give a realistic estimate of the detection limits, for giant telescopes, in standard conditions of seeing.     

Mass-radius curve for extrasolar Earth-like planets and ocean planets

By comparison with the Earth-like planets and the large icy satellites of the Solar System, one can model the internal structure of extrasolar planets. The input parameters are the composition of the star (Fe/Si and Mg/Si), the Mg content of the mantle (Mg# = Mg/[Mg + Fe]), the amount of H₂O and the total mass of the planet. Equation of State (EoS) of the different materials that are likely to be present within such planets have been obtained thanks to recent progress in high-pressure experiments. They are used to compute the planetary radius as a function of the total mass. Based on accretion models and data on planetary differentiation, the internal structure is likely to consist of an iron-rich core, a silicate mantle and an outer silicate crust resulting from magma formation in the mantle. The amount of H₂O and the surface temperature control the possibility for these planets to harbor an ocean. In preparation to the interpretation of the forthcoming data from the CNES led CoRoT (Convection Rotation and Transit) mission and from ground-based observations, this paper investigates the relationship between radius and mass. If H₂O is not an important component (less than 0.1%) of the total mass of the planet, then a relation (R/REarth)=ab(M/MEarth) is calculated with (a,b)=(1,0.306) and (a,b)=(1,0.274) for ¹⁰⁻²MEarth<M<MEarth and MEarth<M<10MEarth, respectively. Calculations for a planet that contains 50% H₂O suggest that the radius would be more than 25% larger than that based on the Earth-like model, with (a,b)=(1.258,0.302) for ¹⁰⁻²MEarth<M<MEarth and (a,b)=(1.262,0.275) for MEarth<M<10MEarth, respectively. For a surface temperature of 300 K, the thickness of the ocean varies from 150 to 50 km for planets 1 to 10 times the Earth's mass, respectively. Application of this algorithm to bodies of the Solar System provides not only a good fit to most terrestrial planets and large icy satellites, but also insights for discussing future observations of exoplanets.     

A 2D study of the effects of the size of a solid core on the equatorial flow in giant planets

Zonal winds simulated in two-dimensional computer models of turbulent convection in the equatorial plane of giant planets have greater surface amplitudes for cases with smaller solid cores, and therefore larger buoyancy driving, all other properties being equal. This differential rotation in radius is maintained by the convergence of angular momentum flux, which occurs because of the convective flow that develops due to the effects of planetary rotation and density stratification. The superposition of the convective flow and the stronger zonal flow produces wave-like, instead of cellular convection.     

The maximal runaway temperature of Earth-like planets

In Simpson’s (Simpson, G.C. [1927]. Mem. R. Meteorol. Soc. II (16), 69–95) classical derivation of the temperature of the Earth in the semi-gray model, the surface temperature diverges as the fourth root of the thermal radiation’s optical depth. No resolution to this apparent paradox was yet obtained under the strict semi-gray approximation. Using this approximation and a simplified approach, we study the saturation of the runaway greenhouse effect.First we generalize the problem of the semi-gray model to cases in which a non-negligible fraction of the stellar radiation falls on the long-wavelength range, and/or that the planetary long-wavelength emission penetrates into the transparent short wavelength domain of the absorption.Second, applying the most general assumptions and independently of any particular properties of an absorber, we show that the greenhouse effect saturates and that any Earth-like planet has a maximal temperature which depends on the type of and distance to its main-sequence star, its albedo and the primary atmospheric components which determine the cutoff frequency below which the atmosphere is optically thick. For example, a hypothetical convection-less planet similar to Venus, that is optically thin in the visible, could have at most a surface temperature of 1200–1300 K irrespective of the nature of the greenhouse gas.We show that two primary mechanisms are responsible for the saturation of the runaway greenhouse effect, depending on the value of λ_cut, the wavelength above which the atmosphere becomes optically thick. Unless λ_cut is small and resides in the optical region, saturation is achieved by radiating the thermal flux of the planet through the short wavelength tail of the thermal distribution. This has an interesting observational implication, the radiation from such a planet should be skewed towards the NIR. Otherwise, saturation takes place by radiating through windows in the FIR.     

The origin and nature of Neptune-like planets orbiting close to solar type stars

The sample of known exoplanets is strongly biased to masses larger than the ones of the giant gaseous planets of the Solar System. Recently, the discovery of two extrasolar planets of considerably lower masses around the nearby Stars GJ 436 and ρ Cancri was reported. They are like our outermost icy giants, Uranus and Neptune, but in contrast, these new planets are orbiting at only some hundredth of the Earth-Sun distance from their host stars, raising several new questions about their origin and constitution. Here we report numerical simulations of planetary accretion that show, for the first time through N-body integrations that the formation of compact systems of Neptune-like planets close to the hosts stars could be a common by-product of planetary formation. We found a regime of planetary accretion, in which orbital migration accumulates protoplanets in a narrow region around the inner edge of the nebula, where they collide each other giving rise to Neptune-like planets. Our results suggest that, if a protoplanetary solar environment is common in the Galaxy, the discovery of a vast population of this sort of ‘hot cores’ should be expected in the near future.     

The deep wind structure of the giant planets: Results from an anelastic general circulation model

The giant gas planets have hot convective interiors, and therefore a common assumption is that these deep atmospheres are close to a barotropic state. Here we show using a new anelastic general circulation model that baroclinic vorticity contributions are not negligible, and drive the system away from an isentropic and therefore barotropic state. The motion is still aligned with the direction of the axis of rotation as in a barotropic rotating fluid, but the wind structure has a vertical shear with stronger winds in the atmosphere than in the interior. This shear is associated with baroclinic compressibility effects. Most previous convection models of giant planets have used the Boussinesq approximation, which assumes the density is constant in depth; however, Jupiter's actual density varies by four orders of magnitude through its deep molecular envelope. We therefore developed a new general circulation model (based on the MITgcm) that is anelastic and thereby incorporates this density variation. The model's geometry is a full 3D sphere down to a small inner core. It is nonhydrostatic, uses an equation of state suitable for hydrogen-helium mixtures (SCVH), and is driven by an internal heating profile. We demonstrate the effect of compressibility by comparing anelastic and Boussinesq cases. The simulations develop a mean state that is geostrophic and hydrostatic including the often neglected, but significant, vertical Coriolis contribution. This leads to modification of the standard thermal wind relation for a deep compressible atmosphere. The interior flow organizes in large cyclonically rotating columnar eddies parallel to the rotation axis, which drive upgradient angular momentum eddy fluxes, generating the observed equatorial superrotation. Heat fluxes align with the axis of rotation, and provide a mechanism for the transport of heat poleward, which can cause the observed flat meridional emission. We address the issue of over-forcing which is common in such convection models and analyze the dependence of our results on this; showing that the vertical wind structure is not very sensitive to the Rayleigh number. We also study the effect of rotation, showing how the transition from a rapidly to a slowly rotating system affects the dynamics.     

A new family of planets? “Ocean-Planets”

A new family of planets is considered which is in between the rocky terrestrial planets and the gaseous giants, “Ocean-Planets.” We present the possible formation, composition and internal structure of these putative planets. We consider their oceans, as well as their possible Exobiology interest. These exoplanets should be detectable by Space missions such as Eddington, Kepler, and possibly COROT (launch scheduled in 2006). They have a density lower than that of rocky planets. Their rather large radius would make them attractive targets for exoplanet spectroscopic missions such as Darwin/TPF, all the more because a robust biosignature appears to exist.     

Detecting the glint of starlight on the oceans of distant planets

We propose that astronomers will be eventually be able to discriminate between extrasolar Earth-like planets with surface oceans and those without using the shape of phase light curves in the visible and near-IR spectrum. We model the visible light curves of planets having Earth-like surfaces, seasons, and optically-thin atmospheres with idealized diffuse-scattering clouds. We show that planets partially covered by water will appear measurably brighter near crescent phase (relative to Lambertian planets) because of the efficient specular reflection (“glint”) of starlight incident on their surfaces at a highly oblique angle. Planets on orbits within 30° of edge-on orientation (50% of all planets) will show pronounced glint over a sizeable range of orbital longitudes, from quadrature to crescent, all outside the glare of their parent stars. Also, water-covered planets will appear darker than a Lambertian disk near full illumination. Finally, we show that planets with a mixed land/water surface will polarize the reflected signal by as much as 30-70%. These results suggest several new ways of directly identifying water on distant planets.     

Orbital stability of systems of closely-spaced planets

An investigation of the stability of systems of 1 M_⊕ (Earth-mass) bodies orbiting a Sun-like star has been conducted for virtual times reaching 10 billion years. For the majority of the tests, a symplectic integrator with a fixed timestep of between 1 and 10 days was employed; however, smaller timesteps and a Bulirsch-Stoer integrator were also selectively utilized to increase confidence in the results. In most cases, the planets were started on initially coplanar, circular orbits, and the longitudinal initial positions of neighboring planets were widely separated. The ratio of the semimajor axes of consecutive planets in each system was approximately uniform (so the spacing between consecutive planets increased slowly in terms of distance from the star). The stability time for a system was taken to be the time at which the orbits of two or more planets crossed. Our results show that, for a given class of system (e.g., three 1 M_⊕ planets), orbit crossing times vary with planetary spacing approximately as a power law over a wide range of separation in semimajor axis. Chaos tests indicate that deviations from this power law persist for changed initial longitudes and also for small but non-trivial changes in orbital spacing. We find that the stability time increases more rapidly at large initial orbital separations than the power-law dependence predicted from moderate initial orbital separations. Systems of five planets are less stable than systems of three planets for a specified semimajor axis spacing. Furthermore, systems of less massive planets can be packed more closely, being about as stable as 1 M_⊕ planets when the radial separation between planets is scaled using the mutual Hill radius. Finally, systems with retrograde planets can be packed substantially more closely than prograde systems with equal numbers of planets.     

Chemistry of atmospheres formed during accretion of the Earth and other terrestrial planets

We used chemical equilibrium and chemical kinetic calculations to model chemistry of the volatiles released by heating different types of carbonaceous, ordinary and enstatite chondritic material as a function of temperature and pressure. Our results predict the composition of atmospheres formed by outgassing during accretion of the Earth and other terrestrial planets. Outgassing of CI and CM carbonaceous chondritic material produces H₂O-rich (steam) atmospheres in agreement with the results of impact experiments. However, outgassing of other types of chondritic material produces atmospheres dominated by other gases. Outgassing of ordinary (H, L, LL) and high iron enstatite (EH) chondritic material yields H₂-rich atmospheres with CO and H₂O being the second and third most abundant gases. Outgassing of low iron enstatite (EL) chondritic material gives a CO-rich atmosphere with H₂, CO₂, and H₂O being the next most abundant gases. Outgassing of CV carbonaceous chondritic material gives a CO₂-rich atmosphere with H₂O being the second most abundant gas. Our results predict that the atmospheres formed during accretion of the Earth and Mars were probably H₂-rich unless the accreted material was dominantly CI and CM carbonaceous chondritic material. We also predict significant amounts of S, P, Cl, F, Na, and K in accretionary atmospheres at high temperatures (1500-2500 K). Finally, our results may be useful for interpreting spectroscopic observations of accreting extrasolar terrestrial planets.