ONIOM模型重组能分解求助

Doraemon10001

wal 发表于 2026-6-16 09:33
我个人是H和L一起优化的。你可以根据你的目的自行取舍

你好，我用您说的方法去尝试了，但是体系太大（已经减少到1230多个原子了），发现底层用PM7基态优化都算不动，感觉底层还是要用力场算，力场似乎加不了TD,您有什么办法吗。

wal

Doraemon10001 发表于 2026-6-17 09:44
你好，我用您说的方法去尝试了，但是体系太大（已经减少到1230多个原子了），发现底层用PM7基态优化都算 ...

这样的话我就没招了。等其他人来解答吧

Doraemon10001

wal 发表于 2026-6-17 13:30
这样的话我就没招了。等其他人来解答吧

谢谢你，已经帮助我很多了，我就是很好奇文献中有许多都是用这个方法的，用QM/MM都成功分解了重组能，可能是用Gaussian自带的fcht就是分解不了这种情况，我还找了Dushin分解代刷，但都说体系太大分解不了，应该都是不愿意修改阈值，重新编译，可能得用MOMAP这种收费软件了，还是希望能有其他人能解答。

Kalinite

使用静电嵌入的时候L层只以静电荷的形式出现，且你说你在优化的时候冻结了L层，也就是说L层完全可以不需要以原子的形式出现，直接改成背景电荷就行了。另外AIE的体系里H层和L层之间大概也不涉及切断化学键。
考虑以上这点的话，你做的ONIOM实际上和做带背景电荷的QM没有区别，那就没有必要做ONIOM。那么HR计算的问题也能解决了。

Kalinite

Doraemon10001 发表于 2026-6-17 14:23
谢谢你，已经帮助我很多了，我就是很好奇文献中有许多都是用这个方法的，用QM/MM都成功分解了重组能，可 ...

QM/MM分解重组能又是另一回事。静电嵌入的框架下，QM/MM除了极少数in group code可以做解析Hessian（依然昂贵），绝大多数的文章都是只对QM区进行差分得到QM区近似的Hessian，然后再分解重组能。这和QM分解重组能没有区别。