计算化学公社

标题: 芳香体系双自由基特征计算问题 [打印本页]

作者
Author: Sway 时间: 2021-2-26 23:19
标题: 芳香体系双自由基特征计算问题
最近在做一个大共轭芳香性体系的双自由基特征计算，遇到了以下问题：1.想先计算对称破缺态，做了以下输入：
(, 下载次数 Times of downloads: 40)
然后输出为：
(, 下载次数 Times of downloads: 57)
不知道是否得到对称破缺的波函数了呢？

2.根据帖子http://sobereva.com/264进行UNO的分析过程中，轨道的occupied接近2为双占据，1-2之间为自由基的特征，结果我的计算结果如下,occupied的数值也太大了。不知道是哪里出了问题，求指教

(, 下载次数 Times of downloads: 47)

作者
Author: granvia 时间: 2021-2-27 00:00
直接看下g16输出中的Mulliken自旋密度

作者
Author: sobereva 时间: 2021-2-27 00:39
1 如果最终的波函数的<S**2>大于1，或者用Multiwfn绘制自旋密度图发现不是处处为0，就说明是对称破缺单重态。

直接看Gaussian输出的Mulliken原子自旋布居不是很稳妥，对有的体系原子上的alpha单电子和beta单电子数恰好相同，会误判。

2 操作方式有误。先把博文里的例子重复一遍确认操作正确。

作者
Author: Sway 时间: 2021-2-27 14:25

sobereva 发表于 2021-2-27 00:39
1 如果最终的波函数的大于1，或者用Multiwfn绘制自旋密度图发现不是处处为0，就说明是对称破缺单重态。

...

好的，我再确认下2是不是操作错误。谢谢sob老师！！！

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