[交流] 【求助完毕】VASP结构优化时离子步原子是如何移动的？

S(E；F1，F2，...， Fn)=f(r1，r2， ...， rn).

S(E；F1，F2，...， Fn)=
S'(E)=f(R=r1，r2， ...， rn).

R(n=new)=R(o=old)+deltaR

问题是，在S'(E)=f(R)方程具体形式未知的情况下，如何确定deltaR呢？

Fx=-deltaE/deltaRx
deltaE>0,向x轴正方向移动势能增加，原子受力应该指向x负方向,Fx<0；
deltaE<0,向x轴正方向移动势能降低，原子受力应该指向X正方向,Fx>0。
所以，力Fx以x轴正方向为正，负方向为负。其它以此类推。