含过渡金属的团簇几何优化scf不收敛

wyqms · 发表于 Post on 2022-6-28 22:08:32

各位老师好，本人在利用Gaussian16进行PCN-777（MOF）的优化时，遇到了SCF不收敛问题，研读了sob老师的解决SCF不收敛问题的方法 - 思想家公社的门口：量子化学·分子模拟·二次元 (sobereva.com)这篇博文，尝试了其中几种方法后，依旧难以收敛，考虑是不是截取的片段有问题，请各位老师帮忙检查一下！

wzkchem5 · 发表于 Post on 2022-6-28 22:59:19

团簇电荷不对。仔细数一下总体系的电子数，确定正确的体系总电荷以后，用正确的电荷做计算

wyqms · 发表于 Post on 2022-6-28 23:20:18

wzkchem5 发表于 2022-6-28 22:59
团簇电荷不对。仔细数一下总体系的电子数，确定正确的体系总电荷以后，用正确的电荷做计算

老师，您的意思是整个团簇电荷不为0？

wzkchem5 · 发表于 Post on 2022-6-29 00:13:26

wyqms 发表于 2022-6-28 16:20
老师，您的意思是整个团簇电荷不为0？

对。算一下每个金属和配体的电荷，加起来就知道了

wyqms · 发表于 Post on 2022-6-29 11:07:42

wzkchem5 发表于 2022-6-29 00:13
对。算一下每个金属和配体的电荷，加起来就知道了

老师您好，我有个地方不明白，在这个MOF中金属中心的氧原子处于对称位置化合价应该是相同的吧？为什么一个是-2价另一个是-1价，是这个单位分子式标错了还是我理解错了？

wzkchem5 · 发表于 Post on 2022-6-29 14:04:58

wyqms 发表于 2022-6-29 04:07
老师您好，我有个地方不明白，在这个MOF中金属中心的氧原子处于对称位置化合价应该是相同的吧？为什么一 ...

是标错了。不要盲目相信结构式里面画的价态，有的时候结构式只是一个示意图，并不是里面每一个字都是准确的。即便结构式里的价态标的是对的，也要通过化学常识去验算

wyqms · 发表于 Post on 2022-6-30 23:53:40

本帖最后由 wyqms 于 2022-7-1 00:01 编辑

wzkchem5 发表于 2022-6-29 14:04
是标错了。不要盲目相信结构式里面画的价态，有的时候结构式只是一个示意图，并不是里面每一个字都是准确 ...

老师，您好，实在不好意思，还是需要您来解惑，我按照理解数了电荷，整体是2个正电荷，Delta-E在循环了45次后一直都维持在小数点后4位，目前第99次循环的Delta-E= -0.000106736803，RMSDP=9.19D-06，始终无法逼近08，麻烦老师帮我看看是我结构设计有问题还是说本身就很难收敛，而且在尝试了sob老师解决SCF不收敛问题的方法 - 思想家公社的门口：量子化学·分子模拟·二次元 (sobereva.com)中的前4条方案发现每一步相比不加关键词的Delta-E，差别在小数点后9位。这已经困扰我十多天了，麻烦老师解惑了！

wzkchem5 · 发表于 Post on 2022-7-1 01:51:13

wyqms 发表于 2022-6-30 16:53
老师，您好，实在不好意思，还是需要您来解惑，我按照理解数了电荷，整体是2个正电荷，Delta-E在循环了45 ...

电荷还是数错了。如果用金属电荷数、配体电荷数加起来的方法数不清楚，那就一个一个数电子，用总核电荷数减去总电子数的方法计算。

wyqms · 发表于 Post on 2022-7-1 12:05:00

本帖最后由 wyqms 于 2022-7-1 12:06 编辑

wzkchem5 发表于 2022-7-1 01:51
电荷还是数错了。如果用金属电荷数、配体电荷数加起来的方法数不清楚，那就一个一个数电子，用总核电荷数 ...

老师您好，电荷我是按照Zr周围36个电子，C-O键中的氧为周围9个电子，周围三个O-Zr键的O为10个电子计算的电荷，我建立的体系总核电荷数是710，算出来还是2个正电荷，请问老师是我化合价理解错了还是说其他的原因？

wzkchem5 · 发表于 Post on 2022-7-1 14:51:56

wyqms 发表于 2022-7-1 05:05
老师您好，电荷我是按照Zr周围36个电子，C-O键中的氧为周围9个电子，周围三个O-Zr键的O为10个电子计算的 ...

这个结构里有的C-O是单键，有的是双键，此外有的Zr-O键是共价键、有的是配位键，数电子的时候考虑到这两个问题了吗

wyqms · 发表于 Post on 2022-7-1 16:13:44

wzkchem5 发表于 2022-7-1 14:51
这个结构里有的C-O是单键，有的是双键，此外有的Zr-O键是共价键、有的是配位键，数电子的时候考虑到这两 ...

老师，周围有3个Zr-O键的氧是0价吗，应该按照 8个电子计算吗？我算的时候是按照结构式标的-2价算的，是不是这个地方不能用-2价算？

wzkchem5 · 发表于 Post on 2022-7-1 16:16:43

wyqms 发表于 2022-7-1 09:13
老师，周围有3个Zr-O键的氧是0价吗，应该按照 8个电子计算吗？我算的时候是按照结构式标的-2价算的，是不 ...

这里面所有氧都是-2价的

wyqms · 发表于 Post on 2022-7-1 17:25:48

wzkchem5 发表于 2022-7-1 16:16
这里面所有氧都是-2价的

老师，我重新算的是-10，另外，关于计算MK电荷计算的关键词，手册只给了Pop=MK IOp(6/50=1)，但是我在本论坛看到有人用了pop=MK iop(6/33=2,6/41=10,6/42=17)，后面的6/33这个代表的意义试什么，想要计算MK不需要加上6/50=1吗？

wzkchem5 · 发表于 Post on 2022-7-1 17:30:46

wyqms 发表于 2022-7-1 10:25
老师，我重新算的是-10，另外，关于计算MK电荷计算的关键词，手册只给了Pop=MK IOp(6/50=1)，但是我在本 ...

这种问题自己查IOp手册，http://gaussian.com/iops/

wyqms · 发表于 Post on 2022-7-1 20:02:38

wzkchem5 发表于 2022-7-1 17:30
这种问题自己查IOp手册，http://gaussian.com/iops/

多谢老师的指点，我马上去优化，谢谢老师了

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

[新手求助] 含过渡金属的团簇几何优化scf不收敛

浏览过的版块