计算化学公社

标题: 求助双自由基特征值y的计算问题 [打印本页]

作者
Author: 不干正事 时间: 2022-2-15 08:46
标题: 求助双自由基特征值y的计算问题
大家好，对于一个双自由基体系，已经用guess=mix关键词得到了残缺的波函数，然后根据之前得到的一个方法，用图中显示的一个方法计算得到的，感觉不是特别正确。然后查相关文献也看到了一个公式来计算y值，但不清楚公式中的各项需要代入什么？

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

#p UBHandHLYP/6-31G(d) nosymm guess=mix stable=opt pop=no
title
0 1
H 0.0 0.0 0.0
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，后人应避免犯同样的错误。

作者
Author: 不干正事 时间: 2022-2-15 10:29

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

好的，非常感谢

作者
Author: QuantumicGuy 时间: 2022-2-15 12:00
CASSCF计算双自由基以及双自由基特征的计算
http://sobereva.com/264（http://bbs.keinsci.com/thread-322-1-1.html）

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

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