计算化学公社

标题: 求助苯炔自由基B3LYP优化后有虚频 [打印本页]
作者
Author: 哒咩    时间: 2025-2-20 15:43
标题: 求助苯炔自由基B3LYP优化后有虚频
(, 下载次数 Times of downloads: 4) (, 下载次数 Times of downloads: 4) 请问各位老师我用uB3LYP/6-31g（d）优化苯炔自由基并计算频率时有-143的虚频，已经尝试手动调整过结构，但虚频越改越大，这该怎么解决呢？
作者
Author: sobereva    时间: 2025-2-20 16:35
按虚频调结构，调得幅度别太小，并且加上calcfc。若还有问题上传输出文件
作者
Author: 哒咩    时间: 2025-2-20 17:12
老师我按您的方法添加了关键词在opt文件里，但是附件类型限制我上传opt的输出文件，我把opt的输入文件和频率的输出文件已上传，调完结构后虚频反而更大了
作者
Author: 哒咩    时间: 2025-2-20 17:12
sobereva 发表于 2025-2-20 16:35
按虚频调结构，调得幅度别太小，并且加上calcfc。若还有问题上传输出文件

好的老师已上传相关文件
作者
Author: sobereva    时间: 2025-2-20 17:38
得到了无虚频的弯曲的结构
别写没必要的关键词，自找麻烦
作者
Author: Uus/pMeC6H4-/キ    时间: 2025-2-20 17:47
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
复制代码


有关几何优化和波函数稳定的问题，参考这帖二楼。
作者
Author: 哒咩    时间: 2025-2-20 21:39
Uus/pMeC6H4-/キ 发表于 2025-2-20 17:47
Ph-C≡C·应该叫苯乙炔基自由基，“苯炔”（Benzyne）指的是苯环去掉两个氢所成的活泼物种，不能混淆。nosy ...

好的谢谢
作者
Author: 哒咩    时间: 2025-2-20 21:39
sobereva 发表于 2025-2-20 17:38
得到了无虚频的弯曲的结构
别写没必要的关键词，自找麻烦

好的感谢老师
作者
Author: Uus/pMeC6H4-/キ    时间: 2025-2-20 23:26
sobereva 发表于 2025-2-20 17:38
得到了无虚频的弯曲的结构
别写没必要的关键词，自找麻烦

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

作者
Author: wzkchem5    时间: 2025-2-21 05:00
Uus/pMeC6H4-/キ 发表于 2025-2-20 16:26
社长，顺便借此机会问一下，这种体系的优化需要考虑多参考特征么？我也是第一次接触，用UB3LYP-D3(BJ)/de ...

假如要确定炔基到底在不在苯环平面内，有没有弯曲等等，我觉得还是用多参考态方法比较保险。假如只是为了得到一个结构用于算高精度单点能，我觉得DFT问题不大
作者
Author: sobereva    时间: 2025-2-21 06:38
Uus/pMeC6H4-/キ 发表于 2025-2-20 23:26
社长，顺便借此机会问一下，这种体系的优化需要考虑多参考特征么？我也是第一次接触，用UB3LYP-D3(BJ)/de ...

T1无法用来判断某种DFT泛函是否一定不适用，当前只能说HF很不适合
不十分讲究的话，可以几种常用泛函（B3LYP、PBE0、wB97XD等）都做opt freq，结果基本都对上就OK了。特别讲究、当前体系的计算精度是研究的问题的重点的话，应当用多参考方法确认




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3