1H-1，2，4-三唑衍生物爆轟性能的理論計算

馬忠亮，周 華，王建龍，李景滿

（1．中北大學化工與環境學院，山西太原030051；2．四川通達化工有限責任公司，四川達州635000）

1H-1，2，4-三唑衍生物爆轟性能的理論計算

馬忠亮1，周 華1，王建龍1，李景滿2

（1．中北大學化工與環境學院，山西太原030051；2．四川通達化工有限責任公司，四川達州635000）

采用密度泛函理論（DFT）B3LYP方法，在6-31G（d，p）基組水平下對1H-1，2，4-三唑衍生物進行幾何構型全優化，設計了合適的等鍵反應計算1 H-1，2，4-三唑衍生物的生成焓。基于0.001e·bohr－3等電子密度面所包圍的體積空間求得分子平均摩爾體積（V）和理論密度（ρ），運用Kamlet-Jacobs方程估算1H-1，2，4-三唑衍生物的爆速（D）和爆壓（p），結果表明，1H-1，2，4-三唑發生N-取代或引入－N3基團都能顯著提高衍生物的生成焓，－NF2或－ONO2基團能夠增強衍生物的爆轟性能。通過計算弱鍵的化學解離能（BDE），得出部分1H-1，2，4三唑衍生物具有較好的熱穩定性。

有機化學；量子化學；密度泛函理論；1H-1，2，4三唑；生成焓；爆轟性能；熱穩定性

引 言

低感度、高能量密度化合物是含能材料研究領域的重要方向。三唑是含有3個氮原子的五元環化合物，分子結構中的高氮（50.9%）低碳氫含量不僅使其具有較高的密度，而且更容易達到氧平衡，其衍生物有可能作為高能量密度材料［1］。鈍感炸藥5-氨基-3-硝基-1，2，4-三唑（ANTA）和3-硝基-1，2，4-三唑-5-酮（NTO）是1，2，4-三唑重要的衍生物，ANTA可用于火箭推進劑，NTO及其鹽類是具有廣泛應用前景的鈍感含能材料。

近年來，三唑衍生物作為含能材料或含能材料中間體備受關注。Krishna kumar等［2］研究了3，5-二硝基-1，2，4-三唑分子間的相互作用；XU L N等［3］研究了NTO分子間的相互作用；Siker A K等［4］研究了ANTA及其衍生物的熱穩定性；朱佳平等［5］計算了DNMT的爆轟性能；魯亞琳［6］研究了1，2，3-三唑及1，2，4-三唑的前線軌道能及熱力學性質。為了深入研究取代基團對1 H-1，2，4-三唑性能的影響，本研究對1H-1，2，4-三唑衍生物的爆轟性能進行了系統的計算，以期為1 H-1，2，4-三唑的實驗研究提供參考。

1 計算方法

根據取代基種類、位置及數量的不同，將1H-1，2，4-三唑衍生物分為A、B、C、D和E五類（如圖1和表1所示）。運用Gaussian03程序［7］，采用密度泛函理論（DFT）B3LYP［8-9］方法結合6-31G（d，p）基組，對1 H-1，2，4-三唑衍生物進行分子幾何構型優化計算。經振動頻率分析，確定優化構型具有勢能面上極小點（無虛頻），得到分子的總能量。

圖1 1H-1，2，4三唑衍生物的分子結構Fig.1 Molecular structure of 1H-1，2，4-triazole derivatives

表1 1H-1，2，4-三唑衍生物的取代基Table 1 Substituent of 1H-1，2，4-triazole derivatives

等鍵反應廣泛用于計算化合物的生成焓［10-11］。本研究設計如下等鍵反應：

根據式（1）所確定的等鍵反應計算化合物的氣相生成焓，由于含能材料大部分以固相形式存在，在估算爆轟性能時應先計算其固相生成焓。根據蓋斯定律：ΔHf，solid＝ΔHf，gas－ΔHsub，因此，只需計算升華焓，即可得到目標化合物的固相生成焓。Politzer［12-14］計算升華焓的經驗公式為：

式中：a＝11.166×10－4kJ·mol－1A－4；b＝6.900 kJ·mol－1；c＝12.404 kJ·mol－1［15］

運用Monte-Carlo方法由0.001e／Bohr3等電子密度面所包圍的體積求得理論密度，采用公式（3）修正求得其晶體密度：

式中：α＝0.9183；β＝0.0028；γ＝0.0443［16］

在求得固相生成焓和晶體密度的基礎上，采用Kamlet-Jacobs公式［17-18］估算了1 H-1，2，4-三唑衍生物的爆速和爆壓，篩選出爆轟性能較好的化合物。

為了考察化合物的熱穩定性，以鍵解離能（BDE）作為衡量化合物熱穩定性的一個指標，采用公式（6）計算1 H-1，2，4-三唑衍生物中較弱化學鍵的鍵解離能：

2 結果與討論

2.1 生成焓

表1列出7種相關化合物的氣相生成焓，其中1 H-1，2，4-三唑（C2H3N3）、CH4、CH3NH2、CH3NO2和CH3ONO2生成焓的實驗值來自文獻［19-21］。由表1可看出，1H-1，2，4-三唑、CH4、CH3NN2、CH3NO2和CH3ONO2的氣相生成焓計算值與實驗值的相對誤差分別為0.39%、3.77%、0.71%、4.32%、8.89%，表明在G2水平下相關化合物的計算值可信。

根據等鍵反應式（1）計算1 H-1，2，4-三唑衍生物的氣相生成焓、升華焓以及固相生成焓，結果見表3。從表3可以看出，1 H-1，2，4-三唑衍生物的生成焓受取代基種類、位置以及數量的影響。觀察A、B、C系列可見，發生C－取代時引入－N3基團能大幅度提高衍生物的生成焓，引入－NH2、－NF2或－ONO2基團則導致衍生物的生成焓降低。發生N－取代時（D系列）引入－N3、－NO2、－NH2、－NF2或－ONO2基團都能提高衍生物生成焓，引入基團能夠提高衍生物生成焓的順序為：－N3＞－NO2＞－NH2＞－NF2＞－ONO2。

表2 相關化合物的氣相生成焓及相對誤差Table 2 Heats of formation and relative error for the related compounds

表3 1H-1，2，4-三唑衍生物的總能量、零點能、熱校正和生成焓的計算值Table 2 Calculated values of total energies，zero-point energies，thermal correction，heats of formation for the 1H-1，2，4-triazole derivatives

從表3看出，發生C－取代的衍生物，取代基位置對衍生物生成焓的影響并不明顯，而發生N－取代的衍生物生成焓則大幅提高。

1 H-1，2，4-三唑衍生物引入基團數與生成焓的關系見表4。從表4可以看出，隨著C上的H逐漸被－NH2、－NF2或－ONO2取代，衍生物的生成焓越來越小；而被－N3取代時，衍生物的生成焓隨著－N3基團數的增多而增大；當取代基為－NO2時，隨著－NO2基團數的增加，衍生物生成焓呈先減小后增大的規律。對－NH2、－N3和－ONO2三種基團衍生物，比較P-Mono（C3）＋P-Mono（C5）和P-Di（C，C）的值以及P-Mono（C3）＋P-Mono（C5）＋PMono（N）和P-Tri（C，C，N）的值，發現1 H-1，2，4-三唑引入－NH2，－N3和－ONO2基團后其衍生物的生成焓有較好的基團加和性。而當取代基為－NF2基團時由于位阻效應，導致衍生物生成焓不符合基團加和性。當引入－NO2基團時，既有位阻效應又能生成弱氫鍵，同時影響衍生物的生成焓。

表4 1H-1，2，4-三唑引入取代基數目與生成焓的關系Table 4 Relationship of heats of formation with the substituents number join to 1H-1，2，4-triazole

2.2 能量性能

采用Kamlet-Jacobs公式估算爆速（D）和爆壓（p）。從Kamlet-Jacobs公式可以看出，爆速、爆壓與密度和爆熱有關。運用Monte-Carlo方法由0. 001e／Bohr3等電子密度面所包圍的體積求得理論密度，通過修正公式（3）求得其晶體密度，根據得到的固相生成焓計算其爆熱。衍生物的密度、爆熱、氧平衡、爆速和爆壓計算值見表5。

表5可以顯示，取代基的種類、位置以及基團數量對衍生物的密度有一定的影響。除A1和D1外，衍生物的密度都大于1 H-1，2，4-三唑的密度，基團對密度增大的貢獻順序為：－NF2＞－ONO2＞－NO2＞－N3＞－NH2。同一個基團發生C－取代比發生N－取代得到的衍生物的密度大，而在發生C－取代時取代位置對密度的影響卻不大。衍生物密度隨著取代基數量的增加而增大，這種增大的趨勢隨著取代基數量的不斷增大而變緩。

從表5可以看出，－NF2、－ONO2、－NO2和－N3基團發生C－取代或N－取代都能提高1H-1，2，4-三唑衍生物的爆熱。而－NH2基團發生N－取代時1H-1，2，4-三唑衍生物的爆熱增加，但發生C－取代爆熱降低。

表5 1H-1，2，4-三唑衍生物的密度、爆熱、氧平衡、爆速和爆壓計算值Table 5 Predicted densities，heats of detonation，oxygen balance，detonation velocities and detonation pressures of 1H-1，2，4-triazole derivatives

續表5

從表5可以看出，密度和爆熱受取代基團類型、取代位置和取代數量的影響而變化，根據Kamlet-Jacobs公式，密度和爆熱的變化決定了爆速和爆壓的變化。因此，爆速和爆壓也與取代基的類型、取代位置和取代數量有關。－NF2、－NO2、－NO2和－N3基團不但能提高衍生物的密度，同時還能提高其爆熱，因此這些基團的引入可大幅度提高1 H-1，2，-4-三唑衍生物的爆轟性能。

由表5可以看出，C3，C5，E3和E5的爆轟性能高于RDX，但是只有C3和E3的爆轟性能優于HMX。

2.3 熱穩定性

在計算過程中根據Mulliken鍵級選擇相對較弱的鍵，由式（6）計算其化學鍵解離能（BDE）。表6列出1 H-1，2，4-三唑衍生物中弱鍵的解離能和Mulliken鍵級。對于發生C－取代時（A、B和C系列），含－NH2取代基時環內NI－N2鍵的BDE小于C－R鍵，說明衍生物開環反應是誘導反應；相反，含－NO2和－N3基團的衍生物環內N1－N2鍵的BDE大于C3－R鍵，說明在衍生物中取代基的脫落是誘導反應；含－NF2和－ONO2基團的衍生物，N－F和O－NO2鍵的斷裂是解離過程中最容易發生的。而發生N－取代（D和E系列），并且只有－ONO2基團引入到N1位置時，O－NO2鍵最弱；當－NH2、－N3、－NF2或－NO2引入到N1位置時，N1－R鍵的BDE遠小于其他鍵，說明1H-1，2，4-三唑衍生物中只要發生N－取代，化合物中最可能發生斷裂的鍵是N1－R鍵，而－ONO2基團引入到N1位置時，最容易發生斷裂的鍵則是ONO2鍵。

Chung［23］和Su［24］等提出，作為新的高能材料候選物的首要條件是其鍵解離能應超過80k J／mol。從表5數據可知，所有目標化合物的鍵解離能高于80 kJ／mol，但在實際應用中有待進一步實驗驗證。

表6 1H-1，2，4-三唑衍生物中弱鍵的解離能和Mulliken鍵級Table 6 Bond dissociation energies and Mulliken bond orders of the relatively weak bonds in the 1H-1，2，4-triazole derivatives

續表6

3 結 論

（1）運用DFT-B3LYP方法研究了1H-1，2，4-三唑衍生物的生成焓、爆轟性能和熱穩定性。取代基的種類和位置影響衍生物的生成焓，引入－N3基團能顯著增加衍生物的生成焓。對于同類取代基發生N－取代衍生物生成焓值高于發生C－取代衍生物。

（2）－NF2、－ONO2、－NO2和－N3基團不僅能提高衍生物的密度，還能提高其爆熱，有效提高1H-1，2，4-三唑衍生物的爆轟性能。所有衍生物的熱穩定性較好。綜合考慮爆轟性能和熱穩定性，C3、C5、E3和E5有可能成為潛在的高能量密度材料。

［1］ 陽世清，岳守體.國外四嗪四唑類高氮含能材料研究進展［J］.含能材料，2003，11（4）：231-235.

YANG Shi-qing，YUE Shou-ti.Progress in high-nitrogen energetic materials derails derived from tetrazine and tetrazole［J］.Chinese Journal of Energetic Materials，2003，11（4）：231-235.

［2］ Krishna K V，Keresztury G，Sundius T，et al.Hydrogen bonding and molecular vibrations of 3，5-diamino-1，2，4-triazole［J］.Spectrochim Acta A，2005，61：261-267.

［3］ Xu L N，Xiao H M，Fang G Y，et al.Theoretical study on intermolecular interactions of 3-nitro-1，2， 4triazole-5-one-dimers［J］.Acta Chimica Sinica，2005，63（12）：1062-1068.

［4］ Sikder A K，Geetha M，Sarwade D B，et al.Studies on characterization and thermal behaviour of 3-amino-5-nitro-1，2，4-triazole and its derivatives［J］.Hazard，2001，82（1）：1-12.

［5］ 朱佳平，劉勝楠，李永祥.幾種唑類炸藥分子的結構和爆轟性能的密度泛函理論研究［J］.中北大學學報（自然科學版），2011，2：200-206.

ZHU Jia-ping，LIU Sheng-nan，LI Yong-xiang. Structures and detonation performance for several explosive molecule of thiiazoles by density functional theory［J］.Chinese Journal of North University of China（Natural Science Edition），2011，2：200-206.

［6］ Lu Y L，Gong X D，Ju X H.structures and properties of 1，2，3-triazoles and 1，2，4-triazoles［J］.Chinese J Struct Chem，2006，25（5）：582-588.

［7］ Frisch M J，Trucks G W，Schlegel H B.et al.Gaussian 03，Revision B.03.［M］.Pittsburgh：Gaussian Inc，2003.

［8］ Lee C，Yang W，Parr R G.Development of the collsalvetti correlation-energy formula into a functional of the electron density［J］.Phys Rev B，1988，37：785-789.

［9］ Becke A D.Density-functional thermochemistry.Ⅲ. the role of exact exchange［J］.Chem Phys，1993，98：5648-5652.

［10］Zhang X W，Zhu W H，Xiao H M.Comparative theoretical studies of energetic substituted carbon-and ni-trogen-bridged difurazans［J］.Phys Chem A，2010，114：603-612.

［11］Ju X H，Wang X，Bei F L.Substituent effects on heats of formation，group interactions，and detonation properties of polyazidocubanes［J］.Comput Chem，2005，26（12）：1263-1269.

［12］Politzer P，Lane P，Murray J S.Computational characterization of a potential energetic compound：1，3，5，7-Tetranitro-2，4，6，8-tetraazacubane［J］.Cent Eur J Energ Mater，2011，8：39-52.

［13］Politzer P，Murray J S.Some perspectives on estimating detonation properties of C，H，N，O compounds［J］.Cent Eur J Energ Mater，2011，8：209-220.

［14］Politzer P，Murray J S，Grice M E，et al.Calculation of heats of sublimation and solid phase heats of formation［J］.Mol Phys，1997，91：923-928.

［15］Byrd E F C，Rice B M.Improved prediction of heats of formation of energetic materials using quantum chemical calculations［J］.Phys Chem A，2006，110：1005-1013.

［16］Politzer P，Martinez J，Murray JS.et al.An electrostatic interaction correction for improved crystal density predictions［J］.Mol Phys，2009，107：2095-2101.

［17］Kamlet M J，Jacobs S J.Chemistry of detonation.I.A simple method for calculating detonation properties of C，H，N，O explosives［J］.Chem Phys，1968，48：23-35.

［18］Kamlet M J，Ablard J E.Chemistry of detonations. II.buffered equilibria［J］.Chem Phys，1968，48：36-42.

［19］Dean J A.Lange′s Handbook of Chemistry［M］.15th Ed.New York：McGraw-Hill，1999.

［20］David R L.Handbook of chemistry and physics［M］. 84th edn.Boca Raton：CRC Press，2003-2004.

［21］Afeefy H Y，Liebman J F，Stein S E.Neutral Thermochemical Data［M］.Gaithersburg：National Institute of Standards and Technology，2000.

［22］Talawar M B，Sivabalan R，Mukundan T，et al.Environmentally compatible next generation green energetic materials（GEMs）［J］.Hazard Mater，2009，161：589-607.

［23］Chung G S，Schimidt M W，Gordon M S.An ab initio study of potential energy surfaces for N8 isomers［J］. Phys Chem A，2000，104：5647-5650.

［24］Su X F，Cheng X L，Meng C M.et al.Quantum chemical study on nitroimidazole，polynitroimidazole and their methyl derivatives［J］.Hazard Mater，2009，161（1）：551-558.

Theoretical Calculation on Detonation Performance of 1H-1，2，4-Triazole Derivatives

MA Zhong-liang1，ZHOU Hua1，WANG Jian-long1，LI Jing-man2
（1.School of Chemical Engineering and Environment，North University of China，Taiyuan 030051，China；2.Sichuan Tongda Chemical Industry Co.，Ltd，Dazhou Sichuan 635000，China）

The geometries of 1H-1，2，4-triazole derivatives were obtained by density functional theory at DFTB3LYP／6-31G（d，p）level.The heats of formation（HOF）were predicted by designed isodesmic reactions.The average molar volume（V）and theoretical density（ρ）were estimated based on 0.001 e·bohr－3density space.The detonation velocity（D）and pressure（p）of 1H-1，2，4-triazole derivatives were estimated by using the Kamlet-Jacbos equation.It was found that N-substituent in 1H-1，2，4-triazole can improve the HOF values of the derivatives，and－N3is an effective structural unit for improving the HOF values of the derivatives.The－NF2or－ONO2group is very useful for enhancing the detonation performance of these derivatives.The analysis of the bond dissociation energies for several relatively weak bonds suggest that most of the derivatives have good thermal stability.

organic chemistry；quantum chemistry；density functional theory；1H-1，2，4-triazole derivatives；heat of formation；detonation performance；thermal stability

TJ55；O<641 文獻標志碼：A class="emphasis_bold">641 文獻標志碼：A 文章編號：1007-7812（2014）04-0038-07641 文獻標志碼：A

1007-7812（2014）04-0038-07

A 文章編號：1007-7812（2014）04-0038-07

2014-03-10；

2014-05-30

馬忠亮（1967－），男，教授，博導，從事新型發射藥及裝藥研究。