Related papers: Improving the accuracy of simulated chaotic $N$-bo…
A new approach is developed to integrate numerically the equations of motion for systems of interacting rigid polyatomic molecules. With the aid of a leapfrog framework, we directly involve principal angular velocities into the integration,…
We derive variational integrators for stochastic Hamiltonian systems on Lie groups using a discrete version of the stochastic Hamiltonian phase space principle. The structure-preserving properties of the resulting scheme, such as…
We study how inexact nonlinear solvers lead to a loss of exact symplecticity in the Symplectic Euler (SE) and Stormer-Verlet (SV) schemes when applied to general nonseparable Hamiltonian systems. These schemes are implicit and require…
We present a simple algorithm to switch between $N$-body time integrators in a reversible way. We apply it to planetary systems undergoing arbitrarily close encounters and highly eccentric orbits, but the potential applications are broader.…
In this work, the benefits of the phase fitting technique are embedded in high order discrete Lagrangian integrators. The proposed methodology creates integrators with zero phase lag in a test Lagrangian in a similar way used in phase…
Symplectic integrators with long-term preservation of integrals of motion are introduced for the guiding-center model of plasma particles in toroidal magnetic fields of general topology. An efficient transformation to canonical coordinates…
Modified Hamiltonians are used in the field of geometric numerical integration to show that symplectic schemes for Hamiltonian systems are accurate over long times. For nonlinear systems the series defining the modified Hamiltonian usually…
Simulating the evolution of the gravitational N-body problem becomes extremely computationally expensive as N increases since the problem complexity scales quadratically with the number of bodies. We study the use of Artificial Neural…
We apply one of the exactly symplectic integrators, that we call HB15, of \cite{HB15}, along with the Kepler problem solver of \cite{WH15}, to solve planetary system $N$-body problems. We compare the method to Wisdom-Holman methods (WH) in…
We present a simple choice of integration variables that can be used to exploit the near-integrable character of problems in celestial mechanics. The approach is based on the well-known principle of variation of parameters: instead of…
Nonadiabatic behavior of metastable systems modeled by anharmonic Hamiltonians is reproduced by the Fokker-Planck and imaginary time Schrodinger equation scheme with subsequent symplectic integration. Example solutions capture ergodicity…
Symplectic integrators evolve dynamical systems according to modified Hamiltonians whose error terms are also well-defined Hamiltonians. The error of the algorithm is the sum of each error Hamiltonian's perturbation on the exact solution.…
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…
Chaotic hyperbolic dynamical systems enjoy a surprising degree of rigidity, a fact which is well known in the mathematics community but perhaps less so in theoretical physics circles. Low-dimensional hyperbolic systems are either conjugate…
In general relativity, the motion of an extended body moving in a given spacetime can be described by a particle on a (generally non-geodesic) worldline. In first approximation, this worldline is a geodesic of the underlying spacetime, and…
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…
We split the generic conformal mechanical system into a "radial" and an "angular" part, where the latter is defined as the Hamiltonian system on the orbit of the conformal group, with the Casimir function in the role of the Hamiltonian. We…
Harmonic inversion has already been proven to be a powerful tool for the analysis of quantum spectra and the periodic orbit orbit quantization of chaotic systems. The harmonic inversion technique circumvents the convergence problems of the…
We present several methods, which utilize symplectic integration techniques based on two and three part operator splitting, for numerically solving the equations of motion of the disordered, discrete nonlinear Schr\"odinger (DDNLS)…
Symplectic integrators offer many advantages for the numerical solution of Hamiltonian differential equations, including bounded energy error and the preservation of invariant sets. Two of the central Hamiltonian systems encountered in…