Gaussian混合基组输入文件的写法请教

Dave666

昨天问了个fchk转化出错的问题，后来改成纯的大基组计算，发现没问题。另外纯的大中小基组算出来的能量相差不大，且使用orca用混合基组算了下，能量也相近。这几个值之间相差最多就3 hartree。
但是用gaussian混合基组输入计算的话，能量和前面的数值就相差了2000多hartree，算出来轨道也很奇怪。
所以推测应该是gaussian混合基组的输入文件写的有问题，但也不知道哪里出错了。麻烦各位老师指教。其中oldchk读取的是原来在def2tzvp基组下结构优化得到的chk文件。
输入文件如下：
%rwf=235-wB97XD-def2qzvp_def2tzvp-energy.rwf
%oldchk=235-wB97XD-def2tzvp-opt.chk
%mem=150GB
%nprocshared=64
%chk=235-wB97XD-def2qzvp_def2tzvp-energy.chk
#p guess=read maxdisk=12TB geom=allcheck wb97xd gen

Si Sn 0      
def2qzvp
****
C H N O P 0
def2tzvp
****

zjxitcc

（1）def2系列基组对Sn是有赝势的，因此gjf文件应当写genecp，不能写gen。而且坐标底下还需要指定赝势，例如
C H N O P 0
def2tzvp
****
Si Sn 0      
def2qzvp
****
[空行]
Sn 0
def2
[空行]
[空行]

（2）注意不要犯这三种错误，否则计算大概率不合理（一般来说可以正常运行，但计算过程可能会比较离谱）
i. 当前为带赝势的计算，读取的chk文件却是用全电子基组算的；
ii. 当前为全电子基组计算，读取的chk文件却是在带赝势的情况下计算的;
iii. 当前为带def2赝势的计算，读取的chk文件却是在SDD或LANL2下算的。

（3）“相差最多就3 hartree”这句话是比较可怕的，虽然不一定是错误的。因为Hartree是非常大的单位，1 Hartree = 627.51 kcal/mol。如果在Gaussian/ORCA使用严格一致的分子坐标、泛函、基组计算单点能，对于KS-DFT方法俩程序能量区别通常在1e-3 Hartree数量级，如果ORCA各种阈值设置更tight一点，能量区别还可以降到1e-4 Hartree数量级。如果是HF/MP2等波函数方法，在没有RI加速的情况下，俩程序能量区别通常在1e-6 Hartree数量级甚至更低。有一种特殊情况是两程序SCF收敛到了不同的解，能量会差1e-2 Hartree数量级甚至1e-1。所以你说“相差最多就3 hartree”非常容易让读者误解。