计算化学公社

标题: CASSCF求助 [打印本页]

作者
Author: zidu113 时间: 2015-8-27 00:56
标题: CASSCF求助
本帖最后由 zidu113 于 2015-8-27 01:03 编辑

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

   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)可靠？
   望能得到各位老师不吝赐教，万分感激。

作者
Author: jiangning198511 时间: 2015-8-27 08:56
CASSCF的波函数是多个SLATER行列式的组合，每个行列式前有个系数，在CASSCF方法中行列式系数和分子轨道系数同时优化。

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

作者
Author: jiangning198511 时间: 2015-8-27 08:58
但这只是一个标准，还需要看T1诊断等。这个在计算CCSD(T)中可以看到，这个值应该小于0.02(这个值具体记不清楚了，你可以查查相关文献)

作者
Author: zidu113 时间: 2015-8-27 11:53

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

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

。希望再给我这个新手指点指点，谢谢。

作者
Author: sobereva 时间: 2015-8-27 16:35

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

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

作者
Author: zidu113 时间: 2015-8-27 16:54

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

作者
Author: sobereva 时间: 2015-8-27 16:54

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

是

作者
Author: zidu113 时间: 2015-8-27 18:15
谢谢sobereva 老师。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)