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_ip_j+\sum_{i=1}^{n} b_ip_i$, where $a_{ij}=a_{ij}(\textbf{q})$ and $b_i=b_i(\textbf{q})$ are functions of coordinates $\textbf{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.

Symplectic Method Based on Dual Variable Principle

A High Order Symplectic Algorithm Based on Weighted Residual Method

A Symplectic Algorithm For The Stability Analysis Of Nonlinear Parametric Excited Systems

A Symplectic Approximation with Nonlinear Stability and Convergence Analysis for Efficiently Solving Semi-Linear Klein–Gordon Equations

A Symplectic Algorithm for Dynamics of Rigid Body

A Highly Precise Symplectic Direct Integration Method Based on Phase Errors for Hamiltonian Systems

Adjustment of force-gradient operator in symplectic methods

Symplectic finite-difference methods for solving partial differential equations

A variational symplectic scheme based on Simpson's quadrature

Explicit High-Order Noncanonical Symplectic Algorithms for Ideal Two-Fluid Systems

Generalized flow-composed symplectic methods for post-Newtonian Hamiltonian systems

Explicit near-symplectic integrators for post-Newtonian Hamiltonian systems

Conservative Multiquadric Quasi-Interpolation Method for Hamiltonian Wave Equations

Symplectic Runge-Kutta Methods for Hamiltonian Systems Driven by Gaussian Rough Paths

Explicit High-Order Gauge-Independent Symplectic Algorithms for Relativistic Charged Particle Dynamics

Nonstandard Finite Difference Variational Integrators for Multisymplectic Pdes

Symplectic Hamiltonian Finite Element Methods for Semilinear Wave Propagation

High order explicit symplectic integrators for the Discrete Non Linear Schrödinger equation

General Local Energy-Preserving Integrators for Solving Multi-Symplectic Hamiltonian PDEs.

Extended phase-space symplectic-like integrators for coherent post-Newtonian Euler-Lagrange equations