计算化学公社

 找回密码 Forget password
 注册 Register

求助:vasp计算分子与原胞优化无法消除虚频

查看数: 15177 | 评论数: 6 | 收藏 Add to favorites 1
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2021-3-27 23:51

正文摘要:

各位老师好,我是刚开始接触vasp计算的小白。目前我计算了三个模型,分别是晶格各自里面的NO与H2O2分子,还有Fe3O4原胞,分别对上述三个模型进行了opt优化与freq计算,三个模型的opt计算与freq计算的输入文件详见附 ...

回复 Reply

gg302694005 发表于 Post on 2021-3-28 11:37:54
卡开发发 发表于 2021-3-28 11:10
1、暂时也不能判断其他哪些关键字可能造成虚频,只能先进行尝试再根据原理分析,LREAL=T可能会引入一些噪 ...

非常感谢老师耐心细致的解答!
卡开发发 发表于 Post on 2021-3-28 11:10:17
gg302694005 发表于 2021-3-28 10:41
老师,您提的第一点建议,除了设置参数ISPIN=2,把IREAL参数改为默认,还需要设置哪些参数才合理呢。您的 ...

1、暂时也不能判断其他哪些关键字可能造成虚频,只能先进行尝试再根据原理分析,LREAL=T可能会引入一些噪声。
2、VASP的话优化晶格反复拷贝CONTCAR为POSCAR可能对结果会有些许影响。
3、VASP放在盒子里分子有虚频正常,有可能未必消除得了,但不会那么大(如100cm-1以上)。可能的原因是格点太粗糙破缺了分子的旋转不变性,除此之外,人为引入的周期性边界也会使得这种旋转不变性发生破缺,所以有时候这些虚频是难以被消除的,除非网格取得足够密集和盒子取得足够大。

评分 Rate

参与人数
Participants 2
eV +4 收起 理由
Reason
晨晨晨嘛 + 2
sobereva + 2

查看全部评分 View all ratings

gg302694005 发表于 Post on 2021-3-28 10:41:26
卡开发发 发表于 2021-3-28 10:07
1、Fe3O4应该是亚铁磁性的,直接开自旋极化是不足够的,一般情况就别用LREAL=T了。
2、声子计算最好和优化 ...

老师,您提的第一点建议,除了设置参数ISPIN=2,把IREAL参数改为默认,还需要设置哪些参数才合理呢。您的第二点建议,我之前是按照bigbro(大师兄)的教程,如何基于结构优化的INCAR设置的频率计算的INCAR,我之前也按照相关的建议将opt的EDIFF设置为1E-06,但是频率计算结果还是存在虚频。您的第三点建议我马上尝试一下,但是我有一个可能不正确的疑问,如果INCAR的收敛标准不变,那么把CONTCAR变为POSCAR再算一次,结果会不一样吗?最后就是比较邪乎的问题,为什么vasp在空盒子里面算分子,还会出现虚频呢。非常感谢老师的解答
gg302694005 发表于 Post on 2021-3-28 10:33:41
get-it 发表于 2021-3-28 04:40
算分子为什么要用vasp

老师,因为要算吸附,先把分子放盒子里算一遍,然后再算吸附
卡开发发 发表于 Post on 2021-3-28 10:07:56
1、Fe3O4应该是亚铁磁性的,直接开自旋极化是不足够的,一般情况就别用LREAL=T了。
2、声子计算最好和优化结构使用同一精度,一开始优化结构就把SCF收敛标准和构型优化的收敛标准设置严格些。
3、最好多拷贝几次CONTCAR为POSCAR优化。

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
晨晨晨嘛 + 2

查看全部评分 View all ratings

get-it 发表于 Post on 2021-3-28 04:40:15
算分子为什么要用vasp

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-17 20:06 , Processed in 0.176815 second(s), 26 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list