重组能计算，中性分子得电子结构优化会断键问题

LC98 · 发表于 Post on 2025-5-22 15:10:16

我最近在进行重组能计算，发现部分有机小分子，在电荷为0，多重度为1的时候优化结构正常，但是当设置电荷为-1 多重度为2的时候会发生断键，这种断键方式有所不同，断键越离谱，算出来的重组能越大，我认为是不合理的。开始我是在真空下计算，后来我觉得是因为没有考虑极化，所以我加了SMD溶剂模型，但是也会出现同样的情况，后面看文献，相同的分子，文献算出来的值远小于我算的值(他用的B3LYP /6-311+G(d）)，同时我换了理论方法，从开始的B3LYP 6-31G*换成了M062X 6-31G*，但是还是会发生断键，我想请问如何解决这种断键问题，或者断键后结构浮动不要那么大。

wzkchem5 · 发表于 Post on 2025-5-22 15:20:54

有的分子在得到电子之后本来就会断键，这是正常的

LC98 · 发表于 Post on 2025-5-22 15:29:39

wzkchem5 发表于 2025-5-22 15:20
有的分子在得到电子之后本来就会断键，这是正常的

老师是这样的，我看文献里面，相同的分子，看他算出来的重组能比较低，符合sob老师帖子里大致的区间范围，但是我算的有的甚至达到了4个电子伏特以上，所以我认为我算的可能有问题，如果说我进行限制性优化您觉得合理吗？跟完全放开的有比较的价值吗？

wzkchem5 · 发表于 Post on 2025-5-22 15:46:54

LC98 发表于 2025-5-22 15:29
老师是这样的，我看文献里面，相同的分子，看他算出来的重组能比较低，符合sob老师帖子里大致的区间范围 ...

如果你的理论级别和文献严格保持一致，能重复出文献的结果吗

LC98 · 发表于 Post on 2025-5-22 15:57:13

我感觉重复不出来，因为他算出来的这个结构，我认为是不会断键的，并且我初始优化是在B3LYP 6-31G*这个级别下，但是还是会断键，根据后续算出来的值，远大于他的值，他是0.33我大约是4.5eV。并且是这样的，我在其他没有断键或者断键不离谱的分子可以和他得到类似的结果，我用的是M062X def2TZVP em=gd3。所以我觉得是结构差异导致的结果差异。

gezq · 发表于 Post on 2025-5-23 00:43:38

本帖最后由 gezq 于 2025-5-23 00:44 编辑

一个有机分子捕获一个热电子，有时会发生电子贴附解离，断键确实是有可能的

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

[新手求助] 重组能计算，中性分子得电子结构优化会断键问题

浏览过的版块