Adjustment of force-gradient operator in symplectic methods

Abstract

Many force-gradient explicit symplectic integration algorithms have been designed for the Hamiltonian $H=T(\mathbf{p})+V(\mathbf{q})$ with kinetic energy $T(\mathbf{p})={\mathbf{p}}^2/2$ in the existing references. When the force-gradient operator is appropriately adjusted as a new operator, they are still suitable for a class of Hamiltonian problems $H=K(\mathbf{p},\mathbf{q})+V(\mathbf{q})$ with \emph{integrable} part $K(\mathbf{p},\mathbf{q})={\sum }_{i=1}^n{\sum }_{j=1}^n{a}_{ij}{p}_i{p}_j+{\sum }_{i=1}^n{b}_i{p}_i$, where $a_{ij}=a_{ij}(\text{q})$ and $b_i=b_i(\text{q})$ are functions of coordinates $\text{q}$. The newly adjusted operator is not a force-gradient operator but is similar to the momentum-version operator associated to the potential $V$. The newly extended (or adjusted) algorithms are no longer solvers of the original Hamiltonian, but are solvers of slightly modified Hamiltonians. They are explicit symplectic integrators with symmetry or time-reversibility. Numerical tests show that the standard symplectic integrators without the new operator are generally poorer than the corresponding extended methods with the new operator in computational accuracies and efficiencies. The optimized methods have better accuracies than the corresponding non-optimized counterparts. Among the tested symplectic methods, the two extended optimized seven-stage fourth-order methods of Omelyan, Mryglod and Folk exhibit the best numerical performance. As a result, one of the two optimized algorithms is used to study the orbital dynamical features of a modified Hénon-Heiles system and a spring pendulum. These extended integrators allow for integrations in Hamiltonian problems, such as the spiral structure in self-consistent models of rotating galaxies and the spiral arms in galaxies.

arXiv comment21 pages, 11 figures

References in corpus (14)

• A New Test for Chaos
• Lyapunov indices with two nearby trajectories in a curved spacetime
• Transition from Regular to Chaotic Circulation in Magnetized Coronae near Compact Objects
• On the construction of high-order force gradient algorithms for integration of motion in classical and quantum systems
• Pseudo-High-Order Symplectic Integrators
• Explicit methods in extended phase space for inseparable Hamiltonian problems
• Optimized Forest-Ruth- and Suzuki-like algorithms for integration of motion in many-body systems
• Chaotic motion and spiral structure in self-consistent models of rotating galaxies
• Construction of Explicit Symplectic Integrators in General Relativity. I. Schwarzschild Black Holes
• On post-Newtonian orbits and the Galactic-center stars
• Symplectic Integration of Post-Newtonian Equations of Motion with Spin
• Generalized, energy-conserving numerical simulations of particles in general relativity. I. Time-like and null geodesics
• A study of symplectic integrators for planetary system problems: error analysis and comparisons
• Analytical forms of chaotic spiral arms

Cited by in corpus (2)

• Impact of electric charges on chaos in magnetized $Reissner-Nordstr\ddot{o}m$ spacetimes
• Study of chaos in rotating galaxies using extended force-gradient symplectic methods