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.     

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.     

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.     

On the relationship of long-period comets to known and unknown planets

Statistical aspects of the question of the dynamical relationship of long-period comets to giant planets (Saturn, Uranus, Neptune), dwarf planets (Pluto, 2003 EL₆₁, 2003 UB₃₁₃), and five hypothetical planets are investigated. Data for 859 and 888 comets (2005, 2006) with periods P > 200 years are used in the statistics. No comets of the Kreutz, Marsden, Kracht, and Meyer groups are considered. The minimum interorbital distances of comets from the listed planetary bodies are mainly analyzed. Effects testifying to the kinematical relationship of some of the comets to planets have been established by testing the cometary data relative to 67 planes. The distribution of the perihelia of long-period comets is also discussed.     

On the internal energy source of the large planets

Based on spacecraft observations of the giant planets we call attention to an apparent regular mass-luminosity relation for the large planets (giant planets plus Earth). It can provide important clues and constraints on the question of their energy source, including various relevant theories and their difficulties. We describe in some detail a new thermonuclear reaction mechanism, valid for plasma fusions at both super—high and low temperatures, and including the original thermonuclear mechanism as a particular case. It can explain not only the internal energy souce of the large planets, but also many anomalous fusion phenomena observed in the laboratory and in astronomical bodies. Although these appear in various guises, their basic principle is the same, that is, they can all be understood in terms of the new mechanism.     

On the existence of small comets and their interactions with planets

Arguments are presented for a substantial, unexplored population of comets with radii less than 1 km. Known examples confirm this population and extrapolation of any plausible size-distribution function indicates large numbers. However, their accurate numbers, orbital characteristics, and physical properties are unknown. Thus, even though the small comets may be the most frequent cometary bodies impacting the planets, a quantitative evaluation is not currently possible. We advocate an optimized, dedicated search program to characterize this population.Laboratory for Atmospheric and Space Physics, University of Colorado at BoulderLaboratory for Atmospheric and Space Physics, University of Colorado at BoulderLaboratory for Atmospheric and Space Physics, University of Colorado at Boulder     

Obliquity evolution of extrasolar terrestrial planets

We have investigated the obliquity evolution of terrestrial planets in habitable zones (at ∼1 AU) in extrasolar planetary systems, due to tidal interactions with their satellite and host star with wide varieties of satellite-to-planet mass ratio (m/Mp) and initial obliquity (γ₀), through numerical calculations and analytical arguments. The obliquity, the angle between planetary spin axis and its orbit normal, of a terrestrial planet is one of the key factors in determining the planetary surface environments. A recent scenario of terrestrial planet accretion implies that giant impacts of Mars-sized or larger bodies determine the planetary spin and form satellites. Since the giant impacts would be isotropic, tilted spins (sinγ₀∼1) are more likely to be produced than straight ones (sinγ₀∼0). The ratio m/Mp is dependent on the impact parameters and impactors' mass. However, most of previous studies on tidal evolution of the planet-satellite systems have focused on a particular case of the Earth-Moon systems in which m/Mp?0.0125 and γ₀∼10° or the two-body planar problem in which γ₀=0° and stellar torque is neglected. We numerically integrated the evolution of planetary spin and a satellite orbit with various m/Mp (from 0.0025 to 0.05) and γ₀ (from 0° to 180°), taking into account the stellar torques and precessional motions of the spin and the orbit. We start with the spin axis that almost coincides with the satellite orbit normal, assuming that the spin and the satellite are formed by one dominant impact. With initially straight spins, the evolution is similar to that of the Earth-Moon system. The satellite monotonically recedes from the planet until synchronous state between the spin period and the satellite orbital period is realized. The obliquity gradually increases initially but it starts decreasing down to zero as approaching the synchronous state. However, we have found that the evolution with initially tiled spins is completely different. The satellite's orbit migrates outward with almost constant obliquity until the orbit reaches the critical radius ∼10-20 planetary radii, but then the migration is reversed to inward one. At the reversal, the obliquity starts oscillation with large amplitude. The oscillation gradually ceases and the obliquity is reduced to ∼0° during the inward migration. The satellite eventually falls onto the planetary surface or it is captured at the synchronous state at several planetary radii. We found that the character change of precession about total angular momentum vector into that about the planetary orbit normal is responsible for the oscillation with large amplitude and the reversal of migration. With the results of numerical integration and analytical arguments, we divided the m/Mp-γ₀ space into the regions of the qualitatively different evolution. The peculiar tidal evolution with initially tiled spins give deep insights into dynamics of extrasolar planet-satellite systems and discussions of surface environments of the planets.     

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.     

On the possibility of earth-type habitable planets around 47 UMa

We investigate whether Earth-type habitable planets can in principle exist in the planetary system of 47 UMa. The system of 47 UMa consists of two Jupiter-size planets beyond the outer edge of the stellar habitable zone, and thus resembles our own Solar System most closely compared to all exosolar planetary systems discovered so far. Our study of habitability deliberately follows an Earth-based view according to the concept of Franck and colleagues, which assumes the long-term possibility of photosynthetic biomass production under geodynamic conditions. Consequently, a broad variety of climatological, biogeochemical, and geodynamical processes involved in the generation of photosynthesis-driven life conditions is taken into account. The stellar luminosity and the age of the star/planet system are of fundamental importance for planetary habitability. Our study considers different types of planetary continental growth models and takes into account a careful assessment of the stellar parameters. In the event of successful formation and orbital stability, two subjects of intense research, we find that Earth-type habitable planets around 47 UMa are in principle possible! The likelihood of those planets is increased if assumed that 47 UMa is relatively young (?6 Gyr) and has a relatively small stellar luminosity as permitted by the observational range of those parameters.     

Nonuniform cratering of the terrestrial planets

We estimate the impact flux and cratering rate as a function of latitude on the terrestrial planets using a model distribution of planet crossing asteroids and comets [Bottke, W.F., Morbidelli, A., Jedicke, R., Petit, J.-M., Levison, H.F., Michel, P., Metcalfe, T.S., 2002. Icarus 156, 399-433]. After determining the planetary impact probabilities as a function of the relative encounter velocity and encounter inclination, the impact positions are calculated analytically, assuming the projectiles follow hyperbolic paths during the encounter phase. As the source of projectiles is not isotropic, latitudinal variations of the impact flux are predicted: the calculated ratio between the pole and equator is 1.05 for Mercury, 1.00 for Venus, 0.96 for the Earth, 0.90 for the Moon, and 1.14 for Mars over its long-term obliquity variation history. By taking into account the latitudinal dependence of the impact velocity and impact angle, and by using a crater scaling law that depends on the vertical component of the impact velocity, the latitudinal variations of the cratering rate (the number of craters with a given size formed per unit time and unit area) is in general enhanced. With respect to the equator, the polar cratering rate is about 30% larger on Mars and 10% on Mercury, whereas it is 10% less on the Earth and 20% less on the Moon. The cratering rate is found to be uniform on Venus. The relative global impact fluxes on Mercury, Venus, the Earth and Mars are calculated with respect to the Moon, and we find values of 1.9, 1.8, 1.6, and 2.8, respectively. Our results show that the relative shape of the crater size-frequency distribution does not noticeably depend upon latitude for any of the terrestrial bodies in this study. Nevertheless, by neglecting the expected latitudinal variations of the cratering rate, systematic errors of 20-30% in the age of planetary surfaces could exist between equatorial and polar regions when using the crater chronology method.     

Planets and minor planets CCD observations with the CMT and the CMASF

CCD observations of planets and minor planets obtained while surveying the northern and southern hemispheres using the CMT on La Palma (Spain) and the CMASF at El Leoncito (Argentina) are presented.Also presented is the possibility of observing some Solar System bodies with these instruments in the future.     

On the width and shape of the corotation region for low-mass planets

We study the coorbital flow for embedded, low-mass planets. We provide a simple semi-analytic model for the corotation region, which is subsequently compared to high-resolution numerical simulations. The model is used to derive an expression for the half-width of the horseshoe region, x _s, which in the limit of zero softening is given by   x _s/ r _p= 1.68( q / h )^1/2  , where q is the planet to central star mass ratio, h is the disc aspect ratio and   r _p  is the orbital radius. This is in very good agreement with the same quantity measured from simulations. This result is used to show that horseshoe drag is about an order of magnitude larger than the linear corotation torque in the zero-softening limit. Thus, the horseshoe drag, the sign of which depends on the gradient of specific vorticity, is important for estimates of the total torque acting on the planet. We further show that phenomena, such as the Lindblad wakes, with a radial separation from corotation of approximately a pressure scaleheight H can affect x _s, even though for low-mass planets   x _s≪ H   . The effect is to distort streamlines and reduce x _s through the action of a back pressure. This effect is reduced for smaller gravitational softening parameters and planets of higher mass, for which x _s becomes comparable to H .     

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 dynamics of two massive planets on inclined orbits

The significant orbital eccentricities of most giant extrasolar planets may have their origin in the gravitational dynamics of initially unstable multiple planet systems. In this work, we explore the dynamics of two close planets on inclined orbits through both analytical techniques and extensive numerical scattering experiments. We derive a criterion for two equal mass planets on circular inclined orbits to achieve Hill stability, and conclude that significant radial migration and eccentricity pumping of both planets occurs predominantly by 2:1 and 5:3 mean motion resonant interactions. Using Laplace-Lagrange secular theory, we obtain analytical secular solutions for the orbital inclinations and longitudes of ascending nodes, and use those solutions to distinguish between the secular and resonant dynamics which arise in numerical simulations. We also illustrate how encounter maps, typically used to trace the motion of massless particles, may be modified to reproduce the gross instability seen by the numerical integrations. Such a correlation suggests promising future use of such maps to model the dynamics of more coplanar massive planet systems.     

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.     

Effects of mass loss for highly-irradiated giant planets

We present calculations for the evolution and surviving mass of highly-irradiated extrasolar giant planets (EGPs) at orbital semimajor axes ranging from 0.023 to 0.057 AU using a generalized scaled theory for mass loss, together with new surface-condition grids for hot EGPs and a consistent treatment of tidal truncation. Theoretical estimates for the rate of energy-limited hydrogen escape from giant-planet atmospheres differ by two orders of magnitude, when one holds planetary mass, composition, and irradiation constant. Baraffe et al. [Baraffe, I., Selsis, F., Chabrier, G., Barman, T.S., Allard, F., Hauschildt, P.H., Lammer, H., 2004. Astron. Astrophys. 419, L13-L16] predict the highest rate, based on the theory of Lammer et al. [Lammer, H., Selsis, F., Ribas, I., Guinan, E.F., Bauer, S.J., Weiss, W.W., 2003. Astrophys. J. 598, L121-L124]. Scaling the theory of Watson et al. [Watson, A.J., Donahue, T.M., Walker, J.C.G., 1981. Icarus 48, 150-166] to parameters for a highly-irradiated exoplanet, we find an escape rate ∼¹⁰² lower than Baraffe's. With the scaled Watson theory we find modest mass loss, occurring early in the history of a hot EGP. In this theory, mass loss including the effect of Roche-lobe overflow becomes significant primarily for masses below a Saturn mass, for semimajor axes ?0.023 AU. This contrasts with the Baraffe model, where hot EGPs are claimed to be remnants of much more massive bodies, originally several times Jupiter and still losing substantial mass fractions at present.     

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.     

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