Jupiter and Saturn formed in a few million years (ref. 1) from a gas-dominated protoplanetary disk, and were susceptible to gas-driven migration of their orbits on timescales of only ∼100,000 years (ref. 2). Hydrodynamic simulations show that these giant planets can undergo a two-stage, inward-then-outward, migration3,4,5. The terrestrial planets finished accreting much later6, and their characteristics, including Mars’ small mass, are best reproduced by starting from a planetesimal disk with an outer edge at about one astronomical unit from the Sun7,8 (1 AU is the Earth–Sun distance). Here we report simulations of the early Solar System that show how the inward migration of Jupiter to 1.5 AU, and its subsequent outward migration, lead to a planetesimal disk truncated at 1 AU; the terrestrial planets then form from this disk over the next 30–50 million years, with an Earth/Mars mass ratio consistent with observations. Scattering by Jupiter initially empties but then repopulates the asteroid belt, with inner-belt bodies originating between 1 and 3 AU and outer-belt bodies originating between and beyond the giant planets. This explains the significant compositional differences across the asteroid belt. The key aspect missing from previous models of terrestrial planet formation is the substantial radial migration of the giant planets, which suggests that their behaviour is more similar to that inferred for extrasolar planets than previously thought.
We investigate the evolution of grains composed of an ice shell surrounding an olivine core as they pass through a spiral shock in a protoplanetary disk. We use published three-dimensional radiation-hydrodynamics simulations of massive self-gravitating protoplanetary disks to extract the thermodynamics of spiral shocks in the region between 10 and 20 AU from the central star. As the density wave passes, it heats the grains, causing them to lose their ice shell and resulting in a lowering of the grain opacity. In addition, since grains of different sizes will have slightly different temperatures, there will be a migration of ice from hotter grains to cooler ones. The rate of migration depends on the temperature of the background gas, so a grain distribution that is effectively stable for low temperatures can undergo an irreversible change in opacity if the gas is temporarily heated to above ~150 K. We find that the opacity can drop more and at a significantly faster rate throughout the spiral shocks relative to the prediction of the standard dust grain model adopted in hydrodynamical calculations of protoplanetary disks. This would lead to faster gas cooling within spiral arms. We discuss the implications of our results on the susceptibility of disks to fragment into sub-stellar objects at distances of a few tens of astronomical units.
We present the results of a search for white-light flares on ~23,000 cool dwarfs in the Kepler Quarter 1 long cadence data. We have identified 373 flaring stars, some of which flare multiple times during the observation period. We calculate relative flare energies, flare rates, and durations and compare these with the quiescent photometric variability of our sample. We find that M dwarfs tend to flare more frequently but for shorter durations than K dwarfs and that they emit more energy relative to their quiescent luminosity in a given flare than K dwarfs. Stars that are more photometrically variable in quiescence tend to emit relatively more energy during flares, but variability is only weakly correlated with flare frequency. We estimate distances for our sample of flare stars and find that the flaring fraction agrees well with other observations of flare statistics for stars within 300 pc above the Galactic plane. These observations provide a more rounded view of stellar flares by sampling stars that have not been pre-selected by their activity, and are informative for understanding the influence of these flares on planetary habitability.
We apply two equilibrium tide theories – a constant-phase-lag model and a constant-time-lag model – to compute the obliquity evolution of terrestrial planets orbiting in the habitable zones around LMSs. The time for the obliquity to decrease from an Earth-like obliquity of 23.5° to 5°, the “tilt erosion time”, is compared to the traditional insolation habitable zone (IHZ) in the parameter space spanned by the semi-major axis a, the eccentricity e, and the stellar mass Ms. We also compute tidal heating and equilibrium rotation caused by obliquity tides as further constraints on habitability. The Super-Earth Gl581 d and the planet candidate Gl581 g are studied as examples for these tidal processes.
We evaluate the orbital evolution and several plausible origin scenarios for the mutually inclined orbits of υ And c and d. These two planets have orbital elements that oscillate with large amplitudes and lie close to the stability boundary. This configuration, and in particular the observed mutual inclination, demands an explanation. The planetary system may be influenced by a nearby low-mass star, υ And B, which could perturb the planetary orbits, but we find it cannot modify two coplanar orbits into the observed mutual inclination of 30°. However, it could incite ejections or collisions between planetary companions that subsequently raise the mutual inclination to >30°. Our simulated systems with large mutual inclinations tend to be further from the stability boundary than υ And, but we are able to produce similar systems. We conclude that scattering is a plausible mechanism to explain the observed orbits of υ And c and d, but we cannot determine whether the scattering was caused by instabilities among the planets themselves or by perturbations from υ And B. We also develop a procedure to quantitatively compare numerous properties of the observed system to our numerical models. Although we only implement this procedure to υ And, it may be applied to any exoplanetary system.
We provide an overview of stellar variability in the first quarter data from the Kepler mission. The intent of this paper is to examine the entire sample of over 150,000 target stars for periodic behavior in their light curves and relate this to stellar characteristics. This data set constitutes an unprecedented study of stellar variability given its great precision and complete time coverage (with a half hour cadence). Because the full Kepler pipeline is not currently suitable for a study of stellar variability of this sort, we describe our procedures for treating the “raw” pipeline data. About half of the total sample exhibits convincing periodic variability up to two weeks, with amplitudes ranging from differential intensity changes of less than 10–4 up to more than 10%. K and M dwarfs have a greater fraction of period behavior than G dwarfs. The giants in the sample have distinctive quasi-periodic behavior, but are not periodic in the way we define it. Not all periodicities are due to rotation, and the most significant period is not necessarily the rotation period. We discuss properties of the light curves, and in particular look at a sample of very clearly periodic G dwarfs. It is clear that a large number of them do vary because of rotation and starspots, but it will take further analysis to fully exploit this.
The migration of the giant planets due to the scattering of planetesimals causes powerful resonances to move through the asteroid belt and the terrestrial planet region. Exactly when and how the giant planets migrated is not well known. In this paper we present results of an investigation of the formation of the terrestrial planets during and after the migration of the giant planets. The latter is assumed to have occurred immediately after the dissipation of the nebular disk – i.e. “early” with respect to the timing of the late heavy bombardment (LHB). The presumed cause of our modeled early migration of the giant planets is angular mometum transfer between the planets and scattered planetesimals.
Circulations and habitable zones of planets orbiting low-mass stars are investigated. Many of these planets are expected to rotate synchronously relative to their parent stars, thereby raising questions about their surface temperature distributions and habitability. We use a global circulation model to study idealized, synchronously rotating (tidally locked) planets of various rotation periods, with surfaces of all land or all water, but with an Earth-like atmosphere and solar insolation. The dry planets exhibit wide variations in surface temperature: >80 °C on the dayside to <-110 °C on the nightside for the 240-h rotator, for example. The water-covered aquaplanets are warmer and exhibit narrower ranges of surface temperatures, e.g., ∼40 °C to >-60 °C for the 240-h orbiter. They also have a larger habitable area, defined here as the region where average surface temperatures are between 0 °C and 50 °C. This concept has little relevance for either dry or aquaplanets, but might become relevant on a planet with both land area and oceans.
Main sequence M stars pose an interesting problem for astrobiology: their abundance in our galaxy makes them likely targets in the hunt for habitable planets, but their strong chromospheric activity produces high-energy radiation and charged particles that may be detrimental to life. We studied the impact of the 1985 April 12 flare from the M dwarf AD Leonis (AD Leo), simulating the effects from both UV radiation and protons on the atmospheric chemistry of a hypothetical, Earth-like planet located within its habitable zone. Based on observations of solar proton events and the Neupert effect, we estimated a proton flux associated with the flare of 5.9 × 108 protons cm−2 sr−1 s−1 for particles with energies >10 MeV. Then we calculated the abundance of nitrogen oxides produced by the flare by scaling the production of these compounds during a large solar proton event called the Carrington event. The simulations were performed with a 1-D photochemical model coupled to a 1-D radiative/convective model. Our results indicate that the UV radiation emitted during the flare does not produce a significant change in the ozone column depth of the planet. When the action of protons is included, the ozone depletion reaches a maximum of 94% two years after the flare for a planet with no magnetic field. At the peak of the flare, the calculated UV fluxes that reach the surface, in the wavelength ranges that are damaging for life, exceed those received on Earth during less than 100 s. Therefore, flares may not present a direct hazard for life on the surface of an orbiting habitable planet. Given that AD Leo is one of the most magnetically active M dwarfs known, this conclusion should apply to planets around other M dwarfs with lower levels of chromospheric activity. Key Words: M dwarf—Flare—Habitable zone—Planetary atmospheres. Astrobiology 10, 751–771.
Interstellar material is highly processed when subjected to the physical conditions that prevail in the inner regions of protoplanetary disks, the potential birthplace of habitable planets. Polycyclic aromatic hydrocarbons (PAHs) are abundant in the interstellar medium, and they have also been observed in the disks around young stars, with evidence for some modification in the latter. Using a chemical model developed for sooting flames, we have investigated the chemical evolution of PAHs in warm (1000–2000 K) and oxygen-rich (C/O < 1) conditions appropriate for the region where habitable planets may eventually form. Our study focuses on (1) delineating the conditions under which PAHs will react and (2) identifying the key reaction pathways and reaction products characterizing this chemical evolution. We find that reactions with H, OH and O are the main pathways for destroying PAHs over disk timescale at temperatures greater than about 1000 K. In the process, high abundances of C2H2 persist over long timescales due to the kinetic inhibition of reactions that eventually drive the carbon into CO, CO2 and CH4. The thermal destruction of PAHs may thus be the cause of the abundant C2H2 that has been observed in disks. We propose that protoplanetary disks have a ‘soot line’, within which PAHs are irreversibly destroyed via thermally-driven reactions. The soot line will play an important role, analogous to that of the ‘snow line’, in the bulk carbon content of meteorites and habitable planets.
We study the dynamical stability of planetary systems consisting of one hypothetical terrestrial-mass planet (1 or 10 M ⊕) and one massive planet (10 M ⊕-10 M jup). We consider masses and orbits that cover the range of observed planetary system architectures (including non-zero initial eccentricities), determine the stability limit through N-body simulations, and compare it to the analytic Hill stability boundary. We show that for given masses and orbits of a two-planet system, a single parameter, which can be calculated analytically, describes the Lagrange stability boundary (no ejections or exchanges) but diverges significantly from the Hill stability boundary. However, we do find that the actual boundary is fractal, and therefore we also identify a second parameter which demarcates the transition from stable to unstable evolution. We show the portions of the habitable zones (HZs) of ρ CrB, HD 164922, GJ 674, and HD 7924 that can support a terrestrial planet. These analyses clarify the stability boundaries in exoplanetary systems and demonstrate that, for most exoplanetary systems, numerical simulations of the stability of potentially habitable planets are only necessary over a narrow region of the parameter space. Finally, we also identify and provide a catalog of known systems that can host terrestrial planets in their HZs.
The habitable zones (HZs) of main-sequence stars have traditionally been defined as the range of orbits that intercept the appropriate amount of stellar flux to permit surface water on a planet. Terrestrial exoplanets discovered to orbit M stars in these zones, which are close-in due to decreased stellar luminosity, may also undergo significant tidal heating. Tidal heating may span a wide range for terrestrial exoplanets and may significantly affect conditions near the surface. For example, if heating rates on an exoplanet are near or greater than that on Io (where tides drive volcanism that resurfaces the planet at least every 1 Myr) and produce similar surface conditions, then the development of life seems unlikely. On the other hand, if the tidal heating rate is less than the minimum to initiate plate tectonics, then CO2 may not be recycled through subduction, leading to a runaway greenhouse that sterilizes the planet. These two cases represent potential boundaries to habitability and are presented along with the range of the traditional HZ for main-sequence, low-mass stars. We propose a revised HZ that incorporates both stellar insolation and tidal heating. We apply these criteria to GJ 581 d and find that it is in the traditional HZ, but its tidal heating alone may be insufficient for plate tectonics.
To date, no accretion model has succeeded in reproducing all observed constraints in the inner Solar System. These constraints include: (1) the orbits, in particular the small eccentricities, and (2) the masses of the terrestrial planets – Mars’ relatively small mass in particular has not been adequately reproduced in previous simulations; (3) the formation timescales of Earth and Mars, as interpreted from Hf/W isotopes; (4) the bulk structure of the asteroid belt, in particular the lack of an imprint of planetary embryo-sized objects; and (5) Earth’s relatively large water content, assuming that it was delivered in the form of water-rich primitive asteroidal material. Here we present results of 40 high-resolution (N = 1000–2000) dynamical simulations of late-stage planetary accretion with the goal of reproducing these constraints, although neglecting the planet Mercury.
We present N-body simulations of planetary accretion beginning with 1 km radius planetesimals in orbit about a 1 M⊙ star at 0.4 AU. The initial disk of planetesimals contains too many bodies for any current N-body code to integrate; therefore, we model a sample patch of the disk. Although this greatly reduces the number of bodies, we still track in excess of 105 particles. We consider three initial velocity distributions and monitor the growth of the planetesimals. The masses of some particles increase by more than a factor of 100. Additionally, the escape speed of the largest particle grows considerably faster than the velocity dispersion of the particles, suggesting impending runaway growth, although no particle grows large enough to detach itself from the power law size-frequency distribution. These results are in general agreement with previous statistical and analytical results.
The known extrasolar multiple-planet systems share a surprising dynamical attribute: they cluster just beyond the Hill stability boundary. Here we show that the planet-planet scattering model, which naturally explains the observed exoplanet eccentricity distribution, can reproduce the observed distribution of dynamical configurations. We calculated how each of our scattered systems would appear over an appropriate range of viewing geometries; as Hill stability is weakly dependent on the masses, the mass-inclination degeneracy does not significantly affect our results. We consider a wide range of initial planetary mass distributions and find that some are poor fits to the observed systems. In fact, many of our scattering experiments overproduce systems very close to the stability boundary. The distribution of dynamical configurations of two-planet systems may provide better discrimination between scattering models than the distribution of eccentricity. Our results imply that, at least in their inner regions which are weakly affected by gas or planetesimal disks, planetary systems should be “packed,” with no large gaps between planets.
Three planets with minimum masses less than 10 M⊕ orbit the star HD 40307, suggesting these planets may be rocky. However, with only radial velocity data, it is impossible to determine if these planets are rocky or gaseous. Here we exploit various dynamical features of the system in order to assess the physical properties of the planets. Observations allow for circular orbits, but a numerical integration shows that the eccentricities must be at least 10–4. Also, planets b and c are so close to the star that tidal effects are significant. If planet b has tidal parameters similar to the terrestrial planets in the solar system and a remnant eccentricity larger than 10–3, then, going back in time, the system would have been unstable within the lifetime of the star (which we estimate to be 6.1 ± 1.6 Gyr). Moreover, if the eccentricities are that large and the inner planet is rocky, then its tidal heating may be an order of magnitude greater than extremely volcanic Io, on a per unit surface area basis. If planet b is not terrestrial, e.g., Neptune-like, these physical constraints would not apply. This analysis suggests the planets are not terrestrial-like, and are more like our giant planets. In either case, we find that the planets probably formed at larger radii and migrated early-on (via disk interactions) into their current orbits. This study demonstrates how the orbital and dynamical properties of exoplanet systems may be used to constrain the planets’ physical properties.
We study the dynamical stability of an additional, potentially habitable planet in the HD 47186 planetary system. Two planets are currently known in this system: a “hot Neptune” with a period of 4.08 days and a Saturn-mass planet with a period of 3.7 years. Here we consider the possibility that one or more undetected planets exist between the two known planets and possibly within the habitable zone (HZ) in this system. Given the relatively low masses of the known planets, additional planets could have masses $ \rlap{<}{\lower1.0ex\hbox{$\sim $}}10 {\,M_\oplus }$, and hence be terrestrial-like and further improving potential habitability. We perform N-body simulations to identify the stable zone between planets b and c and find that much of the inner HZ can harbor a 10 M ⊕ planet. With the current radial velocity threshold of ~1 m s–1, an additional planet should be detectable if it lies at the inner edge of the habitable zone at 0.8 AU. We also show that the stable zone could contain two additional planets of 10 M ⊕ each if their eccentricities are lower than ~0.3.
The tidal heating of hypothetical rocky (or terrestrial) extrasolar planets spans a wide range of values depending on stellar masses and initial orbits. Tidal heating may be sufficiently large (in many cases, in excess of radiogenic heating) and long-lived to drive plate tectonics, similar to the Earth’s, which may enhance the planet’s habitability. In other cases, excessive tidal heating may result in Io-like planets with violent volcanism, probably rendering them unsuitable for life. On water-rich planets, tidal heating may generate subsurface oceans analogous to Europa’s with similar prospects for habitability. Tidal heating may enhance the outgassing of volatiles, contributing to the formation and replenishment of a planet’s atmosphere. To address these issues, we model the tidal heating and evolution of hypothetical extrasolar terrestrial planets. The results presented here constrain the orbital and physical properties required for planets to be habitable.
Tides raised on a planet by the gravity of its host star can reduce the planet’s orbital semi-major axis and eccentricity. This effect is only relevant for planets orbiting very close to their host stars. The habitable zones of low-mass stars are also close in, and tides can alter the orbits of planets in these locations. We calculate the tidal evolution of hypothetical terrestrial planets around low-mass stars and show that tides can evolve planets past the inner edge of the habitable zone, sometimes in less than 1 billion years. This migration requires large eccentricities (>0.5) and low-mass stars (≲0.35 M⊙). Such migration may have important implications for the evolution of the atmosphere, internal heating, and the Gaia hypothesis. Similarly, a planet that is detected interior to the habitable zone could have been habitable in the past. We consider the past habitability of the recently discovered, ∼5 M⊕ planet, Gliese 581 c. We find that it could have been habitable for reasonable choices of orbital and physical properties as recently as 2 Gyr ago. However, when constraints derived from the additional companions are included, most parameter choices that indicate past habitability require the two inner planets of the system to have crossed their mutual 3:1 mean motion resonance. As this crossing would likely have resulted in resonance capture, which is not observed, we conclude that Gl 581 c was probably never habitable. Astrobiology 8, 557–568.
To date, two planetary systems have been discovered with close-in, terrestrial-mass planets forumla. Many more such discoveries are anticipated in the coming years with radial velocity and transit searches. Here we investigate the different mechanisms that could form ‘hot Earths’ and their observable predictions. Models include: (1) in situ accretion; (2) formation at larger orbital distance followed by inward ‘type 1’ migration; (3) formation from material being ‘shepherded’ inward by a migrating gas giant planet; (4) formation from material being shepherded by moving secular resonances during dispersal of the protoplanetary disc; (5) tidal circularization of eccentric terrestrial planets with close-in perihelion distances and (6) photoevaporative mass-loss of a close-in giant planet. Models 1–4 have been validated in previous work. We show that tidal circularization can form hot Earths, but only for relatively massive planets forumla with very close-in perihelion distances (≲0.025 au), and even then the net inward movement in orbital distance is at most only 0.1–0.15 au. For planets of less than forumla, photoevaporation can remove the planet’s envelope and leave behind the solid core on a Gyr time-scale, but only for planets inside 0.025–0.05 au. Using two quantities that are observable by current and upcoming missions, we show that these models each produce unique signatures, and can be observationally distinguished. These observables are the planetary system architecture (detectable with radial velocities, transits and transit timing) and the bulk composition of transiting close-in terrestrial planets (measured by transits via the planet’s radius).
We report the results of our Hubble Space Telescope (HST) snapshot survey with the ACS HRC PR200L prism, designed to measure the near-UV emission in a sample of nearby M dwarfs. Thirty-three stars were observed, spanning the mass range from 0.1 to 0.6 solar masses (Teff ~ 2200-4000 K) where the UV energy distributions vary widely between active and inactive stars. These observations provide much needed constraints on models of the habitability zone and the atmospheres of possible terrestrial planets orbiting M dwarf hosts and will be useful in refining the target selection for future space missions such as Terrestrial Planet Finder (TPF). We compare our data with a new generation of M dwarf atmospheric models and discuss their implications for the chromospheric energy budget. These NUV data will also be valuable in conjunction with existing optical, FUV, and X-ray data to explore unanswered questions regarding the dynamo generation and magnetic heating in low-mass stars.
Previous research has indicated that high amounts of ozone (O3) and oxygen (O2) may be produced abiotically in atmospheres with high concentrations of CO2. The abiotic production of these two gases, which are also characteristic of photosynthetic life processes, could pose a potential “false-positive” for remote-sensing detection of life on planets around other stars.We show here that such false positives are unlikely on any planet that possesses abundant liquid water, as rainout of oxidized species onto a reduced planetary surface should ensure that atmospheric H2 concentrations remain relatively high, and that O2 and O3 remain low. Our aim is to determine the amount of O3 and O2 formed in a high CO2 atmosphere for a habitable planet without life. We use a photochemical model that considers hydrogen (H2) escape and a detailed hydrogen balance to calculate the O2 and O3 formed on planets with 0.2 of CO2 around the Sun, and 0.02, 0.2 and 2 bars of CO2 around a young Sun-like star with higher UV radiation. The concentrations obtained by the photochemical model were used as input in a radiative transfer model that calculated the spectra of the modeled planets. The O3 and O2 concentrations in the simulated planets are extremely small, and unlikely to produce a detectable signature in the spectra of those planets. We conclude that with a balanced hydrogen budget, and for planets with an active hydrological cycle, abiotic formation of O2 and O3 is unlikely to create a possible false positive for life detection in either the visible/near-infrared or mid-infrared wavelength regimes.
The water content and habitability of terrestrial planets are determined during their final assembly, from perhaps 100 1,000-km “planetary embryos” and a swarm of billions of 1–10-km “planetesimals.” During this process, we assume that water-rich material is accreted by terrestrial planets via impacts of water-rich bodies that originate in the outer asteroid region. We present analysis of water delivery and planetary habitability in five high-resolution simulations containing about 10 times more particles than in previous simulations. These simulations formed 15 terrestrial planets from 0.4 to 2.6 Earth masses, including five planets in the habitable zone. Every planet from each simulation accreted at least the Earth’s current water budget; most accreted several times that amount (assuming no impact depletion). Each planet accreted at least five water-rich embryos and planetesimals from the past 2.5 astronomical units; most accreted 10–20 water-rich bodies. We present a new model for water delivery to terrestrial planets in dynamically calm systems, with low-eccentricity or low-mass giant planets—such systems may be very common in the Galaxy. We suggest that water is accreted in comparable amounts from a few planetary embryos in a “hit or miss” way and from millions of planetesimals in a statistically robust process. Variations in water content are likely to be caused by fluctuations in the number of water-rich embryos accreted, as well as from systematic effects, such as planetary mass and location, and giant planet properties. Key Words: Planetary formation—Water delivery—Extrasolar planets—Cosmochemistry. Astrobiology 7(1), 66–84.
Close-in giant planets (e.g., “hot Jupiters”) are thought to form far from their host stars and migrate inward, through the terrestrial planet zone, via torques with a massive gaseous disk. Here we simulate terrestrial planet growth during and after giant planet migration. Several-Earth-mass planets also form interior to the migrating jovian planet, analogous to recently discovered “hot Earths.” Very-water-rich, Earth-mass planets form from surviving material outside the giant planet’s orbit, often in the habitable zone and with low orbital eccentricities. More than a third of the known systems of giant planets may harbor Earth-like planets.
Coupled one-dimensional photochemical-climate calculations have been performed for hypothetical Earth-like planets around M dwarfs. Visible/near-infrared and thermal-infrared synthetic spectra of these planets were generated to determine which biosignature gases might be observed by a future, space-based telescope. Our star sample included two observed active M dwarfs—AD Leo and GJ 643—and three quiescent model stars. The spectral distribution of these stars in the ultraviolet generates a different photochemistry on these planets. As a result, the biogenic gases CH4, N2O, and CH3Cl have substantially longer lifetimes and higher mixing ratios than on Earth, making them potentially observable by space-based telescopes. On the active M-star planets, an ozone layer similar to Earth’s was developed that resulted in a spectroscopic signature comparable to the terrestrial one. The simultaneous detection of O2 (or O3) and a reduced gas in a planet’s atmosphere has been suggested as strong evidence for life. Planets circling M stars may be good locations to search for such evidence. Key Words: Habitable planets—M dwarfs—Biosignatures—Biomarkers—Terrestrial Planet Finder. Astrobiology 5, 706–725.