计算化学公社

标题: 请各位大佬帮忙给点意见，尝试了卢老师博文中的所有方法，均SCF不收敛 [打印本页]

作者
Author: nianbin 时间: 2021-7-6 15:49
标题: 请各位大佬帮忙给点意见，尝试了卢老师博文中的所有方法，均SCF不收敛
本帖最后由 nianbin 于 2021-7-6 15:51 编辑

我想做一下两个Fe3+离子与有机小分子的配合物的结构优化进而进一步算出binding energy，但是这个结构始终无法SCF收敛，用了卢老师博文http://sobereva.com/61?tdsourcetag=s_pctim_aiomsg，中的解决SCF收敛的办法，仍然无法奏效，请各位老师给点意见，谢谢

作者
Author: wzkchem5 时间: 2021-7-6 16:04
结构不合理。建议找本配位化学的教材看一下，Fe3+是怎么配位的。
不是你随便把两个Fe3+和有机分子放在一起，然后优化就可以的。初始结构必须大体上合理，不能指望如此异想天开的结构，高斯还能给你优化出来。竟然还让两个Fe3+彼此成键，真不知道怎么想的。就算你一点配位化学都没学过，哪怕学过初中物理都知道，电荷同性相斥，结果你硬是把两个高达+3价的离子放得这么近，都成键了，你觉得能合理吗？

作者
Author: wzkchem5 时间: 2021-7-6 16:07
一般而言，现实中不存在的结构，SCF出现不收敛问题的概率会大得多。这是因为，一般的DFT方法本来就是为那些合理的结构而设计的，对于不合理的结构，DFT的基本假设（体系可以用单参考态波函数定性描述）根本不成立，所以第一容易不收敛，第二即使收敛了结果也没有用。

作者
Author: fatpig 时间: 2021-7-6 19:35
能设成singlet基本表示...楼主就是算着玩的，所以不收敛也很正常吧...

作者
Author: nianbin 时间: 2021-7-6 19:44

fatpig 发表于 2021-7-6 19:35
能设成singlet基本表示...楼主就是算着玩的，所以不收敛也很正常吧...

1 3 5 7 都试了

作者
Author: nianbin 时间: 2021-7-6 19:45

wzkchem5 发表于 2021-7-6 16:07
一般而言，现实中不存在的结构，SCF出现不收敛问题的概率会大得多。这是因为，一般的DFT方法本来就是为那些 ...

多谢老师，这个结构是使用molclus联合xtb搜索出来的，应该是有问题，我重新弄一下

作者
Author: k64_cc 时间: 2021-7-6 20:18
xtb不少时候都不对，还是要把化学直觉放第一位呀。

作者
Author: wzkchem5 时间: 2021-7-6 20:24

nianbin 发表于 2021-7-6 12:45
多谢老师，这个结构是使用molclus联合xtb搜索出来的，应该是有问题，我重新弄一下

xtb你设置了电荷了吗？虽然xtb不是特别靠谱，但是不至于让两个带+3价的Fe跑到一起。估计你是没有设置体系总电荷，导致xtb按总电荷为0算的

作者
Author: nianbin 时间: 2021-7-6 20:52

wzkchem5 发表于 2021-7-6 20:24
xtb你设置了电荷了吗？虽然xtb不是特别靠谱，但是不至于让两个带+3价的Fe跑到一起。估计你是没有设置体系 ...

老师好，xtb设置了电荷+6

作者
Author: nianbin 时间: 2021-7-6 20:52

k64_cc 发表于 2021-7-6 20:18
xtb不少时候都不对，还是要把化学直觉放第一位呀。

是的，感谢

作者
Author: wzkchem5 时间: 2021-7-6 22:58

nianbin 发表于 2021-7-6 12:44
1 3 5 7 都试了

这个分子是弱场配体，两个Fe都是高自旋d5，所以应该是11重态。即使更低的自旋态更稳定，也应该从11重态做broken symmetry。
所以还是那句话，要靠化学直觉，不能不动脑子盲试。别人试到7重态就不试了，那是因为别人是单金属配合物，或者强场配合物，或者不成对d电子没那么多。不代表你也是试到7重态就可以停了。

作者
Author: wzkchem5 时间: 2021-7-6 23:01

nianbin 发表于 2021-7-6 13:52
老师好，xtb设置了电荷+6

我用GFN2-xTB给你的初始结构做结构优化，电荷+6，自旋多重度11，发现有一个Fe会跑开。你现在这个结构如果是用xtb找的，肯定说明还是算得有问题，比如电荷没传对，或者自旋多重度没设对。
建议仔细阅读xtb的输出文件（不是molclus的输出文件），确保你所有都设对了，也都正确地传给xtb了。

作者
Author: nianbin 时间: 2021-7-7 00:14

wzkchem5 发表于 2021-7-6 22:58
这个分子是弱场配体，两个Fe都是高自旋d5，所以应该是11重态。即使更低的自旋态更稳定，也应该从11重态做 ...

非常感谢wzkchem5老师

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