计算化学公社

标题: 请教关于用CDA对金属吸附氮气的理论计算表征 [打印本页]

作者
Author: liyaohx 时间: 2020-4-15 19:22
标题: 请教关于用CDA对金属吸附氮气的理论计算表征
做了实验想要用理论去解释为什么有些金属簇可以和氮气反应有些不可以，比如说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

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

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

作者
Author: liyaohx 时间: 2020-4-16 21:18

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

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

欢迎光临计算化学公社 (http://ccc.keinsci.com/)