计算化学公社

标题: 求助M06 int=ultrafine方法可以找到的过渡态，用b3lyp int=ultrafine优化就震荡 [打印本页]

作者
Author: 化小白 时间: 2024-4-14 15:09
标题: 求助M06 int=ultrafine方法可以找到的过渡态，用b3lyp int=ultrafine优化就震荡
老师，您好，我在做Ce(III)催化两个有机物之间的反应，由于实验是Ce(III)催化，那整个反应体系的电荷和自旋多重度就是(0,2)，就很难收敛，就先用(+1,0)找过渡态，然后再用(0,2)优化，但是有个过渡态用M06 int=ultrafine可以找到，走IRC也是想要的，但是用b3lyp int=ultrafine优化就一直震荡，现在用opt(gdiis,maxstep=1或3或5或10,notrust，loose)都试了都震荡，想问一下老师会不会这个过渡态不存在？但是用b3lyp在电荷和自旋多重度为(+1,0)找过渡态也是可以找到，想问一下老师我能不能把M06 int=ultrafine方法找到的结构全部冻结再用b3lyp做优化？

作者
Author: wjc404 时间: 2024-4-14 15:55
本帖最后由 wjc404 于 2024-4-14 16:02 编辑

如果是SCF收敛有困难，且确定是Ce3+催化，可以在几何优化时给Ce用MWB47大核赝势（这样总比一开始假设是Ce4+去优化来得合理）（4F被赝化了，因此自旋多重度设为1），后续算单点再用MWB28小核赝势或全电子基组。

稀土配合物体系B3LYP不是最合适的，PBE0-D3(BJ)以及wB97M-V应该更好。

作者
Author: sobereva 时间: 2024-4-15 04:45
根本不存在(+1,0)，自旋多重度不可能为0
这种体系用B3LYP本来就是下策，还不如M06适合，PBE0-D3(BJ)才是应该优先考虑的
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272（http://bbs.keinsci.com/thread-536-1-1.html）

试maxstep=1这么小的值根本没意义，与其试这个远不如去试诸如maxstep=5,calcfc,recalc=5

作者
Author: 化小白 时间: 2024-4-15 09:47

wjc404 发表于 2024-4-14 15:55
如果是SCF收敛有困难，且确定是Ce3+催化，可以在几何优化时给Ce用MWB47大核赝势（这样总比一开始假设是Ce4+ ...

好的，谢谢老师，我试试

作者
Author: 化小白 时间: 2024-4-15 09:48

sobereva 发表于 2024-4-15 04:45
根本不存在(+1,0)，自旋多重度不可能为0
这种体系用B3LYP本来就是下策，还不如M06适合，PBE0-D3(BJ)才是应 ...

好的，谢谢老师，我马上试试

作者
Author: 化小白 时间: 2024-4-15 10:36

化小白发表于 2024-4-15 09:48
好的，谢谢老师，我马上试试

老师，我这一写recalc=5就报错，说写错了关键词，去掉就能正常运行，这咋办老师

作者
Author: sai77 时间: 2024-4-15 10:43

化小白发表于 2024-4-15 10:36
老师，我这一写recalc=5就报错，说写错了关键词，去掉就能正常运行，这咋办老师

升级到G16

作者
Author: 化小白 时间: 2024-4-15 10:54

化小白发表于 2024-4-15 09:47
好的，谢谢老师，我试试

老师，我对赝势基组的含义了解的太少了，目前只用过SDD和lanl2DZ，不太懂4f被赝化是什么意思，可以理解为是用了MWB47这个基组，对于Ce(III)把电荷和自旋多重度写为（0，1）是合理的吗？可以这样理解吗老师

作者
Author: 化小白 时间: 2024-4-15 10:57

sai77 发表于 2024-4-15 10:43
升级到G16

老师，我们实验室没有G16，那这个方法应该是用不了了，老师还有没有其他办法

作者
Author: sobereva 时间: 2024-4-16 00:56

化小白发表于 2024-4-15 10:54
老师，我对赝势基组的含义了解的太少了，目前只用过SDD和lanl2DZ，不太懂4f被赝化是什么意思，可以理解为 ...

4f电子用赝势代替，不直接被赝势基组描述
赝势的函数形式以及在量子化学程序中定义的方式
http://sobereva.com/188
详解Gaussian中混合基组、自定义基组和赝势基组的输入
http://sobereva.com/60

作者
Author: 化小白 时间: 2024-4-16 15:56

sobereva 发表于 2024-4-16 00:56
4f电子用赝势代替，不直接被赝势基组描述
赝势的函数形式以及在量子化学程序中定义的方式
http://sober ...

老师，那我就不在(0,2)的情况下优化了，直接改为用MWB47基组在(0,1)的时候优化，然后再换小核赝势基组做单点，可是这样优化出来的.log文件依然显示电荷为1呀，感觉跟我在(1,1)时用SDD优化出来的好像一样，这样好像体现不出来Ce(III)是三价，好像还是Ce(IV)+

作者
Author: sobereva 时间: 2024-4-17 04:29

化小白发表于 2024-4-16 15:56
老师，那我就不在(0,2)的情况下优化了，直接改为用MWB47基组在(0,1)的时候优化，然后再换小核赝势基组做 ...

仔细把此文看了
使用Gaussian做镧系金属配合物的量子化学计算
http://sobereva.com/581（http://bbs.keinsci.com/thread-21274-1-1.html）

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