Related papers: JANUS: A bit-wise reversible integrator for N-body…
We present a non-canonically symplectic integration scheme tailored to numerically computing the post-Newtonian motion of a spinning black-hole binary. Using a splitting approach we combine the flows of orbital and spin contributions. In…
REBOUND is a new multi-purpose N-body code which is freely available under an open-source license. It was designed for collisional dynamics such as planetary rings but can also solve the classical N-body problem. It is highly modular and…
We study analytically and experimentally certain symplectic and time-reversible N-body integrators which employ a Kepler solver for each pair-wise interaction, including the method of Hernandez & Bertschinger (2015). Owing to the Kepler…
Simulating physically plausible trajectories toward user-defined goals is a fundamental yet challenging task in fluid dynamics. While particle-based simulators can efficiently reproduce forward dynamics, inverse inference remains difficult,…
We investigate the computational performance of various numerical methods for the integration of the equations of motion and the variational equations for some typical classical many-body models of condensed matter physics: the…
We prove that the dynamical system charaterized by the Hamiltonian $ H = \lambda N \sum_{j}^{N} p_j + \mu \sum_{j,k}^{N} {{(p_j p_k)}^{1\over 2}} \{ cos [ \nu ( q_j - q_k)] \} $ proposed and studied by Calogero [1,2] is equivalent to a…
In this work we introduce a geometric integrator for molecular dynamics simulations of physical systems in the canonical ensemble. In particular, we consider the equations arising from the so-called density dynamics algorithm with any…
Reversible algorithms are algorithms in which each step represents a partial injective function; they are useful for performance optimization in reversible systems. In this study, using Janus, a reversible imperative high-level programming…
We present WHFast, a fast and accurate implementation of a Wisdom-Holman symplectic integrator for long-term orbit integrations of planetary systems. WHFast is significantly faster and conserves energy better than all other Wisdom-Holman…
It is well known that the Hamiltonian of an $n$-dimensional isotropic oscillator admits an $SU(n)$ symmetry, making the system maximally superintegrable. However, the dynamical symmetries of the anisotropic oscillator are much more subtle.…
A special purpose solver, based on the Magnus expansion, well suited for the integration of the linear three neutrino oscillations equations in matter is proposed. The computations are speeded up to two orders of magnitude with respect to a…
Most direct N-body integrations of planetary systems use a symplectic integrator with a fixed timestep. A large timestep is desirable in order to speed up the numerical simulations. However, simulations yield unphysical results if the…
We consider geometric numerical integration algorithms for differential equations evolving on symmetric spaces. The integrators are constructed from canonical operations on the symmetric space, its Lie triple system (LTS), and the…
In early Solar System numerical simulations, where chaos is a primary driver, it is difficult to explore parameter space in a systematic way. In such simulations, stable configurations are hard to come by, and often require special…
Phase Space is the framework best suited for quantizing superintegrable systems, naturally preserving the symmetry algebras of the respective hamiltonian invariants. The power and simplicity of the method is fully illustrated through new…
This paper deals with an implicit Newton-like inertial dynamical system governed by a maximally comonotone inclusion problem in a Hilbert space. Under suitable conditions, we establish not only pointwise estimates and integral estimates for…
We consider a numerical scheme for Hamilton-Jacobi equations based on a direct discretization of the Lax-Oleinik semi-group. We prove that this method is convergent with respect to the time and space stepsizes provided the solution is…
We have introduced self-consistent spin, tidal and dynamical equations of motion into REBOUNDx, a library of additional effects for the popular N-body integrator REBOUND. The equations of motion used are derived from the constant time lag…
We study the influence of relativity on the chaotic properties and dynamical outcomes of an unstable triple system; the Pythagorean three-body problem. To this end, we extend the Brutus N-body code to include Post-Newtonian pairwise terms…
We present a numerical weak-lensing analysis that is fully relativistic and non-perturbative for the scalar part of the gravitational potential and first-order in the vector part, frame dragging. Integrating the photon geodesics backwards…