计算化学公社

标题: 计算Eu@C78（富勒烯C78内嵌镧系金属铕），结构优化l502报错 [打印本页]

作者
Author: gqy 时间: 2022-4-28 15:38
标题: 计算Eu@C78（富勒烯C78内嵌镧系金属铕），结构优化l502报错
本帖最后由 gqy 于 2022-4-28 15:46 编辑

最近在进行Eu@C78的结构优化时，尝试了各种办法（同时参考了老师的文章：解决SCF不收敛问题的方法），但是在输出文件L502模块，一直显示收敛失败。
我们尝试了不同的泛函：PBE1PBE,B3LYP,M062X,TPSSH，C原子用的基组分别有6-31g*和甚至STO-3G，金属Eu一直使用的SDD赝势，SDD赝势基组，在上面各种泛函和基组的组合下，尝试了不同的关键词进行计算。前面计算用的g16，后面同时开始使用g09，但是输出文件一直显示收敛失败。 (, 下载次数 Times of downloads: 8) 文件中附带了计算的输入文件关键词和输出文件的报错部分。图片显示了投入计算的分子结构，还有计算的其中一种方法的输入和输出文件。问题似乎有点穷途末路的感觉。求大神指点迷津，万分感谢。 (, 下载次数 Times of downloads: 49)

作者
Author: wzkchem5 时间: 2022-4-28 15:55
自旋多重度设错了，Eu那么多f电子，自旋多重度不可能是2

作者
Author: gqy 时间: 2022-4-28 16:25

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

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

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

作者
Author: gqy 时间: 2022-4-29 09:38

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.

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

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

作者
Author: gqy 时间: 2022-4-29 20:41

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

好的，感谢您的回复

作者
Author: gqy 时间: 2022-4-29 20:47

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

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

作者
Author: sobereva 时间: 2022-4-30 07:21

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）。

作者
Author: gqy 时间: 2022-5-1 16:22

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

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

作者
Author: 七尺贱 时间: 2022-5-2 10:43
symm=loose，scf=conver=6

作者
Author: gqy 时间: 2022-5-15 15:18
最近终于完成了结构，和单点的计算。感谢以上各位老师的建议。

作者
Author: zjxitcc 时间: 2022-5-15 15:38

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

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

作者
Author: halouhapily 时间: 2022-9-14 17:39

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

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

作者
Author: wzkchem5 时间: 2022-9-14 17:46

halouhapily 发表于 2022-9-14 10:39
老师，关于Eu(lll）体系的自旋多重度的设置（配体都是闭壳层结构），您看我理解的对吗？
Eu(lll)4f电子 ...

对，但小核赝势的情况可以用Hund规则判断，不用试。除非不确定Eu是不是3价，或者除了Eu以外其他地方也有单电子（比如过渡金属、其他稀土、有机自由基），那种情况下可能需要试

作者
Author: halouhapily 时间: 2022-9-14 20:13

wzkchem5 发表于 2022-9-14 17:46
对，但小核赝势的情况可以用Hund规则判断，不用试。除非不确定Eu是不是3价，或者除了Eu以外其他地方也有 ...

好的老师，我明白了，那如果稀土配合物整体显0价呢（Eu离子+3价，配合物-1价*3）？是不是不用单独管Eu离子的最外层电子，整体电荷为0，自旋多重度为1？

作者
Author: wzkchem5 时间: 2022-9-14 20:40

halouhapily 发表于 2022-9-14 13:13
好的老师，我明白了，那如果稀土配合物整体显0价呢（Eu离子+3价，配合物-1价*3）？是不是不用单独管Eu离 ...

自旋多重度和电荷没有任何关系；是否需要考虑外层电子也和电荷没有任何关系，而只和是否用了大核赝势有关。自旋多重度等于自旋向上的单电子数-自旋向下的单电子数+1 （起码当你不开SOC、不用多参考态方法、不做TDDFT时是如此），应该先根据Hund规则来判断自旋向上的单电子数以及自旋向下的单电子数，然后套公式算。不存在“如果配合物整体显x价，那么就不用考虑离子的外层电子/自旋多重度必然是y”的说法

作者
Author: youwillnever 时间: 2022-10-31 20:59

gqy 发表于 2022-5-15 15:18
最近终于完成了结构，和单点的计算。感谢以上各位老师的建议。

请问是如何解决问题的呢？我优化eu3配合物129轮不收敛

自旋多重度应该怎么调整呢？

作者
Author: wzkchem5 时间: 2022-11-1 00:04

youwillnever 发表于 2022-10-31 13:59
请问是如何解决问题的呢？我优化eu3配合物129轮不收敛
自旋多重度应该怎么调整呢？

用Hund规则判断自旋多重度，如果不能唯一判断的话，用Hund规则确定所有可能的自旋多重度（比如你的体系有两个4f7离子，那就只需要试1和15），分别计算，取能量较低的。其中涉及反铁磁耦合的体系（例如我举的例子里面自旋多重度是1的那个计算）需要按http://sobereva.com/82的方法准备片段初猜，光把自旋多重度设对了是不够的

作者
Author: youwillnever 时间: 2022-11-1 12:24

wzkchem5 发表于 2022-11-1 00:04
用Hund规则判断自旋多重度，如果不能唯一判断的话，用Hund规则确定所有可能的自旋多重度（比如你的体系有 ...

好的，我去看看，感谢您的回复

作者
Author: bms 时间: 2024-3-11 19:23

gqy 发表于 2022-5-15 15:18
最近终于完成了结构，和单点的计算。感谢以上各位老师的建议。

你好请问您测试的Eu的自旋多重度是多少呢

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