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

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

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

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

猫小琪 发表于 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'

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

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

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

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

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

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

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

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

其实按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设得非常奇怪，刚入门的计算者最好把不是很熟悉的参数都关掉。

如果没有理由非得用VASP，就绝对不要用VASP做这个任务
拿Gaussian、ORCA比VASP算这种普通分子体系强百倍