spin density和spin population做出来不一样

ysqyxfxrg

用Multiwfn做一些配合物结构的spin density和spin population分析如下：(isosurface 0.003)

结构一：

Atom     Alpha pop.   Beta pop.    Spin pop.     Atomic charge

    1(Au)      9.67754      9.08654      0.59101         0.23592

    2(Au)      9.75340      9.44291      0.31050        -0.19631

结构二：

Atom     Alpha pop.   Beta pop.    Spin pop.     Atomic charge

    1(Au)      9.47686      9.04291      0.43395         0.48023

    2(Au)      9.90894      9.49685      0.41209        -0.40580

想问一下，结构一和结构二的spin population在Au(1)和Au(2)上都是差不多的，为何spin density图上显得只有Au(1)有alpha spin？按照spin density图的样子，应该是Au(1)和Au(2)的spin population相差很大。为啥spin density和spin population趋势不一样？

sobereva

让Multiwfn在图上把原子序号显示出来，别人根本不知道你说的是哪个原子
并且说清楚计算细节，什么级别得到的波函数、在Multiwfn里用什么方法得到的自旋布居

ysqyxfxrg

老师好，用的计算方法如下：
! CAM-B3LYP def2-TZVP def2/J RIJCOSX tightopt SlowConv miniprint nopop
将产生的molden文件放入Multiwfn。

绘制使用的方法是http://sobereva.com/353博文里的，

自旋密度方法：
5    //计算格点数据
5    //自旋密度
3    //用高质量格点

自旋布居方法：
7   //布居分析
5   //Mulliken分析
1

Isosurface 0.003

Atom      Alpha pop.   Beta pop.    Spin pop.     Atomic charge
     1(Au)      9.67754      9.08654      0.59101         0.23592
     2(Au)      9.75340      9.44291      0.31050        -0.19631

Isosurface 0.003

Atom      Alpha pop.   Beta pop.    Spin pop.     Atomic charge
     1(Au)      9.47686      9.04291      0.43395         0.48023
     2(Au)      9.90894      9.49685      0.41209        -0.40580



然后要比较的是两个文件里的两个Au的自旋密度和自旋布居，上图的两个Au1的电子密度图上的电子密度都远大于Au2，但是自旋布居却两个Au差不多，想问一下为啥，谢谢老师