计算化学公社

标题: Gaussian用优化后的结构进行单点能计算时出现502错误合理吗 [打印本页]

作者
Author: abu123 时间: 2023-10-19 15:20
标题: Gaussian用优化后的结构进行单点能计算时出现502错误合理吗
刚接触gaussian计算不久，计算过程中存在一些疑问，希望各位老师能答疑
几何优化收敛，然后用优化后的结构进行单点能计算时出现502错误，显示收敛失败，这种情况合理吗？

另外，几何优化后计算所得的结合能E1与用优化后的结构进行单点能计算(选择更高精度的方法与基组)，然后算的结合能差距较大？这是什么原因呢？合理吗？

作者
Author: 大喵说多谢 时间: 2023-10-19 19:26
http://bbs.keinsci.com/thread-4829-1-1.html
L502的错在这里，你对照输出文件看看哪里出问题了。

作者
Author: abu123 时间: 2023-10-19 20:35

大喵说多谢发表于 2023-10-19 19:26
http://bbs.keinsci.com/thread-4829-1-1.html
L502的错在这里，你对照输出文件看看哪里出问题了。

我是在几何优化后的结构上用ccsd,def2tzvpp计算单点能，然后有很大部分结构显示SCF收敛失败，还有一个是什么错误我没有找出来。错误如下：
DD1Dir will call FoFJK  11 times, MxPair=    914
NAB=  3364 NAA=  1653 NBB=  1653 NumPrc=  4.
FoFJK:  IHMeth= 1 ICntrl=    200 DoSepK=F KAlg= 0 I1Cent= 0 FoldK=F
IRaf= 990000000 NMat= 914 IRICut= 1142 DoRegI=T DoRafI=T ISym2E=-1.
FoFCou: FMM=F IPFlag=          0 FMFlag=    100000 FMFlg1=          0
      NFxFlg=          0 DoJE=F BraDBF=F KetDBF=F FulRan=T
      wScrn=  0.000000 ICntrl=    200 IOpCl=  0 I1Cent=          0 NGrid=          0
      NMat0=  914 NMatS0=    0 NMatT0=  457 NMatD0=  914 NMtDS0= 0 NMtDT0= 0
Petite list used in FoFCou.
这种情况是因为优化的结构存在问题吗？

作者
Author: sobereva 时间: 2023-10-20 07:59

abu123 发表于 2023-10-19 20:35
我是在几何优化后的结构上用ccsd,def2tzvpp计算单点能，然后有很大部分结构显示SCF收敛失败，还有一个是 ...

没事甭用CCSD，毫无意义，精度明显比不上CCSD(T)，耗时也没比CCSD(T)便宜哪去。顶级的双杂化泛函的精度甚至有时都能超过贵得多的CCSD

SCF不收敛该怎么解决怎么解决
解决SCF不收敛问题的方法
http://sobereva.com/61

另外，几何优化后计算所得的结合能E1与用优化后的结构进行单点能计算(选择更高精度的方法与基组)，然后算的结合能差距较大？这是什么原因呢？合理吗？

谁知道你具体用的什么级别。最关键信息不说，问这种问题没意义

作者
Author: abu123 时间: 2023-10-20 11:11

sobereva 发表于 2023-10-20 07:59
没事甭用CCSD，毫无意义，精度明显比不上CCSD(T)，耗时也没比CCSD(T)便宜哪去。顶级的双杂化泛函的精度甚 ...

我是用B3LYP,DEF2TZVP对Ti4Ta4C12进行了几何优化，计算了体系的结合能,之后用MP2,DEF2TZVPP再优化完的结构上进行了单点能的计算，最后算了结合能，然后发现两个结合能相差大概30ev

作者
Author: wzkchem5 时间: 2023-10-20 23:23

abu123 发表于 2023-10-20 04:11
我是用B3LYP,DEF2TZVP对Ti4Ta4C12进行了几何优化，计算了体系的结合能,之后用MP2,DEF2TZVPP再优化完的结 ...

这个体系的HOMO-LUMO gap多少？如果特别小的话，MP2微扰发散，MP2结果不能用。

作者
Author: sobereva 时间: 2023-10-21 02:55

abu123 发表于 2023-10-20 11:11
我是用B3LYP,DEF2TZVP对Ti4Ta4C12进行了几何优化，计算了体系的结合能,之后用MP2,DEF2TZVPP再优化完的结 ...

没事甭用MP2，尤其是算涉及过渡金属的体系
比MP2便宜得多的大量普通泛函比MP2结果还更好。这年头如果没有特殊理由绝对不要用MP2
仔细看
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272（http://bbs.keinsci.com/thread-536-1-1.html）

作者
Author: abu123 时间: 2023-10-21 16:14

wzkchem5 发表于 2023-10-20 23:23
这个体系的HOMO-LUMO gap多少？如果特别小的话，MP2微扰发散，MP2结果不能用。

1.19ev左右

作者
Author: abu123 时间: 2023-10-21 16:14

sobereva 发表于 2023-10-21 02:55
没事甭用MP2，尤其是算涉及过渡金属的体系
比MP2便宜得多的大量普通泛函比MP2结果还更好。这年头如果没 ...

好的，谢谢老师

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