彭紅建 ，武青，周姣連，李小波，謝佑卿

(1. 中南大學 化學化工學院，湖南 長沙，410083；2. 中南大學 材料科學與工程學院，湖南 長沙，410083)

價鍵結(jié)構(gòu)是指固體或分子中原子的價電子結(jié)構(gòu)和原子與原子之間形成的鍵以及兩者之間的關(guān)系。研究材料的價鍵結(jié)構(gòu)有助于從原子結(jié)構(gòu)層次和電子結(jié)構(gòu)層次更深入地了解材料的微觀結(jié)構(gòu)，設計出滿足人們需要的新材料，也是當今材料設計最為活躍的前沿領(lǐng)域之一[1-6]。確定材料的價鍵結(jié)構(gòu)主要有2種：第一原理方法和鍵價法。鍵價法是在 Pauling電價規(guī)則基礎上發(fā)展起來的，采用與量子力學不同的規(guī)則，通過鍵中求價來確定材料的鍵價關(guān)系。本文作者根據(jù)金屬材料系統(tǒng)科學框架[7-11]，從“能”和“形”2個方面確定價電子結(jié)構(gòu)。能形法[12]是在單原子狀態(tài)自洽法的基礎上將三態(tài)雜化擴展為多態(tài)雜化，原子雜化狀態(tài)可連續(xù)地變化，其中能量最低的原子狀態(tài)即為最穩(wěn)定結(jié)構(gòu)的原子狀態(tài)，克服兩態(tài)雜化和三態(tài)雜化多解的問題及基本態(tài)選擇的困難，定出更加精細的原子狀態(tài)。根據(jù)能形法確定金屬Mo的價鍵結(jié)構(gòu)，與第一原理方法和鍵距差等方法進行比較，并研究其隨溫度和壓力的變化。

1 能形法的原理與方法

金屬的價鍵結(jié)構(gòu)包括價電子結(jié)構(gòu)、鍵結(jié)構(gòu)和單鍵半徑如圖1所示。

圖1 BCC體系價鍵結(jié)構(gòu)的示意圖Fig.1 Valence bond structure of BCC-centered cubic crystal

根據(jù)晶格常數(shù)(式(1))和結(jié)合能(式(2))便可計算金屬的晶格常數(shù)和結(jié)合能。

式中：ec為雜化原子的共價電子數(shù)，sc為 s軌道共價電子數(shù)；pc為 p軌道共價電子數(shù)；ef為雜化原子的自由電子數(shù)；a為晶格常數(shù)；s代表不同鍵的類型(1代表最近鄰鍵，2代表第二近鄰鍵，3代表第三近鄰鍵，…)；ec和R分別為雜化原子的共價電子數(shù)和單鍵半徑；rs，Is和ns分別為各類鍵的鍵長、等同鍵數(shù)和鍵價；Gs和Is為與晶格類型相關(guān)的常數(shù)；β按EET理論取值。對于fcc和bcc晶體結(jié)構(gòu)，僅考慮原子的近鄰、第二近鄰和第三近鄰的相互作用。

式中：A表示原子中電子電荷對核電荷的屏蔽常數(shù)，由于原子內(nèi)層電子的屏蔽、電子的庫侖、交換及關(guān)聯(lián)作用對鍵能影響的總效應，它是一個與元素相關(guān)的量，其值通過晶體結(jié)合能的實驗值來定；z是原子外層電子總數(shù)；f 和 f′分別表示共價電子和近自由電子的成鍵能力，它與價電子的關(guān)系為：

任意溫度T下的平衡晶格常數(shù)為：

其中：a293為293 K時的晶格常數(shù)；α為線熱膨脹系數(shù)。

對于隨壓力的變化采用 Rose物態(tài)方程[13]，其函數(shù)表達式采用如下形式：

其中：

B0，V0和E0分別是零溫零壓下體彈性模量、原子體積和結(jié)合能，從式(5)可知壓力與晶格常數(shù)的關(guān)系。

2 金屬Mo的價鍵結(jié)構(gòu)

從金屬Mo的外層電子結(jié)構(gòu)來看，約4.910 4個(平均數(shù))共價電子用來形成具有定域特征的共價鍵，1.089 6個(平均數(shù))近自由電子用來形成具有離域特征金屬鍵，這與金屬導電性和硬度極其相關(guān)；由于原子軌道的d-s-p雜化，約1.650 0個(平均數(shù))的d或者s軌道電子升級到p軌道以形成更強的共價鍵。

從能量分配來看，金屬的結(jié)合能主要來源于近鄰和次近鄰原子的相互作用所形成的共價能，共價電子形成的共價能為634.366 3 kJ/mol，但近自由電子所形成的自由能為23.972 4 kJ/mol，也是一種不可忽視的能量。對于金屬Mo外層價電子數(shù)有很多，但只有一種的能量是最低的，相對結(jié)合能隨sc和pc電子占居數(shù)的變化如圖2所示，因此，可以確定最穩(wěn)定狀態(tài)的sc和pc電子占居數(shù)。

對于Mo的價態(tài)，Pauling將它的價定為6；Raju[15]則認為是ds雜化，將其價態(tài)定為4d55s1，與自由原子相同；能形法將 Mo定為固定價 4.910 4，流動價1.089 6，其總價等于6.000 0；Mulliken方法計算金屬Mo的原子價軌占居數(shù)，外層電子的值為鍵或離域成鍵，又不能計算，因為第一原理是按能量的高低填充電子，因此，可以看出能形法能更精細地確定價電子結(jié)構(gòu)。

圖2 金屬Mo的相對結(jié)合能隨sc和pc電子占居數(shù)的變化Fig.2 Relative cohesive energy of Mo with sc and pc electrons

另外，鍵距差[15](BLD)法確定金屬 Mo價電子結(jié)構(gòu)為啞對電子數(shù)nd(0.761 4)，共價電子數(shù)nc(4)，晶格電子數(shù)nl(1.238 6)，總價電子數(shù)nT(5.238 6)，單鍵半徑R(0.123 78 nm)。BLD法是元素2個極限價態(tài)之間的某一離散價態(tài)，該價態(tài)的晶格常數(shù)與實驗晶格常數(shù)之差小于等于0.005 nm。由能形法所得無素價態(tài)是元素所有可能的價態(tài)中，與實驗晶格常數(shù)和結(jié)合能相差最小的價態(tài)。相比BLD法，能形法增加了“能”的判據(jù)。

對于單鍵半徑，Pauling的單鍵半徑為0.129 6 nm，計算值為0.128 08 nm。比較小的原因主要有2個：一是 Pauling只考慮了原子的近鄰相互作用，能形法不僅考慮近鄰原子的相互作用，而且考慮第二近鄰、第三近鄰原子的相互作用等；二是近自由電子對共價鍵的鍵長不起作用，由于它的存在導致共價電子數(shù)減少。

3 金屬Mo的價鍵結(jié)構(gòu)隨溫度的變化

金屬Mo的價鍵結(jié)構(gòu)隨溫度變化關(guān)系如圖3所示。其中：金屬 Mo價電子結(jié)構(gòu)(dc，sc，pc，ec和 ef)隨溫度變化如圖3(a)所示；單鍵半徑(R)和晶格常數(shù)(r2)隨溫度變化如圖3(b)和3(c)所示；最近鄰鍵長(r1)及鍵價(n1)隨溫度變化如圖3(d)和3(e)所示；結(jié)合能(Ecoh)、相對結(jié)合能(Ecoh/A)及屏蔽常數(shù)A隨溫度變化如圖3(f)，(g)和(h)所示。

圖3 金屬Mo的價鍵結(jié)構(gòu)與溫度的關(guān)系Fig.3 Relationship between valence bond structure of Mo metal and temperature

從圖 3(a)可以看出：隨著溫度升高，共價電子ec(sc，pc，dc) 減少，近自由電子ef增多，但在固相范圍內(nèi)變化不大；圖3(b)和(c)說明了鍵長與溫度的關(guān)系，鍵長隨溫度的升高而增大。圖3(d)和(e)說明了最近鄰鍵長與其鍵價的相互關(guān)系，可見：隨著溫度升高，鍵長增長，鍵價降低；圖3(f)，(g)和(h)說明了鍵能與溫度的關(guān)系，發(fā)現(xiàn)隨著溫度升高，鍵長變長，鍵價減少，即鍵上電子對數(shù)減少，鍵能也降低。這是成鍵電子減少、原子間距變大造成的。

4 金屬Mo的價鍵結(jié)構(gòu)隨壓力的變化

圖4 金屬Mo的價鍵結(jié)構(gòu)與壓力的關(guān)系Fig.4 Relationship between alence bond structure of Mo metal and pressure

根據(jù)式(5)和(7)及Rose在文獻[14]中確定3個參數(shù)B0，V0和E0，可以確定金屬Mo的價鍵結(jié)構(gòu)隨壓力變化如圖4所示。其中：金屬Mo價電子結(jié)構(gòu)(dc, sc, pc, ec和ef)隨壓力變化如圖4(a)所示；單鍵半徑(R)，晶格常數(shù)(r2)隨壓力變化如圖4(b)和4(c)所示；最近鄰鍵長(r1)及鍵價(n1)隨壓力變化如圖 4(d)和 4(e)所示；結(jié)合能(Ecoh)，相對結(jié)合能(Ecoh/A)及屏蔽常數(shù)A隨壓力變化如圖 4(f)，(g)和(h)所示。

從圖 4(a)可以看出：隨著壓力的增大，金屬 Mo的共價電子數(shù)ec(sc，pc，dc)增多，近自由電子數(shù)ef減少；在壓力0～130 GPa范圍內(nèi)，壓力對價電子結(jié)構(gòu)影響顯著，其表現(xiàn)金屬性下降，共價性增加。圖(b)和(c)說明了鍵長與壓力的關(guān)系，單鍵半徑隨壓力的增大而減小。圖4(d)和(e)說明了最近鄰鍵長與其鍵價的相互關(guān)系，可見鍵長變短，鍵價增加，即鍵上電子對數(shù)增加。圖4(f)，(g)和(h)說明了鍵能與壓力的關(guān)系，可見壓力增大，結(jié)合能降低。這是外界做功轉(zhuǎn)化為內(nèi)能的緣故。在整個壓力范圍內(nèi)，結(jié)合能變化還是很大，相對結(jié)合能與壓力呈單調(diào)遞增關(guān)系，這是成鍵電子增加、原子間距變小造成的。

通過結(jié)合能和相對結(jié)合能可以定出屏蔽常數(shù)，屏蔽常數(shù)與壓力呈單調(diào)遞減關(guān)系，因此，溫度和壓力對結(jié)合能的影響機理是不一樣的，結(jié)合能的變化是相對結(jié)合能和屏蔽常數(shù)相互作用的結(jié)果。當溫度變化時，相對結(jié)合能的減少量超過屏蔽常數(shù)的增大量，導致結(jié)合能降低；而壓力變化時，相對結(jié)合能的增加量小于屏蔽常數(shù)的減少量，導致結(jié)合能降低。

5 結(jié)論

(2) 隨著溫度升高，共價電子減少，近自由電子增加，但在固相范圍內(nèi)變化不大；隨著溫度升高，鍵長增長，鍵價降低，鍵能也降低。

(4) 隨著壓力的升高，共價電子增加，近自由電子減少，晶體的共價性增加，金屬性下降，且鍵長縮短，鍵價增大，鍵能降低。

[1] Brown I D. The chemical bond in inorganic chemistry: the valence bond model[M]. New York: Oxford University Press,2002.

[2] Brown I D. Influence of chemical and spatial constraints on the structures of inorganic compounds[J]. Acta Cryst. B, 1997, 53:381-393.

[3] Hunter B A, Howard C J, Kim D J. Bond valence analysis of tetragonal zirconias[J]. J Solid State Chem, 1999, 146(2):363-368.

[4] Swenson J, Adams S. Application of the bond valence method to reverse Monte Carlo produced structural models of superionic glasses[J]. Phys Rev B, 2001, 64(2): 024204(1-10).

[5] Adams S, Swenson J. Bond valence analysis of reverse Monte Carlo produced structural models; a way to understand ion conduction in glasses[J]. J Phys: Condens Matter, 2005, 17(5):S87-S101.

[6] Urusov V S. Semi-empirical groundwork of the bond-valence model[J]. Acta Cryst B, 1995, 51: 641-649.

[7] 陶輝錦, 謝佑卿, 彭紅建, 等. 面心立方和亞穩(wěn)液相金屬 Cu的電子結(jié)構(gòu)和物理性質(zhì)[J]. 中南大學學報: 自然科學版, 2007,38(1): 1-8.TAO Hui-jin, XIE You-qing, PENG Hong-jian, et al. Electronic structures and physical properties of FCC and metastable liquid Cu metals[J]. Journal of Central South University: Science and Technology, 2007, 38(1): 1-8.

[8] PENG Hong-jian, XIE You-qing, NIE Yao-zhuang,Thermodynamic properties of stable and metastable phases of Pt metal[J]. Tans Nonferrous Met Soc China, 2009, 19(2):428-432.

[9] XIE You-qing. The electronic structure and properties of pure iron[J]. Acta Metall Mater, 1994, 42(11): 3705-3715.

[10] 彭紅建, 謝佑卿, 李小波. 單原子狀態(tài)自洽法的算法[J]. 中南大學學報: 自然科學版, 2007, 38(5): 900-905.PENG Hong-jian, XIE You-qing, Li Xiao-bo. Algorithms of single-atom self-consistency method[J]. Journal of Central South University. Science and Technology, 2007, 38(5): 900-905.

[11] LI Xiao-bo, XIE You-qing, NIE Yao-zhuang, et al. Energy and shape method for determining the valence and bond structure of crystals[J]. Physica B, 2007, 391(2): 249-255.

[12] Rose J H, Smith J R, Guinea F, et al. Universal features of the equation of state of metals[J]. Phys Rev B, 1984, 29(6):2963-2969.

[13] Gray D E. American Institute of Physics, American Institute of Physics Handbook[M]. 3rd ed. New York: McGraw-Hill Book Company, 1972: 23-24.

[14] Raju S, Mohandas E, Raghunathan V S. Engel-brewer electron correlation model: A critical discussion and revision of concepts[J]. Mater Trans JIM, 1996, 37(3): 195-202.

[15] 張瑞林. 固體與分子經(jīng)驗電子理論[M]. 長春: 吉林科學技術(shù)出版社, 1993: 255-257.ZHANG Rui-lin. The empirical electronic theory of solids and molecules[M]. Changchun: Jilin Science and Technology Press,1993: 255-257.