计算化学公社

标题: molpro对MoF4O-C2H4做sapt能量分解报错 [打印本页]

作者
Author: 小翟w 时间: 2023-11-20 20:12
标题: molpro对MoF4O-C2H4做sapt能量分解报错
老师，你好，我是做分子间弱相互作用，我的结构是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,

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

作者
Author: 小翟w 时间: 2023-11-22 10:34

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

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

作者
Author: sobereva 时间: 2023-11-22 19:16

小翟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）

作者
Author: 小翟w 时间: 2023-12-2 20:35

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

好的，谢谢老师

作者
Author: zjxitcc 时间: 2023-12-2 21:04
aug-cc-pVDZ基组对Mo元素没有定义，而aug-cc-pVDZ-PP对Mo有定义

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