BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Sean N. Raymond

Sean N. Raymond appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.EPastro-ph.SRastro-phastro-ph.GAastro-ph.IMastro-ph.HEphysics.geo-phphysics.pop-ph

Catalog footprint

What is connected

79works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

79 published item(s)

preprint2025arXiv

Disk fraction among free-floating planetary-mass objects in Upper Scorpius

Free-floating planetary-mass objects (FFPs) have been detected through direct imaging in several young, nearby star-forming regions. The properties of circumstellar disks around these objects may provide a valuable probe into their origin but are currently limited by the small sample sizes explored. We aim to perform a statistical study of the occurrence of circumstellar disks down to the planetary-mass regime. We performed a systematic survey of disks among the population identified in the 5-10 Myr-old Upper Scorpius association (USC), restricted to members outside the younger, embedded Ophiuchus region and with estimated masses below 105 M_Jup. We took advantage of unWISE photometry to search for mid-infrared excesses in the WISE (W1-W2) color. We implemented a Bayesian outlier detection method, which models the photospheric sequence and computes excess probabilities for each object, enabling a statistically sound estimation of disk fractions. We explored disk fractions across an unprecedentedly fine mass grid, reaching down to objects as low as ~6 M_Jup assuming 5 Myr or ~8 M_Jup assuming 10 Myr, thus extending the previous lower boundary of disk fraction studies. Depending on the age, our sample includes between 17 and 40 FFPs. We confirm that the disk fraction steadily rises with decreasing mass and exceeds 30% near the substellar-to-planetary mass boundary at ~13 M_Jup. We find hints of a possible flattening in this trend around 25-45 M_Jup, potentially signaling a transition in the dominant formation processes. This shift in trend should be considered with caution and needs to be confirmed with more sensitive observations. Our results are consistent with the gradual dispersal of disks over time, as disk fractions in Upper Scorpius appear systematically lower than those in younger regions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Chemical Habitability: Supply and Retention of Life's Essential Elements During Planet Formation

Carbon, Hydrogen, Nitrogen, Oxygen, Phosphorus and Sulfur (CHNOPS) play key roles in the origin and proliferation of life on Earth. Given the universality of physics and chemistry, not least the ubiquity of water as a solvent and carbon as a backbone of complex molecules, CHNOPS are likely crucial to most habitable worlds. To help guide and inform the search for potentially habitable and ultimately inhabited environments, we begin by summarizing the CHNOPS budget of various reservoirs on Earth, their role in shaping our biosphere, and their origins in the Solar Nebula. We then synthesize our current understanding of how these elements behave and are distributed in diverse astrophysical settings, tracing their journeys from synthesis in dying stars to molecular clouds, protoplanetary settings, and ultimately temperate rocky planets around main sequence stars. We end by identifying key branching points during this journey, highlighting instances where a forming planets' distribution of CHNOPS can be altered dramatically, and speculating about the consequences for the chemical habitability of these worlds.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Early Solar System instability triggered by dispersal of the gaseous disk

The Solar System's orbital structure is thought to have been sculpted by an episode of dynamical instability among the giant planets. However, the instability trigger and timing have not been clearly established. Hydrodynamical modeling has shown that while the Sun's gaseous protoplanetary disk was present the giant planets migrated into a compact orbital configuration in a chain of resonances. Here we use dynamical simulations to show that the giant planets' instability was likely triggered by the dispersal of the gaseous disk. As the disk evaporated from the inside-out, its inner edge swept successively across and dynamically perturbed each planet's orbit in turn. The associated orbital shift caused a dynamical compression of the exterior part of the system, ultimately triggering instability. The final orbits of our simulated systems match those of the Solar System for a viable range of astrophysical parameters. The giant planet instability therefore took place as the gaseous disk dissipated, constrained by astronomical observations to be a few to ten million years after the birth of the Solar System. Terrestrial planet formation would not complete until after such an early giant planet instability; the growing terrestrial planets may even have been sculpted by its perturbations, explaining the small mass of Mars relative to Earth.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Mathematical encoding within multi-resonant planetary systems as SETI beacons

How might an advanced alien civilization manipulate the orbits within a planetary system to create a durable signpost that communicates its existence? While it is still debated whether such a purposeful advertisement would be prudent and wise, we propose that mean-motion resonances between neighboring planets -- with orbital periods that form integer ratios -- could in principle be used to encode simple sequences that one would not expect to form in nature. In this Letter we build four multi-resonant planetary systems and test their long-term orbital stability. The four systems each contain 6 or 7 planets and consist of: (i) consecutive integers from 1 to 6; (ii) prime numbers from 2 to 11; (iii) the Fibonacci sequence from 1 to 13; and (iv) the Lazy Caterer sequence from 1 to 16. We built each system using N-body simulations with artificial migration forces. We evaluated the stability of each system over the full 10 Gyr integration of the Sun's main sequence phase. We then tested the stability of these systems for an additional 10 Gyr, during and after post-main sequence evolution of the central stars (assumed to be Sun-like) to their final, white dwarf phase. The only system that was destabilized was the consecutive integer sequence (system i). The other three sequences therefore represent potential SETI beacons.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Thermal processing of Jupiter Family Comets during their chaotic orbital evolution

Evidence for cometary activity beyond Jupiter and Saturn's orbits -- such as that observed for Centaurs and long period comets -- suggests that the thermal processing of comet nuclei starts long before they enter the inner Solar System, where they are typically observed and monitored. Such observations raise questions as to the depth of unprocessed material, and whether the activity of JFCs can be representative of any primitive material. Here we model the coupled thermal and dynamical evolution of Jupiter Family Comets (JFCs), from the moment they leave their outer Solar System reservoirs until their ejection into interstellar space. We apply a thermal evolution model to a sample of simulated JFCs obtained from dynamical simulations (arXiv:1706.07447) that successfully reproduce the orbital distribution of observed JFCs. We show that due to the stochastic nature of comet trajectories toward the inner solar system, all simulated JFCs undergo multiple heating episodes resulting in significant modifications of their initial volatile contents. A statistical analysis constrains the extent of such processing. We suggest that primordial condensed hypervolatile ices should be entirely lost from the layers that contribute to cometary activity observed today. Our results demonstrate that understanding the orbital (and thus, heating) history of JFCs is essential when putting observations in a broader context.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Origin and dynamical evolution of the asteroid belt

The asteroid belt was dynamically shaped during and after planet formation. Despite representing a broad ring of stable orbits, the belt contains less than one one-thousandth of an Earth mass. The asteroid orbits are dynamically excited with a wide range in eccentricity and inclination and their compositions are diverse, with a general trend toward dry objects in the inner belt and more water-rich objects in the outer belt. Here we review models of the asteroid belt's origins and dynamical history. The classical view is that the belt was born with several Earth masses in planetesimals, then strongly depleted. However, it is possible that very few planetesimals ever formed in the asteroid region and that the belt's story is one of implantation rather than depletion. A number of processes may have implanted asteroids from different regions of the Solar System, dynamically removed them, and excited their orbits. During the gaseous disk phase these include the effects of giant planet growth and migration and sweeping secular resonances. After the gaseous disk phase these include scattering from resident planetary embryos, chaos in the giant planets' orbits, the giant planet instability, and long-term dynamical evolution. Different global models for Solar System formation imply contrasting dynamical histories of the asteroid belt. Vesta and Ceres may have been implanted from opposite regions of the Solar System -- Ceres from the Jupiter-Saturn region and Vesta from the terrestrial planet region -- and could therefore represent very different formation conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Planetesimal rings as the cause of the Solar System's planetary architecture

Astronomical observations reveal that protoplanetary disks around young stars commonly have ring- and gap-like structures in their dust distributions. These features are associated with pressure bumps trapping dust particles at specific locations, which simulations show are ideal sites for planetesimal formation. Here we show that our Solar System may have formed from rings of planetesimals -- created by pressure bumps -- rather than a continuous disk. We model the gaseous disk phase assuming the existence of pressure bumps near the silicate sublimation line (at $T \sim$1400~K), water snowline (at $T \sim$170~K), and CO-snowline (at $T \sim$30~K). Our simulations show that dust piles up at the bumps and forms up to three rings of planetesimals: a narrow ring near 1~au, a wide ring between $\sim$3-4~au and $\sim$10-20~au, and a distant ring between $\sim$20~au and $\sim$45~au. We use a series of simulations to follow the evolution of the innermost ring and show how it can explain the orbital structure of the inner Solar System and provides a framework to explain the origins of isotopic signatures of Earth, Mars and different classes of meteorites. The central ring contains enough mass to explain the rapid growth of the giant planets' cores. The outermost ring is consistent with dynamical models of Solar System evolution proposing that the early Solar System had a primordial planetesimal disk beyond the current orbit of Uranus.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Planet formation: key mechanisms and global models

Models of planet formation are built on underlying physical processes. In order to make sense of the origin of the planets we must first understand the origin of their building blocks. This review comes in two parts. The first part presents a detailed description of six key mechanisms of planet formation: 1) The structure and evolution of protoplanetary disks 2) The formation of planetesimals 3) Accretion of protoplanets 4) Orbital migration of growing planets 5) Gas accretion and giant planet migration 6) Resonance trapping during planet migration. While this is not a comprehensive list, it includes processes for which our understanding has changed in recent years or for which key uncertainties remain. The second part of this review shows how global models are built out of planet formation processes. We present global models to explain different populations of known planetary systems, including close-in small/low-mass planets (i.e., super-Earths), giant exoplanets, and the Solar System's planets. We discuss the different sources of water on rocky exoplanets, and use cosmochemical measurements to constrain the origin of Earth's water. We point out the successes and failings of different models and how they may be falsified. Finally, we lay out a path for the future trajectory of planet formation studies.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The First Habitable Zone Earth-sized Planet from TESS. I: Validation of the TOI-700 System

We present the discovery and validation of a three-planet system orbiting the nearby (31.1 pc) M2 dwarf star TOI-700 (TIC 150428135). TOI-700 lies in the TESS continuous viewing zone in the Southern Ecliptic Hemisphere; observations spanning 11 sectors reveal three planets with radii ranging from 1 R$_\oplus$ to 2.6 R$_\oplus$ and orbital periods ranging from 9.98 to 37.43 days. Ground-based follow-up combined with diagnostic vetting and validation tests enable us to rule out common astrophysical false-positive scenarios and validate the system of planets. The outermost planet, TOI-700 d, has a radius of $1.19\pm0.11$ R$_\oplus$ and resides in the conservative habitable zone of its host star, where it receives a flux from its star that is approximately 86% of the Earth's insolation. In contrast to some other low-mass stars that host Earth-sized planets in their habitable zones, TOI-700 exhibits low levels of stellar activity, presenting a valuable opportunity to study potentially-rocky planets over a wide range of conditions affecting atmospheric escape. While atmospheric characterization of TOI-700 d with the James Webb Space Telescope (JWST) will be challenging, the larger sub-Neptune, TOI-700 c (R = 2.63 R$_\oplus$), will be an excellent target for JWST and beyond. TESS is scheduled to return to the Southern Hemisphere and observe TOI-700 for an additional 11 sectors in its extended mission, which should provide further constraints on the known planet parameters and searches for additional planets and transit timing variations in the system.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The nature and origins of sub-Neptune size planets

Planets intermediate in size between the Earth and Neptune, and orbiting closer to their host stars than Mercury does the Sun, are the most common type of planet revealed by exoplanet surveys over the last quarter century. Results from NASA's Kepler mission have revealed a bimodality in the radius distribution of these objects, with a relative underabundance of planets between 1.5 and 2.0 $R_{\oplus}$. This bimodality suggests that sub-Neptunes are mostly rocky planets that were born with primary atmospheres a few percent by mass accreted from the protoplanetary nebula. Planets above the radius gap were able to retain their atmospheres ("gas-rich super-Earths"), while planets below the radius gap lost their atmospheres and are stripped cores ("true super-Earths"). The mechanism that drives atmospheric loss for these planets remains an outstanding question, with photoevaporation and core-powered mass loss being the prime candidates. As with the mass-loss mechanism, there are two contenders for the origins of the solids in sub-Neptune planets: the migration model involves the growth and migration of embryos from beyond the ice line, while the drift model involves inward-drifting pebbles that coagulate to form planets close-in. Atmospheric studies have the potential to break degeneracies in interior structure models and place additional constraints on the origins of these planets. However, most atmospheric characterization efforts have been confounded by aerosols. Observations with upcoming facilities are expected to finally reveal the atmospheric compositions of these worlds, which are arguably the first fundamentally new type of planetary object identified from the study of exoplanets.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The origins of nearly coplanar, non-resonant systems of close-in super-Earths

Some systems of close-in "super-Earths" contain five or more planets on non-resonant but compact and nearly coplanar orbits. The Kepler-11 system is an iconic representative of this class of system. It is challenging to explain their origins given that planet-disk interactions are thought to be essential to maintain such a high degree of coplanarity, yet these same interactions invariably cause planets to migrate into chains of mean motion resonances. Here we mine a large dataset of dynamical simulations of super-Earth formation by migration. These simulations match the observed period ratio distribution as long as the vast majority of planet pairs in resonance become dynamically unstable. When instabilities take place resonances are broken during a late phase of giant impacts, and typical surviving systems have planet pairs with significant mutual orbital inclinations. However, a subset of our unstable simulations matches the Kepler-11 system in terms of coplanarity, compactness, planet-multiplicity and non-resonant state. This subset have dynamical instability phases typically much shorter than ordinary systems. Unstable systems may keep a high degree of coplanarity post-instability if planets collide at very low orbital inclinations ($\lesssim1^\circ$) or if collisions promote efficient damping of orbital inclinations. If planetary scattering during the instability takes place at low orbital inclinations ($\text{i}\lesssim1^\circ$), orbital inclinations are barely increased by encounters before planets collide.When planetary scattering pumps orbital inclinations to higher values ($\gtrsim 1^\circ$) planets tend to collide at higher mutual orbital inclinations, but depending on the geometry of collisions mergers' orbital inclinations may be efficiently damped. Each of these formation pathways can produce analogues to the Kepler-11 system.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A record of the final phase of giant planet migration fossilized in the asteroid belt's orbital structure

The asteroid belt is characterized by an extreme low total mass of material on dynamically excited orbits. The Nice Model explains many peculiar qualities of the solar system, including the belt's excited state, by invoking an orbital instability between the outer planets. However, previous studies of the Nice Model's effect on the belt's structure struggle to reproduce the innermost asteroids' orbital inclination distribution. Here, we show how the final phase of giant planet migration sculpts the asteroid belt, in particular its inclination distribution. As interactions with leftover planetesimals cause Saturn to move away from Jupiter, its rate of orbital precession slows as the two planets' mutual interactions weaken. When the planets approach their modern separation, where Jupiter completes just short of five orbits for every two of Saturn's, Jupiter's eccentric forcing on Saturn strengthens. We use numerical simulations to show that the absence of asteroids with orbits that precess between 24-28 arcsec/yr is related to the inclination problem. As Saturn's precession speeds back up, high inclination asteroids are excited on to planet crossing orbits and removed from the inner main belt. Through this process, the asteroid belt's orbital structure is reshaped, leading to markedly improved simulation outcomes.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Solar System Formation in the Context of Extra-Solar Planets

Exoplanet surveys have confirmed one of humanity's (and all teenagers') worst fears: we are weird. If our Solar System were observed with present-day Earth technology -- to put our system and exoplanets on the same footing -- Jupiter is the only planet that would be detectable. The statistics of exo-Jupiters indicate that the Solar System is unusual at the ~1% level among Sun-like stars (or ~0.1% among all stars). But why are we different? Successful formation models for both the Solar System and exoplanet systems rely on two key processes: orbital migration and dynamical instability. Systems of close-in super-Earths or sub-Neptunes require substantial radial inward motion of solids either as drifting mm- to cm-sized pebbles or migrating Earth-mass or larger planetary embryos. We argue that, regardless of their formation mode, the late evolution of super-Earth systems involves migration into chains of mean motion resonances, generally followed by instability when the disk dissipates. This pattern is likely also ubiquitous in giant planet systems. We present three models for inner Solar System formation -- the low-mass asteroid belt, Grand Tack, and Early Instability models -- each invoking a combination of migration and instability. We identify bifurcation points in planetary system formation. We present a series of events to explain why our Solar System is so weird. Jupiter's core must have formed fast enough to quench the growth of Earth's building blocks by blocking the flux of inward-drifting pebbles. The large Jupiter/Saturn mass ratio is rare among giant exoplanets but may be required to maintain Jupiter's wide orbit. The giant planets' instability must have been gentle, with no close encounters between Jupiter and Saturn, also unusual in the larger (exoplanet) context. Our Solar System system is thus the outcome of multiple unusual, but not unheard of, events.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Did Jupiter's core form in the innermost parts of the Sun's protoplanetary disk?

Jupiter's core is generally assumed to have formed beyond the snow line. Here we consider an alternative scenario, that Jupiter's core may have accumulated in the innermost parts of the protoplanetary disk. A growing body of research suggests that small particles ("pebbles") continually drift inward through the disk. If a fraction of drifting pebbles is trapped at the inner edge of the disk a several Earth-mass core can quickly grow. Subsequently, the core may migrate outward beyond the snow line via planet-disk interactions. Of course, to reach the outer Solar System Jupiter's core must traverse the terrestrial planet-forming region. We use N-body simulations including synthetic forces from an underlying gaseous disk to study how the outward migration of Jupiter's core sculpts the terrestrial zone. If the outward migration is fast (Tmig~10^4 years), the core simply migrates past resident planetesimals and planetary embryos. However, if its migration is slower (Tmig~10^5 years) the core removes solids from the inner disk by shepherding objects in mean motion resonances. In many cases the disk interior to 0.5-1 AU is cleared of embryos and most planetesimals. By generating a mass deficit close to the Sun, the outward migration of Jupiter's core may thus explain the absence of terrestrial planets closer than Mercury. Jupiter's migrating core often stimulates the growth of another large (~Earth-mass) core -- that may provide a seed for Saturn's core -- trapped in exterior resonance. The migrating core also may transport a fraction of terrestrial planetesimals, such as the putative parent bodies of iron meteorites, to the asteroid belt.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Is there an exoplanet in the Solar System?

We investigate the prospects for the capture of the proposed Planet 9 from other stars in the Sun's birth cluster. Any capture scenario must satisfy three conditions: the encounter must be more distant than ~150 au to avoid perturbing the Kuiper belt; the other star must have a wide-orbit planet (a>~100au); the planet must be captured onto an appropriate orbit to sculpt the orbital distribution of wide-orbit Solar System bodies. Here we use N-body simulations to show that these criteria may be simultaneously satisfied. In a few percent of slow close encounters in a cluster, bodies are captured onto heliocentric, Planet 9-like orbits. During the ~100 Myr cluster phase, many stars are likely to host planets on highly-eccentric orbits with apastron distances beyond 100 au if Neptune-sized planets are common and susceptible to planet--planet scattering. While the existence of Planet 9 remains unproven, we consider capture from one of the Sun's young brethren a plausible route to explain such an object's orbit. Capture appears to predict a large population of Trans-Neptunian Objects (TNOs) whose orbits are aligned with the captured planet, and we propose that different formation mechanisms will be distinguishable based on their imprint on the distribution of TNOs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Asteroid Belt as a Relic From a Chaotic Early Solar System

The orbital structure of the asteroid belt holds a record of the Solar System's dynamical history. The current belt only contains ${\rm \sim 10^{-3}}$ Earth masses yet the asteroids' orbits are dynamically excited, with a large spread in eccentricity and inclination. In the context of models of terrestrial planet formation, the belt may have been excited by Jupiter's orbital migration. The terrestrial planets can also be reproduced without invoking a migrating Jupiter; however, as it requires a severe mass deficit beyond Earth's orbit, this model systematically under-excites the asteroid belt. Here we show that the orbits of the asteroids may have been excited to their current state if Jupiter and Saturn's early orbits were chaotic. Stochastic variations in the gas giants' orbits cause resonances to continually jump across the main belt and excite the asteroids' orbits on a timescale of tens of millions of years. While hydrodynamical simulations show that the gas giants were likely in mean motion resonance at the end of the gaseous disk phase, small perturbations could have driven them into a chaotic but stable state. The gas giants' current orbits were achieved later, during an instability in the outer Solar System. Although it is well known that the present-day Solar System exhibits chaotic behavior, our results suggest that the early Solar System may also have been chaotic.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The habitability of Proxima Centauri b II. Possible climates and Observability

Radial velocity monitoring has found the signature of a $M \sin i = 1.3$~M$_\oplus$ planet located within the Habitable Zone (HZ) of Proxima Centauri \citep{Anglada16}. Despite a hotter past and an active host star the planet Proxima~b could have retained enough volatiles to sustain surface habitability \citep{Ribas2016}. Here we use a 3D Global Climate Model (GCM) to simulate Proxima b's atmosphere and water cycle for its two likely rotation modes (1:1 and 3:2 spin-orbit resonances) while varying the unconstrained surface water inventory and atmospheric greenhouse effect. We find that a broad range of atmospheric compositions allow surface liquid water. On a tidally-locked planet with sufficient surface water inventory, liquid water is always present, at least in the substellar region. With a non-synchronous rotation, this requires a minimum greenhouse warming ($\sim$10~mbar of CO$_2$ and 1~bar of N$_2$). If the planet is dryer, $\sim$0.5~bar/1.5~bars of CO$_2$ (respectively for asynchronous/synchronous rotation) suffice to prevent the trapping of any arbitrary small water inventory into polar/nightside ice caps. We produce reflection/emission spectra and phase curves for the simulated climates. We find that atmospheric characterization will be possible by direct imaging with forthcoming large telescopes. The angular separation of $7 λ/D$ at 1~$μ$m (with the E-ELT) and a contrast of $\sim$10$^{-7}$ will enable high-resolution spectroscopy and the search for molecular signatures, including H$_2$O, O$_2$, and CO$_2$. The observation of thermal phase curves can be attempted with JWST, thanks to a contrast of $2\times10^{-5}$ at 10~$μ$m. Proxima~b will also be an exceptional target for future IR interferometers. Within a decade it will be possible to image Proxima~b and possibly determine whether this exoplanet's surface is habitable.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The habitability of Proxima Centauri b. I. Irradiation, rotation and volatile inventory from formation to the present

Proxima b is a planet with a minimum mass of 1.3 MEarth orbiting within the habitable zone (HZ) of Proxima Centauri, a very low-mass, active star and the Sun's closest neighbor. Here we investigate a number of factors related to the potential habitability of Proxima b and its ability to maintain liquid water on its surface. We set the stage by estimating the current high-energy irradiance of the planet and show that the planet currently receives 30 times more EUV radiation than Earth and 250 times more X-rays. We compute the time evolution of the star's spectrum, which is essential for modeling the flux received over Proxima b's lifetime. We also show that Proxima b's obliquity is likely null and its spin is either synchronous or in a 3:2 spin-orbit resonance, depending on the planet's eccentricity and level of triaxiality. Next we consider the evolution of Proxima b's water inventory. We use our spectral energy distribution to compute the hydrogen loss from the planet with an improved energy-limited escape formalism. Despite the high level of stellar activity we find that Proxima b is likely to have lost less than an Earth ocean's worth of hydrogen before it reached the HZ 100-200 Myr after its formation. The largest uncertainty in our work is the initial water budget, which is not constrained by planet formation models. We conclude that Proxima b is a viable candidate habitable planet.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Water loss from Earth-sized planets in the habitable zones of ultracool dwarfs: Implications for the planets of TRAPPIST-1

Ultracool dwarfs (UCD; $T_{\rm eff}<\sim3000~$K) cool to settle on the main sequence after $\sim$1 Gyr. For brown dwarfs, this cooling never stops. Their habitable zone (HZ) thus sweeps inward at least during the first Gyr of their lives. Assuming they possess water, planets found in the HZ of UCDs have experienced a runaway greenhouse phase too hot for liquid water prior to entering the HZ. It has been proposed that such planets are desiccated by this hot early phase and enter the HZ as dry worlds. Here we model the water loss during this pre-HZ hot phase taking into account recent upper limits on the XUV emission of UCDs and using 1D radiation-hydrodynamic simulations. We address the whole range of UCDs but also focus on the planets recently found around the $0.08~M_\odot$ dwarf TRAPPIST-1. Despite assumptions maximizing the FUV-photolysis of water and the XUV-driven escape of hydrogen, we find that planets can retain significant amounts of water in the HZ of UCDs, with a sweet spot in the $0.04$-$0.06~M_\odot$ range. We also studied the TRAPPIST-1 system using observed constraints on the XUV-flux. We find that TRAPPIST-1b and c may have lost as much as 15 Earth Oceans and planet d -- which might be inside the HZ -- may have lost less than 1 Earth Ocean. Depending on their initial water contents, they could have enough water to remain habitable. TRAPPIST-1 planets are key targets for atmospheric characterization and could provide strong constraints on the water erosion around UCDs.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A primordial origin for the composition similarity between the Earth and the Moon

Most of the properties of the Earth-Moon system can be explained by a collision between a planetary embryo and the growing Earth late in the accretion process. Simulations show that most of the material that eventually aggregates to form the Moon originates from the impactor. However, analysis of the terrestrial and lunar isotopic composition show them to be highly similar. In contrast, the compositions of other solar system bodies are significantly different than the Earth and Moon. This poses a major challenge to the giant impact scenario since the Moon-forming impactor is then thought to also have differed in composition from the proto-Earth. Here we track the feeding zones of growing planets in a suite of simulations of planetary accretion, in order to measure the composition of Moon-forming impactors. We find that different planets formed in the same simulation have distinct compositions, but the compositions of giant impactors are systematically more similar to the planets they impact. A significant fraction of planet-impactor pairs have virtually identical compositions. Thus, the similarity in composition between the Earth and Moon could be a natural consequence of a late giant impact.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Accretion of Uranus and Neptune from inward-migrating planetary embryos blocked by Jupiter and Saturn

Reproducing Uranus and Neptune remains a challenge for simulations of solar system formation. The ice giants' peculiar obliquities suggest that they both suffered giant collisions during their formation. Thus, there must have been an epoch of accretion dominated by collisions among large planetary embryos in the primordial outer solar system. We test this idea using N-body numerical simulations including the effects of a gaseous protoplanetary disk. One strong constraint is that the masses of the ice giants are very similar -- the Neptune/Uranus mass ratio is $\sim1.18$. We show that similar-size ice giants do indeed form by collisions between planetary embryos beyond Saturn. The fraction of successful simulations varies depending on the initial number of planetary embryos in the system, their individual and total masses. Similar-sized ice giants are consistently reproduced in simulations starting with 5-10 planetary embryos with initial masses of $\sim$3-6 ${\rm M_\oplus}$. We conclude that accretion from a population of planetary embryos is a plausible scenario for the origin of Uranus and Neptune.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Gas giant planets as dynamical barriers to inward-migrating super-Earths

Planets of 1-4 times Earth's size on orbits shorter than 100 days exist around 30-50% of all Sun-like stars. In fact, the Solar System is particularly outstanding in its lack of "hot super-Earths" (or "mini-Neptunes"). These planets -- or their building blocks -- may have formed on wider orbits and migrated inward due to interactions with the gaseous protoplanetary disk. Here, we use a suite of dynamical simulations to show that gas giant planets act as barriers to the inward migration of super-Earths initially placed on more distant orbits. Jupiter's early formation may have prevented Uranus and Neptune (and perhaps Saturn's core) from becoming hot super-Earths. Our model predicts that the populations of hot super-Earth systems and Jupiter-like planets should be anti-correlated: gas giants (especially if they form early) should be rare in systems with many hot super-Earths. Testing this prediction will constitute a crucial assessment of the validity of the migration hypothesis for the origin of close-in super-Earths.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Highly Siderophile Elements in the Earth's Mantle as a Clock for the Moon-forming Impact

According to the generally accepted scenario, the last giant impact on the Earth formed the Moon and initiated the final phase of core formation by melting the Earth's mantle. A key goal of geochemistry is to date this event, but different ages have been proposed. Some argue for an early Moon-forming event, approximately 30 million years (Myr) after the condensation of the first solids in the Solar System, whereas others claim a date later than 50 Myr (and possibly as late as around 100 My) after condensation. Here we show that a Moon-forming event at 40 Myr after condensation, or earlier, is ruled out at a 99.9 per cent confidence level. We use a large number of N-body simulations to demonstrate a relationship between the time of the last giant impact on an Earth-like planet and the amount of mass subsequently added during the era known as Late Accretion. As the last giant impact is delayed, the late-accreted mass decreases in a predictable fashion. This relationship exists within both the classical scenario and the Grand Tack scenario of terrestrial planet formation, and it holds across a wide range of disk conditions. The concentration of highly siderophile elements (HSEs) in Earth's mantle constrains the mass of chondritic material added to Earth during Late Accretion. Using HSE abundance measurements, we determine a Moon-formation age of 95 +/- 32 Myr since condensation. The possibility exists that some late projectiles were differentiated and left an incomplete HSE record in Earth's mantle. Even in this case, various isotopic constraints strongly suggest that the late-accreted mass did not exceed 1 per cent of Earth's mass, and so the HSE clock still robustly limits the timing of the Moon-forming event to significantly later than 40 My after condensation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Long-lived Chaotic Orbital Evolution of Exoplanets in Mean Motion Resonances with Mutual Inclinations

We present N-body simulations of resonant planets with inclined orbits that show chaotically evolving eccentricities and inclinations that can persist for at least 10 Gyr. A wide range of behavior is possible, from fast, low amplitude variations to systems in which eccentricities reach 0.9999 and inclinations 179.9 degrees. While the orbital elements evolve chaotically, at least one resonant argument always librates. We show that the HD 73526, HD 45364 and HD 60532 systems may be in chaotically-evolving resonances. Chaotic evolution is apparent in the 2:1, 3:1 and 3:2 resonances, and for planetary masses from lunar- to Jupiter-mass. In some cases, orbital disruption occurs after several Gyr, implying the mechanism is not rigorously stable, just long-lived relative to the main sequence lifetimes of solar-type stars. Planet-planet scattering appears to yield planets in inclined resonances that evolve chaotically in about 0.5% of cases. These results suggest that 1) approximate methods for identifying unstable orbital architectures may have limited applicability, 2) the observed close-in exoplanets may be produced during the high eccentricity phases induced by inclined resonances, 3) those exoplanets' orbital planes may be misaligned with the host star's spin axis, 4) systems with resonances may be systematically younger than those without, 5) the distribution of period ratios of adjacent planets detected via transit may be skewed due to inclined resonances, and 6) potentially habitable planets in resonances may have dramatically different climatic evolution than the Earth. The GAIA spacecraft is capable of discovering giant planets in these types of orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Mercury-T: A new code to study tidally evolving multi-planet systems. Applications to Kepler-62

A large proportion of observed planetary systems contain several planets in a compact orbital configuration, and often harbor at least one close-in object. These systems are then most likely tidally evolving. We investigate how the effects of planet-planet interactions influence the tidal evolution of planets. We introduce for that purpose a new open-source addition to the Mercury N-body code, Mercury-T, which takes into account tides, general relativity and the effect of rotation-induced flattening in order to simulate the dynamical and tidal evolution of multi-planet systems. It uses a standard equilibrium tidal model, the constant time lag model. Besides, the evolution of the radius of several host bodies has been implemented (brown dwarfs, M-dwarfs of mass $0.1~M_\odot$, Sun-like stars, Jupiter). We validate the new code by comparing its output for one-planet systems to the secular equations results. We find that this code does respect the conservation of total angular momentum. We applied this new tool to the planetary system Kepler-62. We find that tides influence the stability of the system in some cases. We also show that while the four inner planets of the systems are likely to have slow rotation rates and small obliquities, the fifth planet could have a fast rotation rate and a high obliquity. This means that the two habitable zone planets of this system, Kepler-62e ad f are likely to have very different climate features, and this of course would influence their potential at hosting surface liquid water.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Terrestrial Planet Formation Constrained by Mars and the Structure of the Asteroid Belt

Reproducing the large Earth/Mars mass ratio requires a strong mass depletion in solids within the protoplanetary disk between 1 and 3 AU. The Grand Tack model invokes a specific migration history of the giant planets to remove most of the mass initially beyond 1 AU and to dynamically excite the asteroid belt. However, one could also invoke a steep density gradient created by inward drift and pile-up of small particles induced by gas-drag, as has been proposed to explain the formation of close-in super Earths. Here we show that the asteroid belt's orbital excitation provides a crucial constraint against this scenario for the Solar System. We performed a series of simulations of terrestrial planet formation and asteroid belt evolution starting from disks of planetesimals and planetary embryos with various radial density gradients and including Jupiter and Saturn on nearly circular and coplanar orbits. Disks with shallow density gradients reproduce the dynamical excitation of the asteroid belt by gravitational self-stirring but form Mars analogs significantly more massive than the real planet. In contrast, a disk with a surface density gradient proportional to $r^{-5.5}$ reproduces the Earth/Mars mass ratio but leaves the asteroid belt in a dynamical state that is far colder than the real belt. We conclude that no disk profile can simultaneously explain the structure of the terrestrial planets and asteroid belt. The asteroid belt must have been depleted and dynamically excited by a different mechanism such as, for instance, in the Grand Tack scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An Earth-sized Planet in the Habitable Zone of a Cool Star

The quest for Earth-like planets represents a major focus of current exoplanet research. While planets that are Earth-sized and smaller have been detected, these planets reside in orbits that are too close to their host star to allow liquid water on their surface. We present the detection of Kepler-186f, a 1.11+\-0.14 Earth radius planet that is the outermost of five planets - all roughly Earth-sized - that transit a 0.47+\-0.05 Rsun star. The intensity and spectrum of the star's radiation places Kepler-186f in the stellar habitable zone, implying that if Kepler-186f has an Earth-like atmosphere and H2O at its surface, then some of this H2O is likely to be in liquid form.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dynamics of exoplanetary systems, links to their habitability

Our knowledge of planets' orbital dynamics, which was based on Solar System studies, has been challenged by the diversity of exoplanetary systems. Around cool and ultra cool dwarfs, the influence of tides on the orbital and spin evolution of planets can strongly affect their climate and their capacity to host surface liquid water. We illustrate the role of tides and dynamics with the extreme case of planets orbiting around brown dwarfs. In multiple planet systems, the eccentricity is excited by planet-planet interactions. Planets are therefore heated up from the inside by the tidally-induced friction. This process can heat a habitable zone planet to such a level that surface liquid water cannot exist. We also talk about the newly discovered potentially habitable Earth-sized planet Kepler-186f. Given the poorly estimated age of the system, the planet could still be evolving towards synchronization and have a high obliquity or be pseudo-synchronized with a zero obliquity. These two configurations would have a different effect on the climate of this planet.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Formation, tidal evolution and habitability of the Kepler-186 system

The Kepler-186 system consists of five planets orbiting an early-M dwarf. The planets have physical radii of 1.0-1.50 R$_\oplus$ and orbital periods of 4 to 130 days. The $1.1~$R$_\oplus$ Kepler-186f with a period of 130 days is of particular interest. Its insolation of roughly $0.32~S_\odot$places it within the liquid water habitable zone. We present a multi-faceted study of the Kepler-186 system. First, we show that the distribution of planet masses can be roughly reproduced if the planets accreted from a high-surface density disk presumably sculpted by an earlier phase of migration. However, our simulations predict the existence of 1-2 undetected planets between planets e and f. Next, we present a dynamical analysis of the system including the effect of tides. The timescale for tidal evolution is short enough that the four inner planets must have small obliquities and near-synchronous rotation rates. Tidal evolution of Kepler-186f is slow enough that its current spin state depends on a combination of its dissipation rate and the stellar age. Finally, we study the habitability of Kepler-186f with a 1-D climate model. The planet's surface temperature can be raised above 273 K with 0.5-5 bars of CO$_2$, depending on the amount of N$_2$ present. Kepler-186f represents a case study of an Earth-sized planet in the cooler regions of the habitable zone of a cool star.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Hot super-Earths and giant planet cores from different migration histories

Planetary embryos embedded in gaseous protoplanetary disks undergo Type I orbital migration. Migration can be inward or outward depending on the local disk properties but, in general, only planets more massive than several $M_\oplus$ can migrate outward. Here we propose that an embryo's migration history determines whether it becomes a hot super-Earth or the core of a giant planet. Systems of hot super-Earths (or mini-Neptunes) form when embryos migrate inward and pile up at the inner edge of the disk. Giant planet cores form when inward-migrating embryos become massive enough to switch direction and migrate outward. We present simulations of this process using a modified N-body code, starting from a swarm of planetary embryos. Systems of hot super-Earths form in resonant chains with the innermost planet at or interior to the disk inner edge. Resonant chains are disrupted by late dynamical instabilities triggered by the dispersal of the gaseous disk. Giant planet cores migrate outward toward zero-torque zones, which move inward and eventually disappear as the disk disperses. Giant planet cores migrate inward with these zones and are stranded at ~1-5 AU. Our model reproduces several properties of the observed extra-solar planet populations. The frequency of giant planet cores increases strongly when the mass in solids is increased, consistent with the observed giant exoplanet - stellar metallicity correlation. The frequency of hot super-Earths is not a function of stellar metallicity, also in agreement with observations. Our simulations can reproduce the broad characteristics of the observed super-Earth population.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

No universal minimum-mass extrasolar nebula: Evidence against in-situ accretion of systems of hot super-Earths

It has been proposed that the observed systems of hot super-Earths formed in situ from high-mass disks. By fitting a disk profile to the entire population of Kepler planet candidates, Chiang & Laughlin (2013) constructed a "minimum-mass extrasolar nebula" with surface density profile Sigma r^-1.6. Here we use multiple-planet systems to show that it is inconsistent to assume a universal disk profile. Systems with 3-6 low-mass planets (or planet candidates) produce a diversity of minimum-mass disks with surface density profiles ranging from Sigma r^-3.2 to Sigma r^0.5 (5th-95th percentile). By simulating the transit detection of populations of synthetic planetary systems designed to match the properties of observed super-Earth systems, we show that a universal disk profile is statistically excluded at high confidence. Rather, the underlying distribution of minimum-mass disks is characterized by a broad range of surface density slopes. Models of gaseous disks can only explain a narrow range of slopes (roughly between r^0 and r^-1.5). Yet accretion of terrestrial planets in a gas-free environment preserves the initial radial distribution of building blocks. The known systems of hot super-Earths must therefore not represent the structure of their parent gas disks and can not have predominantly formed in situ. We instead interpret the diversity of disk slopes as the imprint of a process that re-arranged the solids relative to the gas in the inner parts of protoplanetary disks. A plausible mechanism is inward type 1 migration of Mars- to Earth-mass planetary embryos, perhaps followed by a final assembly phase.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Orbital Dynamics of Multi-Planet Systems with Eccentricity Diversity

Since exoplanets were detected using the radial velocity method, they have revealed a diverse distribution of orbital configurations. Amongst these are planets in highly eccentric orbits (e > 0.5). Most of these systems consist of a single planet but several have been found to also contain a longer period planet in a near-circular orbit. Here we use the latest Keplerian orbital solutions to investigate four known systems which exhibit this extreme eccentricity diversity; HD 37605, HD 74156, HD 163607, and HD 168443. We place limits on the presence of additional planets in these systems based on the radial velocity residuals. We show that the two known planets in each system exchange angular momentum through secular oscillations of their eccentricities. We calculate the amplitude and timescale for these eccentricity oscillations and associated periastron precession. We further demonstrate the effect of mutual orbital inclinations on the amplitude of high-frequency eccentricity oscillations. Finally, we discuss the implications of these oscillations in the context of possible origin scenarios for unequal eccentricities.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Planetesimal driven migration as an explanation for observations of high levels of warm, exozodiacal dust

High levels of exozodiacal dust have been observed in the inner regions of a large fraction of main sequence stars. Given the short lifetime of the observed small dust grains, these 'exozodis' are difficult to explain, especially for old (>100 Myr) stars. The exozodiacal dust may be observed as excess emission in the mid-infrared, or using interferometry. We hypothesise that exozodi are sustained by planetesimals scattered by planets inwards from an outer planetesimal belt, where collision timescales are long. In this work, we use N-body simulations to show that the outwards migration of a planet into a belt, driven by the scattering of planetesimals, can increase, or sustain, the rate at which planetesimals are scattered from the outer belt to the exozodi region. We hypothesise that this increase is sufficient to sustain the observed exozodi on Gyr timescales. No correlation between observations of an outer belt and an exozodi is required for this scenario to work, as the outer belt may be too faint to detect. If planetesimal driven migration does explain the observed exozodi, this work suggests that the presence of an exozodi indicates the presence of outer planets and a planetesimal belt.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Rossiter-McLaughlin Observations of 55 Cnc e

We present Rossiter-McLaughlin observations of the transiting super-Earth 55 Cnc e collected during six transit events between January 2012 and November 2013 with HARPS and HARPS-N. We detect no radial-velocity signal above 35 cm/s (3-sigma) and confine the stellar v sin i to 0.2 +/- 0.5 km/s. The star appears to be a very slow rotator, producing a very low amplitude Rossiter-McLaughlin effect. Given such a low amplitude, the Rossiter-McLaughlin effect of 55 Cnc e is undetected in our data, and any spin-orbit angle of the system remains possible. We also performed Doppler tomography and reach a similar conclusion. Our results offer a glimpse of the capacity of future instrumentation to study low amplitude Rossiter-McLaughlin effects produced by super-Earths.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Terrestrial Planet Formation at Home and Abroad

We review the state of the field of terrestrial planet formation with the goal of understanding the formation of the inner Solar System and low-mass exoplanets. We review the dynamics and timescales of accretion from planetesimals to planetary embryos and from embryos to terrestrial planets. We discuss radial mixing and water delivery, planetary spins and the importance of parameters regarding the disk and embryo properties. Next, we connect accretion models to exoplanets. We first explain why the observed hot Super Earths probably formed by in situ accretion or inward migration. We show how terrestrial planet formation is altered in systems with gas giants by the mechanisms of giant planet migration and dynamical instabilities. Standard models of terrestrial accretion fail to reproduce the inner Solar System. The "Grand Tack" model solves this problem using ideas first developed to explain the giant exoplanets. Finally, we discuss whether most terrestrial planet systems form in the same way as ours, and highlight the key ingredients missing in the current generation of simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Terrestrial Planet Formation in the Presence of Migrating Super-earths

Super-Earths with orbital periods less than 100 days are extremely abundant around Sun-like stars. It is unlikely that these planets formed at their current locations. Rather, they likely formed at large distances from the star and subsequently migrated inward. Here we use N-body simulations to study the effect of super-Earths on the accretion of rocky planets. In our simulations, one or more super-Earths migrates inward through a disk of planetary embryos and planetesimals embedded in a gaseous disk. We tested a wide range of migration speeds and configurations. Fast-migrating super-Earths ($τ_{mig} \sim$0.01-0.1 Myr) only have a modest effect on the protoplanetary embryos and planetesimals. Sufficient material survives to form rocky, Earth-like planets on orbits exterior to the super-Earths'. In contrast, slowly migrating super-Earths shepherd rocky material interior to their orbits and strongly deplete the terrestrial planet-forming zone. In this situation any Earth-sized planets in the habitable zone are extremely volatile-rich and are therefore probably not Earth-like.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Grand Tack model: a critical review

The `Grand Tack' model proposes that the inner Solar System was sculpted by the giant planets' orbital migration in the gaseous protoplanetary disk. Jupiter first migrated inward then Jupiter and Saturn migrated back outward together. If Jupiter's turnaround or "tack" point was at ~1.5 AU the inner disk of terrestrial building blocks would have been truncated at ~1 AU, naturally producing the terrestrial planets' masses and spacing. During the gas giants' migration the asteroid belt is severely depleted but repopulated by distinct planetesimal reservoirs that can be associated with the present-day S and C types. The giant planets' orbits are consistent with the later evolution of the outer Solar System. Here we confront common criticisms of the Grand Tack model. We show that some uncertainties remain regarding the Tack mechanism itself; the most critical unknown is the timing and rate of gas accretion onto Saturn and Jupiter. Current isotopic and compositional measurements of Solar System bodies -- including the D/H ratios of Saturn's satellites -- do not refute the model. We discuss how alternate models for the formation of the terrestrial planets each suffer from an internal inconsistency and/or place a strong and very specific requirement on the properties of the protoplanetary disk. We conclude that the Grand Tack model remains viable and consistent with our current understanding of planet formation. Nonetheless, we encourage additional tests of the Grand Tack as well as the construction of alternate models.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Tidal evolution in multiple planet systems: application to Kepler-62 and Kepler-186

A large number of observed exoplanets are part of multiple planet systems. Most of these systems are sufficiently close-in to be tidally evolving. In such systems, there is a competition between the excitation caused by planet-planet interactions and tidal damping. Using as an example two multiple planet systems, which host planets in the surface liquid water habitable zone (HZ): Kepler-62 and Kepler-186, we show the importance and effect of both planetary and stellar tides on the dynamical evolution of planets and on the climate of the HZ planets.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Vega's hot dust from icy planetesimals scattered inward by an outward-migrating planetary system

Vega has been shown to host multiple dust populations, including both hot exo-zodiacal dust at sub-AU radii and a cold debris disk extending beyond 100 AU. We use dynamical simulations to show how Vega's hot dust can be created by long-range gravitational scattering of planetesimals from its cold outer regions. Planetesimals are scattered progressively inward by a system of 5-7 planets from 30-60 AU to very close-in. In successful simulations the outermost planets are typically Neptune-mass. The back-reaction of planetesimal scattering causes these planets to migrate outward and continually interact with fresh planetesimals, replenishing the source of scattered bodies. The most favorable cases for producing Vega's exo-zodi have negative radial mass gradients, with sub-Saturn- to Jupiter-mass inner planets at 5-10 AU and outer planets of 2.5 to 20 Earth masses. The mechanism fails if a Jupiter-sized planet exists beyond ~15 AU because the planet preferentially ejects planetesimals before they can reach the inner system. Direct-imaging planet searches can therefore directly test this mechanism.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Very Wide Binary Stars as the Primary Source of Stellar Collisions in the Galaxy

We present numerical simulations modeling the orbital evolution of very wide binaries, pairs of stars separated by over ~1000 AU. Due to perturbations from other passing stars and the Milky Way's tide, the orbits of very wide binary stars occasionally become extremely eccentric, which forces close encounters between the companion stars (Kaib et al. 2013). We show that this process causes a stellar collision between very wide binary companion stars once every 1000-7500 years on average in the Milky Way. One of the main uncertainties in this collision rate is the amount of energy dissipated by dynamic tides during close (but not collisional) periastron passages. This dissipation presents a dynamical barrier to stellar collisions and can instead transform very wide binaries into close or contact binaries. However, for any plausible tidal dissipation model, very wide binary stars are an unrealized, and potentially the dominant, source of stellar collisions in our Galaxy. Such collisions should occur throughout the thin disk of the Milky Way. Stellar collisions within very wide binaries should yield a small population of single, Li-depleted, rapidly rotating massive stars.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Water Delivery and Giant Impacts in the 'Grand Tack' Scenario

A new model for terrestrial planet formation (Hansen 2009, Walsh et al. 2011) has explored accretion in a truncated protoplanetary disk, and found that such a configuration is able to reproduce the distribution of mass among the planets in the Solar System, especially the Earth/Mars mass ratio, which earlier simulations have generally not been able to match. Walsh et al. tested a possible mechanism to truncate the disk--a two-stage, inward-then-outward migration of Jupiter and Saturn, as found in numerous hydrodynamical simulations of giant planet formation. In addition to truncating the disk and producing a more realistic Earth/Mars mass ratio, the migration of the giant planets also populates the asteroid belt with two distinct populations of bodies--the inner belt is filled by bodies originating inside of 3 AU, and the outer belt is filled with bodies originating from between and beyond the giant planets (which are hereafter referred to as `primitive' bodies). We find here that the planets will accrete on order 1-2% of their total mass from primitive planetesimals scattered onto planet-crossing orbits during the formation of the planets. For an assumed value of 10% for the water mass fraction of the primitive planetesimals, this model delivers a total amount of water comparable to that estimated to be on the Earth today. While the radial distribution of the planetary masses and the dynamical excitation of their orbits are a good match to the observed system, we find that the last giant impact is typically earlier than 20 Myr, and a substantial amount of mass is accreted after that event. However, 5 of the 27 planets larger than half an Earth mass formed in all simulations do experience large late impacts and subsequent accretion consistent with the dating of the Moon-forming impact and the estimated amount of mass accreted by Earth following that event.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamical and collisional constraints on a stochastic late veneer on the terrestrial planets

Given their tendency to be incorporated into the core during differentiation, the highly-siderophile elements (HSEs) in Earth's mantle are thought to have been accreted as a `late veneer' after the end of the giant impact phase. Bottke et al (2010) proposed that the large Earth-to-Moon HSE abundance ratio can be explained if the late veneer was characterized by large (D = 1000-4000km) impactors. Here we simulate the evolution of the terrestrial planets during a stochastic late veneer phase from the end of accretion until the start of the late heavy bombardment ~500 Myr later. We show that a late veneer population of 0.05 Earth masses dominated by large (D > 1000km) bodies naturally delivers a ~0.01 Earth mass veneer to Earth, consistent with constraints. The eccentricities and inclinations of the terrestrial planets are excited by close encounters with the largest late veneer bodies. We find the best agreement with their post-veneer orbits if either a) the terrestrial planets' pre-veneer angular momentum deficit AMD_0 was less than half of the current one AMD_now, or b) AMD_0 <= AMD_now and the veneer was limited to smaller (D_max <= 2000km) bodies. Veneer impacts on Venus, Earth and Mars were mostly accretionary but on Mercury and the Moon they were mostly erosive. In ~20% of simulations an energetic impact occurred that could have removed >25% of Mercury's mass, thereby increasing its iron mass fraction. We show that, due to the erosive nature of larger impacts, the Moon cannot accrete any material from objects larger than 500-1000km. The large Earth-to-Moon HSE abundance ratio is explained if the late veneer included large impactors (D >= 500-1000km) regardless of their size distribution, as long as most of Earth's veneer came from large bodies. The spin angular momenta imparted by stochastic late veneer impacts was far in excess of the current ones for Mercury and Venus.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Habitable Planets Eclipsing Brown Dwarfs: Strategies for Detection and Characterization

Given the very close proximity of their habitable zones, brown dwarfs represent high-value targets in the search for nearby transiting habitable planets that may be suitable for follow-up occultation spectroscopy. In this paper we develop search strategies to find habitable planets transiting brown dwarfs depending on their maximum habitable orbital period (PHZ out). Habitable planets with PHZ out shorter than the useful duration of a night (e.g. 8-10 hrs) can be screened with 100 percent completeness from a single location and in a single night (near-IR). More luminous brown dwarfs require continuous monitoring for longer duration, e.g. from space or from a longitude-distributed network (one test scheduling achieved - 3 telescopes, 13.5 contiguous hours). Using a simulated survey of the 21 closest known brown dwarfs (within 7 pc) we find that the probability of detecting at least one transiting habitable planet is between 4.5 +5.6-1.4 and 56 +31-13 percent, depending on our assumptions. We calculate that brown dwarfs within 5-10 pc are characterizable for potential biosignatures with a 6.5 m space telescope using approx. 1 percent of a 5-year mission's life-time spread over a contiguous segment only one fifth to one tenth of this duration.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Making systems of Super Earths by inward migration of planetary embryos

Using N-body simulations with planet-disk interactions, we present a mechanism capable of forming compact systems of hot super Earths such as Kepler 11. Recent studies show that outward migration is common in the inner parts of radiative disks. However we show that two processes naturally tip the balance in favor of inward migration. First the corotation torque is too weak to generate outward migration for planetary embryos less massive than $4\mearth$. Second, system of multiple embryos generate sustained non-zero eccentricities that damp the corotation torque and again favor inward migration. Migration and accretion of planetary embryos in realistic disks naturally produce super Earths in resonant chains near the disk inner edge. Their compact configuration is similar to the observed systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Planetary System Disruption by Galactic Perturbations to Wide Binary Stars

Nearly half of the exoplanets found within binary star systems reside in very wide binaries with average stellar separations beyond 1,000 AU (1 AU being the Earth-Sun distance), yet the influence of such distant binary companions on planetary evolution remains largely unstudied. Unlike their tighter counterparts, the stellar orbits of wide binaries continually change under the influence of the Galactic tide and impulses from other passing stars. Here we report numerical simulations demonstrating that the variable nature of wide binary star orbits dramatically reshapes the planetary systems they host, typically Gyrs after formation. Contrary to previous understanding, wide binary companions may often strongly perturb planetary systems, triggering planetary ejections and exciting orbital eccentricities of surviving planets. Indeed, observed exoplanet eccentricities offer evidence of this; giant exoplanet orbits within wide binaries are statistically more eccentric than those around isolated stars. Both eccentricity distributions are well-reproduced when we assume isolated stars and wide binaries host similar planetary systems whose outermost giant planets are scattered beyond ~10 AU from their parent stars via early internal instabilities. Consequently, our results suggest that although wide binaries eventually truncate their planetary systems, most isolated giant exoplanet systems harbor additional distant, still undetected planets.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Secular Behavior of Exoplanets: Self-Consistency and Comparisons with the Planet-Planet Scattering Hypothesis

If mutual gravitational scattering among exoplanets occurs, then it may produce unique orbital properties. For example, two-planet systems that lie near the boundary between circulation and libration of their periapses could result if planet-planet scattering ejected a former third planet quickly, leaving one planet on an eccentric orbit and the other on a circular orbit. We first improve upon previous work that examined the apsidal behavior of known multiplanet systems by doubling the sample size and including observational uncertainties. This analysis recovers previous results that demonstrated that many systems lay on the apsidal boundary between libration and circulation. We then performed over 12,000 three-dimensional N-body simulations of hypothetical three-body systems that are unstable, but stabilize to two-body systems after an ejection. Using these synthetic two-planet systems, we test the planet-planet scattering hypothesis by comparing their apsidal behavior, over a range of viewing angles, to that of the observed systems and find that they are statistically consistent regardless of the multiplicity of the observed systems. Finally, we combine our results with previous studies to show that, from the sampled cases, the most likely planetary mass function prior to planet-planet scattering follows a power law with index -1.1. We find that this pre-scattering mass function predicts a mutual inclination frequency distribution that follows an exponential function with an index between -0.06 and -0.1.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Solar System Moons as Analogs for Compact Exoplanetary Systems

The field of exoplanetary science has experienced a recent surge of new systems that is largely due to the precision photometry provided by the Kepler mission. The latest discoveries have included compact planetary systems in which the orbits of the planets all lie relatively close to the host star, which presents interesting challenges in terms of formation and dynamical evolution. The compact exoplanetary systems are analogous to the moons orbiting the giant planets in our Solar System, in terms of their relative sizes and semi-major axes. We present a study that quantifies the scaled sizes and separations of the Solar System moons with respect to their hosts. We perform a similar study for a large sample of confirmed Kepler planets in multi-planet systems. We show that a comparison between the two samples leads to a similar correlation between their scaled sizes and separation distributions. The different gradients of the correlations may be indicative of differences in the formation and/or long-term dynamics of moon and planetary systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Long-Term Dynamical Evolution of Planetary Systems

This chapter concerns the long-term dynamical evolution of planetary systems from both theoretical and observational perspectives. We begin by discussing the planet-planet interactions that take place within our own Solar System. We then describe such interactions in more tightly-packed planetary systems. As planet-planet interactions build up, some systems become dynamically unstable, leading to strong encounters and ultimately either ejections or collisions of planets. After discussing the basic physical processes involved, we consider how these interactions apply to extrasolar planetary systems and explore the constraints provided by observed systems. The presence of a residual planetesimal disc can lead to planetary migration and hence cause instabilities induced by resonance crossing; however, such discs can also stabilise planetary systems. The crowded birth environment of a planetary system can have a significant impact: close encounters and binary companions can act to destabilise systems, or sculpt their properties. In the case of binaries, the Kozai mechanism can place planets on extremely eccentric orbits which may later circularise to produce hot Jupiters.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Tidal dissipation and eccentricity pumping: Implications for the depth of the secondary eclipse of 55 Cnc e

We use the super Earth 55 Cnc e as a case study to address an observable effect of tidal heating. We investigate whether planet-planet interactions can force the eccentricity of this planet to a level affecting the eclipse depth observed with Spitzer. Using the constant time lag tidal model, we first calculate the observed planet flux as a function of albedo and eccentricity, for different tidal dissipation constants and for two extreme cases: a planet with no heat redistribution and a planet with full heat redistribution. We derive the values of albedo and eccentricity that match the observed transit depth. We then perform N-body simulations of the planetary system including tides and General Relativity to follow the evolution of the eccentricity of planet e. We compare the range of eccentricities given by the simulations with the eccentricities required to alter the eclipse depth. We find that the eccentricity of planet e can be large enough to contribute at a measurable level to the thermal emission measured with Spitzer. This affects the constraints on the albedo of the planet, which can be as high as 0.9 (instead of 0.55 when ignoring tidal heating). We also derive a maximum value for the eccentricity of planet e directly from the eclipse depth: e<0.015 assuming Earth's dissipation constant. Transiting exoplanets in multiple planet systems are exceptional targets for testing tidal models because their tidal luminosity may be observable. Future multi-wavelengths observations of eclipse depth and phase curves (for instance with EChO and JWST) should allow us to better resolve the temperature map of these planets and break the degeneracy between albedo and tidal heating that remains for single band observations. In addition, an accurate determination of the eccentricity will make it possible to constrain the dissipation rate of the planet and to probe its internal structure.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A low mass for Mars from Jupiter's early gas-driven migration

Jupiter and Saturn formed in a few million years (Haisch et al. 2001) from a gas-dominated protoplanetary disk, and were susceptible to gas-driven migration of their orbits on timescales of only ~100,000 years (Armitage 2007). Hydrodynamic simulations show that these giant planets can undergo a two-stage, inward-then-outward, migration (Masset & Snellgrove 2001, Morbidelli & Crida 2007, Pierens & Nelson 2008). The terrestrial planets finished accreting much later (Klein et al. 2009), 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 Sun (Wetherill 1978, Hansen 2009) (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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Building Terrestrial Planets

This paper reviews our current understanding of terrestrial planets formation. The focus is on computer simulations of the dynamical aspects of the accretion process. Throughout the chapter, we combine the results of these theoretical models with geochemical, cosmochemical and chronological constraints, in order to outline a comprehensive scenario of the early evolution of our Solar System. Given that the giant planets formed first in the protoplanetary disk, we stress the sensitive dependence of the terrestrial planet accretion process on the orbital architecture of the giant planets and on their evolution. This suggests a great diversity among the terrestrial planets populations in extrasolar systems. Issues such as the cause for the different masses and accretion timescales between Mars and the Earth and the origin of water (and other volatiles) on our planet are discussed at depth.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Debris disks as signposts of terrestrial planet formation. II Dependence of exoplanet architectures on giant planet and disk properties

We present models for the formation of terrestrial planets, and the collisional evolution of debris disks, in planetary systems that contain multiple unstable gas giants. We previously showed that the dynamics of the giant planets introduces a correlation between the presence of terrestrial planets and debris disks. Here we present new simulations that show that this connection is qualitatively robust to changes in: the mass distribution of the giant planets, the width and mass distribution of the outer planetesimal disk, and the presence of gas in the disk. We discuss how variations in these parameters affect the evolution. Systems with equal-mass giant planets undergo the most violent instabilities, and these destroy both terrestrial planets and the outer planetesimal disks that produce debris disks. In contrast, systems with low-mass giant planets efficiently produce both terrestrial planets and debris disks. A large fraction of systems with low-mass outermost giant planets have stable gaps between these planets that are frequently populated by planetesimals. Planetesimal belts between outer giant planets may affect debris disk SEDs. If Earth-mass seeds are present in outer planetesimal disks, the disks radially spread to colder temperatures. We argue that this may explain the very low frequency of > 1 Gyr-old solar-type stars with observed 24 micron excesses. Among the (limited) set of configurations explored, the best candidates for hosting terrestrial planets at ~1 AU are stars older than 0.1-1 Gyr with bright debris disks at 70 micron but with no currently-known giant planets. These systems combine evidence for rocky building blocks, with giant planet properties least likely to undergo destructive dynamical evolution. We predict an anti-correlation between debris disks and eccentric giant planets, and a positive correlation between debris disks and terrestrial planets.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Effect of the stellar spin history on the tidal evolution of close-in planets

We investigate how the evolution of the stellar spin rate affects, and is affected by, planets in close orbits, via star-planet tidal interactions. To do this, we used a standard equilibrium tidal model to compute the orbital evolution of single planets orbiting both Sun-like stars and 0.1 M\odot M-dwarfs. We tested two stellar spin evolution profiles, one with fast initial rotation (P=1.2 day) and one with slow initial rotation (P=8 day). We tested the effect of varying the stellar and planetary dissipation and the planet's mass and initial orbital radius. Conclusions: Tidal evolution allows to differentiate the early behaviors of extremely close-in planets orbiting either a rapidly rotating star or a slowly rotating star. The early spin-up of the star allows the close-in planets around fast rotators to survive the early evolution. For planets around M-dwarfs, surviving the early evolution means surviving on Gyr timescales whereas for Sun-like stars the spin-down brings about late mergers of Jupiter planets. In light of this study, we can say that differentiating between one spin evolution from another given the present position of planets can be very tricky. Unless we can observe some markers of former evolution it is nearly impossible to distinguish the two very different spin profiles, let alone intermediate spin profiles. Though some conclusions can still be drawn from statistical distributions of planets around fully convective M-dwarfs. However, if the tidal evolution brings about a merger late in its history it can also entail a noticeable acceleration of the star in late ages, so that it is possible to have old stars that spin rapidly. This raises the question of better constraining the age of stars.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Mini-Oort clouds: Compact isotropic planetesimal clouds from planet-planet scattering

Starting from planetary systems with three giant planets and an outer disk of planetesimals, we use dynamical simulations to show how dynamical instabilities can transform planetesimal disks into 100-1000 AU-scale isotropic clouds. The instabilities involve a phase of planet-planet scattering that concludes with the ejection of one or more planets and the inward-scattering of the surviving gas giant(s) to remove them from direct dynamical contact with the planetesimals. "Mini-Oort clouds" are thus formed from scattered planetesimals whose orbits are frozen by the abrupt disappearance of the perturbing giant planet. Although the planetesimal orbits are virtually isotropic, the surviving giant planets tend to have modest inclinations (typically ~10 degrees) with respect to the initial orbital plane. The collisional lifetimes of mini-Oort clouds are long (10 Myr to >10 Gyr) and there is a window of ~100 Myr or longer during which they produce spherical clouds of potentially observable dust at 70 microns. If the formation channel for hot Jupiters commonly involves planetary close encounters, we predict a correlation between this subset of extrasolar planetary systems and mini-Oort clouds.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Planet-Planet Scattering Alone Cannot Explain the Free-Floating Planet Population

Recent gravitational microlensing observations predict a vast population of free-floating giant planets that outnumbers main sequence stars almost twofold. A frequently-invoked mechanism for generating this population is a dynamical instability that incites planet-planet scattering and the ejection of one or more planets in isolated main sequence planetary systems. Here, we demonstrate that this process alone probably cannot represent the sole source of these galactic wanderers. By using straightforward quantitative arguments and N-body simulations, we argue that the observed number of exoplanets exceeds the plausible number of ejected planets per system from scattering. Thus, other potential sources of free-floaters, such as planetary stripping in stellar clusters and post-main-sequence ejection, must be considered.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Primordial origins of Earth's carbon

In this chapter we review the astrophysical origins of Earth's carbon, starting from the products of the Big Bang and culminating with the Earth's formation. We review the measured compositions of different primitive objects including comets, various classes of meteorites and interstellar dust particles. We discuss the composition of the Solar Nebula, especially with regards to the distribution of volatiles such as carbon. We discuss dynamical models of planetary formation from planetesimals and planetary embryos, and the timescale for volatile delivery to the growing Earth from different sources. Finally, we review Earth's carbon reservoirs. Throughout the chapter we highlight open questions related to planet formation, meteoritics, and geochemistry.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Compositional Diversity of Extrasolar Terrestrial Planets: II. Migration Simulations

Prior work has found that a variety of terrestrial planetary compositions are expected to occur within known extrasolar planetary systems. However, such studies ignored the effects of giant planet migration, which is thought to be very common in extra-solar systems. Here we present calculations of the compositions of terrestrial planets that formed in dynamical simulations incorporating varying degrees of giant planet migration. We used chemical equilibrium models of the solid material present in the disks of five known planetary host stars: the Sun, GJ 777, HD4203, HD19994 and HD213240. Giant planet migration has a strong effect on the compositions of simulated terrestrial planets as the migration results large-scale mixing between terrestrial planet building blocks that condensed at a range of temperatures. This mixing acts to 1) increase the typical abundance of Mg-rich silicates in the terrestrial planets feeding zones and thus increase the frequency of planets with Earth-like compositions compared with simulations with static giant planet orbits; and 2) drastically increase the efficiency of the delivery of hydrous phases (water and serpentine) to terrestrial planets and thus produce water worlds and/or wet Earths. Our results demonstrate that although a wide variety of terrestrial planet compositions can still be produced, planets with Earth-like compositions should be common within extrasolar planetary systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

55 Cancri: A Coplanar Planetary System that is Likely Misaligned with its Star

Although the 55 Cnc system contains multiple, closely packed planets that are presumably in a coplanar configuration, we use numerical simulations to demonstrate that they are likely to be highly inclined to their parent star's spin axis. Due to perturbations from its distant binary companion, this planetary system precesses like a rigid body about its parent star. Consequently, the parent star's spin axis and the planetary orbit normal likely diverged long ago. Because only the projected separation of the binary is known, we study this effect statistically, assuming an isotropic distribution for wide binary orbits. We find that the most likely projected spin-orbit angle is ~50 degrees, with a ~30% chance of a retrograde configuration. Transit observations of the innermost planet - 55 Cnc e - may be used to verify these findings via the Rossiter-McLaughlin effect. 55 Cancri may thus represent a new class of planetary systems with well-ordered, coplanar orbits that are inclined with respect to the stellar equator.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

55 Cancri: Stellar Astrophysical Parameters, a Planet in the Habitable Zone, and Implications for the Radius of a Transiting Super-Earth

The bright star 55 Cancri is known to host five planets, including a transiting super-Earth. The study presented here yields directly determined values for 55 Cnc's stellar astrophysical parameters based on improved interferometry: $R=0.943 \pm 0.010 R_{\odot}$, $T_{\rm EFF} = 5196 \pm 24$ K. We use isochrone fitting to determine 55 Cnc's age to be 10.2 $\pm$ 2.5 Gyr, implying a stellar mass of $0.905 \pm 0.015 M_{\odot}$. Our analysis of the location and extent of the system's habitable zone (0.67--1.32 AU) shows that planet f, with period $\sim$ 260 days and $M \sin i = 0.155 M_{Jupiter}$, spends the majority of the duration of its elliptical orbit in the circumstellar habitable zone. Though planet f is too massive to harbor liquid water on any planetary surface, we elaborate on the potential of alternative low-mass objects in planet f's vicinity: a large moon, and a low-mass planet on a dynamically stable orbit within the habitable zone. Finally, our direct value for 55 Cancri's stellar radius allows for a model-independent calculation of the physical diameter of the transiting super-Earth 55 Cnc e ($\sim 2.05 \pm 0.15 R_{\earth}$), which, depending on the planetary mass assumed, implies a bulk density of 0.76 $ρ_{\earth}$ or 1.07 $ρ_{\earth}$.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Debris disks as signposts of terrestrial planet formation

Circumstantial evidence suggests that most known extra-solar planetary systems are survivors of violent dynamical instabilities. Here we explore how giant planet instabilities affect the formation and survival of terrestrial planets. We simulate planetary system evolution around Sun-like stars from initial conditions that comprise: an inner disk of planetesimals and planetary embryos, three giant planets at Jupiter-Saturn distances, and a massive outer planetesimal disk. We then calculate dust production rates and debris disk SEDs assuming that each planetesimal particle represents an ensemble of smaller bodies in collisional equilibrium. We predict a strong correlation between the presence of terrestrial planets and debris disks, mediated by the giant planets. Strong giant planet instabilities destroy all rocky material - including fully-formed terrestrial planets if the instabilities occur late - along with the icy planetesimals. Stable or weakly unstable systems allow terrestrial planets to accrete and significant dust to be produced in their outer regions. Stars older than ~100 Myr with bright cold dust emission (at ~70 microns) signpost the dynamically calm environments conducive to efficient terrestrial accretion. We predict that while the typical eccentricities of terrestrial planets are small, there should exist a novel class of terrestrial planet system whose single planet undergoes large amplitude oscillations in eccentricity and inclination. By scaling to the observed semimajor axis distribution of giant exoplanets, we estimate that terrestrial exoplanets in the same systems should be a few times more abundant at 0.5 AU than giant or terrestrial exoplanets at 1 AU. Finally, we discuss the Solar System, which appears to be unusual in combining a rich terrestrial planet system with a low dust content.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The debris disk - terrestrial planet connection

The eccentric orbits of the known extrasolar giant planets provide evidence that most planet-forming environments undergo violent dynamical instabilities. Here, we numerically simulate the impact of giant planet instabilities on planetary systems as a whole. We find that populations of inner rocky and outer icy bodies are both shaped by the giant planet dynamics and are naturally correlated. Strong instabilities -- those with very eccentric surviving giant planets -- completely clear out their inner and outer regions. In contrast, systems with stable or low-mass giant planets form terrestrial planets in their inner regions and outer icy bodies produce dust that is observable as debris disks at mid-infrared wavelengths. Fifteen to twenty percent of old stars are observed to have bright debris disks (at wavelengths of ~70 microns) and we predict that these signpost dynamically calm environments that should contain terrestrial planets.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Tidal evolution of planets around brown dwarfs

The tidal evolution of planets orbiting brown dwarfs (BDs) presents an interesting case study because BDs' terrestrial planet forming region is located extremely close-in. In fact, the habitable zones of BDs range from roughly 0.001 to 0.03 AU and for the lowest-mass BDs are located interior to the Roche limit. In contrast with stars, BDs spin up as they age. Thus, the corotation distance moves inward. This has important implications for the tidal evolution of planets around BDs. We used a standard equilibrium tidal model to compute the orbital evolution of a large ensemble of planet-BD systems. We tested the effect of numerous parameters such as the initial semi-major axis and eccentricity, the rotation period of the BD, the masses of both the BD and planet, and the tidal dissipation factors. We find that all planets that form at or beyond the corotation distance and with initial eccentricities smaller than \sim 0.1 are repelled from the BD. Some planets initially interior to corotation can survive if their inward tidal evolution is slower than the BD's spin evolution, but most initially close-in planets fall onto the BD. We find that the most important parameter for the tidal evolution is the initial orbital distance with respect to the corotation distance. Some planets can survive in the habitable zone for Gyr timescales, although in many cases the habitable zone moves inward past the planet's orbit in just tens to hundreds of Myr. Surviving planets can have orbital periods of less than 10 days (as small as 10 hrs), so they could be observable by transit.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Two phase, inward-then-outward migration of Jupiter and Saturn in the gaseous Solar Nebula

It has recently been shown that the terrestrial planets and asteroid belt can be reproduced if the giant planets underwent an inward-then-outward migration (the "Grand Tack"; Walsh et al 2011). Inward migration occurs when Jupiter opens a gap and type II migrates inward. The planets "tack" and migrate outward when Saturn reaches the gap-opening mass and is caught in the 3:2 resonance with Jupiter. The aim is to test the viability of the Grand Tack model and to study the dynamical evolution of Jupiter and Saturn during their growth from 10 Earth masses cores. We have performed numerical simulations using a grid-based hydrodynamical code. Most of our simulations assume an isothermal equation of state for the disk but a subset use a fully-radiative version of the code. For an isothermal disk the two phase migration of Jupiter and Saturn is very robust and independent of the mass-growth history of these planets provided the disk is cool enough. For a radiative disk the we find some outcomes with two phase migrations and others with more complicated behavior. We construct a simple, 1-D model of an evolving viscous disk to calculate the evolution of the disk's radiative properties: the disk transitions from radiative to isothermal from its outermost regions inward in time. We show that a two-phase migration is a natural outcome at late times even under the limiting assumption that isothermal conditions are required. Thus, our simulations provide strong support for the Grand Tack scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Generalized Milankovitch Cycles and Longterm Climatic Habitability

Although the Earth's orbit is never far from circular, terrestrial planets around other stars might experience substantial changes in eccentricity that could lead to climate changes, including possible "phase transitions" such as the snowball transition (or its opposite). There is evidence that Earth has gone through at least one globally frozen, "snowball" state in the last billion years, which it is thought to have exited after several million years because global ice-cover shut off the carbonate-silicate cycle, thereby allowing greenhouse gases to build up to sufficient concentration to melt the ice. Due to the positive feedback caused by the high albedo of snow and ice, susceptibility to falling into snowball states might be a generic feature of water-rich planets with the capacity to host life. This paper has two main thrusts. First, we revisit one-dimensional energy balance climate models as tools for probing possible climates of exoplanets, investigate the dimensional scaling of such models, and introduce a simple algorithm to treat the melting of the ice layer on a globally-frozen planet. We show that if a terrestrial planet undergoes Milankovitch-like oscillations of eccentricity that are of great enough magnitude, it could melt out of a snowball state. Second, we examine the kinds of variations of eccentricity that a terrestrial planet might experience due to the gravitational influence of a giant companion. We show that a giant planet on a sufficiently eccentric orbit can excite extreme eccentricity oscillations in the orbit of a habitable terrestrial planet. More generally, these two results demonstrate that the longterm habitability (and astronomical observables) of a terrestrial planet can depend on the detailed architecture of the planetary system in which it resides.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Habitable Climates: The Influence of Eccentricity

In the outer regions of the habitable zone, the risk of transitioning into a globally frozen "snowball" state poses a threat to the habitability of planets with the capacity to host water-based life. We use a one-dimensional energy balance climate model (EBM) to examine how obliquity, spin rate, orbital eccentricity, and ocean coverage might influence the onset of such a snowball state. For an exoplanet, these parameters may be strikingly different from the values observed for Earth. Since, for constant semimajor axis, the annual mean stellar irradiation scales with (1-e^2)^(-1/2), one might expect the greatest habitable semimajor axis (for fixed atmospheric composition) to scale as (1-e^2)^(-1/4). We find that this standard ansatz provides a reasonable lower bound on the outer boundary of the habitable zone, but the influence of obliquity and ocean fraction can be profound in the context of planets on eccentric orbits. For planets with eccentricity 0.5, our EBM suggests that the greatest habitable semimajor axis can vary by more than 0.8 AU (78%!) depending on obliquity, with higher obliquity worlds generally more stable against snowball transitions. One might also expect that the long winter at an eccentric planet's apoastron would render it more susceptible to global freezing. Our models suggest that this is not a significant risk for Earth-like planets around Sun-like stars since such planets are buffered by the thermal inertia provided by oceans covering at least 10% of their surface. Since planets on eccentric orbits spend much of their year particularly far from the star, such worlds might turn out to be especially good targets for direct observations with missions such as TPF-Darwin. Nevertheless, the extreme temperature variations achieved on highly eccentric exo-Earths raise questions about the adaptability of life to marginally or transiently habitable conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Planet-planet scattering in planetesimal disks II: Predictions for outer extrasolar planetary systems

We develop an idealized dynamical model to predict the typical properties of outer extrasolar planetary systems, at radii beyond 5 AU. Our hypothesis is that dynamical evolution in outer planetary systems is controlled by a combination of planet-planet scattering and planetary interactions with an exterior disk of small bodies ("planetesimals"). Using 5,000 long duration N-body simulations, we follow the evolution of three planets surrounded by a 50 Earth mass primordial planetesimal disk. For large planet masses (above that of Saturn) the influence of the disk is modest, and we recover the observed eccentricity distribution of extrasolar planets (observed primarily at smaller radii). We explain the observed mass dependence of the eccentricity by invoking strong correlations between planet masses in the same system. For lower mass planets we observe diverse dynamical behavior: strong scattering events, sudden jumps in eccentricity due to resonance crossings, and re-circularization of scattered low-mass planets in the disk. We present distributions of the final eccentricity and inclination, and discuss how they vary with planet mass and initial system architecture. We predict a transition to lower eccentricities for low mass planets at radii where disks influence the dynamics. Radial velocity measurements capable of detecting planets with K~5 m/s and periods in excess of 10 years will constrain this regime. We also study the population of resonant and non-resonant multiple planet systems. We show that, among systems with Jupiter-mass planets that avoid close encounters, the planetesimal disk acts as a damping mechanism that frequently populates mean motion resonances. Resonant chains ought to be common among massive planets in outer planetary systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Roles of Tidal Evolution and Evaporative Mass Loss in the Origin of CoRoT-7 b

CoRoT-7 b is the first confirmed rocky exoplanet, but, with an orbital semi-major axis of 0.0172 AU, its origins may be unlike any rocky planet in our solar system. In this study, we consider the roles of tidal evolution and evaporative mass loss in CoRoT-7 b's history, which together have modified the planet's mass and orbit. If CoRoT-7 b has always been a rocky body, evaporation may have driven off almost half its original mass, but the mass loss may depend sensitively on the extent of tidal decay of its orbit. As tides caused CoRoT-7 b's orbit to decay, they brought the planet closer to its host star, thereby enhancing the mass loss rate. Such a large mass loss also suggests the possibility that CoRoT-7 b began as a gas giant planet and had its original atmosphere completely evaporated. In this case, we find that CoRoT-7 b's original mass probably didn't exceed 200 Earth masses (about 2/3 of a Jupiter mass). Tides raised on the host star by the planet may have significantly reduced the orbital semi-major axis, perhaps causing the planet to migrate through mean-motion resonances with the other planet in the system, CoRoT-7 c. The coupling between tidal evolution and mass loss may be important not only for CoRoT-7 b but also for other close-in exoplanets, and future studies of mass loss and orbital evolution may provide insight into the origin and fate of close-in planets, both rocky and gaseous.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Tidal Effects on the Habitability of Exoplanets: The Case of GJ 581 d

Tides may be crucial to the habitability of exoplanets. If such planets form around low-mass stars, then those in the circumstellar habitable zone will be close enough to their host stars to experience strong tidal forces. Tides may result in orbital decay and circularization, evolution toward zero obliquity, a fixed rotation rate (not necessarily synchronous), and substantial internal heating. Due to tidal effects, the range of habitable orbital locations may be quite different from that defined by the traditional concept of a habitable zone (HZ) based on stellar insolation, atmospheric effects, and liquid water on a planet's surface. Tidal heating may make locations within the traditional HZ too hot, while planets outside the traditional zone could be rendered quite habitable due to tides. Here we consider these effects on the exoplanet GJ 581 d.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Al-26 and the formation of the Solar System from a molecular cloud contaminated by Wolf-Rayet winds

In agreement with previous work, we show that the presence of the short-lived radionuclide Al-26 in the early Solar System was unlikely (<2% a priori probability) to be the result of direct introduction of supernova ejecta into the gaseous disk during the Class II stage of protosolar evolution. We also show that any Bondi-Hoyle accretion of contaminated residual gas from the natal star cluster made a negligible contribution to the primordial Al-26 inventory of the Solar System. These results are consistent with the absence of the oxygen isotopic signature expected with any late introduction of supernova ejecta into the protoplanetary disk. Instead, the presence of Al-26 in the oldest Solar System solids (calcium-aluminum-rich inclusions or CAIs) and its apparent uniform distribution with the inferred canonical Al-26/Al-27 ratio of (4.5-5)E-5 support the inheritance of Al-26 from the parent giant molecular cloud. We propose that this radionuclide originated in a prior generation of massive stars that formed in the same host molecular cloud as the Sun and contaminated that cloud by Wolf-Rayet winds. We calculated the Galactic distribution of Al-26/Al-27 ratios that arise from such contamination using the established embedded cluster mass and stellar initial mass functions, published nucleosynthetic yields from the winds of massive stars, and by assuming rapid and uniform mixing into the cloud. Although our model predicts that the majority of stellar systems contain no Al-26 from massive stars, and that the a priori probability that the Al-26/Al-27 ratio will reach or exceed the canonical Solar System value is only ~6%, the maximum in the distribution of non-zero values is close to the canonical ratio.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Building the Terrestrial Planets: Constrained Accretion in the Inner Solar System

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 assume that Jupiter and Saturn are fully-formed at the start of each simulation, and test orbital configurations that are both consistent with and contrary to the "Nice model." We find that a configuration with Jupiter and Saturn on circular orbits forms low-eccentricity terrestrial planets and a water-rich Earth on the correct timescale, but Mars' mass is too large by a factor of 5-10 and embryos are often stranded in the asteroid belt. A configuration with Jupiter and Saturn in their current locations but with slightly higher initial eccentricities (e = 0.07-0.1) produces a small Mars, an embryo-free asteroid belt, and a reasonable Earth analog but rarely allows water delivery to Earth. None of the configurations we tested reproduced all the observed constraints. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

CoRoT-7 b: Super-Earth or Super-Io?

CoRoT-7 b, a planet about 70% larger than the Earth orbiting a Sun-like star, is the first-discovered rocky exoplanet, and hence has been dubbed a "super-Earth". Some initial studies suggested that since the planet is so close to its host star, it receives enough insolation to partially melt its surface. However, these past studies failed to take into consideration the role that tides may play in this system. Even if the planet's eccentricity has always been zero, we show that tidal decay of semi-major axis could have been large enough that the planet formed on a wider orbit which received less insolation. Moreover, CoRoT-7 b could be tidally heated at a rate that dominates its geophysics and drives extreme volcanism. In this case, CoRoT-7 b is a "super-Io" that, like Jupiter's volcanic moon, is dominated by volcanism and rapid resurfacing. Such heating could occur with an eccentricity of just 10^-5. This small value could be driven by CoRoT-7 c if its own eccentricity is larger than ~10^-4. CoRoT-7 b may be the first of a class of planetary super-Ios likely to be revealed by the CoRoT and Kepler spacecraft.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Planet-planet scattering leads to tightly packed planetary systems

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The HD 40307 Planetary System: Super-Earths or Mini-Neptunes?

Three planets with minimum masses less than 10 Earth masses 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 0.0001. 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 0.001, 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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

The Successful Prediction of the Extrasolar Planet HD 74156 d

Most of the first-discovered extrasolar multi-planet systems were found to lie close to dynamically unstable configurations. However a few observed multi-planet systems (e.g. HD 74156) did not show this trait. Those systems could share this property if they contain an additional planet in between those that are known. Previous investigations identified the properties of hypothetical planets that would place these systems near instability. The hypothetical planet in HD 74156 was expected to have a mass about equal to that of Saturn, a semi-major axis between 0.9 and 1.4 AU, and an eccentricity less than 0.2. HD 74156 d, a planet with a mass of 1.3 Saturn masses at 1.04 AU with an eccentricity of 0.25, was recently announced. We have reanalyzed all published data on this system in order to place tighter constraints on the properties of the new planet. We find two possible orbits for this planet, one close to that already identified and another (with a slightly better fit to the data) at ~0.89 AU. We also review the current status of other planet predictions, discuss the observed single planet systems, and suggest other systems which may contain planets in between those that are already known. The confirmation of the existence of HD 74156 d suggests that planet formation is an efficient process, and planetary systems should typically contain many planets.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Formation of Earth-like Planets During and After Giant Planet Migration

Close-in giant planets are thought to have formed in the cold outer regions of planetary systems and migrated inward, passing through the orbital parameter space occupied by the terrestrial planets in our own Solar System. We present dynamical simulations of the effects of a migrating giant planet on a disk of protoplanetary material and the subsequent evolution of the planetary system. We numerically investigate the dynamics of post-migration planetary systems over 200 million years using models with a single migrating giant planet, one migrating and one non-migrating giant planet, and excluding the effects of a gas disk. Material that is shepherded in front of the migrating giant planet by moving mean motion resonances accretes into "hot Earths", but survival of these bodies is strongly dependent on dynamical damping. Furthermore, a significant amount of material scattered outward by the giant planet survives in highly excited orbits; the orbits of these scattered bodies are then damped by gas drag and dynamical friction over the remaining accretion time. In all simulations Earth-mass planets accrete on approximately 100 Myr timescales, often with orbits in the Habitable Zone. These planets range in mass and water content, with both quantities increasing with the presence of a gas disk and decreasing with the presence of an outer giant planet. We use scaling arguments and previous results to derive a simple recipe that constrains which giant planet systems are able to form and harbor Earth-like planets in the Habitable Zone, demonstrating that roughly one third of the known planetary systems are potentially habitable.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Spectroscopic Properties of Cool Stars in the SDSS: An Analysis of Magnetic Activity and a Search for Subdwarfs

We present a spectroscopic analysis of nearly 8000 late-type dwarfs in the Sloan Digital Sky Survey. Using the Halpha emission line as an activity indicator, we investigate the fraction of active stars as a function of spectral type and find a peak near type M8, confirming previous results. In contrast to past findings, we find that not all M7-M8 stars are active. We show that this may be a selection effect of the distance distributions of previous samples, as the active stars appear to be concentrated near the Galactic Plane. We also examine the activity strength (ratio of the luminosity emitted in Halpha to the bolometric luminosity) for each star, and find that the mean activity strength is constant over the range M0-M5 and declines at later types. The decline begins at a slightly earlier spectral type than previously found. We explore the effect that activity has on the broadband photometric colors and find no significant differences between active and inactive stars. We also carry out a search for subdwarfs using spectroscopic metallicity indicators, and find 60 subdwarf candidates. Several of these candidates are near the extreme subdwarf boundary. The spectroscopic subdwarf candidates are redder by \~0.2 magnitudes in g-r compared to disk dwarfs at the same r-i color.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

The formation and habitability of terrestrial planets in the presence of hot jupiters

`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 (Armitage 2003). We present results of simulations of terrestrial planet formation in the presence of hot 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 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 hot jupiters.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Making other Earths: Dynamical Simulations of Terrestrial Planet Formation and Water Delivery

We present results from 42 simulations of late stage planetary accretion, focusing on the delivery of volatiles (primarily water) to the terrestrial planets. Our simulations include both planetary "embryos" (defined as Moon to Mars sized protoplanets) and planetesimals, assuming that the embryos formed via oligarchic growth. We investigate volatile delivery as a function of Jupiter's mass, position and eccentricity, the position of the snow line, and the density (in solids) of the solar nebula. In all simulations, we form 1-4 terrestrial planets inside 2 AU, which vary in mass and volatile content. In 42 simulations we have formed 43 planets between 0.8 and 1.5 AU, including 11 "habitable" planets between 0.9 and 1.1 AU. These planets range from dry worlds to "water worlds" with 100+ oceans of water (1 ocean = 1.5x10^24 g), and vary in mass between 0.23 and 3.85 Earth masses. There is a good deal of stochastic noise in these simulations, but the most important parameter is the planetesimal mass we choose, which reflects the surface density in solids past the snow line. A high density in this region results in the formation of a smaller number of terrestrial planets with larger masses and higher water content, as compared with planets which form in systems with lower densities. We find that an eccentric Jupiter produces drier terrestrial planets with higher eccentricities than a circular one. In cases with Jupiter at 7 AU, we form what we call "super embryos," 1-2 Earth mass protoplanets which can serve as the accretion seeds for 2+ Earth mass planets with large water contents.

0 citations0 reviewsNo code linkedOpen access