GVB计算能量远高于UHF

BinWang

各位老师：
      大家好！我用GAMMES进行GVB计算发现计算出来的能量远远高于UHF的结果，这是什么原因呢，是因为我固定了内层电子轨道的原因吗？我的初始波函数是从UHF读入的（利用MOKIT转化成原始的GVB后通过GAMMES读入）。输入文件如下（MOKIT自动生成）：
    $CONTRL SCFTYP=GVB RUNTYP=ENERGY ICHARG=1 MULT=3 NOSYM=1 ICUT=11
  MAXIT=500 RELWFN=DK ISPHER=1 $END X2C
$SYSTEM MWORDS=1953 $END
$SCF NCO=97 NPAIR=127 NSETO=2 NO(1)=1,1 DIRSCF=.T. CONV=5d-4 $END
$MOFRZ FRZ=.T.
  IFRZ(1)=1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17    #这是固定内层轨道的部分
  IFRZ(18)=18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34
  IFRZ(35)=35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51
  IFRZ(52)=52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68
  IFRZ(69)=69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85
  IFRZ(86)=86,87,88,89,90,91,92,93,94,95,96,97
$END
$GUESS GUESS=MOREAD NORB=981 $END
$DATA
GAMESS inp file produced by MOKIT,nif=981,nbf=981
C1   1
Ir   77.         5.36500002         8.21116298        10.13060499
   S 10
。。。。。。。。。。。。
后面是坐标+基组+波函数，太长省略


UHF和GVB优化出来的能量差距很大，-20659.799516177 Hartree（这个是GAMMES中GVB的结果，全电子基组x2c-svpall） vs -20674.87185862 Hartree(这个是PySCF中UHF的结果，全电子基组x2c-svpall)

hebrewsnabla

mokit 的自动工作流里面，如果开了 x2c，gamess的GVB计算用的是 DKH2, pyscf 用的是 sfx2c ，二者的绝对能量无法比较。如果不开相对论哈密顿，是可以比较的。

BinWang

hebrewsnabla 发表于 2026-1-14 14:53
mokit 的自动工作流里面，如果开了 x2c，gamess的GVB计算用的是 DKH2, pyscf 用的是 sfx2c ，二者的绝对能 ...

非常感谢老师的回复，这种情况下gamess产生的GVB波函数还能顺利读到pyscf里面吗？我有必要将gamess的GVB计算改为sf2c吗？

hebrewsnabla

BinWang 发表于 2026-1-14 15:15
非常感谢老师的回复，这种情况下gamess产生的GVB波函数还能顺利读到pyscf里面吗？能 ...

可以

BinWang

hebrewsnabla 发表于 2026-1-14 15:18
可以

非常感谢老师