第一性原理計算Mg2Ni氫化物的電子結(jié)構(gòu)及其穩(wěn)定性分析*

劉春華歐陽楚英嵇英華

第一性原理計算Mg2Ni氫化物的電子結(jié)構(gòu)及其穩(wěn)定性分析*

劉春華歐陽楚英 嵇英華

(江西師范大學(xué)物理系，江西省光電子重點(diǎn)實(shí)驗室，南昌330022)
(2010年7月8日收到;2010年10月9日收到修改稿)

對Mg2Ni及其氫化物的能量和電子結(jié)構(gòu)進(jìn)行了第一性原理計算，并對Mg2Ni低溫氫化物和高溫氫化物的穩(wěn)定性進(jìn)行了分析．結(jié)果發(fā)現(xiàn):在Mg2Ni中，Mg與Ni存在較強(qiáng)的相互作用．H原子的加入使得Mg的價電子向Ni的3d軌道轉(zhuǎn)移，并在Ni-3d軌道和Ni-4p軌道間形成了帶隙．Ni-4s和H-1 s電子之間的成鍵作用，使得Mg和Ni間相互作用減弱了;從LT-Mg2NiH4中去掉一個H原子比HT-Mg2NiH4中去掉一個H原子需要更多的能量，表明高溫相有利于脫氫反應(yīng)，而低溫相更適合儲氫反應(yīng)．

第一性原理，Mg2Ni，Mg2NiH4，結(jié)構(gòu)穩(wěn)定性

PACS:71．15．Mb，31．15．es，31．10．+z

1.引言

在大力提倡綠色能源的今天，氫作為綠色燃料被大規(guī)模應(yīng)用于各種燃料電池．發(fā)展氫經(jīng)濟(jì)成為大勢所趨，而如何安全且經(jīng)濟(jì)地儲存氫則是其中關(guān)鍵．近年來，許多科研小組從實(shí)驗、理論方面對各種儲氫材料對進(jìn)行了大量的工作，且取得了大量的成果．其中鎂基氫化物以其較輕的質(zhì)量、較高的儲氫量、資源豐富和良好的吸放氫等優(yōu)勢，被視為最具發(fā)展?jié)摿Φ膬洳牧现唬?，2］．其中最簡單的氫化物形式MgH2含氫量達(dá)到7.67 wt%［3］．但其存在著吸/放氫動力學(xué)性能差，放氫溫度在570 K以外，且用作二次化學(xué)電源負(fù)極材料時易腐蝕而導(dǎo)致循環(huán)穩(wěn)定性差，這些缺點(diǎn)使得其在實(shí)際應(yīng)用中受到一定的限制［4］．如何認(rèn)識進(jìn)而改良這些缺點(diǎn)也是科研人員對鎂基儲氫材料研究的關(guān)鍵點(diǎn)．

近年來，科研人員對各種各樣的鎂基合金進(jìn)行了大量的研究，以求改進(jìn)其吸附氫的動力學(xué)性能和活躍性能［5，6］．其中Mg2NiH4因其較高的含氫量(3.6 wt%)及較快的吸放氫速度，被認(rèn)為很具發(fā)展?jié)摿Φ牟牧希?］．針對Mg2NiH4的理論研究也已經(jīng)取得不錯的結(jié)果．如García等［8］對高溫相的Mg2NiH4進(jìn)行了大量理論研究，計算得到了氫原子在高溫相中所在的位置，Myers等［9］計算得到了低溫相的Mg2NiH4的晶體結(jié)構(gòu)，討論了其熱生成焓、電子結(jié)構(gòu)及光學(xué)屬性．Mg2NiH4被認(rèn)為是一種帯隙約為1.7—2.0 eV的半導(dǎo)體［10，11］，其氫吸附焓值約為0.70 eV/H2

［12］，比MgH2的氫吸附焓值(0.77 eV/ H2)稍低，因而其放氫溫度仍然須達(dá)到到510 K左右［13］．這么高的放氫溫度在實(shí)際應(yīng)用中就受到一定的限制，如應(yīng)用于聚合物電解質(zhì)膜燃料電池(PEM)這個溫度就不適合了．因此，如何進(jìn)一步改善鎂基儲氫合金的動力學(xué)性能是當(dāng)前儲氫材料研究的重點(diǎn)．

本文通過基于密度泛函理論(DFT)的第一性原理平面波贗勢方法［14］，模擬計算得出了Mg2Ni及其低溫相氫化物(LT-Mg2NiH4)和高溫相的氫化物(HT-Mg2NiH4)的電子結(jié)構(gòu)，并且在各個數(shù)據(jù)層面比對了高溫和低溫相氫化物的穩(wěn)定性及其影響因素．

2.計算模型與方法

作為鎂基儲氫合金之一的Mg2Ni具有良好的儲氫能力，其在吸氫時首先會形成結(jié)構(gòu)與Mg2Ni基本一致的α-H0.3Mg2Ni，只是晶胞體積會發(fā)生微小的膨脹［9］．在高溫和高壓下再形成Mg2NiH4，不同條件下可以形成高溫相氫化物(HT-Mg2NiH4)和低溫相氫化物(LT-Mg2NiH4)．我們對Mg2Ni及其氫化物進(jìn)行了充分的結(jié)構(gòu)優(yōu)化，并建立計算模型，分別得出Mg2Ni，HT-Mg2NiH4及LT-Mg2NiH4的晶格參數(shù)，計算結(jié)果見表1．

Mg2Ni的原子結(jié)構(gòu)為六方晶系(如圖1(a))，其原胞具有高對稱的D46，空間群為P6222 (No.180)［15］．晶格參數(shù)為a=5.219，c=13.293［16］．每個原胞包含12個鎂原子和6個鎳原子，可以分為平行于xy平面6個層面．每隔一層鎂原子中有2個鎳原子分別占據(jù)著3 b和3 d位，同時每隔一層鎳原子中有4個鎂原子分別占據(jù)著6 f和6 i位．

在室溫條件下(298 K)，Mg2NiH4為低溫相(LT-Mg2NiH4)，在X射線和中子衍射實(shí)驗中顯示單斜晶體結(jié)構(gòu)．(如圖1(b))其單位晶胞空間群為C2/ c(No.15)，晶格常數(shù)a=14.343，b=6.4038，c= 6.4830．每個晶胞包含16個鎂原子、8個鎳原子和32個氫原子．其中氫原子(8 f)具四種相互獨(dú)立的空間對稱性，且在鎳原子(8 f)周圍顯四面體方位排布．

當(dāng)溫度在553 K時，Mg2NiH4為高溫相(HTMg2NiH4)，在中子衍射中顯現(xiàn)面心立方結(jié)構(gòu)［17，18］．(如圖1(b))其原胞空間群為Fm3m (No.255)．晶格常數(shù)a=6.507．每個原胞包含8個鎂原子、4個鎳原子和16個氫原子．金屬原子顯反螢石結(jié)構(gòu)排布．在衍射實(shí)驗中，可以測出鎂原子和鎳原子分別占據(jù)著8 c位和4a位，而由于NiH4原子團(tuán)在Fm3 m群中的不確定性及鎂和鎳在中子散射中的長度關(guān)系，從而通過衍射方法不能確定氫原子的具體位置．García等通過計算氫原子不同排布的HT-Mg2NiH4的能量，發(fā)現(xiàn)其在費(fèi)米能級周圍的能帶分布跟氫原子在鎳原子周圍的排布情況有很大的相關(guān)系［19］．當(dāng)氫原子在鎳原子周圍顯四面體方位排布時，體系能量為最低，且會在費(fèi)米能級周圍產(chǎn)生個虛假的能帶結(jié)構(gòu)．甚至氫原子顯正四面體排布時，能使得體系顯示半導(dǎo)體屬性．自從Mg2Ni氫化物被認(rèn)為具有半導(dǎo)體屬性，我們對HTMg2NiH4所進(jìn)行的計算都是把氫原子看作是顯正四面體排布的，可以近似認(rèn)為氫原子占據(jù)96 j位［20］．

另外，為了便于低溫相和高溫相的Mg2NiH4實(shí)驗數(shù)據(jù)更好的進(jìn)行對比，我們對HT-Mg2NiH4建立了如圖1(c)所示的超胞，從而使得實(shí)驗計算中的HT-Mg2NiH4也具有16個鎂原子、8個鎳原子和32個氫原子．

本文采用基于密度泛函理論的第一性原理方法和超軟贗勢平面波方法［13］對Mg2Ni及其氫化物的幾何結(jié)構(gòu)和電子結(jié)構(gòu)進(jìn)行了理論研究．文章所有的計算都是采用VASP(Vienna an initio simulation package)程序包進(jìn)行計算的［21］．計算中電子結(jié)構(gòu)用密度泛函理論和廣義梯度近似描述［22］．贗勢由平面波綴加贗勢(PAW)來描述，電子之間的交換關(guān)聯(lián)能用Perdew-Wang(PW91)的廣義梯度近似(GGA)［23］，這種近似被廣泛應(yīng)用于密度泛函理論的第一性原理計算中［23，24］．計算過程中對Mg2Ni的高溫氫化物(HT-Mg2NiH4)和低溫氫化物(LTMg2NiH4)使用PAW-GGA法及第一性原理贗勢法，對Kohn-Sham方程泛函進(jìn)行自洽求解［25］．在對Mg2Ni，HT-Mg2NiH4和LT-Mg2NiH4進(jìn)行靜態(tài)計算之前，對所有的晶格參數(shù)和原子位置都進(jìn)行了充分的弛豫，直至每個晶胞上的能量收斂于10－5eV，受力收斂于0.05 eV/．自洽迭代過程中的簡約布里淵區(qū)分別取為7×7×3(Mg2Ni)，2×4×4(LTMg2NiH4)和5×5×5(HT-Mg2NiH4)．對費(fèi)米能采用高斯展寬，展開寬度為0.2 eV，平面波截斷能取為520 eV．文章所有計算結(jié)果都未考慮自旋的影響．

3.結(jié)果與討論

3.1.電子結(jié)構(gòu)分析

圖2是吸氫前的Mg2Ni的總態(tài)密度和軌道投影態(tài)密度，費(fèi)米能級處設(shè)為零eV．從圖中我們可以看出，價電子的能級范圍大約從Fermi能級以下7.5 eV一直延續(xù)到Fermi能級0eV處．Ni-3 d電子局域性較強(qiáng)，主要集中在靠近Fermi能級的－3 eV到0 eV．而Mg的3 s和2 p電子的色散比較大，能量范圍從－7.5 eV一直延伸到0 eV．從圖中，我們可以看出，Mg2Ni電子結(jié)構(gòu)中存在較強(qiáng)的Mg和Ni的軌道雜化，并因此而導(dǎo)致Mg和Ni間有較強(qiáng)的相互作用．另外，Mg的s電子的主體部分能級范圍明顯比Mg的p電子能級的主體部分的能級分布偏低，Mg和Ni之間的相互作用主要存在Ni-3 d軌道和Mg-p軌道間．

圖3 是吸氫以后LT-Mg2NiH4和HT-Mg2NiH4的軌道投影態(tài)密度和總態(tài)密度分布．從圖中我們可以看出，由于H原子的介入，價電子的能級范圍延伸到了Fermi能級以下9.5 eV的地方．這主要是由于H原子能級相對于Ni-3 d和Mg-3s電子能級都較低造成的．非常有意思的是，由于H的介入，材料的電子結(jié)構(gòu)由金屬性轉(zhuǎn)變成了絕緣體，其帶隙寬度大約為1.4 eV．這和實(shí)驗觀察到的結(jié)果相一致［10，11］．這種電子結(jié)構(gòu)的轉(zhuǎn)變，伴隨著Ni和Mg之間的電子轉(zhuǎn)移．圖4給出了Mg2Ni吸附H前后的電荷密度．從圖中可以明顯看出，吸附H前，Mg原子周圍存在一定的電子云，而在吸附H后，Mg上面電子密度幾乎為0．這種電子轉(zhuǎn)移發(fā)生后，對體系電子結(jié)構(gòu)產(chǎn)生了明顯的變化:Ni-3 d從Mg原子上轉(zhuǎn)移過來得到了2個電子，電子構(gòu)型變成了Ni-3d10，從而Ni的3 d軌道變成滿帶．同時，Ni的4s電子和H的s電子之間形成比較強(qiáng)的共價鍵，其成鍵態(tài)能級在Ni-3 d軌道下方．這樣，總體上體系的帶隙由Ni-3d滿帶和Ni-4p空帶之間形成．另外，由于H-1s和Ni-4s電子之間形成的成鍵態(tài)的能級發(fā)生了下降，在它和Ni-3 d電子態(tài)之間也形成了一個明顯的能隙．

相較于Mg2Ni，氫化后的態(tài)密度中顯示的Mg和Ni軌道雜化在有些區(qū)域明顯減弱，這主要由于Mg轉(zhuǎn)移電子到Ni-3d軌道上后，Ni-3d10電子構(gòu)型的局域性變得更強(qiáng)．Ni—H鍵的形成，也減弱了Mg和Ni之間的相互作用．比較LT-Mg2NiH4和HT-Mg2NiH4的態(tài)密度分布，發(fā)現(xiàn)兩者成鍵峰的分布的能量范圍有稍許區(qū)別，但總體上兩者總態(tài)密度圖形狀非常類似，都存在兩處能隙，說明H在高溫和低溫相里面，電子之間的相互作用和成鍵沒有明顯的差別．仔細(xì)分析原子間投影態(tài)密度，發(fā)現(xiàn)HT-Mg2NiH4比LTMg2NiH4的能隙要更加大些，這說明低溫相Ni-3 d電子的局域性相對被減弱了．而這種Ni-3 d電子局域性的減弱，是由于Ni-3 d電子和Mg的s電子之間存在一定的相互作用．如圖3(b)所示，在－2.5到－3 eV左右，Mg的s電子有一個比較明顯的峰存在(如圖中*號標(biāo)記)．

電荷密度分布圖可以直觀地看出原子間的相互作用．圖4展示了Mg2Ni，LT-Mg2NiH4及HTMg2NiH4的等電荷密度分布．從圖4(a)中我們可以發(fā)現(xiàn)Mg上面明顯的有一定的電子分布，這表明在Mg2Ni中Mg原子的電子沒有完全轉(zhuǎn)移到Ni的d軌道上，并且這些電子和Ni原子間存在一定的相互作用，這與態(tài)密度圖中得到的結(jié)論相一致．從圖4(b)和(c)中我們可以發(fā)現(xiàn)，在LT/HT-Mg2NiH4中，Ni原子與H原子間電子云有較大的重合，而Mg原子上基本上沒有電荷分布，這說明Ni原子與H間的相互作用明顯比Mg原子與H原子間的相互作用強(qiáng)．如前面對電子結(jié)構(gòu)分析，其主要原因在于Ni-4 s和H-1 s之間形成了比較強(qiáng)烈的共價鍵作用．Takahashi和Jasen等也發(fā)現(xiàn)在LT/HT-Mg2NiH4中Ni與H的成鍵相互作用遠(yuǎn)強(qiáng)于Mg與H之間的相互作用，這與我們的計算結(jié)果相一致［26］．

3.2.LT/HT-Mg2NiH4的相對穩(wěn)定性以及對儲氫性能的影響

從上面對Mg2Ni及其氫化物態(tài)密度的分析，我們可以看出，由于H的介入，導(dǎo)致Mg和Ni之間的相互作用減弱，而Ni和H之間形成了比較強(qiáng)的相互作用．相對于MgH2儲氫，Mg2Ni合金的優(yōu)勢在于Ni和H的相互作用，使得整體上來說，Mg2Ni合金氫化物的脫氫性能更好，進(jìn)而降低脫H溫度．另外，由于在低溫相中，Mg的s軌道和Ni的d軌道之間存在的比較弱的雜化，導(dǎo)致LT-Mg2NiH4的穩(wěn)定性更好，從而高溫相的脫H性能更好．但是我們分析出LT-Mg2NiH4比HT-Mg2NiH4更加穩(wěn)定的結(jié)論只是從比較兩者態(tài)密度曲線分布得出的，并沒有從宏觀數(shù)據(jù)上進(jìn)行比較．為了更進(jìn)一步比較兩者哪種更容易脫氫，我們分別計算了從LT-Mg2NiH4和HTMg2NiH4的Ni原子周圍去掉一個氫原子所需要的能量，我們可以用如下方程式進(jìn)行計算:

ΔE表示從Mg2NiH4中脫掉一個H原子所需要的能量，Etot(Mg16Ni8H32)和Etot(Mg16Ni8H31)分別表示Mg16Ni8H32晶體的總能量及脫掉一個氫原子后的總能量，Etot(H)表示單個氫原子的總能量．從上述定義，我們計算出ΔELT，ΔEHT的值分別為1.20 eV和0.64 eV，如表2所示．我們從中可以看出從LTMg2NiH4脫掉一個氫原子所需的能量比從HTMg2NiH4脫掉一個氫原子所需的能量要高，也就是說HT-Mg2NiH4要更加容易脫掉氫原子，這和實(shí)驗所得的高溫相適合脫氫和低溫相適合儲氫的結(jié)論一致．從LT到HT的過程是個穩(wěn)定性下降的過程，這跟從態(tài)密度圖中分析所得的結(jié)果一致．

4.結(jié)論

本文我們從態(tài)密度分布、電荷密度分布及能量等方面對Mg2Ni及其低溫和高溫的Mg2NiH4的結(jié)構(gòu)進(jìn)行了研究．所有結(jié)果都是基于密度泛函理論的第一性原理計算得出的．計算中弛豫得到的晶格常數(shù)和原子位置與從X射線和中子衍射實(shí)驗得到的結(jié)果符合較好．在Mg2Ni中，Mg和Ni的電子結(jié)構(gòu)存在較強(qiáng)的相互作用．H原子的加入使得Mg的價電子向Ni的3 d軌道轉(zhuǎn)移，并在Ni-3 d軌道和Ni-4 p軌道間形成了帶隙．Ni-4s和H-1 s電子之間的成鍵作用，使得Mg和Ni間相互作用減弱．HT-Mg2NiH4比LT-Mg2NiH4的能隙要更加大些，這說明低溫相Ni-3d電子的局域性相對被減弱了．而這種Ni-3 d電子局域性的減弱，是由于Ni-3d電子和Mg的s電子之間存在一定的相互作用．從LT-Mg2NiH4中去掉一個H原子比HT-Mg2NiH4中去掉一個H原子需要更多的能量，表明高溫相有利于脫氫反應(yīng)，而低溫相更適合儲氫反應(yīng)．

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

［2］Reilly J J，Wiswall R H 1968 Inorg Chem．7 2254

［3］Han X Q，Jiang L J，Liu Y F 2010 Acta Phys．Sin．59 4542(in Chinese)［韓曉琴、蔣利娟、劉玉芳2010物理學(xué)報59 4542］

［4］García G N，Abriata J P 1999 Phys．Rev．B 59 11746

［5］Soubeyroux J L，F(xiàn)ruchart D，Mikou A，Pezat M，Dar-riet B1984 Mater．Res．Bull 19 1119

［6］Zolliker P，Yvon K，Jorgensen J D，Rotella F J 1986 Inorg．Chem．25 3590

［7］Zeng K J，Klassen T，Oelerich W，Bormann R 1999 J．Alloys Compound．28 3151

［8］Myer W R，Wang L W，Richardson T J，Rubin M D 2002 Appl．Phys．91 4879

［9］Blomqvist H，R nnebro E，Noréus D，Kuji T 2002 Alloys Compound．268 330

［10］Post M L，Murray J 1987 Less-Common Met．134 15

［11］Liao M S，Zhang Q E 1997 Inorg．Chem．36 396

［12］Myers W R，Wang L W，Richardson T J，Rubin M D 2002 Appl．Phys．91 4879

［13］Richardson T J，Slack J L，Armitage R D，Kostecki R，F(xiàn)arangis B，Rubin M D 2001 Appl．Phys．78 3047

［14］Ouyang C Y，Xiong Z H，Ouyang Q Z，Liu G D，Ye Z Q，Lei M S 2006 Chin．Phys．15 1585

［15］Simiˇcic＇M V，Zdujic＇M，Dimitrijevic＇R，Nikolic＇-Bujanpvic＇Lj，Popovic＇N H 2006 J．Power Sources 158 730

［16］Zolliker P，Yvon K，Jorgensen J D，Rotella F J 1986 Inorg．Chem．25 3590

［17］Yvon K，Stucki J 1981 Inorg．Chem．20 2776

［18］Noréus D，Olsson L G 1983 J．Chem．Phys．78 2419

［19］García G N，Abriata J P，Sofo J O 1999 Phys．Rev．B 59 11746

［20］García G N，Abriata J P，Sofo J O 1999 Phys．Rev．B 59 11754

［21］Kresse G，Hafner J 1993 Phys．Rev．B 47 558

［22］Blochl P E 1994 Phys．Rev．B 50 17953

［23］Perdew J P，Chevary J A 1992 Phys．Rev．B 46 6671

［24］Zhong Z Y，Nie Z X，Du Y L，Ouyang C Y，Shi S Q，Lei M S 2009 Chin．Phys．18 2492

［25］Perdew J P 1992 Phys．Rev．B 46 6671

［26］Takahashi Y，Yukawa H，Morinaga M 1996 J．Alloys Compound．242 98

PACS:71．15．Mb，31．15．es，31．10．+z

*Project supported by the National Natural Science Foundation of China(Grant Nos．11064004 and 10864002)and Natural Science Foundation of Jiangxi Province，China(Grant Nos．2010 GZW0028)．

E-mail:chunh_liu@163．com

First principles investigation of electronic structures and stabilities of Mg2Ni and its complex hydrides*

Liu Chun-HuaOuyang Chu-Ying Ji Ying-Hua
(Department of Physics of Jiangxi Normal University，Key Laboratory of Photo Electronic＆Telecommunication of Jiangxi Province，Nanchang 330022，China)
(Received 8 July 2010;revised manuscript received 9 October 2010)

First-principle calculations are carried out to investigate the energies and the electronic structures of Mg2Ni alloy and its high/low temperature hydrides．Results show that strong bonding interaction between Mg and Ni atoms exists in the Mg2Ni alloy．With H inserted，the charge transfer occurs from Mg to Ni-3d，which gives rise to a band gap between Ni-3 d and Ni-4 p orbitals．The interaction between Mg and Ni in LT/HT-Mg2NiH4hydrides is weakened due to the covalent bonding between Ni-4 s and H-1s．Therefore，it neads more energy to remove H atoms from LT-Mg2NiH4than from HTMg2NiH4，indicating that HT-Mg2NiH4is more suitable for de-hydrogenation reaction while the LT-Mg2NiH4is suitable for hydrogenation reaction．

first-principles，Mg2Ni，Mg2NiH4，stability of structures

*國家自然科學(xué)基金(批準(zhǔn)號:11064004和10864002)和江西省自然科學(xué)基金(批準(zhǔn)號:2010 GZW0028)資助的課題．

E-mail:chunh_liu@163．com