计算化学公社

标题: 计算含C N O H S F原子的分子做优化+振动分析计算，出现Gaussian L1002 error [打印本页]

作者
Author: Shawn_ 时间: yesterday 17:18
标题: 计算含C N O H S F原子的分子做优化+振动分析计算，出现Gaussian L1002 error
本帖最后由 Shawn_ 于 2025-8-14 21:01 编辑

Gaussian16 C01 计算含C N O H S F原子的分子做优化+振动分析计算【关键字：# opt freq b3lyp/6-311g(d,p) em=gd3bj，出现L1002错误，请问该如何解决呢？

******************************************
%nprocshared=8
Will use up to 8 processors via shared memory.
%mem=30GB
%chk=/home/shawn/workspace/Zh-F\zh-f.chk
-------------------------------------
...
**** Warning!!: The smallest alpha delta epsilon is  0.69219171D-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
G2DrvN: will do 29 centers at a time, making 8 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.
Calling FoFCou, ICntrl=  3107 FMM=T I1Cent= 0 AccDes= 0.00D+00.
Calling FoFCou, ICntrl=  3107 FMM=T I1Cent= 0 AccDes= 0.00D+00.
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=11111111111111111111111111111111111111111111111111
      IDoAtm=11111111111111111111111111111111111111111111111111
      IDoAtm=1111111111111111111111111
      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 678 degrees of freedom in the 1st order CPHF.  IDoFFX=6 NUNeed=    3.
678 vectors produced by pass  0 Test12= 1.11D-12 1.00D-09 XBig12= 9.11D+03 8.42D+01.
AX will form 48 AO Fock derivatives at one time.
OrtVc1:  Ph=1 IOff=    0 IPass=80 DotMx1= 1.68+278
OrtVc1:  Ph=1 M=  1107216 NPass=80 Test1= 1.11D-12 Small= 1.11D-06 VSmall= 1.00D-12
OrtVc1 failed #1.
Error termination via Lnk1e in /home/shawn/gaussian/g16/l1002.exe at Thu Aug 14 15:11:07 2025.
Job cpu time:    3 days 11 hours 41 minutes  1.4 seconds.
Elapsed time:    0 days 10 hours 28 minutes 40.1 seconds.
File lengths (MBytes):  RWF= 150996 Int=    0 D2E=    0 Chk= 594 Scr=    1
另外其中的cpu time与实际时间不符，实际程序运行时间大约为一天半，但程序最后的cpu time显示为三天，请问这又是什么原因呢？
谢谢各位老师。

作者
Author: sobereva 时间: yesterday 20:25
这叫优化+振动分析计算，不叫优化计算，描述必须准确，你的问题完全是振动分析步骤出现的

225原子体系指望用8核做这么贵的任务，完全是开玩笑。至少去找个>30核、不少于200 GB可用内存的机子跑

CPU Time是所有核跑当前任务耗时的总和。

作者
Author: Shawn_ 时间: yesterday 21:06

sobereva 发表于 2025-8-14 20:25
这叫优化+振动分析计算，不叫优化计算，描述必须准确，你的问题完全是振动分析步骤出现的

225原子体系指 ...

非常感谢sob老师的答疑，已修改帖子内容。我的机器是16核 96G，因为用同样的配置计算完一个200原子体系（相同），耗时2天，觉得可以接受，所以才用相同的配置计算这个体系。因为组里条件不足所以才尝试使用个人电脑计算。

作者
Author: sobereva 时间: yesterday 22:07

Shawn_ 发表于 2025-8-14 21:06
非常感谢sob老师的答疑，已修改帖子内容。我的机器是16核 96G，因为用同样的配置计算完一个200原子体系（ ...

本来此任务就特别耗时和耗内存，有多少资源就全都用上，给16核、90GB，还不行就换机子，或者基组降到6-31G*

作者
Author: Shawn_ 时间: yesterday 23:39

sobereva 发表于 2025-8-14 22:07
本来此任务就特别耗时和耗内存，有多少资源就全都用上，给16核、90GB，还不行就换机子，或者基组降到6-31 ...

谢谢sob老师，我再尝试一次

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