【計算論文精讀】Appl. Surf. Sci.:DFT研究小分子在表面的吸附和活化規律 2023年12月8日 上午9:29 ? 計算 ? 閱讀 21 研究背景 氮氧化物會引起一系列的環境問題,如酸雨、煙霧和光化學煙霧。氨選擇性催化還原(NH3-SCR)因其相對成熟、高效的技術,已成為控制氮氧化物排放的主流技術。催化劑設計是NH3-SCR反應更為高效的重要舉措。因此,尋找新型、高效、環保的低溫NH3-SCR催化劑具有重要的現實意義。實驗證明了Fe3O4(111)表面對NH3-SCR反應的催化作用,但仍缺乏系統的原子水平研究。 武漢理工大學謝峻林和方德等人采用第一性原理計算方法研究了NH3、NOx和O2與Fe3O4(111)表面的相互作用,這一研究的目的是為了更好地理解NH3-SCR反應中小分子在Fe3O4(111)表面上的吸附和活化規律。 計算方法 本文采用Materials Studio軟件包的CASTEP模塊對其催化性能及電子性質進行理論研究,計算基于廣義梯度近似(GGA)下的PBE泛函計算電子交換相關性,為了避免層間的相互作用,將Z方向的真空空間設置為15 ?。截斷能設置為400 eV,將能量、最大力和最大位移的收斂公差分別設為1× 10-5 eV/atom、0.03 eV·?-1和0.001 ?。 結果與討論 Fe3O4(111)表面是自然生長的主要晶面,其催化活性高于其他表面。根據終端表面原子類型的不同,沿Fe3O4(111)方向有兩種主要的表面類型,可以表示為Fetet 端面和Feoct -tet端面。如圖1所示,Fetet表面由四面體Fe位點(Fetet)、孤立O位點(Oiso)、非孤立O位點(Onon-iso)和中空位點(Ohollow)組成。Feoct-tet端面包括最外八面體Fe位點(Feoct)、次外Fetet和Oiso。對于這兩種表面模型,最上面的三層原子被允許松弛,而其他原子層被固定,以確保它們的表面性質可以更準確地模擬。 圖1. Fe3O4(111)的Fetet端面和Feoct-tet端面結構圖 圖2 (a)和(b)顯示了氨吸附在表面的優化平衡構型。氨位于Fetet的上方。可以看出,當氨接近Fetet的頂部時,吸附能為-1.48 eV,是一種穩定的化學吸附能。此時,氨的N原子與Fetet之間形成了一個化學鍵,N-Fe鍵長為2.109 ?。Mulliken電荷可以解釋分子中電子的分布和轉移。電荷表示分子中每個原子以不同構型所擁有的電荷量,而電子的變化表示分子總電荷的增益或損失。可以觀察到,在構型(a)中,一小部分電子(0.04 e)從氨轉移到催化劑端面,這說明氨和Fetet之間存在電子相互作用。 Ohollow與H原子之間存在三方弱相互作用,使NH3以平衡構型(b)所示的形式吸附在表面,吸附能為-0.12 eV。它接近物理吸附,表明平衡構型不太穩定。然而,在三個Oiso的影響下,所有的N-H鍵長都被拉長。與其他平衡構型相比,NH3可以從表面獲得更多的電子(-0.10 e),兩者之間的電子相互作用更強,但這種平衡構型的吸附效果最差。 圖2. Fe3O4(111)的Fetet -端面吸附NH3的構型 圖3 (A)-(C)顯示了氨吸附在Feoct-tet端接表面的平衡構型。圖3 (A)和(B)中氨垂直吸附在Feoct和Fetet上的吸附能分別為?0.70 eV和?0.27 eV,而氨中的鍵長和鍵角只有極小的變化。這表明,當氨垂直吸附在氨端端表面時,它們之間只有微弱的相互作用。此外,N-Fe鍵長分別為2.125 ?和2.214 ?。這可能是由于亞外層吸附位點的空間位阻大于最外層,導致吸附能更小,N-Fe鍵長越長。這是最有利的吸附方法,即在氨、Feoct和Oiso之間生成N-Fe化學鍵和O-H氫鍵。 圖3. Fe3O4(111)的Feoct-tet端面吸附NH3的結構 在Fetet和Feoct-tet表面吸附NH3前后,相關原子或分子最穩定構型的PDOS如圖4所示。這些結果表明,兩種表面都對NH3在吸附位點上有穩定作用。進一步分析和比較兩個表面吸附位點的PDOS,可以更深入地了解NH3與表面的相互作用。這些共振峰表明NH3的N原子與Fetet(或Feoct)之間形成了N- Fe鍵。共振峰的形狀、數量和強度表明,NH3在Fetet端部表面更為穩定。N 2p和Fe 3d原子軌道的PDOS峰的大量重疊,揭示了N – Fe鍵存在雜化。此外,NH3非成鍵軌道的峰強度減弱,峰形狀變寬,證明主要是NH3孤對電子與表面有強相互作用。 圖4. 吸附NH3的態密度變化 隨后研究NO在Fe3O4(111)的Fetet端面和Feoct-tet端面的吸附位置。Fetet端面吸附時,NO在表面不同位置的吸附能在0.20 eV到1.32 eV之間,說明在Fetet端部表面既存在物理吸附又存在化學吸附。當NO吸附在表面的Ohollo上時,吸附能均為-0.20 eV,為強物理吸附。圖5中NO以N原子為頭吸附在Fetet上的平衡構型(a)是最穩定的。其吸附能最高,為-1.32 eV。有0.25 e的電子從表面轉移到NO。對于Feoct-tet端面,N原子對Feoct和Fetet)的吸附能分別為-2.23 eV和-0.16 eV。NO與吸附位點之間的距離分別為1.757 ?和1.798 ?。相反,對于圖6中的平衡構型(C)和(D),當NO以O原子為頭在Feoct和Fetet上時,吸附能分別只有-0.39 eV和-0.05 eV,相對較低。在上述四種平衡構型中,無論哪種形式的NO靠近Feoct-tet端表面,Feoct上的吸附能都大于Fetet上的吸附能,這主要是空間位阻效應造成的。 圖5. Fe3O4(111)的Fetet 端面吸附NO的構型 圖6. Fe3O4(111)的Feoct-tet端面吸附NO的結構 對于最佳平衡構型的PDOS進行分析,見圖7。與吸附前相比,主要貢獻的分子軌道對應的峰值位置在能量較高的方向上略有偏移。這些結果都表明NO在Fetet上被激活,其能量也有所增加。對比圖7 (B)和(C),可以發現8個主共振峰。此外,主要受Fetet的3d和4p原子軌道的影響,費米能級附近的NO π*反鍵軌道和費米能級以上的分子軌道的能量被拉高并完全分散在Fetet原子軌道中,這證明Fetet對NO的穩定作用。NO被Feoct端面激活,同時也受到該位點相對較強的穩定性的影響。更重要的是,Feoct很好地平衡了兩者的作用,使吸附能更高。 圖7. 吸附NO的態密度變化 隨后研究O2在其表面的吸附平衡構型。如圖8 (a)-(d)所示,Fetet端部表面的四個不同位置對O2的吸附。結構優化后,Oiso、Onon-iso和Ohollow均能吸附O2,吸附能約為-1.20 eV。O2在Fetet上形成O-Fe化學鍵,進一步加強了吸附效果,提高吸附能。平衡構型(a)中的O2受表面feet的影響,通過吸附位點從表面獲得0.18 e電子,O2的O-O鍵長被拉長,說明兩者之間存在較強的相互作用。其他三種平衡構型幾乎沒有變化,說明Fetet端部表面可能存在許多表面吸附氧。 對于Feoct-tet端部表面,分別研究了O2接近Feoct、Fetet和Oiso時的構型。O2水平吸附在Fetet上,兩個O原子同時與Fetet結合,形成三角形吸附氧結構,吸附能最高(-2.54 eV)。Mulliken電荷表明,分別有0.32 e和0.33 e電子從表面轉移到兩個O原子上,O2共獲得0.65 e電子。O2與Feoct相互作用強烈,形成兩個O-Fe化學鍵,鍵長分別為2.014 ?和2.020 ?。O-O鍵長被拉長至1.381 ?,表明O2被激活。 圖8. Fe3O4(111)的Fetet 端面吸附O2的構型 圖9. Fe3O4(111)的Feoct-tet端面吸附O2的結構 吸附氧氣前后端PDOS見圖10。O2的價層電子數比NO多1個,分子軌道有一定相似性。DOS圖中?25.91 eV和?13.86 eV的峰值可歸因于O2的σ 2s和σ*2s軌道。當O2吸附在Fetet上時,O2所有已占軌道對應的峰都發生了分裂,能量略有降低。結果表明,該位點對O2既有穩定作用,又有活化作用。對于Feoct表面,也可以看到,所有已占軌道對應的峰都被分裂,這表明O2的分子軌道被輕微破壞。而成鍵軌道和反鍵軌道分別朝著能量增加和能量減少的方向運動,趨向于O原子本身的原子軌道能量,這表明O2進一步被激活。對比發現,Feoct端部表面上穩定的三角吸附氧結構是非常重要的。 圖10. 吸附O2的態密度變化 NH3– SCR反應中NH3吸附在催化劑表面活性位點上形成-NH2的脫氫過程稱為NH3活化。這是一個重要的基本反應,使用CASTEP中的LST/QST算法在Fe3O4(111)表面上搜索NH3氧化脫氫的過渡態(TS),結果如圖11所示。NH3首先吸附在Fe3O4(111)表面,獲得最穩定的構型,分別釋放1.48 eV和1.19 eV的能量。在Fetet表面,在最靠近吸附位點的晶格氧的影響下,NH3發生氧化脫氫反應。其中一個氫與晶格中的氧結合在表面形成羥基,然后自身轉化為-NH2。該反應的能壘較高(1.49 eV),與吸附能相當。可能的原因是NH3在Fetet表面的穩定吸附構型是垂直的,晶格氧對遠處NH3 中H的影響有限。 相反,對于Feoct-tet端部表面,由于穩定的水平吸附構型,NH3的氫與最近的晶格氧之間存在氫鍵相互作用,使N-H鍵在一定程度上被削弱和激活。NH3氧化脫氫反應能壘(0.79 eV)降低。NH3氧化脫氫的整個反應過程是吸熱的,但反應熱足以被吸附時釋放的能量所覆蓋,說明在Fe3O4(111)表面發生的基元反應是一個自發的過程。 圖11. NH3活化步驟的過渡態結果 結論與展望 本文探討了Fe3O4(111)兩個暴露表面上NH3、NO和O2分子的吸附規律,由于水平吸附構型,NH3明顯受到最近晶格氧的影響,這使得氧化脫氫反應更容易發生。這為在原子水平上理解其表面上的NH3– SCR反應提供依據。 文獻信息 Hou, S., Fang, D., Jin, Q., Ye, Y., Li, Q., He, F., & Xie, J. (2022). First principle study on the interactions of NH3, NOx and O2 with Fe3O4 (1 1 1) surfaces.?Applied Surface Science,?605, 154645. https://doi.org/10.1016/j.apsusc.2022.154645 原創文章,作者:計算搬磚工程師,如若轉載,請注明來源華算科技,注明出處:http://www.zzhhcy.com/index.php/2023/12/08/5c2a1a72cb/ 贊 (0) 0 生成海報 相關推薦 崔小強等:應用機器學習快速篩選過渡金屬單原子負載N摻雜γ相石墨炔用于一氧化碳電還原C1產物的理論研究 2024年1月9日 Nature 子刊 | 化學家和機器人都可以讀懂,用于機器人合成可重復性的通用化學編程語言 2024年1月18日 Angew.:調控金屬間的電荷轉移促進水氧化電催化 2023年11月23日 復旦大學董曉麗研究員團隊綜述:路在何方之未來低溫電池電解液 2024年4月10日 各種尺寸、形狀都適用,圖卷積神經網絡探索金屬納米粒子的電化學穩定性 2024年3月13日 【DFT+實驗】山東大學Angew.:近100%選擇性!BP/BWO S-型異質結助力苯胺氧化和CO2光轉化 2024年2月28日