張 暢，孫小偉，宋 婷，田俊紅，劉子江

（蘭州交通大學(xué)數(shù)理學(xué)院, 甘肅 蘭州 730070）

超硬材料是指維氏硬度大于40 GPa 的材料[1-2]。這些材料具有優(yōu)異的性能，如高硬度、良好的導(dǎo)熱性和不可壓縮性、穩(wěn)定的折射率和化學(xué)性質(zhì)，廣泛應(yīng)用于工業(yè)和日常生活中，如磨料涂層、切割和拋光工具等[3-5]。2005 年，Kaner 等[6]提出，過渡金屬具有較高的價(jià)電子密度和較強(qiáng)的不可壓縮性，輕元素之間的共價(jià)鍵為結(jié)構(gòu)提供了較高的剪切模量，過渡金屬硼化物、碳化物和氮化物具有較高的硬度，甚至可能達(dá)到超硬材料的標(biāo)準(zhǔn)。鉻硼化物是一種耐磨材料，人們從實(shí)驗(yàn)和理論上對(duì)其特性進(jìn)行了廣泛的研究[7-8]。鉻硼化物中，通過實(shí)驗(yàn)證實(shí)其結(jié)構(gòu)有CrB4、CrB2、Cr2B3、Cr3B4、CrB、Cr5B3和Cr2B[9-13]。其中，CrB4因其高硬度備受關(guān)注，例如：Niu 等[14]計(jì)算出空間群為Pnnm的CrB4的硬度為48 GPa，而后，Zhang 等[15]利用第一性原理方法系統(tǒng)地研究了CrB4的硬度隨壓力的變化規(guī)律。其余鉻硼化物普遍具有較低的硬度[16-20]，限制了其工業(yè)應(yīng)用。摻雜是一種提高材料硬度的行之有效的方法。Dovale-Farelo 等[21]采用第一性原理方法研究了不同濃度下Cr1-xMoxB2的機(jī)械性能，結(jié)果表明，CrB2的硬度為17 GPa，不滿足硬質(zhì)材料的標(biāo)準(zhǔn)，而Cr3MoB8、Cr2MoB6、CrMoB4、CrMo2B6和CrMo3B8的硬度均超過20 GPa，為硬質(zhì)材料。在鉻硼化物體系中，Cr2B3的力學(xué)性質(zhì)仍缺少系統(tǒng)的理論計(jì)算。1985 年，Okada 等[22]通過高溫Al-Cr-B 熔體制備了一種新型硼化鉻晶體，即Cr2B3，隨后采用單晶衍射法對(duì)Cr2B3的結(jié)構(gòu)進(jìn)行了研究[11]；1996 年，該團(tuán)隊(duì)發(fā)現(xiàn)Cr2B3結(jié)晶呈正交晶系，空間群為Cmcm，在(010)平面上測(cè)得的硬度為(22.4±1.7) GPa[17]。由此可見，Cr2B3的硬度僅達(dá)到硬質(zhì)材料的標(biāo)準(zhǔn)。通過元素?fù)诫s方式提高Cr2B3的硬度，將有助于拓展其工業(yè)應(yīng)用范圍。此外，在實(shí)際應(yīng)用中，往往會(huì)遇到高壓這一特殊工況。因此，研究摻雜Cr2B3在高壓下的結(jié)構(gòu)和力學(xué)性質(zhì)尤為重要。

本工作首先將Cr2B3中的部分Cr 元素分別替換為Sc、Ti、V 和Zr 元素，考慮兩種不同的替換方式，然后篩選出更優(yōu)的替換位置；而后，利用第一性原理方法，計(jì)算摻雜結(jié)構(gòu)的晶格常數(shù)、彈性常數(shù)、態(tài)密度、電子局域函數(shù)等；最后，基于電子結(jié)構(gòu)討論過渡金屬元素Sc、Ti、V、Zr 的摻雜對(duì)Cr2B3維氏硬度和彈性模量的影響。

1 計(jì)算方法

采用密度泛函理論框架下的贗勢(shì)平面波方法，運(yùn)用CASTEP 軟件包[23]對(duì)Cr2B3及其摻雜化合物的結(jié)構(gòu)進(jìn)行優(yōu)化。實(shí)驗(yàn)中證實(shí)Cr2B3為正交晶系結(jié)構(gòu)[22]，晶胞中共包含20 個(gè)原子，考慮50%的摻雜濃度，則其摻雜化合物的化學(xué)式可表示為CrMB3（M=Sc, Ti, V, Zr）。電子與電子之間的交換關(guān)聯(lián)能采用Perdew 等[24]修正的PBESOL 形式的廣義梯度近似（generalized gradient approximation，GGA），離子實(shí)與價(jià)電子之間的相互作用勢(shì)用超軟贗勢(shì)描述[25]。Cr 原子和B 原子的電子組態(tài)分別為 3d54s1和2s22p1，摻雜原子 Sc、Ti、V 和 Zr 的電子組態(tài)分別為3d14s2、3d24s2、3d34s2和4d25s2。倒易空間布里淵區(qū)k點(diǎn)采用Monkhorst-Pack[26]方法選取，積分網(wǎng)格為14×2×14，平面波基函數(shù)的截?cái)嗄芰咳?50 eV。迭代過程中系統(tǒng)能量的收斂標(biāo)準(zhǔn)為5×10-6eV/atom，作用在晶胞中每個(gè)原子上的力小于0.01 eV/?，晶胞應(yīng)力偏差低于0.02 GPa，最大位移的收斂公差小于5×10-4?。

2 結(jié)果與討論

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

Cr2B3晶胞包含4 個(gè)Cr2B3單元，其空間群為Cmcm，晶格常數(shù)a= 3.026 4 ?，b= 18.115 0 ?，c= 2.954 2 ?。該結(jié)構(gòu)共包含8 個(gè)Cr 原子、12 個(gè)B 原子，其中：Cr 原子有2 種占據(jù)方式，Cr1 和Cr2 占據(jù)的Wyckoff 位置分別為4c (0, 0.427 66, 0.25)和4c (0, -0.296 30, 0.25)；B 原子有3 種占據(jù)方式，即4c (0, 0.023 60, 0.25)、4c (0, 0.118 60, 0.25)和4c (0, -0.171 30, 0.25)，晶體結(jié)構(gòu)如圖1(a)所示。將該結(jié)構(gòu)作為母體結(jié)構(gòu)，對(duì)Cr2B3進(jìn)行摻雜。用 Sc、Ti、V 和Zr 分別替換Cr2B3晶格中的Cr1 和Cr2，并對(duì)其進(jìn)行幾何優(yōu)化，計(jì)算結(jié)果如表1 所示。除Sc 替換Cr2 位置外，其他替換情況下計(jì)算結(jié)果均可達(dá)收斂標(biāo)準(zhǔn)。Ti、V、Zr 摻雜元素替換Cr1 位置時(shí)具有更低的形成能和生成焓，由此可見，Cr1 位置為更優(yōu)的替換位置。因此，將Cr2B3晶格中的Cr1 原子分別替換為過渡金屬原子 Sc、Ti、V 和Zr，獲得了摻雜結(jié)構(gòu)CrMB3（M = Sc, Ti, V, Zr），如圖1(b)所示。表1 也給出了優(yōu)化后的晶格常數(shù)，與實(shí)驗(yàn)值[11]比較發(fā)現(xiàn)，二者之間的差異很小，表明本研究選取的計(jì)算參數(shù)是可靠的。同時(shí)，摻雜結(jié)構(gòu)的空間群未發(fā)生變化。由此可知，原子的替換并未對(duì)Cr2B3的晶格結(jié)構(gòu)產(chǎn)生明顯的影響。

圖1 (a) Cr2B3 和 (b) CrMB3（M = Sc, Ti, V, Zr）的晶體結(jié)構(gòu)（藍(lán)色、綠色和粉色小球分別代表 Cr、B 和過渡金屬原子）Fig. 1 Crystal structures of (a) Cr2B3 and (b) CrMB3 (M=Sc, Ti, V, Zr), where the blue, green and pink spheres represent Cr, B and transition metal atoms, respectively

為研究摻雜元素對(duì)Cr2B3彈性性質(zhì)的影響，首先要判斷摻雜結(jié)構(gòu)的穩(wěn)定性。固體的結(jié)構(gòu)穩(wěn)定性可通過熱力學(xué)穩(wěn)定性、力學(xué)穩(wěn)定性和動(dòng)力學(xué)穩(wěn)定性判斷。

對(duì)于母體結(jié)構(gòu)Cr2B3，其熱力學(xué)穩(wěn)定性可通過生成焓判斷

式中：ΔH(Cr2B3) 為化合物 Cr2B3的生成焓，H(Cr2B3)為Cr2B3的焓值，H(Cr)和H(B)分別為Cr 和B 作為單質(zhì)存在的焓值。對(duì)于摻雜結(jié)構(gòu)，則需要計(jì)算其形成能Ef[27]

式中：Ef為化合物 Cr2B3的形成能，E(CrMB3)為Cr2B3中摻雜Sc、Ti、V 和Zr 原子的總能量，E(Cr2B3)為Cr2B3的能量值，n為摻雜的原子數(shù)，μ為相應(yīng)原子的化學(xué)勢(shì)。 ΔH和Ef的計(jì)算結(jié)果也列于表1 中?？梢钥闯?，在零壓下，Cr2B3的ΔH為負(fù)值，通過替換 Cr1 位置得到的摻雜結(jié)構(gòu)CrMB3（M = Sc, Ti, V, Zr）的ΔH和Ef均為負(fù)值，表明它們都具備熱力學(xué)穩(wěn)定性。

表1 零壓下Cr2B3 及摻雜結(jié)構(gòu)CrMB3（M=Sc, Ti, V, Zr）的晶格常數(shù)、形成焓及摻雜結(jié)構(gòu)的形成能Table 1 Lattice constants, formation enthalpy of Cr2B3 and CrMB3 (M = Sc, Ti, V, Zr), and impurity formation energy of CrMB3 (M=Sc, Ti, V, Zr) at zero pressure

本研究將通過計(jì)算上述結(jié)構(gòu)的彈性常數(shù)來判斷它們的力學(xué)穩(wěn)定性。上述所有結(jié)構(gòu)都屬于正交晶系，該晶系具有9 個(gè)獨(dú)立的彈性常數(shù)，即C11、C22、C33、C44、C55、C66、C12、C13和C23。根據(jù)Born 穩(wěn)定性準(zhǔn)則[28]，有

圖2 和圖3 給出了根據(jù)應(yīng)力-應(yīng)變關(guān)系計(jì)算得到的Cr2B3及其摻雜結(jié)構(gòu)在不同壓力下的彈性常數(shù)。上述所有結(jié)構(gòu)在 0～150 GPa 的壓力范圍內(nèi)均滿足正交晶系的力學(xué)穩(wěn)定性條件。將摻雜結(jié)構(gòu)的彈性常數(shù)與 Cr2B3比較發(fā)現(xiàn)，V 摻雜可使Cr2B3的彈性常數(shù)C11、C22和C33均增大。C11、C22和C33是表征材料在a、b和c軸抗形變能力的參數(shù)，其值越大，抗形變能力越強(qiáng)。圖2 顯示，隨著壓力的升高，Cr2B3晶體在a、b和c軸上的抗形變能力逐漸增強(qiáng)。C44作為剪切彈性常數(shù)，其值越大，材料的剪切變形抗力越強(qiáng)。由圖3 可知，CrTiB3和CrVB3的C44高于其他幾種結(jié)構(gòu)，這兩種結(jié)構(gòu)對(duì)應(yīng)的剪切模量也高于其他幾種結(jié)構(gòu)，因此，CrTiB3和CrVB3具有更強(qiáng)的剪切變形抗力。

圖3 CrMB3（M=Sc, Ti, V, Zr）的彈性常數(shù)隨壓力的變化Fig. 3 Pressure dependence of the elastic constants for CrMB3 (M=Sc, Ti, V, Zr)

為判斷Cr2B3的動(dòng)力學(xué)穩(wěn)定性，通過有限位移法分別計(jì)算了Cr2B3在零壓和150 GPa 高壓下的聲子色散曲線，如圖4 所示。零壓下，Cr2B3晶體的聲子色散關(guān)系在整個(gè)布里淵區(qū)內(nèi)未出現(xiàn)虛頻，表明Cr2B3在零壓下滿足動(dòng)力學(xué)穩(wěn)定性；相比于零壓，Cr2B3在150 GPa 下的聲子色散曲線的頻率有一定的提高，但其在X附近表現(xiàn)出虛頻，因而動(dòng)力學(xué)不穩(wěn)定。圖5 為CrScB3、CrTiB3、CrVB3和 CrZrB3在零壓下的聲子色散曲線，其中 CrVB3在G點(diǎn)出現(xiàn)微小虛頻，CrScB3、CrTiB3和 CrZrB3出現(xiàn)了明顯虛頻，均不滿足動(dòng)力學(xué)穩(wěn)定性。

圖4 Cr2B3 在零壓(a)和 150 GPa (b)下的聲子色散曲線Fig. 4 Phonon-dispersion curves of Cr2B3 at (a) 0 GPa and (b) 150 GPa

2.2 彈性性質(zhì)

根據(jù)彈性常數(shù)，可獲得材料的體積彈性模量、剪切模量以及楊氏模量等力學(xué)參量。體積彈性模量B和剪切模量G采用Voigt-Reuss-Hill 近似法計(jì)算[29]， 楊氏模量E可由B和G推導(dǎo)得出，計(jì)算公式如下

式中：BV和BR為Voigt 和Reuss 模型的體積彈性模量，GV和GR分別為Voigt 剪切模量和Reuss 剪切模量。圖6 為Cr2B3和CrMB3（M=Sc, Ti, V, Zr）的體積彈性模量、剪切模量和楊氏模量隨壓力的變化?？梢钥闯觯S著壓力增加，Cr2B3及其摻雜結(jié)構(gòu)的體積彈性模量、剪切模量、楊氏模量均增大。Cr2B3在整個(gè)壓力范圍內(nèi)的體積彈性模量比其他4 種結(jié)構(gòu)更大。相比于零壓下的Cr2B3，零壓下CrVB3和CrTiB3的剪切模量分別提高14.3%和16.2%，楊氏模量分別提高8.2%和12.0%。

圖6 Cr2B3 和 CrMB3（M=Sc, Ti, V, Zr）的體積彈性模量(a)、剪切模量(b)和楊氏模量(c)隨壓力的變化Fig. 6 Pressure dependence of (a) the bulk moduli, (b) the shear moduli, (c) the Young’s moduli for Cr2B3 and CrMB3 (M=Sc, Ti, V, Zr)

硬度是材料抵抗永久變形能力的量度，是材料的重要力學(xué)參數(shù)。Cr2B3及其摻雜結(jié)構(gòu)的理論維氏硬度采用半經(jīng)驗(yàn)?zāi)Ｐ蚚30]計(jì)算

圖7 為 Cr2B3及其摻雜結(jié)構(gòu)的維氏硬度隨壓力的變化曲線。可以看出，在研究的整個(gè)壓力范圍內(nèi)，除壓力超120 GPa 的Zr 摻雜外，Sc、Ti、V 和Zr 元素的摻雜均使Cr2B3的維氏硬度有所提高。值得注意的是，Ti 摻雜使Cr2B3的硬度從26.3 GPa 提高至40.2 GPa，提高52.9%，達(dá)到超硬材料的標(biāo)準(zhǔn)。隨著壓力的增加，Cr2B3的硬度呈線性增加，但其硬度仍較低，150 GPa 下僅為28.3 GPa，CrZrB3的硬度總體呈下降趨勢(shì)， CrVB3的硬度在37 GPa 上下小幅波動(dòng)，而CrScB3和CrTiB3的硬度隨著壓力的增加發(fā)生明顯的波動(dòng)。在CrScB3的硬度隨壓力的變化曲線中，100 GPa 時(shí)出現(xiàn)明顯的波動(dòng)，計(jì)算得到該點(diǎn)的形成能Ef為0.083 295 eV，為正值，不滿足熱力學(xué)穩(wěn)定性，因此出現(xiàn)了力學(xué)性質(zhì)異常。CrTiB3在零壓下的維氏硬度明顯高于其他壓力下的硬度，達(dá)到超硬材料的標(biāo)準(zhǔn)，且在整個(gè)壓力范圍內(nèi)，CrTiB3的硬度始終高于另外4 種結(jié)構(gòu)。總體看來，CrVB3同時(shí)具備較高的硬度以及較小的壓力依賴性，是理想的硬質(zhì)材料。

Cr2B3的硬度隨壓力變化呈現(xiàn)出明顯的規(guī)律性，因此針對(duì)Cr2B3在不同壓力下的B-B 鍵的鍵長及布居數(shù)做進(jìn)一步計(jì)算，結(jié)果如表2 所示。通常，鍵長越短，材料的硬度越高。表2 顯示，隨著壓力的升高，Cr2B3的3 類B-B 鍵，即B1-B1、B1-B2 和B2-B3 的鍵長均縮短。此外，還可通過布居數(shù)判斷B-B 共價(jià)鍵強(qiáng)弱（通常，正值表示共價(jià)鍵，正值越大，共價(jià)性越強(qiáng)），進(jìn)而分析硬度隨壓力變化的原因。Cr2B3的B-B 鍵的布居數(shù)隨壓力升高均增大，表明B-B 鍵的共價(jià)性隨壓力的升高而增強(qiáng)，與Cr2B3的硬度隨壓力的變化趨勢(shì)一致。

表2 Cr2B3 中B-B 鍵的鍵長以及布居數(shù)隨壓力的變化Table 2 Pressure dependence of B-B bond length and population for Cr2B3

過渡金屬元素的摻雜雖提高了Cr2B3的硬度，但同時(shí)使材料的韌性有所降低。根據(jù)Pugh 規(guī)則[31]，可以判斷出固體的韌脆性。一般說來，當(dāng)B/G小于1.75 時(shí)，材料表現(xiàn)出脆性，反之則表現(xiàn)出韌性。圖8給出了Cr2B3和CrMB3（M=Sc, Ti, V, Zr）的B/G隨壓力的變化。零壓下，Cr2B3及其摻雜結(jié)構(gòu)的B/G均小于1.75，表現(xiàn)出脆性，而摻雜結(jié)構(gòu)的B/G比Cr2B3的更小，表現(xiàn)出更明顯的脆性；Cr2B3和CrZrB3分別在65.0 和93.7 GPa 由脆性轉(zhuǎn)變?yōu)轫g性。

圖8 Cr2B3 及 CrMB3（M=Sc, Ti, V, Zr）的B/G 隨壓力的變化Fig. 8 Pressure dependence of the B/G for Cr2B3 and CrMB3(M=Sc, Ti, V, Zr)

圖9、圖10、圖11 和圖12 為零壓下Cr2B3及其摻雜結(jié)構(gòu)的體積彈性模量和楊氏模量在xy、xz和yz面的投影圖。在yz平面上，Cr2B3及其摻雜結(jié)構(gòu)的體積彈性模量的二維投影近似為圓形，表明材料的體積彈性模量在yz平面上呈各向同性；在xy平面和xz平面上，Cr2B3及其摻雜結(jié)構(gòu)的體積彈性模量和楊氏模量的二維投影偏離圓形，表明材料的體積彈性模量和楊氏模量在xy和xz平面上呈各向異性。摻雜結(jié)構(gòu)與Cr2B3的體積彈性模量大致相同，說明摻雜對(duì)體積彈性模量的各向異性沒有明顯的影響。CrTiB3和CrVB3的楊氏模量在xz平面上的投影相比于Cr2B3僅出現(xiàn)微小變化，說明在xz平面上V 和Ti 元素的摻雜對(duì)楊氏模量各向異性的影響較小。

圖9 零壓下 Cr2B3 的體積彈性模量的三維表示(a)及其在 xy、xz 和 yz 面上的投影(b)Fig. 9 (a) 3D representation and (b) 2D projections on xy, xz and yz planes of the bulk modulus for Cr2B3 at 0 GPa

圖10 零壓下 Cr2B3 的楊氏模量的三維表示(a)及其在 xy、xz 和 yz 面上的投影(b)Fig. 10 (a) 3D representation and (b) 2D projections on xy, xz and yz planes of the Young’s modulus for Cr2B3 at 0 GPa

圖11 零壓下 CrMB3（M=Sc, Ti, V, Zr）的體積彈性模量在 xy、xz 和 yz 面的投影Fig. 11 2D projections of the bulk modulus for CrMB3 (M=Sc, Ti, V, Zr) on xy, xz and yz planes at 0 GPa

圖12 零壓下 CrMB3（M=Sc, Ti, V, Zr）的楊氏模量在 xy，xz 和 yz 面的投影Fig. 12 2D projections of the Young’s modulus for CrMB3 (M=Sc, Ti, V, Zr) on xy, xz and yz planes at 0 GPa

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

電荷密度和態(tài)密度（density of states，DOS）可以用于揭示結(jié)構(gòu)的電子性質(zhì)。計(jì)算了零壓下 Cr2B3及其摻雜結(jié)構(gòu)的總態(tài)密度和分波態(tài)密度，如圖13 和圖14 所示，其中虛線所處位置為費(fèi)米能級(jí)。在費(fèi)米能級(jí)上，各結(jié)構(gòu)的態(tài)密度值均大于零，表明所有結(jié)構(gòu)都具有明顯的金屬性。通過態(tài)密度分析可以看出：對(duì)于Cr2B3，在-5.0～-3.4 eV 區(qū)間，B-2p 和Cr-3d 具有相似形狀的峰，出現(xiàn)雜化現(xiàn)象，說明Cr 原子與B 原子之間的化學(xué)鍵存在共價(jià)鍵成分，而在8.9～12.3 eV 區(qū)間，B-2p 與B-2s 發(fā)生了強(qiáng)烈的雜化，使得Cr2B3具有較高的硬度；對(duì)于CrScB3，在2.3～3.2 eV 區(qū)間，B-2p 與Cr-3d 發(fā)生雜化；對(duì)于CrTiB3，在-3.8～-2.5 eV 區(qū)間，B-2p 與Ti-3d 發(fā)生雜化；對(duì)于CrVB3，在-4.4～-3.0 eV 區(qū)間，V-3d 與B-2p 發(fā)生雜化；對(duì)于CrZrB3，在-3.4～-2.3 eV和1.5～3.0 eV 區(qū)間，B-2p 與Cr-3d 以及Zr-4d 出現(xiàn)了形狀相似的峰，說明此處出現(xiàn)雜化現(xiàn)象。

圖13 零壓下 Cr2B3 的總態(tài)密度以及分波態(tài)密度Fig. 13 Total density of states and partial density of states for Cr2B3 at 0 GPa

圖14 零壓下 CrMB3（M=Sc, Ti, V, Zr）的總態(tài)密度以及分波態(tài)密度Fig. 14 Total density of states and partial density of states for CrMB3 (M=Sc, Ti, V, Zr) at 0 GPa

金屬與非金屬之間的B―M 共價(jià)鍵有利于提高剪切模量，形成的M―M 金屬鍵將降低材料的剪切模量。為了揭示化學(xué)鍵合的性質(zhì)與維氏硬度的關(guān)系，計(jì)算了Cr2B3和CrMB3（M=Sc, Ti, V, Zr）的電子局域函數(shù)（electronic local functions，ELF）[32]，并選取了(100)平面和( 0 01)平面進(jìn)行分析，結(jié)果如圖15 所示。圖15 中電荷高度局域位置用紅色顯示，表明該處存在強(qiáng)的共價(jià)鍵。對(duì)比Cr2B3及其摻雜結(jié)構(gòu)在(100)面上的電荷局域分布情況，可以看出摻雜結(jié)構(gòu)B-B 之間的紅色區(qū)域相較于Cr2B3均有所擴(kuò)大。在( 0 01)晶面上，Cr2B3的電荷主要集中在B-B 之間的連線上，而摻雜結(jié)構(gòu)的B 原子被紅色區(qū)域包圍，代表著電子集中在B 原子的周圍。材料的硬度往往與B-B 共價(jià)鍵的強(qiáng)弱相關(guān)。摻雜結(jié)構(gòu)具有更強(qiáng)的B-B共價(jià)鍵，因而擁有更高的硬度。

圖15 零壓下(a) Cr2B3、(b) CrScB3、(c) CrTiB3、(d) CrVB3、(e) CrZrB3 在(100)平面的 ELF 以及零壓下(f) Cr2B3、(g) CrScB3、(h) CrTiB3、(i) CrVB3 和(j) CrZrB3 在( 001)平面的 ELFFig. 15 Electronic local functions contours for (a) Cr2B3, (b) CrScB3, (c) CrTiB3, (d) CrVB3 and (e) CrZrB3 in plane(100) at 0 GPa, and electronic local functions contours for (f) Cr2B3, (g) CrScB3, (h) CrTiB3,(i) CrVB3 and (j) CrZrB3 in plane ( 001) at 0 GPa

3 結(jié) 論

利用第一性原理方法，對(duì)Cr2B3的生成焓，CrMB3（M=Sc, Ti, V, Zr）的形成能，以及他們的彈性常數(shù)、彈性模量和電子結(jié)構(gòu)進(jìn)行了詳細(xì)的計(jì)算和討論。結(jié)果表明，在零壓下 Sc、Ti、V 和 Zr 元素的摻雜均使Cr2B3的維氏硬度有所提高。其中，CrTiB3的硬度達(dá)到40.2 GPa，為潛在的超硬材料。在 0～150 GPa壓力范圍內(nèi)，Cr2B3、CrVB3和CrZrB3的硬度隨壓力的變化趨勢(shì)呈現(xiàn)一定的規(guī)律性：Cr2B3的硬度隨壓力的增加而單調(diào)增加；CrVB3的硬度受壓力的影響很小，且相較于Cr2B3具有更高的硬度；CrZrB3的硬度隨壓力的增加，總體呈下降趨勢(shì)。Cr2B3的B-B 鍵長隨壓力增大而縮短，同時(shí)布居數(shù)增大，表明Cr2B3的B-B 共價(jià)鍵增強(qiáng)，導(dǎo)致Cr2B3的硬度隨壓力的增加而增大。