单原子对称性破缺态如何做自旋布居分析

ggdh

以对称性破缺单重态碳卡宾为例
这个碳上有两个单电子，一个自旋向上，一个自旋向下，但是这两个电子没有配对，处于不同轨道上。
如果做普通的自旋布居分析，那这个碳上的大概是0左右。
如何用自旋布居来表示这种单个碳上不配对的情况呢？
研究了下sob大大的博文：
http://sobereva.com/353
觉得用Hirshfeld/Becke盆分别对正负自旋密度进行积分似乎可以。
但是博文中给出的例子，没有分别对正负密度进行积分。不知如何才能实现？
或者是否还有其他更好的表征手段？

sobereva

在function.f90中的function userfunc里面合适的位置插入这一段
if (iuserfunc==700) then
tmpval=fspindens(x,y,z,1)
if (tmpval>0) userfunc=tmpval
end if
代表，如果settings.ini里的iuserfunc设成了700，那么利用15-1-100来积分的话，积分的都是每个原子上自旋密度为正的部分。类似地，可以添加代码积分自旋密度为负的部分。

ggdh

sobereva 发表于 2018-1-6 07:27
在function.f90中的function userfunc里面合适的位置插入这一段
if (iuserfunc==700) then
tmpval=fspind ...

谢谢sob
把spinnatural orbital 弄出来后发现了另外一种方法，可以不用重新编译Multiwfn。
首先用Multiwfn打开SNO的wfn文件，用：

1. 6 26 轨道编号

复制代码

把所有本征值为负（beta自旋）的轨道占据数改为0（这些轨道排在最后，比如一共900个轨道，），基本上700-900号轨道是负的。
然后用：

1. 15 1 1

复制代码

把alpha电子自旋的布居弄出来。
再用同样的方法对把beta电子自旋的布居弄出来就行了。