计算化学公社

标题: IBRION=1 or 2和DOSCAR之间有关系吗 [打印本页]

作者
Author: gur0bi 时间: 2022-2-2 12:37
标题: IBRION=1 or 2和DOSCAR之间有关系吗
研究生刚刚接触VASP，想求助大家一个问题，我目前有几组geometry optimizaiton的计算结果，正好包括IBRION=1和两组IBRION=2的计算。其中在IBRION=1的设置下，DOSCAR几乎是空的没有dos信息（如附图），而IBRION=2则是生成了正常的DOSCAR。其他INCAR中的设置都是一样的，很疑惑为什么这其中会有联系？大家有遇到过类似的情况吗？

作者
Author: 卡开发发 时间: 2022-2-2 22:05
要不然两组数据打包或者给个网盘链接过来看看？我不确定你IBRION=1的构型优化是否收敛。

作者
Author: gur0bi 时间: 2022-2-3 13:32

卡开发发发表于 2022-2-2 22:05
要不然两组数据打包或者给个网盘链接过来看看？我不确定你IBRION=1的构型优化是否收敛。

十分感谢，上传了IBRION1和IBRION2的两个压缩包，各包括两组计算，其中IBRION1计算中的两个STEP2计算结果里的DOSCAR有问题
链接：https://pan.baidu.com/s/1qnATa63COdr5v7NlPGGXdw
提取码：h3wm

作者
Author: 卡开发发 时间: 2022-2-3 19:49

gur0bi 发表于 2022-2-3 13:32
十分感谢，上传了IBRION1和IBRION2的两个压缩包，各包括两组计算，其中IBRION1计算中的两个STEP2计算结果 ...

我翻了一下辕马，在main.F：

ELSE IF (DYN%IBRION/=1) THEN
...
WRITE(16, ... #16是DOSCAR的io number
CALL XML_DOS(... #XML_DOS是将dos写入xml的子程序
...
ENDIF

复制代码

看上去IBRION=1确实不能输出DOS信息，但我没有深究前后其他程序结构以及不知道VASP这么设计的原因是什么。目前来说，建议上还是跑个单点能自洽来获得DOS比较合理。

作者
Author: gur0bi 时间: 2022-2-4 05:50

卡开发发发表于 2022-2-3 19:49
我翻了一下辕马，在main.F：

看上去IBRION=1确实不能输出DOS信息，但我没有深究前后其他程序结构以及 ...

懂了，感谢指点！

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