BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Derek C. Richardson

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

ResearcherUnclaimed source record
astro-ph.EPcond-mat.softastro-ph.GAastro-ph.HEastro-ph.IMastro-ph.SRcond-mat.mtrl-sciphysics.comp-phphysics.geo-ph

Catalog footprint

What is connected

20works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

20 published item(s)

preprint2026arXiv

A change of the rotation period of asteroid (65803) Didymos caused by the DART impact

On 26 September 2022, the NASA Double Asteroid Redirection Test (DART) spacecraft impacted Dimorphos, the secondary component of the binary asteroid (65803)~Didymos. This experiment tested the Kinetic Impactor technology for diverting dangerous asteroids. Due to the impact, the binary system's angular momentum has changed, resulting in a significant change in the orbital period of Dimorphos. Precise values of the pre- and post-impact orbital periods were derived from a large set of photometric light curves measured for the Didymos-Dimorphos system during six apparitions from 2003 to 2023. We used these data to detect a possible change in the rotation period of the primary as a consequence of the impact. We analyzed the binary system's light curves using the binary asteroid light curve decomposition method. We selected parts of the light curves covering orbital phases outside mutual events, which represent the primary rotational light curves. We applied the light curve inversion method to construct a convex shape model of Didymos and determine its rotation period before and after the impact. These two periods were treated as independent free parameters of the modeling. We found a value of 2.2603891 +/- 0.0000002 h for the pre-impact period and 2.260440 +/- 0.000008 h for the post-impact period. Their difference 0.18 +/- 0.03 s is small yet significant, indicating that the rotation of Didymos became slower after the DART impact. The most plausible physical explanation is Didymos's post-impact reshaping, making its shape more oblate.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Dynamical Evolution of the Didymos-Dimorphos Binary Asteroid as Rubble Piles following the DART Impact

Previous efforts have modeled the Didymos system as two irregularly shaped rigid bodies, although it is likely that one or both components are in fact rubble piles. Here, we relax the rigid-body assumption to quantify how this affects the spin and orbital dynamics of the system following the DART impact. Given known fundamental differences between our simulation codes, we find that faster rigid-body simulations produce nearly the same result as rubble-pile models in scenarios with a moderate value for the momentum enhancement factor, $β$ ($β{\sim}3$) and an ellipsoidal secondary. This indicates that the rigid-body approach is likely adequate for propagating the post-impact dynamics necessary to meet DART Mission requirements. Although, if Dimorphos has a highly-irregular shape or structure, or if $β$ is unexpectedly large, then rubble-pile effects may become important. If Dimorphos's orbit and spin state are sufficiently excited, then surface particle motion is also possible. However, these simulations are limited in their resolution and range of material parameters, so they serve as a demonstration of principle, and Future work is required to fully understand the likelihood and magnitude of surface motion.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Experimentally measuring rolling and sliding in three-dimensional dense granular packings

We experimentally measure a three-dimensional (3D) granular system's reversibility under cyclic compression. We image the grains using a refractive-index-matched fluid, then analyze the images using the artificial intelligence of variational autoencoders. These techniques allow us to track all the grains' translations and 3D rotations with accuracy sufficient to infer sliding and rolling displacements. Our observations reveal unique roles played by 3D rotational motions in granular flows. We find that rotations and contact-point motion dominate the dynamics in the bulk, far from the perturbation's source. Furthermore, we determine that 3D rotations are irreversible under cyclic compression. Consequently, contact-point sliding, which is dissipative, accumulates throughout the cycle. Using numerical simulations whose accuracy our experiment supports, we discover that much of the dissipation occurs in the bulk, where grains rotate more than they translate. Our observations suggest that the analysis of 3D rotations is needed for understanding granular materials' unique and powerful ability to absorb and dissipate energy.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Predictions for the Dynamical States of the Didymos System before and after the Planned DART Impact

NASA's Double Asteroid Redirection Test (DART) spacecraft is planned to impact the natural satellite of (65803) Didymos, Dimorphos, around 23:14 UTC on 26 September 2022, causing a reduction in its orbital period that will be measurable with ground-based observations. This test of kinetic impactor technology will provide the first estimate of the momentum transfer enhancement factor $β$ at a realistic scale, wherein ejecta from the impact provides an additional deflection to the target. Earth-based observations, the LICIACube spacecraft (to be detached from DART prior to impact), and ESA's follow-up Hera mission to launch in 2024, will provide additional characterization of the deflection test. Together Hera and DART comprise the Asteroid Impact and Deflection Assessment (AIDA) cooperation between NASA and ESA. Here the predicted dynamical states of the binary system upon arrival and after impact are presented. The assumed dynamically relaxed state of the system will be excited by the impact, leading to an increase in eccentricity and slight tilt of the orbit together with enhanced libration of Dimorphos with amplitude dependent on the currently poorly known target shape. Free rotation around the moon's long axis may also be triggered and the orbital period will experience variations from seconds to minutes over timescales of days to months. Shape change of either body due to cratering or mass wasting triggered by crater formation and ejecta may affect $β$ but can be constrained through additional measurements. Both BYORP and gravity tides may cause measurable orbital changes on the timescale of Hera's rendezvous.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Rotation-induced granular motion on the secondary component of binary asteroids: Application to the DART impact on Dimorphos

NASA's Double Asteroid Redirection Test (DART) mission will kinetically impact Dimorphos, the secondary component of the Didymos binary asteroid system, which will excite Dimorphos's dynamical state and lead to significant libration about the synchronous state and possibly chaotic non-principal axis rotation. Although this particular outcome is human caused, many other secondary components of binary systems are also prone to such exotic spin states. For a satellite in an excited spin state, the time-varying tidal and rotational environment can lead to significant surface accelerations. Depending on the circumstances, this mechanism may drive granular motion on the surface of the secondary. We modeled the dynamical evolution of a Didymos-like binary asteroid system using a fully coupled, three-dimensional simulation code. Then, we computed the time-varying gravitational and rotational accelerations felt over the entire surface resulting from the secondary's perturbed dynamical state. We find that an excited spin and orbit can induce large changes in the effective surface slope, potentially triggering granular motion and surface refreshment. However, for the case of the DART impact, this effect is highly dependent on many unknowns, such as Dimorphos's detailed shape, bulk density, surface geology, and the momentum transferred. Aside from the Didymos system and the DART mission, this effect also has important implications for binary systems in general.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Surface Properties of Near-Sun Asteroids

Near-Earth Asteroids (NEAs) with small perihelion distances reach sub-solar temperatures of > 1000 K. They are hypothesized to undergo "super-catastrophic" disruption, potentially caused by near-Sun processes such as thermal cracking, spin-up, meteoroid impacts, and subsurface volatile release; all of which are likely to cause surface alteration, which may change the spectral slope of the surface. We attempted to observe 35 of the 53 known near-Sun asteroids with q < 0.15 au from January 2017 to March 2020 to search for trends related to near-Sun processes. We report the optical colors and spectral slopes of 22 objects that we successfully observed and the measured rotation periods for three objects. We find the distribution of colors to be overall bluer than the color distribution of NEAs, though there is large overlap. We attribute large scatter to unknown dynamical histories and compositions for individual objects, as well as competing surface altering processes. We also investigated potential correlations between colors and other properties (e.g., perihelion distance, Tisserand parameter, rotation period), and searched for evidence of activity. Finally, we have compiled all known physical and dynamical properties of these objects, including probabilistic source regions and dwell times with q < 0.15 au.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Critical spin periods of sub-km-sized cohesive rubble-pile asteroids: dependencies on material parameters

In this work, we employ a soft-sphere discrete element method with a cohesion implementation to model the dynamical process of sub-km-sized cohesive rubble piles under continuous spinup. The dependencies of critical spin periods $T_c$ on several material parameters for oblate rubble piles with different bulk diameters $D$ are explored. Our numerical simulations show that both the increase of interparticle cohesion and particle shape parameter in our model can strengthen the bodies, especially for the smaller ones. In addition, we find there exists some critical diameter $D_{cri,ρ}$ at which the variation trend of $T_c$ with the bulk density $ρ$ reverses. Though a greater static friction coefficient $μ_S$ can strengthen the body, this effect attains a minimum at a critical diameter $D_{cri,ϕ}$ close to $D_{cri,ρ}$. The continuum theory (analytical method) is used for comparison and two equivalent critical diameters are obtained. The numerical results were fitted with the analytical method and the ratio of the interparticle cohesion $c$ to the bulk cohesion $C$ is estimated to be roughly 88.3. We find this ratio keeps constant for different $c$ and $ρ$, while it strongly depends on the friction angle $ϕ$. Also, our numerical results further show that the dependency of $T_c$ on $ϕ$ is opposite from that predicted by the continuum theory when $D$ < $D_{cri,ϕ}$. Finally, we find that the two critical diameters happen to be close to the diameter when the mean normal stress of the body equals zero, which is the separation between the compressive regime and the tensile regime.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Superparticle Method for Simulating Collisions

For problems in astrophysics, planetary science and beyond, numerical simulations are often limited to simulating fewer particles than in the real system. To model collisions, the simulated particles (aka superparticles) need to be inflated to represent a collectively large collisional cross section of real particles. Here we develop a superparticle-based method that replicates the kinetic energy loss during real-world collisions, implement it in an $N$-body code and test it. The tests provide interesting insights into dynamics of self gravitating collisional systems. They show how particle systems evolve over several free fall timescales to form central concentrations and equilibrated outer shells. The superparticle method can be extended to account for the accretional growth of objects during inelastic mergers.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Hektor - an exceptional D-type family among Jovian Trojans

In this work, we analyze Jovian Trojans in the space of suitable resonant elements and we identify clusters of possible collisional origin by two independent methods: the hierarchical clustering and a so-called "randombox". Compared to our previous work (Brož and Rozehnal 2011), we study a twice larger sample. Apart from Eurybates, Ennomos and $1996\,\rm RJ$ families, we have found three more clusters --- namely families around asteroids (20961)~Arkesilaos, (624)~Hektor in the $L_4$ libration zone and (247341)~$2001\,\rm UV_{209}$ in $L_5$. The families fulfill our stringent criteria, i.e. a high statistical significance, an albedo homogeneity and a steeper size-frequency distribution than that of background. In order to understand their nature, we simulate their long term collisional evolution with the Boulder code (Morbidelli et al. 2009) and dynamical evolution using a modified SWIFT integrator (Levison and Duncan, 1994). Within the framework of our evolutionary model, we were able to constrain the the age of the Hektor family to be either 1 to 4 Gyr or, less likely, 0.1 to 2.5 Gyr, depending on initial impact geometry. Since (624) Hektor itself seems to be a bilobed--shape body with a satellite (Marchis et al. 2014), i.e. an exceptional object, we address its association with the D--type family and we demonstrate that the moon and family could be created during a single impact event. We simulated the cratering event using a Smoothed Particle Hydrodynamics (SPH, Benz and Asphaug, 1994). This is also the first case of a family associated with a D--type parent body.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Role of the Kozai-Lidov Mechanism in Black Hole Binary Mergers in Galactic Centers

In order to understand the rate of merger of stellar-mass black hole binaries (BHBs) by gravitational wave (GW) emission it is important to determine the major pathways to merger. We use numerical simulations to explore the evolution of BHBs inside the radius of influence of supermassive black holes (SMBHs) in galactic centers. In this region the evolution of binaries is dominated by perturbations from the central SMBH. In particular, as first pointed out by Antonini and Perets, the Kozai-Lidov (KL) mechanism trades relative inclination of the BHB to the SMBH for eccentricity of the BHB, and for some orientations can bring the BHB to an eccentricity near unity. At very high eccentricities, GW emission from the BHB can become efficient, causing the members of the BHB to coalesce. We use a novel combination of two N-body codes to follow this evolution. We are forced to simulate small systems to follow the behavior accurately. We have completed 400 simulations that range from $\sim$ 300 stars around a $10^{3}$ M$_{\odot}$ black hole to $\sim$ 4500 stars around a $10^{4}$ M$_{\odot}$ black hole. These simulations are the first to follow the internal orbit of a binary near a SMBH while also following the changes to its external orbit self-consistently. We find that this mechanism could produce mergers at a maximum rate per volume of $\sim 100$ Gpc$^{-3}$ yr$^{-1}$ or considerably less if the inclination oscillations of the binary remain constant as the BHB inclination to the SMBH changes, or if the binary black hole fraction is small.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Collisional Formation and Modeling of Asteroid Families

In the last decade, thanks to the development of sophisticated numerical codes, major breakthroughs have been achieved in our understanding of the formation of asteroid families by catastrophic disruption of large parent bodies. In this review, we describe numerical simulations of asteroid collisions that reproduced the main properties of families, accounting for both the fragmentation of an asteroid at the time of impact and the subsequent gravitational interactions of the generated fragments. The simulations demonstrate that the catastrophic disruption of bodies larger than a few hundred meters in diameter leads to the formation of large aggregates due to gravitational reaccumulation of smaller fragments, which helps explain the presence of large members within asteroid families. Thus, for the first time, numerical simulations successfully reproduced the sizes and ejection velocities of members of representative families. Moreover, the simulations provide constraints on the family dynamical histories and on the possible internal structure of family members and their parent bodies.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Coupled Spin and Shape Evolution of Small Rubble-Pile Asteroids: Self-Limitation of the YORP Effect

We present the first self-consistent simulations of the coupled spin-shape evolution of small gravitational aggregates under the influence of the YORP effect. Because of YORP's sensitivity to surface topography, even small centrifugally driven reconfigurations of aggregates can alter the YORP torque dramatically, resulting in spin evolution that can differ qualitatively from the rigid-body prediction. One third of our simulations follow a simple evolution described as a modified YORP cycle. Two-thirds exhibit one or more of three distinct behaviors---stochastic YORP, self-governed YORP, and stagnating YORP---which together result in YORP self-limitation. Self-limitation confines rotation rates of evolving aggregates to far narrower ranges than those expected in the classical YORP cycle, greatly prolonging the times over which objects can preserve their sense of rotation. Simulated objects are initially randomly packed, disordered aggregates of identical spheres in rotating equilibrium, with low internal angles of friction. Their shape evolution is characterized by rearrangement of the entire body, including the deep interior. They do not evolve to axisymmetric top shapes with equatorial ridges. Mass loss occurs in one-third of the simulations, typically in small amounts from the ends of a prolate-triaxial body. We conjecture that YORP self-limitation may inhibit formation of top-shapes, binaries, or both, by restricting the amount of angular momentum that can be imparted to a deformable body. Stochastic YORP, in particular, will affect the evolution of collisional families whose orbits drift apart under the influence of Yarkovsky forces, in observable ways.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Low-speed Impact Simulations into Regolith in Support of Asteroid Sampling Mechanism Design I.: Comparison with 1-g Experiments

This study is carried out in the framework of sample-return missions to asteroids that use a low-speed projectile as the primary component of its sampling mechanism (e.g., JAXA's Hayabusa and Hayabusa2 missions). We perform numerical simulations of such impacts into granular materials using different projectile shapes under Earth's gravity. We then compare the amounts of ejected mass obtained in our simulations against what was found in experiments that used similar setups, which allows us to validate our numerical approach. For the targets, we consider 2 different monodisperse grain-diameter sizes: 5 mm and 3 mm. The impact speed of the projectile is 11 m s$^{-1}$ directed downward, perpendicular to the surface of the targets. Using an implementation of the soft-sphere discrete element method (SSDEM) in the $N$-Body gravity tree code PKDGRAV, previously validated in the context of low-speed impacts into sintered glass bead agglomerates, we find a noticeable dependence of the amount of ejected mass on the projectile shape. As found in experiments, in the case of the larger target grain size (5 mm), a conically shaped projectile ejects a greater amount of mass than do projectiles of other shapes, including disks and spheres. We then find that numerically the results are sensitive to the normal coefficient of restitution of the grains, especially for impacts into targets comprised of smaller grains (3 mm). We also find that static friction plays a more important role for impacts into targets comprised of the larger grains. As a preliminary demonstration, one of these considered setups is simulated in a microgravity environment. As expected, a reduction in gravity increases both the amount of ejected mass and the timescale of the impact process. Our methodology is also adaptable to the conditions of sampling mechanisms included in specific mission designs.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Numerical predictions of surface effects during the 2029 close approach of asteroid 99942 Apophis

Asteroid (99942) Apophis' close approach in 2029 will be one of the most significant small-body encounter events in the near future and offers a good opportunity for in situ exploration to determine the asteroid's surface properties and measure any tidal effects that might alter its regolith configuration. Resurfacing mechanics has become a new focus for asteroid researchers due to its important implications for interpreting surface observations, including space weathering effects. This paper provides a prediction for the tidal effects during the 2029 encounter, with an emphasis on whether surface refreshing due to regolith movement will occur. The potential shape modification of the object due to the tidal encounter is first confirmed to be negligibly small with systematic simulations, thus only the external perturbations are taken into account for this work (despite this, seismic shaking induced by shifting blocks might still play a weak role and we will look into this mechanism in future work). A two-stage approach is developed to model the responses of asteroid surface particles (the regolith) based on the soft-sphere implementation of the parallel N-body gravity tree code pkdgrav. A full-body model of Apophis is sent past the Earth on the predicted trajectory to generate the data of all forces acting at a target point on the surface. A sandpile constructed in the local frame is then used to approximate the regolith materials; all the forces the sandpile feels during the encounter are imposed as external perturbations to mimic the regolith's behavior in the full scenario. The local mechanical environment on the asteroid surface is represented in detail, leading to an estimation of the change in global surface environment due to the encounter. Typical patterns of perturbation are presented that depend on the asteroid orientation and sense of rotation at perigee.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Brazil-nut effect and its application to asteroids

Out of the handful of asteroids that have been imaged, some have distributions of blocks that are not easily explained. In this paper, we investigate the possibility that seismic shaking leads to the size sorting of particles in asteroids. In particular, we focus on the so-called Brazil Nut Effect (BNE) that separates large particles from small ones under vibrations. We study the BNE over a wide range of parameters by using the N-body code PKDGRAV, and find that the effect is largely insensitive to the coefficients of restitution, but sensitive to friction constants and oscillation speeds. Agreeing with the previous results, we find that convection drives the BNE, where the intruder rises to the top of the particle bed. For the wide-cylinder case, we also observe a "whale" effect, where the intruder follows the convective current and does not stay at the surface. We show that the non-dimensional critical conditions for the BNE agree well with previous studies. We also show that the BNE is scalable for low-gravity environments and that the rise speed of an intruder is proportional to the square root of the gravitational acceleration. Finally, we apply the critical conditions to observed asteroids, and find that the critical oscillation speeds are comparable to the seismic oscillation speeds that are expected from non-destructive impacts.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Numerical simulations of granular dynamics II. Particle dynamics in a shaken granular material

Surfaces of planets and small bodies of our Solar System are often covered by a layer of granular material that can range from a fine regolith to a gravel-like structure of varying depths. Therefore, the dynamics of granular materials are involved in many events occurring during planetary and small-body evolution thus contributing to their geological properties. We demonstrate that the new adaptation of the parallel N-body hard-sphere code pkdgrav has the capability to model accurately the key features of the collective motion of bidisperse granular materials in a dense regime as a result of shaking. As a stringent test of the numerical code we investigate the complex collective ordering and motion of granular material by direct comparison with laboratory experiments. We demonstrate that, as experimentally observed, the scale of the collective motion increases with increasing small-particle additive concentration. We then extend our investigations to assess how self-gravity and external gravity affect collective motion. In our reduced-gravity simulations both the gravitational conditions and the frequency of the vibrations roughly match the conditions on asteroids subjected to seismic shaking, though real regolith is likely to be much more heterogeneous and less ordered than in our idealised simulations. We also show that collective motion can occur in a granular material under a wide range of inter-particle gravity conditions and in the absence of an external gravitational field. These investigations demonstrate the great interest of being able to simulate conditions that are to relevant planetary science yet unreachable by Earth-based laboratory experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Numerical simulations of granular dynamics. I. Hard-sphere discrete element method and tests

We present a new particle-based (discrete element) numerical method for the simulation of granular dynamics, with application to motions of particles on small solar system body and planetary surfaces. The method employs the parallel N-body tree code pkdgrav to search for collisions and compute particle trajectories. Collisions are treated as instantaneous point-contact events between rigid spheres. Particle confinement is achieved by combining arbitrary combinations of four provided wall primitives, namely infinite plane, finite disk, infinite cylinder, and finite cylinder, and degenerate cases of these. Various wall movements, including translation, oscillation, and rotation, are supported. We provide full derivations of collision prediction and resolution equations for all geometries and motions. Several tests of the method are described, including a model granular "atmosphere" that achieves correct energy equipartition, and a series of tumbler simulations that show the expected transition from tumbling to centrifuging as a function of rotation rate.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

An analytic parameterization of self-gravity wakes in Saturn's rings, with application to occultations and propellers

We have developed a semi-analytic method of parameterizing N-body simulations of self-gravity wakes in Saturn's rings, describing their properties by means of only 6 numbers: 3 optical depths and 3 weighting factors. These numbers are obtained using a density-estimation procedure that finds the frequencies of various values of local density within a simulated ring patch. Application of our parameterization to a suite of N-body simulations implies that the distribution of local optical depths is trimodal, rather than bimodal as previous authors have assumed. Rings dominated by self-gravity wakes appear to be mostly empty space. The implications of this result for the analysis of occultation data are more conceptual than practical. The only adjustment needed is that the model parameter $τ_{gap}$ should be interpreted as representing the area-weighted average optical depth within the gaps (or inter-wake regions). The most significant consequence of our results applies to the question of why "propeller" structures observed in the mid-A ring are seen as relative-bright features, even though the most prominent features of simulated propellers are regions of relatively low density. We find preliminary quantitative support for the hypothesis that propellers would be bright if they involve a local and temporary disruption of self-gravity wakes, flooding the region with more "photometrically active" material (i.e., material that can contribute to the rings' local optical depth) even though the overall density is lower, and thus raising their apparent brightnesses in agreement with observations. We suggest that this mechanism be tested by future detailed numerical models.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Fragment properties at the catastrophic disruption threshold: The effect of the parent body's internal structure

Numerical simulations of asteroid break-ups, including both the fragmentation of the parent body and the gravitational interactions between the fragments, have allowed us to reproduce successfully the main properties of asteroid families formed in different regimes of impact energy, starting from a non-porous parent body. In this paper, using the same approach, we concentrate on a single regime of impact energy, the so-called catastrophic threshold usually designated by Q*D, which results in the escape of half of the target's mass. Thanks to our recent implementation of a model of fragmentation of porous materials, we can characterize Q*D for both porous and non-porous targets with a wide range of diameters. We can then analyze the potential influence of porosity on the value of Q*D, and by computing the gravitational phase of the collision in the gravity regime, we can characterize the collisional outcome in terms of the fragment size and ejection speed distributions, which are the main outcome properties used by collisional models to study the evolutions of the different populations of small bodies. We also check the dependency of Q*D on the impact speed of the projectile. In the strength regime, which corresponds to target sizes below a few hundreds of meters, we find that porous targets are more difficult to disrupt than non-porous ones. In the gravity regime, the outcome is controlled purely by gravity and porosity in the case of porous targets. In the case of non-porous targets, the outcome also depends on strength. We then propose some power-law relationships between Q*D and both target's size and impact speed that can be used in collisional evolution models.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

N-Body Simulations of Growth from 1 km Planetesimals at 0.4 AU

We present N-body simulations of planetary accretion beginning with 1 km radius planetesimals in orbit about a 1 solar mass star at 0.4 AU. The initial disk of planetesimals contains too many bodies for any current N-body code to integrate; therefore, we model a sample patch of the disk. Although this greatly reduces the number of bodies, we still track in excess of 10^5 particles. We consider three initial velocity distributions and monitor the growth of the planetesimals. The masses of some particles increase by more than a factor of 100. Additionally, the escape speed of the largest particle grows considerably faster than the velocity dispersion of the particles, suggesting impending runaway growth, although no particle grows large enough to detach itself from the power law size-frequency distribution. These results are in general agreement with previous statistical and analytical results. We compute rotation rates by assuming conservation of angular momentum around the center of mass at impact and that merged planetesimals relax to spherical shapes. At the end of our simulations, the majority of bodies that have undergone at least one merger are rotating faster than the breakup frequency. This implies that the assumption of completely inelastic collisions (perfect accretion), which is made in most simulations of planetary growth at sizes 1 km and above, is inappropriate. Our simulations reveal that, subsequent to the number of particles in the patch having been decreased by mergers to half its initial value, the presence of larger bodies in neighboring regions of the disk may limit the validity of simulations employing the patch approximation.

0 citations0 reviewsNo code linkedOpen access