请问win版Gaussian16模拟报错，但输出文件中未见Err怎么办？

hhann

使用win版Gaussian16进行DFT模拟，都是跑到同一位置后出现Link died，但输出文件中未见Err。大致在跑L1110处错误，求助请帮忙看看是什么问题？后续又用不同电脑用G09跑相同文件，出错位置与G16版不同。输入文件及输出文件均已放至附件。
以下为输入文件部分
%nprocshared=8
%mem=40GB
%chk=D:\sihangaussian\EDOT0204.chk
# opt freq b3lyp/6-311g(d,p) em=gd3bj


EDOT0204


0 1
C                 -9.25860000   -1.54470000   11.45980000
C                 -9.61950000   -2.32630000   10.41650000
C                 -8.81720000   -2.60940000    9.36510000
C                 -7.54300000   -2.18440000    9.48870000
C                 -7.18260000   -1.39610000   10.52700000
C                 -8.03000000   -0.98900000   11.49960000
C                 -6.44060000   -2.39000000    8.74630000
C                 -5.42020000   -1.71170000    9.29780000
S                 -5.76670000   -0.97410000   10.49900000
C                -10.30550000   -1.50190000   12.30360000
C                -11.31350000   -2.21350000   11.77090000
S                -11.01140000   -2.81600000   10.48590000
N                 -5.94220000   -2.96680000    7.72260000
C                 -4.70970000   -2.70560000    7.57480000
C                 -4.32240000   -1.89900000    8.56030000
C                 -3.67860000   -2.95630000    6.77040000
C                 -2.60320000   -2.30910000    7.26560000
S                 -2.91650000   -1.52120000    8.48360000
O                 -7.60510000   -0.20830000   12.54410000
O                 -9.22470000   -3.43810000    8.35100000
C                 -7.57970000    1.17360000   12.27070000
C                 -6.91360000    1.92330000   13.42950000
C                 -9.97080000   -2.80360000    7.33800000
C                -10.52850000   -3.85000000    6.36700000
N                -10.75540000   -1.06720000   13.41670000
C                -11.95220000   -1.43230000   13.62790000
C                -12.35960000   -2.16860000   12.59940000
C                -12.93840000   -1.33990000   14.51800000
C                -14.01960000   -2.02010000   14.06610000
S                -13.73020000   -2.63460000   12.75260000
C                -15.15090000   -2.05760000   14.80460000
C                -16.39890000   -2.57970000   14.73830000
C                -17.37020000   -2.44140000   15.69130000
C                -18.45890000   -3.10730000   15.24220000
C                -18.16490000   -3.63650000   14.05170000
C                -16.89390000   -3.32850000   13.70800000
C                -10.08310000   -0.24950000   14.43720000
C                 -6.58550000   -3.86220000    6.74970000
C                 -6.35400000   -5.35320000    7.06430000
C                 -6.97830000   -5.78310000    8.40090000
C                 -6.88620000   -6.24090000    5.92620000
C                 -9.48780000   -1.08750000   15.58530000
C                 -8.38060000   -2.04290000   15.11370000
C                 -8.95370000   -0.17640000   16.70400000
C                -11.46470000   -4.84010000    7.08020000
C                -11.26500000   -3.16000000    5.20650000
C                 -5.46120000    1.46290000   13.63950000
C                 -6.96140000    3.44010000   13.17920000
C                -19.69600000   -3.32970000   15.71900000
C                -20.58960000   -4.05340000   15.02000000
C                -20.26120000   -4.57300000   13.82160000
C                -19.03070000   -4.36230000   13.32880000



以下为输出文件最后
Rotational constants (GHZ):           0.0183201           0.0007353           0.0007242
Standard basis: 6-311G(d,p) (5D, 7F)
There are  3695 symmetry adapted cartesian basis functions of A   symmetry.
There are  3534 symmetry adapted basis functions of A   symmetry.
  3534 basis functions,  6030 primitive gaussians,  3695 cartesian basis functions
   597 alpha electrons      597 beta electrons
       nuclear repulsion energy     34385.9302406618 Hartrees.
NAtoms=  255 NActive=  255 NUniq=  255 SFac= 1.00D+00 NAtFMM=   60 NAOKFM=T Big=T
Integral buffers will be    131072 words long.
Raffenetti 2 integral format.
Two-electron integral symmetry is turned on.
Nuclear repulsion after empirical dispersion term =    34385.1920048234 Hartrees.
One-electron integrals computed using PRISM.
NBasis=  3534 RedAO= T EigKep=  7.15D-05  NBF=  3534
NBsUse=  3534 1.00D-06 EigRej= -1.00D+00 NBFU=  3534
Initial guess from the checkpoint file:  "D:\sihangaussian\EDOT0204.chk"
B after Tr=     0.000000    0.000000    0.000000
         Rot=    1.000000    0.000000    0.000000    0.000000 Ang=   0.00 deg.
Requested convergence on RMS density matrix=1.00D-08 within 128 cycles.
Requested convergence on MAX density matrix=1.00D-06.
Requested convergence on             energy=1.00D-06.
No special actions if energy rises.
SCF Done:  E(RB3LYP) =  -10252.6477631     A.U. after    1 cycles
            NFock=  1  Conv=0.41D-08     -V/T= 2.0030
DoSCS=F DFT=T ScalE2(SS,OS)=  1.000000  1.000000
Range of M.O.s used for correlation:     1  3534
NBasis=  3534 NAE=   597 NBE=   597 NFC=     0 NFV=     0
NROrb=   3534 NOA=   597 NOB=   597 NVA=  2937 NVB=  2937


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




**** Warning!!: The smallest alpha delta epsilon is  0.70748254D-01


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

Hugo_314cat

简单看了一下，真是好大的分子啊。输入文件没什么问题，只是你给的配置真的够吗？

hhann

Hugo_314cat 发表于 2025-2-10 11:21
简单看了一下，真是好大的分子啊。输入文件没什么问题，只是你给的配置真的够吗？

是的，大约有250多个原子。我的电脑是i9，16核。曾试过给到16核，但不久后就死机了，不知是否是给的太大的原因。曾咨询过别人用win版G09跑相似大小分子，大约跑了1个多月也是可以跑出来的。但我这个跑几天就link died

Hugo_314cat

hhann 发表于 2025-2-10 11:35
是的，大约有250多个原子。我的电脑是i9，16核。曾试过给到16核，但不久后就死机了，不知是否是给的太大 ...

搞个服务器吧，个人认为这种用小电脑+win跑纯纯折磨人。就算是花了一个多月他能跑出来，我觉得也很不值得。

lemon_electron

自己电脑跑这个的话基本也干不了别的事情了吧， 搞个服务器吧

Santz

报错是l502.exe Convergence failure -- run terminated.

hhann

Hugo_314cat 发表于 2025-2-10 12:01
搞个服务器吧，个人认为这种用小电脑+win跑纯纯折磨人。就算是花了一个多月他能跑出来，我觉得也很不值得 ...

好的，感谢您的意见。因为条件原因可能需要外包，所以想先排除一下输入文件是否无误，以防浪费钱财

savorw

250多个原子用个中等基组6-31G(d)跑就差不多了。还是windows机器，还用6-311g(d,p) 的基组。

hhann

Santz 发表于 2025-2-10 14:05
报错是l502.exe Convergence failure -- run terminated.

谢谢您的帮助。还想冒昧问一下，这是因为不收敛的原因吗？是我基组选的不对还是结构问题呢？我在Gaussian计算前用chem3d的MM2进行了一下结构优化，不知是否与此相关？

sobereva

认真看
解决SCF不收敛问题的方法
http://sobereva.com/61

200原子以上的DFT的opt freq，没个像样的双路服务器根本甭考虑。拿个人计算机计算纯属瞎折腾
花费2、3天的任务已经属于很高耗时的任务了，通常来说需要更多耗时要么说明机子太烂要么说明计算模型严重不合理。一个任务算一个月纯属胡搞瞎搞。windows版Gaussian算这个完全是开玩笑