2022, Vol. 36
2. 山西浙大新材料与化工研究院, 山西 太原 030024
2. Shanxi-Zheda Institute of Advanced Materials and Chemical Engineering, Taiyuan 030024, China
丙烷(C3H8)催化氧化脱氢是有效利用天然气和页岩气的重要手段, 二氧化碳(CO2)作为弱氧化剂参与催化反应, 这对于能源利用和解决全球变暖问题具有重要的研究价值[1-2]. 近年来, CO2氧化C3H8反应逐渐成为当前研究的热点[3-4], 其中C3H8的第一个C－H键活化是氧化脱氢过程的关键步骤[1, 5-6]. C3H8初步氧化可以通过C(1)－H或C(2)－H键的断裂进行, 生成的中间体可以再继续裂解. 目前大部分研究表明C3H8的C(2)－H键断裂更具优势, 而且, C3H8的C－H键断裂生成丙基为速率控制步骤, 在C3H8氧化脱氢过程中起着决定性作用[7-9].
目前, 推动C3H8氧化脱氢的关键在于催化剂的设计. 研究表明, 氮化硼(BN)边缘能有效的吸附并选择性活化C3H8的C－H键[10-11], 促进C3H8的氧化脱氢, 但因其大带隙[12]对CO2的吸附活化并不敏感, 而Baei[13]和Qu[14]等研究发现BN双环团簇可以有效地捕获并激活CO2分子. 另外, Yang等[15-18]发现在BN团簇内掺杂少量过渡金属可以调节B-p轨道的占据, 在催化反应中表现出良好的催化活性. 研究表明, 相对于纯的BN材料, 有金属颗粒掺杂的BN类催化体系展示出相对更好的催化性能, 将其用于CO2氧化C3H8反应, 或可同时促进C3H8的C－H键断裂脱氢与CO2的活化加氢.
在目前广泛研究的C3H8脱氢催化剂中, 镍基催化剂可作为C3H8氧化脱氢的有效催化剂, 如介孔NiO在450 ℃时, 催化C3H8氧化的转化率可达25.8%, C3H6选择性为51.3%, 且72 h内保持稳定[19], 但镍基催化剂的NiOx易被还原, 促进C3H8的C－C键断裂同时抑制C－H的活化[20]. 双金属合金相对于孤立的金属基催化剂在催化C3H8选择性脱氢方面则具有更为优异的C－H活化性能, 如Ni/Mo摩尔比为3/1的NiMo/Al2O3催化剂在570 ℃下的C3H6选择性为63.1%[21]; Fe3Ni(111)表面亦更倾向于C3H8的C－H键氧化裂解[22]; PtCu合金催化剂在催化丙烷直接脱氢(PDH)过程中, C3H8的C－C键断裂的能垒(2.0 eV)远高于C－H键(1.54 eV), 表明该体系更利于C－H键活化[23]. Cu作为Ni的邻位金属, 具有与其相似的电子分布状态, 且Cu基催化剂在CO2活化方面有着不可替代的优势, Ni-Cu双金属间的相互作用或能改善活性中心的化学环境, 进而显著提高催化剂的活性[24]. 因此, 为了同时活化C3H8和CO2, 助力C3H8氧化脱氢, NiCu双金属合金催化体系或可作为一种不错的选择. 但Ni和Cu两种金属的含量占比作为影响催化性能的重要因素, 仍有待探究.
在这项工作中, 所涉及到的催化剂模型均为双金属共掺杂的硼氮类富勒烯. 利用VASP软件理论计算了C3H8和CO2在催化剂NixCuy-B24N28(x+y=4, x=1、2、3、4)表面相应的吸附能和反应活化能, 并探讨了CO2氧化C3H8反应速控步骤的吸附以及氧化机理.
1 计算细节及催化剂构型 1.1 计算细节DFT计算均使用Vienna Ab-initio Simulation Package模拟软件包(VASP)进行[25-27], 所用模型为非周期性簇模型, 同时在Z轴方向设置1.5 nm的真空度, 构型优化过程中所有原子均处于驰豫状态且不考虑自旋极化. 运用投影仪增强波(PAW)[28]描述离子核与电子之间的相互作用, 通过广义梯度近似(GGA)[29]处理交换相关函数, 截断能设定为270 eV, 布里渊区k点采样为1×1×1 [30]. 采用CI-NEB方法[31]确定过渡态(TS), 之后进一步计算虚频, 有且仅有一个虚频证实过渡态计算正确.
吸附能定义如下:
| $ \mathrm{E}_{\text {abs }}=\mathrm{E}_{\text {absorbate/ ( catalyst })}-\mathrm{E}_{\text {catalyst }}-\mathrm{E}_{\text {absorbate }} $ | (1) |
其中Eabsorbate/(catalyst)、Ecatalyst、Eabsorbate分别意味着吸附物种吸附在催化剂上的总能量、催化剂本身的能量以及吸附物种本身的能量.
对于在催化剂NixCuy-B24N28表面上发生的A+B→C+D之类的反应, 通过以下公式计算反应热(ΔH)和活化能垒(Ea):
| $ \varDelta H=\mathrm{E}_{(\mathrm{C}+\mathrm{D}) /(\text { catalyst })}-\mathrm{E}_{(\mathrm{A}+\mathrm{B}) /(\text { catalyst })} $ | (2) |
| $ \mathrm{E}_{\mathrm{a}}=\mathrm{E}_{\mathrm{TS} /(\text { catalyst }}-\mathrm{E}_{(\mathrm{A}+\mathrm{B}) /(\text { catalyst })} $ | (3) |
其中, E(A+B)/(catalyst)表示催化剂表面共吸附的反应物即A和B的总能量, E(C+D)/(catalyst)表示产物即被吸附的C和D的总能量, ETS/(catalyst)是反应过渡态的能量.
p带中心的公式如下所示[32]:
| $ \varepsilon \mathrm{p}=\frac{\int_{+\infty}^{-\infty} \mathrm{EP}(\mathrm{E}) \mathrm{dE}}{\int_{+\infty}^{-\infty} \mathrm{P}(\mathrm{E}) \mathrm{dE}} $ | (4) |
其中, E代表能量, P(E)代表p电子密度.
1.2 催化剂构型我们所用催化剂为双金属共掺杂的硼氮类富勒烯, 这是一个(BN)28笼状模型, 通过电弧熔炼合成在其顶部以4个金属原子取代掺杂4个B原子的结构[33-35]. 以NiCu3-B24N28为例, 模型的主视图和俯视图如图 1所示. 深蓝色为氮原子, 粉色为硼原子, 黄色为铜原子, 浅蓝色为镍原子.
|
图 1 NiCu3-B24N28催化剂的主视图和俯视图 Fig.1 Front view and top view of the NiCu3-B24N28 catalyst |
(BN) 28富勒烯是由28个硼原子和28个氮原子组成的四元环与六元环依次相连的稳定笼结构. 我们计算了金属于不同位点取代掺杂(BN)28簇模型的形成能, 搭建了如图 1所示的以4个金属粒子取代4个B原子的稳定催化剂模型, 之后确定了C3H8、CO2及其中间体的吸附态, 包括CH3CH2CH3、CH2CH2CH3、CH3CHCH3、OCO、OCHO和HOCO. 各物种在NiCu3-B24N28表面的稳定吸附构型图如图 2所示, 而在其它NixCuy-B24N28表面的吸附构型与在NiCu3-B24N28表面类似, 与之对应的吸附能列于表 1.
|
图 2 C3H8及相关物种在NiCu3-B24N28表面的稳定构型图 Fig.2 The stable configuration of C3H8 and related species on NiCu3-B24N28 surface |
|
表 1 C3H8及相关物种在NixCuy-B24N28表面的吸附能(eV) Table 1 Adsorption energy(eV) of C3H8 and related species on the surface of NixCuy-B24N28 |
C3Hx(x = 7~8)的吸附: C3H8吸附在H1位点的上方与Cu-N-Cu组成的三元环相互作用, 其在NixCuy-B24N28表面吸附最稳定的结构是C－C键, 几乎与催化剂表面的Cu－N键平行, 且为物理吸附. CH2CH2CH3通过C3H8伯碳的C－H键断裂, 它的最稳定吸附状态是通过伯碳吸附在T1位点, 通过表 1可以看出, 由于失去一个H, 其在催化剂表面的吸附能比C3H8大得多, 是化学吸附. CH3CHCH3由C3H8仲碳的C－H键断裂生成, 通过仲碳吸附在催化剂表面的T1位点, CH3CHCH3与CH2CH2CH3类似, 都是化学吸附.
HxCO2(x = 0~1)的吸附: CO2吸附在T2位点的上方, 其在NixCuy-B24N28表面吸附最稳定的结构是CO2的C原子位于T1位点N原子的上方, 且整个CO2分子与催化剂表面呈现15°到45°的倾斜角度, 是物理吸附. OCO的C原子被C3H8脱掉的H原子攻击生成OCHO中间体, 通过O原子与催化剂表面N原子之间的相互作用, 化学吸附在NixCuy-B24N28表面. OCO通过与C3H8脱掉的H原子结合生成HOCO, 中间体HOCO与OCHO类似, 都是化学吸附.
2.2 C3H8和CO2在NixCuy-B24N28表面的初步氧化过程利用VASP软件计算了NixCuy-B24N28表面CO2氧化C3H8反应速控步骤的反应热和活化能. 在NiCu3-B24N28催化剂表面, CO2氧化C3H8反应速控步骤的初态(IS)、过渡态(TS)、终态(FS)的构型如图 3所示. 其它NixCuy-B24N28表面与NiCu3-B24N28的IS、TS、FS的构型类似, 相应的反应热和活化能的计算结果在表 2列出.
|
图 3 NiCu3-B24N28表面CO2氧化C3H8速控步骤的IS、TS、FS的稳定构型图 Fig.3 The stable configuration of IS, TS and FS of the rate-controlling step of the oxidation of C3H8 by CO2 on the NiCu3-B24N28 surface |
|
|
表 2 NixCuy-B24N28表面CO2氧化C3H8速控步骤的相关反应热和活化能 Table 2 Related reaction heat and activation energy of the rate-controlling step of the oxidation of C3H8 by CO2 on the surface of NixCuy-B24N28 |
R1: C3H8+CO2→CH2CH2CH3+HOCO. C3H8在NiCu3-B24N28表面的初始氧化始于C3H8和CO2分子的吸附. C3H8的C－C键与表面平行, CO2与表面呈45°倾斜角, 作为反应物的初始构型. C3H8的初始氧化有六条路径, 分别是C(1)－O(1)、C(1)－O(2)、C(1)－C、C(2)－O(1)、C(2)－O(2) 和C(2)－C键结合路径. 该反应表示的是C3H8和CO2开始初步氧化, 经过C3H8伯碳脱氢与CO2的O结合, 即通过C(1)－O(1)路径生成CH2CH2CH3和HOCO. CH2CH2CH3和HOCO通过C原子分别与催化剂表面的T1、T2位点成键, 该反应放热-2.09 eV, 活化能为1.73 eV.
R2: C3H8+CO2→CH2CH2CH3+OCOH. C3H8伯碳脱氢也可以和CO2的另一个活性O结合, 即通过C(1)－O(2)路径生成CH2CH2CH3和OCOH, 该反应放热-2.10 eV, 需克服1.85 eV的能垒. 其活化能相对于C(1)－O(1)路径较高, 因此不支持该路径优先反应.
R3: C3H8+CO2→CH2CH2CH3+OCHO. C3H8伯碳脱氢还可以与CO2的C结合, 即通过C(1)－C路径跨越1.68 eV的能垒, 生成CH2CH2CH3和OCHO. 可见C(1)－C路径所需的活化能略低于C(1)－O(1)路径, 两者之间存在较激烈的竞争(1.68、1.73 eV), 但前者较容易.
R4: C3H8+CO2→CH3CHCH3+HOCO. C3H8还可以通过仲碳脱氢与氧化剂CO2的O结合, 该过程中反应放热为-2.19 eV, H原子需克服1.96 eV的能垒, 通过C(2)－O(1)路径生成CH3CHCH3和HOCO. 将其与伯碳脱氢的3条路径对比, 发现R4路径发生反应的概率较低.
R5: C3H8+CO2→CH3CHCH3+OCOH. C3H8通过仲碳脱氢经TS5形成CH3CHCH3和OCOH的活化能为1.42 eV, 反应放热-2.09 eV. 根据表 2中的活化能数据, 可以发现在NiCu3-B24N28表面, C(2)－O(2)路径反应所需活化能最低, 且低于同路径下的Ni2Cu2-B24N28、Ni3Cu-B24N28以及Ni4-B24N28表面, 因此可作为最有可能发生的路径, 这与之前的研究结果一致[36].
R6: C3H8+CO2→CH3CHCH3+OCHO. C3H8仲碳脱氢的C(2)－C路径也可能发生, 该过程放热-0.10 eV, 活化能为2.37 eV. 生成CH3CHCH3和OCHO的活化能高于生成CH3CHCH3和OCOH的活化能, 说明C(2)－O(2)路径更有利.
通过图 4可以发现, NixCuy-B24N28催化CO2氧化C3H8反应中C3H8的第一个C－H键断裂的活化能随着Cu占比的增大逐渐降低, 这说明催化活性与催化体系中掺杂金属原子的相对数量有一定的标度关系. 在NixCuy-B24N28催化剂表面, C3H8的C－H键的活化以及CO2的C－O键的活化可能与Cu驱使的双金属向BN层的电子转移数量有关.
|
图 4 NixCuy-B24N28表面C3H8第一个C－H键活化关系图 Fig.4 The relationship between the activation of the first C－H bond for C3H8 on the surface of NixCuy-B24N28 |
我们计算了NixCuy-B24N28表面C(2)－O(2)路径中丙烷C(2)原子及催化剂N原子的p轨道的局域态密度(PDOS)及其p带中心, 图 5中0点为费米能级, 不同颜色的曲线分别代表C原子(黑色)和N原子(红色)的PDOS, 垂直实线分别表示其p带中心. 由图 5(a)-(d)可知, 在双金属掺杂催化体系中, C3H7中间体C原子与吸附位点N原子的p轨道的重叠程度逐渐增大, 即C原子与N原子p带中心差值分别为NiCu3-B24N28(0.52 eV)、Ni2Cu2-B24N28(0.82 eV)、Ni3Cu-B24N28(0.96 eV)和Ni4-B24N28(1.06 eV). 结果表明, 随着Cu含量占比的逐渐升高, 丙烷C原子与催化剂表面N原子之间的电子转移量逐渐增大, 使得催化剂对丙烷C－H键的活化作用逐渐增强, 反应所需的活化能逐渐降低. 该分析结果与我们所得活化能变化规律相匹配.
|
图 5 NixCuy-B24N28表面C原子与N原子p轨道的PDOS及p带中心 Fig.5 The PDOS and p-band center of p orbital of carbon atom and nitrogen atom on the surface of NixCuy-B24N28 |
通过VASP软件计算了CO2氧化C3H8速控步骤中6条可能路径下的反应热和活化能, 研究了C3H8和CO2在NixCuy-B24N28(x+y=4, x=1、2、3、4)催化剂表面吸附和速控步骤的反应机理. 根据计算结果得出以下结论: (1) C3H8和CO2在NixCuy-B24N28表面是物理吸附, C3H8+CO2→CH3CHCH3+OCOH是最有利的路径. (2)在NixCuy-B24N28(x+y=4, x=1、2、3、4)催化体系中, Cu含量的增加有助于提高CO2氧化C3H8速控步骤的活性, 其在不同催化剂表面的活化能顺序是NiCu3-B24N28(1.42 eV)、Ni2Cu2-B24N28(1.57 eV)、Ni3Cu-B24N28(1.62 eV)、Ni4-B24N28(1.75 eV).
| [1] |
a. Zhang Tao(张涛), Liu Zhi-cheng(刘志成), Yang Wei-min(杨为民). Progress in the catalytic conversion of light alkanes with carbon dioxide(低碳烷烃与二氧化碳催化转化研究进展)[J]. Sci Sinchim(中国科学: 化学), 2021, 51(2): 154-164. b. Zhang Jin-he(张金赫), Zhang Chang-wu(张昌武), Wang Xin-ge(王新鸽), et al. Self-hydrolysis-driven syn-thesis of nanocrystalline aluminum nitride supported platinum-zinc catalyst for propane dehydrogenation(自水解驱动PtZn/纳米AlN催化剂的合成及其丙烷脱氢性能研究)[J]. J Mol Catal(China)(分子催化), 2021, 35(5): 395-403. c. Gao Xi-ran(高希然), Dilinur·Ali(迪丽努尔), Fang Ya-ping(方亚平), et al. FTIR studies of the metal-sup-port interactions on Au/HZSM-5: Acidity and dehydr-ogenation of new Si-O (H) –Au species(Au/HZSM-5催化剂上Si-O(H)-Au的酸性和脱氢性能研究)[J]. J Mol Catal(China)(分子催化), 2021, 35(3): 226-234. |
| [2] |
Tian Hai-feng(田海锋), Liao Jian-kang(廖建康), Cha Fei(查飞), et al. Thermodynamic simulation of propane coupling with carbon dioxide to propylene(丙烷与二氧化碳耦合制丙烯的热力学模拟研究)[J]. J North-west Normal Univer (Nat Sci)(西北师范大学学报(自然科学版)), 2019, 55(2): 94–99.
|
| [3] |
Atanga M A, Rezaei F, Jawad A, et al. Oxidative dehydro-genation of propane to propylene with carbon dioxide[J]. Appl Catal B-Environ, 2018, 220: 429–445.
DOI:10.1016/j.apcatb.2017.08.052 |
| [4] |
Kang J, Czaja A D, Guliants V V. Carbon dioxide as feedstock in selective oxidation of propane[J]. Eur J Inorg Chem, 2017, 2017(40): 4757–4762.
DOI:10.1002/ejic.201701049 |
| [5] |
Coperet , Christophe . C－H bond activation and organ-ometallic intermediates on isolated metal centers on oxide surfaces[J]. Chem Rev, 2010, 110(2): 656–680.
DOI:10.1021/cr900122p |
| [6] |
Zhao Z J, Chiu C C, Gong J. Molecular understandings on the activation of light hydrocarbons over heterogeneous catalysts[J]. Chem Sci, 2015, 6(8): 4403–4425.
DOI:10.1039/C5SC01227A |
| [7] |
Xu Zhi-kang(徐志康), Huang Jia-lu(黄佳露), Wang Ting-hai(王廷海), et al. Advances in catalysts for prop-ane dehydrogenation to propylene(丙烷脱氢制丙烯催化剂的研究进展)[J]. Chem Indus Eng Prog(化工进展), 2021, 40(4): 1893–1916.
|
| [8] |
Si C W, Lian Z, Olanrele S O, et al. Revealing the origin of the reactivity of metal-free boron nitride catalysts in oxidative dehydrogenation of propane[J]. Appl Surf Sci, 2020, 519: 146241.
DOI:10.1016/j.apsusc.2020.146241 |
| [9] |
Zhao Zhi-bing(赵志冰). A DFT study on the mechanism of BN catalyzed propane oxidative dehydrogenation(BN催化丙烷氧化脱氢反应机理的DFT研究)[D]. Dalian(大连): Dalian University of Technology(大连理工大学), 2018.
|
| [10] |
Wang L C, Wang Y Y, Zhang R, et al. Edge-zctivating CO2-mediated ethylbenzene dehydrogenation by a hierarchical porous BN catalyst[J]. ACS Catal, 2020, 10(12): 6697–6706.
DOI:10.1021/acscatal.0c00070 |
| [11] |
Rajbanshi B, Saha S, Fricke C, et al. Oxidative dehydr-ogenation of propane on the oxygen adsorbed edges of boron nitride nanoribbons[J]. Catal Sci Technol, 2020, 10(15): 5181–5195.
DOI:10.1039/D0CY01031F |
| [12] |
Arenal R, Ferrari A C, Reich S, et al. Raman spectro-scopy of single-wall boron nitride nanotubes[J]. Nano Lett, 2006, 6(8): 1812–1816.
DOI:10.1021/nl0602544 |
| [13] |
Baei M T, Peyghan A A, Bagheri Z. A DFT study on CO2 interaction with a BN nano-cage[J]. B Korean Chem Soc, 2012, 33(10): 3338–3342.
DOI:10.5012/bkcs.2012.33.10.3338 |
| [14] |
Qu M, Qin G, Du A, et al. B80 fullerene: A promising metal-free photocatalyst for efficient conversion of CO2 to HCOOH[J]. J Phys Chem C, 2019, 123(39): 24193–24199.
DOI:10.1021/acs.jpcc.9b07562 |
| [15] |
Yang X W, Zhou S, Huang S L, et al. Metal-encapsulated boron nitride nanocages for solar-driven nitrogen fixation[J]. J Phys Chem C, 2020, 124(43): 23798–23806.
DOI:10.1021/acs.jpcc.0c07968 |
| [16] |
Nigam S, Majumder C. CO oxidation by BN fullerene cage: Effect of impurity on the chemical reactivity[J]. ACS Nano, 2008, 2(7): 1422–1428.
DOI:10.1021/nn8001455 |
| [17] |
Wen S H, Deng W Q, Han K L. Endohedral BN metall-ofullerene M@B36N36 complex as promising hydrogen storage materials[J]. J Phys Chem C, 2008, 112(32): 12195–12200.
DOI:10.1021/jp801893f |
| [18] |
Zhao J J, Du Q Y, Zhou S, et al. Endohedrally doped cage clusters[J]. Chem Rev, 2020, 120(17): 9021–9163.
DOI:10.1021/acs.chemrev.9b00651 |
| [19] |
Li J H, Wang C C, Huang C J, et al. Mesoporous nickel oxides as effective catalysts for oxidative dehydrogenation of propane to propene[J]. Appl Catal A-Gen, 2010, 382(1): 99–105.
DOI:10.1016/j.apcata.2010.04.034 |
| [20] |
Saelee T, Namuangruk S, Kungwan N, et al. Theoretical insight into catalytic propane dehydrogenation on Ni(111)[J]. J Phys Chem C, 2018, 122(26): 14678–14690.
DOI:10.1021/acs.jpcc.8b03939 |
| [21] |
Wang T, Cui X J, Winther K T, et al. Theory-aided discovery of metallic catalysts for selective propane dehydrogenation to propylene[J]. ACS Catal, 2021, 11(10): 6290–6297.
DOI:10.1021/acscatal.0c05711 |
| [22] |
Gomez E, Kattel A S, Yan B, et al. Combining CO2 redu-ction with propane oxidative dehydrogenation over bimetallic catalysts[J]. Nat Commun, 2018, 9(1): 1398–1404.
DOI:10.1038/s41467-018-03793-w |
| [23] |
Sun S J, Sun G D, Pei C L, et al. Origin of performances of Pt/Cu single-atom alloy catalysts for propane dehydr-ogenation[J]. J Phys Chem C, 2021, 125(34): 18708–18716.
DOI:10.1021/acs.jpcc.1c04295 |
| [24] |
Tian Ya-ni(田亚妮), Zhao Fang(赵芳), Wang Chang-zhen(王长真), et al. Effect of Cu additive on the hydrog-enation of 1, 4-butynedio over Ni-Cu/SiO2 bimetallic catalyst(Cu对Ni-Cu/SiO2催化1, 4-丁炔二醇加氢性能的影响)[J]. J Mol Catal(China)(分子催化), 2019, 33(2): 132–139.
|
| [25] |
Hafner J. Ab-initio simulations of materials using VASP: Density-functional theory and beyond[J]. J Comput Chem, 2008, 29(13): 2044–2078.
DOI:10.1002/jcc.21057 |
| [26] |
Sun G, Kürti J, Rajczy P, et al. Performance of the vienna ab initio simulation package (VASP) in chemical applications[J]. J Mol Struc-Theochem, 2015, 624(1/3): 37–45.
|
| [27] |
Kresse G, Hafner J. Ab initio molecular-dynamics simu-lation of the liquid-metal-amorphous-semiconductor transition in germanium[J]. Phys Rev B, 1994, 49(20): 14251–14269.
DOI:10.1103/PhysRevB.49.14251 |
| [28] |
Kresse G, Joubert D. From ultrasoft pseudopotentials to the projector augmented-wave method[J]. Phys Rev B, 1999, 59(3): 1758–1775.
DOI:10.1103/PhysRevB.59.1758 |
| [29] |
Perdew J P, Burke K, Ernzerhof M. Generalized gradient approximation made simple[J]. Phys Rev Lett, 1998, 77(18): 3865–3868.
|
| [30] |
Wang W N, Ren R P, Lv Y K. DFT study on the mecha-nism of methanol to methyl formate on the M@C16B8 surface[J]. Mater Today Commun, 2021, 26: 102090.
DOI:10.1016/j.mtcomm.2021.102090 |
| [31] |
Rusydi F, Madinah R, Puspitasari I, et al. Teaching reaction kinetics through isomerization cases with the basis of density-functional calculations[J]. Biochem Mol Biol Edu, 2021, 49(2): 216–227.
DOI:10.1002/bmb.21433 |
| [32] |
Dickens C F, Montoya J H, Kulkarni A R, et al. An electronic structure descriptor for oxygen reactivity at metal and metal-oxide surfaces[J]. Surf Sci, 2019, 681: 122–129.
DOI:10.1016/j.susc.2018.11.019 |
| [33] |
Feng L, Lu Y H, Kong J, et al. Theoretical studies on the structure and properties of BN clusters (BN)n and endohedral metallo-BN clusters M@(BN)n[J]. Comput Theor Chem, 2011, 964(1/3): 56–64.
|
| [34] |
Oku T, Suganuma K. High-resolution electron micro-scopy and structural optimization of c, b n and fe@b n clusters[J]. Diam Relat Mater, 2001, 10(3/7): 1205–1209.
|
| [35] |
Li R Y, Wang Y. Modification of boron nitride nanocages by titanium doping results unexpectedly in exohedral complexes[J]. Nat Commun, 2019, 10(1): 4908–4916.
DOI:10.1038/s41467-019-12877-0 |
| [36] |
Liu Y M, Xie S H, Cao Y, et al. Synthesis of novel cage-like mesoporous vanadosilicate and its efficient performance for oxidation dehydrogenation of propane[J]. J Phys Chem C, 2010, 114(13): 5941–5946.
DOI:10.1021/jp909749w |