计算化学公社

 找回密码 Forget password
 注册 Register
Views: 12974|回复 Reply: 9
打印 Print 上一主题 Last thread 下一主题 Next thread

[VASP] 过渡态有多个虚频,尝试了很多办法都解决不了

[复制链接 Copy URL]

28

帖子

0

威望

181

eV
积分
209

Level 3 能力者

{VASP!!!!!!!!!!!!!!!!!表面催化!!!!!!}要疯了,为啥我无论怎么优化,最后频率验证总是有多个虚频啊,而且是没法忽略的那种。试了以下方法:
1,IBRION=1,有的优化了无数步,力很小了。有的收敛了,一算频率。。。2个然后提高收敛精度,还是不行
2,dimer,最后的确优化到只剩一个虚频,但是结构变了,并不是我想要的结构了(大家用dimer的时候没有遇到这种情况吗!!!最后的结构显然不合理)


我要放弃了,是不是有多个虚频就不是过渡态啊,频率这种东西也太难搞了,我直接吐了



另外,为解决这个问题,主要参考了刘大佬的贴子:http://bbs.keinsci.com/thread-13065-1-1.html
针对其提出的过渡态优化不收敛依次尝试以下方法做了以下改进
1,避免不合理的初始结构:我虽然用的nebmake,但所有结构我都检查并将不合理结构手动调整了。
2,CI-NEB插点的数目不合理,emmm,我觉得还行,因为我的距离都差不多,基本每个基元步骤都插4个点,一次CINEB就差不多了。
3,先用NEB或者CINEB进行粗收敛,再用DIMER/CINEB精收敛,NEB+DIMER/CINEB+DIMER都尝试了,最大的问题是DIMER之后结构变得不合理,举个例子CH3的分解,原本CINEB结束H原子(最终落在红圈的地方)在催化剂原子的顶位(侧顶位),dimer完,好家伙,直接给我优化到穴位了,这沃妮马




最后附上这狗日的频率,我本来已经很宽容了,好多基元反应100以下的虚频我都忽略了,可是还是有很多小瘪犊子不行

微信图片_20201015104450.png (23.84 KB, 下载次数 Times of downloads: 28)

微信图片_20201015104450.png

微信图片_20201015103925.png (139.14 KB, 下载次数 Times of downloads: 17)

微信图片_20201015103925.png

微信图片_20201015103918.png (186.2 KB, 下载次数 Times of downloads: 24)

微信图片_20201015103918.png

115

帖子

0

威望

1367

eV
积分
1482

Level 4 (黑子)

2#
发表于 Post on 2020-10-15 14:47:44 | 只看该作者 Only view this author
正在跑dimer,第三步curvature还是正值的我正在瑟瑟发抖(T_T)

20

帖子

0

威望

111

eV
积分
131

Level 2 能力者

3#
发表于 Post on 2023-3-13 16:56:32 | 只看该作者 Only view this author
请问没有虚频怎么办呢?提高了电子和力精度也不行,但是能垒数值和路径是和文献一致的

10

帖子

0

威望

171

eV
积分
181

Level 3 能力者

4#
发表于 Post on 2023-6-27 15:56:20 | 只看该作者 Only view this author
请问楼主是怎么解决的

28

帖子

0

威望

281

eV
积分
309

Level 3 能力者

5#
发表于 Post on 2023-8-2 21:20:33 | 只看该作者 Only view this author
qmlearner 发表于 2020-10-15 14:47
正在跑dimer,第三步curvature还是正值的我正在瑟瑟发抖(T_T)

很正常

28

帖子

0

威望

281

eV
积分
309

Level 3 能力者

6#
发表于 Post on 2023-8-2 21:21:04 | 只看该作者 Only view this author
忽略100以下的虚频,这样会导致零点能不准确

3

帖子

0

威望

295

eV
积分
298

Level 3 能力者

7#
发表于 Post on 2023-8-2 21:35:04 | 只看该作者 Only view this author
可以尝试一下qvasp程序,里面有一个命令qvasp -zc 可以消除过渡态虚频。

28

帖子

0

威望

281

eV
积分
309

Level 3 能力者

8#
发表于 Post on 2023-8-5 10:39:47 | 只看该作者 Only view this author
zhang1999501 发表于 2023-8-2 21:35
可以尝试一下qvasp程序,里面有一个命令qvasp -zc 可以消除过渡态虚频。

结构优化虚频可以消除吗

3

帖子

0

威望

295

eV
积分
298

Level 3 能力者

9#
发表于 Post on 2023-8-5 17:43:55 | 只看该作者 Only view this author
xmuqyj 发表于 2023-8-5 10:39
结构优化虚频可以消除吗

手册中写的是消除过渡态的小虚频,确实有效,但是否可以消除结构优化的虚频,你可以咨询qvasp作者易文才。

211

帖子

0

威望

1466

eV
积分
1677

Level 5 (御坂)

10#
发表于 Post on 2023-8-6 13:56:11 | 只看该作者 Only view this author
xmuqyj 发表于 2023-8-2 21:21
忽略100以下的虚频,这样会导致零点能不准确

您好,如果优化出【同种反应途径,成键角度大的】100以下的唯一虚频,当提升泛函基组进行溶剂化时发现零点能不准了,但是我们也定位到【同种反应途径不同成键角度的过渡态,角度小的其虚频在200-500之间】,同类反应途径随泛函基组的改变没有影响随角度变化变化的趋势。但不同的反应性途径,随着泛函基组的提升,反应的优势路径发生了改变。这种情况需要将我的【大成键角度】过渡态虚频提高才可以吗?

本版积分规则 Credits rule

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

GMT+8, 2026-2-21 13:51 , Processed in 0.189411 second(s), 29 queries , Gzip On.

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