关于G16中G4方法bug的一个解决办法

我本是个娃娃

关于G16中G4方法bug的一个解决办法
文/我本是个娃娃

摘要：在G16(A02, B01)中，使用G4方法计算存在不收敛的情况，使得整体计算结果错误。而G09中，又无法实现G16独有的G4(SP)和G4(noopt)功能。据此，本文提出一个解决这个bug的方案。经过测试表明，方案不仅可以解决这个问题，并使得最终的计算结果符合标准，并且可以灵活修改，完成Gaussian程序内置G4方法所无法完成的计算任务。

G4.md (7.58 KB, 下载次数 Times of downloads: 43)

2022-11-13 13:10 上传 Uploaded

点击下载
Click to download

G4.zip (68 KB, 下载次数 Times of downloads: 45)

2022-11-13 13:10 上传 Uploaded

点击下载
Click to download

mfdsrax2

为什么一篇文章只有标题作者和摘要，没有正文

冰释之川

mfdsrax2 发表于 2022-11-14 09:06
为什么一篇文章只有标题作者和摘要，没有正文

正文当然是最关键的脚本与安装包啦

糖糖糖豆9988

请教一下老师，我利用“ #p ub3lyp 6-31G(2df,p) guess=mix stable=opt nosymm”和“ #p G4=SP guess=read geom=allcheck nosymm”计算双自由基的单点能，out文件末尾如图，能量结果是异常的。out文件中只有最后一个SCF的<S**2>= 0.0000，其余都是接近1的，这是不是您提到的“G4方法计算存在不收敛的情况，使得整体计算结果错误”的原因呢？或者说明G4不适合计算这个双自由基的单点能呢？谢谢老师

我本是个娃娃

糖糖糖豆9988 发表于 2024-4-11 09:23
请教一下老师，我利用“ #p ub3lyp 6-31G(2df,p) guess=mix stable=opt nosymm”和“ #p G4=SP guess=read  ...

软件是用的G09还是G16？
如果是G16，建议改成G09试试。

阿豪0605

您好楼主，请问怎么看这一项DE(HF)=-29532.446718 # 注意这一项是否出错，
或是怎么查看G4的输出文件是有问题
我用文献的坐标和G4的方法来验证能垒，但是过渡态的能量较高，得出的正向和逆向能垒
都是负的，输出文件也显示Normal termination，下面是C6H12G4输出的能量值
Temperature=              298.150000 Pressure=                      1.000000
E(ZPE)=                     0.176579 E(Thermal)=                    0.187662
E(QCISD(T))=             -439.665635 E(Empiric)=                   -0.187769
E(CCSD(T))=                -0.187769 E(Empiric)=                   -0.025401
DE(Plus)=                  -0.368915 DE(2DF)=                      -0.573841
E(Delta-G3XP)=             -0.046009 DE(HF)=                     -440.690991
G4(0 K)=                 -440.679908 G4 Energy=                  -440.678963
G4 Enthalpy=             -440.734306 G4 Free Energy=

我本是个娃娃

阿豪0605 发表于 2026-1-13 16:19
您好楼主，请问怎么看这一项DE(HF)=-29532.446718 # 注意这一项是否出错，
或是怎么查看G4的输出文件是有 ...

1. 尝试使用G16的其它版本，或者用G09；
2. G4方法默认是优化到极小点，所以要添加opt=(ts, noeigen, calcfc)

阿豪0605

我本是个娃娃 发表于 2026-1-15 19:54
1. 尝试使用G16的其它版本，或者用G09；
2. G4方法默认是优化到极小点，所以要添加opt=(ts, noeigen, ca ...

你好楼主，下面是我用关键词# G4=SP，只想用G4求单点能，最后一行这个值G4=-235.1646122是吗？
DE(HF)=-235.164612这个值和好像，不加SP关键词，也是最后一行G4和DE(HF)一样，这是G16  A03版本的问题嘛
E(QCISD(T))=             -234.432148 E(Empiric)=                   -0.123617
E(CCSD(T))=                -0.123617 E(Empiric)=                   -0.011386
DE(Plus)=                  -0.226330 DE(2DF)=                      -0.347862
E(Delta-G3XP)=             -0.023269 DE(HF)=                     -235.164612
G4 Energy=      
1\1\GINC-C08R1N24\Mixed\G4\G4\C6H11(2)\LIUZHAO0605\17-Jan-2026\0\\# G4
=SP\\SPECIES NAME/COMMENTS: C 20/05/2015 17:53:49\\0,2\C,0,0.048717,-0
.748207,1.257251\C,0,0.048717,0.786799,1.278344\C,0,0.56593,1.349416,0
.\C,0,0.048717,0.786799,-1.278344\C,0,0.048717,-0.748207,-1.257251\C,0
,-0.643498,-1.27126,0.\H,0,-0.99587,1.121557,1.434002\H,0,0.619949,1.1
64809,2.133589\H,0,1.087208,-1.10787,1.270246\H,0,-0.438731,-1.136955,
2.159114\H,0,-0.99587,1.121557,-1.434002\H,0,0.619949,1.164809,-2.1335
89\H,0,-0.438731,-1.136955,-2.159114\H,0,1.087208,-1.10787,-1.270246\H
,0,-1.693246,-0.939144,0.\H,0,-0.658274,-2.367612,0.\H,0,1.1026,2.2916
35,0.\\Version=ES64L-G16RevA.03\State=2-A'\MP2/GTBas1=-234.336055\MP4/
GTBas1=-234.4288372\CCSD(T)/G3Bas1=-234.4321478\MP2/GTBas2=-234.346435
\MP4/GTBas2=-234.440223\MP2/GTBas3=-234.5494768\MP4/GTBas3=-234.655166
8\HF/GTLargeXP=-233.6519002\MP2/GTLargeXP=-234.9077191\HF/GFHFB1=-233.
6704581\HF/GFHFB2=-233.6742465\G4=-235.1646122