Related papers: High order symplectic integrators for perturbed Ha…
It is known that the standard and the inverted harmonic oscillator are different. Replacing thus of {\omega} by i{\omega} in the regular oscillator is necessary going to give the inverted oscillator H^{r}. This replacement would lead to…
This article considers non-relativistic charged particle dynamics in both static and non-static electromagnetic fields, which are governed by nonseparable, possibly time-dependent Hamiltonians. For the first time, explicit symplectic…
Hamilton's equations of motion form a fundamental framework in various branches of physics, including astronomy, quantum mechanics, particle physics, and climate science. Classical numerical solvers are typically employed to compute the…
Recently, a whole class of evergy-preserving integrators has been derived for the numerical solution of Hamiltonian problems. In the mainstream of this research, we have defined a new family of symplectic integrators depending on a real…
The superintegrability of four Hamiltonians $\tilde{H_r} = \lambda\, H_r$, $r=a,b,c,d$, where $H_r$ are known Hamiltonians and $\lambda$ is a certain function defined on the configuration space and depending of a parameter $\kappa$, is…
This article reviews some integrators particularly suitable for the numerical resolution of differential equations on a large time interval. Symplectic integrators are presented. Their stability on exponentially large time is shown through…
We consider a Hamiltonian system on the symplectic space $({\mathbb{R}}^{2n}, dy\wedge dx)$ with a real-analytic Hamiltonian $H : {\mathbb{R}}^{2n}\to {\mathbb{R}}$. We assume that the system has a non-degenerate equilibrium position at the…
We follow up on our previous works which presented a possible approach for deriving symplectic schemes for a certain class of highly oscillatory Hamiltonian systems. The approach considers the Hamilton-Jacobi form of the equations of…
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…
An accurate and efficient method dealing with the few-body dynamics is important for simulating collisional N-body systems like star clusters and to follow the formation and evolution of compact binaries. We describe such a method which…
Most numerical integration algorithms are not designed specifically for Hamiltonian systems and do not respect their characteristic properties, which include the preservation of phase space volume with time. This can lead to spurious…
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 methods that preserve geometric invariants of the system, such as energy, momentum or the symplectic form, are called geometric integrators. In this paper we present a method to construct symplectic-momentum integrators for…
We develop Lie-Poisson integrators for general Hamiltonian systems on $\mathbf{R}^{3}$ equipped with the rigid body bracket. The method uses symplectic realisation of $\mathbf{R}^{3}$ on $T^{*}\mathbf{R}^{2}$ and application of symplectic…
The dynamic equation of mass point in rotating coordinates is governed by Coriolis and centrifugal force, besides a corotating potential relative to frame. Such a system is no longer a canonical Hamiltonian system so that the construction…
Long-term stability studies of nonlinear Hamiltonian systems require symplectic integration algorithms which are both fast and accurate. In this paper, we study a symplectic integration method wherein the symplectic map representing the…
A variational and perturbative treatment is provided for a family of generalized spiked harmonic oscillator Hamiltonians H = -(d/dx)^2 + B x^2 + A/x^2 + lambda/x^alpha, where B > 0, A >= 0, and alpha and lambda denote two real positive…
We present a family of multistep integrators based on the Adams-Bashforth methods. These schemes can be constructed for arbitrary convergence order with arbitrary step size variation. The step size can differ between different subdomains of…
In this paper, we are concerned with the construction and analysis of a new class of methods obtained as double jump compositions with complex coefficients and projection on the real axis. It is shown in particular that the new integrators…
In a recent work of Wu, Wang, Sun and Liu, a second-order explicit symplectic integrator was proposed for the integrable Kerr spacetime geometry. It is still suited for simulating the nonintegrable dynamics of charged particles moving…