Related papers: Should $N$-body integrators be symplectic everywhe…
Since they were introduced in the 1990s, Lie group integrators have become a method of choice in many application areas. These include multibody dynamics, shape analysis, data science, image registration and biophysical simulations. Two…
In a previous paper, second- and fourth-order explicit symplectic integrators were designed for a Hamiltonian of the Schwarzschild black hole. Following this work, we continue to trace the possibility of the construction of explicit…
Numerical simulations of self-gravitating systems are generally based on N-body codes, which solve the equations of motion of a large number of interacting particles. This approach suffers from poor statistical sampling in regions of low…
The conservation of energy, linear momentum and angular momentum are important drivers for our physical understanding of the evolution of the Universe. These quantities are also conserved in Newton's laws of motion under gravity…
We describe an algorithm for long-term planetary orbit integrations, including the dominant post-Newtonian effects, that employs individual timesteps for each planet. The algorithm is symplectic and exhibits short-term errors that are…
We implement several symplectic integrators, which are based on two part splitting, for studying the chaotic behavior of one- and two-dimensional disordered Klein-Gordon lattices with many degrees of freedom and investigate their numerical…
We introduce a novel numerical method to integrate partial differential equations representing the Hamiltonian dynamics of field theories. It is a multi-symplectic integrator that locally conserves the stress-energy tensor with an excellent…
Accurate $N$-body simulations of multiple systems such as binaries and triples are essential for understanding the formation and evolution of interacting binaries and binary mergers, including gravitational wave sources, blue stragglers and…
Astrophysical Challenges which demand the solution of the one million (or more) gravitating body problem are briefly discussed for the fields of cosmology, galactic nuclei and globular star clusters. Results from the classical three-body…
In this paper we describe a new algorithm for the long-term numerical integration of the two-body problem, in which two particles interact under a Newtonian gravitational potential. Although analytical solutions exist in the unperturbed and…
In this paper we study the performance of a symplectic numerical integrator based on the splitting method. This method is applied to a subtle problem i.e. higher order resonance of the elastic pendulum. In order to numerically study the…
A wide variety of outstanding problems in astrophysics involve the motion of a large number of particles ($N\gtrsim 10^{6}$) under the force of gravity. These include the global evolution of globular clusters, tidal disruptions of stars by…
Two documents have recently been posted on the arXiv describing a numerical integration algorithm: "symplectic orbit/spin tracking code for all-electric storage rings" [1] and some computational results therefrom [2]. This note comments…
We provide topological obstructions to the existence of orbit cylinders of symmetric orbits, for mechanical systems preserved by antisymplectic involutions (e.g. the restricted three-body problem). Such cylinders induce continuous paths…
Periodic and quasi-periodic orbits of the $n$-body problem are critical points of the action functional constrained to the Sobolev space of symmetric loops. Variational methods yield collisionless orbits provided the group of symmetries…
We review the implementation of individual particle time-stepping for N-body dynamics. We present a class of integrators derived from second order Hamiltonian splitting. In contrast to the usual implementation of individual time-stepping,…
An energetically balanced, implicit integrator for non-hydrostatic vertical atmospheric dynamics on the sphere is presented. The integrator allows for the exact balance of energy exchanges in space and time for vertical atmospheric motions…
This article is a survey of classical and quantum completely integrable systems from the viewpoint of local ``phase space'' analysis. It advocates the use of normal forms and shows how to get global information from glueing local pieces.…
We consider the numerical integration of the matrix Hill's equation. Parametric resonances can appear and this property is of great interest in many different physical applications. Usually, the Hill's equations originate from a Hamiltonian…
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…