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.

Metals observed in the atmospheres of white dwarfs suggest that many have recently accreted planetary bodies. In some cases, the compositions observed suggest the accretion of material dominantly from the core (or the mantle) of a differentiated planetary body. Collisions between differentiated exoplanetesimals produce such fragments. In this work, we take advantage of the large numbers of white dwarfs where at least one siderophile (core-loving) and one lithophile (rock-loving) species have been detected to assess how commonly exoplanetesimals differentiate. We utilise N-body simulations that track the fate of core and mantle material during the collisional evolution of planetary systems to show that most remnants of differentiated planetesimals retain core fractions similar to their parents, whilst some are extremely core-rich or mantle-rich. Comparison with the white dwarf data for calcium and iron indicates that the data are consistent with a model in which $66^{+4}_{-6}\%$ have accreted the remnants of differentiated planetesimals, whilst $31^{+5}_{-5}\%$ have Ca/Fe abundances altered by the effects of heating (although the former can be as high as $100\%$, if heating is ignored). These conclusions assume pollution by a single body and that collisional evolution retains similar features across diverse planetary systems. These results imply that both collisions and differentiation are key processes in exoplanetary systems. We highlight the need for a larger sample of polluted white dwarfs with precisely determined metal abundances to better understand the process of differentiation in exoplanetary systems.

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.

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.

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.

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.

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.

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.

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.

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 $\mu$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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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)

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.

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.

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.

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.

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.

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.

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.