芳香体系双自由基特征计算问题

Sway · 发表于 Post on 2021-2-26 23:19:33

最近在做一个大共轭芳香性体系的双自由基特征计算，遇到了以下问题：1.想先计算对称破缺态，做了以下输入：

然后输出为：

不知道是否得到对称破缺的波函数了呢？

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

granvia · 发表于 Post on 2021-2-27 00:00:12

直接看下g16输出中的Mulliken自旋密度

sobereva · 发表于 Post on 2021-2-27 00:39:28

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

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

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

Sway · 发表于 Post on 2021-2-27 14:25:55

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

...

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] 芳香体系双自由基特征计算问题

浏览过的版块