计算化学公社

标题: 请问win版Gaussian16模拟报错，但输出文件中未见Err怎么办？ [打印本页]

作者
Author: hhann 时间: 2025-2-10 11:11
标题: 请问win版Gaussian16模拟报错，但输出文件中未见Err怎么办？
使用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

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

作者
Author: hhann 时间: 2025-2-10 11:35

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

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

作者
Author: Hugo_314cat 时间: 2025-2-10 12:01

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

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

作者
Author: lemon_electron 时间: 2025-2-10 13:34
自己电脑跑这个的话基本也干不了别的事情了吧，

搞个服务器吧

作者
Author: Santz 时间: 2025-2-10 14:05
报错是l502.exe Convergence failure -- run terminated. (, 下载次数 Times of downloads: 5)

作者
Author: hhann 时间: 2025-2-10 14:41

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

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

作者
Author: savorw 时间: 2025-2-10 14:48
250多个原子用个中等基组6-31G(d)跑就差不多了。还是windows机器，还用6-311g(d,p) 的基组。

作者
Author: hhann 时间: 2025-2-10 14:50

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

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

作者
Author: sobereva 时间: 2025-2-11 00:12
认真看
解决SCF不收敛问题的方法
http://sobereva.com/61

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

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