陳宏善陳華君

H2在MgO團簇吸附的從頭計算研究*

陳宏善陳華君

(西北師范大學物理與電子工程學院，蘭州730070)
(2010年12月22日收到;2011年1月27日收到修改稿)

根據(jù)最穩(wěn)定幻數(shù)MgO團簇的結構特點，對H2在巖鹽和管狀結構(MgO)9，12表面的吸附性質(zhì)進行了從頭計算研究．結果表明H2可以在處于團簇不同位置的Mg正離子或氧負離子上發(fā)生物理吸附;在Mg離子上H2以側(cè)位方式吸附并向Mg轉(zhuǎn)移電子，在O離子上H2以端位方式吸附并被明顯極化．吸附的穩(wěn)定性主要依賴于吸附離子的位置即配位數(shù)，Mg/O離子的配位數(shù)越低則吸附越強;在配位數(shù)相同時，H2在Mg離子上的吸附強于在O離子上的吸附．隨吸附位置的不同，吸附能在0.03—0.08 eV．計算結果與最近H2在MgO塊體表面吸附的實驗結果一致．

團簇，MgO，H2吸附，DFT

PACS:36．40.－c，31．15．A－，68．43．Bc

1.引言

能源在現(xiàn)代社會發(fā)展中具有基礎地位，不斷增長的能源需求與化石燃料的有限儲量也是當今社會的主要矛盾之一．氫的儲量豐富、燃燒值高、清潔無污染，是化石燃料的最佳替代品．但氣、液、固態(tài)儲氫的安全性及化學儲氫的經(jīng)濟性阻礙了氫能的實際應用，因此對各種材料儲氫性能的研究具有重要意義［1—6］．團簇具有豐富的物理化學性質(zhì)和極大的表面/體積比，團簇儲氫有可能解決氫的儲運問題．H2在離子型固體表面的吸附很早就引起了人們的注意，對于H2在晶體氧化鎂表面的吸附，曾利用氦原子散射或中子散射進行了研究［7—9］．最近實驗表明H2可以在氧化鎂(100)表面發(fā)生物理吸附，H2以平行于表面的方式吸附在Mg離子之上，吸附率在0.2—0.8之間時吸附位呈c(2×2)有序結構，而吸附率在0.8—1.05之間時呈六邊形結構［9］．理論研究方面，Hermansson曾利用分子動力學結合從頭計算方法對H2在MgO(111)面上的吸附進行了模擬［10］．最近，Zhang等人用密度泛函方法對H2在MgO(001)表面的吸附進行了計算［11］，Dawoud還利用Monte Carlo方法和微擾理論對D2在(001)面上的吸附進行了研究［12］．對于H2在MgO團簇上的吸附，僅有Morokuma等人利用(MgO)n=4，6，8模擬氧化鎂晶體的面、棱、角，用從頭計算方法研究了H2在其上的解離問題［13，14］，而他們的計算并沒有考慮MgO團簇本身的穩(wěn)定結構特點．實驗［15—17］和理論［16—23］研究表明，(MgO)n在n為3的倍數(shù)時特別穩(wěn)定，而這些幻數(shù)團簇的最穩(wěn)定結構是巖鹽結構或六元管狀結構．本文利用從頭計算方法研究了H2在巖鹽結構(MgO)9，12和管狀(MgO)12的吸附規(guī)律．

2 ．計算方法

本文計算是在密度泛函理論B3 LYP［24］框架下進行的，即電子交換能采用HF和Becke三參數(shù)混合形式，電子相關能的定域部分使用函數(shù)VWNⅢ［25］，而非定域部分使用LYP的泛函形式［26］．計算使用雙ζ基組6-31 G．計算過程首先對(MgO)n團簇的結構進行優(yōu)化，對各種吸附方式進行非限制完全優(yōu)化的基礎上，進一步進行頻率分析并消去可能存在的虛頻．為了比較，還利用二級微擾理論MP2/6-311G對H2在(MgO)9上的吸附進行了計算．所有計算使用Gaussian03程序完成［27］．

3 ．結果討論

3．1.H2在巖鹽結構(MgO)9表面的吸附

考慮團簇的不同面、邊、角，并考慮H2吸附在單個原子之上、不同原子之間以及H2的不同取向，設計了51種不同的吸附模式，對各結構進行完全優(yōu)化及頻率計算，對存在的虛頻進行簡正模式分析并消去虛頻以得到穩(wěn)定結構．結果表明，比較穩(wěn)定的吸附方式是H2吸附在單個Mg或O離子之上．在Mg離子上方，H2以側(cè)位方式吸附，即分子的方向與分子到Mg離子的連線垂直，兩個H原子處于近似等價的位置;而在O離子上方，H2以端位方式吸附，即分子的方向與到O離子的連線平行．圖1中給出了H2在(MgO)9上最穩(wěn)定的6個吸附方式．表1中列出了H2分子的吸附距離、吸附能以及兩個H原子上的Mulligan電荷．H2吸附在Mg離子上時，給出的是Mg到H2分子中心的距離，H2吸附在O離子上時，給出的是O到較近H原子的距離;吸附能被定義為吸附體系的能量與孤立團簇和H2分子的能量之差．H2分子在MgO團簇上的吸附僅為弱的物理吸附，吸附距離在2.4—3.6之間，吸附能主要與Mg/O離子的位置，即它們的配位數(shù)相關．當吸附在3配位的角上時，吸附能較大，吸附在5配位的面上時，吸附能最小;Mg，O離子配位數(shù)相同時，H2在Mg離子上的吸附比在O離子上的吸附更強．電荷布居數(shù)表明，H2吸附在Mg離子上時，H2分子向Mg離子轉(zhuǎn)移電子，H2吸附在O離子上時，H2分子被極化．表1還給出了在MP2/6-311 G水平計算得到的吸附能，表明兩種方法給出的結果非常接近．

3.2.H2在巖鹽結構(MgO)12表面的吸附

結合(MgO)12的結構特點及H2在(MgO)9上的吸附性質(zhì)，設計了H2在(MgO)12上吸附的三十余種不同方式，對各個結構進行完全優(yōu)化并消去虛頻得到穩(wěn)定的吸附結構．圖1中給出了最穩(wěn)定的9種吸附方式．同樣，吸附能的大小主要依賴于吸附位置，吸附離子的配位數(shù)越小吸附能越大．與(MgO)9上的各對應吸附方式比較，H2在(MgO)12上的吸附明顯增強．在(MgO)12上，也存在H2吸附在Mg/O離子對上方的穩(wěn)定結構，H2分子的方向大約平行于同O離子的連線而垂直于同Mg離子的連線，但這種吸附結構與H2吸附在Mg或O離子正上方(亦是穩(wěn)定的吸附結構)的穩(wěn)定性差別很?。?/p>

我們還優(yōu)化了多個H2分子在(MgO)12上的吸附結構，吸附方式的選擇考慮了角、邊、面上的不同離子并把H2分子限制在團簇的同一側(cè)．圖1中給出了3個這樣的吸附結構，優(yōu)化表明當多個H2分子同時吸附在團簇表面時，吸附結構與單個分子吸附時變化很小．當兩個H2吸附在Mg3c和O4c時，吸附能0.113 eV略小于兩個H2分子分別吸附在對應位置的吸附能之和(0.085+0.039=0.124)．當兩個H2同時吸附在O3c和Mg5c時，吸附能0.084 eV比兩個H2分子單獨吸附在對應位置的吸附能(0.069+ 0.034=0.103)小20%．而當三個H2分子吸附在同一面上的Mg3c，O4c和O5c/Mg5c時，吸附能0.119 eV比三個H2分別吸附的能量之和(0.154 eV)小23%．這表明由于H2在MgO團簇表面僅存在較弱的物理吸附，多個H2分子可以同時形成穩(wěn)定的吸附，吸附的強度較單個分子吸附時有所減弱，但減弱的程度不大．

3.3.H2在管狀(MgO)12表面的吸附

考慮管狀(MgO)12的對稱性及結構特點，設計了H2分子在管外各個不等價位置的吸附方式并進行結構優(yōu)化和頻率分析，圖1給出了最穩(wěn)定的4種吸附方式．與在巖鹽結構上的吸附類似，H2分子以側(cè)位方式吸附在Mg離子上方或以端位方式吸附在O離子上方，且吸附能與在巖鹽結構上具有相同的配位數(shù)的Mg/O離子上的吸附能接近．

對于H2在MgO表面的吸附，只有Hermansson等人的計算給出了解離吸附［10］．他們利用分子動力學對MgO層狀模型進行模擬退火，結合Hartree-Fock CO-LCAO計算得到H2在(111)面上形成化學吸附，吸附能為0.6 eV．Zhang等人用密度泛函方法對H2在MgO(001)表面的吸附進行了研究［11］，表明解離吸附的能壘為2.2 eV，因此只能發(fā)生物理吸附，而在Mg，O離子及橋位上形成物理吸附的能量分別為0.054，0.015和0.053 eV．Dawoud利用Monte Carlo方法結合微擾理論對D2在(001)面上的研究也表明僅發(fā)生物理吸附［12］，隨吸附量的增加，D2在(001)面上形成p(2×2)，p(4×2)等結構，而在Mg，O離子上的吸附能為0.016和0.012 eV．對于H2在(MgO)n=4，6，8上的吸附，Morokuma利用從頭計算研究表明解離過程的能壘在2.3 eV以上［13，14］，對H2的物理吸附，他們僅給出了在(MgO)4上的結果，不同精度的計算表明H2在Mg， O離子上的吸附能約為0.08 eV，在Mg離子上的吸附略強于在O離子上的吸附．

3．4.H2吸附在(MgO)n團簇上的分子軌道特征

對單個Mg，O原子，原子軌道的能級順序為Mg的2 s，2 p，O的2s，2 p，Mg的3s，H2分子的軌道能介于O的2p和Mg的3s之間．Mg，O原子形成(MgO)n團簇后，Mg的3 s電子向O的2 p空軌道轉(zhuǎn)移，O的軌道能級被抬高，它們形成的能帶位于約－14—－10 eV之間，而最低空軌道具有很小的正的軌道能，因此(MgO)n團簇的能帶具有典型的絕緣體特征．當H2分子吸附在MgO團簇上時，其軌道能降低到約－16 eV，位于O的2 p和Mg的3s形成的能帶之下;由于團簇的價帶是滿的，它們不能與H2有效的形成分子軌道，因此吸附為弱的物理吸附．圖2給出了H2分子吸附在(MgO)9上時有H參與形成的典型軌道．我們看到由H2為主形成的軌道保持了自由H2分子的軌道特征，同時H的軌道也以很小的分量參與了一些團簇的軌道．當H2吸附在Mg離子上時，由于Mg參與形成分子軌道的3 s是空的(即形成正離子)，H2分子以側(cè)位方位吸附在離子上方并向Mg離子轉(zhuǎn)移電子;當H2吸附在O離子上時，由于O的2 p軌道是占滿的并在H2的附近形成了較大的電荷分布，H2分子以端位方式吸附在O離子上方并被極化．

3.5.吸附H2的脫附能量曲線

最后我們考察了吸附在巖鹽結構(MgO)12上的H2分子的脫附過程．通過控制H2分子到吸附位(Mg/O離子)的距離進行部分結構優(yōu)化得出脫附過程的能量變化．圖3給出了吸附在Mg3c及O4c/Mg3c的H2在脫附過程的勢能曲線．正如一般的物理吸附所具有的普遍特征，隨著H2分子吸附距離的增大，體系的能量單調(diào)升高，并不存在能壘．

4 ．結論

設計了H2在巖鹽結構(MgO)9及巖鹽和管狀(MgO)12團簇上的各種可能吸附方式，對吸附性質(zhì)進行了第一性原理優(yōu)化計算．結果表明H2在MgO團簇上存在物理吸附，且能夠被吸附在不同位置的Mg正離子和O負離子之上．吸附的穩(wěn)定性主要決定于Mg/O離子的位置，即他們的配位數(shù)，當H2被吸附在3配位的邊角位置時吸附最強，吸附在4配位的棱上時吸附強度減弱，吸附在5配位的面上位置時強度最弱;當Mg/O離子所處位置相同時，H2在Mg離子上的吸附更加穩(wěn)定．當H2被吸附在Mg離子上時，H2的方向垂直于到Mg離子的方向，兩個H近似等價并向Mg離子轉(zhuǎn)移電荷，而在O離子上方H2以端位方式吸附并被明顯極化．H2在MgO團簇的吸附能介于0.03—0.08 eV，當在同一團簇吸附多個H2分子時，吸附強度會有所減弱．

［1］Schlapbach L，Züttel A 2001 Nature 414 353

［2］Coontz R，Hanson B 2004 Science 305 957

［3］Zhao Y F，Kim Y H，Dillon A C，Heben M J，Zhang S B 2005 Phys．Rev．Lett．94 155504

［4］Zhou J J，Chen Y G，Wu Z L，Zheng X，F(xiàn)ang Y C，Gao T 2009 Acta Phys．Sin．58 4853(in Chinese)［周晶晶、陳云貴、吳朝玲、鄭欣、房玉超、高濤2009物理學報58 4853］

［5］Ni M Y，Wang X L，Zeng Z 2009 Chin．Phys．B 18 357

［6］Liu X Y，Wang C Y，Tang Y J，Sun W G，Wu W D 2010 Chin．Phys．B 19 036103

［7］Degenhardt D，Lauter H J，Haensel R 1987 Jpn．J．Appl．Phys．26 341

［8］Skofronick J G，Toennies J P，Traeger F，Weiss H 2003 Phys．Rev．B 67 035413

［9］Larese J Z，Arnold T，F(xiàn)razier L，Hinde R J，Ramirez－Cuersta A J 2008 Phys．Rev．Lett．101 165302

［10］Hermansson K，Baudin M，Ensing B，Alfredsson M，Wojcik M 1998 J．Chem．Phys．109 7515

［11］Wu G X，Zhang J Y，Wu Y Q，Li Q，Chou K C，Bao X H 2009 J．Alloys Compd．480 788

［12］Dawoud J N，Jack D B 2009 Appl．Surf．Sci．256 1443

［13］Sawabe K，Koga N，Morokuma K，Lwasawa Y 1992 J．Chem． Phys．97 6871

［14］Sawabe K，Koga N，Morokuma K，Lwasawa Y 1994 J．Chem．Phys．101 4819

［15］Ziemann P J，Castleman A W Jr．1991 J．Chem．Phys．94 718

［16］Saunders W A，1989 Z．Phys．D:At．Mol．Clust．12 601

［17］Katakuse I，Ichihara T，Ito H，Hirai M 1990 Mass Spectr．4 16

［18］Fu X H，Yin Y，Zhang L，Ye H 2009 Acta Phys．Sin．58 5007 (in Chinese)［傅興海、尹伊、張磊、葉輝2009物理學報58 5007］

［19］Dong R B，Chen X S，Wang X F，Lu W 2008 J．Chem．Phys．129 044705

［20］Carrasco J，Illas F，Bromley S T 2007 Phys．Rev．Lett．99 235502

［21］Jain A，Kumar V，Sluiter M，Kawazoe Y 2006 Comput．Mater．Sci．36 171

［22］Bawa F，Panas I 2002 Phys．Chem．Chem．Phys．4 103

［23］Calvo F 2003 Phys．Rev．B 67 131403

［24］Becke A D 1993 J．Chem．Phys．98 5648

［25］Vosko S H，Wilk L，Nusair M 1980 Can．J．Phys．58 1200

［26］Lee C，Yang W，Parr R G 1988 Phys．Rev．B 37 785

［27］Frisch，M J，Trucks G W，Schlegel H B et al Gaussian 03．Revision E．01．Wallingford CT:Gaussian Inc．2004．

PACS:36．40．－c，31．15．A－，68．43．Bc

*Project supported by the National Natural Science Foundation of China(Grant No．20873102)and the Northwest Normal University(Grant No．NWNU-KJCXGC-03-62)．

E-mail:chenhs@nwnu．edu．cn

Adsorption of H2on MgO clusters studied by ab initio method*

Chen Hong-ShanChen Hua-Jun
(College of Physics and Electronic Engineering，Northwest Normal University，Lanzhou 730070，China)
(Received 22 December 2010;revised manuscript received 27 January 2011)

Based on the most stable structures of the magic MgO clusters，the adsorption of H2onto the rocksalt and hexagonal tube(MgO)9，12is studied using first principal method．The results show that physical adsorption can be formed on Mg+qor O－qions at different sites of the clusters．On the top of the Mg+qions，H2is adsorbed in the side-on manner and donates electrons to the ions．On the top of the O－qions，H2is adsorbed in the end-on manner and the molecule is polarized．The strength of the adsorption depends mainly on the location of the Mg/O ions，and the smaller the coordination number of the ions，the stronger the adsorption of H2is．For the Mg/O ions with the same coordination number，the adsorption formed onto the Mg+qis stronger．According to the Mg or O ions and their locations，the adsorption energy ranges from 0.03 eV to 0.08 eV．

cluster，MgO，H2adsorption，DFT

*國家自然科學基金(批準號:20873102)和西北師范大學科技創(chuàng)新工程(批準號:NWNU-KJCXGC03-62)資助的課題．

E-mail:chenhs@nwnu．edu．cn