Higher-order force gradient symplectic algorithms

被引:21
作者
Chin, SA [1 ]
Kidwell, DW [1 ]
机构
[1] Texas A&M Univ, Dept Phys, Ctr Theoret Phys, College Stn, TX 77843 USA
基金
美国国家科学基金会;
关键词
D O I
10.1103/PhysRevE.62.8746
中图分类号
O35 [流体力学]; O53 [等离子体物理学];
学科分类号
070204 ; 080103 ; 080704 ;
摘要
We show that a recently discovered fourth order symplectic algorithm, which requires one evaluation of force gradient in addition to three evaluations of the force, when iterated to higher order, yielded algorithms that are far superior to similarly iterated higher order algorithms based on the standard Forest-Ruth algorithm. We gauge the accuracy of each algorithm by comparing the step-size independent error functions associated with energy conservation and the rotation of the Laplace-Runge-Lenz vector when solving a highly eccentric Kepler problem. For orders 6, 8, 10, and 12, the new algorithms are approximately a factor of 10(3), 10(4), 10(4), and 10(5) better.
引用
收藏
页码:8746 / 8752
页数:7
相关论文
共 14 条