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.

Towards Using Feedback Mechanisms for Numerical Integration of Hamiltonian Systems: A Spring-Mass Benchmark

Feedback Integrators for Mechanical Systems with Holonomic Constraints

Feedback Stabilization of Quasi Nonintegrable Hamiltonian Systems by Using Lyapunov Exponent

A Novel Method for Evaluating the Averaged Drift and Diffusion Coefficients of High DOF Quasi-Non-integrable Hamiltonian Systems

Feedback Integrators for Nonholonomic Mechanical Systems

Hamiltonian Matching for Symplectic Neural Integrators

A new type of energy-preserving integrators for quasi-bi-Hamiltonian systems

Optimal Nonlinear Feedback Control of Quasi-Hamiltonian Systems

Numerical comparisons among some methods for Hamiltonian problems

Explicit near-symplectic integrators for post-Newtonian Hamiltonian systems

Numerical integration methods for large-scale biophysical simulations

Geometric integrators and the Hamiltonian Monte Carlo method

On the numerical integration of the Fokker-Planck equation driven by a mechanical force and the Bismut-Elworthy-Li formula

Adjustment of force-gradient operator in symplectic methods

Simulating Nonholonomic Dynamics

Quasi-Lagrangian equations and its energy-conservative numerical integration for nonlinear dynamic systems

Dynamics and non-integrability of the double spring pendulum

Numerical methods for a periodic optimal control problem

Efficiency analysis of numerical integrations for finite element substructure in real-time hybrid simulation

Implicit-Explicit simulation of Mass-Spring-Charge Systems

An energetically balanced, quasi-Newton integrator for non-hydrostatic vertical atmospheric dynamics