计算化学公社

标题: 请教OUTCAR中力的收敛输出 [打印本页]

作者
Author: 菩城 时间: 2017-10-16 18:24
标题: 请教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？
请各位前辈解答，谢谢！

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

作者
Author: 菩城 时间: 2017-10-17 08:05

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

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

欢迎光临计算化学公社 (http://ccc.keinsci.com/)