王猛杰, 常春然, 徐抗震, 馬海霞, 趙鳳起

(1. 西北大學(xué)化工學(xué)院，陜西 西安 710069；2. 西安交通大學(xué)化學(xué)工程與技術(shù)學(xué)院，陜西 西安 710049；3. 西安近代化學(xué)研究所燃燒與爆炸技術(shù)重點實驗室，陜西 西安 710065)

?

鐵酸系列復(fù)合金屬氧化物鋁熱反應(yīng)焓的理論研究

王猛杰1, 常春然2, 徐抗震1, 馬海霞1, 趙鳳起3

(1. 西北大學(xué)化工學(xué)院，陜西 西安 710069；2. 西安交通大學(xué)化學(xué)工程與技術(shù)學(xué)院，陜西 西安 710049；3. 西安近代化學(xué)研究所燃燒與爆炸技術(shù)重點實驗室，陜西 西安 710065)

運用密度泛函理論(DFT)，在Material Studio程序包的GGA-PBE、 GGA-BLYP、 GGA-PW91泛函結(jié)合DNP基組水平上，計算了7種常見金屬氧化物的鋁熱反應(yīng)焓。將計算值與實驗值進行對比，確定出最佳計算基組水平。在此基礎(chǔ)上，對6種鐵酸系列復(fù)合金屬氧化物(CuFe2O4、NiFe2O4、CoFe2O4、MgFe2O4、ZnFe2O4、MnFe2O4)的鋁熱反應(yīng)焓進行了理論計算，并運用蓋斯定律導(dǎo)出其標準摩爾生成焓。結(jié)果表明，GGA-PBE/DNP方法計算的標準摩爾反應(yīng)焓精度高、誤差小，平均誤差7.072kJ/mol；在GGA-PBE/DNP水平下，6種復(fù)合金屬氧化物的鋁熱反應(yīng)焓分別為-3695.02、-3388.53、-3380.13、-841.06、-3142.57和-2738.40kJ/mol，與等量物理混合金屬氧化物的鋁熱反應(yīng)焓相差不大。6種復(fù)合金屬氧化物標準摩爾生成焓為-992.96、-1092.12、-1090.13、-1431.13、-1185.15和-1311.78kJ/mol。

量子化學(xué)；密度泛函理論;DFT；鋁熱反應(yīng)；復(fù)合金屬氧化物

引 言

復(fù)合金屬氧化物具有耐高溫、化學(xué)穩(wěn)定性好、硬度高等特點，在功能材料、電極材料、氣敏材料、磁性材料、催化材料等[1-6]方面具有巨大的應(yīng)用價值。近年來，一些學(xué)者將復(fù)合金屬氧化物代替金屬氧化物作為燃燒催化劑應(yīng)用到推進劑領(lǐng)域，發(fā)現(xiàn)其優(yōu)越的性能不僅能極大地提高推進劑燃速，減小壓力指數(shù)，而且能使燃燒更加穩(wěn)定[7-8]。

超級鋁熱劑在燃燒彈、微型推進器、含能表面涂層、納米焊接和推進劑等方面應(yīng)用廣泛[9-11]。復(fù)合金屬氧化物與鋁的鋁熱反應(yīng)是一類重要反應(yīng)，然而由于復(fù)合氧化物優(yōu)異的穩(wěn)定性，很難用實驗方法測得其鋁熱反應(yīng)焓，因此，其基礎(chǔ)熱力學(xué)數(shù)據(jù)比較匱乏。目前，對復(fù)合金屬氧化物與鋁的鋁熱反應(yīng)報道較少。

本研究設(shè)計7種常見氧化物的鋁熱反應(yīng)，在3種泛函水平下計算了其鋁熱反應(yīng)焓，篩選出GGA-PBE/DNP作為研究的適宜方法和基組，在該水平下，進一步對復(fù)合金屬氧化物進行理論計算研究。此外，依照得到的復(fù)合金屬氧化物與鋁的反應(yīng)熱，通過蓋斯定律，得到其在0K下的生成焓，并依照基?；舴蚬降玫狡湓?98.15K下的標準摩爾生成焓，以期為進一步研究復(fù)合金屬氧化物的熱性能提供有效的熱力學(xué)參數(shù)。

1 計算方法

1.1 晶胞的構(gòu)建

根據(jù)ICSD數(shù)據(jù)庫，構(gòu)建了氧化鐵(Fe2O3)、氧化銅(CuO)、氧化鎳(NiO)、二氧化鉛(PbO2)、三氧化鎢(WO3)、三氧化鉬(MoO3)、三氧化二鉻(Cr2O3)和鋁等物質(zhì)的晶胞，其中Al2O3和Fe2O3兩種金屬氧化物的晶胞結(jié)構(gòu)如圖1所示。由圖1可見,Al2O3和Fe2O3的空間結(jié)構(gòu)非常相近，空間構(gòu)型均屬于R-3CH,空間群編號為167(晶格常數(shù)：a=b=0.475nm;c=1.297nm)，每個Al2O3晶胞中帶有12個Al3+和18個O2-。

圖1 Al2O3和 Fe2O3的晶胞構(gòu)型Fig.1 Unit cell configurations of Al2O3 and Fe2O3

鐵酸銅(CuFe2O4)、鐵酸鎳(NiFe2O4)等復(fù)合金屬氧化物屬于類鈣鈦礦型復(fù)合氧化物，其典型的晶胞結(jié)構(gòu)見圖2。從圖2可以看出，兩者都屬于AB2O4型類鈣鈦礦型復(fù)合氧化物，空間構(gòu)型為FD-3MS，屬于立方晶系11，空間群編號227。尖晶石[12]晶胞可以劃分成8個小的立方單位，分別由4個A型和4個B型小單位拼在一起。每個A型、B型小單位都有4個O2-離子，晶胞中O2-的個數(shù)是32個。Cu2+(Ni2+)處于A型小單位的中心及一半的頂點及B型小單位一半的頂點上，晶胞中Cu2+(Ni2+)的數(shù)目是8個。Cu2+(Ni2+)呈四配位，即占據(jù)O2-密堆積中的四面體空隙。每個B型小單位中有4個Fe3+，晶胞中Fe3+的個數(shù)是16個。Fe3+呈八配位，即占據(jù)O2-密堆積中的八面體空隙。每個鐵原子分別與4個氧原子相連，每個銅(鎳)原子分別與6個氧原子相連，構(gòu)成尖晶石型化合物。

圖2 鐵酸銅的晶胞結(jié)構(gòu)Fig.2 The unit cell structure of CuFe2O4

1.2 計算方法

運用密度泛函理論(DFT)方法，在GGA-PBE/DNP、GGA-BLYP/DNP、GGA-PW91/DNP水平下計算了相關(guān)氧化物和金屬在0K下完美晶型的總能量(E)，計算結(jié)果列于表1。通過與實驗值進行對比擬合，確定出計算值與實驗值相關(guān)性最好的基組水平。晶體模擬均在Material Studio軟件包的DMol3模塊下進行，晶胞中k-point 網(wǎng)格參數(shù)為9×9×9，軌道截止能量為0.50nm，SCF收斂值為10×e-6。

表1 不同泛函水平下各個分子的總能量

注: 1hatree=2625.5kJ/mol

1.3 鋁熱反應(yīng)焓計算

以Fe2O3的鋁熱反應(yīng)為例，計算其在0K下的摩爾反應(yīng)焓(ΔrHm)。據(jù)能量守恒定律: ΔH=Q

Fe2O3+2Al=2Fe+Al2O3

(1)

EFe2O3+EAl=EFe+EAl2O3+Q

(2)

ΔH=(Ef,p-Ef,R) × 2625.5

(3)

ΔH298K=ΔH0K+ΔHT

(4)

式中:Ef,R和Ef,p分別為0K完美晶型下反應(yīng)物和生成物的總能量。ΔHT為溫度校正項，Ef,R=EFe2O3+2EAl,Ef,p=EAl2O3+2EFe,其值可以通過基希霍夫公式計算得到。根據(jù)式(2)、(3)、(4)，可得到Fe2O3的鋁熱反應(yīng)焓。常見金屬氧化物的鋁熱反應(yīng)焓可按上述計算方法得到。

對于復(fù)合金屬氧化物，以計算CuFe2O4的鋁熱反應(yīng)焓為例：

Fe2O3+2Al=2Fe+Al2O3

(5)

3CuO+2Al=Al2O3+3Cu

(6)

3Fe2O3+8Al+3CuO=6Fe+4Al2O3+3Cu

(7)

3CuFe2O4+8Al=6Fe+4Al2O3+3Cu

(8)

反應(yīng)(7)為反應(yīng)(5)與反應(yīng)(6)的物理混合體系，且反應(yīng)(7)的反應(yīng)物所含離子數(shù)與反應(yīng)(8)相等；根據(jù)能量守恒定律：反應(yīng)(7)和反應(yīng)(8)的鋁熱反應(yīng)焓應(yīng)相等，即ΔH(7)=ΔH(8)。又，ΔH(7)=3ΔH(5)+ΔH(6)。根據(jù)反應(yīng)(5)與(6)的鋁熱反應(yīng)焓，可求得ΔH(7)。而Al、Fe、Al2O3和Cu的鋁熱反應(yīng)焓(E)見表1，同時存在下列關(guān)系式：

ΔE=Ep-ER

(9)

EP=6EFe+4EAl2O3+3ECu

(10)

ER=3ECuFe2O4+8EAl

(11)

ΔH(8)=ΔE×2625.5kJ/mol

(12)

式中:EP和ER分別為反應(yīng)物和生成物的總能量;ΔH(8)為0K下復(fù)合金屬氧化物與鋁的摩爾反應(yīng)焓(ΔrHm)。

1.4 標準摩爾生成焓的計算

根據(jù)蓋斯定律：

(13)

2 結(jié)果與討論

2.1 鋁熱反應(yīng)焓誤差分析

在PBE、PW91和BLYP泛函水平下，計算得到7種常見氧化物的鋁熱反應(yīng)焓，結(jié)果見表2。由表2可知，GGA-PBE/DNP與GGA-BLYP/DNP條件下的計算值與實驗值的線性相關(guān)系數(shù)均大于0.9990，兩種泛函沒有較大的區(qū)別，都可以用于計算總能量。GGA-BLYP/DNP水平下，計算所得的最大偏差和最小偏差分別是-27.020kJ/mol和-11.485kJ/mol。GGA-PBE/DNP水平下，7種氧化物的標準摩爾生成焓與實驗值的均方根誤差最小。而GGA-PW91/DNP泛函下的計算值與實驗值偏差較大，原因主要有兩點：首先不同的泛函對應(yīng)的研究對象不同，GGA-PW91是一類沒有修正過的泛函，它不適用于金屬的研究計算；其次，對于研究對象，一部分晶體是軟件中自帶的，一部分是從ICSD數(shù)據(jù)庫導(dǎo)入的cif文件計算所得，還有一部分是查詢所得的晶胞參數(shù)自建晶體所得，導(dǎo)致得到的晶體存在一定差距。因此，GGA-PBE計算所得數(shù)據(jù)誤差最小?？梢姡瑢τ谶M一步計算二元復(fù)合金屬氧化物的鋁熱反應(yīng)焓，GGA-PBE/DNP水平的精確度較高。

表2 7種金屬氧化物鋁熱反應(yīng)焓計算值與實驗值的比較

注：ΔrHm為鋁熱反應(yīng)焓；δ為絕對偏差，為計算值與實驗值的差值;R為計算值與實驗值的線性相關(guān)系數(shù)；RMSE為 均方根誤差；ME為平均偏差。

從表2還可以看出，GGA-PBE/DNP水平下CuO、PbO2、MoO3、NiO等7種氧化物鋁熱反應(yīng)焓的計算值與實驗值平均誤差相差不大。這表明系統(tǒng)誤差與晶胞的大小無關(guān)，因此GGA-PBE/DNP可用于計算復(fù)合金屬氧化物這類晶胞體積較大的物質(zhì)。

2.2 復(fù)合金屬氧化物的鋁熱反應(yīng)焓

表3為GGA-PBE/DNP水平下計算得到的鐵酸系列復(fù)合金屬氧化物的鋁熱反應(yīng)焓。所有反應(yīng)值均是在0K下計算得到。由于MgO與Al不發(fā)生反應(yīng)，所以MgFe2O4/Al的物理復(fù)合物中只有Fe2O3與Al反應(yīng)；其余5種物理混合物的鋁熱反應(yīng)焓與復(fù)合氧化物的鋁熱反應(yīng)焓相比，相對偏差分別是0.59%、0.91%、2.93%，0.22%和0.328%。由此更近一步證實了鐵酸系列復(fù)合金屬氧化物與其相應(yīng)的物理混合物能量相差不大，因此可以運用其數(shù)據(jù)計算復(fù)合金屬氧化物的標準摩爾生成焓。

表3 GGA-PBE/DNP水平下復(fù)合金屬氧化物的鋁熱反應(yīng)焓

注：1～6號物理混合物分別為CuO與Fe2O3，NiO與Fe2O3，CoO與Fe2O3，MgO與Fe2O3，ZnO與Fe2O3，MnO與Fe2O3。

2.3 復(fù)合金屬氧化物的熱性能

在0K下，通過表3中CuFe2O4鋁熱反應(yīng)焓，由公式(13)計算得到CuFe2O4的標準摩爾生成焓。其他鐵酸系列復(fù)合金屬氧化物的標準摩爾生成焓可根據(jù)計算CuFe2O4方法得到，結(jié)果列于表4。

表4 GGA-PBE/DNP水平下復(fù)合金屬氧化物的標準摩爾生成焓

圖3 GGA-PBE/DNP水平下計算得到的6種復(fù)合金屬氧化物的標準摩爾生成焓計算值與文獻值的擬合曲線Fig.3 Fitted curves of the calculated and reference values of the standard molar enthalpies of formation of six kinds of metal composite oxides at GGA-PBE/DNP level

3 結(jié) 論

(1)通過比較GGA-PBE/DNP、GGA-BLYP/DNP、GGA-PW91/DNP三種方法計算所得常見金屬氧化物和Al的標準摩爾反應(yīng)焓與實驗值，得出GGA-PBE/DNP方法具有精度高、誤差小的特點,可估算復(fù)合金屬氧化物與Al的標準摩爾反應(yīng)焓。

(2)6種復(fù)合金屬氧化物的鋁熱反應(yīng)焓與其等量物理混合物的鋁熱反應(yīng)焓的差值分別為21.89、-30.6、-96.39、5.21、-7.05和9.02kJ/mol，表明等量物理混合物與其相應(yīng)的復(fù)合金屬氧化物所帶的能量相差不大。

(3)6種復(fù)合金屬氧化物中，最穩(wěn)定的是MgFe2O4，最不穩(wěn)定的是CuFe2O4。GGA-PBE/DNP方法計算所得CuFe2O4、NiFe2O4、CoFe2O4、MgFe2O4、ZnFe2O4和MnFe2O4的標準摩爾生成焓分別為-992.62、-1092.12、-1090.13、-1431.13、-1185.15和-1311.78kJ/mol。

[1] Aakash，Choubey R，Das D，et al. Effect of doping of manganese ions on the structural and magnetic properties of nickel ferrite[J]. Journal of Alloys and Compounds，2016，668(6)：33-39.

[2] 石曉波，傅海萍，李春根. 鐵酸鎳納米微粒的固相合成及其催化性能[J]. 應(yīng)用化學(xué)，2003，20(1)：55-58. SHI Xiao-bo，F(xiàn)U Hai-ping，LI Chun-gen. Solid state synthesis and catalytic activity of nano-NiFe2O4[J]. Chinese Journal of Applied Chemistry，2003，20(1)：55-58.

[3] WEI T T，ZHANG Y，XU K Z，et al. Catalytic action of nano Bi2WO6on thermal decompositions of AP，RDX，HMX and combustion of NG/NC propellant[J]. RCS Advance，2015，5(86)：70323-70328.

[4] ZHANG Y，WEI T T，XU K Z，et al. Catalytic decomposition action of hollow CuFe2O4nanospheres on RDX and FOX-7[J]. RCS Advance，2015，5(92)：75630-75635.

[5] Maqusood A，Mohd J A，Hisham A A，et al. Copper ferrite nanoparticle-induced cytotoxicity and oxidative stress in human breast cancer MCF-7 cells[J]. Colloids and Surfaces B，2016，142:46-54.

[6] 王虹，王軍利，李翠清，等. ACo2O4/HZSM-5催化劑上N2O的直接分解[J]. 物理化學(xué)學(xué)報，2010，26(10)：2739-2744.WANG Hong，WANG Jun-li，LI CUI-qing，et al. Decomposition of N2O on ACo2O4/HZSM-5 catalysts[J]. Acta Physico-Chimica Sinica,2010，26(10)：2739-2744.[7] 安亭，趙鳳起，裴慶，等. 超級鋁熱劑的制備、表征及其燃燒催化作用[J]. 無機化學(xué)學(xué)報，2011，27(2)：231-238.

AN Ting，ZHAO Feng-qi，PEI Qing，et al. Preparation，characterization and combustion catalytic activity of nanopartical super thermites[J]. Chinese Journal of Inorganic Chemistry，2011，27(2)：231-238.

[8] 安亭，趙鳳起，裴慶，等. 超級鋁熱劑對雙基推進劑燃燒性能的影響[J]. 固體火箭技術(shù)，2012，35(6)：773-778.

AN Ting, ZHAO Feng-qi, PEI Qing, et al. Effects of super thermites on combustion properties of double-base propellant[J]. Journal of Solid Rocket Technology，2012，35(6)：773-778.

[9] 馮曉軍，王曉峰，徐洪濤，等. Al粉對炸藥爆炸加速的影響[J]. 火炸藥學(xué)報，2014，37(5)：25-32.

FENG Xiao-jun，WANG Xiao-feng,XU Hong-tao，et al. Influence of Al powder on the explosion acceleration ability for explosive[J]. Chinese Journal of Explosives and Propellants(Huozhayao Xuebao)，2014，37(5)：25-32.

[10] 王軍，張文超，沈瑞琪，等. 納米鋁熱劑的研究進展[J]. 火炸藥學(xué)報，2014，37(4)：1-8.

WANG Jun，ZHANG Wen-chao，SHEN Rui-qi，et al. Research of progress on nano thermite[J]. Chinese Journal of Explosives & Propellants(Huozhayao Xuebao)，2014，37(4)：1-8.

[11] Zhao N N，He C C，Liu J B，et al. Dependence of catalytic properties of Al/Fe2O3thermites on morphology of Fe2O3particles in combustion reactions[J]. Journal of Solid State Chemistry，2014，219(11)，67-73.

[12] 熊熙然， 張玲玲，蒼大強. 不同晶體結(jié)構(gòu)對鐵酸銅發(fā)射率影響的研究[J]. 冶金能源，2013，32(1)：15-19.

XIONG Xi-ran，ZHANG Ling-ling，CANG Da-qiang. Study on different crystal structure impact on infrared emissivity of CuFe2O4[J]. Energy for Metallurgical Industry，2013，32(1)：15-19.

[13] Ye Da-lun，Hu Jian-hua. Handbook of Thermodynamic Data for Practical Inorganic Materials [M]. Beijing：Metallurgical Industry Press，2002.

[14] Lide D R. CRC Handbook of Chemistry and Physics [M]. Boca Raton:CRC Press，2004.

Theoretical Study of Thermite Reaction Enthalpy for Metal Composite Oxides of Ferrite Series

WANG Meng-jie1,CHANG Chun-ran2,XU Kang-zhen1,MA Hai-xia1,ZHAO Feng-qi3

(1.School of Chemical Engineering, Northwest University, Xi′an 710069, China; 2. School of Chemical Engineering and Technology, Xi′an Jiaotong University, Xi′an 710049, China; 3. Science and Technology on Combustion Laboratory,Xi′an Modern Chemistry Research Institute, Xi′an 710065, China)

By using the density functional theory (DFT), the thermite reaction enthalpies of 7 kinds of common metal oxides were calculated with GGA-PBE, GGA-BLYP, GGA-PW91 functional in conjuction with DNP basis set in the Material Studio package. The calculated values were compared with the experimental ones to determine the optimum calculation basis set level. On this basis, the thermite reaction enthalpies of six kinds of metal composite oxides of ferrite series (CuFe2O4, NiFe2O4, CoFe2O4, MgFe2O4, ZnFe2O4and MnFe2O4) were theoretically calculated, and Hess′s law was employed to calculate their standard molar enthalpies of formation. The results show that the claclulation method based on the GGA-PBE functional in conjuction with DNP basis set is most precise, and the mean error (ME) of common thermite reactions between calculated and experimental enthalpies is 7.072kJ/mol. At GGA-PBE/DNP level, thermite reaction enthalpies of the six compounds are -3695.02,-3388.53, -3380.13,-841.06, -3142.57 and -2738.40 kJ/mol, respectively,which is little different with that of the equivalent physical mixed metal oxides. The standard molar enthalpies of formation of six substances are -992.62, -1092.12, -1090.13, -1431.13, -1185.15 and -1311.78 kJ/mol, respectively.

quantum chemistry;density functional theory;DFT; thermite reaction; metal composite oxides

10.14077/j.issn.1007-7812.2016.05.007

2016-07-15；

2016-09-06

國家自然科學(xué)基金(No.21673178;No.21241003)；中國博士后科學(xué)基金資助 (No.2014M552480)

王猛杰(1992-)，男，碩士，從事含能材料的理論計算研究。E-mail: 492048485@qq.com

徐抗震(1976-)，男，博士，教授，從事含能材料研究。 E-mail: xukz@nwu.edu.cn

TJ55;TQ203

A

1007-7812(2016)05-0046-05