A family of symplectic integrators: Stability, accuracy, and molecular dynamics applications
A family of symplectic integrators: Stability, accuracy, and molecular dynamics applications
复制标题
DOI:
10.1137/s1064827595282350
复制
发表时间:
1997-01-01
影响因子:
3.1
通讯作者:
Schlick, T
中科院分区:
文献类型:
--
作者:
Skeel, RD;Zhang, GH;Schlick, T
The following integration methods for special second-order ordinary differential equations are studied: leapfrog, implicit midpoint, trapezoid, Stormer-Verlet, and Cowell-Numerov. We show that all are members, or equivalent to members, of a one-parameter family of schemes. Some methods have more than one common form, and we discuss a systematic enumeration of these forms. We also present a stability and accuracy analysis based on the idea of ''modified equations'' and a proof of symplecticness. It follows that Cowell-Numerov and ''LIM2'' (a method proposed by Zhang and Schlick) are symplectic. A different interpretation of the values used by these integrators leads to higher accuracy and better energy conservation. Hence, we suggest that the straightforward analysis of energy conservation is misleading.