从发帖到现在几年过去了，这期间时不时有朋友从帖子这边私信过来问如何解决收敛的问题。因为我不是时常在线的状态，所以总是不能及时回复各位道友的信息（拱手）。先把之前解决收敛我的计算流程和一些参数设置放在这条回复里面，后面有问题的朋友也欢迎联系。在这块儿有经验的朋友也欢迎分享一些经验。

当时计算内嵌Eu元素的时候，确实很收敛，开始改变了很多关键词(请仔细参考：http://sobereva.com/61)，进行尝试。而镧系锕系配合物的计算甚至更困难，由于f轨道近简并导致SCF比过渡金属配合物还更难收敛(请仔细参考：http://sobereva.com/581)。最后是先用小的基组初步优化，（以及用SCF=conver=N关键词改变收敛限，我先设置到N为6或者5，如：#P PBE1PBE/genecp pseudo=read opt freq SCF=(Fermi,vshift=500,conver=6,maxcyc=180) IOp(5/13=1) 进行粗略优化）提取出优化后的坐标，再用大一点的基组，（这时收敛限为正常收敛）尝试不同的关键词看哪一个收敛更好。
#P PBE1PBE/genecp pseudo=read opt freq SCF=(vshift=450,novaracc,maxcyc=180)  

坐标

C 0                                                               
6-31G*                                                              
****                                                               
Eu 0                                                                 
SDD                                                            
****                                                               
                                                                  
Eu 0                                                                 
SDD

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

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

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

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

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

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

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

请问是如何解决问题的呢？我优化eu3配合物129轮不收敛  
自旋多重度应该怎么调整呢？

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

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

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

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

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

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

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

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

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

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

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

symm=loose，scf=conver=6

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

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

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