请教关于用CDA对金属吸附氮气的理论计算表征

liyaohx

做了实验想要用理论去解释为什么有些金属簇可以和氮气反应有些不可以，比如说FeTi2C2-和FeTi2C3-，这两个只差一个碳原子，都是形成Ti-N键。想着用mutiwen的拓扑分析，找到了金属氮键的BCP点去分析，但是发现不能解释实验（不能反应的电子密度比能反应的簇还大），静电势试了也不行，最后想着用CDA去分析一下氮分子向金属空轨道填电子的情况来分析，但是我发现我得到的值都很小（最大的也就是33号轨道），有点懵，不知道是不是我用的方法不对：
                     **** Result for alpha electrons ****

    Orb.      Occ.          d           b        d - b          r
       1    1.000000   -0.000000    0.000000   -0.000000    0.000000
       2    1.000000   -0.000000   -0.000000    0.000000    0.000000
       3    1.000000   -0.000000   -0.000004    0.000004   -0.000001
       4    1.000000   -0.000000    0.000017   -0.000017    0.000000
       5    1.000000   -0.000001   -0.000002    0.000002   -0.000000
       6    1.000000   -0.000000    0.000000   -0.000000   -0.000000
       7    1.000000   -0.000000    0.000000   -0.000000   -0.000000
       8    1.000000   -0.000027   -0.000785    0.000758   -0.000175
       9    1.000000   -0.000001   -0.000020    0.000019    0.000000
      10    1.000000   -0.000000   -0.000001    0.000001   -0.000000
      11    1.000000   -0.000000   -0.000002    0.000002   -0.000000
      12    1.000000   -0.000006    0.000000   -0.000006    0.000009
      13    1.000000   -0.000003   -0.000005    0.000002    0.000014
      14    1.000000   -0.000000    0.000001   -0.000001   -0.000000
      15    1.000000   -0.000000    0.000000   -0.000000   -0.000000
      16    1.000000    0.000496   -0.000002    0.000498    0.000015
      17    1.000000   -0.000128   -0.000005   -0.000123    0.000004
      18    1.000000   -0.000002    0.000046   -0.000048   -0.000000
      19    1.000000   -0.000001   -0.000007    0.000006    0.000001
      20    1.000000   -0.000003    0.000058   -0.000061    0.000006
      21    1.000000   -0.000003   -0.000392    0.000390   -0.000045
      22    1.000000   -0.000096   -0.006564    0.006468   -0.001505
      23    1.000000   -0.000000   -0.000031    0.000030   -0.000000
      24    1.000000   -0.000000    0.000000   -0.000000    0.000000
      25    1.000000   -0.000001   -0.000064    0.000063   -0.000002
      26    1.000000   -0.000014    0.000080   -0.000094    0.000754
      27    1.000000   -0.000060    0.000415   -0.000475    0.003129
      28    1.000000    0.000000   -0.000204    0.000205    0.000027
      29    1.000000    0.000000   -0.000057    0.000057    0.000008
      30    1.000000   -0.000012   -0.000064    0.000052    0.000853
      31    1.000000   -0.000059    0.000200   -0.000259    0.001743
      32    1.000000    0.010869   -0.000105    0.010974    0.000221
      33    1.000000    0.031681   -0.013312    0.044992    0.019131
      34    1.000000    0.001106   -0.007309    0.008415   -0.027906
      35    1.000000   -0.000201   -0.001319    0.001118    0.000355
      36    1.000000    0.001818    0.000188    0.001630    0.003706
      37    1.000000    0.004064   -0.000687    0.004751    0.000819
      38    1.000000    0.013497    0.001751    0.011746   -0.006674
      39    1.000000    0.000089   -0.000541    0.000629   -0.001703
      40    1.000000    0.000003   -0.000090    0.000093   -0.000025
      41    1.000000    0.001232   -0.001909    0.003141   -0.007000
      42    1.000000   -0.001682    0.000554   -0.002236   -0.006989
      43    1.000000   -0.000455    0.001457   -0.001912    0.001211
      44    1.000000   -0.000041    0.000608   -0.000650   -0.000095
      45    1.000000    0.001343    0.002413   -0.001070   -0.005988
      46    1.000000    0.001625   -0.001961    0.003586   -0.028966
      47    1.000000   -0.000204    0.002938   -0.003141   -0.002503
      48    1.000000    0.000108    0.007963   -0.007855   -0.008933
      49    1.000000    0.000002    0.000429   -0.000427   -0.000830
      50    1.000000   -0.002759    0.017341   -0.020100   -0.000183
      51    0.000000    0.000000    0.000000    0.000000    0.000000

sobereva

我不知道你说的反应具体指什么反应，计算过势垒没有？应当根据机理多方面考虑，未必光是考察个Ti-N键强度就能解释的。

本来金属和N2之间电荷转移量就不会很大，更何况拆分成一堆轨道的贡献，其中贡献最大的那个只有你加粗标注的数量级是非常正常的。

liyaohx

sobereva 发表于 2020-4-16 13:13
我不知道你说的反应具体指什么反应，计算过势垒没有？应当根据机理多方面考虑，未必光是考察个Ti-N键强度就 ...

老师，我说的反应是金属簇对氮气吸附生成金属氮化物这一步的，不过老师你倒是提醒了我