捐赠本论坛 (Donate this site)
开启辅助访问

计算化学公社

 找回密码 Forget password
 注册 Register
计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 量子化学 (Quantum Chemistry) › 求助关于G16中叔丁氧自由基对利多卡因碳氢键的反应能， ...

求助关于G16中叔丁氧自由基对利多卡因碳氢键的反应能，先结构优化，但SCF收敛失败

查看数: 181 | 评论数: 10 | 收藏 Add to favorites 0
关灯 | 提示：支持键盘翻页<-左 右->
帖子模式
    组图打开中，请稍候......
kkkkkkkkkk
发布时间: 2026-1-14 12:14

正文摘要:

本帖最后由 kkkkkkkkkk 于 2026-1-14 12:19 编辑 gjf文件代码 %nprocshared=1 %mem=6GB %chk=ts-6.chk # stable=opt uwb97xd/6-311++g(d,p) ts 0 2 C           ...

回复 Reply

kkkkkkkkkk 发表于 Post on 2026-1-15 10:45:11
sobereva 发表于 2026-1-15 02:03
gdiis和Cartesian不兼容，下文明确说了
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/ ...

好的好的！谢谢！
北大-陶豫 发表于 Post on 2026-1-15 10:38:35
sobereva 发表于 2026-1-15 02:01
PM7自身已经带了色散校正项，没法再加D3校正（即便程序不报错，也未必结果有意义，自己也不知道实际算的 ...

谢谢指正！我平时都只把它当作粗略优化，优化得差不多就行，都没关注过这些
sobereva 发表于 Post on 2026-1-15 02:03:47
gdiis和Cartesian不兼容，下文明确说了
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

甭写多余的noraman，下文说了
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331（http://bbs.keinsci.com/thread-3460-1-1.html）

弥散加的完全多余，尤其是氢的弥散。仔细看下文
谈谈弥散函数和“月份”基组
http://sobereva.com/119
sobereva 发表于 Post on 2026-1-15 02:01:31
北大-陶豫 发表于 2026-1-14 13:26
当前结构太烂，容易导致电子结构诡异，难以收敛，再加上你基组还用弥散函数，更是雪上加霜。
好好看：谈谈 ...

PM7自身已经带了色散校正项，没法再加D3校正（即便程序不报错，也未必结果有意义，自己也不知道实际算的是什么）
可以直接用PM7，或者PM6D3关键词
kkkkkkkkkk 发表于 Post on 2026-1-14 14:28:48
BinWang 发表于 2026-1-14 13:03
看了一眼你的结构，初猜太烂，重新构建模型吧

谢谢！
kkkkkkkkkk 发表于 Post on 2026-1-14 14:28:24
BinWang 发表于 2026-1-14 12:54
难收敛的话，先减小基组，生成一个过渡态，之后将该过渡态作为初猜用大基组重新找过渡态。

好的好的！谢谢！
kkkkkkkkkk 发表于 Post on 2026-1-14 13:57:28
北大-陶豫 发表于 2026-1-14 13:26
当前结构太烂，容易导致电子结构诡异，难以收敛，再加上你基组还用弥散函数，更是雪上加霜。
好好看：谈谈 ...

谢谢！
北大-陶豫 发表于 Post on 2026-1-14 13:26:45
当前结构太烂，容易导致电子结构诡异，难以收敛，再加上你基组还用弥散函数，更是雪上加霜。
好好看：谈谈量子化学中基组的选择 http://sobereva.com/336
刚搭出来的结构烂不要紧，用 PM7D3 等方法粗略优化一下即可。关键词：# opt PM7 Em=GD3
我给你这个体系用 PM7D3 粗略优化了一下，能量下降了 415.38 kcal，结构也正常多了，优化好的结构如下，你可以直接拿来用。
  1. C                  3.69016700   -1.01660300    0.55845800
  2. C                  4.03405700   -0.15834300   -0.47966300
  3. C                  3.10688100    0.14721100   -1.47283400
  4. C                  1.82611700   -0.40347300   -1.43769800
  5. C                  1.48913800   -1.28397200   -0.39225000
  6. C                  2.41741800   -1.59168700    0.61531000
  7. C                  2.07763600   -2.49695500    1.74697100
  8. N                  0.19640900   -1.89738500   -0.38555200
  9. C                  0.83118600   -0.02474000   -2.47826000
  10. C                 -0.73217500   -1.59695700    0.59582000
  11. C                 -2.12003600   -2.22503100    0.41919200
  12. O                 -0.44609600   -0.91406400    1.55832000
  13. N                 -3.05320600   -1.34227400   -0.29281400
  14. C                 -3.92153900   -0.51625800    0.56102000
  15. C                 -2.63729100   -0.88160800   -1.51599400
  16. C                 -3.24092300    0.65465700    1.27006400
  17. C                 -3.00665200    0.43440700   -2.07353900
  18. H                  4.41571700   -1.24366500    1.33613200
  19. H                  5.02960500    0.27986600   -0.51483600
  20. H                  3.37960600    0.82793500   -2.27566500
  21. H                  1.58908700   -1.92728400    2.55921300
  22. H                  2.95993500   -2.98714600    2.17550700
  23. H                  1.37370800   -3.28826900    1.45349700
  24. H                 -0.11982300   -2.25363200   -1.28143100
  25. H                  1.23940000    0.67422100   -3.22010000
  26. H                  0.44074200   -0.88781600   -3.02946300
  27. H                 -0.03156500    0.49784300   -2.01368900
  28. H                 -2.05653300   -3.19138000   -0.13118700
  29. H                 -2.52114800   -2.46955000    1.43104300
  30. H                 -4.39458600   -1.19073600    1.31195400
  31. H                 -4.75338400   -0.13949700   -0.08030400
  32. H                 -0.47726200    0.79810200    0.45714000
  33. H                 -2.03372800   -1.54059900   -2.11781000
  34. H                 -2.73936700    1.32433900    0.55756700
  35. H                 -2.47473000    0.32045600    1.98045800
  36. H                 -3.96557300    1.25510300    1.82734500
  37. H                 -4.06864600    0.67847100   -1.93472000
  38. H                 -2.42086400    1.24200900   -1.59243100
  39. H                 -2.80220700    0.48973200   -3.15093800
  40. C                  0.34977600    2.60788700    0.55402100
  41. C                 -0.59337000    3.37944400    1.48729600
  42. C                  1.49191800    1.95377400    1.33414600
  43. O                 -0.42255700    1.60624700   -0.08990900
  44. C                  0.86468700    3.47947800   -0.59585300
  45. H                 -1.47162200    3.73522200    0.93762100
  46. H                 -0.09679100    4.24264900    1.93509000
  47. H                 -0.95544200    2.73579600    2.29602600
  48. H                  2.20895600    1.47512900    0.65481500
  49. H                  2.03511800    2.67642500    1.94540400
  50. H                  1.11937900    1.16011900    1.99693900
  51. H                  1.47231900    4.30778900   -0.22668300
  52. H                  1.47376400    2.88224000   -1.28560200
  53. H                  0.03076900    3.88889800   -1.17684900

复制代码

BinWang 发表于 Post on 2026-1-14 13:03:43
看了一眼你的结构，初猜太烂，重新构建模型吧
BinWang 发表于 Post on 2026-1-14 12:54:25
难收敛的话，先减小基组，生成一个过渡态，之后将该过渡态作为初猜用大基组重新找过渡态。

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

GMT+8, 2026-1-24 04:22 , Processed in 0.174910 second(s), 25 queries , Gzip On.

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