2022, Vol. 36
2. 长江师范学院 电子信息工程学院 超常配位键工程与新材料技术重庆市重点实验室, 重庆 408100
2. Key Laboratory of Extraordinary Bond Engineering and Advanced Materials Technology of Chongqing, School of Electronic Information Engineering, Yangtze Normal University, Chongqing 408100, China
臭氧（O3）的催化分解是环境污染控制领域的重要反应. 一方面O3本身是一种污染物, 会导致光化学烟雾, 影响人体肺功能, 引起中枢神经障碍[1-2], 而催化分解是去除O3的有效方法[3-7]之一. 另一方面, 在催化O3深度氧化和等离子体-催化技术降解有机物过程中, 高压放电产生的O3在催化剂表面解离成氧化能力更强的活性O原子和OH自由基, 从而对反应物和放电过程中未被完全降解的产物进行深度氧化, 达到去除有机物的目的[8-9]. 近年来, 人们对催化分解O3反应进行了广泛研究, 过渡金属氧化物由于成本低、毒性小等优点, 是进行O3分解的重要催化剂之一, 逐渐成为研究热点[10-11]. 氧化铜（CuO）广泛应用于催化O3分解过程, 例如SiO2负载金属氧化物分解O3的效率依次为: Cu/SiO2 > Co/SiO2 > Ni/SiO2 > Fe/SiO2 > Mn/SiO2[12]. CuO、CuO/沸石等催化剂能够有效催化O3氧化空气中的甲醛[13-14].Cu2O-CuO-Cu（OH）2分级纳米结构在室温下表现出100%的O3分解效率, 并具有高稳定性和抗水性, 纳米结构的表面形态和化学成分在催化O3分解中起着重要作用, 纳米结构中的离域空穴能促进中间产物氧阴离子和催化剂表面之间的电子转移[15].
O3在过渡金属氧化物表面吸附时, 表面可以作为一个能隙较小的电子供体, 电子从金属离子转移到O3, 并分解为一个O2分子和一个O原子, O原子可与另一个O3分子反应形成自由氧分子和超氧自由基O2-[16]. Wu等[17]利用原位红外检测到O3在CaO、ZnO、γ-Al2O3、CuO和α-Fe2O3表面吸附和分解后的物种有O3-、O和O2-. 在Al2O3（110）和Cu2O（111）表面, O3可自发地在Al和Cu离子上解离成氧分子和О原子, 在Al2O3（110）表面的解离能为-1.02 eV[18-19]. 在MnO2表面, O3末端的O原子通过氢键吸附在表面酸位点上, 另一末端O原子插入氧空位位点, 通过电子转移形成氧分子, 留下O2-物种, 同时酸位点再生; O3末端的O原子形成的O2-发生作用, 形成O22-同时释放氧气, 最后O22-脱附形成氧气, 氧空位点再生[20]. Oison等[21]对比了通过GGA+U和LDA+U计算得到的O3在CuO（111）表面的吸附能, LDA+U得到的吸附能的值偏大, GGA+U计算得到O3和表面之间离子型相互作用更多, 共价型相互作用更少.
目前, 尽管近年来人们对催化分解O3反应进行了一定的实验研究和理论计算, 但是O3分子在催化剂表面的分解路径尚不清晰, 对于氧化还原反应过程中催化剂氧空位和表面吸附氧的生成机理的研究还不够充分. 因此, 我们基于GGA+U计算, 研究O3在完整和具有氧空位的CuO（111）表面吸附的吸附位、吸附能、吸附结构和电子结构特征, 比较了O3在完整表面和具有氧空位的表面分解的路径和能垒, 分析该过程中氧空位和表面吸附氧的产生, 以及对O3吸附和分解的影响, 以便深入理解催化分解O3的机理和O3分解对于催化剂氧化还原性能的影响.
1 计算方法及模型计算采用VASP软件包[22-23], 基于广义梯度近似（GGA-PBE）[24]和投影增强波（PAW）[25]方法, 截断动能为520 eV. 优化后的CuO的晶格参数为a=0.47 nm、b=0.34 nm、c=0.51 nm和β=99.54°, 与之前的计算结果基本一致[26]. 根据以往的研究, CuO（111）表面具有较低的表面能, 是最稳定的表面[27]. 采用p（4×2）超胞建模, 三层O-Cu-O结构, z方向上增加1.5 nm的真空层[28-29].
在GGA+U计算中, 设置Cu的Hubbard参数U=7.5 eV, 可以很好地处理Cu 3d轨道的强电子相关性[30]. 表面结构弛豫收敛的判定标准为Hellmann-Feynman力小于0.2 eV/nm[27, 31]. 分别使用2×2×1和7×7×1的Monkhorst-Pack网格进行弛豫和态密度计算. 采用CI-NEB方法搜索过渡态[32].
吸附能（Eads, kJ/mol）和差分电荷密度（CDD, Δρ, 1000 e/（nm）3）的计算方法如下:
| $ E_{\mathrm{ads}}=E_{\mathrm{O}_{3}+\mathrm{CuO}}-E_{\mathrm{CuO}}-E_{\mathrm{O}_{3}} $ | (1) |
其中, EO3 + CuO是吸附O3后系统的总能, ECuO和EO3分别为CuO和气相O3分子的总能.
| $ \Delta \rho=\rho_{\mathrm{O}_{3}+\mathrm{CuO}}-\rho_{\mathrm{CuO}}-\rho_{\mathrm{O}_{3}} $ | (2) |
其中, ρO3 + CuO是吸附O3后系统的电荷密度; ρCuO和ρO3分别为CuO和O3分子的电荷密度.
2 结果和讨论 2.1 O3在洁净表面的吸附O3在CuO（111）的表面吸附形成4种比较稳定的O3构型, 如图 1（a）所示. 表 1给出了O3在CuO（111）表面的吸附能. 当O3的两端的O原子分别与CuO表面的两个Cu原子成键形成桥式结构时, 吸附能较大, 构型bri（1）和bri（2）的吸附能分别为-0.77和-1.22 eV. 当O3一端的O原子为自由端时, 吸附能较小, 构型top（1）和top（2）的吸附能分别为-0.75和-0.28 eV. bri（1）中O3的O－O键分别为0.13和0.14 nm, top（1）中的O－O键分别为0.13和0.14 nm, 此时O3趋于分解成O*和O2*; 而bri（2）中O3的O－O键由0.13对称地拉长至0.14 nm, 说明O3的O和表面的Cu具有较强的相互作用, 形成的Cu－O键长较短（均为0.20 nm）, 因此该构型的吸附能最大.
|
图 1 O3吸附在CuO（111）表面的结构侧视图 Fig.1 Side view of the configurations of O3 adsorbed on CuO(111) surface |
|
表 1 O3在CuO（111）表面的吸附能 Table 1 Adsorption energy of O3 on CuO(111) surface |
为进一步分析成键情况, 计算了构型bri（2）的差分电荷密度（CDD）图. 如图 2（a）所示, Cu吸附位的电荷密度明显减小, O3末端两个O-Cu原子下方的电荷密度显著增强, 说明形成了Cu－O离子键. O3中间的O-O附近电荷密度也有所增强, 说明电荷从CuO表面转移到O3. O3的3个O原子之间的区域电荷密度减小, 表示O3中的O－O键减弱, 与O－O键长增长相符. 图 2（b）所示为构型bri（2）中O3末端和O3中间O原子的投影态密度（PDOS）图, s轨道和p轨道分别用实线和虚线表示. O-Cu和O-O的s轨道均在-14 eV处有强峰, 在-7 eV处有弱峰; p轨道在-8.56、-8.31、-6.82、-6.57、-0.58和-0.09 eV处交叠, 说明O-Cu和O-O之间仍有很强的共价键. 由于Cu－O键的形成, 在-2.70、-1.83、-1.59和-1.34 eV处有强峰, 而O-O在该处无峰, 说明O-Cu和Cu形成了很强的离子键.
|
图 2 （a）构型bri（2）和ovbri（1）的CDD图; （b）构型bri（2）和（c）ovbri（1）的投影态密度（黄色和蓝色部分分别表示电荷密度增加和减少了50 e/（nm）3; E=0 eV为费米能级; O-Cu表示O3中与Cu成键的O, O-O表示O3中间的O） Fig.2 (a) CDD plotof bri (2) and ovbri (1) with an isosurface level of 50 e/(nm)3; charge accumulation is in yellow and charge depletion is in blue. PDOS of O-Cu and O-O in (b) bri (2) and (c) ovbri (1), where O-Cu represents the O at one edge of O3 that bound with Cu and O-O represents the O in the middle of O3. The energy at E= 0 eV represents the Fermi energy. |
氧空位是金属氧化物表面普遍存在的一种表面缺陷, 在氧化还原反应中起重要作用, 因此, 研究了O3在具有氧空位的CuO表面的吸附. 图 1（b）是O3在含有氧空位的CuO（111）表面Cu位上的3种吸附构型ovbri（1）、ovbri（2）和ovbri（3）, 均为桥式构型, 不再形成O3一端的O原子为自由端的构型. O3在含有氧空位的CuO（111）表面的吸附能远高于在完整CuO（111）表面的吸附能, Cu－O键长变短, 说明Cu和O-Cu的相互作用变强. 其中ovbri（3）吸附结构的吸附能最大, 为-2.95 eV, O3一端的O原子填补了氧空位, 另一端的O原子与表面两个Cu形成桥键, O－O键分别拉长至0.14和0.15 nm.
由图 2（a）构型ovbri（1）的差分电荷密度图可以看出, 与完整CuO表面的O3吸附相比, O3在具有氧空位的CuO表面吸附时, 电荷转移明显增强, Cu附近的电荷密度减弱程度和O3的O原子附近电荷密度增强程度均大于完整表面, 说明O3分子和CuO表面存在强烈的相互作用. 从图 2（c）构型ovbri（1）的PDOS图可以看到, O-Cu和O-O在费米能级附近的态密度减弱并向低能量偏移, O-Cu和O-O的p轨道在-8.63、-6.90、-0.52和-0.22 eV处交叠, 能量低于完整表面, 说明O3在具有氧空位的表面吸附的结构更稳定.
2.3 O3在CuO表面的分解路径和表面氧物种的生成反应路径终态吸附的O2*和O*是能够引发氧化反应的表面活性氧, 表 2列出了这些表面活性氧的价态和键长. 吸附的O*价态为-0.55和-1.61; 完整表面吸附的O2*价态分别为-0.04、-0.21、-0.18和-0.15, O－O键长分别为0.12、0.13、0.13和0.13 nm, 略高于自由氧分子的O－O键长（0.12 nm）, 氧空位处吸附的O2*价态为-0.10和-1.03, O－O键长为0.12和0.15 nm, 可见O3分解后形成的表面吸附O2*和O*外层电子配位均不饱和, 形成了超氧物种[33], 相对于自由O2更容易参与氧化反应.
|
|
表 2 反应路径终态表面活性氧的价态和键长 Table 2 Valence states and bondlengths of surface active oxygen in the final state of reaction path |
根据O3在CuO催化剂上吸附的初始结构, 提出构型bri（1）、bri（2）、ovbri（2）和ovbri（3）分解并形成表面活性氧的合理途径. O3在CuO（111）表面分解的活化能垒图, 过渡态和终态的几何结构如图 3所示. 对于O3在完整CuO（111）表面的分解, 反应从构型bri（1）和bri（2）开始, O3两端的O原子与CuO表面的两个Cu原子形成桥式结构, 随着O3的O原子与CuO表面的Cu原子相互作用, 电子从Cu原子转移到O原子上, O3与CuO结合形成中间体, 即过渡态TS1和TS2, 最终O3解离成O*和O2*, 反应的最高能垒分别为0.58（路径bir（1））和0.52 eV（路径bir（2））. O2*是弱吸附分子, 从表面脱附分别需要0.42（路径bir（1））和0.52 eV（路径bir（2））的能量, 脱附后形成自由氧分子. O*和CuO表面的晶格氧结合, 形成另一个O2*, 此O2*脱附需要较高的能量2.06 eV, 脱附后形成含有氧空位的表面. 该过程中, 氧空位形成的能垒最高; 若不考虑氧空位的形成, 以bir（2）为初始结构的反应路径能垒较低, 为0.52 eV, 是反应可能发生的路径.
|
图 3 O3在CuO表面的分解路径和活化能垒图 Fig.3 O3 decomposition path and activation energy barrier on CuO surface |
O3在具有氧空位表面的分解从构型ovbri（1）和ovbir（3）开始, 在ovbir（1）中, O3中的两个O与表面的两个Cu分别成键, 形成较为对称的桥式结构, 吸附能为-1.66 eV. 在ovbir（3）中, O3一端的O原子与CuO表面的两个Cu原子形成桥式结构, 另一端的O原子填补了氧空位, 吸附能为2.95 eV. 吸附后的O3经TSl和TS2两个过渡态形成O2*, 另一个O变成CuO表面的晶格氧, 反应能垒分别为0.30（路径ovbir（1））和0.12 eV（路径ovbir（3））, O2*从反应表面脱附后, 留下完整的CuO表面, 脱附所需能量分别为0.27（路径ovbir（1））和0.52 eV（路径ovbir（3））. 由此可见, 以ovbir（3）为初始结构的反应路径能垒较低, 是反应可能发生的路径, 但是该路径最终O2*脱附的能垒比ovbir（3）路径高. 同时, 在具有氧空位表面, O3分解的反应能垒低于完整表面. 氧空位的形成提高了吸附能, 使O3分子更容易吸附在CuO表面, 从而促使催化反应发生; 同时, 氧空位的形成降低了反应能垒, 利于反应进行, 因此在具有氧空位的表面O3更易发生分解.
3 结论采用GGA+U计算详细探讨了O3在CuO（111）表面的吸附和解离过程. O3分子在完整CuO表面吸附形成4种稳定构型, 3种构型为化学吸附, 吸附能最高为-1.22 eV, 1种构型为物理吸附, 吸附能为-0.22 eV. O3在具有氧空位的CuO表面吸附形成3种稳定构型, 均为化学吸附, 吸附能最高为-2.95 eV, 显著高于完整表面的吸附能. O3吸附后, Cu吸附位的电荷密度减小, O3中的O原子附近的电荷密度显著增强, 电荷从CuO表面转移到O3, 并形成Cu－O离子键. O3在完整CuO（111）表面的分解从桥式构型bri（1）和bri（2）开始, 通过反应能垒分别为0.58和0.52 eV的过渡态, 解离成O*和O2*, O2*的脱附能分别为0.42和0.52 eV. O*和CuO表面的晶格氧结合形成的O2*脱附后形成氧空位, 脱附能为2.06 eV. O3在具有氧空位表面的反应从构型ovbri（1）和ovbir（3）开始, 经反应能垒分别为0.30和0.12 eV的过渡态解离形成O*和O2*, O*变成CuO表面的晶格氧, O2*从CuO表面脱附后形成完整CuO表面, 脱附能分别为0.27和0.52 eV. 氧空位的形成增大了吸附能, 降低了反应能垒, 使O3分子更容易吸附在CuO表面, 加快了O3的催化分解. O3分解后形成了超氧物种, 提高了表面的氧化活性.
| [1] |
Tao L, Zhang Z, Chen P, et al. Thin-felt Al-fiber-structured Pd-Co-MnOx/Al2O3 catalyst with high mois-ture resistance for high-throughput O3 decomposition[J]. Appl Surf Sci, 2019, 481: 802–810.
DOI:10.1016/j.apsusc.2019.03.134 |
| [2] |
Zhang Yi(张溢), Shi Jin(石琎), Zhang Zhi-xiang(张志翔), et al. Study on simultaneous removal of formald-ehyde and ozone by manganese oxygen octahedral molecular sieve at room temperature(锰氧八面体分子筛室温同时去除甲醛和臭氧的性能研究)[J]. J Mol Catal(China)(分子催化), 2019, 33(3): 201–208.
|
| [3] |
Tao L, Zhao G, Chen P, et al. High‐performance Co‐MnOx composite oxide catalyst structured onto Alfiberfelt for high‐throughput O3 decomposition[J]. Chem-CatChem, 2019, 11(3): 1131–1142.
|
| [4] |
Lin F, Wang Z, Zhang Z, et al. Flue gas treatment with ozone oxidation: An overview on NOx, organic pollutants, and mercury[J]. Chem Eng J, 2020, 382: 123030.
DOI:10.1016/j.cej.2019.123030 |
| [5] |
Ding Jing-ya(丁静亚), Xi Yang-chun(席阳春), Zhou Yi(周毅), et al. Efficient decomposition of ozone over Co-modified g-C3N4 catalysts(g-C3N4负载Co及其分解臭氧性能研究)[J]. J Mol Catal(China)(分子催化), 2020, 34(1): 19–27.
|
| [6] |
Kamal M S, Razzak S A, Hossain M M. Catalytic oxidation of volatile organic compounds (VOCs)-A review[J]. Atmos Environ, 2016, 140: 117–134.
DOI:10.1016/j.atmosenv.2016.05.031 |
| [7] |
Chen J, Chen X, Xu W, et al. Homogeneous introduction of CeOy into MnOx-based catalyst for oxidation of aromatic VOCs[J]. Appl Catal B-Environ, 2018, 224: 825–835.
DOI:10.1016/j.apcatb.2017.11.036 |
| [8] |
a. Whitehead J C. Plasma-catalysis: The known knowns, the known unknowns and the unknown unknowns[J]. J Phys D Appl Phys, 2016, 49(24): 243001. b. Ji Zeng-bao(季增宝), Yang Bao-jun(杨宝军), Sun Jian-yu(孙剑宇), et al. Treatment UDMH waste water byO3/UV heterogeneous catalytic method(O3/UV多相催化氧化法处理UDMH废水的研究)[J]. J Mol Catal(China)(分子催化), 2022, 36(2): 99-108. |
| [9] |
Xu X, Wang P, Xu W, et al. Plasma-catalysis of metal loaded SBA-15 for toluene removal: Comparison of continuously introduced and adsorption-discharge plasma system[J]. Chem Eng J, 2016, 283: 276–284.
DOI:10.1016/j.cej.2015.07.050 |
| [10] |
Zhang Zhen-zhen(张珍珍), Li Xin-he(李鑫恒). Recent advances on catalytic oxidation of formaldehyde(基于催化氧化技术去除甲醛的研究进展)[J]. J Mol Catal(China)(分子催化), 2019, 33(4): 382–390.
|
| [11] |
Dong Y, Zhao J, Zhang J Y, et al. Synergy of Mn and Ni enhanced catalytic performance for toluene combustion over Ni-doped α-MnO2 catalysts[J]. Chem Eng J, 2020, 388: 124244.
DOI:10.1016/j.cej.2020.124244 |
| [12] |
Einaga H, Maeda N, Nagai Y. Comparison of catalytic properties of supported metal oxides for benzene oxidation using ozone[J]. Catal Sci Technol, 2015, 5(6): 3147–3158.
DOI:10.1039/C5CY00315F |
| [13] |
Rashidi R, Yousefinejad S, Mokarami H. Catalytic ozon-ation process using CuO/clinoptilolite zeolite for the removal of formaldehyde from the air stream[J]. Int J Environ Sci Technol, 2019, 16(11): 6629–6636.
DOI:10.1007/s13762-018-2059-2 |
| [14] |
Turkay O, Inan H, Dimoglo A. Experimental and theor-etical investigations of CuO-catalyzed ozonation of humic acid[J]. Sep Purif Technol, 2014, 134: 110–116.
DOI:10.1016/j.seppur.2014.07.040 |
| [15] |
Gong S, Wu X, Zhang J, et al. Facile solution synthesis of Cu2O-CuO-Cu(OH)2 hierarchical nanostructures for effective catalytic ozone decomposition[J]. Cryst-EngComm, 2018, 20(22): 3096–3104.
|
| [16] |
Zhang X, Li X, Qin W. Investigation of the catalytic activity for ozonation on the surface of NiO nanoparticles[J]. Chem Phys Lett, 2009, 479(4/6): 310–315.
|
| [17] |
Wu J, Su T, Jiang Y, et al. In situ DRIFTS study of O3 adsorption on CaO, γ-Al2O3, CuO, α-Fe2O3 and ZnOat room temperature for the catalytic ozonation of cinnam-aldehyde[J]. Appl Surf Sci, 2017, 412: 290–305.
DOI:10.1016/j.apsusc.2017.03.237 |
| [18] |
Hu Xiao-hong(呼小红). Theoretical mechanism investi-gation on catalytic ozonation by metal oxide(金属氧化物催化臭氧生成羟基自由基机理理论研究)[D]. Harbin(哈尔滨): Harbin Institute of Technology(哈尔滨工业大学), 2013.
|
| [19] |
Ouali H, Lambert M C, Raymond L, et al. Mechanism of O3 sensing on Cu2O (111) surface: First principle calculations[J]. Appl Surf Sci, 2015, 351: 840–845.
DOI:10.1016/j.apsusc.2015.06.017 |
| [20] |
Liu Y, Yang W, Zhang P, et al. Nitric acid-treated birnessite-type MnO2 : An efficient and hydrophobic material for humid ozone decomposition[J]. Appl Surf Sci, 2018, 442: 640–649.
DOI:10.1016/j.apsusc.2018.02.204 |
| [21] |
Oison V, Ouali H, Lambert M C, et al. Experimental and ab initio study of the O3 detection at the CuO (111) surface[J]. Surf Sci, 2014, 622: 44–50.
DOI:10.1016/j.susc.2013.12.003 |
| [22] |
Kresse G, Furthmüller J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set[J]. Phys Rev B, 1996, 54(16): 11169.
DOI:10.1103/PhysRevB.54.11169 |
| [23] |
Kresse G, Furthmüller J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set[J]. Comp Mater Sci, 1996, 6(1): 15–50.
DOI:10.1016/0927-0256(96)00008-0 |
| [24] |
Perdew J P, Burke K, Ernzerhof M. Generalized gradient approximation made simple[J]. Phys Rev Lett, 1996, 77(18): 3865.
DOI:10.1103/PhysRevLett.77.3865 |
| [25] |
Blöchl P E. Projector augmented-wave method[J]. Phys Rev B, 1994, 50(24): 17953.
DOI:10.1103/PhysRevB.50.17953 |
| [26] |
Moreno J L V, Padama A A B, Kasai H. A density func-tional theory-based study on the dissociation of NO on a CuO (110) surface[J]. CrystEngComm, 2014, 16(11): 2260–2265.
DOI:10.1039/c3ce42103a |
| [27] |
Maimaiti Y, Nolan M, Elliott S D. Reduction mechanisms of the CuO (111) surface through surface oxygen vacancy formation and hydrogen adsorption[J]. Phys Chem Chem Phys, 2014, 16(7): 3036–3046.
DOI:10.1039/c3cp53991a |
| [28] |
Zheng C, Zhao H. The microscopic oxidation mechanism of NH3 on CuO (111): A first-principles study[J]. Fuel Process Technol, 2021, 213: 106712.
DOI:10.1016/j.fuproc.2020.106712 |
| [29] |
Liu L, Gong P, Shao G, et al. Density functional theory plus U study of methanol adsorption and decomposition on CuO surfaces with oxygen vacancy[J]. Aerosol Air Qual Res, 2022, 22(1): 210253.
DOI:10.4209/aaqr.210253 |
| [30] |
Zheng C, Cao J, Zhang Y, et al. Insight into the oxidation mechanism of a Cu-based oxygen carrier (Cu→ Cu2O → CuO) in chemical looping combustion[J]. Energ Fuel, 2020, 34(7): 8718–8725.
DOI:10.1021/acs.energyfuels.0c00941 |
| [31] |
Liu L, Shao G, Gong P, et al. Density functional theory studies on ortho-position adsorption of SO3 at step sites of a CaO surface with SO2 and CO2[J]. Fuel, 2022, 310: 122174.
DOI:10.1016/j.fuel.2021.122174 |
| [32] |
Henkelman G, Uberuaga B P, Jónsson H. A climbing image nudged elastic band method for finding saddle points and minimum energy paths[J]. J Chem Phys, 2000, 113(22): 9901–9904.
DOI:10.1063/1.1329672 |
| [33] |
Cen W, Liu Y, Wu Z, et al. A theoretic insight into the catalytic activity promotion of CeO2 surfaces by Mn doping[J]. Phys Chem Chem Phys, 2012, 14(16): 5769–5777.
DOI:10.1039/c2cp00061j |