摘要
An operator associated with third-order potential derivatives and a force gradient operator corresponding to second-order potential derivatives are used together to design a number of new fourth-order explicit symplectic integrators for the natural splitting of a Hamiltonian into both the kinetic energy with a quadratic form of momenta and the potential energy as a function of position coordinates.Numerical simulations show that some new optimal symplectic algorithms are much better than their non-optimal counterparts in terms of accuracy of energy and position calculations.
An operator associated with third-order potential derivatives and a force gradient operator corresponding to second-order potential derivatives are used together to design a number of new fourth-order explicit symplectic integrators for the natural splitting of a Hamiltonian into both the kinetic energy with a quadratic form of momenta and the potential energy as a function of position coordinates.Numerical simulations show that some new optimal symplectic algorithms are much better than their non-optimal counterparts in terms of accuracy of energy and position calculations.
基金
supported by the National Natural Science Foundation of China(Grant No.10873007)
supported by the Science Foundation of Jiangxi Education Bureau(GJJ09072)
the Program for an Innovative Research Team of Nanchang University