求助，双自由基体系的开壳单重态与双自由基特征、自旋密度的计算问题

ruiwang

    Sob老师以及各位老师好。我最近开始接触双自由基体系的开壳单重态与双自由基特征的计算问题。自己读过了sob老师的相关帖子，http://sobereva.com/82与http://sobereva.com/264。
    我计算的是一个D-A-D型的共轭小分子。双自由基体系的开壳单重态计算流程如下：

1、我先做了单点计算和波函数的稳定性检验。我的关键词如下：

%chk=step-1.chk

#p UM062X/defTZVPP em=GD3 nosymm guess=mix scf(xqc,maxcycle=500) stable=opt

计算结果：在log文件中，我看到了 The wavefunction is already stable. SCF DONE后看到<S**2>=0.4563；

2、读取上一步的波函数，进行结构优化和频率分析。我的关键词如下：

%chk=step-2.chk

#p opt freq UM062X chkbasis em=GD3 nosymm guess=read geom=allcheck

计算结果：这一步没有虚频；

3.读取上一步的波函数和坐标，进行波函数稳定性检验。我的关键词如下：

%chk=step-3.chk

#p UM062X chkbasis em=GD3 nosymm guess=read geom=allcheck stable=opt

计算结果：在log文件中，我看到了 The wavefunction is already stable. SCF DONE后看到<S**2>=0；

我的问题如下：

1、双自由基体系的最终结构和能量是不是Step2优化之后的分子结构、能量；

2、我在计算双自由基特征y和绘制自旋密度图的时候，是不是把应该把上述Step3得到的坐标导出来，用# UM062X/def2TZVPP guess=mix关键词计算体系得到对称破缺波函数，然后得到的chk文件用Multiwfn去分析。（我这样做了，但是计算后log文件里面<S**2>=0，用Multiwfn查看轨道间的重叠积分，结果都是1。自旋密度图也没有分布）；
3、为什么上第三步之后，SCF DONE后看到<S**2>=0。

sobereva

1 是。最后的结构下已经成了闭壳层体系

2 优化完的结构下基态是闭壳层的，双自由基特征不用算也知道是0

参考北京科音中级量子化学培训班（http://www.keinsci.com/KBQC）里讲双自由基体系优化的这页幻灯片