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

lw198

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

哦哦好的

1848820738

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

已经算出来了

虎王

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

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

zjxitcc

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

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

qeheikj

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

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

zjxitcc

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

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

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

qeheikj

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

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

zjxitcc

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

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

qeheikj

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

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

zjxitcc

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

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

qeheikj

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

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

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 18:50
有解决办法，例如：找有g16的课题组合作；购买g16；或者查阅课题组/单位的Gaussian许可证，如果是Site Li ...

好的，谢谢您

caiyan

各位老师我也想问一下，我也是在高斯结构优化的时候使用的是# opt freq b3lyp/6-311g(d,p)方法；显示这个报错 DoSCS=F DFT=T ScalE2(SS,OS)=  1.000000  1.000000
Range of M.O.s used for correlation:     1   906
NBasis=   906 NAE=   156 NBE=   156 NFC=     0 NFV=     0
NROrb=    906 NOA=   156 NOB=   156 NVA=   750 NVB=   750

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

PrsmSu:  requested number of processors reduced to:  27 ShMem   1 Linda.
PrsmSu:  requested number of processors reduced to:   1 ShMem   1 Linda.
Symmetrizing basis deriv contribution to polar:
IMax=3 JMax=2 DiffMx= 0.00D+00
G2DrvN: will do     1 centers at a time, making   72 passes.
Estimated number of processors is:   36
Calling FoFCou, ICntrl=  3107 FMM=T I1Cent=   0 AccDes= 0.00D+00.
Out-of-memory error in routine FMMEnt-BalFmn (IEnd=     191447358 MxCore=     191436081)
Use %mem=225MW to provide the minimum amount of memory required to complete this step.
Error termination via Lnk1e in /home/use/soft/g16/l1110.exe at Wed Apr 22 11:49:30 2026.
Job cpu time:       0 days  0 hours 16 minutes 46.0 seconds.
Elapsed time:       0 days  0 hours  0 minutes 21.3 seconds.
File lengths (MBytes):  RWF=   1100 Int=      0 D2E=      0 Chk=     62 Scr=      1
是什么原因呢？