计算化学公社

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

[新手求助] Gaussian软件opt=cartesian对[Th(H2O)9]4+频率计算后出现虚频,如何消除虚频?

[复制链接 Copy URL]

81

帖子

0

威望

808

eV
积分
889

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 wangzh 于 2025-12-20 12:09 编辑

老师们好,我使用Gaussian软件对[Th(H2O)9]4+进行优化和频率计算,发现这个体系在水溶液中一直优化失败,根据Gaussian 常见报错及解决方法;新手求助报错时的注意事项
http://bbs.keinsci.com/forum.php ... =4829&fromuid=42500进行重新优化,发现这个体系是C3点群,使用opt=cartesian才能优化成功,输入文件是1.gif。但是发现有一个虚频,然后我根据帖子Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法 - 思想家公社的门口:量子化学·分子模拟·二次元进行消虚频,一直失败,请问老师我应该怎么解决这个问题?
以下是我已经进行的消除虚频的方向:
2文件:在Gview里打开Edit-Point Group工具,将tolerance设低一些,使Gview能判断出此分子应有的最高阶点群,然后点Symmetrize,即调整结构使得体系的结构严格满足这样的点群。接下来再用Gaussian进行优化的输入文件。
3.out文件:是在1.out最后的结构基础上进行opt=(calcall,tight) 计算。
4.out文件:是在1.out最后的结构基础上进行opt=(maxstep=5,notrustupdate,cartesian,calcall,tight) 计算。
5文件:是在GaussView观看振动模式的界面中选中虚频模式,选上Manual Displacement并拉动滑条来沿着虚频模式略微调整结构,然后点save structure得到新的结构,再拿这个新结构保存新的输入文件,重新优化和做振动分析。
以上方法都失败了,请问老师我应该怎么办?感谢老师的指导!



1.zip

4.39 MB, 下载次数 Times of downloads: 14

2.zip

2.94 MB, 下载次数 Times of downloads: 3

6万

帖子

99

威望

5万

eV
积分
124662

管理员

公社社长

2#
发表于 Post on 2025-12-19 04:50:12 | 只看该作者 Only view this author
按虚频调结构后,maxstep=3并结合int=superfine
北京科音自然科学研究中心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

81

帖子

0

威望

808

eV
积分
889

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2025-12-20 12:08:06 | 只看该作者 Only view this author
本帖最后由 wangzh 于 2025-12-20 12:10 编辑
sobereva 发表于 2025-12-19 04:50
按虚频调结构后,maxstep=3并结合int=superfine

老师您好,我在GaussView观看振动模式的界面中选中虚频模式,选上Manual Displacement并拉动滑条到1.00来沿着频模式进行调整结构。显示优化失败,输入输出文件在压缩文件2,我进行了上传,请问老师我还应该怎么修改?感谢老师解答

76

帖子

2

威望

3504

eV
积分
3620

Level 5 (御坂)

见习砖工

4#
发表于 Post on 2025-12-21 07:11:33 | 只看该作者 Only view this author
本帖最后由 Satoru 于 2025-12-20 19:48 编辑

int=superfine基础上加GDIIS、NoSymm后可以很快收敛,没有虚频
Th-OH2_9_4 .7z (223.14 KB, 下载次数 Times of downloads: 3)

观察压缩包2的优化轨迹可以发现分子朝向不断变化,可能是分子接近球形使得每步优化后标准朝向都会明显改变,在积分格点中的额外旋转引入了噪声。我个人认为这个时候开nosymm是完全合理的


顺便试了一下,不需要int=superfine也可以得到无虚频的结构
Th-OH2_9_4 _ultrafine.7z (283.55 KB, 下载次数 Times of downloads: 2)

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
cokie + 5 给好心人点赞

查看全部评分 View all ratings

人には定めの星がある。綺羅星、明星、流れ星。己の星は見えずとも、見上げる私は今日限り。99期生、星見純那。掴んで見せます、自分星!

81

帖子

0

威望

808

eV
积分
889

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2025-12-22 23:13:23 | 只看该作者 Only view this author
Satoru 发表于 2025-12-21 07:11
int=superfine基础上加GDIIS、NoSymm后可以很快收敛,没有虚频

非常感谢老师指导,我懂啦

本版积分规则 Credits rule

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

GMT+8, 2026-1-23 13:15 , Processed in 0.263589 second(s), 24 queries , Gzip On.

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