Related papers: Long-term stability of collisionless equilibrium c…
We investigate the connection between collisionless equilibria and the phase-space relation between density $\rho$ and velocity dispersion $\sigma$ found in simulations of dark matter halo formation, $F=\psd \propto r^{-\alpha}$.…
We developed a new direct-tree hybrid N-body algorithm for fully self-consistent N-body simulations of star clusters in their parent galaxies. In such simulations, star clusters need high accuracy, while galaxies need a fast scheme because…
The main problem is to understand and to find periodic symmetric orbits in the $n$-body problem, in the sense of finding methods to prove or compute their existence, and more importantly to describe their qualitative and quantitative…
We assemble and observe a sample of poor galaxy systems that is suitable for testing N-body simulations of hierarchical clustering (Navarro, Frenk, & White 1997; NFW) and other dynamical halo models (e.g., Hernquist 1990). We (1) determine…
Two-body scattering and other discreteness effects are unimportant in cosmological gravitational clustering in most scenarios, since the dark matter has a small particle mass. The collective field should determine evolution: Two-body…
In a Keplerian system, a large number of bodies orbit a central mass. Accretion disks, protoplanetary disks, asteroid belts, and planetary rings are examples. Simulations of these systems require algorithms that are computationally…
Numerical integrations of the Solar System have been carried out for decades. Their results have been used, for example, to determine whether the Solar System is chaotic, whether Mercury's orbit is stable, or to help discern Earth's climate…
Galaxy observations and N-body cosmological simulations produce conflicting dark matter halo density profiles for galaxy central regions. While simulations suggest a cuspy and universal density profile (UDP) of this region, the majority of…
We have implemented a parallel version of the Barnes-Hut 3-D N-body tree algorithm under PVM 3.2.5, adopting an SPMD paradigm. We parallelize the problem by decomposing the physical domain by means of the {\bf Orthogonal Recursive…
We present the first field-level comparison of cosmological N-body simulations, considering various widely used codes: Abacus, CUBEP$^3$M, Enzo, Gadget, Gizmo, PKDGrav, and Ramses. Unlike previous comparisons focused on summary statistics,…
This work discusses the main analogies and differences between the deterministic approach underlying most cosmological N-body simulations and the probabilistic interpretation of the problem that is often considered in mathematics and…
We present an N-body computer code - aimed at studies of galactic dynamics - with a CPU-efficient algorithm for a continuous (i.e. time-dependent) stellar mass-loss. First, we summarize available data on stellar mass-loss and derive the…
The gauge polyvalence of a new numerical code is tested, both in harmonic-coordinate simulations (gauge-waves testbed) and in singularity-avoiding coordinates (simple Black-Hole simulations, either with or without shift). The code is built…
The first attempt at developing a fully self-consistent code coupling dynamics and collisions to study debris discs (Kral, Thebault, Charnoz, 2013) is presented. So far, these two crucial mechanisms were studied separately, with N-body and…
Using the Ricci and scalar curvatures of the configuration manifold of gravitational N-body systems, we study the exponential instability in their trajectories. It is found that the exponentiation time-scale for isotropic Plummer spheres…
We present an efficient method for building equilibrium multi-component galaxies with non-spherical haloes and bulges. The gist of our approach is to tailor the velocity ellipsoid directly to the geometry of the mass distribution. Thus we…
We propose a new numerical procedure to simulate a single dark halo of any size and mass in a hierarchical framework coupling the extended Press-Schechter formalism (EPSF) to N-body simulations. The procedure consists of assigning…
In this work we present cosmological N-body simulations of the Local Universe with initial conditions constrained by the Two-Micron Redshift Survey (2MRS) within a cubic volume of 180 Mpc/h side-length centred at the Local Group. We use a…
Self-gravitating Newtonian systems consisting of a very large number of particles have generally defied attempts to describe them using statistical mechanics. This is paradoxical since many astronomical systems, or simulations thereof,…
This paper presents a study of the use of numerical simulation and Bayesian optimisation techniques to investigate the dynamics of celestial systems. Initially, the study focuses on Lagrange points in restricted three-body systems where a…