Source author record

A. Morbidelli

A. Morbidelli 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

27works
3topics
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

27 published item(s)

preprint2024arXiv

Dynamics of the Beta Pictoris planetary system and its falling evaporating bodies

For decades, the spectral variations of Beta Pictoris have been modelled as the result of the evaporation of exocomets close to the star, termed falling evaporating bodies (FEBs). Resonant perturbations by a giant planet have been proposed to explain the dynamical origin of these stargrazers. The disk is now known to harbour two giant planets, Beta Pic b and c, orbiting the star at 9.9 au and 2.7 au. While the former almost matches the planet formerly suspected, the discovery of the latter complicates the picture. We first question the stability of the two-planet system. Then we investigate the dynamics of a disk of planetesimals orbiting the star with both planets to check the validity of the FEB generation mechanism. Symplectic N-body simulations are used to determine which regions of the planetesimal disk are dynamically stable. Then we focus on regions where disk particles are able to reach high eccentricities thanks to resonant mechanisms. The first result is that the system is dynamically stable. Both planets may temporarily fall in 7:1 mean motion resonance (MMR). Then, simulations reveal that the whole region extending between ~1.5 au and ~25 au is unstable to planetary perturbations. However, a disk below 1.5 au survives, which appears to constitute an active source of FEBs via high-order MMRs with Beta Pic c. Beta Pic b acts as a distant perturber that helps sustain the whole process. These simulations rule out the preceding FEB generation mechanism model, which placed their origin at around 4-5 au. Conversely, FEBs are likely to originate from a region much further in and related to MMRs with Beta Pic c. That mechanism also appears to last longer, as new planetesimals are able to continuously enter the MMRs and evolve towards the FEB state. Subsequently, the physical nature of the FEBs may differ from that previously thought, and presumably may not be icy.

preprint2022arXiv

The Metal-Silicate Partitioning of Carbon During Earth's Accretion and its Distribution in the Early Solar System

Carbon is an essential element for the existence and evolution of life on Earth. Its abundance in Earth's crust and mantle (the Bulk Silicate Earth, BSE) is surprisingly high given that carbon is strongly siderophile (metal-loving) at low pressures and temperatures, and hence should have segregated almost completely into Earth's core during accretion. Estimates of the concentration of carbon in the BSE lie in the range 100-260 ppm and are much higher than expected based on simple models of core-mantle differentiation. Here we show through experiments at the putative conditions of Earth's core formation (49-71 GPa and 3600-4000 K) that carbon is significantly less siderophile at these conditions than at the low pressures ($\leq$13 GPa) and temperatures ($\leq$2500 K) of previous large volume press studies, but at least an order of magnitude more siderophile than proposed recently based on an experimental approach that is similar to ours. Using our new data along with previously published results, we derive a new parameterization of the pressure-temperature dependence of the metal-silicate partitioning of carbon. We apply this parameterization in a model that combines planet formation and core-mantle differentiation that is based on astrophysical N-body accretion simulations. Because differentiated planetesimals were almost completely depleted in carbon due to sublimation at high temperatures, almost all carbon in the BSE was added by the accretion of fully-oxidized carbonaceous chondrite material from the outer solar system. Carbon is added to the mantle continuously throughout accretion and its concentration reaches values within the BSE range (e.g. 140$\pm$40 ppm) at the end of accretion. The corresponding final core and bulk Earth carbon concentrations are 1270$\pm$300 ppm and 495$\pm$125 ppm respectively.

preprint2021arXiv

Migration of Jupiter mass planets in discs with laminar accretion flows

Migration of giant planets in discs with low viscosity has been studied recently. The proportionality between migration speed and the disc's viscosity is broken by the presence of vortices that appear at the edges of the planet-induced gap. Our goal is to investigate vortex-driven migration in low-viscosity discs in the presence of radial advection of gas, as expected from angular momentum removal by magnetised disc winds. We performed three dimensional simulations using the grid-based code FARGOCA. We mimicked the effects of a disc wind by applying a synthetic torque on a surface layer of the disc characterised by a prescribed column density Sigma_A so that it results in a disc accretion rate of 10^-8 Solar masses per year. Discs with this structure are called 'layered' and the layer where the torque is applied is denoted as 'active'. We also consider the case of accretion focussed near the disc midplane to mimic transport properties induced by a large Hall effect or by weak Ohmic diffusion. We observe two migration phases: in the first phase, the migration of the planet is driven by the vortex and is directed inwards. This phase ends when the vortex disappear. Migration depends on the ability of the torque from the planet to block the accretion flow. When the flow is fast and unimpeded (small Sigma_A) migration is very slow. When the accretion flow is completely blocked, migration is faster and the speed is controlled by the rate at which the accretion flow refills the gap behind the migrating planet. The migration speed of a giant planet in a layered protoplanetary disc depends on the thickness of the accreting layer. The lack of large-scale migration apparently experienced by the majority of giant exoplanets can be explained if the accreting layer is sufficiently thin to allow unimpeded accretion through the disc.

preprint2020arXiv

A re-assessment of the Kuiper belt size distribution for sub-kilometer objects, revealing collisional equilibrium at small sizes

We combine several constraints provided by the crater records on Arrokoth and the worlds of the Pluto system to compute the size-frequency distribution (SFD) of the crater production function for craters with diameter D<10km. For this purpose, we use a Kuiper belt objects (KBO) population model calibrated on telescopic surveys, that describes also the evolution of the KBO population during the early Solar System. We further calibrate this model using the crater record on Pluto, Charon and Nix. Using this model, we compute the impact probability on Arrokoth, integrated over the age of the Solar System. This probability is then used together with other observational constraints to determine the slope of the crater-production function on Arrokoth. In addition, we use our Kuiper belt model also to compare the impact rates and velocities of KBOs on Arrokoth with those on Charon, integrated over the crater retention ages of their respective surfaces. This allows us to establish a relationship between the spatial density of sub-km craters on Arrokoth and of D~20km craters on Charon. Together, all these considerations suggest the crater production function on these worlds has a cumulative power law slope of -1.5<q<-1.2. Converted into a projectile SFD slope, we find -1.2<q_KBO<-1.0. These values are close to the cumulative slope of main belt asteroids in the 0.2--2~km range, a population in collisional equilibrium. For KBOs, however, this slope appears to extend from ~2km down to objects a few tens of meters in diameter, as inferred from sub-km craters on Arrokoth. From the measurement of the dust density in the Kuiper belt made by the New Horizons mission, we predict that the SFD of the KBOs becomes steep again below ~10-30m. All these considerations strongly indicate that the size distribution of the KBO population is in collisional equilibrium.

preprint2020arXiv

Migration of Jupiter mass planets in low viscosity discs

Type-II migration of giant planets has a speed proportional to the disc's viscosity for values of the alpha viscosity parameter larger than 1.e-4 . At lower viscosities previous studies, based on 2D simulations have shown that migration can be very chaotic and often characterized by phases of fast migration. The reason is that in low-viscosity discs vortices appear due to the Rossby-wave instability at the edges of the gap opened by the planet. Migration is then determined by vortex-planet interactions. Our aim is to study migration in low viscosity 3D discs. We performed numerical simulations using 2D (including self-gravity) and 3D codes. After selecting disc masses for which self-gravity is not important, 3D simulations without self-gravity can be safely used. In our nominal simulation only numerical viscosity is present. We then performed simulations with prescribed viscosity to assess the threshold below which the new migration processes appear. We show that for alpha viscosity <= 1.e-5 two migration modes are possible which differ from classical Type-II migration, in the sense that they are not proportional to the disc's viscosity. The first occurs when the gap opened by the planet is not very deep. This occurs in 3D simulations and/or when a big vortex forms at the outer edge of the planetary gap, diffusing material into the gap. We call this type of migration "vortex-driven migration". This migration is very slow and cannot continue indefinitely, because eventually the vortex dissolves. The second migration mode occurs when the gap is deep so that the planet's eccentricity grows to a value ~0.2 due to inefficient eccentricity damping by corotation resonances. This second, faster migration mode appears to be typical of 2D models in discs with slower damping of temperature's perturbations.

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.

preprint2020arXiv

Subsolar Al/Si and Mg/Si ratios of non-carbonaceous chondrites reveal planetesimal formation during early condensation in the protoplanetary disk

The Al/Si and Mg/Si ratios in non-carbonaceous chondrites are lower than the solar (i.e., CI-chondritic) values, in sharp contrast to the non-CI carbonaceous meteorites and the Earth, which are enriched in refractory elements and have Mg/Si ratios that are solar or larger. We show that the formation of a first generation of planetesimals during the condensation of refractory elements implies the subsequent formation of residual condensates with strongly sub-solar Al/Si and Mg/Si ratios. The mixing of residual condensates with different amounts of material with solar refractory/Si element ratios explains the Al/Si and Mg/Si values of non-carbonaceous chondrites. To match quantitatively the observed ratios, we find that the first-planetesimals should have accreted when the disk temperature was ~1,330-1,400 K depending on pressure and assuming a solar C/O ratio of the disk. We discuss how this model relates to our current understanding of disk evolution, grain dynamics, and planetesimal formation. We also extend the discussion to moderately volatile elements (e.g., Na), explaining how it may be possible that the depletion of these elements in non-carbonaceous chondrites is correlated with the depletion of refractory elements (e.g., Al). Extending the analysis to Cr, we find evidence for a higher than solar C/O ratio in the protosolar disk's gas when/where condensation from a fractionated gas occurred. Finally, we discuss the possibility that the supra-solar Al/Si and Mg/Si ratios of the Earth are due to the accretion of ~40% of the mass of our planet from the first-generation of refractory-rich planetesimals.

preprint2016arXiv

Circumplanetary disk or circumplanetary envelope?

We present three-dimensional simulations with nested meshes of the dynamics of the gas around a Jupiter mass planet with the JUPITER and FARGOCA codes. We implemented a radiative transfer module into the JUPITER code to account for realistic heating and cooling of the gas. We focus on the circumplanetary gas flow, determining its characteristics at very high resolution ($80\%$ of Jupiter's diameter). In our nominal simulation where the temperature evolves freely by the radiative module and reaches 13000 K at the planet, a circumplanetary envelope was formed filling the entire Roche-lobe. Because of our equation of state is simplified and probably overestimates the temperature, we also performed simulations with limited maximal temperatures in the planet region (1000 K, 1500 K, and 2000 K). In these fixed temperature cases circumplanetary disks (CPDs) were formed. This suggests that the capability to form a circumplanetary disk is not simply linked to the mass of the planet and its ability to open a gap. Instead, the gas temperature at the planet's location, which depends on its accretion history, plays also fundamental role. The CPDs in the simulations are hot and cooling very slowly, they have very steep temperature and density profiles, and are strongly sub-Keplerian. Moreover, the CPDs are fed by a strong vertical influx, which shocks on the CPD surfaces creating a hot and luminous shock-front. In contrast, the pressure supported circumplanetary envelope is characterized by internal convection and almost stalled rotation.

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.

preprint2015arXiv

Fossilized condensation lines in the Solar System protoplanetary disk

The terrestrial planets and the asteroids dominant in the inner asteroid belt are water poor. However, in the protoplanetary disk the temperature should have decreased below water condensation level well before the disk was photoevaporated. Thus, the global water depletion of the inner Solar System is puzling. We show that, even if the inner disk becomes cold, there cannot be direct condensation of water. This is because the snowline moves towards the Sun more slowly than the gas itself. The appearance of ice in a range of heliocentric distances swept by the snowline can only be due to the radial drift of icy particles from the outer disk. However, if a sufficiently massive planet is present, the radial drift of particles is interrupted, because the disk acquires a superKeplerian rotation just outside of the planetary orbit. From this result, we propose that the precursor of Jupiter achieved about 20 Earth masses when the snowline was still around 3 AU. This effectively fossilized the snowline at that location. Although cooling, the disk inside of the Jovian orbit remained ice-depleted because the flow of icy particles from the outer system was intercepted by the planet. This scenario predicts that planetary systems without giant planets should be much more rich in water in their inner regions than our system. We also show that the inner edge of the planetesimal disk at 0.7AU, required in terrestrial planet formation models to explain the small mass of Mercury and the absence of planets inside of its orbit, could be due to the silicate condensation line, fossilized at the end of the phase of streaming instability that generated the planetesimal seeds. Thus, when the disk cooled, silicate particles started to drift inwards of 0.7AU without being sublimated, but they could not be accreted by any pre-existing planetesimals.

preprint2015arXiv

Outwards migration for planets in stellar irradiated 3D discs

For the very first time we present 3D simulations of planets embedded in stellar irradiated discs. It is well known that thermal effects could reverse the direction of planetary migration from inwards to outwards, potentially saving planets in the inner, optically thick parts of the protoplanetary disc. When considering stellar irradiation in addition to viscous friction as a source of heating, the outer disc changes from a shadowed to a flared structure. Using a suited analytical formula it has been shown that in the flared part of the disc the migration is inwards; planets can migrate outwards only in shadowed regions of the disc, { because the radial gradient of entropy is stronger there}. In order to confirm this result numerically, we have computed the total torque acting on planets held on fixed orbits embedded in stellar irradiated 3D discs using the hydrodynamical code FARGOCA. We find qualitatively good agreement between the total torque obtained with numerical simulations and the one predicted by the analytical formula. For large masses (>20 Earth masses) we find quantitative agreement, and we obtain outwards migration regions for planets up to 60 Earth masses in the early stages of accretional discs. We find nevertheless that the agreement with the analytic formula is quite fortuitous because the formula underestimates the size of the horseshoe region; this error is compensated by imperfect estimates of other terms, most likely the cooling rate and saturation.

preprint2015arXiv

The great dichotomy of the Solar System: small terrestrial embryos and massive giant planet cores

The basic structure of the solar system is set by the presence of low-mass terrestrial planets in its inner part and giant planets in its outer part. This is the result of the formation of a system of multiple embryos with approximately the mass of Mars in the inner disk and of a few multi-Earth-mass cores in the outer disk, within the lifetime of the gaseous component of the protoplanetary disk. What was the origin of this dichotomy in the mass distribution of embryos/cores? We show in this paper that the classic processes of runaway and oligarchic growth from a disk of planetesimals cannot explain this dichotomy, even if the original surface density of solids increased at the snowline. Instead, the accretion of drifting pebbles by embryos and cores can explain the dichotomy, provided that some assumptions hold true. We propose that the mass-flow of pebbles is two-times lower and the characteristic size of the pebbles is approximately ten times smaller within the snowline than beyond the snowline (respectively at heliocentric distance $r<r_{ice}$ and $r>r_{ice}$, where $r_{ice}$ is the snowline heliocentric distance), due to ice sublimation and the splitting of icy pebbles into a collection of chondrule-size silicate grains. In this case, objects of original sub-lunar mass would grow at drastically different rates in the two regions of the disk. Within the snowline these bodies would reach approximately the mass of Mars while beyond the snowline they would grow to $\sim 20$ Earth masses. The results may change quantitatively with changes to the assumed parameters, but the establishment of a clear dichotomy in the mass distribution of protoplanets appears robust, provided that there is enough turbulence in the disk to prevent the sedimentation of the silicate grains into a very thin layer.

preprint2014arXiv

Accretion of Jupiter-mass Planets in the Limit of Vanishing Viscosity

In the core-accretion model the nominal runaway gas-accretion phase brings most planets to multiple Jupiter masses. However, known giant planets are predominantly Jupiter-mass bodies. Obtaining longer timescales for gas accretion may require using realistic equations of states, or accounting for the dynamics of the circumplanetary disk (CPD) in low-viscosity regime, or both. Here we explore the second way using global, three-dimensional isothermal hydrodynamical simulations with 8 levels of nested grids around the planet. In our simulations the vertical inflow from the circumstellar disk (CSD) to the CPD determines the shape of the CPD and its accretion rate. Even without prescribed viscosity Jupiter's mass-doubling time is $\sim 10^4$ years, assuming the planet at 5.2 AU and a Minimum Mass Solar Nebula. However, we show that this high accretion rate is due to resolution-dependent numerical viscosity. Furthermore, we consider the scenario of a layered CSD, viscous only in its surface layer, and an inviscid CPD. We identify two planet-accretion mechanisms that are independent of the viscosity in the CPD: (i) the polar inflow -- defined as a part of the vertical inflow with a centrifugal radius smaller than 2 Jupiter-radii and (ii) the torque exerted by the star on the CPD. In the limit of zero effective viscosity, these two mechanisms would produce an accretion rate 40 times smaller than in the simulation.

preprint2014arXiv

Meridional circulation of gas into gaps opened by giant planets in three-dimensional low-viscosity disks

We examine the gas circulation near a gap opened by a giant planet in a protoplanetary disk. We show with high resolution 3D simulations that the gas flows into the gap at high altitude over the mid-plane, at a rate dependent on viscosity. We explain this observation with a simple conceptual model. From this model we derive an estimate of the amount of gas flowing into a gap opened by a planet with Hill radius comparable to the scale-height of a layered disk (i. e. a disk with viscous upper layer and inviscid midplane). Our estimate agrees with modern MRI simulations(Gressel et al., 2013). We conclude that gap opening in a layered disk can not slow down significantly the runaway gas accretion of Saturn to Jupiter-mass planets.

preprint2014arXiv

Migration of Earth-size planets in 3D radiative discs

In this paper, we address the migration of small mass planets in 3D radiative disks. Indeed, migration of small planets is known to be too fast inwards in locally isothermal conditions. However, thermal effects could reverse its direction, potentially saving planets in the inner, optically thick parts of the protoplanetary disc. This effect has been seen for masses larger than 5 Earth masses, but the minimum mass for this to happen has never been probed numerically, although it is of crucial importance for planet formation scenarios. We have extended the hydro-dynamical code FARGO to 3D, with thermal diffusion. With this code, we perform simulations of embedded planets down to 2 Earth masses. For a set of discs parameters for which outward migration has been shown in the range of $[5, 35]$ Earth masses, we find that the transition to inward migration occurs for masses in the range $[3, 5]$ Earth masses. The transition appears to be due to an unexpected phenomenon: the formation of an asymmetric cold and dense finger of gas driven by circulation and libration streamlines. We recover this phenomenon in 2D simulations where we control the cooling effects of the gas through a simple modeling of the energy equation.

preprint2013arXiv

Constraining the cometary flux through the asteroid belt during the late heavy bombardment

In the Nice model, the late heavy bombardment (LHB) is related to an orbital instability of giant planets which causes a fast dynamical dispersion of a transneptunian cometary disk. We study effects produced by these hypothetical cometary projectiles on main-belt asteroids. In particular, we want to check whether the observed collisional families provide a lower or an upper limit for the cometary flux during the LHB. We present an updated list of observed asteroid families as identified in the space of synthetic proper elements by the hierarchical clustering method, colour data, albedo data and dynamical considerations and we estimate their physical parameters. We selected 12 families which may be related to the LHB according to their dynamical ages. We then used collisional models and N-body orbital simulations to gain insight into the long-term dynamical evolution of synthetic LHB families over 4 Gyr. We account for the mutual collisions, the physical disruptions of comets, the Yarkovsky/YORP drift, chaotic diffusion, or possible perturbations by the giant-planet migration. Assuming a "standard" size-frequency distribution of primordial comets, we predict the number of families with parent-body sizes D_PB >= 200 km which seems consistent with observations. However, more than 100 asteroid families with D_PB >= 100 km should be created at the same time which are not observed. This discrepancy can be nevertheless explained by the following processes: i) asteroid families are efficiently destroyed by comminution (via collisional cascade), ii) disruptions of comets below some critical perihelion distance (q <~ 1.5 AU) are common. Given the freedom in the cometary-disruption law, we cannot provide stringent limits on the cometary flux, but we can conclude that the observed distribution of asteroid families does not contradict with a cometary LHB.

preprint2013arXiv

DUst Around NEarby Stars. The survey observational results

Debris discs are a consequence of the planet formation process and constitute the fingerprints of planetesimal systems. Their solar system's counterparts are the asteroid and Edgeworth-Kuiper belts. The DUNES survey aims at detecting extra-solar analogues to the Edgeworth-Kuiper belt around solar-type stars, putting in this way the solar system into context. The survey allows us to address some questions related to the prevalence and properties of planetesimal systems. We used {\it Herschel}/PACS to observe a sample of nearby FGK stars. Data at 100 and 160 $μ$m were obtained, complemented in some cases with observations at 70 $μ$m, and at 250, 350 and 500 $μ$m using SPIRE. The observing strategy was to integrate as deep as possible at 100 $μ$m to detect the stellar photosphere. Debris discs have been detected at a fractional luminosity level down to several times that of the Edgeworth-Kuiper belt. The incidence rate of discs around the DUNES stars is increased from a rate of $\sim$ 12.1% $\pm$ 5% before \emph{Herschel} to $\sim$ 20.2% $\pm$ 2%. A significant fraction ($\sim$ 52%) of the discs are resolved, which represents an enormous step ahead from the previously known resolved discs. Some stars are associated with faint far-IR excesses attributed to a new class of cold discs. Although it cannot be excluded that these excesses are produced by coincidental alignment of background galaxies, statistical arguments suggest that at least some of them are true debris discs. Some discs display peculiar SEDs with spectral indexes in the 70-160$μ$m range steeper than the Rayleigh-Jeans one. An analysis of the debris disc parameters suggests that a decrease might exist of the mean black body radius from the F-type to the K-type stars. In addition, a weak trend is suggested for a correlation of disc sizes and an anticorrelation of disc temperatures with the stellar age.

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.

preprint2013arXiv

Origin of the peculiar eccentricity distribution of the inner cold Kuiper belt

Dawson and Murray-Clay (2012) pointed out that the inner part of the cold population in the Kuiper belt (that with semi major axis a<43.5 AU) has orbital eccentricities significantly smaller than the limit imposed by stability constraints. Here, we confirm their result by looking at the orbital distribution and stability properties in proper element space. We show that the observed distribution could have been produced by the slow sweeping of the 4/7 mean motion resonance with Neptune that accompanied the end of Neptune's migration process. The orbital distribution of the hot Kuiper belt is not significantly affected in this process, for the reasons discussed in the main text. Therefore, the peculiar eccentricity distribution of the inner cold population can not be unequivocally interpreted as evidence that the cold population formed in-situ and was only moderately excited in eccentricity; it can simply be the signature of Neptune's radial motion, starting from a moderately eccentric orbit. We discuss how this agrees with a scenario of giant planet evolution following a dynamical instability and, possibly, with the radial transport of the cold population.

preprint2012arXiv

A Method to Constrain the Size of the Protosolar Nebula

Observations indicate that the gaseous circumstellar disks around young stars vary significantly in size, ranging from tens to thousands of AU. Models of planet formation depend critically upon the properties of these primordial disks, yet in general it is impossible to connect an existing planetary system with an observed disk. We present a method by which we can constrain the size of our own protosolar nebula using the properties of the small body reservoirs in the solar system. In standard planet formation theory, after Jupiter and Saturn formed they scattered a significant number of remnant planetesimals into highly eccentric orbits. In this paper, we show that if there had been a massive, extended protoplanetary disk at that time, then the disk would have excited Kozai oscillations in some of the scattered objects, driving them into high-inclination (i > 50 deg), low-eccentricity orbits (q > 30 AU). The dissipation of the gaseous disk would strand a subset of objects in these high-inclination orbits; orbits that are stable on Gyr time scales. To date, surveys have not detected any Kuiper Belt Objects with orbits consistent with this dynamical mechanism. Using these non-detections by the Deep Ecliptic Survey (DES) and the Palomar Distant Solar System Survey we are able to rule out an extended gaseous protoplanetary disk (R_D > 80 AU) in our solar system at the time of Jupiter's formation. Future deep all sky surveys such as the Large Synoptic Survey Telescope (LSST) will all us to further constrain the size of the protoplanetary disk.

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.

preprint2012arXiv

Dynamical capture in the Pluto-Charon system

This paper explores the possibility that the progenitors of the small satellites of Pluto got captured in the Pluto-Charon system from the massive heliocentric planetesimal disk in which Pluto was originally embedded into. We find that, if the dynamical excitation of the disk is small, temporary capture in the Pluto-Charon system can occur with non-negligible probability, due to the dynamical perturbations exerted by the binary nature of the Pluto-Charon pair. However, the captured objects remain on very elliptic orbits and the typical capture time is only 100 years. In order to explain the origin of the small satellites of Pluto, we conjecture that some of these objects got disrupted during their Pluto-bound phase by a collision with a planetesimal of the disk. This could have generated a debris disk, which damped under internal collisional evolution, until turning itself into an accretional disk that could form small satellites on circular orbits, co-planar with Charon. Unfortunately, we find that objects large enough to carry a sufficient amount of mass to generate the small satellites of Pluto have collisional lifetimes orders of magnitude longer than the capture time. Thus, this scenario cannot explain the origin of the small satellites of Pluto, which remains elusive.

preprint2011arXiv

Did the Hilda collisional family form during the late heavy bombardment?

We model the long-term evolution of the Hilda collisional family located in the 3/2 mean-motion resonance with Jupiter. Its eccentricity distribution evolves mostly due to the Yarkovsky/YORP effect and assuming that: (i) impact disruption was isotropic, and (ii) albedo distribution of small asteroids is the same as for large ones, we can estimate the age of the Hilda family to be $4_{-1}^{+0}\,{\rm Gyr}$. We also calculate collisional activity in the J3/2 region. Our results indicate that current collisional rates are very low for a 200\,km parent body such that the number of expected events over Gyrs is much smaller than one. The large age and the low probability of the collisional disruption lead us to the conclusion that the Hilda family might have been created during the Late Heavy Bombardment when the collisions were much more frequent. The Hilda family may thus serve as a test of orbital behavior of planets during the LHB. We tested the influence of the giant-planet migration on the distribution of the family members. The scenarios that are consistent with the observed Hilda family are those with fast migration time scales $\simeq 0.3\,{\rm Myr}$ to $3\,{\rm Myr}$, because longer time scales produce a family that is depleted and too much spread in eccentricity. Moreover, there is an indication that Jupiter and Saturn were no longer in a compact configuration (with period ratio $P_{\rm S}/P_{\rm J} > 2.09$) at the time when the Hilda family was created.

preprint2011arXiv

The Cratering History of Asteroid (21) Lutetia

The European Space Agency's Rosetta spacecraft passed by the main belt asteroid (21) Lutetia the 10th July 2010. With its ~100km size, Lutetia is one of the largest asteroids ever imaged by a spacecraft. During the flyby, the on-board OSIRIS imaging system acquired spectacular images of Lutetia's northern hemisphere revealing a complex surface scarred by numerous impact craters, reaching the maximum dimension of about 55km. In this paper, we assess the cratering history of the asteroid. For this purpose, we apply current models describing the formation and evolution of main belt asteroids, that provide the rate and velocity distributions of impactors. These models, coupled with appropriate crater scaling laws, allow us to interpret the observed crater size-frequency distribution (SFD) and constrain the cratering history. Thanks to this approach, we derive the crater retention age of several regions on Lutetia, namely the time lapsed since their formation or global surface reset. We also investigate the influence of various factors -like Lutetia's bulk structure and crater obliteration- on the observed crater SFDs and the estimated surface ages. From our analysis, it emerges that Lutetia underwent a complex collisional evolution, involving major local resurfacing events till recent times. The difference in crater density between the youngest and oldest recognized units implies a difference in age of more than a factor of 10. The youngest unit (Beatica) has an estimated age of tens to hundreds of Myr, while the oldest one (Achaia) formed during a period when the bombardment of asteroids was more intense than the current one, presumably around 3.6Gyr ago or older.

preprint2010arXiv

Cold dust around nearby stars (DUNES). First results: A resolved exo-Kuiper belt around the solar-like star zeta^2 Ret

We present the first far-IR observations of the solar-type stars delta Pav, HR 8501, 51 Peg and zeta^2 Ret, taken within the context of the DUNES Herschel Open Time Key Programme (OTKP). This project uses the PACS and SPIRE instruments with the objective of studying infrared excesses due to exo-Kuiper belts around nearby solar-type stars. The observed 100 um fluxes from delta Pav, HR 8501, and 51 Peg agree with the predicted photospheric fluxes, excluding debris disks brighter than Ldust/Lstar ~ 5 x 10^-7 (1 sigma level) around those stars. A flattened, disk-like structure with a semi-major axis of ~ 100 AU in size is detected around zeta^2 Ret. The resolved structure suggests the presence of an eccentric dust ring, which we interpret as an exo-Kuiper belt with Ldust/Lstar ~ 10^-5.

preprint2009arXiv

Heating of near-Earth objects and meteoroids due to close approaches to the Sun

It is known that near-Earth objects (NEOs) during their orbital evolution may often undergo close approaches to the Sun. Indeed it is estimated that up to ~70% of them end their orbital evolution colliding with the Sun. Starting from the present orbital properties, it is possible to compute the most likely past evolution for every NEO, and to trace its distance from the Sun. We find that a large fraction of the population may have experienced in the past frequent close approaches, and thus, as a consequence, a considerable Sun-driven heating, not trivially correlated to the present orbits. The detailed dynamical behaviour, the rotational and the thermal properties of NEOs determine the exact amount of the resulting heating due to the Sun. In the present paper we discuss the general features of the process, providing estimates of the surface temperature reached by NEOs during their evolution. Moreover, we investigate the effects of this process on meteor-size bodies, analyzing possible differences with the NEO population. We also discuss some possible effects of the heating which can be observed through remote sensing by ground-based surveys or space missions.

preprint2009arXiv

Long range outward migration of giant planets, with application to Fomalhaut b

Recent observations of exoplanets by direct imaging, reveal that giant planets orbit at a few dozens to more than a hundred of AU from their central star. The question of the origin of these planets challenges the standard theories of planet formation. We propose a new way of obtaining such far planets, by outward migration of a pair of planets formed in the 10 AU region. Two giant planets in mean motion resonance in a common gap in the protoplanetary disk migrate outwards, if the inner one is significantly more massive than the outer one. Using hydrodynamical simulations, we show that their semi major axes can increase by almost one order of magnitude. In a flared disk, the pair of planets should reach an asymptotic radius. This mechanism could account for the presence of Fomalhaut b ; then, a second, more massive planet, should be orbiting Fomalhaut at about 75 AU.