Related papers: Variational Integrators for the Gravitational N-Bo…
Symplectic integrators have made it possible to study the long-term evolution of planetary systems with direct N-body simulations. In this paper we reassess the accuracy of such simulations by running a convergence test on 20Myr…
We use a systolic N-body algorithm to evaluate the linear stability of the gravitational N-body problem for N up to 1.3 x 10^5, two orders of magnitude greater than in previous experiments. For the first time, a clear ~ln N-dependence of…
This work proposes a suite of numerical techniques to facilitate the design of structure-preserving integrators for nonlinear dynamics. The celebrated LaBudde-Greenspan integrator and various energy-momentum schemes adopt a difference…
Symplectic methods, in particular the Wisdom-Holman map, have revolutionized our ability to model the long-term, conservative dynamics of planetary systems. However, many astrophysically important effects are dissipative. The consequences…
Direct-summation N-body algorithms compute the gravitational interaction between stars in an exact way and have a computational complexity of O(N^2). Performance can be greatly enhanced via the use of special-purpose accelerator boards like…
We revisit the problem of the equations of motion of a system of $N$ self-interacting massive particles (without spins) in the first post-Minkowskian (1PM) approximation of general relativity. We write the equations of motion, gravitational…
Large-scale cosmological simulations are an indispensable tool for modern cosmology. To enable model-space exploration, fast and accurate predictions are critical. In this paper, we show that the performance of such simulations can be…
Most physical systems are modelled by an ordinary or a partial differential equation, like the n-body problem in celestial mechanics. In some cases, for example when studying the long term behaviour of the solar system or for complex…
A new methodology is developed to integrate numerically the equations of motion for classical many-body systems in molecular dynamics simulations. Its distinguishable feature is the possibility to preserve, independently on the size of the…
Symplectic integrators separate a problem into parts that can be solved in isolation, alternately advancing these sub-problems to approximate the evolution of the complete system. Problems with a single, dominant mass can use mixed-variable…
We show that adaptive time stepping in particle accelerator simulation is an enhancement for certain problems. The new algorithm has been implemented in the OPAL (Object Oriented Parallel Accelerator Library) framework, and is compared to…
We present a parallelized primal-dual algorithm for solving constrained convex optimization problems. The algorithm is "block-based," in that vectors of primal and dual variables are partitioned into blocks, each of which is updated only by…
We present a new mixed finite element method for a class of parabolic equations with $p$-Laplacian and nonlinear memory. The applicability, stability and convergence of the method are studied. First, the problem is written in a mixed…
We present and discuss the characteristics and performances, both in term of computational speed and precision, of a numerical code which numerically integrates the equation of motions of N 'particles' interacting via Newtonian gravitation…
We discuss an elementary derivation of variational symmetries and corresponding integrals of motion for the Lagrangian systems depending on acceleration. Providing several examples, we make the manuscript accessible to a wide range of…
We suggest a numerical integration procedure for solving the equations of motion of certain classical spin systems which preserves the underlying symplectic structure of the phase space. Such symplectic integrators have been successfully…
Modeling self-gravity of collisionless fluids (e.g. ensembles of dark matter, stars, black holes, dust, planetary bodies) in simulations is challenging and requires some force softening. It is often desirable to allow softenings to evolve…
We develop a formalism for General Relativistic N-body simulations in the weak field regime, suitable for cosmological applications. The problem is kept tractable by retaining the metric perturbations to first order, the first derivatives…
Numerical solutions to Newton's equations of motion for chaotic self gravitating systems of more than 2 bodies are often regarded to be irreversible. This is due to the exponential growth of errors introduced by the integration scheme and…
To describe the ``slow'' motions of n interacting mass points, we give the most general 4-d non-instantaneous, non-particle symmetric Galilei-invariant variational principle. It involves two-body invariants constructed from particle…