计算化学公社

标题: Gaussian中的ONIOM（Gaussian:XTB）金属体系报错 [打印本页]

作者
Author: lianghuashou 时间: 2023-10-16 18:07
标题: Gaussian中的ONIOM（Gaussian:XTB）金属体系报错
各位老师好，我用Gaussian 中的ONIOM（Gaussian:XTB）计算金属体系，采用的接口是SOB老师的脚本xtb.sh ，但是金属的赝势基组读取出错，请问各位老师如何解决。
我的输入文件如下

%nosave
%nproc=20
%chk=$GAUSS_SCRDIR/$JOB_NAME.chk
%mem=100GB
#p  ONIOM(PBE1PBE/genecp: External="sh /scr1/xtb_6_6_1/xtb_interface/xtb.sh" )=Embedcharge  Opt(nomicro)

Title Card Required

0 1
C             0 14.47300000 14.48100000 8.83400000 H
C             0 12.97800000 14.63400000 8.97400000 H
.
.
Pt             0 14.31600000 14.56800000 12.76200000 H

Pt 0
LANL2DZ
****
H C O N 0
6-31G(d)
****

Pt 0
LANL2DZ

报错信息如下
(, 下载次数 Times of downloads: 3)

作者
Author: sobereva 时间: 2023-10-18 04:10
下文提到
将Gaussian与Grimme的xtb程序联用搜索过渡态、产生IRC、做振动分析
http://sobereva.com/421（http://bbs.keinsci.com/thread-10106-1-1.html）

如果你要用本文的做法计算很重的元素，比如Pt，有一点主要注意，即虽然Gaussian与xtb联用时Gaussian自己并不会去做任何量化计算，但是由于你没写基组，Gaussian程序内部默认当前设的基组是STO-3G，程序会进而判断你当前体系里的所有元素是否对于STO-3G有定义，如果没有就报错。对很重的元素STO-3G都是没有定义的。为了避免这个报错，解决办法就是关键词里写上UGBS，这是一个几乎涵盖整个周期表的基组，因此就不会因为当前的基组对元素没有定义而报错了

欢迎光临计算化学公社 (http://bbs.keinsci.com/)