高壓下鍺化鎂的金屬化相變研究?

王君龍 張林基 劉其軍 陳元正 沈如 何竹唐斌 劉秀茹?

1)(西南交通大學(xué)物理科學(xué)與技術(shù)學(xué)院,材料先進(jìn)技術(shù)教育部重點實驗室,成都 610031)

2)(西北工業(yè)大學(xué)材料學(xué)院,西安 710072)

高壓下鍺化鎂的金屬化相變研究?

王君龍1)張林基1)劉其軍1)陳元正1)沈如1)何竹1)唐斌2)劉秀茹1)?

1)(西南交通大學(xué)物理科學(xué)與技術(shù)學(xué)院,材料先進(jìn)技術(shù)教育部重點實驗室,成都 610031)

2)(西北工業(yè)大學(xué)材料學(xué)院,西安 710072)

(2017年3月10日收到;2017年6月7日收到修改稿)

鍺化鎂是一種窄帶半導(dǎo)體,壓力作用可以使鍺化鎂導(dǎo)帶底與價帶頂?shù)哪芟蹲冃?本文基于第一性原理計算了鍺化鎂在高壓下的能帶結(jié)構(gòu)以及反螢石相(常壓穩(wěn)定相)和反氯鉛礦相(高壓相)的焓值,發(fā)現(xiàn)在7.5 GPa時反螢石結(jié)構(gòu)鍺化鎂導(dǎo)帶底與價帶頂?shù)哪芟堕]合,預(yù)示著半導(dǎo)體相轉(zhuǎn)變?yōu)榻饘傧?計算結(jié)果還預(yù)測在11.0 GPa時鍺化鎂發(fā)生從反螢石結(jié)構(gòu)到反氯鉛礦結(jié)構(gòu)的相變.實驗研究方面,本文采用長條形壓砧在連續(xù)加壓條件下測量了鍺化鎂高壓下的電阻變化,采用金剛石對頂壓砧測量了鍺化鎂的高壓原位拉曼光譜,發(fā)現(xiàn)在8.7 GPa鍺化鎂的電阻出現(xiàn)不連續(xù)變化,9.8 GPa以上鍺化鎂的拉曼振動峰消失.由于金屬相的自由電子濃度高會阻礙激發(fā)光進(jìn)入樣品,進(jìn)而引起拉曼振動峰消失,因此我們推測鍺化鎂在9.8 GPa轉(zhuǎn)變?yōu)榻饘傧?

鍺化鎂,金屬化相變,高壓拉曼,第一性原理計算

1 引 言

金屬間化合物鍺化鎂是一種具有良好熱電性能、強(qiáng)抗腐蝕能力的窄帶半導(dǎo)體[1].元素鍺和鎂的原料豐富、價格較低廉、無毒,鍺化鎂作為一種環(huán)保型功能材料引起越來越多的關(guān)注和研究[2?5].鍺化鎂具有高的塞貝克系數(shù)、低的電阻率,被認(rèn)為是新一代熱電材料之一[1,6].鍺化鎂的帶隙寬度約為0.67 eV,在紅外探測器方面具有潛在應(yīng)用[7,8].壓力可以調(diào)控物質(zhì)的結(jié)構(gòu)、改變物質(zhì)的性質(zhì),在高壓條件下物質(zhì)內(nèi)部分子之間、原子之間的距離被壓縮,化學(xué)鍵、能級也會發(fā)生變化.當(dāng)原子發(fā)生重新排列,物質(zhì)將以新的高壓相存在[9?13].理論研究表明鍺化鎂的能隙在高壓下變小,當(dāng)能隙減小為零、價帶和導(dǎo)帶重合時,鍺化鎂將由半導(dǎo)體相轉(zhuǎn)變?yōu)榻饘傧?即在高壓下會發(fā)生金屬化相變[14,15].Kalarasse和Bennecer[14]用第一性原理全勢線性綴加平面波法預(yù)測反螢石結(jié)構(gòu)(anti-fluorite)的鍺化鎂在6.7 GPa時發(fā)生金屬化相變.Yu等[15]利用第一性原理計算了鍺化鎂在高壓下的結(jié)構(gòu),預(yù)測鍺化鎂在8.71 GPa和33.28 GPa將發(fā)生兩個結(jié)構(gòu)相變,分別是從反螢石結(jié)構(gòu)到反氯鉛礦結(jié)構(gòu)(anticotunnite)、反氯鉛礦結(jié)構(gòu)到Ni2In結(jié)構(gòu)的轉(zhuǎn)變,并且預(yù)測兩個高壓相均是金屬相.為了驗證理論預(yù)言,Li等[16]使用金剛石壓砧(DAC)研究了高壓下鍺化鎂的電阻率和霍爾系數(shù)隨壓力的變化關(guān)系,結(jié)合不同壓力下鍺化鎂的電阻率隨溫度的變化關(guān)系,認(rèn)為在7.4GPa反螢石結(jié)構(gòu)的鍺化鎂發(fā)生了金屬化相變,在9.5 GPa和35.8 GPa電阻率的不連續(xù)變化對應(yīng)于理論預(yù)言的兩個結(jié)構(gòu)相變,即從反螢石結(jié)構(gòu)到反氯鉛礦結(jié)構(gòu)、再到Ni2In結(jié)構(gòu)的轉(zhuǎn)變.為了研究高壓下鍺化鎂的金屬化相變和結(jié)構(gòu)相變,本文采用新型長條形壓砧測量連續(xù)加壓條件下鍺化鎂的電阻隨壓力的變化曲線,利用原位高壓拉曼光譜技術(shù)測量高壓下鍺化鎂的拉曼光譜圖,基于第一性原理計算鍺化鎂高壓下的能帶結(jié)構(gòu)以及反螢石結(jié)構(gòu)和反氯鉛礦結(jié)構(gòu)的焓值,分析金屬化相變的機(jī)理.

2 實驗方法

純度為99.99%的鍺化鎂購于阿法埃莎公司(A lfa Aesar Corporation),我們對原始樣品進(jìn)行了X射線衍射分析(XRD,型號:Panalytical X’pert,Cu Kαradiation).高壓下的電阻測量在兩面頂液壓機(jī)上進(jìn)行,使用的高壓模具是自行設(shè)計的長條形平面對頂壓砧,壓砧實物圖如圖1(a)的內(nèi)插圖所示[17].相比于傳統(tǒng)的圓形砧面對頂壓砧(即Bridgman壓砧),長條形壓砧可在沿中心線的狹長區(qū)域內(nèi)產(chǎn)生較均勻分布的高壓力,有利于進(jìn)行高壓下樣品的電阻測量[17].長條形硬質(zhì)合金壓砧的砧面長20mm、寬5mm,封墊材料為葉臘石(長23mm、寬5.4 mm、厚度0.55 mm).鍺化鎂樣品填充在葉臘石封墊的凹槽內(nèi)(長8 mm、寬1 mm、深0.2 mm).本文采用四探針方法測量鍺化鎂樣品的電阻,樣品組裝示意圖如圖1(b)所示.四個探針為厚度為0.2 mm、寬1 mm的銅箔.采用恒流電源為測量電路提供恒定的電流100μA,使用記錄儀(型號:GRAPHTEC GL900)記錄連續(xù)加壓過程中樣品的電壓及壓機(jī)油壓,I+和I?表示測量電路中電流的方向,U+和U?表示被測試樣電壓的正負(fù)極.利用鉍和碲化鋅的相變對壓機(jī)油壓和樣品腔實際壓力進(jìn)行了標(biāo)定[18?20],結(jié)果如圖1和(1)式所示.根據(jù)(1)式將實驗過程中記錄的壓機(jī)油壓Poil(單位為MPa)換算成樣品腔實際壓力P(單位為GPa),根據(jù)歐姆定律將樣品的電壓換算成電阻,進(jìn)而獲得樣品的電阻隨壓力的變化曲線.

圖1 (網(wǎng)刊彩色)(a)長條形對頂壓砧的實物圖及壓力標(biāo)定曲線;(b)長條形壓砧上四探針法測電阻示意圖Fig.1.(color on line)(a)The picturesof strip opposite anvilsand the p ressure calibration cu rve;(b)scheMatic d iagraMof saMp le asseMb ly in a strip opposite anvil.

高壓下鍺化鎂的原位拉曼光譜測量使用的是DAC和具有長焦距物鏡的拉曼光譜儀(型號:Renishaw,in-Via),激發(fā)光波長為532 nm.DAC砧面直徑為300μm,封墊材料為T 301不銹鋼片,預(yù)壓后在封墊中心打孔(直徑約為110μm),未使用傳壓介質(zhì),樣品腔的壓力由紅寶石熒光法測定[21].

本文的第一性原理計算軟件采用Materials studio中的CASTEP模塊[22,23].計算過程采用周期性邊界條件,將多電子體系用平面波函數(shù)展開表示,采用超軟贗勢(ultrasoft)作為平面波基組[24].電子間相互作用采用廣義梯度近似(GGA)中Perdew-Burke-Ernzerhof交換關(guān)聯(lián)能函數(shù)[25].本文中利用Broyden-Fletcher-Gold farb-Shanno優(yōu)化方法來找到每個壓力下物質(zhì)的局域最穩(wěn)定結(jié)構(gòu)[26],平面截斷能取為500 eV.

3 實驗結(jié)果與討論

圖2(a)給出了鍺化鎂的X射線衍射譜.圖中衍射峰尖銳、強(qiáng)度高,并且沒有明顯的雜質(zhì)峰,說明原始樣品結(jié)晶良好、純度高.原始樣品為反螢石結(jié)構(gòu),空間群Fm-3m,晶格常數(shù)為a=6.377?,與文獻(xiàn)報道的鍺化鎂晶格常數(shù)實驗值及理論值一致[27?29].圖2(b)是鍺化鎂反螢石結(jié)構(gòu)的晶胞模型.

鍺化鎂的體積彈性模量為57.6GPa[8],在高壓下鍺化鎂體積會發(fā)生較大塌縮,因此壓力可能會對它的電子態(tài)密度產(chǎn)生明顯的影響.圖3(a)是鍺化鎂的電阻隨著壓力的變化曲線.加載壓力后,鍺化鎂的電阻開始急劇減小.我們認(rèn)為初始加壓階段電阻的急劇下降主要是疏松樣品的壓緊過程引起的,因為高電阻率的空氣在樣品壓緊過程中會被擠出.在2.0 GPa以上,鍺化鎂的電阻隨著壓力的升高而緩慢減小,這歸因于鍺化鎂導(dǎo)帶底與價帶頂?shù)哪芟峨S著壓力的升高而變小.對于本征半導(dǎo)體,電阻率ρ與能隙Eg的關(guān)系式為ρ～exp[Eg/(2kBT)],其中kB和T分別為玻爾茲曼常數(shù)和溫度[20].室溫下,隨著壓力的增加能隙減小,更多的載流子被激發(fā)到導(dǎo)帶,相應(yīng)地鍺化鎂的電阻率變小.為了更清楚地顯示高壓下低阻值范圍內(nèi)的電阻變化情況,將圖3(a)中5 GPa以上的電阻-壓力曲線用對數(shù)坐標(biāo)表示,如圖3(b)所示.圖3(b)中,在8.7 GPa左右鍺化鎂的電阻出現(xiàn)了不連續(xù)變化,電阻值迅速減小.如引言所述,Kalarasse和Bennecer[14]通過理論計算預(yù)測反螢石結(jié)構(gòu)的鍺化鎂在6.7 GPa能隙閉合,發(fā)生金屬化相變.Yu等[15]利用第一性原理計算了鍺化鎂在高壓下的結(jié)構(gòu),預(yù)測鍺化鎂在8.71 GPa時發(fā)生反螢石結(jié)構(gòu)到反氯鉛礦結(jié)構(gòu)的相變,而反氯鉛礦結(jié)構(gòu)的鍺化鎂為金屬相.因此,無論是反螢石結(jié)構(gòu)的鍺化鎂在高壓下發(fā)生能隙閉合,還是轉(zhuǎn)變?yōu)榻饘傧嗟姆绰茹U礦結(jié)構(gòu),理論計算的結(jié)果都預(yù)言鍺化鎂將在6.7—8.7 GPa轉(zhuǎn)變?yōu)榻饘傧?參照理論預(yù)言,我們推測鍺化鎂的電阻在8.7 GPa急劇下降對應(yīng)于金屬化轉(zhuǎn)變,該實驗結(jié)果比Li等[16]通過金剛石壓砧測量的鍺化鎂電阻不連續(xù)變化壓力值7.8 GPa略高,不同的壓力測量方法可能是造成相變壓力結(jié)果不同的因素之一.卸壓時,鍺化鎂的電阻隨著壓力的下降而升高,說明壓力引起的能隙減小是可逆的,與Li等的實驗結(jié)論一致.

圖2 (網(wǎng)刊彩色)(a)鍺化鎂原始樣品的X射線衍射圖譜;(b)反螢石結(jié)構(gòu)鍺化鎂的晶胞模型Fig.2.(color on line)(a)XRD pattern of Mg2Ge poWder;(b)crystal structu re of anti-fl uorite type Mg2Ge.

圖3 鍺化鎂的電阻隨壓力的變化曲線 (a)采用線性坐標(biāo);(b)采用對數(shù)坐標(biāo)Fig.3.Pressu re-dependent cu rves of electric resistance of Mg2Ge under(a)linear coordinate and(b)logarithMic coord inate.

高壓原位拉曼光譜法是研究高壓下物質(zhì)分子結(jié)構(gòu)變化的有效手段,也可以用來判斷物質(zhì)的金屬化相變.通常,金屬相物質(zhì)具有大量的自由載流子,會阻礙激發(fā)光進(jìn)入試樣,因此不能觀察到與試樣分子發(fā)生非彈性碰撞的拉曼散射光.為了進(jìn)一步研究鍺化鎂的壓致金屬化相變,我們測量了鍺化鎂在高壓下的拉曼光譜.常壓下反螢石結(jié)構(gòu)鍺化鎂的拉曼光譜主要有F2g和F1u兩個振動模式[1,30,31].其中,F2g振動模式的特征峰波數(shù)為255 cm?1,F1u振動模式的特征峰波數(shù)在545—550 cm?1范圍內(nèi)[1,25,26].圖4(a)給出了不同壓力下鍺化鎂的拉曼光譜圖,高壓原位拉曼光譜測量的最高壓力為21.1 GPa.圖4(b)是拉曼主峰即F2g振動模的強(qiáng)度隨壓力的變化趨勢.隨著壓力的升高,鍺化鎂拉曼光譜的強(qiáng)度減弱,這是由于隨著壓力的增加鍺化鎂導(dǎo)帶低和價帶頂?shù)哪芟稖p小,更多的電子被激發(fā)到導(dǎo)帶.自由電子的增多增加了對激發(fā)光的散射,減小了激發(fā)光與分子發(fā)生非彈性碰撞的概率.在9.8 GPa以上,拉曼主峰及拉曼光譜中的其他峰消失,這表明鍺化鎂發(fā)生了金屬化相變,因為金屬化以后大量的自由電子阻礙了激發(fā)光進(jìn)入試樣[32].Morozova等[32]用相同的方法判斷摻鋁的硅化鎂金屬化相變.圖4(c)給出了拉曼主峰的峰位隨壓力的變化趨勢.圖4(c)中,隨著壓力的增加,主峰的位置向高波數(shù)方向移動,在2.6GPa有一個不連續(xù)變化.由于高壓拉曼光譜實驗時未使用任何傳壓介質(zhì),因此主峰峰位在2.6 GPa左右的不連續(xù)變化可能是軸向應(yīng)力引起的[33].綜上,在9.8 GPa以上鍺化鎂拉曼振動峰的消失,我們推測是金屬化相變引起的.

圖4 (網(wǎng)刊彩色)(a)高壓下鍺化鎂的拉曼光譜;(b)F2g振動模拉曼峰的強(qiáng)度隨壓力的變化曲線;(c)F2g振動模拉曼峰的峰位置隨壓力的變化曲線Fig.4.(color on line)(a)RaMan spectra of Mg2Ge under high p ressures;(b)p ressure dependence of peak intensity of F2gphonon Mode;(c)p ressu re dependence of peak position of F2gphonon Mode.

為了研究金屬化相變的機(jī)理,本文采用第一性原理計算了反螢石結(jié)構(gòu)鍺化鎂的能隙隨壓力的變化.首先,通過結(jié)構(gòu)優(yōu)化得到了零壓下鍺化鎂的晶格常數(shù)a=6.439?,與實驗值6.377?接近,表明計算結(jié)果是合理的.鍺化鎂能帶圖譜如圖5(a)所示,鍺化鎂導(dǎo)帶的最小值位于X點,價帶最大值位于G點.零壓時的間接能隙為0.128 eV,文獻(xiàn)報道的能隙寬度實驗值為0.6—0.67 eV[7,34],計算值為0.101—0.166 eV[8,35],說明本文中的計算值與文獻(xiàn)比較接近.圖5(b)和圖5(c)分別為鍺化鎂的總態(tài)密度圖和分波態(tài)密度圖,圖中實線表示0 GPa,虛線表示7.0 GPa.從總的態(tài)密度圖可以看出,在零壓下鍺化鎂的價帶寬度是10.15 eV,比文獻(xiàn)報道的計算值9.98 eV略高[35].低價帶的能量區(qū)間為?10.15—?7.04 eV,高價帶的能量區(qū)間為?5.32 eV到價帶最大值.從分波態(tài)密度圖可以看出,在低價帶區(qū)間主要是Ge-s的貢獻(xiàn),其次是Mg-s和Mg-p的貢獻(xiàn);高價帶區(qū)間主要是Ge-p的貢獻(xiàn).為了研究高壓下鍺化鎂電子運(yùn)動狀態(tài),本文計算了高壓下鍺化鎂的能帶圖譜.隨著壓力的增加,價帶寬度發(fā)生了展寬.圖5(d)給出了導(dǎo)帶底與價帶頂?shù)哪芰坎钪惦S壓力的變化曲線.圖5(d)中,壓力約7.5 GPa時能隙寬度變?yōu)榱?預(yù)示著在該壓力以上電子可以填充到彼此重疊的能帶中去,鍺化鎂轉(zhuǎn)變?yōu)榻饘傧?該相變壓力比Kalarasse和Bennecer[14]報道的能隙閉合時的壓力值6.7 GPa略高,更接近Li等[16]報道的金屬化相變的實驗結(jié)果.

圖5 (網(wǎng)刊彩色)(a)鍺化鎂的能帶結(jié)構(gòu);(b)鍺化鎂的總態(tài)密度圖;(c)鍺化鎂的分波態(tài)密度圖;(d)鍺化鎂能隙隨壓力的變化曲線:圓形為本文數(shù)據(jù),方形為K alarasse等[14]的數(shù)據(jù)Fig.5.(color on line)(a)Band structu re of Mg2Ge;the density of states of Mg2Ge:(b)total,(c)partial;(d)the variation of band gap With p ressu re of Mg2Ge:the circle shoWs the resu lts of this work,the square rep resents the resu lts of K alarasse et al.[14].

圖6 鍺化鎂的反螢石結(jié)構(gòu)和反氯鉛礦結(jié)構(gòu)的焓差隨壓力的變化Fig.6.Pressu re-dependent cu rve of enthalpy diff erence between the anti-fl uorite and the anti-cotunnite structure of Mg2Ge.

通過第一性原理計算方法研究了鍺化鎂反螢石相和反氯鉛礦相的焓差隨壓力的變化,計算結(jié)果如圖6所示.從圖中可以看出,常壓下鍺化鎂最穩(wěn)定的結(jié)構(gòu)是反螢石結(jié)構(gòu),這與實驗結(jié)果相符.當(dāng)壓力達(dá)到11.0 GPa時,反螢石結(jié)構(gòu)將相變到反氯鉛礦結(jié)構(gòu).因此,理論計算結(jié)果表明,隨著壓力的增加,反螢石結(jié)構(gòu)先發(fā)生金屬化,進(jìn)而結(jié)構(gòu)相變?yōu)榉绰茹U礦結(jié)構(gòu),這與Li等[16]報道的金屬化相變的實驗結(jié)果基本一致.Li等使用金剛石壓砧研究了高壓下鍺化鎂的電阻率和霍爾系數(shù)隨壓力的變化關(guān)系,結(jié)合不同壓力下鍺化鎂的電阻率隨溫度的變化關(guān)系,認(rèn)為在7.4 GPa反螢石結(jié)構(gòu)的鍺化鎂發(fā)生了金屬化相變,在9.5 GPa發(fā)生了反螢石結(jié)構(gòu)到反氯鉛礦結(jié)構(gòu)的轉(zhuǎn)變[16].本文采用長條形壓砧和金剛石對頂壓砧分別測量了鍺化鎂高壓下的電阻變化以及原位高壓拉曼光譜,發(fā)現(xiàn)在8.7 GPa鍺化鎂的電阻出現(xiàn)不連續(xù)變化,9.8 GPa以上鍺化鎂的拉曼振動峰消失.結(jié)合理論計算與Li等[16]的實驗報道,我們推測高壓下導(dǎo)帶低和價帶頂?shù)哪芟堕]合是鍺化鎂電阻在8.7 GPa突然下降的原因,9.8 GPa以上鍺化鎂有可能轉(zhuǎn)變?yōu)榉绰茹U礦結(jié)構(gòu),反氯鉛礦結(jié)構(gòu)也是金屬相.

4 結(jié) 論

本文研究了鍺化鎂在高壓下的金屬化相變,基于第一性原理的理論計算表明,隨著壓力的增加,反螢石結(jié)構(gòu)鍺化鎂的導(dǎo)帶底與價帶頂?shù)哪芟吨饾u減小,在7.5 GPa時能隙閉合,預(yù)示著半導(dǎo)體轉(zhuǎn)變?yōu)榻饘傧?通過比較反螢石相和反氯鉛礦相的焓值,預(yù)測在11.0 GPa鍺化鎂將由反螢石結(jié)構(gòu)轉(zhuǎn)變?yōu)榉绰茹U礦結(jié)構(gòu).實驗研究發(fā)現(xiàn)鍺化鎂的電阻在8.7 GPa時發(fā)生了不連續(xù)變化,在9.8 GPa以上鍺化鎂的拉曼振動峰消失,以此推測鍺化鎂發(fā)生了金屬化,因為金屬化以后大量的自由電子阻礙了激發(fā)光進(jìn)入樣品.結(jié)合文獻(xiàn)報道的金剛石壓砧上鍺化鎂的電阻率、霍爾系數(shù)和電阻率隨溫度變化的實驗結(jié)果,即在7.4 GPa反螢石結(jié)構(gòu)鍺化鎂發(fā)生金屬化相變,在9.5 GPa發(fā)生反螢石結(jié)構(gòu)到金屬相的反氯鉛礦結(jié)構(gòu)轉(zhuǎn)變[16],我們推測高壓下能隙閉合是鍺化鎂在8.7 GPa電阻突然下降的原因,9.8 GPa以上鍺化鎂有可能轉(zhuǎn)變?yōu)榻饘傧嗟姆绰茹U礦結(jié)構(gòu).

感謝楊文革研究員在高壓原位拉曼光譜測量方面給予的支持,感謝郭春生教授在理論計算方面的建議.

[1]Tani J,K ido H 2008 CoMp.Mater.Sci.42 531

[2]Chung P L,Whitten WB,Danielson G C 1965 J.Phys.Chem.Solids 26 1753

[3]Guo S D 2016 Eur.Phys.J.B 89 1

[4]Liu H Y,Zhu Z Z,Yang Y 2008 Acta Phys.Sin.57 5182(in Chinese)[劉慧英,朱梓忠,楊勇 2008物理學(xué)報57 5182]

[5]Mao J,K iMH S,Shuai J,Liu Z,He R,SaparaMadu U,T ian F,Liu W,Ren Z 2016 Acta Mater.103 633

[6]Martin J J 1972 J.Phys.Chem.Solids 33 1139

[7]Stella A,Lynch D W1964 J.Phys.Chem.So lids 25 1253

[8]Corkill J L,Cohen ML 1993 Phys.Rev.B 48 17138

[9]Xu J A,Wang Y Y,Xu MH 1980 Acta Phys.Sin.29 1063(in Chinese)[徐濟(jì)安,王彥云,徐敏華1980物理學(xué)報29 1063]

[10]Wang J R,Zhu J,Hao Y J,Ji G F,X iang G,Zou Y C 2014 Acta Phys.Sin.63 186401(in Chinese)[王金榮,朱俊,郝彥軍,姬廣富,向鋼,鄒洋春2014物理學(xué)報63 186401]

[11]Jin C Q,Liu Q Q,Deng Z,Zhang S J,X ing L Y,Zhu J L,Kong P P,Wang X C 2013 Chin.J.H igh Press.Phys.27 473(in Chinese)[靳常青,劉青清,鄧正,張思佳,邢令義,朱金龍,孔盼盼,望賢成2013高壓物理學(xué)報27 473]

[12]Zhang S J,Wang X C,SamMynaiken R,Tse J S,Yang L X,Liu Q Q,Desgreniers S,Yao Y,Liu H Z,Jin C Q 2009 Phys.Rev.B 80 014506

[13]Zhang J L,Zhang S J,Weng H M,Zhang W,Yang L X,Liu Q Q,Feng SM,Wang X C,Yu R C,Cao L Z,Wang L,Yang WG,Liu H Z,Zhao WY,Zhang S C,Dai X,Fang Z,Jin C Q 2011 Proc.Natl.Acad.Sci.USA 108 24

[14]K alarasse F,Bennecer B 2008 J.Phys.Chem.Solids 69 1775

[15]Yu F,Sun J X,Chen T H 2011 Phys.B:Condens.Matter 406 1789

[16]Li Y,Gao Y,Han Y,Liu C,Peng G,Wang Q,Ke F,Ma Y,Gao C 2015 Appl.Phys.Lett.107 142103

[17]Tang F,Chen L Y,Liu X R,Wang J L,Zhang L J,Hong SM2016 Acta Phys.Sin.65 100701(in Chinese)[唐菲,陳麗英,劉秀茹,王君龍,張林基,洪時明2016物理學(xué)報65 100701]

[18]Getting IC 1998 Metro logia 35 119

[19]Ohtani A,MotobayashiM,Onodera A 1980 Phys.Lett.A 75 435

[20]Morozova N V,Ovsyannikov SV,Korobeinikov IV 2014 J.Appl.Phys.115 213705

[21]Mao H K,Xu J A,Bell P M1986 J.Geophys.Res 91 4673

[22]Payne MC,Teter MP,A llan D C 1992 Rev.Modern Phys.64 1045

[23]Segall MD,Lindan P J D,Probert MJ 2002 J.Phys.:Condens.Matter 14 2717

[24]Vanderbilt D 1990 Phys.Rev.B 41 7892

[25]PerdeWJP,Burke K,ErnzerhofM1996 Phys.Rev.Lett.77 3865

[26]Fischer T H,A lMlof J 1992 J.Phys.Chem.96 9768

[27]G rosch G H,Range K J 1996 J.A lloy.CoMpd.235 250

[28]Zhou D,Liu J,Xu S,Peng P 2012 CoMp.Mater.Sci.51 409

[29]Janot R,Cuevas F,Latroche M,Percheron-Guégan A 2006 In terMetallics 14 163

[30]Buchenauer C J,Cardona M1971 Phys.Rev.B 3 2504

[31]Anastassakis E,Perry C H 1971 Phys.Rev.B 4 1251

[32]Morozova N V,Ovsyannikov SV,Korobeinikov IV 2014 J.Appl.Phys.115 213705

[33]Mohiuddin T MG,LoMbardo A,Nair R R 2009 Phys.Rev.B 79 205433

[34]Stella A,Brothers A D,Hopkins R H 1967 Phys.Status So lidi 23 697

[35]Benhelal O,Chahed A,Laksari S 2005 Phys.Status Solidi 242 2022

PACS:62.50.–p,64.70.Kg,85.40.LsDOI:10.7498/aps.66.166201

*Project supported by the National Natural Science Foundation of China(G rant No.11004163)and the FundaMental Research Funds for Central Universities(G rant Nos.2682014ZT 31,2682016CX 065).

?Corresponding author.E-Mail:xrliu@sWjtu.edu.cn

P ressu re-induced Metallization transition in Mg2Ge?

Wang Jun-Long1)Zhang Lin-Ji1)Liu Qi-Jun1)Chen Yuan-Zheng1)Shen Ru1)He Zhu1)Tang Bin2)Liu Xiu-Ru1)?

1)(School of Physical Science and Technology,Key Laboratory of Advanced Technologies ofMaterials,Ministry of Education of China,Southwest Jiaotong University,Chengdu 610031,China;)
2)(School ofMaterials Science and Engineering,Northwestern Polytechnical University,X i’an 710072,China)

10 March 2017;revised Manuscrip t

7 June 2017)

Mg2Ge With anti-fluorite structure at ambient pressure is characterized as a narroWband seMiconductor and increasing p ressure resu lts in a decrease of the gap.In this work,the band structure of anti-fl uorite Mg2Ge under high pressure is studied by fi rst p rincip les calcu lations,which suggests that Mg2Ge becoMesMetallic at 7.5 GPa as a result of band gap closure.The enthalpy diff erence between anti-fl uorite phase and anti-cotunnite phase under high p ressure is calculated by the fi rst-p rincip les p lane-waveMethod Within the pseudopotential and generalized gradient app roxiMation.The results shoWthat Mg2Ge undergoes a phase transition froMthe anti-fl uorite structure to anti-cotunnite structure at 11.0 GPa.Then we investigate experimentally the pressure-induced metallization of Mg2Ge by electric resistance MeasureMent in strip anvil cell and RaMan spectroscopy by diaMond anvil cell.The pressure distribution is homogeneous along the central line of the strip anvil and the p ressure is changed ccontinuously by using a hyd raulically d riven two-anvil press.Raman scattering experiment is performed at pressure up to 21.1 GPa on a back scattered Raman spectroMeter.Thewavelength of excitation laser is 532 nm.No pressure-transMitting is used and pressure is deterMined by the shift of the ruby luMinescence line.It is found that neither a discontinuous change of electrical resistance at 8.7GPa nor Raman vibrationmodes of Mg2Ge appear above 9.8GPa.The disappearance of the Raman vibrationmode is ascribed to theMetallization since the the free carrier concentration rises after Metallization has p revented the laser light froMpenetrating into the saMp le.We coMpare these results With those of resistivity MeasureMents in diaMond anvil cell.Li et al.[2015 Appl.Phys.Lett.107 142103]reported that Mg2Ge becomesmetallic phase at 7.4 GPa and is transforMed into Metallic anti-cotunnite phase at around 9.5 GPa.We specu late that the discontinuous change in electric resistance at 8.7 GPa is ascribed to the gap closure of anti-fluorite phase and Mg2Gemay transforMinto the anti-cotunnite phase above 9.8 GPa.

Mg2Ge,Metallization,Raman scattering under high pressure,fi rst princip les calculation

10.7498/aps.66.166201

?國家自然科學(xué)基金(批準(zhǔn)號:11004163)和中央高?；究蒲袠I(yè)務(wù)費(fèi)專項資金(批準(zhǔn)號:2682014ZT 31,2682016CX 065)資助的課題.

?通信作者.E-Mail:x rliu@sWjtu.edu.cn

?2017中國物理學(xué)會C h inese P hysica l Society

http://Wu lixb.iphy.ac.cn