2025/11/13 by Tatsuuma, Misako, Okuzumi, Satoshi, Kataoka, Akimasa +1
Physics and Astronomy · #Astro and Planetary Science #Astronomical and nuclear sciences #Earth and Planetary Astrophysics (astro-ph.EP) #FOS: Physical sciences #Nuclear physics research studies
paper · doi:10.48550/arxiv.2511.09835
openalex publication_date 2025/11/13 · openalex created_date 2025/11/15 · openalex updated_date 2026/07/28
Modeling the contact surfaces formed by pebble collisions is crucial to understanding the formation process of comets, which are thought to be composed of pebbles. In this paper, we develop a new model to estimate the contact surface radius and the number of contact points as functions of collision velocity, and examine the formation process of comet 67P/Churyumov--Gerasimenko. Our model is based on the compressive strength of dust aggregates obtained from numerical simulations and assumes that all the impact energy of the pebbles is used for their mutual compression. We compare our model with numerical simulations of pebble collisions, in which we prepare the initial pebbles in the form of compressed dust aggregate spheres and measure the contact surface and pebble radii using two- and three-dimensional characteristic radii, respectively. We also apply our model to the formation scenario of comet 67P, whose tensile strength and bulk density have already been estimated in the literature. We find that its low tensile strength points to formation via pebble collisions at velocities below ∼10\mathrm cm s-1 when a microscopic filling factor of pebbles is lower than 0.6, suggesting that inelastic bouncing collisions played a role in damping the collision velocities. By assuming that the pebble collision velocity is determined by the transition velocity between bouncing and sticking, we estimate the pebble radius inside comet 67P to be 130 μm or smaller.