Related papers: Step-size effect in the time-transformed leapfrog …
In this paper, an implicit nonsymplectic exact energy-preserving integrator is specifically designed for a ten-dimensional phase-space conservative Hamiltonian system with five degrees of freedom. It is based on a suitable…
Symplectic integrators are the preferred method of solving conservative $N$-body problems in cosmological, stellar cluster, and planetary system simulations because of their superior error properties and ability to compute orbital…
Symplectic integrators are a foundation to the study of dynamical $N$-body phenomena, at scales ranging from from planetary to cosmological. These integrators preserve the Poincar\'e invariants of Hamiltonian dynamics. The $N$-body…
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…
In this work, we present a symplectic integration scheme to numerically compute space debris motion. Such an integrator is particularly suitable to obtain reliable trajectories of objects lying on high orbits, especially geostationary ones.…
A second order explicit one-step numerical method for the initial value problem of the general ordinary differential equation is proposed. It is obtained by natural modifications of the well-known leapfrog method, which is a second order,…
We suggest a numerical integration procedure for solving the equations of motion of certain classical spin systems which preserves the underlying symplectic structure of the phase space. Such symplectic integrators have been successfully…
Hill's equations are an approximation that is useful in a number of areas of astrophysics including planetary rings and planetesimal disks. We derive a symplectic method for integrating Hill's equations based on a generalized leapfrog. This…
A theory for thermodynamic induction (TI) under isothermal conditions is presented. This includes a treatment of the Helmholtz free energy budget available for a gate variable to utilize towards aiding another variable's approach towards…
The choice of numerical integrator in approximating solutions to dynamic partial differential equations depends on the smallest time-scale of the problem at hand. Large-scale deformations in elastic solids contain both shear waves and bulk…
Coherent or exact equations of motion for a post-Newtonian Lagrangian formalism are the Euler-Lagrange equations without any terms truncated. They naturally conserve energy {and} angular momentum. Doubling the phase-space variables of…
For potential applications in spintronics and quantum computing, it is desirable to place a quantum spin Hall insulator [i.e., a 2D topological insulator (TI)] on a substrate while maintaining a large energy gap. Here, we demonstrate a…
In previous papers, explicit symplectic integrators were designed for nonrotating black holes, such as a Schwarzschild black hole. However, they fail to work in the Kerr spacetime because not all variables can be separable, or not all…
Strong stability preserving (SSP) integrators for initial value ODEs preserve temporal monotonicity solution properties in arbitrary norms. All existing SSP methods, including implicit methods, either require small step sizes or achieve…
Geometric integrators of the Schr\"{o}dinger equation conserve exactly many invariants of the exact solution. Among these integrators, the split-operator algorithm is explicit and easy to implement, but, unfortunately, is restricted to…
In the present article, we investigate the behavior of orbits in a time independent axially symmetric galactic type potential. This dynamical model can be considered to describe the motion in the central parts of a galaxy, for values of…
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…
We prove that a class of A-stable symplectic Runge--Kutta time semidiscretizations (including the Gauss--Legendre methods) applied to a class of semilinear Hamiltonian PDEs which are well-posed on spaces of analytic functions with analytic…
This work addresses the Hamiltonian dynamics of the Kepler problem in a deformed phase space, by considering the equatorial orbit. The recursion operators are constructed and used to compute the integrals of motion. The same investigation…
Numerical evolution of time-dependent differential equations via explicit Runge-Kutta or Taylor methods typically fails to preserve symmetries of a system. It is known that there exists no numerical integration method that in general…