计算化学公社

标题: ONIOM自定义基组求助 [打印本页]

作者
Author: yundove 时间: 2019-3-9 16:10
标题: ONIOM自定义基组求助
本帖最后由 yundove 于 2019-3-9 16:15 编辑

各位老师好，请问在使用ONIOM计算时，怎么进行自定义基组输入？我试算了一个例子，下面是输入文件。
%NProcShared=24
%Mem=12GB
# oniom(qcisd(t)/cc-pvdz:m062x/gen)

Title Card Required

0 1 0 1 0 1
C             0 -4.17102978 0.46247818 0.00000000 L
H             0 -3.81435694 0.96687637 0.87365150 L
H             0 -3.81435694 0.96687637 -0.87365150 L
H             0 -5.24102978 0.46249136 0.00000000 L
C             0 -3.65771407 -0.98945397 0.00000000 L H 8
H             0 -4.01438509 -1.49385152 0.87365262 L
H             0 -4.01438879 -1.49385337 -0.87365004 L
C             0 -2.11771407 -0.98947218 -0.00000377 H
H             0 -1.58453936 -0.06177356 -0.00000576 H
O             0 -1.48421519 -2.07678526 -0.00000452 H

-H    0
S 3 1.00
   34.0613410             0.60251978E-02
   5.1235746             0.45021094E-01
   1.1646626             0.20189726
S 1 1.00
   0.32723041          1.0000000
S 1 1.00
   0.10307241          1.0000000
P 1 1.00
   0.8000000             1.0000000
****
-C    0
S 6 1.00
  13575.3496820             0.22245814352E-03
2035.2333680             0.17232738252E-02
463.22562359          0.89255715314E-02
131.20019598          0.35727984502E-01
   42.853015891          0.11076259931
   15.584185766          0.24295627626
S 2 1.00
   6.2067138508          0.41440263448
   2.5764896527          0.23744968655
S 1 1.00
   0.57696339419       1.0000000
S 1 1.00
   0.22972831358       1.0000000
S 1 1.00
   0.95164440028E-01          1.0000000
S 1 1.00
   3.1721480009E-02          1.0000000
P 4 1.00
   34.697232244          0.53333657805E-02
   7.9582622826          0.35864109092E-01
   2.3780826883          0.14215873329
   0.81433208183       0.34270471845
P 1 1.00
   0.28887547253       1.0000000
P 1 1.00
   0.10056823671       1.0000000
P 1 1.00
   3.3522745570E-02    1.0000000
D 1 1.00
   1.09700000          1.0000000
D 1 1.00
   0.31800000          1.0000000
F 1 1.00
   0.76100000          1.0000000
****
-O    0
S 6 1.00
  27032.3826310             0.21726302465E-03
4052.3871392             0.16838662199E-02
922.32722710          0.87395616265E-02
261.24070989          0.35239968808E-01
   85.354641351          0.11153519115
   31.035035245          0.25588953961
S 2 1.00
   12.260860728          0.39768730901
   4.9987076005          0.24627849430
S 1 1.00
   1.1703108158          1.0000000
S 1 1.00
   0.46474740994       1.0000000
S 1 1.00
   0.18504536357       1.0000000
S 1 1.00
   6.1681787857E-02    1.0000000
P 4 1.00
   63.274954801          0.60685103418E-02
   14.627049379          0.41912575824E-01
   4.4501223456          0.16153841088
   1.5275799647          0.35706951311
P 1 1.00
   0.52935117943       1.0000000
P 1 1.00
   0.17478421270       1.0000000
P 1 1.00
   5.8261404233E-02    1.0000000
D 1 1.00
   2.31400000          1.0000000
D 1 1.00
   0.64500000          1.0000000
F 1 1.00
   1.42800000          1.0000000
****

M06-2X计算使用的基组是ma-TZVP（直接从https://comp.chem.umn.edu/basissets/ma-TZVP.gbs拷贝过来的）
在对高层进行DFT计算时会报错（已附上输出文件），不知道是什么原因，烦请各位老师帮忙解答，谢谢！

作者
Author: sobereva 时间: 2019-3-9 16:18
你当前的计算设定毫无意义
qcisd(t)/cc-pvdz完全是用来当做典型反面教材的计算级别，结果精度还没你本来要用的低层好
就这么点体系还用ONIOM，完全是自找麻烦，还由于对层之间的耦合考虑不充分导致降低精度，完全不是该用ONIOM的场合
直接用和ma-TZVP大小相仿佛的may-cc-pVTZ就完了，省得你自定义基组

作者
Author: yundove 时间: 2019-3-9 16:24

sobereva 发表于 2019-3-9 16:18
你当前的计算设定毫无意义
qcisd(t)/cc-pvdz完全是用来当做典型反例的计算级别，结果精度还没你本来要用的 ...

感谢sob老师的回复！我刚接触ONIOM，想先在小体系下试算一下，然后再用到大分子体系。现在这样的方法设置是有问题，之后会用ccsd(t)/cbs做高层计算。我只是想知道这样设置能不能计算成功。参照您的http://sobereva.com/60感觉这样设置应该是能跑通的，但是为什么计算会报错呢？

作者
Author: sobereva 时间: 2019-3-9 16:30

yundove 发表于 2019-3-9 16:24
感谢sob老师的回复！我刚接触ONIOM，想先在小体系下试算一下，然后再用到大分子体系。现在这样的方法设置 ...

把输入文件作为附件传上来
而且别用G09 A.01这么老的版本，bug一大堆

作者
Author: yundove 时间: 2019-3-9 16:33

sobereva 发表于 2019-3-9 16:30
把输入文件作为附件传上来

已上传，麻烦sob老师了！

作者
Author: sobereva 时间: 2019-3-10 05:57

yundove 发表于 2019-3-9 16:33
已上传，麻烦sob老师了！

普通计算时某些设定能用、设置方式合理，不代表结合ONIOM的时候就能正常用

当前情况应该算是Gaussian的bug，没法正确对低层级别正确读入自定义基组设定。改用质量和大小相仿佛的6-311+G(2df,p)，或者may-cc-pVTZ亦可（后者给更贵点）

作者
Author: yundove 时间: 2019-3-10 13:50

sobereva 发表于 2019-3-10 05:57
普通计算时某些设定能用、设置方式合理，不代表结合ONIOM的时候就能正常用

当前情况应该算是Gaussian ...

好的，非常感谢sob老师！！

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