计算Eu@C78（富勒烯C78内嵌镧系金属铕），结构优化l502报错

gqy · 发表于 Post on 2022-4-28 15:38:43

本帖最后由 gqy 于 2022-4-28 15:46 编辑

最近在进行Eu@C78的结构优化时，尝试了各种办法（同时参考了老师的文章：解决SCF不收敛问题的方法），但是在输出文件L502模块，一直显示收敛失败。
我们尝试了不同的泛函：PBE1PBE,B3LYP,M062X,TPSSH，C原子用的基组分别有6-31g*和甚至STO-3G，金属Eu一直使用的SDD赝势，SDD赝势基组，在上面各种泛函和基组的组合下，尝试了不同的关键词进行计算。前面计算用的g16，后面同时开始使用g09，但是输出文件一直显示收敛失败。

计算Eu@C78的输入输出的部分文件.doc (65.5 KB, 下载次数 Times of downloads: 8) 文件中附带了计算的输入文件关键词和输出文件的报错部分。图片显示了投入计算的分子结构，还有计算的其中一种方法的输入和输出文件。问题似乎有点穷途末路的感觉。求大神指点迷津，万分感谢。

wzkchem5 · 发表于 Post on 2022-4-28 15:55:09

自旋多重度设错了，Eu那么多f电子，自旋多重度不可能是2

gqy · 发表于 Post on 2022-4-28 16:25:45

wzkchem5 发表于 2022-4-28 15:55
自旋多重度设错了，Eu那么多f电子，自旋多重度不可能是2

天呐，惊醒。万分感谢您的提示，这就改试试。

sobereva · 发表于 Post on 2022-4-29 09:18:26

失败前最后一轮迭代的能量变化已经很小了，可以用scf=conver=6先做预优化，等结构收敛得差不多了，再换回默认SCF收敛限。这在《使用Gaussian做镧系金属配合物的量子化学计算》（http://sobereva.com/581）里都说过。

gqy · 发表于 Post on 2022-4-29 09:38:33

sobereva 发表于 2022-4-29 09:18
失败前最后一轮迭代的能量变化已经很小了，可以用scf=conver=6先做预优化，等结构收敛得差不多了，再换回默 ...

万分感谢您的回复。conver=6的方法我们也尝试过，但是在L701报错。

输入文件9
#P PBE1PBE/genecp int=fine guess=read opt freq SCF=(vshift=450,maxcyc=180,conver=6)
C 0
6-31G*
****
Eu 0
SDD
****
Eu 0
SDD
输出log文件9 Link 0=g16
No NMR shielding tensors so no spin-rotation constants.
Leave Link 601 at Wed Apr 27 11:37:01 2022,  MaxMem= 6710886400  cpu:174.7 elap:    6.7
(Enter /data/software/calculate/g16/l701.exe)
SCFChk:  SCF convergence  4.53D-08 required  1.00D-08

SCF Error SCF Error SCF Error SCF Error SCF Error SCF Error SCF Error SCF Error

                                 ERROR!!!!

SCF has not converged.  Gradients and post-SCF results would be GARBAGE!!

SCF Error SCF Error SCF Error SCF Error SCF Error SCF Error SCF Error SCF Error

Error termination via Lnk1e in /data/software/calculate/g16/l701.exe at Wed Apr 27 11:37:01 2022.

ionexchangeC · 发表于 Post on 2022-4-29 14:26:27

总之先把自旋多重度改成8，还是不收敛再试别的方法。不要忘记检查波函数稳定性。

scf · 发表于 Post on 2022-4-29 16:02:45

用scf=xqc试试。另外f电子不大容易成键，Cr2的基态是单重态，但Gd2的基态自旋多重度很高

gqy · 发表于 Post on 2022-4-29 20:41:35

ionexchangeC 发表于 2022-4-29 14:26
总之先把自旋多重度改成8，还是不收敛再试别的方法。不要忘记检查波函数稳定性。

好的，感谢您的回复

gqy · 发表于 Post on 2022-4-29 20:47:58

scf 发表于 2022-4-29 16:02
用scf=xqc试试。另外f电子不大容易成键，Cr2的基态是单重态，但Gd2的基态自旋多重度很高

感谢您的提醒。在开始计算发现不能收敛时，我们也试过scf=QC，当时自旋多重度设置的2，输出文件依旧报错。目前正在尝试改变自旋多重度重新投入计算，希望能得出想要的结果。

sobereva · 发表于 Post on 2022-4-30 07:21:45

gqy 发表于 2022-4-29 09:38
万分感谢您的回复。conver=6的方法我们也尝试过，但是在L701报错。

输入文件9

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

文中明确说了

注2：如果你用G16做几何优化、振动分析等任务，或使用TDDFT、双杂化、后HF等方法，放宽SCF收敛限的时候必须加上特殊的IOp，否则会因为通不过收敛精度检测而遇到报错，参见《令Gaussian 16中SCF未收敛到默认收敛限也能继续做后续计算的方法》（http://sobereva.com/625）。

gqy · 发表于 Post on 2022-5-1 16:22:21

sobereva 发表于 2022-4-30 07:21
解决SCF不收敛问题的方法
http://sobereva.com/61

好的，咱再加上IOp试试。感谢您的提示。

七尺贱 · 发表于 Post on 2022-5-2 10:43:57

symm=loose，scf=conver=6

gqy · 发表于 Post on 2022-5-15 15:18:28

最近终于完成了结构，和单点的计算。感谢以上各位老师的建议。

zjxitcc · 发表于 Post on 2022-5-15 15:38:41

七尺贱发表于 2022-5-2 10:43
symm=loose，scf=conver=6

scf=conver=6只能当做粗优化，后续仍应收敛至正常精度

halouhapily · 发表于 Post on 2022-9-14 17:39:53

wzkchem5 发表于 2022-4-28 15:55
自旋多重度设错了，Eu那么多f电子，自旋多重度不可能是2

老师，关于Eu(lll）体系的自旋多重度的设置（配体都是闭壳层结构），您看我理解的对吗？
Eu(lll)4f电子有6个，如果用大核赝势的话，4f电子被赝势化了，因此可以认为自旋多重度为1。如果用小核赝势的话，基态的自旋多重度需要测试，这样理解对吗

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] 计算Eu@C78（富勒烯C78内嵌镧系金属铕），结构优化l502报错