Adjustment of Force–Gradient Operator in Symplectic Methods

Lina Zhang, Xin Wu and Enwei Liang
Additional contact information
Lina Zhang: School of Physical Science and Technology, Guangxi University, Nanning 530004, China
Xin Wu: School of Physical Science and Technology, Guangxi University, Nanning 530004, China
Enwei Liang: School of Physical Science and Technology, Guangxi University, Nanning 530004, China

Mathematics, 2021, vol. 9, issue 21, 1-21

Abstract: Many force–gradient explicit symplectic integration algorithms have been designed for the Hamiltonian H = T ( p ) + V ( q ) with kinetic energy T ( p ) = p 2 / 2 in the existing references. When a force–gradient operator is appropriately adjusted as a new operator, it is still suitable for a class of Hamiltonian problems H = K ( p , q ) + V ( q ) with integrable part K ( p , q ) = ? i = 1 n ? j = 1 n a i j p i p j + ? i = 1 n b i p i , where a i j = a i j ( q ) and b i = b i ( q ) are functions of coordinates 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.

Keywords: symplectic integration; force gradient; chaos; Hamiltonian systems (search for similar items in EconPapers)
JEL-codes: C (search for similar items in EconPapers)
Date: 2021
References: View complete reference list from CitEc
Citations:

Downloads: (external link)
https://www.mdpi.com/2227-7390/9/21/2718/pdf (application/pdf)
https://www.mdpi.com/2227-7390/9/21/2718/ (text/html)

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:gam:jmathe:v:9:y:2021:i:21:p:2718-:d:665495

Access Statistics for this article

Mathematics is currently edited by Ms. Emma He

More articles in Mathematics from MDPI
Bibliographic data for series maintained by MDPI Indexing Manager ().