Interactions Among Non-Interacting Particles in Planet Formation Simulations
arXiv:2007.11758 · doi:10.3847/2041-8213/aba68d
Abstract
Over the course of the recent decades, -body simulations have become a standard tool for quantifying the gravitational perturbations that ensue in planet-forming disks. Within the context of such simulations, massive non-central bodies are routinely classified into "big" and "small" particles, where big objects interact with all other objects self-consistently, while small bodies interact with big bodies but not with each other. Importantly, this grouping translates to an approximation scheme where the orbital evolution of small bodies is dictated entirely by the dynamics of the big bodies, yielding considerable computational advantages with little added cost in terms of astrophysical accuracy. Here we point out, however, that this scheme can also yield spurious dynamical behaviour, where even in absence of big bodies within a simulation, indirect coupling among small bodies can lead to excitation of the constituent "non-interacting" orbits. We demonstrate this self-stirring by carrying out a sequence of numerical experiments, and confirm that this effect is largely independent of the time-step or the employed integration algorithm. Furthermore, adopting the growth of angular momentum deficit as a proxy for dynamical excitation, we explore its dependence on time, the cumulative mass of the system, as well as the total number of particles present in the simulation. Finally, we examine the degree of such indirect excitation within the context of conventional terrestrial planet formation calculations, and conclude that although some level of caution may be warranted, this effect plays a negligible role in driving the simulated dynamical evolution.
References in corpus (15)
- A low mass for Mars from Jupiter's early gas-driven migration
- REBOUND: An open-source multi-purpose N-body code for collisional dynamics
- IAS15: A fast, adaptive, high-order integrator for gravitational dynamics, accurate to machine precision over a billion orbits
- Building Terrestrial Planets
- Statistical Study of the Early Solar System's Instability with 4, 5 and 6 Giant Planets
- Formation of the Terrestrial Planets from a Narrow Annulus
- Hybrid Symplectic Integrators for Planetary Dynamics
- Neptune's Orbital Migration Was Grainy, Not Smooth
- Lunar and Terrestrial Planet Formation in the Grand Tack Scenario
- Formation of Giant Planet Satellites
- A new inclination instability reshapes Keplerian disks into cones: application to the outer Solar System
- Dynamics of escaping Earth ejecta and their collision probability with different Solar System bodies
- Chaotic Excitation and Tidal Damping in the GJ 876 System
- JANUS: A bit-wise reversible integrator for N-body dynamics
- Appearance of Saturn's F ring azimuthal channels for the anti-alignment configuration between the ring and Prometheus