Gaussian-ORCA-xTB联用进行ONIOM计算

winnerwill

最近论坛上好多联用，挺有意思的。以下欢迎拍砖。
这里试着把三个常用的三个量化软件进行ONIOM计算。
用的是社长原创的接口，这里稍微修改了点，以适应不同情形。
请看 gox.zip (1.91 MB, 下载次数: 64)

2020-8-18 06:30 上传

点击文件名下载附件




试了几个例子，貌似还是可行的。
ONIOM用的不多，不恰当的地方请指出。欢迎试算恰当的实际例子。

用之前先把orca.sh中的orcapath定义为自己的实际路径。
其他照常运行，设置好三个程序正常的运行路径和环境变量。


例子一（Gau-xTB）：
%nproc=1
#P ONIOM(B3LYP/def2SVP em=GD3BJ : External="sh ./xtb.sh") Opt(nomicro)



例子二（Gau-ORCA）：
%nproc=1
#P ONIOM(External="sh ./orca.sh" : uff) Opt(nomicro)


例子三（Gau-ORCA-xTB三者联用）：
%nproc=1
#P ONIOM(External="sh ./orca.sh" : External="sh ./xtb.sh") Opt(nomicro)


例子四（Gaussian非联用）：
#P ONIOM(B3LYP/6-31G(d,p):UFF) Opt

以上仅双层的例子，三层的例子应该也可行。

References
http://bbs.keinsci.com/thread-10141-1-1.html
http://bbs.keinsci.com/thread-10106-1-1.html........

喵星大佬

用ONIOM的情况下调用其他软件计算是用内坐标优化效率比较高，直接调用用的笛卡尔坐标，优化效率很低，貌似以前有人说过这事

winnerwill

喵星大佬 发表于 2020-8-18 23:44
用ONIOM的情况下调用其他软件计算是用内坐标优化效率比较高，直接调用用的笛卡尔坐标，优化效率很低，貌似 ...

Opt里面不指定坐标类型的话，貌似ONIOM默认转为冗余内坐标计算，应该是推荐的。
不过如果结构优化过程中存在变成三点一线的原子，则可能会失败报错。
我试了指定z-matrix，反而优化效率更低，甚至在默认的优化圈数内无法收敛。