计算化学公社
标题: 求助用Gaussian的ONIOM模型计算分子的外重组能 [打印本页]
作者Author: 好多于 时间: 2022-12-13 16:25
标题: 求助用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没有动）。请问各位老师，这样近似计算是否可以？ 除此之外，各位老师是否有更好的方法计算分子的外重组能？
作者Author: sobereva 时间: 2022-12-14 14:44
算得动的话，可以用CP2K结合CDFT在DFT下计算外重组能
作者Author: 好多于 时间: 2022-12-14 15:22
好的，谢谢卢老师
| 欢迎光临 计算化学公社 (http://bbs.keinsci.com/) |
Powered by Discuz! X3.3 |