构造求解一阶微分方程组初值问题的八阶龙格-库塔递推公式,结合描述有机分子运动的一维紧束缚模型,研究PPV原子链中极化子的形成及运动.对碳原子数N=160的PPV原子链,由可控步长八阶龙格-库塔公式求解2N(2N+1)=102720个方程组成的方程组,用Fortran语言编程计算,得到稳定的极化子结构和运动图像；在场强E=1×105 V·cm-1的电场作用下,极化子沿分子链的运动速率约为0.2635(A)·fs-1.计算结果表明,八阶龙格-库塔方程可以有效地用于有机分子链中载流子运动的模拟.
八阶龙格-库塔公式 PPV 极化子 Fortran