将小基组的波函数作为大基组的初猜出错：Error reading in TMP file 53/0 (11024)

退休老干部

我想要比较一下大基组和小基组之间的差距，然而小基组scf能收敛，大基组却不能收敛，于是想要尝试用小基组的波函数作为大基组的初猜，但是却出现了如下的错误：
Guess MOs from SCF MO coefficient file
Reading MOs from coefficient file
Size of previous SCF MO coefficient file: 3806036 (expected); 1984030 (actual).

Q-Chem warning in module /scratch/svnadmin/18824_zgan/qchem/gesman/GuessMan.C, line 722:

Inconsistent size for SCF MO coefficient file.  Please check basis and PURECART setting.


Q-Chem fatal error occurred in module /scratch/svnadmin/18824_zgan/qchem/libmdc/newfileman.C, line 293:

Error reading in TMP file 53/0 (11024)

--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
我的输入文件 sp_0.in：
$molecule
read sp_0.xyz
$end

$rem
exchange  xygjos
basis  6-311+G(2d,2p)
aux_basis rimp2-cc-pvtz
solvent_method  pcm
dft_d   d3_bj
mem_static 600
$end

$solvent
dielectric 78.4
dielectric_infi 1.76
$end

@@@

$molecule
read
$end

$rem
exchange  xygjos
basis  6-311+G(3df,2p)
aux_basis rimp2-cc-pvtz
solvent_method  pcm
scf_guess read
dft_d   d3_bj
mem_static 600
$end

$solvent
dielectric 78.4
dielectric_infi 1.76
$end

------------------------------
sp_0.xyz:
$molecule
-4 1
O                  1.92938953   -1.90925047   -2.78254339
P                  2.64900286   -2.43034408   -1.54747799
O                  2.43886233   -3.91093716   -1.27193686
O                  4.09190754   -1.96561986   -1.42313354
O                  1.82813112   -1.60756652   -0.30189251
P                  1.49853568   -1.97599615    1.21026257
O                  0.23725106   -2.78021549    1.30502853
O                  2.70008081   -2.44844239    1.96729720
O                  1.06974219   -0.48936506    1.76124974
P                  1.80037300    0.92183100    2.03478500
O                  1.44362592    1.35933985    3.41759031
O                  3.22117063    0.91362107    1.58174668
O                  0.92123100    1.89099000    1.06108200
C                  1.54609562    2.56910845   -0.02492435
H                  1.96317558    3.51887273    0.32049516
H                  2.35361575    1.96597574   -0.43911304
C                  0.51293488    2.83635412   -1.09786860
H                  0.96643542    3.45372488   -1.87442493
O                  0.09838476    1.59087133   -1.70755680
C                 -1.26157437    1.32103809   -1.43149513
H                 -1.87154460    1.48514001   -2.32051546
N                 -1.43431492   -0.07220631   -1.08773660
C                 -1.00434651   -0.73549183    0.05179249
H                 -0.48937854   -0.22524703    0.84829559
N                 -1.27993648   -2.01162475    0.02001783
C                 -1.91503054   -2.21637529   -1.19328690
C                 -2.44380531   -3.40005631   -1.78159093
O                 -2.46783980   -4.54806883   -1.34791543
N                 -3.00417121   -3.12047159   -3.04906453
H                 -3.41007854   -3.91959579   -3.51600796
C                 -3.05199930   -1.89443762   -3.65613554
N                 -3.68073004   -1.81902861   -4.85356284
H                 -3.82554740   -2.64930026   -5.40314818
H                 -3.55859717   -0.96219446   -5.36740141
N                 -2.55716185   -0.80803143   -3.10753536
C                 -2.01285757   -1.02043746   -1.89449134
C                 -0.76752159    3.49773352   -0.59499122
O                 -1.31104735    4.30990457   -1.62744169
H                 -2.27120073    4.29913311   -1.51478743
H                 -0.59474281    4.08768072    0.30704911
C                 -1.68554267    2.29695173   -0.33012485
H                 -1.46974354    1.88418432    0.65110956
O                 -3.06299070    2.59945306   -0.47384162
H                 -3.42519681    2.83572189    0.38546117
$end

---------------------------------------------------------------------------
命令： qchem -save -nt 16 sp_0.in sp_0.out save
------------------------------------------------------------------------
PS：我换了一个简单的体系H2+测试了一下，可以顺利进行，不知道问题出在哪里，是我的结构有问题吗？可是看报错似乎是tmp文件的问题

hebrewsnabla

你可以去qchem论坛提问

zjxitcc

基组过于弥散，建议只给O原子用弥散基函数，或许就能收敛了，就不需要这么读取了。

hebrewsnabla

scf_guess read 只支持相同基组的波函数读取；如果要读不同基组的，需要设定basis set projection （搜一下手册就知道了）。

H2+可以读是因为对于H这两个基组是一样的。

另外，XYGJ-OS根本没有D3(BJ)参数。

退休老干部

hebrewsnabla 发表于 2022-8-26 15:09
scf_guess read 只支持相同基组的波函数读取；如果要读不同基组的，需要设定basis set projection （搜一下 ...

感谢回复，这个问题我之前已经在论坛上问过了，他也建议我用基组投影，但是后面没有用基组投影啥的，而是增加了双电子积分的精度，就成功收敛了。至于加D3是我无知，干啥都喜欢往上加，被老师提醒过了