计算化学公社

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

[VASP] 求助 在计算频率的时候出现了比较小的虚频 求教怎么消除虚频

[复制链接 Copy URL]

19

帖子

0

威望

123

eV
积分
142

Level 2 能力者

这是我的INCAR   

PINLV.png (15.8 KB, 下载次数 Times of downloads: 142)

PINLV.png

3

帖子

0

威望

41

eV
积分
44

Level 2 能力者

35#
发表于 Post on 2024-11-25 16:46:35 | 只看该作者 Only view this author
granvia 发表于 2019-5-28 01:22
我认真看了一下那个VASP教程,其实没有问题,因为他计算得到乙醇的频率有3*9=27个值,其中三个是出现虚数 ...

老师,我明白了,因为最后几个频率对应的是平动和转动自由度,所以对应的能量应该远小于振动能级,所以只要虚频在平动和转动自由度上,且其能量相较于振动能量占比很少即可忽略,上述理解不知道是否有问题,谢谢老师

3

帖子

0

威望

41

eV
积分
44

Level 2 能力者

34#
发表于 Post on 2024-11-25 16:02:06 | 只看该作者 Only view this author
granvia 发表于 2019-5-28 01:22
我认真看了一下那个VASP教程,其实没有问题,因为他计算得到乙醇的频率有3*9=27个值,其中三个是出现虚数 ...

老师你好,那这么说是不是不对应振动的虚频无论多大都可以忽略?

1

帖子

0

威望

39

eV
积分
40

Level 2 能力者

33#
发表于 Post on 2024-10-11 16:08:13 | 只看该作者 Only view this author
sobereva 发表于 2019-5-27 04:02
分子计算,50、100的虚频根本不能随随便便忽略,除非对自己的数据根本不负责任。

1 一方面这已经可能意 ...

我的个人看法:
波数越大的振动越接近局域的谐振子,波数越小的振动谐振子的力常数越接近零,越接近自由平动。波数非常接近零的振动自由度,已经接近于理想自由平动模式了,用简谐近似不合理,但直接忽略也不合理。
如果振动模式是由H、N、O、F元素贡献的,且如果用隐式溶剂模型,比如水溶液,计算振动时并没有考虑振动模式与溶剂氢键网络的耦合。这些接近平动的模式实际上不应该存在,而是与氢键网络的振动耦合,(除非该元素与氢键网络没有相互作用),所以把小波数的振动手动改成水的振动波数可能更合理一些。当然这只是推测,没有严格证明。

44

帖子

0

威望

209

eV
积分
253

Level 3 能力者

32#
发表于 Post on 2023-12-23 09:24:32 | 只看该作者 Only view this author
老师您好,我目前遇到了消虚频问题。通过jmol软件查看,我的孤立体系计算也存在虚频。具体数据在我的最新论坛博文里。
我一开始想通过jmol可视化,将结构按照虚频旋转的方向保存进行重新优化。但是,我尝试保存失败了。
后续,我又尝试vaspkit的507程序,但是显示说“no imaginary-freq”。告诉我不存在虚频,与vaspkit论坛里,写的也不相同。
或许老师,有没有什么好的消虚频的方式?

我的体系中设置EDIFFG = -0.03,EDIFF = 1E-7,ENCUT = 500

14

帖子

0

威望

101

eV
积分
115

Level 2 能力者

31#
发表于 Post on 2023-7-17 15:22:15 | 只看该作者 Only view this author
sobereva 发表于 2023-6-27 23:28
得弄清楚造成的原因
诸如优化收敛已经足够精确,但由于CUTOFF太低等数值因素造成Hessian太糙而产生虚频 ...

好的 非常谢谢社长!

6万

帖子

99

威望

5万

eV
积分
120192

管理员

公社社长

30#
发表于 Post on 2023-6-27 23:28:22 | 只看该作者 Only view this author
nosugar 发表于 2023-6-27 20:49
老师您好,想问一下我在计算一个氧气的吸附体系中,结构优化后计算频率,得到三个虚频都大于50cm-1,那这 ...

得弄清楚造成的原因
诸如优化收敛已经足够精确,但由于CUTOFF太低等数值因素造成Hessian太糙而产生虚频,那不是结构自身的问题
再比如优化时对某些部分做了冻结而导致虚频,那也没法说是不稳定。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

14

帖子

0

威望

101

eV
积分
115

Level 2 能力者

29#
发表于 Post on 2023-6-27 20:49:35 | 只看该作者 Only view this author

老师您好,想问一下我在计算一个氧气的吸附体系中,结构优化后计算频率,得到三个虚频都大于50cm-1,那这样的话是证明我的结构不稳定嘛?但是结构优化已经收敛,EDIFFG = -0.02,EDIFF = 1E-5。

9

帖子

0

威望

177

eV
积分
186

Level 3 能力者

28#
发表于 Post on 2023-6-6 09:14:02 | 只看该作者 Only view this author
飞翔的猪 发表于 2019-5-27 20:53
sob老师 ,打扰您了,我想请问一下:vasp进行零点能校正计算的时候,对于虚频为多少才可以忽略。

vasp ...

我最近用VASP计算过渡态,过渡态的虚频一直在十几或者二十几 i cm-1 , 如果按照小于50可以忽略的话,我这过渡态是不是就不存在?如果按照虚频不能随意忽略,我这过渡态就算是找到了。用了几种过渡态的方法和插不同的点,可以相互验证,基本就是这个过渡态结构了。不知道怎么判断了?

6万

帖子

99

威望

5万

eV
积分
120192

管理员

公社社长

27#
发表于 Post on 2022-8-12 06:37:08 | 只看该作者 Only view this author
nemolee 发表于 2022-8-11 22:20
请问Gaussian中计算的小分子这种平动转动的虚频是可以忽略的吗?谢谢。
类似这种:
Low frequencies - ...

可以
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

60

帖子

0

威望

867

eV
积分
927

Level 4 (黑子)

26#
发表于 Post on 2022-8-11 22:20:51 | 只看该作者 Only view this author

请问Gaussian中计算的小分子这种平动转动的虚频是可以忽略的吗?谢谢。
类似这种:
Low frequencies ---   -4.7750   -3.6574   -0.0018    0.0005    0.0017    2.7154
Low frequencies ---   23.8793   31.3131   42.5335

22

帖子

0

威望

119

eV
积分
141

Level 2 能力者

25#
发表于 Post on 2020-1-10 15:25:45 | 只看该作者 Only view this author

收到,非常感谢!

3754

帖子

3

威望

1万

eV
积分
19678

Level 6 (一方通行)

围观吃瓜群众

24#
发表于 Post on 2020-1-8 22:43:24 | 只看该作者 Only view this author
sige 发表于 2020-1-8 15:20
请教卡老师,怎么增加FFT网格密度呢,增大K点吗

动能截断。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

22

帖子

0

威望

119

eV
积分
141

Level 2 能力者

23#
发表于 Post on 2020-1-8 15:20:29 | 只看该作者 Only view this author
卡开发发 发表于 2019-5-27 20:52
其实我觉得比较莫名其妙的地方还是为啥要区分分子振动和固体振动,还要构造两种不同的标准,如果虚频值正好 ...

请教卡老师,怎么增加FFT网格密度呢,增大K点吗

6万

帖子

99

威望

5万

eV
积分
120192

管理员

公社社长

22#
发表于 Post on 2020-1-3 20:36:22 | 只看该作者 Only view this author
xxzj 发表于 2020-1-3 19:02
老师,您好,我对分子结构优化后计算了频率,有一个虚频为-0.66,可以忽略吗?

绝对不能
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

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

GMT+8, 2025-8-18 01:35 , Processed in 0.200498 second(s), 30 queries , Gzip On.

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