计算化学公社

标题: 关于G16中G4方法bug的一个解决办法 [打印本页]

作者
Author: 我本是个娃娃 时间: 2022-11-13 13:10
标题: 关于G16中G4方法bug的一个解决办法
关于G16中G4方法bug的一个解决办法
文/我本是个娃娃

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

(, 下载次数 Times of downloads: 43)

(, 下载次数 Times of downloads: 45)

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

作者
Author: 冰释之川 时间: 2022-11-14 09:57

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

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

作者
Author: 糖糖糖豆9988 时间: 2024-4-11 09:23
请教一下老师，我利用“ #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不适合计算这个双自由基的单点能呢？谢谢老师

作者
Author: 我本是个娃娃 时间: 2024-4-15 10:56

糖糖糖豆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试试。

作者
Author: 阿豪0605 时间: 2026-1-13 16:19
您好楼主，请问怎么看这一项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=

作者
Author: 我本是个娃娃 时间: 2026-1-15 19:54

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

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

作者
Author: 阿豪0605 时间: 6 day ago

我本是个娃娃发表于 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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)