-
Mercury's formation within the Early Instability Scenario
Authors:
Matthew S. Clement,
John E. Chambers,
Nathan A. Kaib,
Sean N. Raymond,
Alan P. Jackson
Abstract:
The inner solar system's modern orbital architecture provides inferences into the epoch of terrestrial planet formation; a ~100 Myr time period of planet growth via collisions with planetesimals and other proto-planets. While classic numerical simulations of this scenario adequately reproduced the correct number of terrestrial worlds, their semi-major axes and approximate formation timescales, the…
▽ More
The inner solar system's modern orbital architecture provides inferences into the epoch of terrestrial planet formation; a ~100 Myr time period of planet growth via collisions with planetesimals and other proto-planets. While classic numerical simulations of this scenario adequately reproduced the correct number of terrestrial worlds, their semi-major axes and approximate formation timescales, they struggled to replicate the Earth-Mars and Venus-Mercury mass ratios. In a series of past independent investigations, we demonstrated that Mars' mass is possibly the result of Jupiter and Saturn's early orbital evolution, while Mercury's diminutive size might be the consequence of a primordial mass deficit in the region. Here, we combine these ideas in a single modeled scenario designed to simultaneously reproduce the formation of all four terrestrial planets and the modern orbits of the giant planets in broad strokes. By evaluating our Mercury analogs' core mass fractions, masses, and orbital offsets from Venus, we favor a scenario where Mercury forms through a series of violent erosive collisions between a number of ~Mercury-mass embryos in the inner part of the terrestrial disk. We also compare cases where the gas giants begin the simulation locked in a compact 3:2 resonant configuration to a more relaxed 2:1 orientation and find the former to be more successful. In 2:1 cases, the entire Mercury-forming region is often depleted due to strong sweeping secular resonances that also tend to overly excite the orbits of Earth and Venus as they grow. While our model is quite successful at replicating Mercury's massive core and dynamically isolated orbit, the planets' low mass remains extremely challenging to match. Finally, we discuss the merits and drawbacks of alternative evolutionary scenarios and initial disk conditions.
△ Less
Submitted 23 January, 2023;
originally announced January 2023.
-
Mercury as the relic of Earth and Venus' outward migration
Authors:
Matthew S. Clement,
Sean N. Raymond,
John E. Chambers
Abstract:
In spite of substantial advancements in simulating planet formation, the planet Mercury's diminutive mass, isolated orbit, and the absence of planets with shorter orbital periods in the solar system continue to befuddle numerical accretion models. Recent studies have shown that, if massive embryos (or even giant planet cores) formed early in the innermost parts of the Sun's gaseous disk, they woul…
▽ More
In spite of substantial advancements in simulating planet formation, the planet Mercury's diminutive mass, isolated orbit, and the absence of planets with shorter orbital periods in the solar system continue to befuddle numerical accretion models. Recent studies have shown that, if massive embryos (or even giant planet cores) formed early in the innermost parts of the Sun's gaseous disk, they would have migrated outward. This migration may have reshaped the surface density profile of terrestrial planet-forming material and generated conditions favorable to the formation of Mercury-like planets. Here, we continue to develop this model with an updated suite of numerical simulations. We favor a scenario where Earth and Venus' progenitor nuclei form closer to the Sun and subsequently sculpt the Mercury-forming region by migrating towards their modern orbits. This rapid formation of ~0.5 Earth-mass cores at ~0.1-0.5 au is consistent with modern high-resolution simulations of planetesimal accretion. In successful realizations, Earth and Venus accrete mostly dry, Enstatite Chondrite-like material as they migrate; thus providing a simple explanation for the masses of all four terrestrial planets, inferred isotopic differences between Earth and Mars, and Mercury's isolated orbit. Furthermore, our models predict that Venus' composition should be similar to the Earth's, and possibly derived from a larger fraction of dry material. Conversely, Mercury analogs in our simulations attain a range of final compositions.
△ Less
Submitted 30 November, 2021;
originally announced December 2021.
-
Exploring tidal obliquity variations with SMERCURY-T
Authors:
Steven M. Kreyche,
Jason W. Barnes,
Billy L. Quarles,
John E. Chambers
Abstract:
We introduce our new code, SMERCURY-T, which is based on existing codes SMERCURY (Lissauer et al. 2012) and Mercury-T (Bolmont et al. 2015). The result is a mixed-variable symplectic N-body integrator that can compute the orbital and spin evolution of a planet within a multi-planet system under the influence of tidal spin torques from its star. We validate our implementation by comparing our exper…
▽ More
We introduce our new code, SMERCURY-T, which is based on existing codes SMERCURY (Lissauer et al. 2012) and Mercury-T (Bolmont et al. 2015). The result is a mixed-variable symplectic N-body integrator that can compute the orbital and spin evolution of a planet within a multi-planet system under the influence of tidal spin torques from its star. We validate our implementation by comparing our experimental results to that of a secular model. As we demonstrate in a series of experiments, SMERCURY-T allows for the study of secular spin-orbit resonance crossings and captures for planets within complex multi-planet systems. These processes can drive a planet's spin state to evolve along vastly different pathways on its road toward tidal equilibrium, as tidal spin torques dampen the planet's spin rate and evolve its obliquity. Additionally, we show the results of a scenario that exemplifies the crossing of a chaotic region that exists as the overlap of two spin-orbit resonances. The test planet experiences violent and chaotic swings in its obliquity until its eventual escape from resonance as it tidally evolves. All of these processes are and have been important over the obliquity evolution of many bodies within the Solar System and beyond, and have implications for planetary climate and habitability. SMERCURY-T is a powerful and versatile tool that allows for further study of these phenomena.
△ Less
Submitted 7 September, 2021;
originally announced September 2021.
-
The early instability scenario: Mars' mass explained by Jupiter's orbit
Authors:
Matthew S. Clement,
Nathan A. Kaib,
Sean N. Raymond,
John E. Chambers
Abstract:
The formation of the solar system's giant planets predated the ultimate epoch of massive impacts that concluded the process of terrestrial planet formation. Following their formation, the giant planets' orbits evolved through an episode of dynamical instability. Several qualities of the solar system have recently been interpreted as evidence of this event transpiring within the first ~100 Myr afte…
▽ More
The formation of the solar system's giant planets predated the ultimate epoch of massive impacts that concluded the process of terrestrial planet formation. Following their formation, the giant planets' orbits evolved through an episode of dynamical instability. Several qualities of the solar system have recently been interpreted as evidence of this event transpiring within the first ~100 Myr after the Sun's birth; around the same time as the final assembly of the inner planets. In a series of recent papers we argued that such an early instability could resolve several problems revealed in classic numerical studies of terrestrial planet formation; namely the small masses of Mars and the asteroid belt. In this paper, we revisit the early instability scenario with a large suite of simulations specifically designed to understand the degree to which Earth and Mars' formation are sensitive to the specific evolution of Jupiter and Saturn's orbits. By deriving our initial terrestrial disks directly from recent high-resolution simulations of planetesimal accretion, our results largely confirm our previous findings regarding the instability's efficiency of truncating the terrestrial disk outside of the Earth-forming region in simulations that best replicate the outer solar system. Moreover, our work validates the primordial 2:1 Jupiter-Saturn resonance within the early instability framework as a viable evolutionary path for the solar system. While our simulations elucidate the fragility of the terrestrial system during the epoch of giant planet migration, many realizations yield outstanding solar system analogs when scrutinized against a number of observational constraints. Finally, we highlight the inability of models to form adequate Mercury-analogs and the low eccentricities of Earth and Venus as the most significant outstanding problems for future numerical studies to resolve.
△ Less
Submitted 9 June, 2021;
originally announced June 2021.
-
Born extra-eccentric: A broad spectrum of primordial configurations of the gas giants that match their present-day orbits
Authors:
Matthew S. Clement,
Rogerio Deienno,
Nathan A. Kaib,
Andre Izidoro,
Sean N. Raymond,
John E. Chambers
Abstract:
In a recent paper we proposed that the giant planets' primordial orbits may have been eccentric (~0.05), and used a suite of dynamical simulations to show outcomes of the giant planet instability that are consistent with their present-day orbits. In this follow-up investigation, we present more comprehensive simulations incorporating superior particle resolution, longer integration times, and elim…
▽ More
In a recent paper we proposed that the giant planets' primordial orbits may have been eccentric (~0.05), and used a suite of dynamical simulations to show outcomes of the giant planet instability that are consistent with their present-day orbits. In this follow-up investigation, we present more comprehensive simulations incorporating superior particle resolution, longer integration times, and eliminating our prior means of artificially forcing instabilities to occur at specified times by shifting a planets' position in its orbit. While we find that the residual phase of planetary migration only minimally alters the the planets' ultimate eccentricities, our work uncovers several intriguing outcomes in realizations where Jupiter and Saturn are born with extremely large eccentricities (~0.10 and ~0.25, respectively). In successful simulations, the planets' orbits damp through interactions with the planetesimal disk prior to the instability, thus loosely replicating the initial conditions considered in our previous work. Our results therefore suggest an even wider range of plausible evolutionary pathways are capable of replicating Jupiter and Saturn's modern orbital architecture.
△ Less
Submitted 23 May, 2021;
originally announced May 2021.
-
Dynamical avenues for Mercury's origin II: in-situ formation in the inner terrestrial disk
Authors:
Matthew S. Clement,
John E. Chambers
Abstract:
Modern terrestrial planet formation models are highly successful at consistently generating planets with masses and orbits analogous to those of Earth and Venus. In stark contrast to classic theoretical predictions and inferred demographics of multi-planet systems of rocky exoplanets, the mass (>10) and orbital period (>2) ratios between Venus and Earth and the neighboring Mercury and Mars are not…
▽ More
Modern terrestrial planet formation models are highly successful at consistently generating planets with masses and orbits analogous to those of Earth and Venus. In stark contrast to classic theoretical predictions and inferred demographics of multi-planet systems of rocky exoplanets, the mass (>10) and orbital period (>2) ratios between Venus and Earth and the neighboring Mercury and Mars are not common outcomes in numerically generated systems. While viable solutions to the small-Mars problem are abundant in the literature, Mercury's peculiar origin remains rather mysterious. In this paper, we investigate the possibility that Mercury formed in a mass-depleted, inner region of the terrestrial disk (a < 0.5 au). This regime is often neglected in terrestrial planet formation models because of the high computational cost of resolving hundreds of short-period objects over ~100 Myr timescales. By testing multiple disk profiles and mass distributions, we identify several promising sets of initial conditions that lead to remarkably successful analog systems. In particular, our most successful simulations consider moderate total masses of Mercury-forming material (0.1-0.25 Earth masses). While larger initial masses tend to yield disproportionate Mercury analogs, smaller values often inhibit the planets' formation as the entire region of material is easily accreted by Venus. Additionally, we find that shallow surface density profiles and larger inventories of small planetesimals moderately improve the likelihood of adequately reproducing Mercury.
△ Less
Submitted 22 April, 2021;
originally announced April 2021.
-
Dynamical avenues for Mercury's origin I: The lone survivor of a primordial generation of short-period proto-planets
Authors:
Matthew S. Clement,
John E. Chambers,
Alan P. Jackson
Abstract:
The absence of planets interior to Mercury continues to puzzle terrestrial planet formation models, particularly when contrasted with the relatively high derived occurrence rates of short-period planets around Sun-like stars. Recent work proposed that the majority of systems hosting hot super-Earths attain their orbital architectures through an epoch of dynamical instability after forming in quasi…
▽ More
The absence of planets interior to Mercury continues to puzzle terrestrial planet formation models, particularly when contrasted with the relatively high derived occurrence rates of short-period planets around Sun-like stars. Recent work proposed that the majority of systems hosting hot super-Earths attain their orbital architectures through an epoch of dynamical instability after forming in quasi-stable, tightly packed configurations. Isotopic evidence seems to suggest that the formation of objects in the super-Earth mass regime is unlikely to have occurred in the solar system as the terrestrial-forming disk is thought to have been significantly mass-deprived starting around 2 Myr after CAI; a consequence of either Jupiter's growth or an intrinsic disk feature. Nevertheless, terrestrial planet formation models and high-resolution investigations of planetesimal dynamics in the gas disk phase occasionally find that quasi-stable proto-planets with masses comparable to that of Mars emerge in the vicinity of Mercury's modern orbit. In this paper, we investigate whether it is possible for a primordial configuration of such objects to be cataclysmically destroyed in a manner that leaves Mercury behind as the sole survivor without disturbing the other terrestrial worlds. We use numerical simulations to show that this scenario is plausible. In many cases, the surviving Mercury analog experiences a series of erosive impacts; thereby boosting its Fe/Si ratio. A caveat of our proposed genesis scenario for Mercury is that Venus typically experiences at least one late giant impact.
△ Less
Submitted 22 April, 2021;
originally announced April 2021.
-
Born eccentric: constraints on Jupiter and Saturn's pre-instability orbits
Authors:
Mattthew S. Clement,
Sean N. Raymond,
Nathan A. Kaib,
Rogerio Deienno,
John E. Chambers,
Andre Izidoro
Abstract:
An episode of dynamical instability is thought to have sculpted the orbital structure of the outer solar system. When modeling this instability, a key constraint comes from Jupiter's fifth eccentric mode (quantified by its amplitude M55), which is an important driver of the solar system's secular evolution. Starting from commonly-assumed near-circular orbits, the present-day giant planets' archite…
▽ More
An episode of dynamical instability is thought to have sculpted the orbital structure of the outer solar system. When modeling this instability, a key constraint comes from Jupiter's fifth eccentric mode (quantified by its amplitude M55), which is an important driver of the solar system's secular evolution. Starting from commonly-assumed near-circular orbits, the present-day giant planets' architecture lies at the limit of numerically generated systems, and M55 is rarely excited to its true value. Here we perform a dynamical analysis of a large batch of artificially triggered instabilities, and test a variety of configurations for the giant planets' primordial orbits. In addition to more standard setups, and motivated by the results of modern hydrodynamical simulations of the giant planets' evolution within the primordial gaseous disk, we consider the possibility that Jupiter and Saturn emerged from the nebular gas locked in 2:1 resonance with non-zero eccentricities. We show that, in such a scenario, the modern Jupiter-Saturn system represents a typical simulation outcome, and M55 is commonly matched. Furthermore, we show that Uranus and Neptune's final orbits are determined by a combination of the mass in the primordial Kuiper belt and that of an ejected ice giant.
△ Less
Submitted 23 September, 2020;
originally announced September 2020.
-
Embryo formation with GPU acceleration: reevaluating the initial conditions for terrestrial accretion
Authors:
Matthew S. Clement,
Nathan A. Kaib,
John E. Chambers
Abstract:
The solar system's terrestrial planets are thought to have accreted over millions of years out of a sea of smaller embryos and planetesimals. Because it is impossible to know the surface density profile for solids and size frequency distribution in the primordial solar nebula, distinguishing between the various proposed evolutionary schemes has been historically difficult. Nearly all previous simu…
▽ More
The solar system's terrestrial planets are thought to have accreted over millions of years out of a sea of smaller embryos and planetesimals. Because it is impossible to know the surface density profile for solids and size frequency distribution in the primordial solar nebula, distinguishing between the various proposed evolutionary schemes has been historically difficult. Nearly all previous simulations of terrestrial planet formation assume that Moon to Mars massed embryos formed throughout the inner solar system during the primordial gas-disk phase. However, validating this assumption through models of embryo accretion is computationally challenging because of the large number of bodies required. Here, we reevaluate this problem with GPU-accelerated, direct N-body simulations of embryo growth starting from r~100 km planetesimals. We find that embryos emerging from the primordial gas phase at a given radial distance already have masses similar to the largest objects at the same semi-major axis in the modern solar system. Thus, Earth and Venus attain ~50% of their modern mass, Mars-massed embryos form in the Mars region, and Ceres-massed objects are prevalent throughout asteroid belt. Consistent with other recent work, our new initial conditions for terrestrial accretion models produce markedly improved solar system analogs when evolved through the giant impact phase of planet formation. However, we still conclude that an additional dynamical mechanism such as giant planet migration is required to prevent Earth-massed Mars analogs from growing.
△ Less
Submitted 7 May, 2020;
originally announced May 2020.
-
Retrograde-rotating exoplanets experience obliquity excitations in an eccentricity-enabled resonance
Authors:
Steven M. Kreyche,
Jason W. Barnes,
Billy L. Quarles,
Jack J. Lissauer,
John E. Chambers,
Matthew M. Hedman
Abstract:
Previous studies have shown that planets that rotate retrograde (backwards with respect to their orbital motion) generally experience less severe obliquity variations than those that rotate prograde (the same direction as their orbital motion). Here we examine retrograde-rotating planets on eccentric orbits and find a previously unknown secular spin-orbit resonance that can drive significant obliq…
▽ More
Previous studies have shown that planets that rotate retrograde (backwards with respect to their orbital motion) generally experience less severe obliquity variations than those that rotate prograde (the same direction as their orbital motion). Here we examine retrograde-rotating planets on eccentric orbits and find a previously unknown secular spin-orbit resonance that can drive significant obliquity variations. This resonance occurs when the frequency of the planet's rotation axis precession becomes commensurate with an orbital eigenfrequency of the planetary system. The planet's eccentricity enables a participating orbital frequency through an interaction in which the apsidal precession of the planet's orbit causes a cyclic nutation of the planet's orbital angular momentum vector. The resulting orbital frequency follows the relationship $f = 2 \dot{\varpi} - \dotΩ$, where $\dot{\varpi}$ and $\dotΩ$ are the rates of the planet's changing longitude of periapsis and ascending node, respectively. We test this mechanism by simulating cases of a simple Earth-Jupiter system, and confirm the predicted resonance. Over the course of 100 Myr, the test Earths with rotation axis precession rates near the predicted resonant frequency experienced pronounced obliquity variations of order $10^\circ$-$30^\circ$. These variations can be significant, and suggest that while retrograde rotation is a stabilizing influence most of the time, retrograde rotators can experience large obliquity variations if they are on eccentric orbits and enter this spin-orbit resonance.
△ Less
Submitted 7 April, 2020; v1 submitted 30 March, 2020;
originally announced March 2020.
-
Dynamical Constraints on Mercury's Collisional Origin
Authors:
Matthew S. Clement,
Nathan A. Kaib,
John E. Chambers
Abstract:
Of the solar system's four terrestrial planets, the origin of Mercury is perhaps the most mysterious. Modern numerical simulations designed to model the dynamics of terrestrial planet formation systematically fail to replicate Mercury; which possesses just 5% the mass of Earth and the highest orbital eccentricity and inclination among the planets. However, Mercury's large iron-rich core and low vo…
▽ More
Of the solar system's four terrestrial planets, the origin of Mercury is perhaps the most mysterious. Modern numerical simulations designed to model the dynamics of terrestrial planet formation systematically fail to replicate Mercury; which possesses just 5% the mass of Earth and the highest orbital eccentricity and inclination among the planets. However, Mercury's large iron-rich core and low volatile inventory stand out among the inner planets, and seem to imply a violent collisional origin. Because most algorithms used for simulating terrestrial accretion do not consider the effects of collisional fragmentation, it has been difficult to test these collisional hypotheses within the larger context of planet formation. Here, we analyze a large suite of terrestrial accretion models that account for the fragmentation of colliding bodies. We find that planets with core mass fractions boosted as a result of repeated hit-and-run collisions are produced in 90% of our simulations. While many of these planets are similar to Mercury in mass, they rarely lie on Mercury-like orbits. Furthermore, we perform an additional batch of simulations designed to specifically test the single giant impact origin scenario. We find less than a 1% probability of simultaneously replicating the Mercury-Venus dynamical spacing and the terrestrial system's degree of orbital excitation after such an event. While dynamical models have made great strides in understanding Mars' low mass, their inability to form accurate Mercury analogs remains a glaring problem.
△ Less
Submitted 3 April, 2019;
originally announced April 2019.
-
The early instability scenario: terrestrial planet formation during the giant planet instability, and the effect of collisional fragmentation
Authors:
Matthew S. Clement,
Nathan A. Kaib,
Sean N. Raymond,
John E. Chambers,
Kevin J. Walsh
Abstract:
The solar system's dynamical state can be explained by an orbital instability among the giant planets. A recent model has proposed that the giant planet instability happened during terrestrial planet formation. This scenario has been shown to match the inner solar system by stunting Mars' growth and preventing planet formation in the asteroid belt. Here we present a large sample of new simulations…
▽ More
The solar system's dynamical state can be explained by an orbital instability among the giant planets. A recent model has proposed that the giant planet instability happened during terrestrial planet formation. This scenario has been shown to match the inner solar system by stunting Mars' growth and preventing planet formation in the asteroid belt. Here we present a large sample of new simulations of the "Early Instability" scenario. We use an N-body integration scheme that accounts for collisional fragmentation, and also perform a large set of control simulations that do not include an early giant planet instability. Since the total particle number decreases slower when collisional fragmentation is accounted for, the growing planets' orbits are damped more strongly via dynamical friction and encounters with small bodies that dissipate angular momentum (eg: hit-and-run impacts). Compared with simulations without collisional fragmentation, our fully evolved systems provide better matches to the solar system's terrestrial planets in terms of their compact mass distribution and dynamically cold orbits. Collisional processes also tend to lengthen the dynamical accretion timescales of Earth analogs, and shorten those of Mars analogs. This yields systems with relative growth timescales more consistent with those inferred from isotopic dating. Accounting for fragmentation is thus supremely important for any successful evolutionary model of the inner solar system.
△ Less
Submitted 18 December, 2018;
originally announced December 2018.
-
Symplectic Integrators: T + V Revisited and Round-Off Reduced
Authors:
John E Chambers
Abstract:
Symplectic integrators separate a problem into parts that can be solved in isolation, alternately advancing these sub-problems to approximate the evolution of the complete system. Problems with a single, dominant mass can use mixed-variable symplectic (MVS) integrators that separate the problem into Keplerian motion of satellites about the primary, and satellite-satellite interactions. Here, we ex…
▽ More
Symplectic integrators separate a problem into parts that can be solved in isolation, alternately advancing these sub-problems to approximate the evolution of the complete system. Problems with a single, dominant mass can use mixed-variable symplectic (MVS) integrators that separate the problem into Keplerian motion of satellites about the primary, and satellite-satellite interactions. Here, we examine T+V algorithms where the problem is separated into kinetic T and potential energy V terms. T+V integrators are typically less efficient than MVS algorithms. This difference is reduced by using different step sizes for primary-satellite and satellite-satellite interactions. The T+V method is improved further using 4th and 6th-order algorithms that include force gradients and symplectic correctors. We describe three 6th-order algorithms, containing 2 or 3 force evaluations per step, that are competitive with MVS in some cases. Round-off errors for T+V integrators can be reduced by several orders of magnitude, at almost no computational cost, using a simple modification that keeps track of accumulated changes in the coordinates and momenta. This makes T+V algorithms desirable for long-term, high-accuracy calculations.
△ Less
Submitted 17 December, 2018;
originally announced December 2018.
-
Pebble Accretion and the Diversity of Planetary Systems
Authors:
J. E. Chambers
Abstract:
I examine the standard model of planet formation, including pebble accretion, using numerical simulations. Planetary embryos large enough to become giant planets do not form beyond the ice line within a typical disk lifetime unless icy pebbles stick at higher speeds than in experiments using rocky pebbles. Systems like the Solar System (small inner planets, giant outer planets) can form if (i) icy…
▽ More
I examine the standard model of planet formation, including pebble accretion, using numerical simulations. Planetary embryos large enough to become giant planets do not form beyond the ice line within a typical disk lifetime unless icy pebbles stick at higher speeds than in experiments using rocky pebbles. Systems like the Solar System (small inner planets, giant outer planets) can form if (i) icy pebbles are stickier than rocky pebbles, and (ii) the planetesimal formation efficiency increases with pebble size, which prevents the formation of massive terrestrial planets. Growth beyond the ice line is dominated by pebble accretion. Most growth occurs early, when the surface density of pebbles is high due to inward drift of pebbles from the outer disk. Growth is much slower after the outer disk is depleted. The outcome is sensitive to the disk radius and turbulence level, which control the lifetime and maximum size of pebbles. The outcome is sensitive to the size of the largest planetesimals since there is a threshold mass for the onset of pebble accretion. The planetesimal formation rate is unimportant provided that some large planetesimals form while pebbles remain abundant. Two outcomes are seen, depending on whether pebble accretion begins while pebbles are still abundant. Either (i) multiple gas giant planets form beyond the ice line, small planets form close to the star, and a Kuiper-belt-like disk of bodies is scattered outwards by the giant planets; or (ii) no giants form and bodies remain Earth-mass or smaller.
△ Less
Submitted 21 April, 2016;
originally announced April 2016.
-
The Frequency of Giant impacts on Earth-like Worlds
Authors:
Elisa V. Quintana,
Thomas Barclay,
William Borucki,
Jason F. Rowe,
John E. Chambers
Abstract:
The late stages of terrestrial planet formation are dominated by giant impacts that collectively influence the growth, composition and habitability of any planets that form. Hitherto, numerical models designed to explore these late stage collisions have been limited by assuming that all collisions lead to perfect accretion, and many of these studies lack the large number of realizations needed to…
▽ More
The late stages of terrestrial planet formation are dominated by giant impacts that collectively influence the growth, composition and habitability of any planets that form. Hitherto, numerical models designed to explore these late stage collisions have been limited by assuming that all collisions lead to perfect accretion, and many of these studies lack the large number of realizations needed to account for the chaotic nature of N-body systems. We improve on these limitations by performing 280 simulations of planet formation around a Sun-like star, half of which used an N-body algorithm that has recently been modified to include fragmentation and hit-and-run (bouncing) collisions. We find that when fragmentation is included, the final planets formed are comparable in terms of mass and number, however their collision histories differ significantly and the accretion time approximately doubles. We explored impacts onto Earth-like planets which we parameterized in terms of their specific impact energies. Only 15 of our 164 Earth-analogs experienced an impact that was energetic enough to strip an entire atmosphere. To strip about half of an atmosphere requires energies comparable to recent models of the Moon-forming giant impact. Almost all Earth-analogs received at least one impact that met this criteria during the 2 Gyr simulations and the median was three giant impacts. The median time of the final giant impact was 43 Myr after the start of the simulations, leading us to conclude that the time-frame of the Moon-forming impact is typical amongst planetary systems around Sun-like stars.
△ Less
Submitted 10 March, 2016; v1 submitted 11 November, 2015;
originally announced November 2015.
-
The Fragility of the Terrestrial Planets During a Giant Planet Instability
Authors:
Nathan A. Kaib,
John E. Chambers
Abstract:
Many features of the outer solar system are replicated in numerical simulations if the giant planets undergo an orbital instability that ejects one or more ice giants. During this instability, Jupiter and Saturn's orbits diverge, crossing their 2:1 mean motion resonance (MMR), and this resonance-crossing can excite the terrestrial planet orbits. Using a large ensemble of simulations of this giant…
▽ More
Many features of the outer solar system are replicated in numerical simulations if the giant planets undergo an orbital instability that ejects one or more ice giants. During this instability, Jupiter and Saturn's orbits diverge, crossing their 2:1 mean motion resonance (MMR), and this resonance-crossing can excite the terrestrial planet orbits. Using a large ensemble of simulations of this giant planet instability, we directly model the evolution of the terrestrial planet orbits during this process, paying special attention to systems that reproduce the basic features of the outer planets. In systems that retain four giant planets and finish with Jupiter and Saturn beyond their 2:1 MMR, we find at least an 85% probability that at least one terrestrial planet is lost. Moreover, systems that manage to retain all four terrestrial planets often finish with terrestrial planet eccentricities and inclinations larger than the observed ones. There is less than a ~5% chance that the terrestrial planet orbits will have a level of excitation comparable to the observed orbits. If we factor in the probability that the outer planetary orbits are well-replicated, we find a probability of 1% or less that the orbital architectures of the inner and outer planets are simultaneously reproduced in the same system. These small probabilities raise the prospect that the giant planet instability occurred before the terrestrial planets had formed. This scenario implies that the giant planet instability is not the source of the Late Heavy Bombardment and that terrestrial planet formation finished with the giant planets in their modern configuration.
△ Less
Submitted 28 October, 2015;
originally announced October 2015.
-
Five New Transit Epochs of the Exoplanet OGLE-TR-111b
Authors:
Sergio Hoyer,
Patricio Rojo,
Mercedes López-Morales,
Rodrigo Díaz,
John E. Chambers,
Dante Minniti
Abstract:
We report five new transit epochs of the extrasolar planet OGLE-TR-111b, observed in the v-HIGH and Bessell I bands with the FORS1 and FORS2 at the ESO Very Large Telescope, between April and May 2008. The new transits have been combined with all previously published transit data for this planet to provide a new Transit Timing Variations (TTVs) analysis of its orbit. We discard TTVs with amplitude…
▽ More
We report five new transit epochs of the extrasolar planet OGLE-TR-111b, observed in the v-HIGH and Bessell I bands with the FORS1 and FORS2 at the ESO Very Large Telescope, between April and May 2008. The new transits have been combined with all previously published transit data for this planet to provide a new Transit Timing Variations (TTVs) analysis of its orbit. We discard TTVs with amplitudes larger than 1.5 minutes over a 4-year observation time baseline, in agreement with the recent result by Adams et al.(2010a). Dynamical simulations fully exclude the presence of additional planets in the system with masses greater than 1.3, 0.4 and 0.5 M_earth at the 3:2, 1:2, 2:1 resonances, respectively. We also place an upper limit of about 30 M_earth on the mass of potential second planets in the region between the 3:2 and 1:2 mean-motion resonances.
△ Less
Submitted 16 February, 2011;
originally announced February 2011.
-
Five Long-period Extrasolar Planets in Eccentric orbits from the Magellan Planet Search Program
Authors:
Pamela Arriagada,
R. Paul Butler,
Dante Minniti,
Mercedes Lopez-Morales,
Stephen A. Shectman,
Fred C. Adams,
Alan P. Boss,
John E. Chambers
Abstract:
Five new planets orbiting G and K dwarfs have emerged from the Magellan velocity survey. These companions are jovian-mass planets in eccentric (e \geq 0.24) intermediate and long-period orbits. HD 86226b orbits a solar metallicity G2 dwarf. The MP sin i mass of the planet is 1.5 MJUP, the semi-major axis is 2.6 AU, and the eccentricity 0.73. HD 129445b orbits a metal rich G6 dwarf. The minimum m…
▽ More
Five new planets orbiting G and K dwarfs have emerged from the Magellan velocity survey. These companions are jovian-mass planets in eccentric (e \geq 0.24) intermediate and long-period orbits. HD 86226b orbits a solar metallicity G2 dwarf. The MP sin i mass of the planet is 1.5 MJUP, the semi-major axis is 2.6 AU, and the eccentricity 0.73. HD 129445b orbits a metal rich G6 dwarf. The minimum mass of the planet is MP sin i =1.6 MJUP, the semi-major axis is 2.9 AU, and the eccentricity 0.70. HD 164604b orbits a K2 dwarf. The MP sin i mass is 2.7 MJUP, semi-major axis is 1.3 AU, and the eccentricity is 0.24. HD 175167b orbits a metal rich G5 star. The MP sin i mass is 7.8 MJUP, the semi-major axis is 2.4 AU, and the eccentricity 0.54. HD 152079b orbits a G6 dwarf. The MP sin i mass of the planet is 3 MJUP, the semi-major axis is 3.2 AU, and the eccentricity is 0.60.
△ Less
Submitted 26 January, 2010; v1 submitted 22 January, 2010;
originally announced January 2010.
-
Low Mass Companions for Five Solar-Type Stars from the Magellan Planet Search Program
Authors:
Dante Minniti,
R. Paul Butler,
Mercedes Lopez-Morales,
Stephen A. Shectman,
Fred C. Adams,
Pamela Arriagada,
Alan P. Boss,
John E. Chambers
Abstract:
We report low mass companions orbiting five Solar-type stars that have emerged from the Magellan precision Doppler velocity survey, with minimum (Msini) masses ranging from 1.2 to 25 Mjup. These nearby target stars range from mildly metal-poor to metal-rich, and appear to have low chromospheric activity. The companions to the brightest two of these stars have previously been reported from the CO…
▽ More
We report low mass companions orbiting five Solar-type stars that have emerged from the Magellan precision Doppler velocity survey, with minimum (Msini) masses ranging from 1.2 to 25 Mjup. These nearby target stars range from mildly metal-poor to metal-rich, and appear to have low chromospheric activity. The companions to the brightest two of these stars have previously been reported from the CORALIE survey. Four of these companions (HD 48265-b, HD 143361-b, HD 28185-b, HD 111232-b) are low-mass Jupiter-like planets in eccentric intermediate and long-period orbits. On the other hand, the companion to HD 43848 appears to be a long period brown dwarf in a very eccentric orbit.
△ Less
Submitted 29 October, 2008;
originally announced October 2008.
-
How eccentric orbital solutions can hide planetary systems in 2:1 resonant orbits
Authors:
Guillem Anglada-Escude,
Mercedes Lopez-Morales,
John E. Chambers
Abstract:
The Doppler technique measures the reflex radial motion of a star induced by the presence of companions and is the most successful method to detect exoplanets. If several planets are present, their signals will appear combined in the radial motion of the star, leading to potential misinterpretations of the data. Specifically, two planets in 2:1 resonant orbits can mimic the signal of a single pl…
▽ More
The Doppler technique measures the reflex radial motion of a star induced by the presence of companions and is the most successful method to detect exoplanets. If several planets are present, their signals will appear combined in the radial motion of the star, leading to potential misinterpretations of the data. Specifically, two planets in 2:1 resonant orbits can mimic the signal of a single planet in an eccentric orbit. We quantify the implications of this statistical degeneracy for a representative sample of the reported single exoplanets with available datasets, finding that 1) around 35 percent of the published eccentric one-planet solutions are statistically indistinguishable from planetary systems in 2:1 orbital resonance, 2) another 40 percent cannot be statistically distinguished from a circular orbital solution and 3) planets with masses comparable to Earth could be hidden in known orbital solutions of eccentric super-Earths and Neptune mass planets.
△ Less
Submitted 25 November, 2009; v1 submitted 8 September, 2008;
originally announced September 2008.
-
Is our Sun a Singleton?
Authors:
D. Malmberg,
M. B. Davies,
J. E. Chambers,
F. De Angeli,
R. P. Church,
D. Mackey,
M. I. Wilkinson
Abstract:
Most stars are formed in a cluster or association, where the number density of stars can be high. This means that a large fraction of initially-single stars will undergo close encounters with other stars and/or exchange into binaries. We describe how such close encounters and exchange encounters can affect the properties of a planetary system around a single star. We define a singleton as a sing…
▽ More
Most stars are formed in a cluster or association, where the number density of stars can be high. This means that a large fraction of initially-single stars will undergo close encounters with other stars and/or exchange into binaries. We describe how such close encounters and exchange encounters can affect the properties of a planetary system around a single star. We define a singleton as a single star which has never suffered close encounters with other stars or spent time within a binary system. It may be that planetary systems similar to our own solar system can only survive around singletons. Close encounters or the presence of a stellar companion will perturb the planetary system, often leaving planets on tighter and more eccentric orbits. Thus planetary systems which initially resembled our own solar system may later more closely resemble some of the observed exoplanet systems.
△ Less
Submitted 19 September, 2007;
originally announced September 2007.
-
N-body integrators for planets in binary star systems
Authors:
John E. Chambers
Abstract:
Symplectic integrators are the tool of choice for many researchers studying dynamical systems because of their good long-term energy conservation properties. For systems with a dominant central mass, symplectic integrators are also highly efficient. In this chapter, I describe the theory of symplectic integrators in terms of Lie series. I show how conventional symplectic algorithms have been ada…
▽ More
Symplectic integrators are the tool of choice for many researchers studying dynamical systems because of their good long-term energy conservation properties. For systems with a dominant central mass, symplectic integrators are also highly efficient. In this chapter, I describe the theory of symplectic integrators in terms of Lie series. I show how conventional symplectic algorithms have been adapted for use in binary-star systems to study problems such as the dynamical stability of multi-planet systems and the accretion of planets from planetesimals. This is achieved by devising new coordinate systems for the wide-binary and close-binary cases separately. I show how the performance of these algorithms can be improved at little extra cost using symplectic correctors. Finally, I discuss drawbacks of these algorithms, in particular in dealing with close encounters with one or both members of the binary, and the prospects for overcoming these problems.
△ Less
Submitted 22 May, 2007;
originally announced May 2007.
-
Terrestrial Planet Formation Around Individual Stars Within Binary Star Systems
Authors:
Elisa V. Quintana,
Fred C. Adams,
Jack J. Lissauer,
John E. Chambers
Abstract:
We calculate herein the late stages of terrestrial planet accumulation around a solar type star that has a binary companion with semimajor axis larger than the terrestrial planet region. We perform more than one hundred simulations to survey binary parameter space and to account for sensitive dependence on initial conditions in these dynamical systems. As expected, sufficiently wide binaries lea…
▽ More
We calculate herein the late stages of terrestrial planet accumulation around a solar type star that has a binary companion with semimajor axis larger than the terrestrial planet region. We perform more than one hundred simulations to survey binary parameter space and to account for sensitive dependence on initial conditions in these dynamical systems. As expected, sufficiently wide binaries leave the planet formation process largely unaffected. As a rough approximation, binary stars with periastron $q_B > 10$ AU have minimal effect on terrestrial planet formation within $\sim 2$ AU of the primary, whereas binary stars with $q_B \la$ 5 AU restrict terrestrial planet formation to within $\sim$ 1 AU of the primary star. Given the observed distribution of binary orbital elements for solar type primaries, we estimate that about 40 -- 50 percent of the binary population is wide enough to allow terrestrial planet formation to take place unimpeded. The large number of simulations allows for us to determine the distribution of results -- the distribution of plausible terrestrial planet systems -- for effectively equivalent starting conditions. We present (rough) distributions for the number of planets, their masses, and their orbital elements.
△ Less
Submitted 9 January, 2007;
originally announced January 2007.
-
The instability of planetary systems in binaries: how the Kozai mechanism leads to strong planet-planet interactions
Authors:
Daniel Malmberg,
Melvyn B. Davies,
John E. Chambers
Abstract:
In this letter we consider the evolution of a planetary system around a star inside a wide binary. We simulate numerically the evolution of the planetary orbits for both co-planar and highly-inclined systems. We find that the Kozai mechanism operates in the latter case. This produces a highly eccentric outer planet whose orbit crosses those of some of the inner planets. Strong planet-planet inte…
▽ More
In this letter we consider the evolution of a planetary system around a star inside a wide binary. We simulate numerically the evolution of the planetary orbits for both co-planar and highly-inclined systems. We find that the Kozai mechanism operates in the latter case. This produces a highly eccentric outer planet whose orbit crosses those of some of the inner planets. Strong planet-planet interactions then follow resulting in the ejection of one or more planets. We note that planetary systems resembling our solar system, formed around single stars in stellar clusters may exchange into binaries and thus will be vulnerable to planet stripping. This process will reduce the number of solar-system like planetary systems, and may produce at least some of the observed extra-solar planets.
△ Less
Submitted 15 January, 2007; v1 submitted 1 December, 2006;
originally announced December 2006.
-
Planet Formation with Migration
Authors:
J. E. Chambers
Abstract:
In the core-accretion model, gas-giant planets form solid cores which then accrete gaseous envelopes. Tidal interactions with disk gas cause a core to undergo inward type-I migration in 10^4 to 10^5 years. Cores must form faster than this to survive. Giant planets clear a gap in the disk and undergo inward type-II migration in <10^6 years if observed disk accretion rates apply to the disk as a w…
▽ More
In the core-accretion model, gas-giant planets form solid cores which then accrete gaseous envelopes. Tidal interactions with disk gas cause a core to undergo inward type-I migration in 10^4 to 10^5 years. Cores must form faster than this to survive. Giant planets clear a gap in the disk and undergo inward type-II migration in <10^6 years if observed disk accretion rates apply to the disk as a whole. Type-II migration times exceed typical disk lifetimes if viscous accretion occurs mainly in the surface layers of disks. Low turbulent viscosities near the midplane may allow planetesimals to form by coagulation of dust grains. The radius r of such planetesimals is unknown. If r<0.5 km, the core formation time is shorter than the type-I migration timescale and cores will survive. Migration is substantial in most cases, leading to a wide range of planetary orbits, consistent with the observed variety of extrasolar systems. When r is of order 100m and midplane alpha is of order 3 times 10^-5, giant planets similar to those in the Solar System can form.
△ Less
Submitted 30 October, 2006;
originally announced October 2006.
-
Short-Term Dynamical Interactions Among Extrasolar Planets
Authors:
Gregory Laughlin,
John E. Chambers
Abstract:
We show that short-term perturbations among massive planets in multiple planet systems can result in radial velocity variations of the central star which differ substantially from velocity variations derived assuming the planets are executing independent Keplerian motions. We discuss two alternate fitting methods which can lead to an improved dynamical description of multiple planet systems. In…
▽ More
We show that short-term perturbations among massive planets in multiple planet systems can result in radial velocity variations of the central star which differ substantially from velocity variations derived assuming the planets are executing independent Keplerian motions. We discuss two alternate fitting methods which can lead to an improved dynamical description of multiple planet systems. In the first method, the osculating orbital elements are determined via a Levenberg-Marquardt minimization scheme driving an N-body integrator. The second method is an improved analytic model in which orbital elements are allowed to vary according to a simple model for resonant interactions between the planets. Both of these methods can determine the true masses for the planets by eliminating the sin(i) degeneracy inherent in fits that assume independent Keplerian motions. We apply our fitting methods to the GJ876 radial velocity data (Marcy et al. 2001), and argue that the mass factors for the two planets are likely in the 1.25-2.0 range
△ Less
Submitted 27 March, 2001; v1 submitted 23 January, 2001;
originally announced January 2001.
-
Pseudo-High-Order Symplectic Integrators
Authors:
J. E. Chambers,
M. A. Murison
Abstract:
Symplectic N-body integrators are widely used to study problems in celestial mechanics. The most popular algorithms are of 2nd and 4th order, requiring 2 and 6 substeps per timestep, respectively. The number of substeps increases rapidly with order in timestep, rendering higher-order methods impractical. However, symplectic integrators are often applied to systems in which perturbations between…
▽ More
Symplectic N-body integrators are widely used to study problems in celestial mechanics. The most popular algorithms are of 2nd and 4th order, requiring 2 and 6 substeps per timestep, respectively. The number of substeps increases rapidly with order in timestep, rendering higher-order methods impractical. However, symplectic integrators are often applied to systems in which perturbations between bodies are a small factor of the force due to a dominant central mass. In this case, it is possible to create optimized symplectic algorithms that require fewer substeps per timestep. This is achieved by only considering error terms of order epsilon, and neglecting those of order epsilon^2, epsilon^3 etc. Here we devise symplectic algorithms with 4 and 6 substeps per step which effectively behave as 4th and 6th-order integrators when epsilon is small. These algorithms are more efficient than the usual 2nd and 4th-order methods when applied to planetary systems.
△ Less
Submitted 14 October, 1999;
originally announced October 1999.