计算化学公社

标题: 求助：vasp计算结构优化耗时过长 [打印本页]

作者
Author: Mr.zhen 时间: 2023-6-5 17:13
标题: 求助：vasp计算结构优化耗时过长
本帖最后由 Mr.zhen 于 2023-6-7 14:38 编辑

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

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

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

作者
Author: Mr.zhen 时间: 2023-6-5 20:20

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

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

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

grep LOOP OUTCAR

复制代码

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

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

作者
Author: Mr.zhen 时间: 2023-6-6 19:22

乐平发表于 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倍，但似乎很多人都自己设定一个其他值，有正有负，大小不一，是需要根据自己研究的体系多算几次慢慢确定吗，具体结果变成什么样算是合适呢？

作者
Author: Stars 时间: 2023-12-15 23:55

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

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

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