Related papers: An N-body approach to modeling debris and ejecta o…
Although planets have been found orbiting binary systems, whether they can survive binary interactions is debated. While the tightest-orbit binaries should host the most dynamically stable and long-lived circumbinary planetary systems, they…
We present a code allowing to evolve three-dimensional self-gravitating collisionless systems with a large number of particles N >= 10^7. FLY (Fast Level-based N-bodY code) is a fully parallel code based on a tree algorithm. It adopts…
A trending astronomical phenomenon to study is the variation in brightness of asteroids, caused by its rotation on its own axis, non-spherical shapes, changes of albedo along its surface and its position relative to the sun. The latter…
During close approaches to planets or stars, the morphological and dynamical properties of rubble-pile small bodies can be modified, and some may catastrophically break up. This phenomenon is of particular interest for understanding the…
When fitting N-body models to astronomical data - including transit times, radial velocity, and astrometric positions at observed times - the derivatives of the model outputs with respect to the initial conditions can help with model…
The development of numerical N -body simulations have allowed to study formation process and evolution of galaxies at different scales. This paper presents the fundamental concepts of N-body systems applied to the cosmological evolution of…
Early detection and tracking of ejecta in the vicinity of small solar system bodies is crucial to guarantee spacecraft safety and support scientific observation. During the visit of active asteroid Bennu, the OSIRIS-REx spacecraft relied on…
Previous research suggests that impacts between planetary embryos and planetesimals during the late stages of planet formation can often determine the percentages of core and mantle material that compose the newly formed planets in a…
We investigate the use of reduced-order modelling to run discrete element simulations at higher speeds. Taking a data-driven approach, we run many offline simulations in advance and train a model to predict the velocity field from the mass…
The "gravitational million-body problem," to model the dynamical evolution of a self-gravitating, collisional N-body system with ~10^6 particles over many relaxation times, remains a major challenge in computational astrophysics.…
We describe the numerical code N-MODY, a parallel particle-mesh code for collisionless N-body simulations in modified Newtonian dynamics (MOND). N-MODY is based on a numerical potential solver in spherical coordinates that solves the…
We present the results of planet formation N-body simulations based on a comprehensive physical model that includes planetary mass growth through mutual embryo collisions and planetesimal/boulder accretion, viscous disc evolution, planetary…
Black hole X-ray binaries in outburst launch discrete, large-scale jet ejections which can propagate to parsec scales. The kinematics of these ejecta appear to be well described by relativistic blast wave models original devised for…
We present a new, simple, fast algorithm to numerically evolve disks of inelastically colliding particles surrounding a central star. Our algorithm adds negligible computational cost to the fastest existing collisionless N-body codes, and…
We present a novel method for including the impact of massive neutrinos in cold dark matter N-body simulations. Our approach is compatible with widely employed Newtonian N-body codes and relies on only three simple modifications. First, we…
In the standard model of terrestrial planet formation, planets are formed through giant impacts of planetary embryos after the dispersal of the protoplanetary gas disc. Traditionally, $N$-body simulations have been used to investigate this…
The decay process of the schematic one-dimensional three-body system is considered. A time-dependent approach is used in combination with a one-dimensional three-body model, which is composed of a heavier core nucleus and two nucleons, with…
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…
The dynamics of planetesimals plays an important role in planet formation, because their velocity distribution sets the growth rate to larger bodies. When planetesimals form in protoplanetary discs, their orbits are nearly circular and…
To more thoroughly study the effects of radiative forces on the orbits of small, astronomical bodies, we introduce the Yarkovsky effect into REBOUNDx, an extensional library for the N-body integrator REBOUND. Two different versions of the…