溫姜霞，吳升聰，王 俊，陶小馬，陳紅梅，歐陽義芳

( 廣西大學(xué)物理科學(xué)與工程技術(shù)學(xué)院，廣西南寧 530004)

0 引言

近年來，針對具有A5B3型化合物的研究越來越多，這是因?yàn)樵擃惢衔锞哂歇?dú)特的物理化學(xué)性能[1-5]，研究表明超過40個金屬間化合物具有A5B3型結(jié)構(gòu)。常見的結(jié)構(gòu)類型有Cr5B3型、W5Si3型和Mn5Si3型這3類，Tao等[5-8]利用基于密度泛函的第一性原理方法對Ta5Si3、W5Si3、V5Si3以及Tl5Te3等化合物的物理化學(xué)性能進(jìn)行模擬計(jì)算研究，獲得了這些化合物的結(jié)構(gòu)穩(wěn)定性、力學(xué)性能、電子結(jié)構(gòu)特性以及熱物理性能和熱電性能等。含堿土金屬Ca的化合物中，Ca5Si3因其有潛在的儲氫性能[3]和熱電性能[9]備受關(guān)注。然而在含Ca的化合物中，Ca5Zn3也具有Cr5B3型結(jié)構(gòu)[1]。2012年，Yang等[10]利用基于密度泛函的第一性原理的方法對Ca-Zn體系的7個化合物進(jìn)行了計(jì)算研究，計(jì)算獲得了形成焓、彈性常數(shù)和電子態(tài)密度等信息。然而，至今還沒有針對Ca5Zn3的更加系統(tǒng)研究的報道。為了能夠獲得更加系統(tǒng)的物理性能，本文將利用基于密度泛函的第一性原理方法對Ca5Zn3化合物進(jìn)行系統(tǒng)計(jì)算研究，獲得該化合物的熱力學(xué)穩(wěn)定性、力學(xué)性能、電子結(jié)構(gòu)特性，并結(jié)合德拜模型對其熱物理性能進(jìn)行研究和預(yù)測。

1 材料與方法

1.1 Ca5Zn3的晶體結(jié)構(gòu)

采用Cr5B3型、W5Si3型和Mn5Si3型等3種結(jié)構(gòu)來討論Ca5Zn3結(jié)構(gòu)穩(wěn)定性。其中Cr5B3型結(jié)構(gòu)是一個四方晶系的結(jié)構(gòu)，Pearson符號為tI32，共有4種不同占位，Cr占據(jù)16l和4c位，B占據(jù)4a和8h位。每個Cr5B3型結(jié)構(gòu)單胞中含有10個Cr原子，6個B原子。其余兩種結(jié)構(gòu)的具體信息可參見文獻(xiàn)[5]。3種晶體結(jié)構(gòu)如圖1所示。

(a) Cr5B3型，(b) W5Si3型，(c) Mn5Si3型;藍(lán)色為Ca原子，灰色為Zn原子(a) Cr5B3 prototype，(b) W5Si3 prototype，(c) Mn5Si3 prototype;the blue atom is Ca，the grey atom is Zn圖1 Ca5Zn3化合物的3種晶體結(jié)構(gòu)Fig.1 Tree crystal structures of Ca5Zn3 compounds

1.2 計(jì)算方法

本文采用基于密度泛函理論(Density Functional Theory，DFT)的第一性原理計(jì)算軟件包Vienna Ab initio Simulation Package (VASP)[11]進(jìn)行計(jì)算。第一性原理計(jì)算方法是投影綴加波方法[12-13]，交換關(guān)聯(lián)勢采用廣義梯度近似的Perdew-Burke-Ernzerhof (PBE)方案[14]。布里淵區(qū)積分采用Monkhorst Pack布點(diǎn)方法[15]，對于Cr5B3和W5Si3結(jié)構(gòu)，K點(diǎn)取11×11×11，而對Mn5Si3結(jié)構(gòu)，K點(diǎn)為9×9×11。在本文計(jì)算中,波函數(shù)的截?cái)嗄茉O(shè)置為500 eV，能量收斂精度為10-6eV。在進(jìn)行性能計(jì)算前，所有結(jié)構(gòu)都進(jìn)行全優(yōu)化，直到原子間的力小于0.01 eV/?為止，以保證計(jì)算的準(zhǔn)確性。有關(guān)形成焓、彈性常數(shù)以及熱物理性質(zhì)的計(jì)算方法參見文獻(xiàn)[5]。

Ca5Zn3化合物形成焓可以用以下的公式進(jìn)行計(jì)算：

ΔH(Ca5Zn3)=Etotal(Ca5Zn3)-5Etotal(Ca)-3Etotal(Zn)，

(1)

其中，Etotal(Ca5Zn3)、Etotal(Ca)和Etotal(Zn)分別是Ca5Zn3合金、Ca以及Zn元素的總能量。

(2)

(3)

(4)

Cr5B3結(jié)構(gòu)類型的Ca5Zn3是四方晶系。四方晶系有6個獨(dú)立彈性常數(shù)，分別為C11,C12,C13,C33,C44和C66。通過下面6個形變矩陣可以求得相應(yīng)彈性常數(shù)(表1)。

表1 四方晶系形變矩陣及相應(yīng)形變能Table 1 Deformation matrix of the tetragonal system and the corresponding deformation energy

德拜溫度作為一個基本的物理參數(shù)，與熱容、彈性常數(shù)等有著密切聯(lián)系。一般來說，德拜溫度越高，化合物熔點(diǎn)越高，意味著原子間的相互作用越強(qiáng)。一般來說，德拜溫度(Θ)可以從平均聲速(um)進(jìn)行估算:

(5)

其中,h和kB分別是普朗克常數(shù)和玻爾茲曼常數(shù)，V是原子平均體積，在多晶體中的平均聲速(um)可以用以下公式計(jì)算:

(6)

其中，ut和ul分別是橫向聲速和縱向聲速，這兩種聲速可以通過以下公式獲得:

(7)

(8)

其中,G為剪切模量，B為體積模量，ρ為質(zhì)量密度。

2 結(jié)果與分析

2.1 結(jié)構(gòu)穩(wěn)定性能

為了獲得Ca5Zn3的結(jié)構(gòu)穩(wěn)定性，考慮了Cr5B3型、W5Si3型和Mn5Si3型3種晶體結(jié)構(gòu)，計(jì)算獲得的晶格常數(shù)、形成焓(表2)，同時分析3種晶體結(jié)構(gòu)的總能隨體積變化的關(guān)系(圖2)。本文計(jì)算獲得的Cr5B3型晶格常數(shù)為7.865 6 ?和15.366 3 ?，Bottcher等[1]的實(shí)驗(yàn)值為 7.963 ?和15.407 ?，其他第一性原理計(jì)算值[10]為 7.802 ?和15.620 ?。本文的晶格常數(shù)計(jì)算值與實(shí)驗(yàn)值之間的誤差為1.22%和0.26%，比Yang等[10]計(jì)算值的誤差(2.02%，1.38%)小，說明本文計(jì)算獲得的晶格常數(shù)更為合理。W5Si3型和Mn5Si3型Ca5Zn3化合物的晶格常數(shù)也作為預(yù)測結(jié)果列于表2。從表2中給出的形成焓可以看出，Cr5B3型Ca5Zn3的形成焓最低為-0.201 eV/atom，其他兩種晶體結(jié)構(gòu)的形成焓分為-0.160 eV/atom和-0.164 eV/atom。從能量的角度，其他兩種晶體結(jié)構(gòu)的形成焓比Cr5B3型的形成焓高出0.04 eV/atom，因此，具有Cr5B3結(jié)構(gòu)的Ca5Zn3是最穩(wěn)定的化合物。從圖2中可以看出，具有Cr5B3型的Ca5Zn3具有最低的總能量，同樣表明Cr5B3型Ca5Zn3是最穩(wěn)定的化合物，這與Bottcher等[1]的實(shí)驗(yàn)結(jié)果是一致的。

表2 Ca5Zn3化合物的晶格常數(shù)和形成焓Table 2 Lattice constant and formation enthalpy of the Ca5Zn3 compounds

圖2 具有3種不同晶體結(jié)構(gòu)的Ca5Zn3化合物總能量與體積的關(guān)系Fig.2 Relationship of total energy and volume of Ca5Zn3 compounds with three different crystal structures

2.2 力學(xué)性能

表3 Ca5Zn3化合物的力學(xué)性能參數(shù)Table 3 Mechanical property parameters of Ca5Zn3 compounds

B/G值一般用來預(yù)測化合物的脆韌性[17]，該比值大于1.75，化合物呈現(xiàn)韌性，比值越大韌性越大；反之，該比值小于1.75，則化合物呈現(xiàn)脆性，且比值越小脆性越大。Ca5Zn3的B/G值為2.61，說明該化合物呈現(xiàn)韌性。為了研究化合物的各向異性，對于四方晶系來說，一般用A1=2C66/(C11-C12)和A2=2C44/(C11-C12)進(jìn)行說明。該數(shù)值為1時，則呈現(xiàn)各向同性，偏離1越遠(yuǎn)則各向異性越明顯。為了更加明顯地呈現(xiàn)該化合物的各向異性，圖3給出了楊氏模量的空間分布圖，由此可以清晰地看出在不同的方向上，楊氏模量的數(shù)值存在明顯的差異。

圖3 Ca5Zn3化合物楊氏模量的空間分布(GPa)Fig.3 Spatial distribution of Young's modulus for Ca5Zn3 compounds (GPa)

2.3 電子結(jié)構(gòu)特性

為了深入了解該化合物原子間的相互作用，本文計(jì)算了該化合物的電子態(tài)密度和電荷密度差(圖4)。由圖4a可以看出，在費(fèi)米面以下的電子態(tài)密度，主要是Ca的s電子和p電子以及Zn的s電子和p電子的貢獻(xiàn)，費(fèi)米面以上的電子態(tài)密度主要是Ca的d電子和Zn的p電子的貢獻(xiàn)。在費(fèi)米面附近，Ca的p電子和d電子跟Zn的p電子發(fā)生雜化效應(yīng)。在費(fèi)米能外態(tài)密度不為零，說明該化合物呈現(xiàn)金屬性。由圖4b和圖4c可知，Ca和Zn之間主要呈現(xiàn)離子鍵特性，Ca失去電子，而Zn則獲得電子，這個和圖4a中的電子態(tài)密度在費(fèi)米面附近發(fā)生雜化效應(yīng)是一致的。由于Ca的電負(fù)性為1.0，Zn的電負(fù)性為1.6，電負(fù)性越大表明越容易獲得電子，由此可見本文的計(jì)算結(jié)果是合理的。

(a)電子態(tài)密度，(b) (001)面，和(c) (110)面(a) electronic density of state,(b) (001) plane and (c) (110) plane圖4 Ca5Zn3的電荷密度差Fig.4 Charge density difference of Ca5Zn3

2.4 熱物理性能

Ca5Zn3化合物的德拜溫度為202 K。平均聲速為2 097 m/s，縱向聲速為3 072 m/s， 橫向聲速為1 901 m/s (表3)。利用德拜模型[18-21]，結(jié)合第一性原理總能量的計(jì)算結(jié)果，計(jì)算了Ca5Zn3化合物在0—1 000 K溫度范圍、0—10 GPa壓強(qiáng)范圍內(nèi)的熱物理性質(zhì)(圖5)。由圖5a可知，化合物的體積隨著溫度的升高而增加，但是隨著壓強(qiáng)的增加則降低。并且隨著壓強(qiáng)的越來越大，體積減小的幅度越來越小，說明原子間的作用力也越來越大。圖5b給出了體積模量隨溫度和壓強(qiáng)的變化關(guān)系，隨著溫度的升高，體積模量有所下降，但是隨著壓強(qiáng)的增加，體積模量也在增加。圖5c表示熱膨脹系數(shù)隨著溫度的升高而增加，隨著壓強(qiáng)的增大而熱膨脹系數(shù)增加的幅度降低。在6 GPa以上時，隨著溫度的上升，熱膨脹系數(shù)增加的幅度極小。表明在高壓下，熱膨脹系數(shù)變化小，體積膨脹不明顯。圖5d展示的是定容熱容隨溫度和壓強(qiáng)的變化關(guān)系，在0—300 K區(qū)間，定容熱容急速增加，在300—400 K區(qū)間，等容熱容的上升趨勢變緩，在400 K以上，已經(jīng)接近等容熱容的極限值，即杜隆珀蒂極限。

圖5 Ca5Zn3的(a)體積，(b)體積模量，(c) 熱膨脹系數(shù)和(d)等容熱容隨溫度和壓強(qiáng)的變化Fig.5 Change of (a) volume,(b) bulk modulus,(c) thermal expansion and (d) heat capacity of Ca5Zn3 with the temperature at various pressures

3 結(jié)論

1)具有Cr5B3型的Ca5Zn3形成焓最低，是最穩(wěn)定的結(jié)構(gòu)，與前人的實(shí)驗(yàn)結(jié)果一致。

2)彈性常數(shù)的計(jì)算值顯示其滿足四方晶系的力學(xué)穩(wěn)定性標(biāo)準(zhǔn)，計(jì)算得到了體積模量、剪切模量、楊氏模量以及B/G值，B/G值顯示該化合物呈現(xiàn)韌性特征。楊氏模量呈現(xiàn)各向異性。

3)電子態(tài)密度顯示Ca的p電子和d電子與Zn的p電子在費(fèi)米能附近發(fā)生雜化效應(yīng)，費(fèi)米面上有態(tài)密度說明該化合物為金屬性。電荷密度差顯示Ca失去電子、Zn得到電子，與這兩種元素的電負(fù)性符合，同時呈現(xiàn)離子鍵特性。

4)體積、體積模量、熱膨脹系數(shù)以及定容熱容隨著溫度和壓強(qiáng)的變化趨勢說明，隨著壓強(qiáng)的增加、體積變小，體積模量變大、熱膨脹系數(shù)變小，定容熱容在低溫下變小，高溫下趨于一致。