Gaussian计算中出现的问题，停在某一步很久了

1848820738

在gaussian优化结构的时候停在一步很长时间，作业显示仍在计算，不知道是什么原因，还望各位老师帮忙解决，以下是log输出文件的最后一部分，chk文件已经两天没有更新了

Error on total polarization charges =  0.06900
SCF Done:  E(RM062X) =  -3278.68563893     A.U. after    1 cycles
            NFock=  1  Conv=0.77D-08     -V/T= 2.0029
SMD-CDS (non-electrostatic) energy       (kcal/mol) =     -20.47
(included in total energy above)
DoSCS=F DFT=T ScalE2(SS,OS)=  1.000000  1.000000
Range of M.O.s used for correlation:     1  1920
NBasis=  1941 NAE=   201 NBE=   201 NFC=     0 NFV=     0
NROrb=   1920 NOA=   201 NOB=   201 NVA=  1719 NVB=  1719

**** Warning!!: The largest alpha MO coefficient is  0.85445555D+02

PrsmSu:  requested number of processors reduced to:   3 ShMem   1 Linda.
Symmetrizing basis deriv contribution to polar:
IMax=3 JMax=2 DiffMx= 0.00D+00
G2DrvN: will do    18 centers at a time, making    6 passes.
Calling FoFCou, ICntrl=  3107 FMM=T I1Cent=   0 AccDes= 0.00D+00.
Calling FoFCou, ICntrl=  3107 FMM=T I1Cent=   0 AccDes= 0.00D+00.
Calling FoFCou, ICntrl=  3107 FMM=T I1Cent=   0 AccDes= 0.00D+00.
Calling FoFCou, ICntrl=  3107 FMM=T I1Cent=   0 AccDes= 0.00D+00.
Calling FoFCou, ICntrl=  3107 FMM=T I1Cent=   0 AccDes= 0.00D+00.
Calling FoFCou, ICntrl=  3107 FMM=T I1Cent=   0 AccDes= 0.00D+00.
NEqPCM:  Using equilibrium solvation (IEInf=0, Eps=   2.3741, EpsInf=   2.2383)
G2PCM: DoFxE=T DoFxN=T DoGrad=T DoDP/DQ/DG/TGxP=FFFF NFrqRd=   0 IEInf=0 SqF1=F DoCFld=F IF1Alg=4.
End of G2Drv F.D. properties file   721 does not exist.
End of G2Drv F.D. properties file   722 does not exist.
End of G2Drv F.D. properties file   788 does not exist.
          IDoAtm=11111111111111111111111111111111111111111111111111
          IDoAtm=11111111111111111111111111111111111111111111111111
          IDoAtm=11111
NEqPCM:  Using equilibrium solvation (IEInf=0, Eps=   2.3741, EpsInf=   2.2383)
          Differentiating once with respect to electric field.
                with respect to dipole field.
          Differentiating once with respect to nuclear coordinates.
CalDSu exits because no D1Ps are significant.
          There are   318 degrees of freedom in the 1st order CPHF.  IDoFFX=6 NUNeed=     3.
    318 vectors produced by pass  0 Test12= 3.46D-13 1.00D-09 XBig12= 5.27D+02 9.04D+00.
AX will form    30 AO Fock derivatives at one time.
    316 vectors produced by pass  1 Test12= 3.46D-13 1.00D-09 XBig12= 7.12D+01 1.13D+00.
    315 vectors produced by pass  2 Test12= 3.46D-13 1.00D-09 XBig12= 1.64D+00 1.18D-01.
    315 vectors produced by pass  3 Test12= 3.46D-13 1.00D-09 XBig12= 2.97D-02 1.23D-02.
    315 vectors produced by pass  4 Test12= 3.46D-13 1.00D-09 XBig12= 2.81D-04 1.14D-03.
    315 vectors produced by pass  5 Test12= 3.46D-13 1.00D-09 XBig12= 2.17D-06 8.33D-05.
    247 vectors produced by pass  6 Test12= 3.46D-13 1.00D-09 XBig12= 1.22D-08 4.70D-06.
    203 vectors produced by pass  7 Test12= 3.46D-13 1.00D-09 XBig12= 5.77D-11 3.57D-07.
    190 vectors produced by pass  8 Test12= 3.46D-13 1.00D-09 XBig12= 3.54D-13 3.09D-08.
    188 vectors produced by pass  9 Test12= 3.46D-13 1.00D-09 XBig12= 3.86D-13 1.82D-08.
    180 vectors produced by pass 10 Test12= 3.46D-13 1.00D-09 XBig12= 5.76D-13 2.33D-08.
    177 vectors produced by pass 11 Test12= 3.46D-13 1.00D-09 XBig12= 3.88D-13 1.77D-08.
    173 vectors produced by pass 12 Test12= 3.46D-13 1.00D-09 XBig12= 3.91D-13 1.46D-08.
    171 vectors produced by pass 13 Test12= 3.46D-13 1.00D-09 XBig12= 4.48D-13 1.61D-08.
    171 vectors produced by pass 14 Test12= 3.46D-13 1.00D-09 XBig12= 4.96D-13 1.74D-08.
    171 vectors produced by pass 15 Test12= 3.46D-13 1.00D-09 XBig12= 4.83D-13 1.60D-08.
    170 vectors produced by pass 16 Test12= 3.46D-13 1.00D-09 XBig12= 5.13D-13 1.55D-08.
    167 vectors produced by pass 17 Test12= 3.46D-13 1.00D-09 XBig12= 4.14D-13 1.28D-08.
    161 vectors produced by pass 18 Test12= 3.46D-13 1.00D-09 XBig12= 4.22D-13 1.54D-08.
    158 vectors produced by pass 19 Test12= 3.46D-13 1.00D-09 XBig12= 3.22D-13 1.27D-08.
    158 vectors produced by pass 20 Test12= 3.46D-13 1.00D-09 XBig12= 3.86D-13 1.28D-08.
    150 vectors produced by pass 21 Test12= 3.46D-13 1.00D-09 XBig12= 3.45D-13 1.60D-08.
    146 vectors produced by pass 22 Test12= 3.46D-13 1.00D-09 XBig12= 2.69D-13 1.15D-08.
    142 vectors produced by pass 23 Test12= 3.46D-13 1.00D-09 XBig12= 4.19D-13 1.31D-08.
    139 vectors produced by pass 24 Test12= 3.46D-13 1.00D-09 XBig12= 2.62D-13 1.23D-08.
    134 vectors produced by pass 25 Test12= 3.46D-13 1.00D-09 XBig12= 2.47D-13 1.05D-08.
    119 vectors produced by pass 26 Test12= 3.46D-13 1.00D-09 XBig12= 2.16D-13 1.08D-08.
    117 vectors produced by pass 27 Test12= 3.46D-13 1.00D-09 XBig12= 1.91D-13 8.42D-09.
    117 vectors produced by pass 28 Test12= 3.46D-13 1.00D-09 XBig12= 1.86D-13 9.97D-09.
    113 vectors produced by pass 29 Test12= 3.46D-13 1.00D-09 XBig12= 1.99D-13 8.78D-09.
    109 vectors produced by pass 30 Test12= 3.46D-13 1.00D-09 XBig12= 1.55D-13 9.85D-09.
    109 vectors produced by pass 31 Test12= 3.46D-13 1.00D-09 XBig12= 2.23D-13 1.10D-08.
    101 vectors produced by pass 32 Test12= 3.46D-13 1.00D-09 XBig12= 1.02D-13 6.55D-09.
    101 vectors produced by pass 33 Test12= 3.46D-13 1.00D-09 XBig12= 1.54D-13 7.94D-09.
    100 vectors produced by pass 34 Test12= 3.46D-13 1.00D-09 XBig12= 1.57D-13 8.09D-09.
    100 vectors produced by pass 35 Test12= 3.46D-13 1.00D-09 XBig12= 1.49D-13 7.78D-09.
    100 vectors produced by pass 36 Test12= 3.46D-13 1.00D-09 XBig12= 1.19D-13 6.92D-09.
    100 vectors produced by pass 37 Test12= 3.46D-13 1.00D-09 XBig12= 1.19D-13 6.20D-09.
    100 vectors produced by pass 38 Test12= 3.46D-13 1.00D-09 XBig12= 1.20D-13 6.67D-09.
     97 vectors produced by pass 39 Test12= 3.46D-13 1.00D-09 XBig12= 1.25D-13 6.84D-09.
     90 vectors produced by pass 40 Test12= 3.46D-13 1.00D-09 XBig12= 1.30D-13 6.97D-09.
     86 vectors produced by pass 41 Test12= 3.46D-13 1.00D-09 XBig12= 7.46D-14 4.89D-09.
     86 vectors produced by pass 42 Test12= 3.46D-13 1.00D-09 XBig12= 9.00D-14 6.43D-09.
     82 vectors produced by pass 43 Test12= 3.46D-13 1.00D-09 XBig12= 8.52D-14 5.66D-09.
     80 vectors produced by pass 44 Test12= 3.46D-13 1.00D-09 XBig12= 6.56D-14 5.14D-09.
     80 vectors produced by pass 45 Test12= 3.46D-13 1.00D-09 XBig12= 1.03D-13 5.96D-09.
     73 vectors produced by pass 46 Test12= 3.46D-13 1.00D-09 XBig12= 5.56D-14 4.44D-09.
     73 vectors produced by pass 47 Test12= 3.46D-13 1.00D-09 XBig12= 6.72D-14 5.41D-09.
     72 vectors produced by pass 48 Test12= 3.46D-13 1.00D-09 XBig12= 6.43D-14 5.00D-09.
     72 vectors produced by pass 49 Test12= 3.46D-13 1.00D-09 XBig12= 7.65D-14 5.44D-09.
     67 vectors produced by pass 50 Test12= 3.46D-13 1.00D-09 XBig12= 6.80D-14 4.41D-09.
     67 vectors produced by pass 51 Test12= 3.46D-13 1.00D-09 XBig12= 5.75D-14 4.58D-09.
     67 vectors produced by pass 52 Test12= 3.46D-13 1.00D-09 XBig12= 5.11D-14 4.19D-09.
     65 vectors produced by pass 53 Test12= 3.46D-13 1.00D-09 XBig12= 6.27D-14 4.77D-09.
     65 vectors produced by pass 54 Test12= 3.46D-13 1.00D-09 XBig12= 8.51D-14 6.23D-09.
     64 vectors produced by pass 55 Test12= 3.46D-13 1.00D-09 XBig12= 7.13D-14 5.03D-09.
     64 vectors produced by pass 56 Test12= 3.46D-13 1.00D-09 XBig12= 5.63D-14 4.29D-09.
     64 vectors produced by pass 57 Test12= 3.46D-13 1.00D-09 XBig12= 7.66D-14 4.91D-09.
     64 vectors produced by pass 58 Test12= 3.46D-13 1.00D-09 XBig12= 6.09D-14 4.45D-09.


因为log文件限制大小，无法上传，只能粘贴一部分

qeheikj

zjxitcc 发表于 2025-9-19 18:50
有解决办法，例如：找有g16的课题组合作；购买g16；或者查阅课题组/单位的Gaussian许可证，如果是Site Li ...

好的，谢谢您

zjxitcc

qeheikj 发表于 2025-9-19 17:33
好吧，组里没有g16。还有什么解决方法吗，或者加快速度的方法呢？

有解决办法，例如：找有g16的课题组合作；购买g16；或者查阅课题组/单位的Gaussian许可证，如果是Site Licence可以一个校区内大家共用。注意，换用其他支持TDDFT解析Hessian的商业量化软件（例如Q-Chem）不在本回答考虑之内，因为那不是解决办法，还不如购买和使用g16。注意，ORCA目前也不支持TDDFT解析Hessian。

没有“加快速度的方法”，Gaussian常见功能 是所有量化软件里最快、最稳健的（不考虑RI加速时，高斯软件几乎是世上最快的），要想更快，只有更新的Gaussian对比更旧的Gaussian，它自己打败过去的自己，例如使用g16可以进行解析Hessian计算。

qeheikj

zjxitcc 发表于 2025-9-19 17:25
g09没有TDDFT解析Hessian，现在是在以数值Hessian在做计算，要产生大量的微小位移结构，算到猴年马月。你 ...

好吧，组里没有g16。还有什么解决方法吗，或者加快速度的方法呢？

zjxitcc

qeheikj 发表于 2025-9-19 16:58
抱歉截图放错了，下面分别是输入与输出文件，麻烦您帮忙看一下

g09没有TDDFT解析Hessian，现在是在以数值Hessian在做计算，要产生大量的微小位移结构，算到猴年马月。你需要使用g16

qeheikj

zjxitcc 发表于 2025-9-19 16:09
你在20L与22L展示的两幅图看起来不是 一组严格对应的输入、输出文件，这很容易给回答者造成误解。当然， ...

抱歉截图放错了，下面分别是输入与输出文件，麻烦您帮忙看一下

zjxitcc

qeheikj 发表于 2025-9-19 16:04
确实，是想解决这个问题，做的是基态优化和频率计算，关键词如图

你在20L与22L展示的两幅图看起来不是 一组严格对应的输入、输出文件，这很容易给回答者造成误解。当然，我也无法回答（它压根不是一个合理的问题）。

qeheikj

zjxitcc 发表于 2025-9-19 10:10
你的问题与波函数稳定性没有直接关系，你可能更想问的是“为什么我这个激发态opt freq任务算了三四天还没 ...

确实，是想解决这个问题，做的是基态优化和频率计算，关键词如图

zjxitcc

qeheikj 发表于 2025-9-19 09:53
老师，怎么看体系是否具有波函数不稳定性呢？我也是opt和freq一起算的，泛函基组b3lyp/6-311g(d)。优化很 ...

你的问题与波函数稳定性没有直接关系，你可能更想问的是“为什么我这个激发态opt freq任务算了三四天还没算完”。举个例子，你问我怎么检验这袋子苹果有没有1斤重，但你实际上正在吃梨，二者风马牛不相及。

另外，如果做的是激发态计算（例如写了TD），应当同时交待激发态计算细节，而不是仅仅给出泛函和基组。当然，如果直接展示gjf文件关键词那一行会更清晰。

qeheikj

zjxitcc 发表于 2021-5-1 10:06
“freq一般两天出不来基本上就是出不来了” 并不是这样的，出现这种情况可能是CPHF难以收敛，或者体系过 ...

老师，怎么看体系是否具有波函数不稳定性呢？我也是opt和freq一起算的，泛函基组b3lyp/6-311g(d)。优化很早就停了，应该是一直在算频率，大概三四天了

zjxitcc

虎王 发表于 2023-4-4 11:34
那您何时会去确认波函数稳定性？

我天天都检验波函数稳定性，我算的体系的电子结构都比较复杂。

虎王

zjxitcc 发表于 2021-5-1 10:06
“freq一般两天出不来基本上就是出不来了” 并不是这样的，出现这种情况可能是CPHF难以收敛，或者体系过 ...

那您何时会去确认波函数稳定性？

1848820738

lw198 发表于 2021-5-1 04:54
这是在做frequency吧，如果长时间停在这里说明很大可能算不出来了。。我之前有在这个步骤持续算了九天还没 ...

已经算出来了

lw198

zjxitcc 发表于 2021-5-1 10:06
“freq一般两天出不来基本上就是出不来了” 并不是这样的，出现这种情况可能是CPHF难以收敛，或者体系过 ...

哦哦好的

zjxitcc

lw198 发表于 2021-5-1 04:54
这是在做frequency吧，如果长时间停在这里说明很大可能算不出来了。。我之前有在这个步骤持续算了九天还没 ...

“freq一般两天出不来基本上就是出不来了” 并不是这样的，出现这种情况可能是CPHF难以收敛，或者体系过大，又或者给的内存过少。高斯是很成熟的量化软件，一般是不会随便就挂掉的（特别是linux下就更稳妥了）。如果由于机器配置不高，实在无法再增大内存，在合适的内存下多等几天总是能算出来的。另外，如果你之前算过的例子有波函数不稳定性，那基于它做freq无意义，自然也就不能拿这种例子来说明freq耗时长短。