计算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) 文件中附带了计算的输入文件关键词和输出文件的报错部分。图片显示了投入计算的分子结构，还有计算的其中一种方法的输入和输出文件。问题似乎有点穷途末路的感觉。求大神指点迷津，万分感谢。

bms · 发表于 Post on 2024-3-11 19:23:23

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

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

youwillnever · 发表于 Post on 2022-11-1 12:24:54

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

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

wzkchem5 · 发表于 Post on 2022-11-1 00:04:04

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

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

youwillnever · 发表于 Post on 2022-10-31 20:59:42

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

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

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

wzkchem5 · 发表于 Post on 2022-9-14 20:40:05

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

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

halouhapily · 发表于 Post on 2022-9-14 20:13:51

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

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

wzkchem5 · 发表于 Post on 2022-9-14 17:46:18

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

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

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。如果用小核赝势的话，基态的自旋多重度需要测试，这样理解对吗

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

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

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

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

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

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

symm=loose，scf=conver=6

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

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

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

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-4-29 20:47:58

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

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

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

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

评分 Rate