含Fe的离子团簇自旋多重度确定求助

wzkchem5

wuyeyangguang 发表于 2021-9-22 04:29
谢谢大佬指点。根据提示做了自旋布居分析，由于计算中加了弥散函数，因此用了Becke方法。闭壳层单重态下 ...

我说的看布居数的意思，不是为了确认总自旋多重度对不对，而是确认波函数对不对。
同一个自旋多重度下可能有不止一个稳定波函数，其中只有一个是对的。而确定一个波函数是不是对的，有以下几种方法：
（1）做布居分析，看每个原子的电荷、自旋布居是否与预期相符；
（2）做自旋密度图、前线轨道图，看是否与预期相符；
（3）把所有可能的波函数都找到，看现在这个波函数是不是能量最低的。
其中（3）最严格，但是一般无法把所有可能的波函数都找全（找全了也很难知道是不是真的找全了），即使能找全也极其费时间。所以推荐（1）和（2），其中（1）更快一些，因为打开输出文件就能看。

wuyeyangguang

wzkchem5 发表于 2021-9-22 15:08
我说的看布居数的意思，不是为了确认总自旋多重度对不对，而是确认波函数对不对。
同一个自旋多重度下可 ...

感谢如此细致的回复。三种情况的原子自旋布居分别是；
Atomic space   Value（自旋多重度11） Value（以闭壳层形式下的1）   Value（自旋极化1）
     1(C )            0.00072589              0.00029870              0.00023641
     2(C )            0.00057724             -0.00010625              0.00030127
     3(C )            0.00057826              0.00030044              0.00005522
     4(H )            0.00026981              0.00003601              0.00009816
     5(H )            0.00030116              0.00016551              0.00013675
     6(H )            0.00028973              0.00011563              0.00008348
     7(N )            0.00059526              0.00005600              0.00042419
     8(C )            0.00216275             -0.00069811             -0.00009637
     9(H )            0.00002779              0.00003042              0.00003833
    10(H )            0.00053614             -0.00045695             -0.00039877
    11(H )            0.00031066             -0.00002379              0.00010866
    12(N )            0.00166461              0.00057274              0.00021315
    13(C )            0.00150417             -0.00062203             -0.00028748
    14(H )            0.00025718             -0.00013545             -0.00007625
    15(H )            0.00009877             -0.00001760              0.00002504
    16(C )            0.00118658              0.00039604              0.00007850
    17(H )            0.00019709              0.00009004             -0.00008130
    18(H )            0.00010536              0.00003980              0.00001158
    19(H )            0.00023553              0.00012174              0.00015112
    20(Cl)            0.26622250             -0.00331347              0.00293561
    21(Cl)            0.26588917             -0.01117699              0.01270492
    22(Cl)            0.27813788             -0.01575011             -0.04100828
    23(Cl)            0.23223735             -0.00658771             -0.04290390
    24(Cl)            0.27154198              0.01759983              0.04046718
    25(Cl)            0.27416173              0.01671542             -0.01101183
    26(Cl)            0.22854347              0.01048725              0.04383617
    27(Cl)            0.26894925              0.00186959             -0.00096227
    28(Fe)            3.95300437             -0.96076845             -0.92984403
    29(Fe)            3.94966517              0.95075992              0.92476451
采用Multiwfn绘制的自旋密度图主要集中在Fe附近。其中11的时候，Fe周围全部为绿色(为正)；对称破缺态的时候两个Fe原子附近都是既有α电子也有β电子，也都集中在Fe附近。奇怪的是按闭壳层方式算，自旋密度也不是处处为0，和对称破缺情况差不多。

wzkchem5

wuyeyangguang 发表于 2021-9-22 11:19
感谢如此细致的回复。三种情况的原子自旋布居分别是；
Atomic space   Value（自旋多重度11） Value（ ...

这个就属于没有收敛到正确的自旋极化态。因为Cl-是弱场配体，FeCl4-是高自旋d5体系，每个FeCl4-都应该有5个自旋，其中一个FeCl4-的自旋布居加起来是+5左右，另一个是-5左右