molpro对MoF4O-C2H4做sapt能量分解报错

小翟w

老师，你好，我是做分子间弱相互作用，我的结构是MoF4O-C2H4，我正在进行的是利用molpro程序进行能量分解（下面是我的输入和输出文件），想问一下为什么使用的基组一直报错
输入文件：
gthresh,energy=1.d-8,orbital=1.d-8,grid=1.d-8
symmetry,nosym
orient,noorient
geometry={
1,Mo1,,0.51187200   -0.00008700   -0.00006800
2,F1,,0.09785600   -1.27297000    1.26759600
3,F2,,0.09975100    1.26533900    1.27572900
4,F3,,0.09766700   -1.26508900   -1.27553100
5,F4,,0.09836300    1.27346400   -1.26723000
6,O1,,2.14700700   -0.00096900   -0.00063800
7,C1,,-2.63728900   -0.66299300    0.00056500
8,C2,,-2.63718900    0.66358700   -0.00019800
9,H1,,-2.64322600   -1.22958200   -0.92096700
10,H2,,-2.64220200   -1.22858100    0.92270400
11,H3,,-2.64229900    1.23019900    0.92131200
12,H4,,-2.64282200    1.22912400   -0.92236200
}

basis={
avdz}
oldnorm=1

ca=2110.2
cb=2120.2
nkerna=1
nkernb=1

{df-hf,basis=avdz; save,2100.2}
edm_hf=energy

dummy,1,2,3,4,5,6
{df-hf,basis=avdz; save,$ca; start,atdens}
ema_hf=energy;
{sapt;monomerA}

dummy,7,8,9,10,11,12
{df-hf,basis=avdz; save,$cb; start,atdens}
emb_hf=energy;
{sapt;monomerB}

{sapt,SAPT_LEVEL=2;intermol,ca=$ca,cb=$cb,icpks=1,fitlevel=3
dfit,basis_coul=avdz,basis_exch=avdz,cfit_scf=3}

eint_hf=(edm_hf-ema_hf-emb_hf)*1000 mH
delta_hf=eint_hf-e1pol-e1ex-e2ind-e2exind

{sapt,SAPT_LEVEL=3;intermol,ca=$ca,cb=$cb,icpks=0,fitlevel=3,nlexfac=0.0
dfit,basis_coul=avdz,basis_exch=avdz,cfit_scf=3}
eint_sapt=e12tot+delta_hf

输出文件：
Variable memory set to    8000000 words,  buffer space   230000 words

ZSYMEL=NOSYM

SETTING BASIS          =    USERDEF
SETTING OLDNORM        =         1.00000000                                 
SETTING CA             =      2110.20000000                                 
SETTING CB             =      2120.20000000                                 
SETTING NKERNA         =         1.00000000                                 
SETTING NKERNB         =         1.00000000                                 


Recomputing integrals since basis changed


Using spherical harmonics

trouble with default basis file; current line:
AVDZ,
complete basis file:
SET,'ORBITAL',
AVDZ,
END,

sobereva

如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并清楚、准确反映出帖子具体内容，避免有任何歧义和含糊性，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题 “molpro sapt能量分解” 改了，以后务必注意

小翟w

sobereva 发表于 2023-11-21 02:15
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此 ...

好的，老师，以后会注意的，那我这个问题可以怎么修改一下

sobereva

小翟w 发表于 2023-11-22 10:34
好的，老师，以后会注意的，那我这个问题可以怎么修改一下

我不用molpro干这个
SAPT用的最多的程序是PSI4
而换做是我我肯定会用sobEDAw研究这种问题
使用sobEDA和sobEDAw方法做非常准确、快速、方便、普适的能量分解分析
http://sobereva.com/685（http://bbs.keinsci.com/thread-39446-1-1.html）

小翟w

sobereva 发表于 2023-11-22 19:16
我不用molpro干这个
SAPT用的最多的程序是PSI4
而换做是我我肯定会用sobEDAw研究这种问题

好的，谢谢老师

zjxitcc

aug-cc-pVDZ基组对Mo元素没有定义，而aug-cc-pVDZ-PP对Mo有定义