求助用Gaussian的ONIOM模型计算分子的外重组能

好多于

请教各位老师，用oniom模型计算分子的外重组能（详见sob老师帖子http://sobereva.com/333），其中进行到第三步时，尝试用-1和opt=readopt方法在ONIOM模型中固定QM部分的坐标，发现无法固定全部的QM部分。（出现All variables have been frozen的报错。）咨询Gaussian客服后，给出micro iteration方法，即：不固定high layer，但是只进行一步优化（maxcyc=1，因此model system不动），并且使用IOp(1/128=1)要求在第一步就完成Micro iteration（low layer分子可以动）。

输入文件关键词如下：

#p opt(maxcyc=1) nosymm iop(1/128=1) oniom...

但是这种情况下结构会显示优化不收敛（因为model system没有动）。请问各位老师，这样近似计算是否可以？ 除此之外，各位老师是否有更好的方法计算分子的外重组能？

sobereva

算得动的话，可以用CP2K结合CDFT在DFT下计算外重组能

好多于

sobereva 发表于 2022-12-14 14:44
算得动的话，可以用CP2K结合CDFT在DFT下计算外重组能

好的，谢谢卢老师