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”非常容易让读者误解。

TYXB521

zjxitcc 发表于 2026-2-27 19:16
（1）def2系列基组对Sn是有赝势的，因此gjf文件应当写genecp，不能写gen。而且坐标底下还需要指定赝势，例 ...

%chk=C:\Users\1\Desktop\PHBAK.chk
%NProcShared=8
%mem=12GB
#p opt=(cartesian,calcfc) freq b3lyp/gen  pseudo=read scrf=(smd,solvent=generic,read)  em=gd3bj

Title Card Required

0 1
C                  0.27153041   -1.62863440    0.00263100
C                  1.66669041   -1.62863440    0.00263100
C                  2.36422841   -0.42088340    0.00263100
C                  1.66657441    0.78762560    0.00143200
C                  0.27174941    0.78754760    0.00095300
C                 -0.42585159   -0.42065840    0.00194900
H                 -0.27822859   -2.58095140    0.00308100
H                  3.46390841   -0.42080340    0.00326500
H                  2.21677441    1.73976860    0.00137300
H                 -1.52545559   -0.42047540    0.00176900
O                 -0.44356461    2.02578195   -0.00028617
C                  2.43624276   -2.96257061    0.00447258
O                  2.75382289   -3.50923077    1.09253328
K                 -0.88924327    2.80278472   -2.53676744
O                  2.79115873   -3.58214747   -1.23450247
K                  4.11451558   -5.87268236   -0.74636670

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

K 0
LANL2DZ

Eps=10.0
EpsInf=2.42
HbondAcidity=0.60
HbondBasicity=0.30
SurfaceTensionAtInterface=40.9
CarbonAromaticity=0.857
ElectronegativeHalogenicity=0.0
老师您好，我这个也是混合基组进行计算的，这几天一直在试错，没在官网找到对应的错误原因，您能帮我看一下是什么问题吗，溶剂模型选用的苯酚

zjxitcc

TYXB521 发表于 2026-6-27 13:04
%chk=C:%users\1\Desktop\PHBAK.chk
%NProcShared=8
%mem=12GB

上传gjf文件和输出文件。不要短时间内重复发帖
http://bbs.keinsci.com/thread-59757-1-1.html