请教OUTCAR中力的收敛输出

菩城

我在OUTCAR中找到三个力的输出，第一个在TOTAL-FORCE下面：
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      4.30570     11.19969      1.73337        -0.017656     -0.021652      0.012046
      7.07248     10.12773      1.80333        -0.067065      0.131452     -0.015895
      1.43307      5.82394      1.70501         0.012386      0.100582      0.002292

------------------------------------------------------------------------------------
不知为何一个原子坐标对应三个力？分别代表三个方向的受力吗？

还有两个力的输出在能量的后面：
  free  energy   TOTEN  =      -314.315861 eV
  energy  without entropy=     -314.315861  energy(sigma->0) =     -314.315861
d Force = 0.2253341E-01[ 0.101E-01, 0.350E-01]  d Energy = 0.2410555E-01-0.157E-02
d Force =-0.3782950E+01[-0.383E+01,-0.374E+01]  d Ewald  =-0.2610905E+01-0.117E+01

不清楚为何会出现两个d Force？如果看力的收敛情况，应该看哪个d Force？
请各位前辈解答，谢谢！

菩城

卡开发发 发表于 2017-10-16 19:56
你说的这两个d force的说明我没有找到明确的来源，如果非要弄懂可能得专门去扒代码，不过最近懒得做这件事 ...

原来如此，非常感谢卡开发发老师！

卡开发发

你说的这两个d force的说明我没有找到明确的来源，如果非要弄懂可能得专门去扒代码，不过最近懒得做这件事。但可以明确的是，vasp对EDIFFG的判断是依据每个原子，也就是∀|Fi|<EDIFFG，i=1,Natom来判断的，也就等同于要求Max(|Fi|)<EDIFFG，其中的|Fi|就是sqrt(Fx^2+Fy^2+Fz^2)，也就是你说的那三个力。