Reproducing the large Earth/Mars mass ratio requires a strong mass depletion in solids within the protoplanetary disc 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 discs of planetesimals and planetary embryos with various radial density gradients and including Jupiter and Saturn on nearly circular and coplanar orbits. Discs with shallow density gradients reproduce the dynamical excitation of the asteroid belt by gravitational self-stirring but form Mars analogues significantly more massive than the real planet. In contrast, a disc with a surface density gradient proportional to r(exp −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 disc 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.
Abstract Reconstructions of the paleoclimate indicate that ancient climatic fluctuations on Earth are often correlated with variations in its orbital elements. However, the chaos inherent in the ...solar system’s orbital evolution prevents numerical simulations from confidently predicting Earth’s past orbital evolution beyond 50–100 Myr. Gravitational interactions among the Sun’s planets and asteroids are believed to set this limiting time horizon, but most prior works approximate the solar system as an isolated system and neglect our surrounding Galaxy. Here we present simulations that include the Sun’s nearby stellar population, and we find that close-passing field stars alter our entire planetary system’s orbital evolution via their gravitational perturbations on the giant planets. This shortens the timespan over which Earth’s orbital evolution can be definitively known by a further ∼10%. In particular, in simulations that include an exceptionally close passage of the Sun-like star HD 7977 2.8 Myr ago, new sequences of Earth’s orbital evolution become possible in epochs before ∼50 Myr ago, which includes the Paleocene–Eocene Thermal Maximum. Thus, simulations predicting Earth’s past orbital evolution before ∼50 Myr ago must consider the additional uncertainty from passing stars, which can open new regimes of past orbital evolution not seen in previous modeling efforts.
The main asteroid belt (MB) is low in mass but dynamically excited. Here we propose a new mechanism to excite the MB during the giant planet (the "Nice model") instability, which is expected to ...feature repeated close encounters between Jupiter and one or more ice giants ("jumping Jupiter" or JJ). We show that, when Jupiter temporarily reaches a high-enough level of excitation, both in eccentricity and inclination, it induces strong forced vectors of eccentricity and inclination across the MB region. Because during the JJ instability Jupiter's orbit "jumps" around, the forced vectors keep changing both in magnitude and phase throughout the whole MB region. The entire cold primordial MB is thus excited as a natural outcome of the JJ instability. The level of such an excitation, however, is typically larger than the current orbital excitation observed in the MB. We show that the subsequent evolution of the solar system is capable of reshaping the resultant overexcited MB to its present-day orbital state, and that a strong mass depletion (∼90%) is associated with the JJ instability phase and its subsequent evolution throughout the age of the solar system.
We have collected transit times for the TRAPPIST-1 system with the Spitzer Space Telescope over four years. We add to these ground-based, HST, and K2 transit-time measurements, and revisit an N-body ...dynamical analysis of the seven-planet system using our complete set of times from which we refine the mass ratios of the planets to the star. We next carry out a photodynamical analysis of the Spitzer light curves to derive the density of the host star and the planet densities. We find that all seven planets' densities may be described with a single rocky mass-radius relation which is depleted in iron relative to Earth, with Fe 21 wt% versus 32 wt% for Earth, and otherwise Earth-like in composition. Alternatively, the planets may have an Earth-like composition but enhanced in light elements, such as a surface water layer or a core-free structure with oxidized iron in the mantle. We measure planet masses to a precision of 3%-5%, equivalent to a radial-velocity (RV) precision of 2.5 cm s−1, or two orders of magnitude more precise than current RV capabilities. We find the eccentricities of the planets are very small, the orbits are extremely coplanar, and the system is stable on 10 Myr timescales. We find evidence of infrequent timing outliers, which we cannot explain with an eighth planet; we instead account for the outliers using a robust likelihood function. We forecast JWST timing observations and speculate on possible implications of the planet densities for the formation, migration, and evolution of the planet system.
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 and Uranus mass ratio is ~1.18. We show that similar-sized 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 five to ten planetary embryos with initial masses of ~3–6 M⊕. We conclude that accretion from a population of planetary embryos is a plausible scenario for the origin of Uranus and Neptune.
Jupiter and Saturn formed in a few million years (ref. 1) from a gas-dominated protoplanetary disk, and were susceptible to gas-driven migration of their orbits on timescales of only ∼100,000 years ...(ref. 2). Hydrodynamic simulations show that these giant planets can undergo a two-stage, inward-then-outward, migration. The terrestrial planets finished accreting much later, 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 (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.
The planetary building blocks that formed in the terrestrial planet region were likely very dry, yet water is comparatively abundant on Earth. Here we review the various mechanisms proposed for the ...origin of water on the terrestrial planets. Various in-situ mechanisms have been suggested, which allow for the incorporation of water into the local planetesimals in the terrestrial planet region or into the planets themselves from local sources, although all of those mechanisms have difficulties. Comets have also been proposed as a source, although there may be problems fitting isotopic constraints, and the delivery efficiency is very low, such that it may be difficult to deliver even a single Earth ocean of water this way. The most promising route for water delivery is the accretion of material from beyond the snow line, similar to carbonaceous chondrites, that is scattered into the terrestrial planet region as the planets are growing. Two main scenarios are discussed in detail. First is the classical scenario in which the giant planets begin roughly in their final locations and the disk of planetesimals and embryos in the terrestrial planet region extends all the way into the outer asteroid belt region. Second is the Grand Tack scenario, where early inward and outward migration of the giant planets implants material from beyond the snow line into the asteroid belt and terrestrial planet region, where it can be accreted by the growing planets. Sufficient water is delivered to the terrestrial planets in both scenarios. While the Grand Tack scenario provides a better fit to most constraints, namely the small mass of Mars, planets may form too fast in the nominal case discussed here. This discrepancy may be reduced as a wider range of initial conditions is explored. Finally, we discuss several more recent models that may have important implications for water delivery to the terrestrial planets.
•Sufficient water is delivered to terrestrial planets in the Grand Tack scenario.•Radial distribution and dynamical excitation are comparable to the Solar System.•Late giant impacts are relatively ...rare with current set of initial conditions.
A new model for terrestrial planet formation (Hansen 2009. Astrophys. J., 703, 1131–1140; Walsh, K.J., et al. 2011. Nature, 2011, 206–209) 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. (Walsh, K.J., et al. 2011. Nature, 2011, 206–209) 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 3AU, and the outer belt is filled with bodies originating from between and beyond the giant planets (which are hereafter referred to as ‘primitive’ bodies).
One implication of the truncation mechanism proposed in Walsh et al. (Walsh, K.J., et al. 2011. Nature, 2011, 206–209) is the scattering of primitive planetesimals onto planet-crossing orbits during the formation of the planets. We find here that the planets will accrete on order 1–2% of their total mass from these bodies. 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. The radial distribution of the planetary masses and the dynamical excitation of their orbits are a good match to the observed system. However, we find that a truncated disk leads to formation timescales more rapid than suggested by radiometric chronometers. In particular, the last giant impact is typically earlier than 20Myr, and a substantial amount of mass is accreted after that event. This is at odds with the dating of the Moon-forming impact and the estimated amount of mass accreted by Earth following 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 those constraints.
Abstract
Evidence for cometary activity beyond Jupiter’s 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 Jupiter-family comets (JFCs) can be representative of any primitive material. Here we model the coupled thermal and dynamical evolution of 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 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.
Radial velocity monitoring has found the signature of a M sin i = 1.3 M(sub ⨁) planet located within the habitable zone (HZ) of Proxima Centauri. Despite a hotter past and an active host star, the ...planet Proxima b could have retained enough volatiles to sustain surface habitability. Here we use a 3D Global Climate Model (GCM) to simulate the atmosphere and water cycle of Proxima b 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. Any low-obliquity, low-eccentricity planet within the HZ of its star should be in one of the climate regimes discussed here. 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 (~10 mbar of CO2 and 1 bar of N2). If the planet is dryer, ~0.5 bar or 1.5 bars of CO2 (for asynchronous or synchronous rotation, respectively) suffice to prevent the trapping of any arbitrary, small water inventory into polar or nightside ice caps. We produce reflection and emission spectra and phase curves for the simulated climates. We find that atmospheric characterization will be possible via direct imaging with forthcoming large telescopes. The angular separation of 7 λ/D at 1 μm (with the E-ELT) and a contrast of ~10(exp -7) will enable high-resolution spectroscopy and the search for molecular signatures, including H2O, O2, and CO2. The observation of thermal phase curves can be attempted with the James Webb Space Telescope, thanks to a contrast of 2 × 10(exp -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 the surface of this exoplanet is habitable.