计算化学公社

标题: 求助计算双自由基特征中自然轨道中占据数出现负值问题 [打印本页]

作者
Author: 不干正事 时间: 2022-2-17 12:50
标题: 求助计算双自由基特征中自然轨道中占据数出现负值问题
大家好，对于一个双自由基体系，想从如下方法计算双自由基特征值：通过自然轨道的n_HOMO和n_LUMO值，代入公式y=1-(2*T)/(1+T^2)，从而求得y值。但是对想要的体系得到对称性破缺波函数后，然后加上关键词guess=(naturalorbitals,read,save,only)后，得到的结果如图显示，出现了负值。不清楚是哪里出现了问题？（由于文件太大，就只放了输入文件和得到开壳层对称残缺波函数的log文件）

作者
Author: 喵星大佬 时间: 2022-2-17 12:52
本帖最后由喵星大佬于 2022-2-17 12:54 编辑

两次用的关键词完全一样么？

作者
Author: 不干正事 时间: 2022-2-17 13:00

喵星大佬发表于 2022-2-17 12:52
两次用的关键词完全一样么？

第一次关键词是#p opt=tight freq ub3lyp/6-311g* guess=(mix,always) ，第二次关键词是#p UB3LYP/6-311G(d) guess=(naturalorbitals,read,save,only)

作者
Author: sobereva 时间: 2022-2-18 07:56
尝试改用Multiwfn产生自然轨道
在Multiwfn中基于fch产生自然轨道的方法与激发态波函数、自旋自然轨道分析实例
http://sobereva.com/403（http://bbs.keinsci.com/thread-9001-1-1.html）

另外注意确保波函数稳定性，产生自然轨道之前做个stable

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