Source author record

R. Brasser

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

22works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

22 published item(s)

preprint2022arXiv

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

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

preprint2022arXiv

Long-term tidal evolution of the TRAPPIST-1 system

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.

preprint2021arXiv

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

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

preprint2021arXiv

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

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.

preprint2021arXiv

Formation of Venus, Earth and Mars: Constrained by isotopes

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.

preprint2020arXiv

No evidence for interstellar planetesimals trapped in the Solar System

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.

preprint2019arXiv

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

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.

preprint2016arXiv

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

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.

preprint2016arXiv

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

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.

preprint2016arXiv

Late veneer and late accretion to the terrestrial planets

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.

preprint2015arXiv

Tilting Saturn without tilting Jupiter: Constraints on giant planet migration

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%).

preprint2014arXiv

A dynamical study on the habitability of terrestrial exoplanets II: The super Earth HD 40307 g

HARPS and it Kepler results indicate that half of solar-type stars host planets with periods P<100 d and masses M < 30 M_E. These super Earth systems are compact and dynamically cold. Here we investigate the stability of the super Earth system around the K-dwarf HD40307. It could host up to six planets, with one in the habitable zone. We analyse the system's stability using numerical simulations from initial conditions within the observational uncertainties. The most stable solution deviates 3.1 sigma from the published value, with planets e and f not in resonance and planets b and c apsidally aligned. We study the habitability of the outer planet through the yearly-averaged insolation and black-body temperature at the pole. Both undergo large variations because of its high eccentricity and are much more intense than on Earth. The insolation variations are precession dominated with periods of 40 kyr and 102 kyr for precession and obliquity if the rotation period is 3 d. A rotation period of about 1.5 d could cause extreme obliquity variations because of capture in a Cassini state. For faster rotation rates the periods converge to 10 kyr and 20 kyr. The large uncertainty in the precession period does not change the overall outcome.

preprint2014arXiv

An Oort cloud origin of the Halley-type comets

The origin of the Halley-type comets (HTCs) is one of the last mysteries of the dynamical evolution of the Solar System. Prior investigation into their origin has focused on two source regions: the Oort cloud and the Scattered Disc. From the former it has been difficult to reproduce the non-isotropic, prograde skew in the inclination distribution of the observed HTCs without invoking a multi-component Oort cloud model and specific fading of the comets. The Scattered Disc origin fares better but suffers from needing an order of magnitude more mass than is currently advocated by theory and observations. Here we revisit the Oort cloud origin and include cometary fading. Our observational sample stems from the JPL catalogue. We only keep comets discovered and observed after 1950 but place no a priori restriction on the maximum perihelion distance of observational completeness. We then numerically evolve half a million comets from the Oort cloud through the realm of the giant planets and keep track of their number of perihelion passages with perihelion distance q<2.5AU, below which the activity is supposed to increase considerably. We can simultaneously fit the HTC inclination and semi-major axis distribution very well with a power law fading function of the form m^-k, where m is the number of perihelion passages with q<2.5 AU and k is the fading index. We match both the inclination and semi-major axis distributions when k~1 and the maximum imposed perihelion distance of the observed sample is q~1.8AU. The value of k is higher than the one obtained for the Long-Period Comets (LPCs), with k~0.7. This increase in k is most likely the result of cometary surface processes. We argue the HTC sample is now most likely complete for q<1.8AU. We calculate that the steady-state number of active HTCs with diameter D>2.3km and q<1.8AU is of the order of 100.

preprint2014arXiv

An updated estimate of the number of Jupiter-family comets using a simple fading law

It has long been hypothesised that the Jupiter-family comets (JFCs) come from the scattered disc, an unstable planetesimal population beyond Neptune. This viewpoint has been widely accepted, but a few issues remain, the most prominent of which are the total number of visible JFCs with a perihelion distance q<2.5 AU and the corresponding number of objects in the scattered disc. In this work we give a robust estimate of the number of visible JFCs with q<2.5 AU and diameter D>2.3 km based on recent observational data. This is combined with numerical simulations that use a simple fading law applied to JFCs that come close to the Sun. For this we numerically evolve thousands of comets from the scattered disc through the realm of the giant planets and keep track of their number of perihelion passages with perihelion distance q<2.5 AU, below which the activity is supposed to increase considerably. We can simultaneously fit the JFC inclination and semi-major axis distribution accurately with a delayed power law fading function and the maximum perihelion distance below which the observational data is complete is q~2.3 AU. From observational data we calculate that a JFC with diameter D=2.3 km has a typical total absolute magnitude H_T=10.8, and the steady-state number of active JFCs with diameter D>2.3 km and q<2.5 AU is of the order of 300 (but with large uncertainties), approximately a factor two higher than earlier estimates. The increased JFC population results in a scattered disc population of 6 billion objects and decreases the observed Oort cloud to scattered disc population ratio to 13, virtually the same as the value of 12 obtained with numerical simulations.

preprint2014arXiv

Reassessing the formation of the Inner Oort cloud in an embedded star cluster II: Probing the inner edge

The detached object Sedna is likely at the inner edge of the Oort cloud, more precisely the inner Oort cloud (IOC). Until recently it was the sole member of this population. The recent discovery of the detached object 2012 VP113 has confirmed that there should be more objects in this region. Three additional IOC candidates with orbits much closer to Neptune have been proposed in the past decade since Sedna's discovery: 2000 CR105, 2004 VN112 and 2010 GB174. Sedna and 2012 VP113 have perhelia near 80 AU and semi-major axes over 250 AU. The latter three have perihelia between 44 AU and 50 AU and semi-major axes between 200 AU and 400 AU. Here we determine whether the latter three objects belong to the IOC or are from the Kuiper Belt's Extended Scattered Disc (ESD) using numerical simulations. We assume that the IOC was formed when the Sun was in its birth cluster. We analyse the evolution of the IOC and the Scattered Disc (SD) during an episode of late giant planet migration. We examine the impact of giant planet migration in the context of four and five planets. We report that the detached objects 2004 VN112 and 2010 GB174 are likely members of the IOC that were placed there while the Sun was in its birth cluster or during an episode of Solar migration in the Galaxy. The origin of 2000 CR105 is ambiguous but it is likely it belongs to the ESD. Based on our simulations we find that the maximum perihelion distance of SD objects is 41 AU when the semi-major axis is higher than 250 AU. Objects closer in are subject to mean motion resonances with Neptune that may raise their perihelia. The five planet model yields the same outcome. We impose a conservative limit and state that all objects with perihelion distance q>45 AU and semi-major axis a>250 AU belong to the inner Oort cloud.

preprint2013arXiv

Constraining the primordial orbits of the Terrestrial Planets

Evidence in the Solar System suggests that the giant planets underwent an epoch of radial migration that was very rapid, with an e-folding timescale shorter than 1~Myr. It is probable that the cause of this migration was that the giant planets experienced an orbital instability that caused them to encounter each other, resulting in radial migration. Several works suggest that this dynamical instability occurred `late', long after all the planets had formed and the solar nebula had dissipated. Assuming that the terrestrial planets had already formed, then their orbits would have been affected by the migration of the giant planets. As a result, how did the orbits of the terrestrial planets change? And can we use this migration to obtain information on the primordial orbits of the terrestrial planets? We directly model a large number of terrestrial planet systems and their response to giant planet migration. We study the change in the Angular Momentum Deficit (AMD) of the terrstrials. We conclude that the primordial AMD should have been lower than ~70\% of the current value, but higher than 10\%. We find that a scenario with five giant planets better satisfies the orbital constraints of the terrestrial planets. We predict that Mars was initially on an eccentric and inclined orbit while the orbits of Mercury, Venus and Earth were more circular and coplanar. The lower primordial dynamical excitement and the peculiar partitioning between planets impose new constraints for terrestrial planet formation simulations.

preprint2013arXiv

Oort Cloud and Scattered Disc formation during a late dynamical instability in the Solar System

One of the outstanding problems of the dynamical evolution of the outer solar system concerns the observed population ratio between the Oort Cloud (OC) and the Scattered Disc (SD): observations suggest that this ratio lies between 100 and 1000 but simulations that produce these two reservoirs simultaneously consistently yield a value of the order of 10. Here we stress that the populations in the OC and SD are inferred from the observed fluxes of new Long Period Comets (LPCs) and Jupiter-family comets (JFCs), brighter than some reference total magnitude. However, the population ratio estimated in the simulations of formation of the SD and OC refers to objects bigger than a given size. There are multiple indications that LPCs are intrinsically brighter than JFCs, i.e. an LPC is smaller than a JFC with the same total absolute magnitude. When taking this into account we revise the SD/JFC population ratio from our simulations relative to Duncan and Levison (1997), and then deduce from the observations that the size-limited population ratio between the OC and the SD is 44 (-34)(+54). From simulations we obtain 12 +- 1 but the agreement cannot be rejected by the null hypothesis.

preprint2012arXiv

A dynamical study on the habitability of terrestrial exoplanets I: Tidally evolved planet-satellite pairs

We investigate the obliquity and spin period of Earth-Moon like systems after 4.5 Gyr of tidal evolution with various satellite masses and initial planetary obliquity and discuss their relations to the habitability of the planet. We find three possible outcomes: either i) the system is still evolving, ii) the system is double synchronous or iii) the satellite has collided with the planet. The transition between case i) and ii) is abrupt and occurs at slightly larger satellite mass ($m_s \sim 0.02m_p$) than the lunar mass. We suggest that cases ii) and iii) are less habitable than case i). Using results from models of giant impacts and satellite accretion, we found that the systems that mimic our own with rotation period $12 < P_p < 48$ h and current planetary obliquity $\varepsilon_p < 40^\circ$ or $\varepsilon_p > 140^\circ$ only represent 14% of the possible outcomes. Elser et al. (2011) conclude that the probability of a terrestrial planet having a heavy satellite is 13%. Combining these results suggests that the probability of ending up with a system such as our own is of the order of 2%.

preprint2012arXiv

Considerations on the accretion of Uranus and Neptune by mutual collisions of planetary embryos in the vicinity of Jupiter and Saturn

Modeling the formation of the ice giants Uranus and Neptune is a long-lasting problem in planetary science. Due to gas-drag, collisional damping, and resonant shepherding, the planetary embryos repel the planetesimals away from their reach and thus they stop growing (Levison et al. 2010). This problem persists independently of whether the accretion took place at the current locations of the ice giants or closer to the Sun. Instead of trying to push the runaway/oligarchic growth of planetary embryos up to 10-15 Earth masses, we envision the possibility that the planetesimal disk could generate a system of planetary embryos of only 1-3 Earth masses. Then we investigate whether these embryos could have collided with each other and grown enough to reach the masses of current Uranus and Neptune. Our results point to two major problems. First, there is typically a large difference in mass between the first and the second most massive core formed and retained beyond Saturn. Second, in many simulations the final planetary system has more than two objects beyond Saturn. The growth of a major planet from a system of embryos requires strong damping of eccentricities and inclinations from the disk of gas. But strong damping also favors embryos and cores to find a stable resonant configuration, so that systems with more than two surviving objects are found. In addition to these problems, in order to have substantial mutual accretion among embryos, it is necessary to assume that the surface density of the gas was several times higher than that of the minimum-mass solar nebula. However this contrasts with the common idea that Uranus and Neptune formed in a gas-starving disk, which is suggested by the relatively small amount of hydrogen and helium contained in the atmospheres of these planets. Only one of our simulations "by chance" successfully reproduced the structure of the outer Solar System.

preprint2011arXiv

An Oort cloud origin for the high-inclination, high-perihelion Centaurs

We analyse the origin of three Centaurs with perihelia in the range 15 AU to 30 AU, inclinations above 70 deg and semi-major axes shorter than 100 AU. Based on long-term numerical simulations we conclude that these objects most likely originate from the Oort cloud rather than the Kuiper Belt or Scattered Disc. We estimate that there are currently between 1 and 200 of these high-inclination, high-perihelion Centaurs with absolute magnitude H<8.

preprint2011arXiv

Reassessing the formation of the inner Oort cloud in an embedded star cluster

We re-examine the formation of the inner Oort comet cloud while the Sun was in its birth cluster with the aid of numerical simulations. This work is a continuation of an earlier study (Brasser et al., 2006) with several substantial modifications. First, the system consisting of stars, planets and comets is treated self-consistently in our N-body simulations, rather than approximating the stellar encounters with the outer Solar System as hyperbolic fly-bys. Second, we have included the expulsion of the cluster gas, a feature that was absent previously. Third, we have used several models for the initial conditions and density profile of the cluster -- either a Hernquist or Plummer potential -- and chose other parameters based on the latest observations of embedded clusters from the literature. {These other parameters result in the stars being on radial orbits and the cluster collapses.} Similar to previous studies, in our simulations the inner Oort cloud is formed from comets being scattered by Jupiter and Saturn and having their pericentres decoupled from the planets by perturbations from the cluster gas and other stars. We find that all inner Oort clouds formed in these clusters have an inner edge ranging from 100 AU to a few hundred AU, and an outer edge at over 100\,000 AU, with little variation in these values for all clusters. All inner Oort clouds formed are consistent with the existence of (90377) Sedna, an inner Oort cloud dwarf planetoid, at the inner edge of the cloud: Sedna tends to be at the innermost 2% for Plummer models, while it is 5% for Hernquist models. We emphasise that the existence of Sedna is a generic outcome. We define a `concentration radius' for the inner Oort cloud and find that its value increases with increasing number of stars in the cluster, ranging from 600 AU to 1500 AU for Hernquist clusters and from 1500 AU to 4000 AU for Plummer clusters...

preprint2011arXiv

Reassessing the origin of Triton

Agnor & Hamilton (2006) demonstrated that the disruption of a binary was an effective mechanism to capture Triton. The subsequent evolution of Triton's post-capture orbit could have proceeded through gravitational tides. The study by Agnor & Hamilton (2006) is repeated in the framework of the Nice model to determine the post-capture orbit of Triton. After capture it is then subjected to tidal evolution. The perturbations from the Sun and the figure of Neptune are included. The perturbations from the Sun acting on Triton cause it to spend a long time in its high-eccentricity phase, usually of the order of 10 Myr, while the typical time to circularise to its current orbit is some 200 Myr. The current orbit of Triton is consistent with an origin through binary capture and tidal evolution, even though the model prefers Triton to be closer to Neptune than it is today. The probability of capturing Triton in this manner is approximately 0.7%. Since the capture of Triton was at most a 50% event -- since only Neptune has one, but Uranus does not -- we deduce that in the primordial trans-Neptunian disc there were 100 binaries with at least one Triton-sized member. Morbidelli et al. (2009) concludes there were some 1000 Triton-sized bodies in the trans-Neptunian proto-planetary disc, so the primordial binary fraction with at least one Triton-sized member is 10%. This value is consistent with theoretical predictions, but at the low end. If Triton was captured at the same time as Neptune's irregular satellites, the far majority of these, including Nereid, would be lost. This suggests either that Triton was captured on an orbit with a small semi-major axis a < 50 R_N (a rare event), or that it was captured before the dynamical instability of the Nice model, or that some other mechanism was at play. The issue of keeping the irregular satellites remains unresolved.