趙定國，郭培民，趙沛

(鋼鐵研究總院 先進鋼鐵流程及材料國家重點實驗室，北京 100081)

化合物生成焓的測定是冶金和材料領域長期需要開展的基礎研究。測定化合物的熱力學標準生成焓的方法很多，同一種化合物的實驗結果也常常有差異，且實驗方法的可靠程度也有待研究[1]。現(xiàn)有的無機化合物中，大部分物相的物性是未知的或不確定的，各種熱力學數(shù)據(jù)手冊及電子數(shù)據(jù)庫也僅針對已知數(shù)據(jù)，更多的熱力學數(shù)據(jù)卻不能得到，因此，通過理論分析和實驗規(guī)律建立的物性估算方法是非常重要的。人們希望能從理論上來計算標準生成焓，這就需要建立必要的模型。20世紀70年代中期，米德瑪(Miedema)和其合作者通過聯(lián)合熱力學概念和電子結構理論建立了化合物標準生成焓估算模型[2]，取得了較好的成果，但對于很多金屬間化合物的估算仍存在較大誤差[3]，如部分鑭系物相誤差達到90 kJ/mol，且對于多原子化合物誤差更大；Ray等[4]建立了擴展的米德瑪模型，考慮了彈性能和結構能，應用于估算Zr-Cu-Ni等三元合金體系；拉蒂默的離子束縛模型也是一種很有代表性的估算模型[5]，但是它對于很多金屬間化合物的計算精度也不高，且無法計算Al系金屬間化合物的有關參數(shù)，屬于半理論半經(jīng)驗模型；國內(nèi)一些學者對一些簡單的金屬間化合物進行了實驗研究和計算[6]，但仍然局限于前人的模型，計算的化合物僅針對某一體系，種類有限。金屬間化合物均是由不同原子以不同的個數(shù)組合出來的，每種原子均有本征熱力學屬性(如標準熵等)，同時金屬之間還存在相互作用的熱力學屬性；而本征熱力學屬性、相互作用熱力學屬性、原子個數(shù)和化合物熱力學相應屬性之間是存在一定關系的，通過對大量的化合物的統(tǒng)計計算，郭培民等[7-9]提出了雙參數(shù)估算模型，已用于估算復合氧化物的標準熵、標準生成焓與比熱容，雙參數(shù)模型還能很好地估算金屬間化合物的標準熵和比熱容[10-11]。為了解決目前金屬間化合物的標準生成焓的估算方法不準確問題，本文作者將利用雙參數(shù)模型對金屬間化合物的標準生成焓進行研究。

1 模型建立

對于金屬間化合物MxNy，M和N為組成該化合物的金屬的符號，x和y分別為它們的系數(shù)。假定MxNy的標準生成焓由2個部分組成：一部分為單質(zhì)金屬的貢獻焓；另一部分為單質(zhì)金屬間的相互作用焓。其表達式為：

由于每種單質(zhì)金屬均含有2個參數(shù)，因此，將此模型稱為“標準生成焓的雙參數(shù)模型”。模型中的參數(shù)需要利用已知二元金屬間化合物的標準生成焓通過回歸方法求得。

雙參數(shù)模型中參數(shù)的求解需要利用已知二元金屬間化合物的標準生成焓，因此，二元金屬間化合物標準生成焓的數(shù)據(jù)是否全面、準確和充分，直接關系到模型的可靠程度。在多種無機物熱力學數(shù)據(jù)手冊中查得109個二元金屬間化合物的標準生成焓[12-15]，涉及Fe，Si，Ca，Al和U等34種金屬，相當于每個金屬有 6.4個數(shù)據(jù)，因此，從統(tǒng)計學角度來看，本文的回歸統(tǒng)計也更具有可靠性。在這些二元金屬間化合物中，鈣化物占10個，鐵化物占8個，硅化物占51個，鋁化物占10個，鈾化物占25個，鎳化物占10個，鈦化物占9個等。

在回歸統(tǒng)計中，需要計算平均標準差S和相關系數(shù) R，用以考察模型的可靠性和估算誤差。S越小，表明估算誤差越?。籖越接近 1，表明模型的可信度越高。S和R的計算式分別為：

式中：yi和 y?i分別為某個二元金屬間化合物的實測標準生成焓和估算標準生成焓；y為所有二元金屬間化合物的標準生成焓平均值。

2 計算結果

通過數(shù)據(jù)分析和模型求解，能夠得到34種金屬的A，A′和常數(shù) D(-1.119 kJ/mol，298 K)，結果見表 1；估算了金屬間化合物的標準生成焓，部分數(shù)據(jù)見表2。利用雙參數(shù)模型，估算的平均誤差為12.47 kJ/mol，相關系數(shù)為0.992。從統(tǒng)計學角度來看，雙參數(shù)模型是可靠的。109個數(shù)據(jù)的標準差S為17.40 kJ/mol，其中誤差大于3S(52.20 kJ/mol)的數(shù)據(jù)只有1個，占總數(shù)據(jù)的0.92%；大于2S(34.80 kJ/mol)的數(shù)據(jù)有6個，占總數(shù)據(jù)的5.51%；大于1S(17.40 kJ/mol)的有31個，占總數(shù)據(jù)的28.44%。圖1所示為標準生成焓誤差的分布。從圖1可見：雙參數(shù)模型109組數(shù)據(jù)的標準生成焓誤差符合正態(tài)分布。

圖1 標準生成焓誤差的分布Fig.1 Error distribution of standard enthalpy

表1 金屬的A和A′Table 1 A and A′ of metal kJ/mol

表2 部分金屬間化合物的標準生成焓Table 2 Part of standard enthalpies of intermetallics kJ/mol

3 各種模型的比較

建立金屬間化合物的估算模型，是用來預測未知物相的標準生成焓，如果模型本身估算的數(shù)據(jù)與已測實驗數(shù)據(jù)相差很大，那么用模型估算未知物相的可靠性就很低。金屬間化合物數(shù)據(jù)眾多，其中硅系化合物個數(shù)最多，其熱力學數(shù)據(jù)經(jīng)過前人的實驗測得，以它作為比較對象具有代表性。

米德瑪及其合作者建立了米德瑪模型用以估算金屬間化合物的標準生成焓：

Miedema等[16-17]考慮的合金化效應及他們用來描述這種效應所使用的2個半經(jīng)驗參數(shù)nws和φ*均屬于電子因素[16-17]，式(1)描述了形成固態(tài)或液態(tài)合金的2種原子參數(shù)的能量影響，第1項是由2種原子的晶?；瘜W勢之差引起的；第2項為在2種不同原子晶粒邊界和威格納-賽茨(Wigner-Seitz)純晶粒邊界形成連續(xù)電子密度所需能量；最后1項是修正項。P和R均是用經(jīng)驗的方法調(diào)整得到的常數(shù)。米德瑪模型能估算至少某一組分為過渡金屬元素形成二元合金的反應熱。

對于金屬間化合物應用離子束縛模型，可以推導出二元金屬間化合物MxNy的標準生成焓為：

式中：xi為化合物中組元i的濃度；Ei和Ci為根據(jù)實驗數(shù)據(jù)和估算模型用最小二乘法確定的參數(shù)。

這里必須指出，離子束縛模型不適用于含鋁的金屬間化合物，否則將產(chǎn)生很大的偏差。

將米德瑪模型和離子束縛模型與雙參數(shù)模型估算的硅系化合物標準生成焓誤差進行了比較，數(shù)據(jù)見表3。令金屬間化合物的實測值與估算值之差的絕對值為絕對誤差(正數(shù))，用 dj表示；絕對誤差與實測絕對值之比的百分數(shù)為相對誤差，用dx表示。即

圖2所示為3種模型的絕對誤差比較。從圖2可看出：米德瑪模型估算的標準生成焓的絕對誤差較大，雙參數(shù)模型標準生成焓的絕對誤差dj大多比離子束縛模型的要小，51種化合物中，有33種金屬間化合物的誤差比離子束縛模型的小。并且離子束縛模型有 6組數(shù)據(jù)絕對誤差大于50 kJ/mol，最大達109.54 kJ/mol，而雙參數(shù)模型最大誤差也只有46.40 kJ/mol，且波動較小。圖3所示為3種模型的相對誤差比較。從圖3可看出：雙參數(shù)模型的相對誤差dx比離子束縛模型的也要小，誤差較大的也很少。

雙參數(shù)模型估算金屬間化合物標準生成焓的標準差為17.17 kJ/mol，而離子束縛模型為31.96 kJ/mol，米德瑪模型為67.09 kJ/mol，雙參數(shù)模型要比離子束縛模型和米德瑪模型的標準差要小很多，指標比離子束縛模型和米德瑪模型的優(yōu)。圖4所示是2種模型偏差的正態(tài)分布。由圖4可見：雙參數(shù)模型偏差的正態(tài)分布呈尖峰形，說明偏差比較集中，都在0附近；而米德瑪模型和離子束縛模型的偏差正態(tài)分布呈禿峰形，說明偏差較大且分散，其中米德瑪模型的分散程度最大。因此，金屬間化合物標準生成焓的雙參數(shù)模型對數(shù)據(jù)的估算比離子束縛模型和米德瑪模型有更高的準確度。

表3 硅系金屬間化合物估算誤差比較Table 3 Estimating error comparison of intermetallics containing silicon kJ/mol

圖2 3種模型的絕對誤差比較Fig.2 Comparisons of absolute error of three models

圖3 3種模型的相對誤差比較Fig.3 Comparisons of relative error of three models

圖4 硅系金屬間化合物標準生成焓誤差分布Fig.4 Error distributions of standard enthalpies of intermetallics containing silicon

4 結論

(1) 建立了金屬間化合物的雙參數(shù)模型，用來估算金屬間化合物的標準生成焓。

(2) 利用雙參數(shù)模型，估算的金屬間化合物的平均誤差為12.47 kJ/mol，相關系數(shù)為0.992，109個數(shù)據(jù)的標準差為17.40 kJ/mol。從統(tǒng)計學角度來看，雙參數(shù)模型是可靠的。

(3) 雙參數(shù)模型估算的標準生成焓要比離子束縛模型的標準差要小，其絕對誤差和相對誤差均比米德瑪模型和離子束縛模型的小，指標比米德瑪模型和離子束縛模型的優(yōu)。

[1] 陳星秋, 嚴新林, 丁學勇, 等. 化合物生成焓：一百年和密度泛函基量子機制的原子模型新時代[J]. 中國稀土學報, 2004,22(增刊): 2-3.CHEN Xing-qiu, YAN Xin-lin, DING Xue-yong, et al.Enthalpies of formation for compounds: a century and comes of age for density functional based computations[J]. Journal of the Chinese Society of Rare Earth, 2004, 22(z1): 2-3.

[2] Latimer W M. Methods of estimating the entropies of solid compounds[J]. J Am Chem Soc, 1951, 73(4): 1480-1482.

[3] 湯振雷, 王為. 計算合金系統(tǒng)熱力學性質(zhì)的 Miedema模型的發(fā)展[J]. 材料導報, 2008, 22(3): 117-118.TANG Zhen-lei, WANG Wei. Review of Miedema model for calculating thermodynamic proper ties of alloys system[J].Materials Review, 2008, 22(3): 117-118.

[4] Ray P K, Akinc M, Kramer M J. Applications of an extended Miedema model for ternary alloys[J]. Journal of Alloys and Compounds, 2010, 489(2): 357-361.

[5] 喬芝郁, 許志宏, 劉洪霖. 冶金和材料計算物理化學[M]. 北京: 冶金工業(yè)出版社, 1999: 65-69.QIAO Zhi-yu, XU Zhi-hong, LIU Hong-lin. Computerized physical chemistry of metallurgy and materials[M]. Beijing:Metallurgical Industry Press, 1999: 65-69.

[6] 周鴻翼, 劉天模, 王金星. Mg-Zn合金的熱力學性質(zhì)計算[J].重慶大學學報, 2006, 29(12): 68-70.ZHOU Hong-yi, LIU Tian-mo, WANG Jin-xing. Calculation of the thermodynamic properties of Mg-Zn alloy[J]. Journal of Chongqing University, 2006, 29(12): 68-70.

[7] 郭培民, 趙沛. 用雙參數(shù)模型估算復合氧化物的標準熵[J].鋼鐵研究學報, 2006, 18(9): 17-20.GUO Pei-min, ZHAO Pei. Estimation of standard entropy of complex oxide using two-parameter model[J]. Journal of Iron and Steel Research, 2006, 18(9): 17-20.

[8] 郭培民, 趙沛. 雙參數(shù)模型預報復合氧化物的生成焓[J]. 鋼鐵研究學報, 2007, 19(5): 25-28.GUO Pei-min, ZHAO Pei. Estimation of standard enthalpies of complex oxide using two-parameter model[J]. Journal of Iron and Steel Research, 2007, 19(5): 25-28.

[9] 郭培民, 趙沛. 雙參數(shù)模型估算復合氧化物比熱的研究[J].鋼鐵研究學報, 2008, 20(8): 5-7.GUO Pei-min, ZHAO Pei. Estimation of specific heat of complex oxide using two-parameter model[J]. Journal of Iron and Steel Research, 2008, 20(8): 5-7.

[10] 趙定國, 郭培民, 趙沛. 金屬間化合物的標準熵估算模型[J].北京科技大學學報, 2010, 32(1): 27-31.ZHAO Ding-guo, GUO Pei-min, Zhao Pei. Estimating model for standard entropy of intermetallics[J]. Journal of University of Science and Technology Beijing, 2010, 32(1): 27-31.

[11] 趙定國, 郭培民, 趙沛. 金屬間化合物的比熱容估算模型[J].北京科技大學學報, 2010, 32(5): 565-569.ZHAO Ding-guo, GUO Pei-min, ZHAO Pei. Estimating model for specific heat of intermetallics[J]. Journal of University of Science and Technology Beijing, 2010, 32(5): 565-569.

[12] 葉大倫, 胡建華. 實用無機物熱力學數(shù)據(jù)手冊[M]. 北京: 冶金工業(yè)出版社, 2002: 10-102.YE Da-lun, HU Jian-hua. Applied thermodynamics data of mineral inorganic substance[M]. Beijing: Metallurgical Industry Press, 2002: 10-102.

[13] David R L. CRC handbook of chemistry and physics 88th edition[M]. Boca Raton: Chemical Rubber Company Press, 2008:10-213.

[14] Dean J A. 蘭氏化學手冊[M]. 15版. 北京: 科學出版社, 2007:125-512.Dean J A. Lange’s chemistry handbook[M]. 15th ed. Beijing:Science Press, 2007: 125-512.

[15] Brain I. 純物質(zhì)熱化學數(shù)據(jù)手冊[M]. 北京: 科學出版社, 2003:1-188.Brain I. Thermochemical data of pure substances[M]. Beijing:Science Press, 2003: 1-188.

[16] 鄭志剛. 合金熱力學性質(zhì)的Miedema理論計算[D]. 南寧: 廣西大學資源與冶金學院, 2005: 6-9.ZHENG Zhi-gang. Calculation of thermodynamic properties of the multi-component alloys with Miedema theory[D]. Nanning:Guangxi University. College of Resources and Environment,2005: 6-9.

[17] WANG Wen-chao, LI Jia-hao, YAN Hua-feng, et al. A thermodynamic model proposed for calculating the standard formation enthalpies of ternary alloy systems[J]. Scripta Materialia, 2007, 56(11): 975-978.