求助：vasp计算结构优化耗时过长

Mr.zhen

vasp新手，在尝试模拟羟基吸附至g-C3N4表面的过程时，发现算到后面需要两天多CONTCAR才会更新一次。计算是通过material studio建模生成的POSCAR，如图，建立的模型为3层，但下面两层使用material studio进行了固定，所以应该是只对第一层的近60个原子进行优化。使用32核进行计算，INCAR如下：
SYSTEM = absorption
ENCUT = 450
ISMEAR = 0
SIGMA = 0.05
EDIFFG = 0.01
NSW = 500
IBRION = 2
POTIM = 0.5
ISIF = 2
EDIFF = 1E-5
ALGO = VeryFast
LCHARG = .FALSE.
NCORE = 8
KPAR = 4  
NELMIN = 6
ISPIN = 2
LORBIT = 11
LREAL = A

不知这么长的耗时是我哪里设置得不当，还是原本就应是耗时很长，希望大家能够指点一二。

Aletyx

固定原子是在POSCAR里面，MS里面固定没有用。此外你的EIFFG=0.01是正确的吗

Mr.zhen

Aletyx 发表于 2023-6-5 20:12
固定原子是在POSCAR里面，MS里面固定没有用。此外你的EIFFG=0.01是正确的吗

MS里固定后导出.xsd文件，然后通过vaspkit可以写出有固定原子的POSCAR。关于INCAR里各参数包括EDIFFG的取值，还在尝试中，一般源于他人的经验

乐平

你可以看看每个电子步的耗时是多少，命令如下：

1. grep LOOP OUTCAR

复制代码

另外，任何时候都要先查手册
https://www.vasp.at/wiki/index.php/EDIFFG

如果 EDIFFG 设置是正值，表示代表弛豫计算过程中相邻两步的总能量变化小于你设置的值时，计算收敛并停止：
如果 EDIFFG 设置为负值，则代表弛豫计算过程中所有受力都小于你设置的绝对值时，计算收敛并停止。这是更方便的设置。

Mr.zhen

乐平 发表于 2023-6-5 22:31
你可以看看每个电子步的耗时是多少，命令如下：

耗时的话是下面这样，算是正常的耗时吗
      LOOP:  cpu time   53.2074: real time   53.2109
      LOOP:  cpu time   50.0015: real time   50.0046
      LOOP:  cpu time   48.7522: real time   48.7546
      LOOP:  cpu time   51.5022: real time   51.5043
      LOOP:  cpu time   50.4074: real time   50.4128
      LOOP:  cpu time   49.8784: real time   49.8799
      LOOP:  cpu time   47.5015: real time   47.5032
      LOOP:  cpu time   41.6458: real time   41.6479
      LOOP:  cpu time   35.3221: real time   35.3237
      LOOP:  cpu time   36.0267: real time   36.0328
     LOOP+:  cpu time  541.6151: real time  541.6489
      LOOP:  cpu time   53.2264: real time   53.2305
      LOOP:  cpu time   47.5270: real time   47.5283
      LOOP:  cpu time   46.1970: real time   46.2001
      LOOP:  cpu time   52.6746: real time   52.6765
      LOOP:  cpu time   53.0061: real time   53.0073
      LOOP:  cpu time   52.4842: real time   52.4860
      LOOP:  cpu time   51.1936: real time   51.1970
      LOOP:  cpu time   45.9841: real time   45.9854
      LOOP:  cpu time   41.7399: real time   41.7420

EDIFFG之前看手册知道是作为收敛停止的标准，但是具体应该设多大还不确定，手册上说默认是EDIFF的10倍，但似乎很多人都自己设定一个其他值，有正有负，大小不一，是需要根据自己研究的体系多算几次慢慢确定吗，具体结果变成什么样算是合适呢？

Stars

Mr.zhen 发表于 2023-6-6 19:22
耗时的话是下面这样，算是正常的耗时吗
      LOOP:  cpu time   53.2074: real time   53.2109
       ...

EDIFF本人强烈建议设置严格一些，至少1E-5，可以提高到1E-6，电子步严格可以让受力更准确，对结构优化有利。
EDIFFG建议用受力为标准，就是-多少，个人意见是-0.02