Related papers: An N-body approach to modeling debris and ejecta o…
We report on our study of asteroidal breakups, i.e. fragmentations of targets, subsequent gravitational reaccumulation and formation of small asteroid families. We focused on parent bodies with diameters $D_{\rm pb} = 10$ km. Simulations…
We present a novel way of modeling common envelope evolution in binary and few-body systems. We consider the common envelope inspiral as driven by a drag force with a power-law dependence in relative distance and velocity. The orbital…
We describe an algorithm for constructing N-body realisations of equilibrium stellar systems. The algorithm complements existing orbit-based modelling techniques using linear programming or other optimization algorithms. The equilibria are…
We introduce a numerical method to integrate tidal effects on collisional systems, using any definition of the external potential as a function of space and time. Rather than using a linearisation of the tidal field, this new method follows…
Since its first formulation in 1983, Milgromian dynamics (MOND) has been very successful in predicting the gravitational potential of galaxies from the distribution of baryons alone, including general scaling relations and detailed rotation…
The general consensus in the N-body community is that statistical results of an ensemble of collisional N-body simulations are accurate, even though individual simulations are not. A way to test this hypothesis is to make a direct…
We present a new C++ code for collisional N-body simulations of star clusters. The code uses the Hermite fourth-order scheme with block time steps, for advancing the particles in time, while the forces and neighboring particles are computed…
We introduce a GPU-accelerated hybrid hydro/N-body code (Enzo-N) designed to address the challenges of concurrently simulating star clusters and their parent galaxies. This task has been exceedingly challenging, primarily due to the…
We investigated numerically the dust growth driven by Brownian motion in a proto-planetary disc around a solar-type young stellar object. This process is considered as the first stage in the transformation of the initially micron-sized…
Recent exoplanet observations have revealed a diversity of exoplanetary systems, which suggests the ubiquity of radial planetary migration. One powerful known mechanism of planetary migration is planetesimal-driven migration (PDM), which…
The Double Asteroid Redirection Test (DART) is a NASA mission intended to crash a projectile on Dimorphos, the secondary component of the binary (65803) Didymos system, to study its orbit deflection. As a consequence of the impact, a dust…
The gravitational $N$-body problem is a nearly universal problem in astrophysics which, despite its deceptive simplicity, still presents a significant computational challenge. For collisional systems such as dense star clusters, the need to…
(Abridged) We describe a new iterative approach for the realization of equilibrium N-body systems for given density distributions. Our method uses elements of Schwarzschild's technique and of the made-to-measure method, but is based on a…
PlanetPack, initially released in 2013, is a command-line software aimed to facilitate exoplanets detection, characterization, and basic dynamical $N$-body simulations. This paper presents the third major release of PlanetPack that…
A model is presented of a ballistic, collisionless, steady state population of ejecta launched at randomly distributed times and velocities and moving under constant gravity above the surface of an airless planetary body. Within a…
Given that secular perturbations in a binary system not only excite high orbital eccentricities but also alter the planetary orbital inclination, the classical Keplerian orbital model is no longer applicable for orbital retrieval. The…
The applicability of the highly idealized secondary infall model to `realistic' initial conditions is investigated. The collapse of proto-halos seeded by $3\sigma$ density perturbations to an Einstein--de Sitter universe is studied here for…
We present the reduced basis method as a tool for developing emulators for equations with tunable parameters within the context of the nuclear many-body problem. The method uses a basis expansion informed by a set of solutions for a few…
We present a new N-body code, gevolution, for the evolution of large scale structure in the Universe. Our code is based on a weak field expansion of General Relativity and calculates all six metric degrees of freedom in Poisson gauge.…
We present the public release of EXP, a basis function expansion C++ library and Python package for running N-body galactic simulations and dynamical discovery. EXP grew out of the need for methodology that seamlessly connects theoretical…