Errors, chaos and the collisionless limit

Abstract

We simultaneously study the dynamics of the growth of errors and the question of the faithfulness of simulations of $N$-body systems. The errors are quantified through the numerical reversibility of small-$N$ spherical systems, and by comparing fixed-timestep runs with different stepsizes. The errors add randomly, before exponential divergence sets in, with exponentiation rate virtually independent of $N$, but scale saturating as $\sim 1/\sqrt{N}$, in line with theoretical estimates presented. In a third phase, the growth rate is initially driven by multiplicative enhancement of errors, as in the exponential stage. It is then qualitatively different for the phase space variables and mean field conserved quantities (energy and momentum); for the former, the errors grow systematically through phase mixing, for the latter they grow diffusively. For energy, the $N$-variation of the `relaxation time' of error growth follows the $N$-scaling of two-body relaxation. This is also true for angular momentum in the fixed stepsize runs, although the associated error threshold is higher and the relaxation time smaller. Due to shrinking saturation scales, the information loss associated with the exponential instability decreases with $N$ and the dynamical entropy vanishes at any finite resolution as $N\to \infty$. A distribution function depending on the integrals of motion in the smooth potential is decreasingly affected. In this sense there is convergence to the collisionless limit, despite the persistence of exponential instability on infinitesimal scales. Nevertheless, the slow $N$-variation in its saturation points to the slowness of the convergence.

arXiv commentRevised version, to appear in MNRAS. Error growth with symplectic leapfrog algorithm examined and compared with fixed stepsize standard Runge Kutta (Section 3.5) and Section 4 largely rewritten for clarity. Main conclusions unchanged

References in corpus (17)

• The Gravitational Million-Body Problem
• Statistics of N-body Simulations I. Equal Masess Before Core Collapse
• N-body simulations of gravitational dynamics
• A Quantification of the Butterfly Effect in Cosmological Simulations and Implications for Galaxy Scaling Relations
• Chaos and Variance in Galaxy Formation
• Stochasticity in N-body Simulations of Disc Galaxies
• Instability of the Gravitational N-Body Problem in the Large-N Limit
• Chaos and the continuum limit in the gravitational N-body problem. I. Integrable potentials
• Smooth potential chaos and N-body simulations
• Chaos and the continuum limit in the gravitational N-body problem II. Nonintegrable potentials
• Time-symmetric integration in astrophysics
• Towards time symmetric N-body integration
• On the Onset of Stochasticity in $\Lambda$CDM Cosmological Simulations
• Shadowing high-dimensional Hamiltonian systems: the gravitational n-body problem
• Regular and chaotic motion in softened gravitational systems
• Particle number dependence in the non-linear evolution of N-body self-gravitating systems
• Symplectic Integration without Roundoff Error

Cited by in corpus (8)

• Large-scale dark matter simulations
• Numerical convergence of hydrodynamical simulations of galaxy formation: the abundance and internal structure of galaxies and their cold dark matter haloes
• $N$-body chaos and the continuum limit in numerical simulations of self-gravitating systems, revisited
• Discreteness effects, $N-$body chaos and the onset of radial-orbit instability
• Symplectic coarse graining approach to the dynamics of spherical self-gravitating systems
• Partial suppression of chaos in relativistic three-body problems
• A disturbance in the force. How force fluctuations hinder dynamical friction and induce core stalling
• Chaos in violent relaxation dynamics. Disentangling micro- and macro-chaos in numerical experiments of dissipationless collapse