Related papers: Symplectic fourth-order maps for the collisional N…
It is well known that symplectic integrators lose their near energy preservation properties when variable step sizes are used. The most common approach to combine adaptive step sizes and symplectic integrators involves the Poincar\'e…
The two-dimensional n-body problem of classical mechanics is a non-integrable Hamiltonian system for n > 2. Traditional numerical integration algorithms, which are polynomials in the time step, typically lead to systematic drifts in the…
We present a method to calculate many-body states of interacting carriers in million atom quantum nanostructures based on atomistic tight-binding calculations and a combination of iterative selection of configurations and perturbation…
Symplectic integration methods based on operator splitting are well established in many branches of science. For Hamiltonian systems which split in more than two parts, symplectic methods of higher order have been studied in detail only for…
Explicit symplectic integrators have been important tools for accurate and efficient approximations of mechanical systems with separable Hamiltonians. For the first time, the article proposes for arbitrary Hamiltonians similar integrators,…
Cosmological N-body simulations are done on massively parallel computers. This necessitates the use of simple time integrators, and, additionally, of mesh-grid approximations of the potentials. Recently, Adamek et al. (2015);…
It has previously been shown that varying the numerical timestep during a symplectic orbital integration leads to a random walk in energy and angular momentum, destroying the phase space-conserving property of symplectic integrators. Here…
Efficient fourth order symplectic integrators are proposed for numerical integration of separable Hamiltonian systems H(p,q)=T(p)+V(q). Symmetric splitting coefficients with five to nine stages are obtained by higher order decomposition of…
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 propose implicit integrators for solving stiff differential equations on unit spheres. Our approach extends the standard backward Euler and Crank-Nicolson methods in Cartesian space by incorporating the geometric constraint inherent to…
We study systems of few two-component fermions interacting via short-range interactions within a harmonic-oscillator trap. The dominant interactions, which are two-body, are organized according to the number of derivatives and defined in a…
We prove the existence of periodic solutions of the N=(n+1)-body problem starting with n bodies whose reduced motion is close to a non-degenerate central configuration and replacing one of them by the center of mass of a pair of bodies…
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…
We study the structure of solutions to linear programming formulations for the traveling salesperson problem (TSP). We perform a detailed analysis of the support of the subtour elimination linear programming relaxation, which leads to…
We construct an explicit reversible symplectic integrator for the planar 3-body problem with zero angular momentum. We start with a Hamiltonian of the planar 3-body problem that is globally regularised and fully symmetry reduced. This…
Elegant integration schemes of second and fourth order for simulations of rigid body systems are presented which treat translational and rotational motion on the same footing. This is made possible by a recent implementation of the exact…
Whether rigid or compliant, contact interactions are inherent to robot motions, enabling them to move or manipulate things. Contact interactions result from complex physical phenomena, that can be mathematically cast as Nonlinear…
We present a class of symplectic integrators adapted for the integration of perturbed Hamiltonian systems of the form $H=A+\epsilon B$. We give a constructive proof that for all integer $p$, there exists an integrator with positive steps…
A new family of methods involving complex coefficients for the numerical integration of differential equations is presented and analyzed. They are constructed as linear combinations of symmetric-conjugate compositions obtained from a basic…
Symplectic integrators are widely used for long-term integration of conservative astrophysical problems due to their ability to preserve the constants of motion; however, they cannot in general be applied in the presence of nonconservative…