-
Formation of super-Mercuries via giant impacts
Authors:
Jingyao Dou,
Philip J. Carter,
Zoë M. Leinhardt
Abstract:
During the final stage of planetary formation, different formation pathways of planetary embryos could significantly influence the observed variations in planetary densities. Of the approximately 5,000 exoplanets identified to date, a notable subset exhibit core fractions reminiscent of Mercury, potentially a consequence of high-velocity giant impacts. In order to better understand the influence o…
▽ More
During the final stage of planetary formation, different formation pathways of planetary embryos could significantly influence the observed variations in planetary densities. Of the approximately 5,000 exoplanets identified to date, a notable subset exhibit core fractions reminiscent of Mercury, potentially a consequence of high-velocity giant impacts. In order to better understand the influence of such collisions on planetary formation and compositional evolution, we conducted an extensive set of smoothed particle hydrodynamics giant impact simulations between two-layered rocky bodies. These simulations spanned a broad range of impact velocities from one to eleven times the mutual escape velocity. We derived novel scaling laws that estimate the mass and core mass fraction of the largest post-collision remnants. Our findings indicate that the extent of core vaporization markedly influences mantle stripping efficiency at low impact angles. We delineate the distinct roles played by two mechanisms -- kinetic momentum transfer and vaporization-induced ejection -- in mantle stripping. Our research suggests that collisional outcomes for multi-layered planets are more complex than those for undifferentiated planetesimal impacts. Thus, a single universal law may not encompass all collision processes. We found a significant decrease in the mantle stripping efficiency as the impact angle increases. To form a 5 M$_{\oplus}$ super-Mercury at $45^{\circ}$, an impact velocity over 200 km s$^{-1}$ is required. This poses a challenge to the formation of super-Mercuries through a single giant impact, implying that their formation would either favor relatively low-angle single impacts or multiple
△ Less
Submitted 6 March, 2024;
originally announced March 2024.
-
Aerosols are not Spherical Cows: Using Discrete Dipole Approximation to Model the Properties of Fractal Particles
Authors:
Matt G. Lodge,
Hannah R. Wakeford,
Zoe M. Leinhardt
Abstract:
The optical properties of particulate-matter aerosols, within the context of exoplanet and brown dwarf atmospheres, are compared using three different models: Mie theory, Modified Mean Field (MMF) Theory, and Discrete Dipole Approximation (DDA). Previous results have demonstrated that fractal haze particles (MMF and DDA) absorb much less long-wavelength radiation than their spherical counterparts…
▽ More
The optical properties of particulate-matter aerosols, within the context of exoplanet and brown dwarf atmospheres, are compared using three different models: Mie theory, Modified Mean Field (MMF) Theory, and Discrete Dipole Approximation (DDA). Previous results have demonstrated that fractal haze particles (MMF and DDA) absorb much less long-wavelength radiation than their spherical counterparts (Mie), however it is shown here that the opposite can also be true if a more varying refractive index profile is used. Additionally, it is demonstrated that absorption and scattering cross-sections, as well as the asymmetry parameter, are underestimated if Mie theory is used. Although DDA can be used to obtain more accurate results, it is known to be much more computationally intensive; to avoid this, the use of low-resolution aerosol models is explored, which could dramatically speed up the process of obtaining accurate computations of optical cross-sections within a certain parameter space. The validity of DDA is probed for wavelengths of interest for observations of aerosols within exoplanet and brown dwarf atmospheres (0.2 to 15 micrometres). Finally, novel code is presented to compare the results of Mie, MMF and DDA theories (CORAL: Comparison Of Radiative AnaLyses), as well as to increase and decrease the resolution of DDA shape files accordingly (SPHERIFY). Both codes can be applied to a range of other interesting astrophysical environments in addition to exoplanet atmospheres, for example dust grains within protoplanetary disks.
△ Less
Submitted 4 December, 2023;
originally announced December 2023.
-
Post-Giant Impact Planetesimals Sustaining Extreme Debris Disks
Authors:
Lewis Watt,
Zoë M. Leinhardt,
Philip J. Carter
Abstract:
Extreme debris disks can show short term behaviour through the evolution and clearing of small grains produced in giant impacts, and potentially a longer period of variability caused by a planetesimal population formed from giant impact ejecta. In this paper, we present results of numerical simulations to explain how a planetesimal populated disk can supply an observed extreme debris disk with sma…
▽ More
Extreme debris disks can show short term behaviour through the evolution and clearing of small grains produced in giant impacts, and potentially a longer period of variability caused by a planetesimal population formed from giant impact ejecta. In this paper, we present results of numerical simulations to explain how a planetesimal populated disk can supply an observed extreme debris disk with small grains. We simulated a sample of giant impacts from which we form a planetesimal population. We then use the $N$-body code {\sc Rebound} to evolve the planetesimals spatially and collisionally. We adopt a simplistic collision criteria in which we define destructive collisions to be between planetesimals with a mutual impact velocity that exceeds two times the catastrophic disruption threshold, $V^*$. We find that for some configurations, a planetesimal populated disk can produce a substantial amount of dust to sustain an observable disk. The semi-major axis at which the giant impact occurs changes the mass added to the observed disk substantially while the orientation of the impact has less of an effect. We determine how the collision rate at the collision point changes over time and show that changes in semi-major axis and orientation only change the initial collision rate of the disk. Collision rates across all disks evolve at a similar rate.
△ Less
Submitted 19 November, 2023;
originally announced November 2023.
-
A super-massive Neptune-sized planet
Authors:
L. Naponiello,
L. Mancini,
A. Sozzetti,
A. S. Bonomo,
A. Morbidelli,
J. Dou,
L. Zeng,
Z. M. Leinhardt,
K. Biazzo,
P. Cubillos,
M. Pinamonti,
D. Locci,
A. Maggio,
M. Damasso,
A. F. Lanza,
J. J. Lissauer,
A. Bignamini,
W. Boschin,
L. G. Bouma,
P. J. Carter,
D. R. Ciardi,
K. A. Collins,
R. Cosentino,
I. Crossfield,
S. Desidera
, et al. (33 additional authors not shown)
Abstract:
Neptune-sized planets exhibit a wide range of compositions and densities, depending onf cators related to their formation and evolution history, such as the distance from their host stars and atmospheric escape processes. They can vary from relatively low-density planets with thick hydrogen-helium atmospheres to higher-density planets with a substantial amount of water or a rocky interior with a t…
▽ More
Neptune-sized planets exhibit a wide range of compositions and densities, depending onf cators related to their formation and evolution history, such as the distance from their host stars and atmospheric escape processes. They can vary from relatively low-density planets with thick hydrogen-helium atmospheres to higher-density planets with a substantial amount of water or a rocky interior with a thinner atmosphere, such as HD 95338 b, TOI-849 b and TOI-2196 b. The discovery of exoplanets in the hot-Neptune desert, a region close to the host stars with a deficit of Neptune-sized planets, provides insights into the formation and evolution of planetary systems, including the existence of this region itself. Here we show observations of the transiting planet TOI-1853 b, which has a radius of 3.46 +- 0.08 Earth radii and orbits a dwarf star every 1.24 days. This planet has a mass of 73.2 +- 2.7 Earth masses, almost twice that of any other Neptune-sized planet known so far, and a density of 9.7 +- 0.8 grams per cubic centimetre. These values place TOI-1853 b in the middle of the Neptunian desert and imply that heavy elements dominate its mass. The properties of TOI-1853 b present a puzzle for conventional theories of planetary formation and evolution, and could be the result of several proto-planet collisions or the final state of an initially high-eccentricity planet that migrated closer to its parent star.
△ Less
Submitted 4 September, 2023;
originally announced September 2023.
-
Isolating the extreme debris disc signature -- explorations of eccentric extreme debris discs formed by giant impacts
Authors:
Thomas Lewis,
Lewis Watt,
Zoë M Leinhardt
Abstract:
In this work, we used N-body simulations and a radiative transfer package to model the evolution of eccentric debris discs produced by giant impacts between planetary embryos. This included how the morphology and infrared emission of these discs varied with embryo eccentricity and collision true anomaly. We found that eccentric discs inherit the eccentric properties of the centre of mass orbit of…
▽ More
In this work, we used N-body simulations and a radiative transfer package to model the evolution of eccentric debris discs produced by giant impacts between planetary embryos. This included how the morphology and infrared emission of these discs varied with embryo eccentricity and collision true anomaly. We found that eccentric discs inherit the eccentric properties of the centre of mass orbit of the two colliding embryos. However, the orientation of the collision with the respect to this orbit plays a key role in determining how closely the disc material resembles the centre of mass orbit. Additionally, we found that increased eccentricity acted to suppress the formation of certain short-term variations in the disc emission depending on the collision position. These short-term variations have been associated with an observational phenomenon called extreme debris discs. Short-term variability has been suggested as a potential signature for giant impacts.
△ Less
Submitted 9 January, 2023;
originally announced January 2023.
-
Atmosphere Loss in Oblique Super-Earth Impacts
Authors:
Thomas R. Denman,
Zoe M. Leinhardt,
Philip J. Carter
Abstract:
Using smoothed particle hydrodynamics we model giant impacts of Super-Earth mass rocky planets between an atmosphere-less projectile and an atmosphere-rich target. In this work we present results from head-on to grazing collisions. The results of the simulations fall into two broad categories: 1) one main post-collision remnant containing material from target and projectile; 2) two main post-colli…
▽ More
Using smoothed particle hydrodynamics we model giant impacts of Super-Earth mass rocky planets between an atmosphere-less projectile and an atmosphere-rich target. In this work we present results from head-on to grazing collisions. The results of the simulations fall into two broad categories: 1) one main post-collision remnant containing material from target and projectile; 2) two main post-collision remnants resulting from `erosive hit-and-run' collisions. All collisions removed at least some of the target atmosphere, in contrast to the idealised hit-and-run definition in which the target mass is unchanged. We find that the boundary between `hit-and-run' collisions and collisions that result in the projectile and target accreting/merging to be strongly correlated with the mutual escape velocity at the predicted point of closest approach. Our work shows that it is very unlikely for a single giant impact to remove all of the atmosphere. For all the atmosphere to be removed, head-on impacts require roughly the energy of catastrophic disruption (i.e. permanent ejection of half the total system mass) and result in significant erosion of the mantle. We show that higher impact angle collisions, which are more common, are less efficient at atmosphere removal than head-on collisions. Therefore, single collisions that remove all the atmosphere without substantially disrupting the planet are not expected during planet formation.
△ Less
Submitted 31 March, 2022;
originally announced March 2022.
-
Atmosphere loss in planet-planet collisions
Authors:
Thomas R. Denman,
Zoe M. Leinhardt,
Phil J. Carter,
Christoph Mordasini
Abstract:
Many of the planets discovered by the Kepler satellite are close orbiting Super-Earths or Mini-Neptunes. Such objects exhibit a wide spread of densities for similar masses. One possible explanation for this density spread is giant collisions stripping planets of their atmospheres. In this paper we present the results from a series of smoothed particle hydrodynamics (SPH) simulations of head-on col…
▽ More
Many of the planets discovered by the Kepler satellite are close orbiting Super-Earths or Mini-Neptunes. Such objects exhibit a wide spread of densities for similar masses. One possible explanation for this density spread is giant collisions stripping planets of their atmospheres. In this paper we present the results from a series of smoothed particle hydrodynamics (SPH) simulations of head-on collisions of planets with significant atmospheres and bare projectiles without atmospheres. Collisions between planets can have sufficient energy to remove substantial fractions of the mass from the target planet. We find the fraction of mass lost splits into two regimes -- at low impact energies only the outer layers are ejected corresponding to atmosphere dominated loss, at higher energies material deeper in the potential is excavated resulting in significant core and mantle loss. Mass removal is less efficient in the atmosphere loss dominated regime compared to the core and mantle loss regime, due to the higher compressibility of atmosphere relative to core and mantle. We find roughly twenty per cent atmosphere remains at the transition between the two regimes. We find that the specific energy of this transition scales linearly with the ratio of projectile to target mass for all projectile-target mass ratios measured. The fraction of atmosphere lost is well approximated by a quadratic in terms of the ratio of specific energy and transition energy. We provide algorithms for the incorporation of our scaling law into future numerical studies.
△ Less
Submitted 2 June, 2020;
originally announced June 2020.
-
Extreme Debris Disk Variability -- Exploring the Diverse Outcomes of Large Asteroid Impacts During the Era of Terrestrial Planet Formation
Authors:
Kate Y. L. Su,
Alan P. Jackson,
Andras Gaspar,
George H. Rieke,
Ruobing Dong,
Johan Olofsson,
G. M. Kennedy,
Zoe M. Leinhardt,
Renu Malhotra,
Michael Hammer,
Huan Y. A. Meng,
W. Rujopakarn,
Joseph E. Rodriguez,
Joshua Pepper,
D. E. Reichart,
David James,
Keivan G. Stassun
Abstract:
The most dramatic phases of terrestrial planet formation are thought to be oligarchic and chaotic growth, on timescales of up to 100-200 Myr, when violent impacts occur between large planetesimals of sizes up to proto-planets. Such events are marked by the production of large amounts of debris as has been observed in some exceptionally bright and young debris disks (termed extreme debris disks). H…
▽ More
The most dramatic phases of terrestrial planet formation are thought to be oligarchic and chaotic growth, on timescales of up to 100-200 Myr, when violent impacts occur between large planetesimals of sizes up to proto-planets. Such events are marked by the production of large amounts of debris as has been observed in some exceptionally bright and young debris disks (termed extreme debris disks). Here we report five years of Spitzer measurements of such systems around two young solar-type stars: ID8 and P1121. The short-term (weekly to monthly) and long-term (yearly) disk variability is consistent with the aftermaths of large impacts involving large asteroid-size bodies. We demonstrate that an impact-produced clump of optically thick dust, under the influence of the dynamical and viewing geometry effects, can produce short-term modulation in the disk light curves. The long-term disk flux variation is related to the collisional evolution within the impact-produced fragments once released into a circumstellar orbit. The time-variable behavior observed in the P1121 system is consistent with a hypervelocity impact prior to 2012 that produced vapor condensates as the dominant impact product. Two distinct short-term modulations in the ID8 system argue for two violent impacts at different times and locations. Its long-term variation is consistent with the collisional evolution of two different populations of impact-produced debris dominated by either vapor condensates or escaping boulders. The bright, variable emission from the dust produced in large impacts from extreme debris disks provides a unique opportunity to study violent events during the era of terrestrial planet formation.
△ Less
Submitted 19 April, 2019; v1 submitted 25 March, 2019;
originally announced March 2019.
-
A giant impact as the likely origin of different twins in the Kepler-107 exoplanet system
Authors:
Aldo S. Bonomo,
Li Zeng,
Mario Damasso,
Zoë M. Leinhardt,
Anders B. Justesen,
Eric Lopez,
Mikkel N. Lund,
Luca Malavolta,
Victor Silva Aguirre,
Lars A. Buchhave,
Enrico Corsaro,
Thomas Denman,
Mercedes Lopez-Morales,
Sean M. Mills,
Annelies Mortier,
Ken Rice,
Alessandro Sozzetti,
Andrew Vanderburg,
Laura Affer,
Torben Arentoft,
Mansour Benbakoura,
François Bouchy,
Jørgen Christensen-Dalsgaard,
Andrew Collier Cameron,
Rosario Cosentino
, et al. (27 additional authors not shown)
Abstract:
Measures of exoplanet bulk densities indicate that small exoplanets with radius less than 3 Earth radii ($R_\oplus$) range from low-density sub-Neptunes containing volatile elements to higher density rocky planets with Earth-like or iron-rich (Mercury-like) compositions. Such astonishing diversity in observed small exoplanet compositions may be the product of different initial conditions of the pl…
▽ More
Measures of exoplanet bulk densities indicate that small exoplanets with radius less than 3 Earth radii ($R_\oplus$) range from low-density sub-Neptunes containing volatile elements to higher density rocky planets with Earth-like or iron-rich (Mercury-like) compositions. Such astonishing diversity in observed small exoplanet compositions may be the product of different initial conditions of the planet-formation process and/or different evolutionary paths that altered the planetary properties after formation. Planet evolution may be especially affected by either photoevaporative mass loss induced by high stellar X-ray and extreme ultraviolet (XUV) flux or giant impacts. Although there is some evidence for the former, there are no unambiguous findings so far about the occurrence of giant impacts in an exoplanet system. Here, we characterize the two innermost planets of the compact and near-resonant system Kepler-107. We show that they have nearly identical radii (about $1.5-1.6~R_\oplus$), but the outer planet Kepler-107c is more than twice as dense (about $12.6~\rm g\,cm^{-3}$) as the innermost Kepler-107b (about $5.3~\rm g\,cm^{-3}$). In consequence, Kepler-107c must have a larger iron core fraction than Kepler-107b. This imbalance cannot be explained by the stellar XUV irradiation, which would conversely make the more-irradiated and less-massive planet Kepler-107b denser than Kepler-107c. Instead, the dissimilar densities are consistent with a giant impact event on Kepler-107c that would have stripped off part of its silicate mantle. This hypothesis is supported by theoretical predictions from collisional mantle stripping, which match the mass and radius of Kepler-107c.
△ Less
Submitted 4 February, 2019;
originally announced February 2019.
-
The origin of the Moon within a terrestrial synestia
Authors:
Simon J. Lock,
Sarah T. Stewart,
Michail I. Petaev,
Zoe M. Leinhardt,
Mia T. Mace,
Stein B. Jacobsen,
Matija Ćuk
Abstract:
The giant impact hypothesis remains the leading theory for lunar origin. However, current models struggle to explain the Moon's composition and isotopic similarity with Earth. Here we present a new lunar origin model. High-energy, high-angular momentum giant impacts can create a post-impact structure that exceeds the corotation limit (CoRoL), which defines the hottest thermal state and angular mom…
▽ More
The giant impact hypothesis remains the leading theory for lunar origin. However, current models struggle to explain the Moon's composition and isotopic similarity with Earth. Here we present a new lunar origin model. High-energy, high-angular momentum giant impacts can create a post-impact structure that exceeds the corotation limit (CoRoL), which defines the hottest thermal state and angular momentum possible for a corotating body. In a typical super-CoRoL body, traditional definitions of mantle, atmosphere and disk are not appropriate, and the body forms a new type of planetary structure, named a synestia. Using simulations of cooling synestias combined with dynamic, thermodynamic and geochemical calculations, we show that satellite formation from a synestia can produce the main features of our Moon. We find that cooling drives mixing of the structure, and condensation generates moonlets that orbit within the synestia, surrounded by tens of bars of bulk silicate Earth (BSE) vapor. The moonlets and growing moon are heated by the vapor until the first major element (Si) begins to vaporize and buffer the temperature. Moonlets equilibrate with BSE vapor at the temperature of silicate vaporization and the pressure of the structure, establishing the lunar isotopic composition and pattern of moderately volatile elements. Eventually, the cooling synestia recedes within the lunar orbit, terminating the main stage of lunar accretion. Our model shifts the paradigm for lunar origin from specifying a certain impact scenario to achieving a Moon-forming synestia. Giant impacts that produce potential Moon-forming synestias were common at the end of terrestrial planet formation.
△ Less
Submitted 27 February, 2018;
originally announced February 2018.
-
Collisional stripping of planetary crusts
Authors:
Philip J. Carter,
Zoë M. Leinhardt,
Tim Elliott,
Sarah T. Stewart,
Michael J. Walter
Abstract:
Geochemical studies of planetary accretion and evolution have invoked various degrees of collisional erosion to explain differences in bulk composition between planets and chondrites. Here we undertake a full, dynamical evaluation of 'crustal stripping' during accretion and its key geochemical consequences. We present smoothed particle hydrodynamics simulations of collisions between differentiated…
▽ More
Geochemical studies of planetary accretion and evolution have invoked various degrees of collisional erosion to explain differences in bulk composition between planets and chondrites. Here we undertake a full, dynamical evaluation of 'crustal stripping' during accretion and its key geochemical consequences. We present smoothed particle hydrodynamics simulations of collisions between differentiated rocky planetesimals and planetary embryos. We find that the crust is preferentially lost relative to the mantle during impacts, and we have developed a scaling law that approximates the mass of crust that remains in the largest remnant. Using this scaling law and a recent set of N-body simulations, we have estimated the maximum effect of crustal stripping on incompatible element abundances during the accretion of planetary embryos. We find that on average one third of the initial crust is stripped from embryos as they accrete, which leads to a reduction of ~20% in the budgets of the heat producing elements if the stripped crust does not reaccrete. Erosion of crusts can lead to non-chondritic ratios of incompatible elements, but the magnitude of this effect depends sensitively on the details of the crust-forming melting process. The Lu/Hf system is fractionated for a wide range of crustal formation scenarios. Using eucrites (the products of planetesimal silicate melting, thought to represent the crust of Vesta) as a guide to the Lu/Hf of planetesimal crust partially lost during accretion, we predict the Earth could evolve to a superchondritic 176-Hf/177-Hf (3-5 parts per ten thousand) at present day. Such values are in keeping with compositional estimates of the bulk Earth. Stripping of planetary crusts during accretion can lead to detectable changes in bulk composition of lithophile elements, but the fractionation is relatively subtle, and sensitive to the efficiency of reaccretion.
△ Less
Submitted 7 December, 2017;
originally announced December 2017.
-
Explaining the variability of WD 1145+017 with simulations of asteroid tidal disruption
Authors:
Dimitri Veras,
Philip J. Carter,
Zoe M. Leinhardt,
Boris T. Gaensicke
Abstract:
Post-main-sequence planetary science has been galvanised by the striking variability, depth and shape of the photometric transit curves due to objects orbiting white dwarf WD 1145+017, a star which also hosts a dusty debris disc and circumstellar gas, and displays strong metal atmospheric pollution. However, the physical properties of the likely asteroid which is discharging disintegrating fragmen…
▽ More
Post-main-sequence planetary science has been galvanised by the striking variability, depth and shape of the photometric transit curves due to objects orbiting white dwarf WD 1145+017, a star which also hosts a dusty debris disc and circumstellar gas, and displays strong metal atmospheric pollution. However, the physical properties of the likely asteroid which is discharging disintegrating fragments remain largely unconstrained from the observations. This process has not yet been modelled numerically. Here, we use the N-body code PKDGRAV to compute dissipation properties for asteroids of different spins, densities, masses, and eccentricities. We simulate both homogeneous and differentiated asteroids, for up to two years, and find that the disruption timescale is strongly dependent on density and eccentricity, but weakly dependent on mass and spin. We find that primarily rocky differentiated bodies with moderate (~3-4 g/cm^3) bulk densities on near-circular (e <~ 0.1) orbits can remain intact while occasionally shedding mass from their mantles. These results suggest that the asteroid orbiting WD 1145+017 is differentiated, resides just outside of the Roche radius for bulk density but just inside the Roche radius for mantle density, and is more akin physically to an asteroid like Vesta instead of one like Itokawa.
△ Less
Submitted 21 October, 2016;
originally announced October 2016.
-
Modelling circumbinary protoplanetary disks II. Gas disk feedback on planetesimal dynamical and collisional evolution in the circumbinary systems Kepler-16 and 34
Authors:
S. Lines,
Z. M. Leinhardt,
C. Baruteau,
S. -J. Paardekooper,
P. J. Carter
Abstract:
Aims. We investigate the feasibility of planetesimal growth in circumbinary protoplanetary disks around the observed systems Kepler- 16 and Kepler-34 under the gravitational influence of a precessing eccentric gas disk. Methods. We embed the results of our previous hydrodynamical simulations of protoplanetary disks around binaries into an N-body code to perform 3D, high-resolution, inter-particle…
▽ More
Aims. We investigate the feasibility of planetesimal growth in circumbinary protoplanetary disks around the observed systems Kepler- 16 and Kepler-34 under the gravitational influence of a precessing eccentric gas disk. Methods. We embed the results of our previous hydrodynamical simulations of protoplanetary disks around binaries into an N-body code to perform 3D, high-resolution, inter-particle gravity-enabled simulations of planetesimal growth and dynamics that include the gravitational force imparted by the gas. Results. Including the full, precessing asymmetric gas disk generates high eccentricity orbits for planetesimals orbiting at the edge of the circumbinary cavity, where the gas surface density and eccentricity have their largest values. The gas disk is able to efficiently align planetesimal pericenters in some regions leading to phased, non-interacting orbits. Outside of these areas eccentric planetesimal orbits become misaligned and overlap leading to crossing orbits and high relative velocities during planetesimal collisions. This can lead to an increase in the number of erosive collisions that far outweighs the number of collisions that result in growth. Gravitational focusing from the static axisymmetric gas disk is weak and does not significantly alter collision outcomes from the gas free case. Conclusions. Due to asymmetries in the gas disk, planetesimals are strongly perturbed onto highly eccentric orbits. Where planetesimals orbits are not well aligned, orbit crossings lead to an increase in the number of erosive collisions. This makes it difficult for sustained planetesimal accretion to occur at the location of Kepler-16b and Kepler-34b and we therefore rule out in-situ growth. This adds further support to our initial suggestions that most circumbinary planets should form further out in the disk and migrate inwards.
△ Less
Submitted 6 April, 2016;
originally announced April 2016.
-
Hiding in the Shadows II: Collisional Dust as Exoplanet Markers
Authors:
Jack Dobinson,
Zoe M. Leinhardt,
Stefan Lines,
Philip J. Carter,
Sarah E. Dodson-Robinson,
Nick A. Teanby
Abstract:
Observations of the youngest planets ($\sim$1-10 Myr for a transitional disk) will increase the accuracy of our planet formation models. Unfortunately, observations of such planets are challenging and time-consuming to undertake even in ideal circumstances. Therefore, we propose the determination of a set of markers that can pre-select promising exoplanet-hosting candidate disks. To this end, N-bo…
▽ More
Observations of the youngest planets ($\sim$1-10 Myr for a transitional disk) will increase the accuracy of our planet formation models. Unfortunately, observations of such planets are challenging and time-consuming to undertake even in ideal circumstances. Therefore, we propose the determination of a set of markers that can pre-select promising exoplanet-hosting candidate disks. To this end, N-body simulations were conducted to investigate the effect of an embedded Jupiter mass planet on the dynamics of the surrounding planetesimal disk and the resulting creation of second generation collisional dust. We use a new collision model that allows fragmentation and erosion of planetesimals, and dust-sized fragments are simulated in a post process step including non-gravitational forces due to stellar radiation and a gaseous protoplanetary disk. Synthetic images from our numerical simulations show a bright double ring at 850 $μ$m for a low eccentricity planet, whereas a high eccentricity planet would produce a characteristic inner ring with asymmetries in the disk. In the presence of first generation primordial dust these markers would be difficult to detect far from the orbit of the embedded planet, but would be detectable inside a gap of planetary origin in a transitional disk.
△ Less
Submitted 14 January, 2016; v1 submitted 13 January, 2016;
originally announced January 2016.
-
Erosive Hit-and-Run Impact Events: Debris Unbound
Authors:
Gal Sarid,
Sarah T. Stewart,
Zoe M. Leinhardt
Abstract:
Erosive collisions among planetary embryos in the inner solar system can lead to multiple remnant bodies, varied in mass, composition and residual velocity. Some of the smaller, unbound debris may become available to seed the main asteroid belt. The makeup of these collisionally produced bodies is different from the canonical chondritic composition, in terms of rock/iron ratio and may contain furt…
▽ More
Erosive collisions among planetary embryos in the inner solar system can lead to multiple remnant bodies, varied in mass, composition and residual velocity. Some of the smaller, unbound debris may become available to seed the main asteroid belt. The makeup of these collisionally produced bodies is different from the canonical chondritic composition, in terms of rock/iron ratio and may contain further shock-processed material. Having some of the material in the asteroid belt owe its origin from collisions of larger planetary bodies may help in explaining some of the diversity and oddities in composition of different asteroid groups.
△ Less
Submitted 31 October, 2016; v1 submitted 16 November, 2015;
originally announced November 2015.
-
Compositional evolution during rocky protoplanet accretion
Authors:
Philip J. Carter,
Zoë M. Leinhardt,
Tim Elliott,
Michael J. Walter,
Sarah T. Stewart
Abstract:
The Earth appears non-chondritic in its abundances of refractory lithophile elements, posing a significant problem for our understanding of its formation and evolution. It has been suggested that this non-chondritic composition may be explained by collisional erosion of differentiated planetesimals of originally chondritic composition. In this work, we present N-body simulations of terrestrial pla…
▽ More
The Earth appears non-chondritic in its abundances of refractory lithophile elements, posing a significant problem for our understanding of its formation and evolution. It has been suggested that this non-chondritic composition may be explained by collisional erosion of differentiated planetesimals of originally chondritic composition. In this work, we present N-body simulations of terrestrial planet formation that track the growth of planetary embryos from planetesimals. We simulate evolution through the runaway and oligarchic growth phases under the Grand Tack model and in the absence of giant planets. These simulations include a state-of-the-art collision model which allows multiple collision outcomes, such as accretion, erosion, and bouncing events, that enables tracking of the evolving core mass fraction of accreting planetesimals. We show that the embryos grown during this intermediate stage of planet formation exhibit a range of core mass fractions, and that with significant dynamical excitation, enough mantle can be stripped from growing embryos to account for the Earth's non-chondritic Fe/Mg ratio. We also find that there is a large diversity in the composition of remnant planetesimals, with both iron-rich and silicate-rich fragments produced via collisions.
△ Less
Submitted 24 September, 2015;
originally announced September 2015.
-
Modelling circumbinary protoplanetary disks: I. Fluid simulations of the Kepler-16 and 34 systems
Authors:
S. Lines,
Z. M. Leinhardt,
C. Baruteau,
S. -J. Paardekooper,
P. J. Carter
Abstract:
The Kepler mission's discovery of a number of circumbinary planets orbiting close (a_p < 1.1 au) to the stellar binary raises questions as to how these planets could have formed given the intense gravitational perturbations the dual stars impart on the disk. The gas component of circumbinary protoplanetary disks is perturbed in a similar manner to the solid, planetesimal dominated counterpart, alt…
▽ More
The Kepler mission's discovery of a number of circumbinary planets orbiting close (a_p < 1.1 au) to the stellar binary raises questions as to how these planets could have formed given the intense gravitational perturbations the dual stars impart on the disk. The gas component of circumbinary protoplanetary disks is perturbed in a similar manner to the solid, planetesimal dominated counterpart, although the mechanism by which disk eccentricity originates differs. This is the first work of a series that aims to investigate the conditions for planet formation in circumbinary protoplanetary disks. We present a number of hydrodynamical simulations that explore the response of gas disks around two observed binary systems: Kepler-16 and Kepler-34. We probe the importance of disk viscosity, aspect-ratio, inner boundary condition, initial surface density gradient, and self-gravity on the dynamical evolution of the disk, as well as its quasi steady-state profile. We find there is a strong influence of binary type on the mean disk eccentricity, e_d, leading to e_d = 0.02 - 0.08 for Kepler-16 and e_d = 0.10 - 0.15 in Kepler-34. The value of alpha-viscosity has little influence on the disk, but we find a strong increase in mean disk eccentricity with increasing aspect-ratio due to wave propagation effects. The choice of inner boundary condition only has a small effect on the surface density and eccentricity of the disk. Our primary finding is that including disk self-gravity has little impact on the evolution or final state of the disk for disks with masses less than 12.5 times that of the minimum-mass solar nebula. This finding contrasts with the results of self-gravity relevance in circumprimary disks, where its inclusion is found to be an important factor in describing the disk evolution.
△ Less
Submitted 24 July, 2015;
originally announced July 2015.
-
Numerically Predicted Indirect Signatures of Terrestrial Planet Formation
Authors:
Zoë M. Leinhardt,
Jack Dobinson,
Philip J. Carter,
Stefan Lines
Abstract:
The intermediate phases of planet formation are not directly observable due to lack of emission from planetesimals. Planet formation is, however, a dynamically active process resulting in collisions between the evolving planetesimals and the production of dust. Thus, indirect observation of planet formation may indeed be possible in the near future. In this paper we present synthetic observations…
▽ More
The intermediate phases of planet formation are not directly observable due to lack of emission from planetesimals. Planet formation is, however, a dynamically active process resulting in collisions between the evolving planetesimals and the production of dust. Thus, indirect observation of planet formation may indeed be possible in the near future. In this paper we present synthetic observations based on numerical N-body simulations of the intermediate phase of planet formation including a state-of-the-art collision model, EDACM, which allows multiple collision outcomes, such as, accretion, erosion, and bouncing events. We show that the formation of planetary embryos may be indirectly observable by a fully functioning ALMA telescope if the surface area involved in planetesimal evolution is sufficiently large and/or the amount of dust produced in the collisions is sufficiently high in mass.
△ Less
Submitted 17 June, 2015;
originally announced June 2015.
-
Formation of planetary debris discs around white dwarfs II: Shrinking extremely eccentric collisionless rings
Authors:
Dimitri Veras,
Zoe M. Leinhardt,
Siegfried Eggl,
Boris T. Gaensicke
Abstract:
The formation channel of the tens of compact debris discs which orbit white dwarfs (WDs) at a distance of one Solar radius remains unknown. Asteroids that survive the giant branch stellar phases beyond a few au are assumed to be dynamically thrust towards the WD and tidally disrupted within its Roche radius, generating extremely eccentric (e>0.98) rings. Here, we establish that WD radiation compre…
▽ More
The formation channel of the tens of compact debris discs which orbit white dwarfs (WDs) at a distance of one Solar radius remains unknown. Asteroids that survive the giant branch stellar phases beyond a few au are assumed to be dynamically thrust towards the WD and tidally disrupted within its Roche radius, generating extremely eccentric (e>0.98) rings. Here, we establish that WD radiation compresses and circularizes the orbits of super-micron to cm-sized ring constituents to entirely within the WD's Roche radius. We derive a closed algebraic formula which well-approximates the shrinking time as a function of WD cooling age, the physical properties of the star and the physical and orbital properties of the ring particles. The shrinking timescale increases with both particle size and cooling age, yielding age-dependent WD debris disc size distributions.
△ Less
Submitted 22 May, 2015;
originally announced May 2015.
-
A Collisional Origin to Earth's Non-chondritic Composition?
Authors:
Amy Bonsor,
Zoë M. Leinhardt,
Philip J. Carter,
Tim Elliott,
Michael J. Walter,
Sarah T. Stewart
Abstract:
Several lines of evidence indicate a non-chondritic composition for Bulk Earth. If Earth formed from the accretion of chondritic material, its non-chondritic composition, in particular the super-chondritic 142Nd/144Nd and low Mg/Fe ratios, might be explained by the collisional erosion of differentiated planetesimals during its formation. In this work we use an N-body code, that includes a state-of…
▽ More
Several lines of evidence indicate a non-chondritic composition for Bulk Earth. If Earth formed from the accretion of chondritic material, its non-chondritic composition, in particular the super-chondritic 142Nd/144Nd and low Mg/Fe ratios, might be explained by the collisional erosion of differentiated planetesimals during its formation. In this work we use an N-body code, that includes a state-of-the-art collision model, to follow the formation of protoplanets, similar to proto-Earth, from differentiated planetesimals (> 100 km) up to isolation mass (> 0.16 M_Earth). Collisions between differentiated bodies have the potential to change the core-mantle ratio of the accreted protoplanets. We show that sufficient mantle material can be stripped from the colliding bodies during runaway and oligarchic growth, such that the final protoplanets could have Mg/Fe and Si/Fe ratios similar to that of bulk Earth, but only if Earth is an extreme case and the core is assumed to contain 10% silicon by mass. This may indicate an important role for collisional differentiation during the giant impact phase if Earth formed from chondritic material.
△ Less
Submitted 13 October, 2014;
originally announced October 2014.
-
Formation of planetary debris discs around white dwarfs I: Tidal disruption of an extremely eccentric asteroid
Authors:
Dimitri Veras,
Zoe M. Leinhardt,
Amy Bonsor,
Boris T. Gaensicke
Abstract:
25%-50% of all white dwarfs (WDs) host observable and dynamically active remnant planetary systems based on the presence of close-in circumstellar dust and gas and photospheric metal pollution. Currently-accepted theoretical explanations for the origin of this matter include asteroids that survive the star's giant branch evolution at au-scale distances and are subsequently perturbed onto WD-grazin…
▽ More
25%-50% of all white dwarfs (WDs) host observable and dynamically active remnant planetary systems based on the presence of close-in circumstellar dust and gas and photospheric metal pollution. Currently-accepted theoretical explanations for the origin of this matter include asteroids that survive the star's giant branch evolution at au-scale distances and are subsequently perturbed onto WD-grazing orbits following stellar mass loss. In this work we investigate the tidal disruption of these highly-eccentric (e > 0.98) asteroids as they approach and tidally disrupt around the WD. We analytically compute the disruption timescale and compare the result with fully self-consistent numerical simulations of rubble piles by using the N-body code PKDGRAV. We find that this timescale is highly dependent on the orbit's pericentre and largely independent of its semimajor axis. We establish that spherical asteroids readily break up and form highly eccentric collisionless rings, which do not accrete onto the WD without additional forces such as radiation or sublimation. This finding highlights the critical importance of such forces in the physics of WD planetary systems.
△ Less
Submitted 8 September, 2014;
originally announced September 2014.
-
Forming Circumbinary Planets: N-body Simulations of Kepler-34
Authors:
Stefan Lines,
Zoe M. Leinhardt,
Sijme-Jan Paardekooper,
Clement Baruteau,
Philippe Thebault
Abstract:
Observations of circumbinary planets orbiting very close to the central stars have shown that planet formation may occur in a very hostile environment, where the gravitational pull from the binary should be very strong on the primordial protoplanetary disk. Elevated impact velocities and orbit crossings from eccentricity oscillations are the primary contributors towards high energy, potentially de…
▽ More
Observations of circumbinary planets orbiting very close to the central stars have shown that planet formation may occur in a very hostile environment, where the gravitational pull from the binary should be very strong on the primordial protoplanetary disk. Elevated impact velocities and orbit crossings from eccentricity oscillations are the primary contributors towards high energy, potentially destructive collisions that inhibit the growth of aspiring planets. In this work, we conduct high resolution, inter-particle gravity enabled N-body simulations to investigate the feasibility of planetesimal growth in the Kepler-34 system. We improve upon previous work by including planetesimal disk self-gravity and an extensive collision model to accurately handle inter-planetesimal interactions. We find that super-catastrophic erosion events are the dominant mechanism up to and including the orbital radius of Kepler-34(AB)b, making in-situ growth unlikely. It is more plausible that Kepler-34(AB)b migrated from a region beyond 1.5 AU. Based on the conclusions that we have made for Kepler-34 it seems likely that all of the currently known circumbinary planets have also migrated significantly from their formation location with the possible exception of Kepler-47(AB)c.
△ Less
Submitted 3 February, 2014;
originally announced February 2014.
-
Hiding in the Shadows: Searching for Planets in Pre--transitional and Transitional Disks
Authors:
Jack Dobinson,
Zoë M. Leinhardt,
Sarah E. Dodson-Robinson,
Nick A. Teanby
Abstract:
Transitional and pre--transitional disks can be explained by a number of mechanisms. This work aims to find a single observationally detectable marker that would imply a planetary origin for the gap and, therefore, indirectly indicate the presence of a young planet. N-body simulations were conducted to investigate the effect of an embedded planet of one Jupiter mass on the production of instantane…
▽ More
Transitional and pre--transitional disks can be explained by a number of mechanisms. This work aims to find a single observationally detectable marker that would imply a planetary origin for the gap and, therefore, indirectly indicate the presence of a young planet. N-body simulations were conducted to investigate the effect of an embedded planet of one Jupiter mass on the production of instantaneous collisional dust derived from a background planetesimal disk. Our new model allows us to predict the dust distribution and resulting observable markers with greater accuracy than previous work. Dynamical influences from a planet on a circular orbit are shown to enhance dust production in the disk interior and exterior to the planet orbit while removing planetesimals from the the orbit itself creating a clearly defined gap. In the case of an eccentric planet the gap opened by the planet is not as clear as the circular case but there is a detectable asymmetry in the dust disk.
△ Less
Submitted 26 September, 2013; v1 submitted 25 September, 2013;
originally announced September 2013.
-
How not to build Tatooine: the difficulty of in situ formation of circumbinary planets Kepler 16b, Kepler 34b and Kepler 35b
Authors:
Sijme-Jan Paardekooper,
Zoe M. Leinhardt,
Philippe Thebault,
Clement Baruteau
Abstract:
We study planetesimal evolution in circumbinary disks, focusing on the three systems Kepler 16, 34 and 35 where planets have been discovered recently. We show that for circumbinary planetesimals, in addition to secular forcing, eccentricities evolve on a dynamical timescale, which leads to orbital crossings even in the presence of gas drag. This makes the current locations of the circumbinary Kepl…
▽ More
We study planetesimal evolution in circumbinary disks, focusing on the three systems Kepler 16, 34 and 35 where planets have been discovered recently. We show that for circumbinary planetesimals, in addition to secular forcing, eccentricities evolve on a dynamical timescale, which leads to orbital crossings even in the presence of gas drag. This makes the current locations of the circumbinary Kepler planets hostile to planetesimal accretion. We then present results from simulations including planetesimal formation and dust accretion, and show that even in the most favourable case of 100% efficient dust accretion, in situ growth starting from planetesimals smaller than ~10 km is difficult for Kepler 16b, Kepler 34b and Kepler 35b. These planets were likely assembled further out in the disk, and migrated inward to their current location.
△ Less
Submitted 15 June, 2012;
originally announced June 2012.
-
Tidal disruption of satellites and formation of narrow rings
Authors:
Zoe M. Leinhardt,
Gordon I. Ogilvie,
Henrik N. Latter,
Eiichiro Kokubo
Abstract:
In this paper we investigate the formation of narrow planetary rings such as those found around Uranus and Saturn through the tidal disruption of a weak, gravitationally bound satellite that migrates within its Roche limit. Using $N$-body simulations, we study the behaviour of rubble piles placed on circular orbits at different distances from a central planet. We consider both homogeneous satellit…
▽ More
In this paper we investigate the formation of narrow planetary rings such as those found around Uranus and Saturn through the tidal disruption of a weak, gravitationally bound satellite that migrates within its Roche limit. Using $N$-body simulations, we study the behaviour of rubble piles placed on circular orbits at different distances from a central planet. We consider both homogeneous satellites and differentiated bodies containing a denser core. We show that the Roche limit for a rubble pile is closer to the planet than for a fluid body of the same mean density. The Roche limit for a differentiated body is also closer to the planet than for a homogeneous satellite of the same mean density. Within its Roche limit, a homogeneous satellite totally disrupts and forms a narrow ring. The initial stages of the disruption are similar to the evolution of a viscous fluid ellipsoid, which can be computed semi-analytically. On the other hand, when a differentiated satellite is just within the Roche limit only the mantle is disrupted. This process is similar to Roche-lobe overflow in interacting binary stars and produces two narrow rings on either side of a remnant satellite. We argue that the Uranian rings, and possibly their shepherd satellites, could have been formed through the tidal disruption of a number of protomoons that were formed inside the corotation radius of Uranus and migrated slowly inwards as a result of tidal interaction with the planet.
△ Less
Submitted 12 June, 2012;
originally announced June 2012.
-
Collisions between Gravity-Dominated Bodies: 2. The Diversity of Impact Outcomes during the End Stage of Planet Formation
Authors:
S. T. Stewart,
Z. M. Leinhardt
Abstract:
Numerical simulations of the stochastic end stage of planet formation typically begin with a population of embryos and planetesimals that grow into planets by merging. We analyzed the impact parameters of collisions leading to the growth of terrestrial planets from recent $N$-body simulations that assumed perfect merging and calculated more realistic outcomes using a new analytic collision physics…
▽ More
Numerical simulations of the stochastic end stage of planet formation typically begin with a population of embryos and planetesimals that grow into planets by merging. We analyzed the impact parameters of collisions leading to the growth of terrestrial planets from recent $N$-body simulations that assumed perfect merging and calculated more realistic outcomes using a new analytic collision physics model. We find that collision outcomes are diverse and span all possible regimes: hit-and-run, merging, partial accretion, partial erosion, and catastrophic disruption. The primary outcomes of giant impacts between planetary embryos are approximately evenly split between partial accretion, graze-and-merge, and hit-and-run events. To explore the cumulative effects of more realistic collision outcomes, we modeled the growth of individual planets with a Monte Carlo technique using the distribution of impact parameters from $N$-body simulations. We find that fewer planets reached masses $>0.7 M_{\rm Earth}$ using the collision physics model compared to simulations that assumed every collision results in perfect merging. For final planets with masses $>0.7 M_{\rm Earth}$, 60% are enriched in their core-to-mantle mass fraction by >10% compared to the initial embryo composition. Fragmentation during planet formation produces significant debris ($\sim15$% of the final mass) and occurs primarily by erosion of the smaller body in partial accretion and hit-and-run events. In partial accretion events, the target body grows by preferentially accreting the iron core of the projectile and the escaping fragments are derived primarily from the silicate mantles of both bodies. Thus, the bulk composition of a planet can evolve via stochastic giant impacts.
△ Less
Submitted 2 May, 2012; v1 submitted 21 September, 2011;
originally announced September 2011.
-
Collisions Between Gravity-Dominated Bodies: 1. Outcome Regimes and Scaling Laws
Authors:
Zoë M. Leinhardt,
Sarah T. Stewart
Abstract:
Collisions are the core agent of planet formation. In this work, we derive an analytic description of the dynamical outcome for any collision between gravity-dominated bodies. We conduct high-resolution simulations of collisions between planetesimals; the results are used to isolate the effects of different impact parameters on collision outcome. During growth from planetesimals to planets, collis…
▽ More
Collisions are the core agent of planet formation. In this work, we derive an analytic description of the dynamical outcome for any collision between gravity-dominated bodies. We conduct high-resolution simulations of collisions between planetesimals; the results are used to isolate the effects of different impact parameters on collision outcome. During growth from planetesimals to planets, collision outcomes span multiple regimes: cratering, merging, disruption, super-catastrophic disruption, and hit-and-run events. We derive equations (scaling laws) to demarcate the transition between collision regimes and to describe the size and velocity distributions of the post-collision bodies. The scaling laws are used to calculate maps of collision outcomes as a function of mass ratio, impact angle, and impact velocity, and we discuss the implications of the probability of each collision regime during planet formation.
The analytic collision model presented in this work will significantly improve the physics of collisions in numerical simulations of planet formation and collisional evolution. (abstract abridged)
△ Less
Submitted 5 January, 2012; v1 submitted 29 June, 2011;
originally announced June 2011.
-
The Formation of the Collisional Family around the Dwarf Planet Haumea
Authors:
Zoe M. Leinhardt,
Robert A. Marcus,
Sarah T. Stewart
Abstract:
Haumea, a rapidly rotating elongated dwarf planet (~ 1500 km in diameter), has two satellites and is associated with a "family" of several smaller Kuiper Belt objects (KBOs) in similar orbits. All members of the Haumea system share a water ice spectral feature that is distinct from all other KBOs. The relative velocities between the Haumea family members are too small to have formed by catastrop…
▽ More
Haumea, a rapidly rotating elongated dwarf planet (~ 1500 km in diameter), has two satellites and is associated with a "family" of several smaller Kuiper Belt objects (KBOs) in similar orbits. All members of the Haumea system share a water ice spectral feature that is distinct from all other KBOs. The relative velocities between the Haumea family members are too small to have formed by catastrophic disruption of a large precursor body, which is the process that formed families around much smaller asteroids in the Main Belt. Here we show that all of the unusual characteristics of the Haumea system are explained by a novel type of giant collision: a graze-and-merge impact between two comparably sized bodies. The grazing encounter imparted the high angular momentum that spun off fragments from the icy crust of the elongated merged body. The fragments became satellites and family members. Giant collision outcomes are extremely sensitive to the impact parameters. Compared to the Main Belt, the largest bodies in the Kuiper Belt are more massive and experience slower velocity collisions; hence, outcomes of giant collisions are dramatically different between the inner and outer solar system. The dwarf planets in the Kuiper Belt record an unexpectedly large number of giant collisions, requiring a special dynamical event at the end of solar system formation.
△ Less
Submitted 30 March, 2010;
originally announced March 2010.
-
Planetesimal collisions in binary systems
Authors:
S. -J. Paardekooper,
Z. M. Leinhardt
Abstract:
We study the collisional evolution of km-sized planetesimals in tight binary star systems to investigate whether accretion towards protoplanets can proceed despite the strong gravitational perturbations from the secondary star. The orbits of planetesimals are numerically integrated in two dimensions under the influence of the two stars and gas drag. The masses and orbits of the planetesimals are…
▽ More
We study the collisional evolution of km-sized planetesimals in tight binary star systems to investigate whether accretion towards protoplanets can proceed despite the strong gravitational perturbations from the secondary star. The orbits of planetesimals are numerically integrated in two dimensions under the influence of the two stars and gas drag. The masses and orbits of the planetesimals are allowed to evolve due to collisions with other planetesimals and accretion of collisional debris. In addition, the mass in debris can evolve due to planetesimal-planetesimal collisions and the creation of new planetesimals. We show that it is possible in principle for km-sized planetesimals to grow by two orders of magnitude in size if the efficiency of planetesimal formation is relatively low. We discuss the limitations of our two-dimensional approach.
△ Less
Submitted 18 January, 2010;
originally announced January 2010.
-
The Validity of the Super-Particle Approximation during Planetesimal Formation
Authors:
Hanno Rein,
Geoffroy Lesur,
Zoe M. Leinhardt
Abstract:
The formation mechanism of planetesimals in protoplanetary discs is hotly debated. Currently, the favoured model involves the accumulation of meter-sized objects within a turbulent disc, followed by a phase of gravitational instability. At best one can simulate a few million particles numerically as opposed to the several trillion meter-sized particles expected in a real protoplanetary disc. The…
▽ More
The formation mechanism of planetesimals in protoplanetary discs is hotly debated. Currently, the favoured model involves the accumulation of meter-sized objects within a turbulent disc, followed by a phase of gravitational instability. At best one can simulate a few million particles numerically as opposed to the several trillion meter-sized particles expected in a real protoplanetary disc. Therefore, single particles are often used as super-particles to represent a distribution of many smaller particles. It is assumed that small scale phenomena do not play a role and particle collisions are not modeled. The super-particle approximation can only be valid in a collisionless or strongly collisional system, however, in many recent numerical simulations this is not the case.
In this work we present new results from numerical simulations of planetesimal formation via gravitational instability. A scaled system is studied that does not require the use of super-particles. We find that the scaled particles can be used to model the initial phases of clumping if the properties of the scaled particles are chosen such that all important timescales in the system are equivalent to what is expected in a real protoplanetary disc. Constraints are given for the number of particles needed in order to achieve numerical convergence.
We compare this new method to the standard super-particle approach. We find that the super-particle approach produces unreliable results that depend on artifacts such as the gravitational softening in both the requirement for gravitational collapse and the resulting clump statistics. Our results show that short range interactions (collisions) have to be modelled properly.
△ Less
Submitted 8 January, 2010; v1 submitted 31 December, 2009;
originally announced January 2010.
-
Planetesimals to Protoplanets II: Effect of Debris on Terrestrial Planet Formation
Authors:
Z. M. Leinhardt,
D. C. Richardson,
G. Lufkin,
J. Haseltine
Abstract:
In this paper we extend our numerical method for simulating terrestrial planet formation from Leinhardt and Richardson (2005) to include dynamical friction from the unresolved debris component. In the previous work we implemented a rubble pile planetesimal collision model into direct N-body simulations of terrestrial planet formation. The new collision model treated both accretion and erosion of…
▽ More
In this paper we extend our numerical method for simulating terrestrial planet formation from Leinhardt and Richardson (2005) to include dynamical friction from the unresolved debris component. In the previous work we implemented a rubble pile planetesimal collision model into direct N-body simulations of terrestrial planet formation. The new collision model treated both accretion and erosion of planetesimals but did not include dynamical friction from debris particles smaller than the resolution limit for the simulation. By extending our numerical model to include dynamical friction from the unresolved debris, we can simulate the dynamical effect of debris produced during collisions and can also investigate the effect of initial debris mass on terrestrial planet formation. We find that significant initial debris mass, 10% or more of the total disk mass, changes the mode of planetesimal growth. Specifically, planetesimals in this situation do not go through a runaway growth phase. Instead they grow concurrently, similar to oligarchic growth. In addition to including the dynamical friction from the unresolved debris, we have implemented particle tracking as a proxy for monitoring compositional mixing. Although there is much less mixing due to collisions and gravitational scattering when dynamical friction of the background debris is included, there is significant inward migration of the largest protoplanets in the most extreme initial conditions.
△ Less
Submitted 13 March, 2009;
originally announced March 2009.
-
Full Numerical Simulations of Catastrophic Small Body Collisions
Authors:
Zoe M. Leinhardt,
Sarah T. Stewart
Abstract:
The outcome of collisions between small icy bodies, such as Kuiper belt objects, is poorly understood and yet a critical component of the evolution of the trans-Neptunian region. The expected physical properties of outer solar system materials (high porosity, mixed ice-rock composition, and low material strength) pose significant computational challenges. We present results from catastrophic sma…
▽ More
The outcome of collisions between small icy bodies, such as Kuiper belt objects, is poorly understood and yet a critical component of the evolution of the trans-Neptunian region. The expected physical properties of outer solar system materials (high porosity, mixed ice-rock composition, and low material strength) pose significant computational challenges. We present results from catastrophic small body collisions using a new hybrid hydrocode to $N$-body code computational technique. This method allows detailed modeling of shock propagation and material modification as well as gravitational reaccumulation. Here, we consider a wide range of material strengths to span the possible range of Kuiper belt objects. We find that the shear strength of the target is important in determining the collision outcome for 2 to 50-km radius bodies, which are traditionally thought to be in a pure gravity regime. The catastrophic disruption and dispersal criteria, $Q_D^*$, can vary by up to a factor of three between strong crystalline and weak aggregate materials. The material within the largest reaccumulated remnants experiences a wide range of shock pressures. The dispersal and reaccumulation process results in the material on the surfaces of the largest remnants having experienced a wider range of shock pressures compared to material in the interior. Hence, depending on the initial structure and composition, the surface materials on large, reaccumulated bodies in the outer solar system may exhibit complex spectral and albedo variations. Finally, we present revised catastrophic disruption criteria for a range of impact velocities and material strengths for outer solar system bodies.
△ Less
Submitted 2 November, 2008;
originally announced November 2008.
-
Physical effects of collisions in the Kuiper belt
Authors:
Zoe M. Leinhardt,
Sarah T. Stewart,
Peter H. Schultz
Abstract:
Collisions are a major modification process over the history of the Kuiper Belt. Recent work illuminates the complex array of possible outcomes of individual collisions onto porous, volatile bodies. The cumulative effects of such collisions on the surface features, composition, and internal structure of Kuiper Belt Objects are not yet known. In this chapter, we present the current state of knowl…
▽ More
Collisions are a major modification process over the history of the Kuiper Belt. Recent work illuminates the complex array of possible outcomes of individual collisions onto porous, volatile bodies. The cumulative effects of such collisions on the surface features, composition, and internal structure of Kuiper Belt Objects are not yet known. In this chapter, we present the current state of knowledge of the physics of cratering and disruptive collisions in KBO analog materials. We summarize the evidence for a rich collisional history in the Kuiper Belt and present the range possible physical modifications on individual objects. The question of how well present day bodies represent primordial planetesimals can be addressed through future studies of the coupled physical and collisional evolution of Kuiper Belt Objects.
△ Less
Submitted 27 May, 2007;
originally announced May 2007.
-
Direct N-body Simulations of Rubble Pile Collisions
Authors:
Z. M. Leinhardt,
D. C. Richardson,
T. Quinn
Abstract:
There is increasing evidence that many km-sized bodies in the Solar System are piles of rubble bound together by gravity. We present results from a project to map the parameter space of collisions between km-sized spherical rubble piles. The results will assist in parameterization of collision outcomes for Solar System formation models and give insight into fragmentation scaling laws. We use a d…
▽ More
There is increasing evidence that many km-sized bodies in the Solar System are piles of rubble bound together by gravity. We present results from a project to map the parameter space of collisions between km-sized spherical rubble piles. The results will assist in parameterization of collision outcomes for Solar System formation models and give insight into fragmentation scaling laws. We use a direct numerical method to evolve the positions and velocities of the rubble pile particles under the constraints of gravity and physical collisions. We test the dependence of the collision outcomes on impact parameter and speed, impactor spin, mass ratio, and coefficient of restitution. Speeds are kept low (< 10 m/s, appropriate for dynamically cool systems such as the primordial disk during early planet formation) so that the maximum strain on the component material does not exceed the crushing strength. We compare our results with analytic estimates and hydrocode simulations. Off-axis collisions can result in fast-spinning elongated remnants or contact binaries while fast collisions result in smaller fragments overall. Clumping of debris escaping from the remnant can occur, leading to the formation of smaller rubble piles. In the cases we tested, less than 2% of the system mass ends up orbiting the remnant. Initial spin can reduce or enhance collision outcomes, depending on the relative orientation of the spin and orbital angular momenta. We derive a relationship between impact speed and angle for critical dispersal of mass in the system. We find that our rubble piles are relatively easy to disperse, even at low impact speed, suggesting that greater dissipation is required if rubble piles are the true progenitors of protoplanets.
△ Less
Submitted 19 August, 1999;
originally announced August 1999.