Showing 1–50 of 67 results for author: Brasser, R

Forming the Trappist-1 system in two steps during the recession of the disc inner edge

Authors: Gabriele Pichierri, Alessandro Morbidelli, Konstantin Batygin, Ramon Brasser

Abstract: Trappist-1 hosts 7 planets where period ratios of neighbouring pairs are close to the 8:5, 5:3, 3:2, 3:2, 4:3, and 3:2 ratios in increasing distance from the star. The Laplace angles associated with neighbouring triplets are observed to be librating, proving the resonant nature of the system. This compact, resonant configuration is a manifest sign of disc-driven migration; however, the preferred o… ▽ More Trappist-1 hosts 7 planets where period ratios of neighbouring pairs are close to the 8:5, 5:3, 3:2, 3:2, 4:3, and 3:2 ratios in increasing distance from the star. The Laplace angles associated with neighbouring triplets are observed to be librating, proving the resonant nature of the system. This compact, resonant configuration is a manifest sign of disc-driven migration; however, the preferred outcome of such evolution is the establishment of first-order resonances, not the high-order resonances observed in the inner system. Here, we explain the observed orbital configuration in a model that is largely independent on the specific disc migration and orbital circularisation efficiencies. Together with migration, the two key elements of our model are: i) the inner border of the protoplanetary disc receded with time; and ii) the system was initially separated in two sub-systems. Specifically, the inner b, c d and e planets were initially placed in a 3:2 resonance chain and then evolved to the 8:5 -- 5:3 commensurability between planets b, c and d under the effect of the recession of the inner edge of the disc, whereas the outer planets migrated to the inner edge at a later time, establishing the remaining resonances. Our results pivot on the dynamical role of the presently unobservable recession of the inner edge of protoplanetary discs. They also reveal the role of recurring phases of convergent migration followed by resonant repulsion with associated orbital circularisation when resonant chains interact with migration barriers. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: To appear in Nature Astronomy

Giant planet formation in the solar system

Authors: Anuja Raorane, Ramon Brasser, Soko Matsumura, Tommy Chi Ho Lau, Man Hoi Lee, Audrey Bouvier

Abstract: The formation history of Jupiter has been of interest due to its ability to shape the solar system's history. Yet little attention has been paid to the formation and growth of Saturn and the other giant planets. Here, we explore the implications of the simplest disc and pebble accretion model with steady-state accretion on the formation of giant planets in the solar system through N-body simulatio… ▽ More The formation history of Jupiter has been of interest due to its ability to shape the solar system's history. Yet little attention has been paid to the formation and growth of Saturn and the other giant planets. Here, we explore the implications of the simplest disc and pebble accretion model with steady-state accretion on the formation of giant planets in the solar system through N-body simulations. We conducted a statistical survey of different disc parameters and initial conditions of the protoplanetary disc to establish which combination best reproduces the present outer solar system. We examined the effect of the initial planetesimal disc mass, the number of planetesimals and their size-frequency distribution slope, pebble accretion prescription, and sticking efficiency on the likelihood of forming gas giants and their orbital distribution. The results reveal that the accretion sticking efficiency is the most sensitive parameter for controlling the final masses and number of giant planets. We have been unable to replicate the formation of all three types of giant planets in the solar system in a single simulation. The probability distribution of the final location of the giant planets is approximately constant in $\log r$, suggesting there is a slight preference for formation closer to the Sun but no preference for more massive planets to form closer. The eccentricity distribution has a higher mean for more massive planets, indicating that systems with more massive planets are more violent. The formation timescales of the cores of the gas giants are distinct, suggesting that they formed sequentially. △ Less

Submitted 30 July, 2024; v1 submitted 26 March, 2024; originally announced March 2024.

Comments: 17 pages, 18 figures, 1 table

Can the giant planets of the Solar System form via pebble accretion in a smooth protoplanetary disc?

Authors: Tommy Chi Ho Lau, Man Hoi Lee, Ramon Brasser, Soko Matsumura

Abstract: Prevailing $N$-body planet formation models typically start with lunar-mass embryos and show a general trend of rapid migration of massive planetary cores to the inner Solar System in the absence of a migration trap. This setup cannot capture the evolution from a planetesimal to embryo, which is crucial to the final architecture of the system. We aim to model planet formation with planet migration… ▽ More Prevailing $N$-body planet formation models typically start with lunar-mass embryos and show a general trend of rapid migration of massive planetary cores to the inner Solar System in the absence of a migration trap. This setup cannot capture the evolution from a planetesimal to embryo, which is crucial to the final architecture of the system. We aim to model planet formation with planet migration starting with planetesimals of $\sim10^{-6}$ -- $10^{-4}M_\oplus$ and reproduce the giant planets of the Solar System. We simulated a population of 1,000 -- 5,000 planetesimals in a smooth protoplanetary disc, which was evolved under the effects of their mutual gravity, pebble accretion, gas accretion, and planet migration, employing the parallelized $N$-body code SyMBAp. We find that the dynamical interactions among growing planetesimals are vigorous and can halt pebble accretion for excited bodies. While a set of results without planet migration produces one to two gas giants and one to two ice giants beyond 6 au, massive planetary cores readily move to the inner Solar System once planet migration is in effect. Dynamical heating is important in a planetesimal disc and the reduced pebble encounter time should be considered in similar models. Planet migration remains a challenge to form cold giant planets in a smooth protoplanetary disc, which suggests an alternative mechanism is required to stop them at wide orbits. △ Less

Submitted 25 March, 2024; v1 submitted 10 January, 2024; originally announced January 2024.

Comments: 17 pages, 14 figures, replaced with published version on A&A

Speeding up the GENGA N-body integrator on consumer-grade graphics cards

Authors: R. Brasser, S. L. Grimm, P. Hatalova, J. G. Stadel

Abstract: GPU computing is popular due to the calculation potential of a single card. The N-body integrator GENGA is built to for this, but it suffers a performance penalty on consumer-grade GPUs due to their truncated double precision (FP64) performance. We aim to speed up GENGA on consumer-grade cards by harvesting their high single-precision performance (FP32). We modified GENGA to be able to compute the… ▽ More GPU computing is popular due to the calculation potential of a single card. The N-body integrator GENGA is built to for this, but it suffers a performance penalty on consumer-grade GPUs due to their truncated double precision (FP64) performance. We aim to speed up GENGA on consumer-grade cards by harvesting their high single-precision performance (FP32). We modified GENGA to be able to compute the long-distance forces between bodies in FP32 precision and tested this with 5 experiments. We ran simulations with similar initial conditions of 6600 planetesimals in both FP32 and FP64 precision. We also ran simulations that i) began with a mixture of planetesimals and planetary embryos, ii) planetesimal-driven giant planet migration, and iii) terrestrial planet formation with a gas disc. Second, we ran the same simulation beginning with 40 000 planetesimals using both FP32 and FP64 precision forces on a variety of consumer-grade and Tesla GPUs to measure the performance boost of FP32 computing. There are no statistical differences when running in FP32 or FP64 precision that can be attributed to the force prescription rather than stochastic effects. The uncertainties in energy are almost identical when using both precisions. However, the uncertainty in the angular momentum using FP32 rather than FP64 precision long-range forces is about two orders of magnitude greater, but still very low. Running the simulations in single precision on consumer-grade cards decreases running time by a factor of three and becomes within a factor of three of a Tesla A100 GPU. Additional tuning speeds up the simulation by a factor of two across all types of cards. The option to compute the long-range forces in single precision in GENGA when using consumer-grade GPUs dramatically improves performance at a little penalty to accuracy. There is an additional environmental benefit because it reduces energy usage. △ Less

Submitted 15 September, 2023; originally announced September 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

Saturn's ancient regular satellites

Authors: Emily. W. Wong, Ramon Brasser, Stephanie. C. Werner, Michelle. R. Kirchoff

Abstract: Are Saturn's regular satellites young or old? And how old are Enceladus' cratered plains? To answer these questions we computed model surface ages of the most heavily cratered terrains on Saturn's regular icy satellites using new high-resolution outer Solar System evolution simulations, and coupled with improved estimates of the trans-Neptunian objects populations. The output of the simulations al… ▽ More Are Saturn's regular satellites young or old? And how old are Enceladus' cratered plains? To answer these questions we computed model surface ages of the most heavily cratered terrains on Saturn's regular icy satellites using new high-resolution outer Solar System evolution simulations, and coupled with improved estimates of the trans-Neptunian objects populations. The output of the simulations allowed us to construct a model impact chronology onto Saturn which automatically applies to its regular satellites. We used crater densities and our impact chronology onto Saturn to construct model impact-crater isochrons, i.e., the scaling of the satellite crater production function with time. The surface ages derived for the cratered plains on Mimas, Enceladus, Tethys, Dione and Rhea range from 4.1 Ga to 4.4 Ga, with the surfaces of Mimas and Enceladus roughly 200 Myr younger than those of the outer three satellites. Uncertainties in these ages are less than 300 Myr. The calculated model surface ages of these satellites are consistent over as much as two orders of magnitude in the observed crater diameter. The similarity of the crater production function amongst all satellites suggests that they were bombarded by a single impactor source. This work supports the idea that Saturn's regular satellites are ancient, and has implications for their formation and their tidal evolution. △ Less

Submitted 7 September, 2023; originally announced September 2023.

Comments: Accepted for publication in Icarus. 30 pages, 11 figures; including Appendix and Appendix Figures

Forming rocky exoplanets around K-dwarf stars

Authors: P. Hatalova, R. Brasser, E. Mamonova, S. C. Werner

Abstract: How multiple close-in super-Earths form around stars with masses lower than that of the Sun is still an open issue. Several recent modeling studies have focused on planet formation around M-dwarf stars, but so far no studies have focused specifically on K dwarfs, which are of particular interest in the search for extraterrestrial life. We aim to reproduce the currently known population of close-in… ▽ More How multiple close-in super-Earths form around stars with masses lower than that of the Sun is still an open issue. Several recent modeling studies have focused on planet formation around M-dwarf stars, but so far no studies have focused specifically on K dwarfs, which are of particular interest in the search for extraterrestrial life. We aim to reproduce the currently known population of close-in super-Earths observed around K-dwarf stars and their system characteristics. We performed 48 high-resolution N-body simulations of planet formation via planetesimal accretion using the existing GENGA software running on GPUs. In the simulations we varied the initial disk mass and the solid and gas surface density profiles. Each simulation began with 12000 bodies with radii of between 200 and 2000 km around two different stars, with masses of 0.6 and 0.8 $M_{\odot}$. Most simulations ran for 20 Myr, with several simulations extended to 40 or 100 Myr. The mass distributions for the planets with masses between 2 and 12 $M_\oplus$ show a strong preference for planets with masses $M_p<6$ $M_\oplus$ and a lesser preference for planets with larger masses, whereas the mass distribution for the observed sample increases almost linearly. However, we managed to reproduce the main characteristics and architectures of the known planetary systems and produce mostly long-term angular-momentum-deficit-stable, nonresonant systems, but we require an initial disk mass of 15 $M_\oplus$ or higher and a gas surface density value at 1 AU of 1500 g cm$^{-2}$ or higher. Our simulations also produce many low-mass planets with $M<2$ $M_\oplus$, which are not yet found in the observed population, probably due to the observational biases. The final systems contain only a small number of planets, which could possibly accrete substantial amounts of gas, and these formed after the gas had mostly dissipated. △ Less

Submitted 11 August, 2023; originally announced August 2023.

Comments: 27 pages, 29 figures, to be published in Astronomy & Astrophysics

Journal ref: A&A 676, A131 (2023)

Terrestrial planet formation from a ring

Authors: J. M. Y. Woo, A. Morbidelli, S. L. Grimm, J. Stadel, R. Brasser

Abstract: It has been long proposed that, if all the terrestrial planets form within a tiny ring of solid material at around 1 AU, the concentrated mass-distance distribution of the current system can be reproduced. Recent planetesimal formation models also support this idea. In this study, we revisit the ring model by performing a number of high-resolution N-body simulations for 10 Myr of a ring of self-in… ▽ More It has been long proposed that, if all the terrestrial planets form within a tiny ring of solid material at around 1 AU, the concentrated mass-distance distribution of the current system can be reproduced. Recent planetesimal formation models also support this idea. In this study, we revisit the ring model by performing a number of high-resolution N-body simulations for 10 Myr of a ring of self-interacting planetesimals, with various radial distributions of the gas disc. We found that even if all the planetesimals form at ~1 AU in a minimum mass solar nebula-like disc, the system tends to spread radially as accretion proceeds, resulting in a system of planetary embryos lacking mass-concentration at ~1 AU. Modifying the surface density of the gas disc into a concave shape with a peak at ~1 AU helps to maintain mass concentrated at ~1 AU and solve the radial dispersion problem. We further propose that such a disc should be short lived (<= 1 Myr) and with a shallower radial gradient in the innermost region (< 1 AU) than previously proposed to prevent a too-rapid growth of Earth. Future studies should extend to ~100 Myr the most promising simulations and address in a self-consistent manner the evolution of the asteroid belt and its role in the formation of the terrestrial planets. △ Less

Submitted 27 February, 2023; originally announced February 2023.

Comments: Accepted by Icarus

Long-term tidal evolution of the TRAPPIST-1 system

Authors: R. Brasser, G. Pichierri, V. Dobos, A. C. Barr

Abstract: The ultracool M-dwarf star TRAPPIST-1 is surrounded by seven planets configured in a resonant chain. Transit-timing variations have shown that the planets are caught in multiple three-body resonances and that their orbits are slightly eccentric, probably caused by resonant forcing. The current values of the eccentricities could be a remnant from their formation. Here we run numerical simulations u… ▽ More The ultracool M-dwarf star TRAPPIST-1 is surrounded by seven planets configured in a resonant chain. Transit-timing variations have shown that the planets are caught in multiple three-body resonances and that their orbits are slightly eccentric, probably caused by resonant forcing. The current values of the eccentricities could be a remnant from their formation. Here we run numerical simulations using fictitious forces of trapping the fully-grown planets in resonances as they migrated in the gas disc, followed by numerical simulations detailing their tidal evolution. For a reduced disc scale height $h\sim 0.03$--0.05, the eccentricities of the planets upon capture in resonance are higher than their current values by factors of a few. We show that the current eccentricities and spacing of planets d to h are natural outcomes of coupled tidal evolution wherein the planets simultaneously damp their eccentricities and separate due to their resonant interaction. We further show that the planets evolve along a set of equilibrium curves in semimajor axis--eccentricity phase space that are defined by the resonances, and that conserve angular momentum. As such, the current 8:5--5:3--(3:2)$^2$--4:3--3:2 resonant configuration cannot be reproduced from a primordial (3:2)$^4$--4:3--3:2 resonant configuration from tidal dissipation in the planets alone. We use our simulations to constrain the long-term tidal parameters $k_2/Q$ for planets b to e, which are in the range $10^{-3}$ to $10^{-2}$, and show that these are mostly consistent with those obtained from interior modelling following reasonable assumptions. △ Less

Submitted 12 July, 2022; originally announced July 2022.

Comments: Accepted in Monthly Notices of the Royal Astronomical Society

Evidence of a primordial isotopic gradient in the inner region of the solar protoplanetary disc

Authors: J. Mah, R. Brasser, J. M. Y. Woo, A. Bouvier, S. J. Mojzsis

Abstract: Not only do the sampled terrestrial worlds (Earth, Mars, and asteroid 4 Vesta) differ in their mass-independent (nucleosynthetic) isotopic compositions of many elements (e.g. $\varepsilon^{48}$Ca, $\varepsilon^{50}$Ti, $\varepsilon^{54}$Cr, $\varepsilon^{92}$Mo), the magnitudes of some of these isotopic anomalies also appear to correlate with heliocentric distance. While the isotopic differences b… ▽ More Not only do the sampled terrestrial worlds (Earth, Mars, and asteroid 4 Vesta) differ in their mass-independent (nucleosynthetic) isotopic compositions of many elements (e.g. $\varepsilon^{48}$Ca, $\varepsilon^{50}$Ti, $\varepsilon^{54}$Cr, $\varepsilon^{92}$Mo), the magnitudes of some of these isotopic anomalies also appear to correlate with heliocentric distance. While the isotopic differences between the Earth and Mars may be readily accounted for by the accretion of mostly local materials in distinct regions of the protoplanetary disc, it is unclear whether this also applies to asteroid Vesta. Here we analysed the available data from our numerical simulation database to determine the formation location of Vesta in the framework of three planet-formation models: classical, Grand Tack, and Depleted Disc. We find that Vesta has a high probability of forming locally in the asteroid belt in models where material mixing in the inner disc is limited; this limited mixing is implied by the isotopic differences between the Earth and Mars. Based on our results, we propose several criteria to explain the apparent correlation between the different nucleosynthetic isotopic compositions of the Earth, Mars, and Vesta: (1) these planetary bodies accreted their building blocks in different regions of the disc, (2) the inner disc is characterised by an isotopic gradient, and (3) the isotopic gradient was preserved during the formation of these planetary bodies and was not diluted by material mixing in the disc (e.g. via giant planet migration). △ Less

Submitted 8 March, 2022; v1 submitted 22 February, 2022; originally announced February 2022.

Journal ref: A&A 660, A36 (2022)

GENGA II: GPU planetary N-body simulations with non-Newtonian forces and high number of particles

Authors: Simon L. Grimm, Joachim G. Stadel, Ramon Brasser, Matthias M. M. Meier, Christoph Mordasini

Abstract: We present recent updates and improvements of the graphical processing unit (GPU) N-body code GENGA. Modern state-of-the-art simulations of planet formation require the use of a very high number of particles to accurately resolve planetary growth and to quantify the effect of dynamical friction. At present the practical upper limit is in the range of 30,000 - 60,000 fully interactive particles; po… ▽ More We present recent updates and improvements of the graphical processing unit (GPU) N-body code GENGA. Modern state-of-the-art simulations of planet formation require the use of a very high number of particles to accurately resolve planetary growth and to quantify the effect of dynamical friction. At present the practical upper limit is in the range of 30,000 - 60,000 fully interactive particles; possibly a little more on the latest GPU devices. While the original hybrid symplectic integration method has difficulties to scale up to these numbers, we have improved the integration method by i) introducing higher level changeover functions and ii) code improvements to better use the most recent GPU hardware efficiently for such large simulations. We added treatments of non-Newtonian forces such as general relativity, tidal interaction, rotational deformation, the Yarkovsky effect, and Poynting-Robertson drag, as well as a new model to treat virtual collisions of small bodies in the solar system. We added new tools to GENGA, such as semi-active test particles that feel more massive bodies but not each other, a more accurate collision handling and a real-time openGL visualization. We present example simulations, including a 1.5 billion year terrestrial planet formation simulation that initially started with 65,536 particles, a 3.5 billion year simulation without gas giants starting with 32,768 particles, the evolution of asteroid fragments in the solar system, and the planetesimal accretion of a growing Jupiter simulation. GENGA runs on modern NVIDIA and AMD GPUs. △ Less

Submitted 8 June, 2022; v1 submitted 24 January, 2022; originally announced January 2022.

Comments: Accepted for publication in ApJ, 30 pages, 19 figures

Effects of pebble accretion on the growth and composition of planetesimals in the inner Solar System

Authors: J. Mah, R. Brasser, A. Bouvier, S. J. Mojzsis

Abstract: Recent work has shown that aside from the classical view of collisions by increasingly massive planetesimals, the accretion of mm- to m-sized 'pebbles' can also reproduce the mass-orbit distribution of the terrestrial planets. Here, we perform N-body simulations to study the effects of pebble accretion onto growing planetesimals of different diameters located in the inner Solar System. The simulat… ▽ More Recent work has shown that aside from the classical view of collisions by increasingly massive planetesimals, the accretion of mm- to m-sized 'pebbles' can also reproduce the mass-orbit distribution of the terrestrial planets. Here, we perform N-body simulations to study the effects of pebble accretion onto growing planetesimals of different diameters located in the inner Solar System. The simulations are run to occur during the lifetime of the gas disc while also simultaneously taking Jupiter's growth into account. We find that pebble accretion can increase the mass in the solid disc by at least a few times its initial mass with reasonable assumptions that pebbles fragment to smaller-sized grains at the snow line and that gas-disc-induced orbital migration effects are in force. Such a large contribution in mass by pebbles would seem to imply that the isotopic composition of the inner Solar System should be similar to the pebble source (i.e. outer Solar System). This implication appears to violate the observed nucleosynthetic isotopic dichotomy of the sampled Solar System. Thus, pebble accretion played little or no role in terrestrial planet formation. △ Less

Submitted 22 December, 2021; originally announced December 2021.

Comments: 18 pages, 11 figures, accepted for publication in MNRAS

The terrestrial planet formation paradox inferred from high-resolution N-body simulations

Authors: Jason Man Yin Woo, Ramon Brasser, Simon L. Grimm, Miles L. Timpe, Joachim Stadel

Abstract: Recent improvements to GPU hardware and the symplectic N-body code GENGA allow for unprecedented resolution in simulations of planet formation. In this paper, we report results from high-resolution N-body simulations of terrestrial planet formation that are mostly direct continuation of our previous 10 Myr simulations (Woo et al. 2021a) until 150 Myr. By assuming that Jupiter and Saturn have alway… ▽ More Recent improvements to GPU hardware and the symplectic N-body code GENGA allow for unprecedented resolution in simulations of planet formation. In this paper, we report results from high-resolution N-body simulations of terrestrial planet formation that are mostly direct continuation of our previous 10 Myr simulations (Woo et al. 2021a) until 150 Myr. By assuming that Jupiter and Saturn have always maintained their current eccentric orbits (EJS), we are able to achieve a reasonably good match to the current inner solar system architecture. However, due to the strong radial mixing that occurs in the EJS scenario, it has difficulties in explaining the known isotopic differences between bodies in the inner solar system, most notably between Earth and Mars. On the other hand, assuming initially circular orbits for Jupiter and Saturn (CJS) can reproduce the observed low degree of radial mixing in the inner solar system, while failing to reproduce the current architecture of the inner solar system. These outcomes suggest a possible paradox between dynamical structure and cosmochemical data for the terrestrial planets within the classical formation scenario. △ Less

Submitted 8 September, 2021; originally announced September 2021.

N-body simulations of planet formation via pebble accretion II. How various giant planets form

Authors: Soko Matsumura, Ramon Brasser, Shigeru Ida

Abstract: Aims. The connection between initial disc conditions and final orbital and physical properties of planets is not well-understood. In this paper, we numerically study the formation of planetary systems via pebble accretion and investigate the effects of disc properties such as masses, dissipation timescales, and metallicities on planet formation outcomes. Methods. We improved the N-body code SyMBA… ▽ More Aims. The connection between initial disc conditions and final orbital and physical properties of planets is not well-understood. In this paper, we numerically study the formation of planetary systems via pebble accretion and investigate the effects of disc properties such as masses, dissipation timescales, and metallicities on planet formation outcomes. Methods. We improved the N-body code SyMBA that was modified by taking account of new planet-disc interaction models and type II migration. We adopted the 'two-alpha disc' model to mimic the effects of both the standard disc turbulence and the mass accretion driven by the magnetic disc wind. Results. We successfully reproduced the overall distribution trends of semi-major axes, eccentricities, and planetary masses of extrasolar giant planets. We find that, when planet formation happens fast enough, giant planets are fully grown (Jupiter mass or higher) and are distributed widely across the disc. On the other hand, when planet formation is limited by the disc's dissipation, discs generally form low-mass cold Jupiters (CJs). Our simulations also naturally explain why hot Jupiters (HJs) tend to be alone and how the observed eccentricity-metallicity trends arise. The low-metallicity discs tend to form nearly circular and coplanar HJs in situ, because planet formation is slower than high-metallicity discs, and thus protoplanetary cores migrate significantly before gas accretion. The high-metallicity discs, on the other hand, generate HJs in situ or via tidal circularisation of eccentric orbits. Both pathways usually involve dynamical instabilities, and thus HJs tend to have broader eccentricity and inclination distributions. When giant planets with very wide orbits ('super-cold Jupiters') are formed, we find that they often belong to metal-rich stars, have eccentric orbits, and tend to have (~80%) companions interior to their orbits. △ Less

Submitted 15 April, 2021; originally announced April 2021.

Comments: 29 pages, 16 figures, 2 tables, A&A (in press)

Journal ref: A&A 650, A116 (2021)

Mars' formation can constrain the primordial orbits of the gas giants

Authors: Jason Man Yin Woo, Joachim Stadel, Simon Grimm, Ramon Brasser

Abstract: Recent high precision meteoritic data infers that Mars finished its accretion rapidly within 10 Myr of the beginning of the Solar system and had an accretion zone that did not entirely overlap with the Earth's. Here we present a detailed study of the accretion zone of planetary embryos from high resolution simulations of planetesimals in a disc. We found that all simulations with Jupiter and Satur… ▽ More Recent high precision meteoritic data infers that Mars finished its accretion rapidly within 10 Myr of the beginning of the Solar system and had an accretion zone that did not entirely overlap with the Earth's. Here we present a detailed study of the accretion zone of planetary embryos from high resolution simulations of planetesimals in a disc. We found that all simulations with Jupiter and Saturn on their current eccentric orbits (EJS) result in a similar accretion zone between fast-forming Mars and Earth region embryos. Assuming more circular orbits for Jupiter and Saturn (CJS), on the other hand, has a significantly higher chance of forming Mars with an accretion zone not entirely dominated by Earth and Venus region embryos, however CJS in general forms Mars slower than in EJS. By further quantifying the degree of overlap between accretion zones of embryos in different regions with the average overlap coefficient (OVL), we found that the OVL of CJS shows a better match with the OVL from a chondritic isotopic mixing model of Earth and Mars, which indicates that the giant planets are likely to have resided on more circular orbits than today during gas disc dissipation, matching their suggested pre-instability orbits. More samples, including those from Mercury and Venus, could potentially confirm this hypothesis. △ Less

Submitted 7 April, 2021; originally announced April 2021.

A new estimate for the age of highly-siderophile element retention in the lunar mantle from late accretion

Authors: R. Brasser, S. J. Mojzsis, S. C. Werner, O. Abramov

Abstract: Subsequent to the Moon's formation, late accretion to the terrestrial planets strongly modified the physical and chemical nature of silicate crusts and mantles. This alteration came in the form of melting through impacts, as well as the belated addition of volatiles and the highly siderophile elements (HSEs). Current debate centres on whether the lunar HSE record is representative of its whole lat… ▽ More Subsequent to the Moon's formation, late accretion to the terrestrial planets strongly modified the physical and chemical nature of silicate crusts and mantles. This alteration came in the form of melting through impacts, as well as the belated addition of volatiles and the highly siderophile elements (HSEs). Current debate centres on whether the lunar HSE record is representative of its whole late accretion history or alternatively that these were only retained in the mantle and crust after a particular time, and if so, when. Here we employ improved Monte Carlo impact simulations of late accretion onto the Moon and Mars and present an updated chronology based on new dynamical simulations of leftover planetesimals and the E-belt. We take into account the inefficient retention of colliding material. We compute the crater and basin densities on the Moon and Mars, the largest objects to strike these planets and the amount of material they accreted. Outputs are used to infer the mass in leftover planetesimals at a particular time period, which is then compared to the lunar HSE abundance. From this estimate we calculate a preferred lunar HSE retention age of ca. 4450 Ma which means that the modelled lunar mantle HSE abundances trace almost all of lunar late accretion. Based on our results, the surface ages of the lunar highlands are at least 4370 Ma. We find that the mass of leftover planetesimals with diameters Di<300 km at 4500 Ma that best fits the crater chronology is approximately 2x10^{-3} Earth mass (ME) while the mass of the E-belt was fixed at 4.5x10^{-4} ME. We also find that a leftover planetesimal mass in excess of 0.01 ME results in a lunar HSE retention age younger than major episodes of lunar differentiation and crust formation, which in turn violates geochemical constraints for the timing and intensity of late accretion to the Earth (Mojzsis et al., 2019). △ Less

Submitted 17 February, 2021; originally announced February 2021.

Comments: Accepted in Icarus after 5 rounds of review

Formation of Venus, Earth and Mars: Constrained by isotopes

Authors: H. Lammer, R. Brasser, A. Johansen, M. Scherf, M. Leitzinger

Abstract: We discuss the current state of knowledge of terrestrial planet formation from the aspects of different planet formation models and isotopic data from 182Hf-182W, U-Pb, lithophile-siderophile elements, 48Ca/44Ca isotope samples from planetary building blocks, 36Ar/38Ar, 20Ne/22Ne, 36Ar/22Ne isotope ratios in Venus' and Earth's atmospheres, the expected solar 3He abundance in Earth's deep mantle an… ▽ More We discuss the current state of knowledge of terrestrial planet formation from the aspects of different planet formation models and isotopic data from 182Hf-182W, U-Pb, lithophile-siderophile elements, 48Ca/44Ca isotope samples from planetary building blocks, 36Ar/38Ar, 20Ne/22Ne, 36Ar/22Ne isotope ratios in Venus' and Earth's atmospheres, the expected solar 3He abundance in Earth's deep mantle and Earth's D/H sea water ratios that shed light on the accretion time of the early protoplanets. Accretion scenarios that can explain the different isotope ratios, including a Moon-forming event after ca. 50 Myr, support the theory that the bulk of Earth's mass (>80%) most likely accreted within 10-30 Myr. From a combined analysis of the before mentioned isotopes, one finds that proto-Earth accreted 0.5-0.6 MEarth within the first ~4-5 Myr, the approximate lifetime of the protoplanetary disk. For Venus, the available atmospheric noble gas data are too uncertain for constraining the planet's accretion scenario accurately. However, from the available Ar and Ne isotope measurements, one finds that proto-Venus could have grown to 0.85-1.0 MVenus before the disk dissipated. Classical terrestrial planet formation models have struggled to grow large planetary embryos quickly from the tiniest materials within the typical lifetime of protoplanetary disks. Pebble accretion could solve this long-standing time scale controversy. Pebble accretion and streaming instabilities produce large planetesimals that grow into Mars-sized and larger planetary embryos during this early accretion phase. The later stage of accretion can be explained well with the Grand-Tack, annulus or depleted disk models. The relative roles of pebble accretion and planetesimal accretion/giant impacts are poorly understood and should be investigated with N-body simulations that include pebbles and multiple protoplanets. △ Less

Submitted 11 February, 2021; originally announced February 2021.

Comments: 49 pages, 7 figures. This is a preprint of an article published in Space Science Reviews. The final authenticated version can be found online at: https://doi.org/10.1007/s11214-020-00778-4

Journal ref: Space Sci Rev 217, 7 (2021)

Growing Mars fast: High-resolution GPU simulations of embryo formation

Authors: Jason Man Yin Woo, Simon L. Grimm, Ramon Brasser, Joachim Stadel

Abstract: Recent high precision meteoritic data improve constraints on the formation timescale and bulk composition of the terrestrial planets. High resolution N-body simulations allow direct comparison of embryo growth timescale and accretion zones to these constraints. In this paper, we present results of high resolution simulations for embryo formation from a disc of up to 41,000 fully-self gravitating p… ▽ More Recent high precision meteoritic data improve constraints on the formation timescale and bulk composition of the terrestrial planets. High resolution N-body simulations allow direct comparison of embryo growth timescale and accretion zones to these constraints. In this paper, we present results of high resolution simulations for embryo formation from a disc of up to 41,000 fully-self gravitating planetesimals with the GPU-based N-body code GENGA. Our results indicate that the growth of embryos are highly dependent on the initial conditions. More massive initial planetesimals, a shorter gas disc decay timescale and initially eccentric Jupiter and Saturn (EJS) all lead to faster growth of embryos. Asteroid belt material can thereby be implanted into the terrestrial planet region via sweeping secular resonances. This could possibly explain the rapid growth of Mars within 10 Myr inferred from its Hf-W chronology. The sweeping secular resonance almost completely clears the asteroid belt and deposits this material in the Mercury-Venus region, altering the composition of embryos there. This could result in embryos in the Mercury-Venus region accreting an unexpectedly high mass fraction from beyond 2 AU. Changing the initial orbits of Jupiter and Saturn to more circular (CJS) or assuming embryos formed in a gas free environment removes the sweeping secular resonance effect and thus greatly decreases material accreted from beyond 2 AU for Mercury-Venus region embryos. We therefore propose that rock samples from Mercury and Venus could aid greatly in deducing the condition and lifetime of the initial protoplanetary gas disc during planetesimal and embryo formation, as well as the initial architecture of the giant planets. △ Less

Submitted 22 January, 2021; originally announced January 2021.

Comments: Accepted for publication in Icarus

The Formation of Bilobate Comet Shapes through Sublimative Torques

Authors: Taylor K. Safrit, Jordan K. Steckloff, Amanda S. Bosh, David Nesvorny, Kevin Walsh, Ramon Brasser, David A. Minton

Abstract: Recent spacecraft and radar observations have found that ~70 percent of short-period comet nuclei, mostly Jupiter-family comets (JFCs), have bilobate shapes (two masses connected by a narrow neck). This is in stark contrast to the shapes of asteroids of similar sizes, of which ~14% are bilobate. This suggests that a process or mechanism unique to comets is producing these shapes. Here we show that… ▽ More Recent spacecraft and radar observations have found that ~70 percent of short-period comet nuclei, mostly Jupiter-family comets (JFCs), have bilobate shapes (two masses connected by a narrow neck). This is in stark contrast to the shapes of asteroids of similar sizes, of which ~14% are bilobate. This suggests that a process or mechanism unique to comets is producing these shapes. Here we show that the bilobate shapes of JFC nuclei are a natural byproduct of sublimative activity during their dynamical migration from their trans-Neptunian reservoir, through the Centaur population, and into the Jupiter family. We model the torques resulting from volatile sublimation during this dynamical migration and find that they tend to spin up these nuclei to disruption. Once disrupted, the rubble pile-like material properties of comet nuclei (tensile strengths of ~1-10 Pa and internal friction angles of ~35$^\circ$) cause them to reform as bilobate objects. We find that JFCs likely experienced rotational disruption events prior to entering the Jupiter family, which could explain the prevalence of bilobate shapes. These results suggest that the bilobate shapes of observed comets developed recently in their history (within the past ~1-10 Myr), rather than during solar system formation or collisions during planet migration and residency in the trans-Neptunian population. △ Less

Submitted 2 November, 2020; originally announced November 2020.

Isotopically distinct terrestrial planets via local accretion

Authors: Jingyi Mah, Ramon Brasser

Abstract: Combining isotopic constraints from meteorite data with dynamical models of planet formation proves to be advantageous in identifying the best model for terrestrial planet formation. Prior studies have shown that the probability of reproducing the distinct isotopic compositions of the Earth and Mars for both classical and Grand Tack models is very low. In the framework of the Grand Tack model, for… ▽ More Combining isotopic constraints from meteorite data with dynamical models of planet formation proves to be advantageous in identifying the best model for terrestrial planet formation. Prior studies have shown that the probability of reproducing the distinct isotopic compositions of the Earth and Mars for both classical and Grand Tack models is very low. In the framework of the Grand Tack model, for Mars to be isotopically different from the Earth, it had to form under very specific conditions. Here, we subjected a fairly new and unexplored model--the depleted disc model--to the test. It presupposes that the region in the inner protoplanetary disc from Mars' orbit and beyond is depleted in mass such that Mars is left with insufficient material to grow to a larger size. Our aim is to test the whether the distinct isotopic compositions of the Earth and Mars are a natural outcome of this model. We found that the terrestrial planets accrete material mostly locally and have feeding zones that are sufficiently distinct. The Earth and Mars, and by extension, Venus, can have distinct isotopic compositions if there is an isotopic gradient in the terrestrial planet region of the protoplanetary disc. Our results suggest that the material in the inner Solar System most likely did not undergo substantial mixing that homogenised the potential isotopic gradient, in contrast to the Grand Tack model where the feeding zones of the terrestrial planets are nearly identical due to the mixing of material by Jupiter's migration. △ Less

Submitted 20 August, 2020; originally announced August 2020.

Comments: 18 pages, 11 figures, 2 tables, accepted for publication in Icarus

Modification of the composition and density of Mercury from late accretion

Authors: Ryuki Hyodo, Hidenori Genda, Ramon Brasser

Abstract: Late accretion is a process that strongly modulated surface geomorphic and geochemical features of Mercury. Yet, the fate of the impactors and their effects on Mercury's surface through the bombardment epoch are not clear. Using Monte-Carlo and analytical approaches of cratering impacts, we investigate the physical and thermodynamical outcomes of late accretion on Mercury. Considering the uncertai… ▽ More Late accretion is a process that strongly modulated surface geomorphic and geochemical features of Mercury. Yet, the fate of the impactors and their effects on Mercury's surface through the bombardment epoch are not clear. Using Monte-Carlo and analytical approaches of cratering impacts, we investigate the physical and thermodynamical outcomes of late accretion on Mercury. Considering the uncertainties in late accretion, we develop scaling laws for the following parameters as a function of impact velocity and total mass of late accretion: (1) depth of crustal erosion, (2) the degree of resurfacing, and (3) mass accreted from impactor material. Existing dynamical models indicate that Mercury experienced an intense impact bombardment (a total mass of $\sim 8 \times 10^{18} - 8 \times 10^{20}$ kg with a typical impact velocity of $30-40$ km s$^{-1}$) after $4.5$ Ga. With this, we find that late accretion could remove 50 m to 10 km of the early crust of Mercury, but the change to its core-to-mantle ratio is negligible. Alternatively, the mantles of putative differentiated planetesimals in the early solar system could be more easily removed and their respective core fraction increased, if Mercury ultimately accreted from such objects. Although the cratering is notable for erasing the older geological surface records on Mercury, we show that $\sim 40-50$ wt% of the impactor's exogenic materials, including the volatile-bearing materials, can be heterogeneously implanted on Mercury's surface as a late veneer ($>3\times 10^{18}-1.6 \times 10^{19}$ kg in total). About half of the accreted impactor's materials are vaporized, and the rest is completely melted upon the impact. We expect that the further interplay between our theoretical results and forthcoming surface observations of Mercury, including the BepiColombo mission, will lead us to a better understanding of Mercury's origin and evolution. △ Less

Submitted 19 August, 2020; originally announced August 2020.

Comments: 19 pages, 7 figures. Accepted for publication in Icarus

No evidence for interstellar planetesimals trapped in the Solar System

Authors: A. Morbidelli, K. Batygin, R. Brasser, S. Raymond

Abstract: In two recent papers published in MNRAS, Namouni and Morais (2018, 2020) claimed evidence for the interstellar origin of some small Solar System bodies, including i) objects in retrograde co-orbital motion with the giant planets, and ii) the highly-inclined Centaurs. Here, we discuss the flaws of those papers that invalidate the authors' conclusions. Numerical simulations backwards in time are not… ▽ More In two recent papers published in MNRAS, Namouni and Morais (2018, 2020) claimed evidence for the interstellar origin of some small Solar System bodies, including i) objects in retrograde co-orbital motion with the giant planets, and ii) the highly-inclined Centaurs. Here, we discuss the flaws of those papers that invalidate the authors' conclusions. Numerical simulations backwards in time are not representative of the past evolution of real bodies. Instead, these simulations are only useful as a means to quantify the short dynamical lifetime of the considered bodies and the fast decay of their population. In light of this fast decay, if the observed bodies were the survivors of populations of objects captured from interstellar space in the early Solar System, these populations should have been implausibly large (e.g. about 10 times the current main asteroid belt population for the retrograde coorbital of Jupiter). More likely, the observed objects are just transient members of a population that is maintained in quasi-steady state by a continuous flux of objects from some parent reservoir in the distant Solar System. We identify in the Halley type comets and the Oort cloud the most likely sources of retrograde coorbitals and highly-inclined Centaurs. △ Less

Submitted 8 June, 2020; originally announced June 2020.

Comments: in press in MNRAS

A new and simple prescription for planet orbital migration and eccentricity damping by planet-disc interactions based on dynamical friction

Authors: Shigeru Ida, Takayuki Muto, Soko Matsumura, Ramon Brasser

Abstract: During planet formation gravitational interaction between a planetary embryo and the protoplanetary gas disc causes orbital migration of the planetary embryo, which plays an important role in shaping the final planetary system. While migration sometimes occurs in the supersonic regime, wherein the relative velocity between the planetary embryo and the gas is higher than the sound speed, migration… ▽ More During planet formation gravitational interaction between a planetary embryo and the protoplanetary gas disc causes orbital migration of the planetary embryo, which plays an important role in shaping the final planetary system. While migration sometimes occurs in the supersonic regime, wherein the relative velocity between the planetary embryo and the gas is higher than the sound speed, migration prescriptions proposed thus far describing the planet-disc interaction force and the timescales of orbital change in the supersonic regime are inconsistent with one another. Here we discuss the details of existing prescriptions in the literature and derive a new simple and intuitive formulation for planet-disc interactions based on dynamical friction that can be applied in both supersonic and subsonic cases. While the existing prescriptions assume particular disc models, ours include the explicit dependence on the disc parameters; hence it can be applied to discs with any radial surface density and temperature dependence (except for the local variations with radial scales less than the disc scale height). Our prescription will reduce the uncertainty originating from different literature formulations of planet migration and will be an important tool to study planet accretion processes, especially when studying the formation of close-in low-mass planets that are commonly found in exoplanetary systems. △ Less

Submitted 6 May, 2020; v1 submitted 16 April, 2020; originally announced April 2020.

Comments: 10 pages, 1 figure, accepted for publication in MNRAS; typos corrected, the reference list was completed

Impact bombardment chronology of the terrestrial planets from 4.5 Ga to 3.5 Ga

Authors: R. Brasser, S. C. Werner, S. J. Mojzsis

Abstract: Subsequent to the Moon's formation, late accretion to the terrestrial planets modified their silicate crusts and mantles. We combine dynamical N-body and Monte Carlo simulations to determine impact probabilities, impact velocities, and expected mass augmentation onto the terrestrial planets from three sources: planetesimals left over from primary accretion, asteroids from the E-belt, and comets ar… ▽ More Subsequent to the Moon's formation, late accretion to the terrestrial planets modified their silicate crusts and mantles. We combine dynamical N-body and Monte Carlo simulations to determine impact probabilities, impact velocities, and expected mass augmentation onto the terrestrial planets from three sources: planetesimals left over from primary accretion, asteroids from the E-belt, and comets arriving from the outer Solar System. We present estimates of the amount of cometary material striking the terrestrial planets in an early episode of planetesimal-driven giant planet migration (Mojzsis et al., 2019). The Moon and Mars suffer proportionally higher cometary accretion than Venus and Earth. The background mass addition from small leftover planetesimals to Earth and Mars is far less than independent estimates based on their mantle abundances of highly-siderophile elements. This supports the theory that both planets were struck by single large bodies that delivered most of their mass augmentation since primary accretion. We present fits to the impact chronologies valid from 4500 Ma to ca. 3700 Ma. The lunar timeline thus obtained is at odds with both the calibrated Neukum (Neukum et al., 2001) and Werner (Werner et al., 2014; Werner, 2019) chronologies. For Mars, the match with its calibrated Werner chronology is no better; by increasing the mass of the E-belt by a factor of four the dynamical lunar and martian chronologies are in line with that of Werner (2019) and match constraints from the current population of Hungaria asteroids. Yet, neither of our dynamical timelines fit well with that of Neukum. The dynamical lunar and martian chronologies are also different from each other. Consequently, the usual extrapolation of such chronologies from one planetary body to the other is technically inappropriate. △ Less

Submitted 24 October, 2019; originally announced October 2019.

Comments: Accepted in Icarus

When did Life Likely Emerge on Earth in an RNA-First Process?

Authors: Steven A. Benner, Elizabeth A. Bell, Elisa Biondi, Ramon Brasser, Thomas Carell, Hyo-Joong Kim, Stephen J. Mojzsis, Arthur Omran, Matthew A. Pasek, Dustin Trail

Abstract: The widespread presence of ribonucleic acid (RNA) catalysts and cofactors in Earth's biosphere today suggests that RNA was the first biopolymer to support Darwinian evolution. However, most "path-hypotheses" to generate building blocks for RNA require reduced nitrogen-containing compounds not made in useful amounts in the CO2-N2-H2O atmospheres of the Hadean. We review models for Earth's impact hi… ▽ More The widespread presence of ribonucleic acid (RNA) catalysts and cofactors in Earth's biosphere today suggests that RNA was the first biopolymer to support Darwinian evolution. However, most "path-hypotheses" to generate building blocks for RNA require reduced nitrogen-containing compounds not made in useful amounts in the CO2-N2-H2O atmospheres of the Hadean. We review models for Earth's impact history that invoke a single ~10^23 kg impactor (Moneta) to account for measured amounts of platinum, gold, and other siderophilic ("iron-loving") elements on the Earth and Moon. If it were the last sterilizing impactor, Moneta would have reduced the atmosphere but not its mantle, opening a "window of opportunity" for RNA synthesis, a period when RNA precursors rained from the atmosphere to land holding oxidized minerals that stabilize advanced RNA precursors and RNA. Surprisingly, this combination of physics, geology, and chemistry suggests a time when RNA formation was most probable, ~120 +/- 100 million years after Moneta's impact, or ~4.36 +/- 0.1 billion years ago. Uncertainties in this time are driven by uncertainties in rates of productive atmosphere loss and amounts of sub-aerial land. △ Less

Submitted 29 August, 2019; originally announced August 2019.

Comments: Paper accepted to ChemSystemsChem [https://onlinelibrary.wiley.com/journal/25704206] August 22, 2019

Journal ref: 19. Benner, S.A., Bell, E.A., Biondi, E., Brasser, R., Carell, T., Kim, H-J., Mojzsis, S.J., Omran, A., Pasek, M.A., and Trail, D. (2020) When did Life Likely Emerge on Earth in an RNA-First Process? ChemSystemsChem 2, e1900035

Mars in the aftermath of a colossal impact

Authors: Jason Man Yin Woo, Hidenori Genda, Ramon Brasser, Stephen J. Mojzsis

Abstract: The abundance of highly siderophile elements (HSEs) inferred for Mars' mantle from martian meteorites implies a Late Veneer (LV) mass addition of ~0.8 wt% with broadly chondritic composition. Late accretion to Mars by a differentiated Ceres-sized (~1000 km diameter) object can account for part of the requisite LV mass, and geochronological constraints suggests that this must have occurred no later… ▽ More The abundance of highly siderophile elements (HSEs) inferred for Mars' mantle from martian meteorites implies a Late Veneer (LV) mass addition of ~0.8 wt% with broadly chondritic composition. Late accretion to Mars by a differentiated Ceres-sized (~1000 km diameter) object can account for part of the requisite LV mass, and geochronological constraints suggests that this must have occurred no later than ca. 4480 Ma. Here, we analyze the outcome of the hypothetical LV giant impact to Mars with smoothed particle hydrodynamics simulations together with analytical theory. Results show that, in general about 50% of the impactor's metallic core shatters into ~10m fragments that subsequently fragment into sub-mm metallic hail at re-accretion. This returns a promising delivery of HSEs into martian mantle compared to either a head-on and hit-and-run collision; in both cases,<10% of impactor's core materials are fragmented and finally embedded in the martian mantle. Isotopic evidence from martian meteorites, and interpretations from atmospheric mapping data show that a global surface water reservoir could be present during the early Noachian (before ca. 4100 Ma). The millimeter-sized metal hail could thus react with a martian hydrosphere to generate ~3 bars of H2, which is adequate to act as a greenhouse and keep early Mars warm. Yet, we also find that this atmosphere is transient. It typically survives shorter than 3 Myr based on the expected extreme ultraviolet (EUV) flux of the early Sun; if the Sun was a slow rotator an accordingly weaker EUV flux could extend this lifetime to >10 Myr. A dense pre-Noachian CO2 atmosphere should lower the escape efficiency of hydrogen by IR emission. A more detailed hydrodynamic atmospheric model of this early hydrogen atmosphere is warranted to examine its effect on pre-Noachian Mars. △ Less

Submitted 20 June, 2019; originally announced June 2019.

Journal ref: Icarus 333 (2019) 87-95

The tidal parameters of TRAPPIST-1 b and c

Authors: R. Brasser, A. C. Barr, V. Dobos

Abstract: The TRAPPIST-1 planetary system consists of seven planets within 0.05 au of each other, five of which are in a multi-resonant chain. {These resonances suggest the system formed via planet migration; subsequent tidal evolution has damped away most of the initial eccentricities. We used dynamical N-body simulations to estimate how long it takes for the multi-resonant configuration that arises during… ▽ More The TRAPPIST-1 planetary system consists of seven planets within 0.05 au of each other, five of which are in a multi-resonant chain. {These resonances suggest the system formed via planet migration; subsequent tidal evolution has damped away most of the initial eccentricities. We used dynamical N-body simulations to estimate how long it takes for the multi-resonant configuration that arises during planet formation to break. From there we use secular theory to pose limits on the tidal parameters of planets b and c. We calibrate our results against multi-layered interior models constructed to fit the masses and radii of the planets, from which the tidal parameters are computed independently.} The dynamical simulations show that the planets typically go unstable 30 Myr after their formation. {Assuming synchronous rotation throughout} we compute $\frac{k_2}{Q} \gtrsim 2\times 10^{-4}$ for planet b and $\frac{k_2}{Q} \gtrsim 10^{-3}$ for planet c. Interior models yield $(0.075-0.37) \times 10^{-4}$ for TRAPPIST-1 b and $(0.4-2)\times 10^{-4}$ for TRAPPIST-1 c. The agreement between the {dynamical and interior} models is not too strong, but is still useful to constrain the dynamical history of the system. We suggest that this two-pronged approach could be of further use in other multi-resonant systems if the planet's orbital and interior parameters are sufficiently well known. △ Less

Submitted 1 May, 2019; originally announced May 2019.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society

The origin of the cratering asymmetry on Triton

Authors: Jingyi Mah, Ramon Brasser

Abstract: The source of impactors capable of producing the cratering asymmetry on Triton is open to debate. While hyperbolic comets are the likely source of impactors, the theoretically predicted crater distribution due to impacts from comets is incompatible with the observed crater distribution on Triton. In this study, we adopted a dynamical approach to investigate the crater distribution produced by thre… ▽ More The source of impactors capable of producing the cratering asymmetry on Triton is open to debate. While hyperbolic comets are the likely source of impactors, the theoretically predicted crater distribution due to impacts from comets is incompatible with the observed crater distribution on Triton. In this study, we adopted a dynamical approach to investigate the crater distribution produced by three sources of impactors: Neptune-centric impactors on prograde, low-inclination orbits, Neptune-centric impactors on high-inclination orbits such as its irregular satellites, and heliocentric impactors with hyperbolic orbits. We employ N-body simulations which take into account the dynamics of the impacts. We found that none of the impactor sources alone are able to reconcile with the crater distribution on Triton. The crater distribution produced by a low-inclination Neptune-centric source is narrowly concentrated on Triton's leading hemisphere whereas the heliocentric impactors produce too many craters on the trailing hemisphere compared to the observed crater distribution. All three impactor populations produce craters on the trailing hemisphere of Triton. We derived the impact distribution of a mixed population of impactors using a Monte Carlo approach and found that if the population consists of mainly heliocentric impactors (70%), the resulting impact distribution fits the observed crater distribution on Triton's leading hemisphere well. △ Less

Submitted 16 April, 2019; originally announced April 2019.

Comments: 7 pages, 8 figures, 1 table, accepted for publication in MNRAS

Onset of giant planet migration before 4480 million years ago

Authors: Stephen J. Mojzsis, Ramon Brasser, Nigel M. Kelly, Oleg Abramov, Stephanie C. Werner

Abstract: Immediately after their formation, the terrestrial planets experienced intense impact bombardment by comets, leftover planetesimals from primary accretion, and asteroids. This temporal interval in solar system evolution, termed late accretion, thermally and chemically modified solid planetary surfaces and may have impeded the emergence of life on the Hadean Earth. The sources and tempo of late acc… ▽ More Immediately after their formation, the terrestrial planets experienced intense impact bombardment by comets, leftover planetesimals from primary accretion, and asteroids. This temporal interval in solar system evolution, termed late accretion, thermally and chemically modified solid planetary surfaces and may have impeded the emergence of life on the Hadean Earth. The sources and tempo of late accretion are, however, vague. Here, we present a timeline that relates variably retentive radiometric ages from asteroidal meteorites, to new dynamical models of late accretion that invokes giant planet migration. Reconciliation of the geochronological data with dynamical models shows that giant planet migration immediately leads to an intense 30 Myr influx of comets to the entire solar system. The absence of whole-sale crustal reset ages after 4450 Ma for the most resilient chronometers from Earth, Moon, Mars, Vesta and various meteorite parent bodies confines the onset of giant planet migration to no later than ca. 4480 Ma. Waning impacts from planetesimals, asteroids (and a minor cometary component) continue to strike the inner planets through a protracted monotonic decline in impactor flux; this is in agreement with predictions from crater chronology. Amended global 3-D thermal analytical bombardment models derived from our new impact mass-production functions show that persistent niches for prebiotic chemistry on the early Hadean Earth could endure late accretion for at least the last 4400 Myr. △ Less

Submitted 21 March, 2019; originally announced March 2019.

Comments: Main text: 46564 characters with spaces/7549 words Tables: 3 Figures:7 References: 115

Journal ref: The Astrophysical Journal, 881:44 (13pp), 2019 August 10

Impact bombardment on the regular satellites of Jupiter and Uranus during an episode of giant planet migration

Authors: E. W. Wong, R. Brasser, S. C. Werner

Abstract: The intensity and effects of early impact bombardment on the major satellites of the giant planets during an episode of giant planet migration is still poorly known. We use a combination of dynamical N-body and Monte Carlo simulations to determine impact probabilities, impact velocities, and expected masses that collide with these satellites to determine the chronology of impacts during the migrat… ▽ More The intensity and effects of early impact bombardment on the major satellites of the giant planets during an episode of giant planet migration is still poorly known. We use a combination of dynamical N-body and Monte Carlo simulations to determine impact probabilities, impact velocities, and expected masses that collide with these satellites to determine the chronology of impacts during the migration. Volatile loss through bombardment is typically 20% for Miranda, a few percents for the larger Uranian satellites and negligible for the Galilean satellites. Due to its small size and the high impact velocity there is a >99% chance that Miranda suffered a catastrophic impact that shattered the satellite. Subsequent re-accretion from a circum-Uranian ring could account for its peculiar surface morphology and low density. The probability to destroy Ariel and Umbriel is 15% and 1% for Titania and Oberon. Approximately 90% of the mass in planetesimals that passes through the Jovian and Uranian satellite systems (about $4 {\rm \ M_{\oplus}}$ and $2 {\rm \ M_{\oplus}}$ respectively) does so in about 15 Myr. This extremely rapid and intense bombardment causes repeated local crustal melting on all satellites. The combination of these effects results in an entirely different impact chronology than that of the inner solar system. We conclude that the simple extrapolation of the lunar chronology to the outer solar system satellites is not correct. The tail end (after 25 Myr) of the chronology function has an e-folding time of 100 Myr at Jupiter, but follows a cumulative Weibull distribution at Uranus, making direct comparisons between the gas and ice giant planets difficult. Based on our results the surfaces of the Uranian satellites, Callisto, and possibly Ganymede, are all about the same age, and are roughly 150 Myr younger than the timing of the dynamical instability. △ Less

Submitted 12 November, 2018; originally announced November 2018.

Comments: 24 pages, 5 figures, accepted for publication in Earth and Planetary Science Letters

How planetary growth outperforms migration

Authors: Anders Johansen, Shigeru Ida, Ramon Brasser

Abstract: Planetary migration is a major challenge for planet formation theories. The speed of Type I migration is proportional to the mass of a protoplanet, while the final decade of growth of a pebble-accreting planetary core takes place at a rate that scales with the mass to the two-thirds power. This results in planetary growth tracks (i.e., the evolution of a protoplanet's mass versus its distance from… ▽ More Planetary migration is a major challenge for planet formation theories. The speed of Type I migration is proportional to the mass of a protoplanet, while the final decade of growth of a pebble-accreting planetary core takes place at a rate that scales with the mass to the two-thirds power. This results in planetary growth tracks (i.e., the evolution of a protoplanet's mass versus its distance from the star) that become increasingly horizontal (migration-dominated) with rising mass of the protoplanet. It has been shown recently that the migration torque on a protoplanet is reduced proportional to the relative height of the gas gap carved by the growing planet. Here we show from 1-D simulations of planet-disc interaction that the mass at which a planet carves a 50% gap is approximately 2.3 times the pebble isolation mass. Our measurements of the pebble isolation mass from 1-D simulations match published 3-D results relatively well, except at very low viscosities where the 3-D pebble isolation mass is significantly higher, possibly due to gap edge instabilities not captured in 1-D. The pebble isolation mass demarks the transition from pebble accretion to gas accretion. Gas accretion to form gas-giant planets therefore takes place over a few astronomical units of migration after reaching first the pebble isolation mass and, shortly after, the 50% gap mass. Our results demonstrate how planetary growth can outperform migration, both during core accretion and during gas accretion, even when the Stokes number of the pebbles is small, St~0.01, and the pebble-to-gas flux ratio in the protoplanetary disc is in the nominal range of 0.01-0.02. We find that planetary growth is very rapid in the first million years of the protoplanetary disc and that the probability for forming gas-giant planets increases with the initial size of the protoplanetary disc and with decreasing turbulent diffusion. △ Less

Submitted 10 January, 2019; v1 submitted 1 November, 2018; originally announced November 2018.

Comments: Accepted for publication in Astronomy & Astrophysics

Enhanced constraints on the interior composition and structure of terrestrial exoplanets

Authors: Haiyang S. Wang, Fan Liu, Trevor R. Ireland, Ramon Brasser, David Yong, Charles H. Lineweaver

Abstract: Exoplanet interior modelling usually makes the assumption that the elemental abundances of a planet are identical to those of its host star. Host stellar abundances are good proxies of planetary abundances, but only for refractory elements. This is particularly true for terrestrial planets, as evidenced by the relative differences in bulk chemical composition between the Sun and the Earth and othe… ▽ More Exoplanet interior modelling usually makes the assumption that the elemental abundances of a planet are identical to those of its host star. Host stellar abundances are good proxies of planetary abundances, but only for refractory elements. This is particularly true for terrestrial planets, as evidenced by the relative differences in bulk chemical composition between the Sun and the Earth and other inner solar system bodies. The elemental abundances of a planet host star must therefore be devolatilised in order to correctly represent the bulk chemical composition of its terrestrial planets. Furthermore, nickel and light elements make an important contribution alongside iron to the core of terrestrial planets. We therefore adopt an extended chemical network of the core, constrained by an Fe/Ni ratio of 18 $\pm$ 4 (by number). By applying these constraints to the Sun, our modelling reproduces the composition of the mantle and core, as well as the core mass fraction of the Earth. We also apply our modelling to four exoplanet host stars with precisely measured elemental abundances: Kepler-10, Kepler-20, Kepler-21 and Kepler-100. If these stars would also host terrestrial planets in their habitable zone, we find that such planets orbiting Kepler-21 would be the most Earth-like, while those orbiting Kepler-10 would be the least. To assess the similarity of a rocky exoplanet to the Earth in terms of interior composition and structure, high-precision host stellar abundances are critical. Our modelling implies that abundance uncertainties should be better than $\sim$ 0.04 dex for such an assessment to be made. △ Less

Submitted 3 December, 2018; v1 submitted 10 October, 2018; originally announced October 2018.

Comments: Updated based on the version of record: updated references, acknowledgements, footnotes and appendix (with the addition of a URL for the access to the source codes); 14 pages, 6 figures, 4 tables; published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 482, Issue 2, 11 January 2019, Pages 2222-2233

Plausible home stars of the interstellar object 'Oumuamua found in Gaia DR2

Authors: C. A. L. Bailer-Jones, D. Farnocchia, K. J. Meech, R. Brasser, M. Micheli, S. Chakrabarti, M. W. Buie, O. R. Hainaut

Abstract: The first detected interstellar object 'Oumuamua that passed within 0.25au of the Sun on 2017 September 9 was presumably ejected from a stellar system. We use its newly determined non-Keplerian trajectory together with the reconstructed Galactic orbits of 7 million stars from Gaia DR2 to identify past close encounters. Such an "encounter" could reveal the home system from which 'Oumuamua was eject… ▽ More The first detected interstellar object 'Oumuamua that passed within 0.25au of the Sun on 2017 September 9 was presumably ejected from a stellar system. We use its newly determined non-Keplerian trajectory together with the reconstructed Galactic orbits of 7 million stars from Gaia DR2 to identify past close encounters. Such an "encounter" could reveal the home system from which 'Oumuamua was ejected. The closest encounter, at 0.60pc (0.53-0.67pc, 90% confidence interval), was with the M2.5 dwarf HIP 3757 at a relative velocity of 24.7km/s, 1Myr ago. A more distant encounter (1.6pc) but with a lower encounter (ejection) velocity of 10.7km/s was with the G5 dwarf HD 292249, 3.8Myr ago. Two more stars have encounter distances and velocities intermediate to these. The encounter parameters are similar across six different non-gravitational trajectories for 'Oumuamua. Ejection of 'Oumuamua by scattering from a giant planet in one of the systems is plausible, but requires a rather unlikely configuration to achieve the high velocities found. A binary star system is more likely to produce the observed velocities. None of the four home candidates have published exoplanets or are known to be binaries. Given that the 7 million stars in Gaia DR2 with 6D phase space information is just a small fraction of all stars for which we can eventually reconstruct orbits, it is a priori unlikely that our current search would find 'Oumuamua's home star system. As 'Oumuamua is expected to pass within 1pc of about 20 stars and brown dwarfs every Myr, the plausibility of a home system depends also on an appropriate (low) encounter velocity. △ Less

Submitted 24 September, 2018; originally announced September 2018.

Comments: Accepted to The Astronomical Journal

Orbital evolution of Saturn's mid-sized moons and the tidal heating of Enceladus

Authors: Ayano Nakajima, Shigeru Ida, Jun Kimura, Ramon Brasser

Abstract: The formation and orbital evolution of Saturn's inner mid-sized moons are still debated. The most puzzling aspects are 1) how the Tethys-Dione pair and the Mimas-Enceladus pair passed through their strong 3:2 mean-motion resonances during the tidal orbital evolution, and 2) the current strong heat flow from Enceladus, which is a few orders of magnitude higher than the tidal energy dissipation caus… ▽ More The formation and orbital evolution of Saturn's inner mid-sized moons are still debated. The most puzzling aspects are 1) how the Tethys-Dione pair and the Mimas-Enceladus pair passed through their strong 3:2 mean-motion resonances during the tidal orbital evolution, and 2) the current strong heat flow from Enceladus, which is a few orders of magnitude higher than the tidal energy dissipation caused by the present orbital eccentricity of Enceladus. Here we perform N-body simulations of the moons' orbital evolution from various initial conditions -- assuming that the moons were formed from Saturn's hypothetical massive ring -- and investigate possible paths to solve the above difficulties. If the influence of the rings is neglected, we find that the Tethys-Dione pair cannot avoid becoming trapped in the mean-motion resonances as they recede from Saturn, and that the Tethys-Enceladus pair cannot avoid collisions after the resonance trapping, in case Saturn's quality factor is smaller than 15,000. These findings are inconsistent with the current orbital configuration. However, taking into account both the eccentricity excitation and the orbital expansion caused by the ring torque, we find that these resonance captures are avoided. With the high eccentricity pumped up by the torque, Enceladus passes through the resonances with Tethys, and the Dione-Tethys pair passes through their 2:1 and possibly the 3:2 resonance as well. After Enceladus resides beyond the 2:1 resonance with the outer ring edge, the eccentricity can be tidally damped. While this is a promising path of evolution, in most runs, Enceladus collides with Tethys by the excited eccentricity. There is a hint that a ring mass decrease could avoid the collision. The tidal heat due to the eccentricity excitation by the ring torque is stored in the moons and slowly radiated away. It may account for the current high heat flow. △ Less

Submitted 3 September, 2018; originally announced September 2018.

Comments: 37 pages, 7 figures, accepted for publication in Icarus

Trapping low-mass planets at the inner edge of the protostellar disc

Authors: R. Brasser, S. Matsumura, T. Muto, S. Ida

Abstract: The formation of multiple close-in low-mass exoplanets is still a mystery. The challenge is to build a system wherein the outermost planet is beyond 0.2 AU from the star. Here we investigate how the prescription for type I planet migration affects the ability to trap multiple planets in a resonant chain near the inner edge of the protostellar disc. A sharp edge modelled as a hyperbolic tangent fun… ▽ More The formation of multiple close-in low-mass exoplanets is still a mystery. The challenge is to build a system wherein the outermost planet is beyond 0.2 AU from the star. Here we investigate how the prescription for type I planet migration affects the ability to trap multiple planets in a resonant chain near the inner edge of the protostellar disc. A sharp edge modelled as a hyperbolic tangent function coupled with supersonic corrections to the classical type I migration torques results in the innermost planets being pushed inside the cavity through resonant interaction with farther planets because migration is starward at slightly supersonic eccentricities. Planets below a few Earth masses are generally trapped in a resonant chain with the outermost planet near the disc edge, but long-term stability is not guaranteed. For more massive planets the migration is so fast that the eccentricity of the innermost resonant pair is excited to highly supersonic levels due to decreased damping on the innermost planet as it is pushed inside the cavity; collisions frequently occur and the system consists one or two intermediate-mass planets residing closer to the star than the disc's inner edge. We found a neat pileup of resonant planets outside the disc edge only if the corotation torque does not rapidly diminish at high eccentricity. We call for detailed studies on planet migration near the disc's inner edge, which is still uncertain, and for an improved understanding of eccentricity damping and disc torques in the supersonic regime. △ Less

Submitted 14 August, 2018; originally announced August 2018.

Comments: Accepted for publication in ApJ Letters

The curious case of Mars formation

Authors: Jason Man Yin Woo, Ramon Brasser, Soko Matsumura, Stephen J. Mojzsis, Shigeru Ida

Abstract: Dynamical models of planet formation coupled with cosmochemical data from martian meteorites show that Mars' isotopic composition is distinct from that of Earth. Reconciliation of formation models with meteorite data require that Mars grew further from the Sun than its present position. Here, we evaluate this compositional difference in more detail by comparing output from two $N$-body planet form… ▽ More Dynamical models of planet formation coupled with cosmochemical data from martian meteorites show that Mars' isotopic composition is distinct from that of Earth. Reconciliation of formation models with meteorite data require that Mars grew further from the Sun than its present position. Here, we evaluate this compositional difference in more detail by comparing output from two $N$-body planet formation models. The first of these planet formation models simulates what is termed the "Classical" case wherein Jupiter and Saturn are kept in their current orbits. We compare these results with another model based on the "Grand Tack", in which Jupiter and Saturn migrate through the primordial asteroid belt. Our estimate of the average fraction of chondrite assembled into Earth and Mars assumes that the initial solid disk consists of only sources of enstatite chondrite composition in the inner region, and ordinary chondrite in the outer region. Results of these analyses show that both models tend to yield Earth and Mars analogues whose accretion zones overlap. The Classical case fares better in forming Mars with its documented composition (29% to 68% enstatite chondrite plus 32% to 67% ordinary chondrite) though the Mars analogues are generally too massive. However, if we include the restriction of mass on the Mars analogues, the Classical model does not work better. We also further calculate the isotopic composition of $^{17} \rm O$, $^{50} \rm Ti$, $^{54} \rm Cr$, $^{142} \rm Nd$, $^{64} \rm Ni$, and $^{92} \rm Mo$ in the martian mantle from the Grand Tack simulations. We find that it is possible to match the calculated isotopic composition of all the above elements in Mars' mantle with their measured values, but the resulting uncertainties are too large to place good restriction on the early dynamical evolution and birth place of Mars. △ Less

Submitted 12 June, 2018; v1 submitted 31 May, 2018; originally announced June 2018.

Comments: 14 pages, 8 figures, presented in the 2017 DPS meeting, 2018 Solar system symposium in Sapporo and 2018 AOGS annual meeting, Accepted for publishing in A&A

Journal ref: A&A 617, A17 (2018)

Jupiter's influence on the building blocks of Mars and Earth

Authors: R. Brasser, N. Dauphas, S. J. Mojzsis

Abstract: Radiometric dating indicates that Mars accreted in the first ~4 Myr of solar system formation, which coincides with the formation and possible migration of Jupiter. While nebular gas from the protoplanetary disk was still present, Jupiter may have migrated inwards and tacked at 1.5 AU in a 3:2 resonance with Saturn. This migration excited planetary building blocks in the inner solar system, result… ▽ More Radiometric dating indicates that Mars accreted in the first ~4 Myr of solar system formation, which coincides with the formation and possible migration of Jupiter. While nebular gas from the protoplanetary disk was still present, Jupiter may have migrated inwards and tacked at 1.5 AU in a 3:2 resonance with Saturn. This migration excited planetary building blocks in the inner solar system, resulting in extensive mixing and planetesimal removal. Here we evaluate the plausible nature of Mars' building blocks, focusing in particular on how its growth was influenced by the formation and migration of Jupiter. We use a combination of dynamical simulations and an isotopic mixing model that traces the accretion of elements with different affinities for metal. Dynamical simulations show that Jupiter's migration causes the late stages of Earth's and Mars' accretion to be dominated by EC-type (enstatite chondrite) material due to the loss of OC (ordinary chondrite) planetesimals. Our analysis of available isotopic data for SNC meteorites shows that Mars consists of approximately 68%$^{+0}_{-39}$ EC + 32%$^{+35}_{-0}$ OC by mass (2σ). The large uncertainties indicate that isotopic analyses of martian samples are for the most part too imprecise to definitely test model predictions; in particular it remains uncertain whether or not Mars accreted predominantly EC material in the latter stages of its formation history. Dynamical simulations also provide no definitive constraint on Mars' accretion history due to the great variety of dynamical pathways that the martian embryo exhibits. The present work calls for new measurements of isotopic anomalies in SNC meteorites targeting siderophile elements (most notably Ni, Mo and Ru) to constrain Mars' accretion history and its formation location. △ Less

Submitted 23 May, 2018; originally announced May 2018.

Comments: Accepted for publication in GRL

Asteroid Family Associations of Active Asteroids

Authors: Henry H. Hsieh, Bojan Novakovic, Yoonyoung Kim, Ramon Brasser

Abstract: We report on the results of a systematic search for associated asteroid families for all active asteroids known to date. We find that 10 out of 12 main-belt comets (MBCs) and 5 out of 7 disrupted asteroids are linked with known or candidate families, rates that have ~0.1% and ~6% probabilities, respectively, of occurring by chance, given an overall family association rate of 37% for asteroids in t… ▽ More We report on the results of a systematic search for associated asteroid families for all active asteroids known to date. We find that 10 out of 12 main-belt comets (MBCs) and 5 out of 7 disrupted asteroids are linked with known or candidate families, rates that have ~0.1% and ~6% probabilities, respectively, of occurring by chance, given an overall family association rate of 37% for asteroids in the inner solar system. We find previously unidentified family associations between 238P/Read and the candidate Gorchakov family, 311P/PANSTARRS and the candidate Behrens family, 324P/La Sagra and the Alauda family, 354P/LINEAR and the Baptistina family, P/2013 R3-B (Catalina-PANSTARRS) and the Mandragora family, P/2015 X6 (PANSTARRS) and the Aeolia family, P/2016 G1 (PANSTARRS) and the Adeona family, and P/2016 J1-A/B (PANSTARRS) and the Theobalda family. All MBCs with family associations belong to families that contain asteroids with primitive taxonomic classifications and low average reported albedos (pV_avg < 0.10), while disrupted asteroids with family associations belong to families that contain asteroids that span wider ranges of taxonomic types and average reported albedos (0.06 < pV_avg < 0.25). These findings are consistent with MBC activity being closely correlated to composition (i.e., whether an object is likely to contain ice), while disrupted asteroid activity is not as sensitive to composition. Given our results, we describe a sequence of processes by which the formation of young asteroid families could lead to the production of present-day MBCs. △ Less

Submitted 3 January, 2018; originally announced January 2018.

Comments: 28 pages, 30 figures, accepted by AJ

Thermal effects of late accretion to the crust and mantle of Mercury

Authors: Stephen J. Mojzsis, Oleg Abramov, Elizabeth A. Frank, Ramon Brasser

Abstract: Impact bombardment on Mercury in the solar system late accretion phase (ca. 4.4 to 3.8 Ga) caused considerable mechanical, chemical and thermal reworking of its silicate reservoirs (crust and mantle). Depending on the frequency, size and velocity of such impactors, effects included regional and global scale crustal melting, and thermal perturbations of the mercurian mantle. We use a 3D transient h… ▽ More Impact bombardment on Mercury in the solar system late accretion phase (ca. 4.4 to 3.8 Ga) caused considerable mechanical, chemical and thermal reworking of its silicate reservoirs (crust and mantle). Depending on the frequency, size and velocity of such impactors, effects included regional and global scale crustal melting, and thermal perturbations of the mercurian mantle. We use a 3D transient heating model to test the effects of two bombardment scenarios on early (pre Tolstojan) Mercury mantle and crust. Results show that rare impacts by the largest (greater than 100 km diameter) bodies deliver sufficient heat to the shallow mercurian mantle producing high-temperature ultra-magnesian (komatiitic s.s.) melts. Impact heating leading to effusive (flood) volcanism can account for eponymous High-Magnesium Region (HMR) observed during the MErcury Surface, Space Environment, GEochemistry Ranging (MESSENGER) mission. We find that late accretion to Mercury induced volumetrically significant crustal melting (less than 58 vol.percent), mantle heating and melt production, which, combined with extensive resurfacing (less than 100 percent), also explains why its oldest cratering record was effectively erased, consistent with crater-counting statistics. △ Less

Submitted 14 November, 2017; originally announced November 2017.

Comments: Accepted for publication in Earth and Planetary Science Letters

The terrestrial late veneer from core disruption of a lunar-sized impactor

Authors: H. Genda, R. Brasser, S. J. Mojzsis

Abstract: Overabundances in highly siderophile elements (HSEs) of Earth's mantle can be explained by conveyance from a singular, immense (3000 km in a diameter) "Late Veneer" impactor of chondritic composition, subsequent to lunar formation and terrestrial core-closure. Such rocky objects of approximately lunar mass (about 0.01 M_E) ought to be differentiated, such that nearly all of their HSE payload is se… ▽ More Overabundances in highly siderophile elements (HSEs) of Earth's mantle can be explained by conveyance from a singular, immense (3000 km in a diameter) "Late Veneer" impactor of chondritic composition, subsequent to lunar formation and terrestrial core-closure. Such rocky objects of approximately lunar mass (about 0.01 M_E) ought to be differentiated, such that nearly all of their HSE payload is sequestered into iron cores. Here, we analyze the mechanical and chemical fate of the core of such a Late Veneer impactor, and trace how its HSEs are suspended - and thus pollute - the mantle. For the statistically most-likely oblique collision (about 45degree), the impactor's core elongates and thereafter disintegrates into a metallic hail of small particles (about 10 m). Some strike the orbiting Moon as sesquinary impactors, but most re-accrete to Earth as secondaries with further fragmentation. We show that a single oblique impactor provides an adequate amount of HSEs to the primordial terrestrial silicate reservoirs via oxidation of (<m-sized) metal particles with a hydrous, pre-impact, early Hadean Earth. △ Less

Submitted 21 September, 2017; originally announced September 2017.

Comments: Accepted for publication in Earth and Planetary Science Letters, 17 pages, 4 figures

The language of exoplanet ranking metrics needs to change

Authors: Elizabeth Tasker, Joshua Tan, Kevin Heng, Stephen Kane, David Spiegel, Ramon Brasser, Andrew Casey, Steven Desch, Caroline Dorn, Christine Houser, John Hernlund, Marine Lasbleis, Matthieu Laneuville, Anne-Sophie Libert, Lena Noack, Cayman Unterborn, June Wicks

Abstract: We have found many Earth-sized worlds but we have no way of determining if their surfaces are Earth-like. This makes it impossible to quantitatively compare habitability, and pretending we can risks damaging the field. We have found many Earth-sized worlds but we have no way of determining if their surfaces are Earth-like. This makes it impossible to quantitatively compare habitability, and pretending we can risks damaging the field. △ Less

Submitted 3 August, 2017; originally announced August 2017.

Comments: Published in Nature Astronomy, Comment, 02/2017

Journal ref: Nature Astronomy, Volume 1, id. 0042 (2017)

A colossal impact enriched Mars' mantle with noble metals

Authors: R. Brasser, S. J. Mojzsis

Abstract: Once the terrestrial planets had mostly completed their assembly, bombardment continued by planetesimals left-over from accretion. Highly siderophile element (HSE) abundances in Mars' mantle imply its late accretion supplement was 0.8 wt.%; Earth and the Moon obtained an additional 0.7 wt.% and 0.02 wt.%, respectively. The disproportionately high Earth/Moon accretion ratio is explicable by stochas… ▽ More Once the terrestrial planets had mostly completed their assembly, bombardment continued by planetesimals left-over from accretion. Highly siderophile element (HSE) abundances in Mars' mantle imply its late accretion supplement was 0.8 wt.%; Earth and the Moon obtained an additional 0.7 wt.% and 0.02 wt.%, respectively. The disproportionately high Earth/Moon accretion ratio is explicable by stochastic addition of a few remaining Ceres-sized bodies that preferentially targeted Earth. Here we show that Mars' late accretion budget also requires a colossal impact, a plausible visible remnant of which is the hemispheric dichotomy. The addition of sufficient HSEs to the martian mantle entails an impactor of at least 1200 km in diameter to have struck Mars before ca. 4430 Ma, by which time crust formation was well underway. Thus, the dichotomy could be one of the oldest geophysical features of the martian crust. Ejected debris could be the source material for its satellites. △ Less

Submitted 6 June, 2017; originally announced June 2017.

Comments: Accepted for publication in Geophysical Research Letters

N-body simulations of planet formation via pebble accretion I: First Results

Authors: Soko Matsumura, Ramon Brasser, Shigeru Ida

Abstract: Context. Planet formation with pebbles has been proposed to solve a couple of long-standing issues in the classical formation model. Some sophisticated simulations have been done to confirm the efficiency of pebble accretion. However, there has not been any global N-body simulations that compare the outcomes of planet formation via pebble accretion with observed extrasolar planetary systems. Aims.… ▽ More Context. Planet formation with pebbles has been proposed to solve a couple of long-standing issues in the classical formation model. Some sophisticated simulations have been done to confirm the efficiency of pebble accretion. However, there has not been any global N-body simulations that compare the outcomes of planet formation via pebble accretion with observed extrasolar planetary systems. Aims. In this paper, we study the effects of a range of initial parameters of planet formation via pebble accretion, and present the first results of our simulations. Methods. We incorporate the pebble accretion model by Ida et al. (2016) in the N-body code SyMBA (Duncan et al. 1998), along with the effects of gas accretion, eccentricity and inclination damping and planet migration in the disc. Results. We confirm that pebble accretion leads to a variety of planetary systems, but have difficulty in reproducing observed properties of exoplanetary systems, such as planetary mass, semimajor axis, and eccentricity distributions. The main reason behind this is a too-efficient type I migration, which sensitively depends on the disc model. However, our simulations also lead to a few interesting predictions. First, we find that formation efficiencies of planets depend on the stellar metallicities, not only for giant planets, but also for Earths (Es) and Super-Earths (SEs). The dependency for Es/SEs is subtle. Although higher metallicity environments lead to faster formation of a larger number of Es/SEs, they also tend to be lost later via dynamical instability. Second, our results indicate that a wide range of bulk densities observed for Es and SEs is a natural consequence of dynamical evolution of planetary systems. Third, the ejection trend of our simulations suggest that one free-floating E/SE may be expected for two smaller-mass planets. △ Less

Submitted 14 July, 2017; v1 submitted 11 May, 2017; originally announced May 2017.

Comments: Accepted for publication in A&A, 21 pages, 15 figures

Saving super-Earths: Interplay between pebble accretion and type I migration

Authors: R. Brasser, B. Bitsch, S. Matsumura

Abstract: Overcoming type I migration and preventing low-mass planets from spiralling into the central star is a long-studied topic. It is well known that outward migration is possible in viscous-heated discs relatively close to the central star because the entropy gradient can be sufficiently steep that the positive corotation torque overcomes the negative Lindblad torque. Yet efficiently trapping planets… ▽ More Overcoming type I migration and preventing low-mass planets from spiralling into the central star is a long-studied topic. It is well known that outward migration is possible in viscous-heated discs relatively close to the central star because the entropy gradient can be sufficiently steep that the positive corotation torque overcomes the negative Lindblad torque. Yet efficiently trapping planets in this region remains elusive. Here we study disc conditions that yield outward migration for low-mass planets under specific planet migration prescriptions. Using a steady-state disc model with a constant alpha-viscosity, outward migration is only possible when the negative temperature gradient exceeds ~0.87. We derive an implicit relation for the maximum mass at which outward migration is possible as a function of viscosity and disc scale height. We apply these criteria, using a simple power-law disc model, to planets that have reached their pebble isolation mass after an episode of rapid accretion. It is possible to trap planets with the pebble isolation mass farther than the inner edge of the disc provided that alpha_crit >~ 0.004 for discs older than 1 Myr. In very young discs the high temperature causes the planets to grow to masses exceeding the maximum for outward migration. As the disc evolves, these more massive planets often reach the central star, generally only towards the end of the disc's lifetime. Saving super-Earths is therefore a delicate interplay between disc viscosity, the opacity profile and the temperature gradient in the viscously heated inner disc. △ Less

Submitted 5 April, 2017; originally announced April 2017.

Comments: Accepted for publication in The Astronomical Journal

The cool and distant formation of Mars

Authors: R. Brasser, S. J. Mojzsis, S. Matsumura, S. Ida

Abstract: With approximately one ninth of Earth's mass, Mars is widely considered to be a stranded planetary embryo that never became a fully-grown planet. A currently popular planet formation theory predicts that Mars formed near Earth and Venus and was subsequently scattered outwards to its present location. In such a scenario, the compositions of the three planets are expected to be similar to each other… ▽ More With approximately one ninth of Earth's mass, Mars is widely considered to be a stranded planetary embryo that never became a fully-grown planet. A currently popular planet formation theory predicts that Mars formed near Earth and Venus and was subsequently scattered outwards to its present location. In such a scenario, the compositions of the three planets are expected to be similar to each other. However, bulk elemental and isotopic data for martian meteorites demonstrate that key aspects of Mars' composition are markedly different from that of Earth. This suggests that Mars formed outside of the terrestrial feeding zone during primary accretion. It is therefore probable that Mars always remained significantly farther from the Sun than Earth; its growth was stunted early and its mass remained relatively low. Here we identify a potential dynamical pathway that forms Mars in the asteroid belt and keeps it outside of Earth's accretion zone while at the same time accounting for strict age and compositional constraints, as well as mass differences. Our uncommon pathway (approximately 2% probability) is based on the Grand Tack scenario of terrestrial planet formation, in which the radial migration by Jupiter gravitationally sculpts the planetesimal disc at Mars' current location. We conclude that Mars' formation requires a specific dynamical pathway, while this is less valid for Earth and Venus. We further predict that} Mars' volatile budget is most likely different from Earth's and that Venus formed close enough to our planet that it is expected to have a nearly identical composition from common building blocks. △ Less

Submitted 1 April, 2017; originally announced April 2017.

Comments: Accepted in Earth and Planetary Science Letters

The Structure of the Distant Kuiper Belt in a Nice Model Scenario

Authors: Rosemary E. Pike, Samantha Lawler, Ramon Brasser, Cory J. Shankman, Mike Alexandersen, J. J. Kavelaars

Abstract: This work explores the orbital distribution of minor bodies in the outer Solar System emplaced as a result of a Nice model migration from the simulations of Brasser & Morbidelli (2013). This planetary migration scatters a planetesimal disk from between 29-34 AU and emplaces a population of objects into the Kuiper belt region. From the 2:1 Neptune resonance and outward, the test particles analyzed… ▽ More This work explores the orbital distribution of minor bodies in the outer Solar System emplaced as a result of a Nice model migration from the simulations of Brasser & Morbidelli (2013). This planetary migration scatters a planetesimal disk from between 29-34 AU and emplaces a population of objects into the Kuiper belt region. From the 2:1 Neptune resonance and outward, the test particles analyzed populate the outer resonances with orbital distributions consistent with trans-Neptunian object (TNO) detections in semi-major axis, inclination, and eccentricity, while capture into the closest resonances is too efficient. The relative populations of the simulated scattering objects and resonant objects in the 3:1 and 4:1 resonances are also consistent with observed populations based on debiased TNO surveys, but the 5:1 resonance is severely underpopulated compared to population estimates from survey results. Scattering emplacement results in the expected orbital distribution for the majority of the TNO populations, however the origin of the large observed population in the 5:1 resonance remains unexplained. △ Less

Submitted 24 January, 2017; originally announced January 2017.

Comments: 15 pages, 5 figures

How primordial is the structure of comet 67P/C-G? Combined collisional and dynamical models suggest a late formation

Authors: M. Jutzi, W. Benz, A. Toliou, A. Morbidelli, R. Brasser

Abstract: There is an active debate about whether the properties of comets as observed today are primordial or, alternatively, if they are a result of collisional evolution or other processes. We investigate the effects of collisions on a comet with a structure like 67P/C-G. We develop scaling laws for the critical specific impact energies required for a significant shape alteration. These are then used in… ▽ More There is an active debate about whether the properties of comets as observed today are primordial or, alternatively, if they are a result of collisional evolution or other processes. We investigate the effects of collisions on a comet with a structure like 67P/C-G. We develop scaling laws for the critical specific impact energies required for a significant shape alteration. These are then used in simulations of the combined dynamical and collisional evolution of comets in order to study the survival probability of a primordially formed object with a shape like 67P/C-G. The effects of impacts on comet 67P/C-G are studied using a SPH shock physics code. The resulting critical specific impact energy defines a minimal projectile size which is used to compute the number of shape-changing collisions in a set of dynamical simulations. These simulations follow the dispersion of the trans-Neptunian disk during the giant planet instability, the formation of a scattered disk, and produce 87 objects that penetrate into the inner solar system with orbits consistent with the observed JFC population. The collisional evolution before the giant planet instability is not considered here. Hence, our study is conservative in its estimation of the number of collisions. We find that in any scenario considered here, comet 67P/C-G would have experienced a significant number of shape-changing collisions, if it formed primordially. This is also the case for generic bi-lobe shapes. Our study also shows that impact heating is very localized and that collisionally processed bodies can still have a high porosity. Our study indicates that the observed bi-lobe structure of comet 67P/C-G may not be primordial, but might have originated in a rather recent event, possibly within the last 1 Gy. This may be the case for any kilometer-sized two-component cometary nuclei. △ Less

Submitted 8 November, 2016; originally announced November 2016.

Comments: Astronomy & Astrophysics, accepted pending minor revisions

Journal ref: A&A 597, A61 (2017)

Late veneer and late accretion to the terrestrial planets

Authors: R. Brasser, S. J. Mojzsis, S. C. Werner, S. Matsumura, S. Ida

Abstract: It is generally accepted that silicate-metal (`rocky') planet formation relies on coagulation from a mixture of sub-Mars sized planetary embryos and (smaller) planetesimals that dynamically emerge from the evolving circum-solar disc in the first few million years of our Solar System. Once the planets have, for the most part, assembled after a giant impact phase, they continue to be bombarded by a… ▽ More It is generally accepted that silicate-metal (`rocky') planet formation relies on coagulation from a mixture of sub-Mars sized planetary embryos and (smaller) planetesimals that dynamically emerge from the evolving circum-solar disc in the first few million years of our Solar System. Once the planets have, for the most part, assembled after a giant impact phase, they continue to be bombarded by a multitude of planetesimals left over from accretion. Here we place limits on the mass and evolution of these planetesimals based on constraints from the highly siderophile element (HSE) budget of the Moon. Outcomes from a combination of N-body and Monte Carlo simulations of planet formation lead us to four key conclusions about the nature of this early epoch. First, matching the terrestrial to lunar HSE ratio requires either that the late veneer on Earth consisted of a single lunar-size impactor striking the Earth before 4.45 Ga, or that it originated from the impact that created the Moon. An added complication is that analysis of lunar samples indicates the Moon does not preserve convincing evidence for a late veneer like Earth. Second, the expected chondritic veneer component on Mars is 0.06 weight percent. Third, the flux of terrestrial impactors must have been low ( <=10^(-6) M_earth/Myr) to avoid wholesale melting of Earth's crust after 4.4~Ga, and to simultaneously match the number of observed lunar basins. This conclusion leads to an Hadean eon which is more clement than assumed previously. Last, after the terrestrial planets had fully formed, the mass in remnant planetesimals was ~10^(-3) M_earth, lower by at least an order of magnitude than most previous models suggest. Our dynamically and geochemically self-consistent scenario requires that future N-body simulations of rocky planet formation either directly incorporate collisional grinding or rely on pebble accretion. △ Less

Submitted 6 September, 2016; originally announced September 2016.

Comments: Accepted for publication in Earth and Planetary Science Letters

Analysis of terrestrial planet formation by the Grand Tack model: System architecture and tack location

Authors: R. Brasser, S. Matsumura, S. Ida, S. J. Mojzsis, S. C. Werner

Abstract: The Grand Tack model of terrestrial planet formation has emerged in recent years as the premier scenario used to account for several observed features of the inner solar system. It relies on early migration of the giant planets to gravitationally sculpt and mix the planetesimal disc down to ~1 AU, after which the terrestrial planets accrete from material left in a narrow circum-solar annulus. Here… ▽ More The Grand Tack model of terrestrial planet formation has emerged in recent years as the premier scenario used to account for several observed features of the inner solar system. It relies on early migration of the giant planets to gravitationally sculpt and mix the planetesimal disc down to ~1 AU, after which the terrestrial planets accrete from material left in a narrow circum-solar annulus. Here we have investigated how the model fares under a range of initial conditions and migration course-change (`tack') locations. We have run a large number of N-body simulations with a tack location of 1.5 AU and 2 AU and tested initial conditions using equal mass planetary embryos and a semi-analytical approach to oligarchic growth. We make use of a recent model of the protosolar disc that takes account of viscous heating, include the full effect of type 1 migration, and employ a realistic mass-radius relation for the growing terrestrial planets. Results show that the canonical tack location of Jupiter at 1.5 AU is inconsistent with the most massive planet residing at 1 AU at greater than 95% confidence. This favours a tack farther out at 2 AU for the disc model and parameters employed. Of the different initial conditions, we find that the oligarchic case is capable of statistically reproducing the orbital architecture and mass distribution of the terrestrial planets, while the equal mass embryo case is not. △ Less

Submitted 3 March, 2016; originally announced March 2016.

Comments: Accepted for publication in The Astrophysical Journal

Effects of Dynamical Evolution of Giant Planets on the Delivery of Atmophile Elements During Terrestrial Planet Formation

Authors: Soko Matsumura, Ramon Brasser, Shigeru Ida

Abstract: Recent observations started revealing the compositions of protostellar discs and planets beyond the Solar System. In this paper, we explore how the compositions of terrestrial planets are affected by dynamical evolution of giant planets. We estimate the initial compositions of building blocks of these rocky planets by using a simple condensation model, and numerically study the compositions of pla… ▽ More Recent observations started revealing the compositions of protostellar discs and planets beyond the Solar System. In this paper, we explore how the compositions of terrestrial planets are affected by dynamical evolution of giant planets. We estimate the initial compositions of building blocks of these rocky planets by using a simple condensation model, and numerically study the compositions of planets formed in a few different formation models of the Solar System. We find that the abundances of refractory and moderately volatile elements are nearly independent of formation models, and that all the models could reproduce the abundances of these elements of the Earth. The abundances of atmophile elements, on the other hand, depend on the scattering rate of icy planetesimals into the inner disc as well as the mixing rate of the inner planetesimal disc. For the classical formation model, neither of these mechanisms are efficient and the accretion of atmophile elements during the final assembly of terrestrial planets appears to be difficult. For the Grand Tack model, both of these mechanisms are efficient, which leads to a relatively uniform accretion of atmophile elements in the inner disc. It is also possible to have a "hybrid" scenario where the mixing is not very efficient but the scattering is efficient. The abundances of atmophile elements in this case increases with orbital radii. Such a scenario may occur in some of the extrasolar planetary systems which are not accompanied by giant planets or those without strong perturbations from giants. We also confirm that the Grand Tack scenario leads to the distribution of asteroid analogues where rocky planetesimals tend to exist interior to icy ones, and show that their overall compositions are consistent with S-type and C-type chondrites, respectively. △ Less

Submitted 27 December, 2015; originally announced December 2015.

Comments: 21 pages, 16 figures, 6 tables, accepted for publication in ApJ

Tilting Saturn without tilting Jupiter: Constraints on giant planet migration

Authors: R. Brasser, Man Hoi Lee

Abstract: The migration and encounter histories of the giant planets in our Solar System can be constrained by the obliquities of Jupiter and Saturn. We have performed secular simulations with imposed migration and N-body simulations with planetesimals to study the expected obliquity distribution of migrating planets with initial conditions resembling those of the smooth migration model, the resonant Nice m… ▽ More The migration and encounter histories of the giant planets in our Solar System can be constrained by the obliquities of Jupiter and Saturn. We have performed secular simulations with imposed migration and N-body simulations with planetesimals to study the expected obliquity distribution of migrating planets with initial conditions resembling those of the smooth migration model, the resonant Nice model and two models with five giant planets initially in resonance (one compact and one loose configuration). For smooth migration, the secular spin-orbit resonance mechanism can tilt Saturn's spin axis to the current obliquity if the product of the migration time scale and the orbital inclinations is sufficiently large (exceeding 30 Myr deg). For the resonant Nice model with imposed migration, it is difficult to reproduce today's obliquity values, because the compactness of the initial system raises the frequency that tilts Saturn above the spin precession frequency of Jupiter, causing a Jupiter spin-orbit resonance crossing. Migration time scales sufficiently long to tilt Saturn generally suffice to tilt Jupiter more than is observed. The full N-body simulations tell a somewhat different story, with Jupiter generally being tilted as often as Saturn, but on average having a higher obliquity. The main obstacle is the final orbital spacing of the giant planets, coupled with the tail of Neptune's migration. The resonant Nice case is barely able to simultaneously reproduce the {orbital and spin} properties of the giant planets, with a probability ~0.15%. The loose five planet model is unable to match all our constraints (probability <0.08%). The compact five planet model has the highest chance of matching the orbital and obliquity constraints simultaneously (probability ~0.3%). △ Less

Submitted 22 September, 2015; originally announced September 2015.

Comments: Accepted for publication in the Astronomical Journal