Paper
Exocomets of $β$ Pictoris II: Two dynamical families of exocomets simulated with REBOUND
Authors
K. P. Jaworska, H. J. Hoeijmakers
Abstract
We investigate the dynamical evolution of particles in the $β$ Pic system to determine likely formation pathways to the present-day observed exocomet populations. We aim to relate these results to similar studies recently carried out since the discovery of the inner planet $β$ Pic c. We simulate the $β$ Pic system using the non-symplectic adaptive N-body integrator IAS15 in REBOUND. We seed the system with over 100,000 mass-less test particles that evolve for 25 Myr, and adopt initial conditions and a particle distribution that closely matches similar simulations in recent literature. Using IAS15, REBOUND resolves close-encounters between test particles and the two gas giants in the system, which is crucial for understanding aspects of the dynamical evolution. Planet-disk interactions rapidly clear most of the system within 35 AU apart from a region within the orbit of $β$ Pic c, and a region between 20 and 25 AU. After 10 Myr, exocomets can be sourced continuously from these regions, as well as from the inner edge of the region beyond ~35 AU where particles are stable on longer timescales. From the region interior to $β$ Pic c, the exocomets are formed by excitation via mean-motion resonance with $β$ Pic c, obtaining a narrow distribution of radial velocities, consistent with spectroscopic observations. Particles initialized in the outer system may enter onto stargrazing orbits due to disruption by the two gas giants, causing a wider radial velocity distribution, and we propose that this population corresponds to a second dynamical family previously observed via spectroscopy. These particles typically undergo chaotic dynamical evolution for $10^2$ to $10^3$ years after passing the water sublimation limit at ~8 AU until reaching the sublimation distance of calcium near 0.4 AU, implying that the two families of exocomets may have different volatile contents.
Metadata
Related papers
Cosmic Shear in Effective Field Theory at Two-Loop Order: Revisiting $S_8$ in Dark Energy Survey Data
Shi-Fan Chen, Joseph DeRose, Mikhail M. Ivanov, Oliver H. E. Philcox • 2026-03-30
Stop Probing, Start Coding: Why Linear Probes and Sparse Autoencoders Fail at Compositional Generalisation
Vitória Barin Pacela, Shruti Joshi, Isabela Camacho, Simon Lacoste-Julien, Da... • 2026-03-30
SNID-SAGE: A Modern Framework for Interactive Supernova Classification and Spectral Analysis
Fiorenzo Stoppa, Stephen J. Smartt • 2026-03-30
Acoustic-to-articulatory Inversion of the Complete Vocal Tract from RT-MRI with Various Audio Embeddings and Dataset Sizes
Sofiane Azzouz, Pierre-André Vuissoz, Yves Laprie • 2026-03-30
Rotating black hole shadows in metric-affine bumblebee gravity
Jose R. Nascimento, Ana R. M. Oliveira, Albert Yu. Petrov, Paulo J. Porfírio,... • 2026-03-30
Raw Data (Debug)
{
"raw_xml": "<entry>\n <id>http://arxiv.org/abs/2603.05600v1</id>\n <title>Exocomets of $β$ Pictoris II: Two dynamical families of exocomets simulated with REBOUND</title>\n <updated>2026-03-05T19:00:47Z</updated>\n <link href='https://arxiv.org/abs/2603.05600v1' rel='alternate' type='text/html'/>\n <link href='https://arxiv.org/pdf/2603.05600v1' rel='related' title='pdf' type='application/pdf'/>\n <summary>We investigate the dynamical evolution of particles in the $β$ Pic system to determine likely formation pathways to the present-day observed exocomet populations. We aim to relate these results to similar studies recently carried out since the discovery of the inner planet $β$ Pic c. We simulate the $β$ Pic system using the non-symplectic adaptive N-body integrator IAS15 in REBOUND. We seed the system with over 100,000 mass-less test particles that evolve for 25 Myr, and adopt initial conditions and a particle distribution that closely matches similar simulations in recent literature. Using IAS15, REBOUND resolves close-encounters between test particles and the two gas giants in the system, which is crucial for understanding aspects of the dynamical evolution. Planet-disk interactions rapidly clear most of the system within 35 AU apart from a region within the orbit of $β$ Pic c, and a region between 20 and 25 AU. After 10 Myr, exocomets can be sourced continuously from these regions, as well as from the inner edge of the region beyond ~35 AU where particles are stable on longer timescales. From the region interior to $β$ Pic c, the exocomets are formed by excitation via mean-motion resonance with $β$ Pic c, obtaining a narrow distribution of radial velocities, consistent with spectroscopic observations. Particles initialized in the outer system may enter onto stargrazing orbits due to disruption by the two gas giants, causing a wider radial velocity distribution, and we propose that this population corresponds to a second dynamical family previously observed via spectroscopy. These particles typically undergo chaotic dynamical evolution for $10^2$ to $10^3$ years after passing the water sublimation limit at ~8 AU until reaching the sublimation distance of calcium near 0.4 AU, implying that the two families of exocomets may have different volatile contents.</summary>\n <category scheme='http://arxiv.org/schemas/atom' term='astro-ph.EP'/>\n <published>2026-03-05T19:00:47Z</published>\n <arxiv:comment>12 pages, 10 figures, accepted for publication in A&A</arxiv:comment>\n <arxiv:primary_category term='astro-ph.EP'/>\n <author>\n <name>K. P. Jaworska</name>\n </author>\n <author>\n <name>H. J. Hoeijmakers</name>\n </author>\n </entry>"
}