Symplectic integrators with potential derivatives to third order  

作　　者：Wei Sun Xin Wu Guo-Qing Huang 

机构地区：[1]Institute of Astronomy ＆ Department of Physics, Nanchang University, Nanchang 330031, China

出　　处：《Research in Astronomy and Astrophysics》2011年第3期353-368,共16页天文和天体物理学研究（英文版）

基　　金：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

摘　　要：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.

关 键 词：celestial mechanics—methods：numerical 

分 类 号：O636.1[理学—高分子化学] O175.8[理学—化学]