计算化学公社

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

[新手求助] 求助苯炔自由基B3LYP优化后有虚频

[复制链接 Copy URL]

16

帖子

0

威望

99

eV
积分
115

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
请问各位老师我用uB3LYP/6-31g(d)优化苯炔自由基并计算频率时有-143的虚频,已经尝试手动调整过结构,但虚频越改越大,这该怎么解决呢?

6万

帖子

99

威望

5万

eV
积分
120158

管理员

公社社长

2#
发表于 Post on 2025-2-20 16:35:27 | 只看该作者 Only view this author
按虚频调结构,调得幅度别太小,并且加上calcfc。若还有问题上传输出文件
北京科音自然科学研究中心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

16

帖子

0

威望

99

eV
积分
115

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2025-2-20 17:12:02 | 只看该作者 Only view this author
老师我按您的方法添加了关键词在opt文件里,但是附件类型限制我上传opt的输出文件,我把opt的输入文件和频率的输出文件已上传,调完结构后虚频反而更大了

202502201710559017..png (73.32 KB, 下载次数 Times of downloads: 4)

202502201710559017..png

2a-radical-calcfc-freq.log

123.39 KB, 下载次数 Times of downloads: 4

2a-radical-calcfc.gjf

928 Bytes, 下载次数 Times of downloads: 3

16

帖子

0

威望

99

eV
积分
115

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2025-2-20 17:12:48 | 只看该作者 Only view this author
sobereva 发表于 2025-2-20 16:35
按虚频调结构,调得幅度别太小,并且加上calcfc。若还有问题上传输出文件

好的老师已上传相关文件

6万

帖子

99

威望

5万

eV
积分
120158

管理员

公社社长

5#
发表于 Post on 2025-2-20 17:38:32 | 只看该作者 Only view this author
得到了无虚频的弯曲的结构
别写没必要的关键词,自找麻烦

optfreq.gjf

1006 Bytes, 下载次数 Times of downloads: 11

optfreq.out

415.34 KB, 下载次数 Times of downloads: 8

北京科音自然科学研究中心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

894

帖子

3

威望

1685

eV
积分
2639

Level 5 (御坂)

傻傻的木瓜

6#
发表于 Post on 2025-2-20 17:47:06 | 只看该作者 Only view this author
Ph-C≡C·应该叫苯乙炔基自由基,“苯炔”(Benzyne)指的是苯环去掉两个氢所成的活泼物种,不能混淆。nosymm、g09defaults、只有6-31G*的genecp都不是普通的几何优化必需的关键词,若有特意使用的理由(比如为了计算键能之类的量,从金属炔基配合物里取出结构优化)最好也说清楚。

仅供参考,如果直接在GaussView画一个苯甲腈再把氮换成碳,用UB3LYP-D3(BJ)/6-311G(d,p)优化并振动分析的结果会有一个34.47i的虚频,振动方向使C(苯环)-C≡C(末端)角度在苯环平面内偏离180度。微调角度可以优化出另一个C(苯环)-C≡C(末端)约170.1度、炔基仍与苯环共面的无虚频结构,但是波函数不稳定。用stable=opt找到该级别下稳定波函数再重新优化,最终C(苯环)-C≡C(末端)约168度、炔基翘出苯环平面,无虚频且波函数稳定。坐标如下。

  1. C                  1.42594700   -1.21512300   -0.01283100
  2. C                  0.04383700   -1.22188900    0.02057600
  3. C                 -0.67775600    0.00029000    0.04106800
  4. C                  0.04425300    1.22215300    0.02053700
  5. C                  1.42634700    1.21483900   -0.01283100
  6. C                  2.11934700   -0.00026700   -0.02939600
  7. H                  1.97304800   -2.14986600   -0.02754700
  8. H                 -0.51325800   -2.14971900    0.03180400
  9. H                 -0.51242100    2.15023500    0.03175200
  10. H                  1.97377600    2.14938700   -0.02757100
  11. H                  3.20267900   -0.00045800   -0.05701700
  12. C                 -2.07378500    0.00021400    0.09209300
  13. C                 -3.32882700   -0.00014700   -0.11112000
复制代码


有关几何优化和波函数稳定的问题,参考这帖二楼。
√546=23.36664289109

16

帖子

0

威望

99

eV
积分
115

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2025-2-20 21:39:32 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2025-2-20 17:47
Ph-C≡C·应该叫苯乙炔基自由基,“苯炔”(Benzyne)指的是苯环去掉两个氢所成的活泼物种,不能混淆。nosy ...

好的谢谢

16

帖子

0

威望

99

eV
积分
115

Level 2 能力者

8#
 楼主 Author| 发表于 Post on 2025-2-20 21:39:53 | 只看该作者 Only view this author
sobereva 发表于 2025-2-20 17:38
得到了无虚频的弯曲的结构
别写没必要的关键词,自找麻烦

好的感谢老师

894

帖子

3

威望

1685

eV
积分
2639

Level 5 (御坂)

傻傻的木瓜

9#
发表于 Post on 2025-2-20 23:26:13 | 只看该作者 Only view this author
sobereva 发表于 2025-2-20 17:38
得到了无虚频的弯曲的结构
别写没必要的关键词,自找麻烦

社长,顺便借此机会问一下,这种体系的优化需要考虑多参考特征么?我也是第一次接触,用UB3LYP-D3(BJ)/def2-TZVP重新优化得到无虚频且波函数稳定的结构,然后用CCSD(T1diag)/def2TZVP关键词做冻核计算:
PhCC_t1diag.log (92.02 KB, 下载次数 Times of downloads: 2)
结果UHF的<S**2>= 1.3266,CCSD迭代收敛后T1 Diagnostic =     0.03969207,这个T1诊断值已经比http://bbs.keinsci.com/thread-5460-1-1.html的不少体系高了。

√546=23.36664289109

1万

帖子

0

威望

8999

eV
积分
20755

Level 6 (一方通行)

10#
发表于 Post on 2025-2-21 05:00:02 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2025-2-20 16:26
社长,顺便借此机会问一下,这种体系的优化需要考虑多参考特征么?我也是第一次接触,用UB3LYP-D3(BJ)/de ...

假如要确定炔基到底在不在苯环平面内,有没有弯曲等等,我觉得还是用多参考态方法比较保险。假如只是为了得到一个结构用于算高精度单点能,我觉得DFT问题不大
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

6万

帖子

99

威望

5万

eV
积分
120158

管理员

公社社长

11#
发表于 Post on 2025-2-21 06:38:03 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2025-2-20 23:26
社长,顺便借此机会问一下,这种体系的优化需要考虑多参考特征么?我也是第一次接触,用UB3LYP-D3(BJ)/de ...

T1无法用来判断某种DFT泛函是否一定不适用,当前只能说HF很不适合
不十分讲究的话,可以几种常用泛函(B3LYP、PBE0、wB97XD等)都做opt freq,结果基本都对上就OK了。特别讲究、当前体系的计算精度是研究的问题的重点的话,应当用多参考方法确认
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2025-8-16 09:04 , Processed in 0.189431 second(s), 23 queries , Gzip On.

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