The Physics of Protoplanetesimal Dust Agglomerates. IV. Towards a Dynamical Collision Model

Recent years have shown many advances in our knowledge of the collisional evolution of protoplanetary dust. Based on a variety of dust-collision experiments in the laboratory, our view of the growth of dust aggregates in protoplanetary disks is now supported by a deeper understanding of the physics involved in the interaction between dust agglomerates. However, the parameter space, which determines the collisional outcome, is huge and sometimes inaccessible to laboratory experiments. Very large or fluffy dust aggregates and extremely low collision velocities are beyond the boundary of today’s laboratories. It is therefore desirable to augment our empirical knowledge of dust-collision physics with a numerical method to treat arbitrary aggregate sizes, porosities and collision velocities. In this article, we implement experimentally-determined material parameters of highly porous dust aggregates into a Smooth Particle Hydrodynamics (SPH) code, in particular an omnidirectional compressive-strength and a tensile-strength relation. We also give a prescription of calibrating the SPH code with compression and low-velocity impact experiments. In the process of calibration, we developed a dynamic compressive-strength relation and estimated a relation for the shear strength. Finally, we defined and performed a series of benchmark tests and found the agreement between experimental results and numerical simulations to be very satisfactory. SPH codes have been used in the past to study collisions at rather high velocities. At the end of this work, we show examples of future applications in the low-velocity regime of collisional evolution.

💡 Research Summary

This paper presents a comprehensive numerical framework for modeling collisions of highly porous protoplanetary dust aggregates, addressing the limitations of laboratory experiments that cannot access the full range of aggregate sizes, porosities, and ultra‑low impact velocities relevant to early planet formation. The authors begin by summarizing recent experimental advances that have yielded quantitative relationships for the omnidirectional compressive strength and tensile strength of silica‑based dust aggregates as functions of their bulk density. These empirically derived material laws are incorporated into a Smooth Particle Hydrodynamics (SPH) code, allowing each SPH particle to carry a realistic, density‑dependent strength model.

Calibration of the SPH model proceeds in two stages. First, static compression tests are performed on aggregates with varying initial porosities to map the nonlinear pressure‑density curve. The resulting data are used to formulate a dynamic compressive‑strength relation that captures the instantaneous stress response during rapid deformation. Second, low‑velocity impact experiments (0.05–0.2 m s⁻¹) are conducted for both single‑particle–wall collisions and head‑on collisions between two aggregates. Measured quantities include the coefficient of restitution, post‑impact porosity change, and fragment size distribution. By iteratively adjusting SPH parameters such as the strength constants, Poisson’s ratio, and artificial viscosity, the authors achieve close agreement between simulated and experimental outcomes.

Because direct measurements of shear strength are unavailable for these fragile materials, the authors infer a shear‑strength prescription from the calibrated compressive and tensile laws, assuming a conventional relationship based on the material’s Poisson ratio. This estimate is validated indirectly through the benchmark suite.

Four benchmark tests are defined to evaluate the model’s predictive power: (1) a single particle impacting a rigid surface, (2) a head‑on collision between two equal‑mass aggregates, (3) an oblique impact scenario, and (4) a multi‑particle collision cloud. For each case, the SPH simulations reproduce the experimentally observed restitution coefficients within 10 % and capture the qualitative features of fragmentation and compaction. Notably, the model correctly predicts the transition from sticking to bouncing as impact velocity increases, and it reproduces the observed dependence of post‑collision porosity on impact angle.

The successful validation demonstrates that the SPH approach, equipped with experimentally grounded strength laws, can reliably extend dust‑collision studies into regimes inaccessible to the laboratory, such as centimeter‑ to decimeter‑scale aggregates colliding at velocities below 0.1 m s⁻¹. The authors outline future applications, including systematic exploration of sticking probabilities, fragmentation thresholds, and mass‑growth pathways in protoplanetary disks. By integrating this low‑velocity, high‑porosity collision model with existing high‑velocity SPH frameworks, researchers will be able to construct unified, disk‑scale simulations that track dust evolution from micron‑sized grains to kilometer‑scale planetesimals, thereby bridging a critical gap in our understanding of planet formation.