劉春華歐陽楚英嵇英華

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

劉春華歐陽楚英 嵇英華

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

對(duì)Mg2Ni及其氫化物的能量和電子結(jié)構(gòu)進(jìn)行了第一性原理計(jì)算，并對(duì)Mg2Ni低溫氫化物和高溫氫化物的穩(wěn)定性進(jìn)行了分析．結(jié)果發(fā)現(xiàn):在Mg2Ni中，Mg與Ni存在較強(qiáng)的相互作用．H原子的加入使得Mg的價(jià)電子向Ni的3d軌道轉(zhuǎn)移，并在Ni-3d軌道和Ni-4p軌道間形成了帶隙．Ni-4s和H-1 s電子之間的成鍵作用，使得Mg和Ni間相互作用減弱了;從LT-Mg2NiH4中去掉一個(gè)H原子比HT-Mg2NiH4中去掉一個(gè)H原子需要更多的能量，表明高溫相有利于脫氫反應(yīng)，而低溫相更適合儲(chǔ)氫反應(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ì)地儲(chǔ)存氫則是其中關(guān)鍵．近年來，許多科研小組從實(shí)驗(yàn)、理論方面對(duì)各種儲(chǔ)氫材料對(duì)進(jìn)行了大量的工作，且取得了大量的成果．其中鎂基氫化物以其較輕的質(zhì)量、較高的儲(chǔ)氫量、資源豐富和良好的吸放氫等優(yōu)勢，被視為最具發(fā)展?jié)摿Φ膬?chǔ)氫材料之一［1，2］．其中最簡單的氫化物形式MgH2含氫量達(dá)到7.67 wt%［3］．但其存在著吸/放氫動(dòng)力學(xué)性能差，放氫溫度在570 K以外，且用作二次化學(xué)電源負(fù)極材料時(shí)易腐蝕而導(dǎo)致循環(huán)穩(wěn)定性差，這些缺點(diǎn)使得其在實(shí)際應(yīng)用中受到一定的限制［4］．如何認(rèn)識(shí)進(jìn)而改良這些缺點(diǎn)也是科研人員對(duì)鎂基儲(chǔ)氫材料研究的關(guān)鍵點(diǎn)．

近年來，科研人員對(duì)各種各樣的鎂基合金進(jìn)行了大量的研究，以求改進(jìn)其吸附氫的動(dòng)力學(xué)性能和活躍性能［5，6］．其中Mg2NiH4因其較高的含氫量(3.6 wt%)及較快的吸放氫速度，被認(rèn)為很具發(fā)展?jié)摿Φ牟牧希?］．針對(duì)Mg2NiH4的理論研究也已經(jīng)取得不錯(cuò)的結(jié)果．如García等［8］對(duì)高溫相的Mg2NiH4進(jìn)行了大量理論研究，計(jì)算得到了氫原子在高溫相中所在的位置，Myers等［9］計(jì)算得到了低溫相的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)這個(gè)溫度就不適合了．因此，如何進(jìn)一步改善鎂基儲(chǔ)氫合金的動(dòng)力學(xué)性能是當(dāng)前儲(chǔ)氫材料研究的重點(diǎn)．

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

2.計(jì)算模型與方法

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

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

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

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

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

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

3.結(jié)果與討論

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

圖2是吸氫前的Mg2Ni的總態(tài)密度和軌道投影態(tài)密度，費(fèi)米能級(jí)處設(shè)為零eV．從圖中我們可以看出，價(jià)電子的能級(jí)范圍大約從Fermi能級(jí)以下7.5 eV一直延續(xù)到Fermi能級(jí)0eV處．Ni-3 d電子局域性較強(qiáng)，主要集中在靠近Fermi能級(jí)的－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電子的主體部分能級(jí)范圍明顯比Mg的p電子能級(jí)的主體部分的能級(jí)分布偏低，Mg和Ni之間的相互作用主要存在Ni-3 d軌道和Mg-p軌道間．

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

相較于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電子的局域性相對(duì)被減弱了．而這種Ni-3 d電子局域性的減弱，是由于Ni-3 d電子和Mg的s電子之間存在一定的相互作用．如圖3(b)所示，在－2.5到－3 eV左右，Mg的s電子有一個(gè)比較明顯的峰存在(如圖中*號(hào)標(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)．如前面對(duì)電子結(jié)構(gòu)分析，其主要原因在于Ni-4 s和H-1 s之間形成了比較強(qiáng)烈的共價(jià)鍵作用．Takahashi和Jasen等也發(fā)現(xiàn)在LT/HT-Mg2NiH4中Ni與H的成鍵相互作用遠(yuǎn)強(qiáng)于Mg與H之間的相互作用，這與我們的計(jì)算結(jié)果相一致［26］．

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

從上面對(duì)Mg2Ni及其氫化物態(tài)密度的分析，我們可以看出，由于H的介入，導(dǎo)致Mg和Ni之間的相互作用減弱，而Ni和H之間形成了比較強(qiáng)的相互作用．相對(duì)于MgH2儲(chǔ)氫，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)一步比較兩者哪種更容易脫氫，我們分別計(jì)算了從LT-Mg2NiH4和HTMg2NiH4的Ni原子周圍去掉一個(gè)氫原子所需要的能量，我們可以用如下方程式進(jìn)行計(jì)算:

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

4.結(jié)論

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

E-mail:chunh_liu@163．com