CASSCF求助

zidu113

各位老师：
       大家好。又要麻烦各位老师了。
       我近来看文献，看见以下两段类似的内容：

      1. In order to test further the validity of the one-electron approximation, we made single-point CASSCF (4,22)/6-311 (4 active electrons and 22 active molecular orbitals) calculations using CCSD(T)/6-311 geometries.We found that the Hartree–Fock configuration is dominant in all three anions:
CHF= 0:953 (B3- ), CHF=0.956 (AL3- ), and CHF=0.953 (Ga3-). Thus, the one-electron picture of the bonding is expected to be valid.
   2.The HF wave function was dominant in the CASSCF(7,12)/6-311+G* expansions for B3 (CHF= 0.8926 among the 108 900 configurations) and in the CASSCF(8,12)/6-311+G* for B3- (CHF=0.9153 among the 122 760 configurations), ensuring that the RCCSD(T) calculations were reliable.
      对CASSCF很陌生。实在不知道这两段文字中的CHF代表什么意思 ？为什么作者这么做能确保 RCCSD(T)可靠？
      望能得到各位老师不吝赐教，万分感激。

jiangning198511

CASSCF的波函数是多个SLATER行列式的组合，每个行列式前有个系数，在CASSCF方法中行列式系数和分子轨道系数同时优化。

如果说有个系数非常大，接近0.95的话，说明体系只需要一个行列式就可以描述的比较好，这时可以说明体系的多参考效应比较小，用RCCSD(T)是不错的选择

jiangning198511

但这只是一个标准，还需要看T1诊断等。这个在计算CCSD(T)中可以看到，这个值应该小于0.02(这个值具体记不清楚了，你可以查查相关文献)

zidu113

jiangning198511 发表于 2015-8-27 08:58
但这只是一个标准，还需要看T1诊断等。这个在计算CCSD(T)中可以看到，这个值应该小于0.02(这个值具体记不清 ...

奥，太好了，非常感谢您的回复。我基本明白了。但是还是想请问，这个CASSCF行列式系数CHF在输出文件的什么位置？我怎么找不到呀，。希望再给我这个新手指点指点，谢谢。

sobereva

zidu113 发表于 2015-8-27 11:53
奥，太好了，非常感谢您的回复。我基本明白了。但是还是想请问，这个CASSCF行列式系数CHF在输出文件的什 ...

CHF是指CASSCF多行列式波函数当中的HF行列式的系数(coefficient of HF)。
量化程序计算CASSCF时会输出构成它的大量行列式当中系数相对较大的那些的系数，从中去找对应于HF的那个就行了，也就是按照能量最低方式占据的那个行列式。

zidu113

sobereva 发表于 2015-8-27 16:35
CHF是指CASSCF多行列式波函数当中的HF行列式的系数(coefficient of HF)。
量化程序计算CASSCF时会输出 ...

谢谢，sobereva 老师的耐心及时回复，但还是不太明白，不知这种行列式在输出文件里是什么形式的,请问所谓的行列式是以下这种形式吗？coefficient of HF是红字部分吗？谢谢sobereva老师。
Full Convergence on CI vector
ITN= 25 MaxIt= 64 E=    -74.0145479719 DE=-1.56D-06 Ac谢谢c= 1.00D-08 Lan=  9
In Davids. NSec=    122760
IMaxHF,MaxMI,MSekIn,IFactt,WhenSc,nroot    5  200    2    2    2    1
          (  1)     Eigenvalue  -74.0145479719   
(       1) 0.9153434 (   13861)-0.0866325 (    1128)-0.0866325 (      91)-0.0713950 (      78)-0.0713950 (    1084) 0.0711702 (   13698)-0.0711702
(    1037)-0.0711702 (   13863)-0.0711701 (    1081)-0.0682345 (   14028)-0.0682344 (    1047) 0.0645888 (   13708)-0.0645888 (   13873) 0.0645888
(    1094) 0.0645888 (       6)-0.0609289 (       3)-0.0609288 (      81)-0.0587268 (      68) 0.0587267 (   14196)-0.0463225 (    1225)-0.0463225
(   40756)-0.0448955 (   14077) 0.0401595 (    3404)-0.0365900 (   20707) 0.0365900 (   21737)-0.0365900 (    4006)-0.0365900 (     326)-0.0364451
(     172) 0.0364451 (   13907)-0.0354572 (    1321) 0.0318647 (   14873) 0.0318647 (      73) 0.0314887 (      86)-0.0314887 (    1185)-0.0309748
(   14037) 0.0309748 (     142) 0.0303383 (   13747)-0.0303176 (   14582)-0.0303175 (     153)-0.0289818 (      21)-0.0283938 (    1138) 0.0279988
(   14039) 0.0279988 (   14206)-0.0279988 (    1187) 0.0279988 (     862) 0.0276204 (   15577)-0.0272056 (    1954)-0.0272055 (    1098)-0.0272026
(   13712)-0.0272026 (
Davidson Disk Diagonalization is being used:
MVecIn=      1 MSekIn=    1 NVConv=    1 MConv=    1 WhenSc=  0 IFact=    1 FinIt=F
NRoot=     1 Conv= 0.00D+00 MaxIt=      122760 Dim=      122760 Max sub=   200
Iteration     1 internal iteration     1 dimension is     1

sobereva

zidu113 发表于 2015-8-27 16:54
谢谢，sobereva 老师的耐心及时回复，但还是不太明白，不知这种行列式在输出文件里是什么形式的,请问所谓 ...

是

zidu113

谢谢sobereva 老师。