计算化学公社

标题: 求助：VASP优化H3AsO3(亚砷酸)分子SCF不收敛 [打印本页]

作者
Author: 猫小琪 时间: 2021-7-15 11:37
标题: 求助：VASP优化H3AsO3(亚砷酸)分子SCF不收敛
本帖最后由猫小琪于 2021-7-15 11:37 编辑

求助，最近在用VASP优化H3AsO3(亚砷酸)分子，设置了15*15*15大小的盒子，H3AsO3已经过粗优化，采用Gamma点，Ismear=0 ，SIGMA 从0.01测试到0.5都不收敛。每次都是第一个粒子步收敛，第二个离子步scf发散了。下面是输入和输出文件，求各位指点。 (, 下载次数 Times of downloads: 5)

(, 下载次数 Times of downloads: 4)

(, 下载次数 Times of downloads: 6) (, 下载次数 Times of downloads: 5)

(, 下载次数 Times of downloads: 12)

作者
Author: sobereva 时间: 2021-7-18 01:37
如果没有理由非得用VASP，就绝对不要用VASP做这个任务
拿Gaussian、ORCA比VASP算这种普通分子体系强百倍

作者
Author: 乘风万里 时间: 2021-7-19 11:08
本帖最后由乘风万里于 2021-7-19 11:21 编辑

其实按sob说的，如果不是做界面催化这种，没必要用VASP算分子。单算分子的话Gaussian、ORCA算得又快又好。
我看你的INCAR没几个参数，分析一下。
1. 首先你的初始结构就不是很合理，且EDIFF设的太粗糙，第一步算的能量都没合理自洽。我看你算的第一个离子步，自洽完自由能+40多eV，很离谱，所以导致算的Hellmann-Feynman力不准，产生的后续构型更加不合理，根本无法收敛。你把电子自洽精度设大点试试，比如EDIFF = 1E-6，分子电子自洽算得挺快的应该。或者检查下初始结构，用xtb这样的半经验软件优化一下分子构型，我怀疑你之前的粗优化就出问题了。
2. 优化结构过程中NELM设那么大干嘛啊，VASP默认值是60，manual强调说40就够用了，可你设成了1000！你这构型本来就不合理，电子自洽1000步，那不是在白算么。正常的是合理的EDIFF下，自洽60步不收敛，马上变成新的构型继续自洽。不合理的结构优化，多算离子步，少算电子步，按VASP默认参数来不好么？
3. 并行参数NCORE, NPAR设得非常奇怪，刚入门的计算者最好把不是很熟悉的参数都关掉。

作者
Author: 猫小琪 时间: 2021-7-20 09:40

sobereva 发表于 2021-7-18 01:37
如果没有理由非得用VASP，就绝对不要用VASP做这个任务
拿Gaussian、ORCA比VASP算这种普通分子体系强百倍

回复社长，这个是计算As在表面吸附的，所以必须用VASP

作者
Author: 卡开发发 时间: 2021-7-20 10:08

猫小琪发表于 2021-7-20 09:40
回复社长，这个是计算As在表面吸附的，所以必须用VASP

我这边按照你的输入测试下来没啥问题（仅调节了NPAR和NCORE），建议确认一下你们编译的版本是否没有问题。

作者
Author: 猫小琪 时间: 2021-7-20 11:04

乘风万里发表于 2021-7-19 11:08
其实按sob说的，如果不是做界面催化这种，没必要用VASP算分子。单算分子的话Gaussian、ORCA算得又快又好。
...

非常感谢您的回复，我又测试了小的EDIFF，并且用B3LYP/def2-TZVPP优化了分子，调整了POTIM，还是不行。虽然结构收敛，但是分子已经散架了。

作者
Author: 乘风万里 时间: 2021-7-20 11:10

卡开发发发表于 2021-7-20 10:08
我这边按照你的输入测试下来没啥问题（仅调节了NPAR和NCORE），建议确认一下你们编译的版本是否没有问题 ...

大佬，的确是这样的，我自己跑了一下，5个离子步，10s就跑完了，每个电子步不超过10（符合分子自洽认知），可他的OSZICAR中电子步都是跑了1000步，额....这个帖子我还真以为是要讨论收敛问题呢，一顿分析。其是他的POTCAR都错完了，你看他OUTCAR里的记录... (, 下载次数 Times of downloads: 25)

作者
Author: 卡开发发 时间: 2021-7-20 11:20

乘风万里发表于 2021-7-20 11:10
大佬，的确是这样的，我自己跑了一下，5个离子步，10s就跑完了，每个电子步不超过10（符合分子自洽认知） ...

那就难怪了。所以还是要仔细点。

作者
Author: 乘风万里 时间: 2021-7-20 11:27

猫小琪发表于 2021-7-20 11:04
非常感谢您的回复，我又测试了小的EDIFF，并且用B3LYP/def2-TZVPP优化了分子，调整了POTIM，还是不行。虽 ...

跑VASP开始的时候，请注意屏幕的stdout有没有显示POTCAR报错：
'WARNING: type information on POSCAR and POTCAR are incompatible
POTCAR overwrites the type information in POSCAR'

作者
Author: 猫小琪 时间: 2021-7-20 11:49

卡开发发发表于 2021-7-20 11:20
那就难怪了。所以还是要仔细点。

非常感谢大神的测试。最近CP2K用多了，还以为VASP多写几个元素也没问题。再次感谢，我这边计算也没问题了。

作者
Author: 猫小琪 时间: 2021-7-20 11:49

乘风万里发表于 2021-7-20 11:27
跑VASP开始的时候，请注意屏幕的stdout有没有显示POTCAR报错：
'WARNING: type information on POSCAR a ...

非常感谢。是我搞混了CP2K和VASP，我这边计算没问题了。

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