求助双自由基特征值y的计算问题

不干正事

大家好，对于一个双自由基体系，已经用guess=mix关键词得到了残缺的波函数，然后根据之前得到的一个方法，用图中显示的一个方法计算得到的，感觉不是特别正确。然后查相关文献也看到了一个公式来计算y值，但不清楚公式中的各项需要代入什么？

zjxitcc

（1）不是“残缺的波函数”，是对称破缺波函数，或自旋极化波函数；
（2）你截图中并不是完整的步骤，也不是最简步骤，因此没法确定你算对了。这里以键解离的H2举个例子，输入文件如下

1. #p UBHandHLYP/6-31G(d) nosymm guess=mix stable=opt pop=no
   
2. 
   
3. title
   
4. 
   
5. 0 1
   
6. H   0.0   0.0   0.0
   
7. H   0.0   0.0   2.0
   

复制代码

到输出文件里找到Natural Orbital Coefficients，底下Eigenvalues处，找到HONO和LUNO的占据数，代入公式算y值即可。如果你一眼看不出HONO和LUNO在哪里，就理解为HOMO、LUMO序号上的轨道的占据数（提示：对于一个偶数电子的体系，HOMO序号 = 电子数/2）。另外，你这贴出的“文献计算y值的方法”公式与http://sobereva.com/264其中的公式是等价的，自己推导一下便知。

注意：HONO是最高占据自然轨道，LUNO是最低未占据自然轨道，与HOMO、LUMO概念是不一样的。一些实验文章不懂概念、不懂计算，该用HONO和LUNO的时候却用成了HOMO、LUMO，后人应避免犯同样的错误。

不干正事

zjxitcc 发表于 2022-2-15 09:52
（1）不是“残缺的波函数”，是对称破缺波函数，或自旋极化波函数；
（2）你截图中并不是完整的步骤，也不 ...

好的，非常感谢

QuantumicGuy

CASSCF计算双自由基以及双自由基特征的计算
http://sobereva.com/264（http://bbs.keinsci.com/thread-322-1-1.html）

sobereva

如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题“双自由基特征值y的计算”改了，以后务必注意