- Research Article
- Open Access
Thermogravimetric analysis, kinetic study, and pyrolysis–GC/MS analysis of 1,1ʹ-azobis-1,2,3-triazole and 4,4ʹ-azobis-1,2,4-triazole
Chemistry Central Journalvolume 12, Article number: 22 (2018)
In general, the greater the number of directly linked nitrogen atoms in a molecule, the better its energetic performance, while the stability will be accordingly lower. But 1,1ʹ-azobis-1,2,3-triazole (1) and 4,4ʹ-azobis-1,2,4-triazole (2) show remarkable properties, such as high enthalpies of formation, high melting points, and relatively high stabilities. In order to rationalize this unexpected behavior of the two compounds, it is necessary to study their thermal decompositions and pyrolyses. Although a great deal of research has been focused on the synthesis and characterization of energetic materials with 1 and 2 as the backbone, a complete report on their fundamental thermodynamic parameters and thermal decomposition properties has not been published.
Thermogravimetric–differential scanning calorimetry were used to obtain the thermal decomposition data of the title compounds. Kissinger and Ozawa–Doyle methods, the two selected non-isothermal methods, are presented for analysis of the solid-state kinetic data. Pyrolysis–gas chromatography/mass spectrometry was used to study the pyrolysis process of the title compounds.
The DSC curves show that the thermal decompositions of 1 and 2 are at different heating rates involved a single exothermic process. The TG curves provide insight into the total weight losses from the compounds associated with this process. At different pyrolysis temperatures, the compositions and types of the pyrolysis products differ greatly and the pyrolysis reaction at 500 °C is more thorough than 400 °C.
Apparent activation energies (E) and pre-exponential factors (lnA/s−1) are 291.4 kJ mol−1 and 75.53 for 1; 396.2 kJ mol−1 and 80.98 for 2 (Kissinger). The values of E are 284.5 kJ mol−1 for 1 and 386.1 kJ mol−1 for 2 (Ozawa–Doyle). The critical temperature of thermal explosion (T b ) is evaluated as 187.01 °C for 1 and 282.78 °C for 2. The title compounds were broken into small fragment ions under the pyrolysis conditions, which then might undergo a multitude of collisions and numerous other reactions, resulting in the formation of C2N2 (m/z 52), etc., before being analyzed by the GC/MS system.
Triazoles are a class of typical nitrogen-rich compounds, which have been widely used in novel energetic materials, medicine, catalysis, and other fields [1,2,3,4]. Over several decades, many studies have shown that azobis-triazole compounds have good energetic properties. This is due to their structure with multiple nitrogen atoms linked directly, along with many C–N and N–N single or double bonds in the molecule, which have good energetic properties [5,6,7,8,9,10,11]. Indeed, compared with a single triazole ring, energetic properties such as energy density, heat of formation, detonation velocity, and detonation pressure, can be significantly improved [8, 12,13,14,15,16].
In general, the greater the number of directly linked nitrogen atoms in a molecule, the better its energetic performance, but stability will be accordingly lower [5, 12, 13, 17,18,19,20,21,22]. The title compounds 1,1ʹ-azobis-1,2,3-triazole (1) and 4,4ʹ-azobis-1,2,4-triazole (2) show remarkable properties, such as high enthalpies of formation, high melting points, and relatively high stabilities [5, 13, 14, 17]. In order to rationalize this unexpected behavior of the two compounds, it is necessary to study their thermal decompositions and pyrolyses. Although a great deal of research has been focused on the synthesis and characterization of energetic materials with 1 and 2 as the backbone [5, 13, 17, 19], a complete report on their fundamental thermodynamic parameters and thermal decomposition properties has not been published.
The stability of the title compounds can be quantitatively determined by studying their thermodynamic properties, such as apparent activation energy (E) and pre-exponential factor (A) of thermal decomposition, and the critical temperature of thermal explosion (T b ).
There are many methods for analyzing non-isothermal solid-state kinetic data from TG and DSC [23,24,25,26]. These methods can be divided into two types: model-fitting and model-free methods, as summarized in Table 1 . Model-fitting methods have been widely used for solid-state reactions because of their ability to directly determine the kinetic parameters from a single TG measurement. However, these methods suffer from several shortcomings, such as their inability to uniquely determine the reaction model, and that application of model-fitting methods to non-isothermal data gives higher values of the kinetic parameters. Conversely, model-free methods require several kinetic curves to perform the analysis. Calculations of several curves at different heating rates are performed on the same value of conversion, which allows calculation of the activation energy at each conversion point. Results from model-free methods tend to be more reliable and reasonable than those from model-fitting methods, especially when dealing with non-thermal data .
Thermal pyrolysis refers to chemical decomposition caused by heat, when the thermal energy applied to the sample exceeding the chemical bond energy of the molecules. PY–GC/MS has been widely used to investigate decomposition processes and pyrolysis products . Analytical pyrolysis allows the thermal breakdown of a molecule into smaller fragments, which are then selected and analyzed by the GC/MS system, providing insight into the decomposition of the sample. The PY–GC/MS can study the thermal pyrolysis of an energetic material and provides information about the nature of the explosion reaction products by GC/MS analysis, which give the information about the reaction process. These in turn can be used to evaluate whether the used explosives are environmental friendly from the identification of the explosive residue species. In addition, the identification of the explosive residue species can also be used in military and counterterrorism practice.
Compound 1 and 2 can be prepared by oxidation of the N–NH2 moieties in 1-amino-1,2,3-triazole and 4-amino-1,2,4-triazole, respectively, with sodium dichloroisocyanurate (SDIC). In this study, the thermal decomposition processes of 1 and 2 have been investigated by dynamic TG–DSC under nitrogen atmosphere at different heating rates, and by PY–GC/MS under helium atmosphere at set temperatures of 400 and 500 °C, respectively. Kinetic parameters have been obtained by two model-free methods, the Kissinger method and the Ozawa–Doyle method combined with a kinetic compensation effect. The data reported herein are expected to be of broad interest to researchers engaged in the study and applications of 1 and 2.
Materials and methodology
1,1ʹ-Azobis-1,2,3-triazole (1) and 4,4ʹ-azobis-1,2,4-triazole (2) were prepared according to literature procedures [13, 14]. The sample purities were > 0.99 (w/w). IR spectra were recorded from solid samples in KBr pellets on a Bruker Tensor 27 spectrometer. Elemental analyses were consistent with the theoretical compositions. The samples were further purified by drying in vacuo at 60 °C for approximately 4 h. The structures of the compounds are shown in Fig. 1.
Synthesis of 1,1ʹ-azobis-1,2,3-triazole (1) 
1-Amino-1,2,3-triazole (1.26 g, 15 mmol) was dissolved in CH3CN (40 mL). The solution was cooled to – 5 to 0 °C and a solution of SDIC (3.33 g, 15 mmol) in water (10 mL) and CH3COOH (5 mL) was added dropwise. The reaction mixture was further stirred at 0 °C for 30 min. It was then neutralized with NaHCO3 and filtered. The filtrate was concentrated to afford the product, which was obtained as a slightly yellow solid after recrystallization from acetone. 1H NMR (400 MHz, [D6]DMSO, 25 °C, TMS): δ = 9.17 (2H, d), 8.21 ppm (2H, d); 13C NMR (100 MHz, [D6]DMSO, 25 °C, TMS): δ = 134.8, 118.0 ppm; IR (KBr): ν = 3128, 1625, 1482, 1320, 1225, 1171 cm−1; MS: m/z 165.0 [M + H]+; elemental analysis calcd. (%) for C4H4N8 (164): C, 29.27; H, 2.44; N, 68.29; found: C, 29.45; H, 2.32; N, 68.23.
Synthesis of 4,4ʹ-azobis-1,2,4-triazole (2) 
Acetic acid (5 mL) was added to a solution of SDIC (5.09 g, 23 mmol) in water (40 mL) with vigorous stirring at 30 °C. After 1 h, the mixture was cooled to 5 °C and a solution of 4-amino-1,2,4-triazole (2.07 g, 25 mmol) in water (10 mL) was added. The reaction mixture was vigorously stirred at 15 °C for 1 h. It was then cooled, and the precipitate that formed was collected by filtration and washed with water at 60 °C. After drying in vacuo, a white product was obtained. 1H NMR (400 MHz; D2O, 25 °C, TMS): δ = 7.22 ppm (4H, s); 13C NMR (100 MHz; D2O, 25 °C, TMS): δ = 138.57 ppm; IR (KBr): ν = 3114, 1493, 1368, 1317, 1180 cm−1; MS: m/z 164 (M+); elemental analysis calcd. (%) for C4H4N8 (164): C, 29.27; H, 2.44; N, 68.29; found: C, 29.32; H, 2.48; N, 68.2.
The thermal decompositions of 1 and 2 under flowing N2 were investigated using a thermogravimetric analyzer (Netzsch STA 449 C; Selb, Germany) and a differential scanning calorimeter (DSC Q2000; New Castle, USA). The TG conditions were as follows: sample mass, ca. 1.0 mg; heating rates, 5, 7, 10, 13, and 20 °C min−1 for 1, and 5, 7, 10, 15, and 20 °C min−1 for 2; atmosphere, N2 (flow rate 30 mL min−1); temperature range, 20–400 °C. The DSC conditions were as follows: heating rates, 5, 7, 10, 13, and 20 °C min−1 for 1, and 5, 7, 10, 15, and 20 °C min−1 for 2; atmosphere, N2 (flow rate 30 mL min−1); temperature range, 20–400 °C. All TG–DSC data were analyzed using Proteus Analysis software.
Thermal pyrolyses of 1 and 2 were investigated using a pyrolysis (EGA/PY-3030D, Fukushima-ken, Japan) gas chromatography–mass spectrometry (QP2010-Ultra, Kyoto, Japan) instrument (PY–GC/MS). The PY conditions were as follows: pyrolysis temperatures, 400 and 500 °C; pyrolysis time, 1 min; injection port temperature, 300 °C. The GC/MS conditions were as follows: capillary chromatographic column, ZB-5HT (30 mm × 0.25 mm × 0.25 μm); heating program: 50 °C, holding for 3 min, then heating to 300 °C at a rate of 10 K min−1, holding to 5 min; injection port temperature, 300 °C; split injection; split ratio 100:1; carrier gas (high-purity helium, 99.9999%) flow rate 1.0 mL min−1; collector temperature, 280 °C; ion source temperature, 250 °C; ion source scan mode, full scan (40–100 m/z).
By using these methods, the kinetic parameters of a solid-state reaction can be obtained without knowing the reaction mechanism.
In previous work [29,30,31,32], Kissinger method was widely used to determine the activation energies with the reaction process that occur under linear heating rate conditions. Although this method has some limitations, it is acceptable when an isoconversional method was used to back up the veracity of the Kissinger method . In addition, since the decomposition reaction process of the title compounds would be very complex, the values of E α of the title compounds obtained by the two isoconversional methods namely Kissinger–Akahira–Sunose  and Starink  method vary greatly and are disorder in the given range of α as 0.05–0.95. And many other types of kinetic methods  were tried to back up the veracity of the Kissinger method, but all of the results are unsatisfactorily. According to the literature [29, 36,37,38,39], an applicable method namely Ozawa–Doyle method was commonly used to back up the Kissinger method in the kinetic calculation of the energetic materials for its’ acceptable result.
In 1957, Kissinger  first introduced a model-free non-isothermal method, which can help researchers to evaluate kinetic parameters without the need of calculating E for each conversion value of the solid-state reaction. This method is described as follows:
In this equation, T P is the peak temperature of the DSC curve. The apparent activation energy (E) and pre-exponential factor (A) can be obtained from the slope −/(RT P ) and intercept ln(AR/E) respectively, of an ln (β/T P 2 ) versus 1/T P plot.
The Ozawa–Doyle method [41, 42] is simple and applicable to reactions that cannot be analyzed by other methods. It has been widely used to determine the apparent activation energy (E) alongside the Kissinger method. The Ozawa–Doyle equation is as follows:
where T P is the peak temperature of the DSC curve and C is a constant. The apparent activation energy (E) can be obtained from a plot of logarithm of heating rates. A plot of logβ versus 1/T P expresses a linear function with an intercept of 0.4567E/R. The calculated E of this method is independent of the mechanism of thermal decomposition.
Results and discussion
The TG and DSC curves of 1 and 2 at different heating rates under N2 atmosphere are shown in Figs. 2 and 3, respectively. The DSC curves show that the thermal decompositions of 1 and 2 at different heating rates involve a single exothermic process. In this process, the molecules of 1 and 2 were broken into smaller pieces by cleavage of the N=N bond linking the two triazole rings and opening of their triazole rings. At the same time, the TG curves provided insight into the total weight losses of the compounds associated with this process.
Characteristic temperatures at different heating rates in the TG–DSC curves of 1 and 2 are shown in Table 2. From Table 2, it can be seen that at 5 °C min−1, thermal decompositions of 1 and 2 started at 176.1 and 281.7 °C, respectively. At higher heating rates, the initial temperature (T 0 ), the extrapolated onset temperature (T e ), and the peak temperature (T P ) of the DSC curves shifted from 176.1, 188.5, and 189.7 °C at 5 °C min−1 to 189.8, 199.7, and 198.0 °C at 20 °C min−1, respectively, for 1. For 2, T 0 , T e , and T P shifted from 281.7, 288.8, and 312.7 °C at 5 °C min−1 to 302.1, 315.1, and 322.7 °C at 20 °C min−1, respectively. These data show that with increased heating rate, the values of T 0 , T e , and T P increase. This behavior can be attributed to heat transfer between the sample and the instrument.
From the thermogravimetric analysis results, we can calculate the kinetic parameters according to the model-free methods. The activation energy (E) and pre-exponential factor (A) were obtained using the Kissinger and Ozawa–Doyle methods.
From the original data of the exothermic peak temperature measured at five different heating rates of 5, 7, 10, 13, and 20 °C min−1 for 1, and 5, 7, 10, 15, and 20 °C min−1 for 2, the apparent activation energies E K and E O , the pre-exponential factors A K , and the linear coefficients r K and r O were determined, as shown in Table 2.
From Table 2, it can be seen that the apparent activation energies (E) for 1 and 2 obtained by the Kissinger method are very close to the values obtained by the Ozawa–Doyle method. The minor differences are assumed to stem from limitations of the method itself and errors in calculation. Moreover, the absolute values of the linear correlation coefficients (r) for 1 and 2 in Table 2 are close to 1, which indicates that the kinetic parameters were obtained with high accuracy.
The Arrhenius equations of the title compounds can be expressed as follows (E is the average of E K and E O ):
where b, c, and d are coefficients.
where E O is the value of E obtained by the Ozawa–Doyle method.
Evidently, the values of apparent activation energy (E), extrapolated onset temperature (T e ), and critical temperature of thermal explosion (T b ) for 2 are consistently higher than those for 1, indicating greater thermodynamic stability of the former. Comparing the values of T b with those for other common energetic compounds: CL-20 (202.07 °C), HMX (267.43 °C), RDX (209.32 °C), NTO (265.53 °C), ENTO (227.44 °C), KNTO (226.32 °C) , ZTO (282.21 °C), ATO (299.64 °C), GZTO·H2O (237.74 °C) , and KZTO·H2O (275.08 °C) , the thermodynamic stability sequence of these compounds can be expressed as: 1 < CL-20 < RDX < KNTO ≈ ENTO < GZTO·H2O < NTO ≈ HMX < KZTO·H2O < ZTO ≈ 2 < ATO.
Thermal pyrolysis analysis
Pyrolysis–gas chromatography–mass spectrometry (PY–GC/MS) can be used to qualitatively analyze pyrolysis products. The pyrolysis chamber was heated to the preset temperature, then the sample was added, after 3 s, a fast heating process was performed and the pyrolysis process was carried out. Fragments were then separated by the GC column and their structures were identified by the MS system. For unknown compounds, one can obtain important information, such as their composition, microstructure, and so on. For known compounds, one can determine the pyrolysis products, and thereby infer the pyrolysis reaction pathways of the compound. This method has many advantages, including a very small injection volume, suitability for a broad range of samples, rapid analysis, and good reproducibility.
Because the molecules of 1 and 2 contain many C–N/N–N single and double bonds, their critical temperatures of thermal explosion (T b ) are below 300 °C. Hence, the explosion reaction must happen during the pyrolysis process, and the actual temperature at the reaction center may briefly reach thousands of degrees Celsius. In the pyrolysis of 1 and 2, following the initial decomposition of the reactive molecule, the pyrolysis products could undergo a multitude of collisions and numerous other reaction processes prior to collection and analysis by the GC/MS system. Hence, investigation of the rapid explosion reaction is a very difficult task, and the mechanistic interpretation inferred by us in this work should be placed in the context of the difficulty in truly isolating microscopic pathways of the explosion reaction.
The pyrolysis–mass spectrometric traces of 1 and 2 at 400 and 500 °C are shown in Figs. 4, 5, 6, 7, respectively. From these figures, it can be seen that for different pyrolysis temperatures, the compositions and types of the pyrolysis products differ greatly. At 400 °C, the retention times of the first three pyrolysis products into the chromatogram were similar to those at 500 °C, but the number of pyrolysis products was fewer at 500 °C, indicating that the pyrolysis reaction at 500 °C is more thorough than 400 °C.
We could estimate the relative contents of the pyrolysis products from a qualitative comparison of the masses of all components using the peak area normalization method. The structures and the relative contents of the pyrolysis products from 1 and 2 are shown in Tables 3 and 4, respectively.
From Table 3, it is clear that with increasing pyrolysis temperature, the numbers of different pyrolysis products were significantly reduced, and the structures were also much simpler. At 400 °C, there were about 24 species among the pyrolysis products of 1, and the differences between the proportions of the different types of pyrolysis products were relatively large. From the pyrolysis results, it is clear that it is difficult to create all of the fragments through direct cleavage of the starting molecules. Hence, we assume that in the explosion reaction, numerous molecules extensively broken into small fragment ions, as shown in Fig. 8, and these fragments underwent secondary reactions, such as coupling, rearrangement, addition, and elimination of hydrogen atoms, prior to detection by the GC/MS system.
The earliest species observed in the pyrolysis at 400 °C were those with m/z 52, m/z 51, and m/z 53, and these species were the sole products at 500 °C. Among these, that at m/z 52 may have arisen from coupling of two m/z 27 (CH=N) fragments by C–C single-bond formed with the elimination of two hydrogen atoms. In the same way, the fragments with m/z 51 and m/z 53 could have been created from that of m/z 26 (CH=CH) coupling with that of m/z 27 (CH=N) through the formation of a C–C single bond, the former with the concurrent elimination reaction of two hydrogen atoms. The larger molecular mass of 103 may have been created by hydrogenation of the fragment with m/z 96. All of the other fragment ions can reasonably be formed from the small fragment ions (as shown in Fig. 8) originating from the explosion reaction of 1, through coupling and rearrangement reactions, sometimes with the concurrent addition or elimination of hydrogen atoms, as shown in Fig. 9.
From Table 4, the same conclusion can be reached. In the explosion reaction process of 2, numerous molecules were extensively broken into much smaller fragment ions, as shown in Fig. 10. In the same way, these fragments could undergo a multitude of secondary reactions prior to detection by the GC/MS system.
The principal pyrolysis products of 2 were also those with m/z 52, m/z 51, and m/z 53, where that at m/z 52 can be created by coupling of two fragments with m/z 27 (CH=N) through C–C single-bond formation with the elimination of two hydrogen atoms. In the same way, those at m/z 51 and m/z 53 can be created by coupling of two fragments with m/z 13 (CH) through C=C double-bond formation creating a fragment with m/z 26 (HC=CH), followed by coupling with a fragment of m/z 27 (CH=N) through formation of a C–C single bond, the former with the concurrent elimination of two hydrogen atoms. The fragment with m/z 81 may have been created as a trimer of that with m/z 27 (CH=N), and that at m/z 141 may have been created by a fragment with m/z 96 coupling with a fragment with m/z 41 (N–N=CH) and the addition of hydrogen atoms, as shown in Fig. 11.
Nevertheless, some fragments are yet unaccounted for, such as those with m/z 77 (Table 3, entries 8, 9, 11), m/z 154 (entries 21, 22), m/z 256 (entry 23), and m/z 175 (entry 24) for 1, and m/z 77 (Table 4, entry 8) for 2.
Experimental kinetic studies on the thermal decomposition processes of two typical nitrogen-rich energetic materials (1 and 2) were described, in which kinetic parameters, namely the apparent activation energy (E) and pre-exponential factor (lnA), were determined by the Kissinger and Ozawa–Doyle methods.
By the Kissinger method, values of E as 291.4 and 396.2 kJ mol−1 were obtained for 1 and 2, respectively, with lnA(s−1) values of 75.53 and 80.98, respectively. By the Ozawa–Doyle method, the values of E were 284.5 and 386.1 kJ mol−1 for 1 and 2, respectively, showing good agreement. The linear correlation coefficients (r) were close to 1, validating the results.
The critical temperatures of thermal explosion (T b ) were determined as 187.01 °C for 1 and 282.78 °C for 2. From the values of E, A, and T b , 2 is clearly more thermodynamically stable than 1. Critical temperatures of thermal explosion follow the sequence: 1 < CL-20 < RDX < KNTO ≈ ENTO < GZTO·H2O < NTO ≈ HMX < KZTO·H2O < ZTO ≈ 2 < ATO.
By PY–GC/MS, thermal pyrolyses of 1 and 2 at 400 and 500 °C generated a greater number of species. By analysis of the possible structures of the pyrolysis products, some conclusions about the pyrolysis pathways of 1 and 2 were drawn. The fragments detected by GC/MS following the pyrolyses of 1 and 2 were likely due to numerous secondary reactions, such as coupling, rearrangement, and addition or elimination of hydrogen atoms, of the smaller ion fragments derived from the explosion reactions (Additional file 1).
Kumar D, Imler GH, Parrish DA, Shreeve JM (2017) A highly stable and insensitive fused triazole–triazine explosive (TTX). Chem Eur J 23:1743–1747
Piercey DG, Chavez DE, Scott BL, Imler GH, Parrish DA (2016) An energetic triazolo-1,2,4-triazine and its N-oxide. Angew Chem Int Ed 55:15315–15318
Wang XS, Huang BS, Liu XY, Zhan P (2016) Discovery of bioactive molecule from CuAAC click-chemistry-based combinatorial libraries. Drug Discov Today 21:118–132
Li PZ, Wang XJ, Liu J, Lin JS, Zou RQ, Zhao YL (2016) A triazole-containing metal-organic framework as a highly effective and substrate size-dependent catalyst for CO2 conversion. J Am Chem Soc 138:2142–2145
Singh RP, Gao HX, Meshri DT, Shreeve JM (2007) Nitrogen-rich heterocycles. Struct Bond 125:35–83
Klapötke TM, Sabaté CM (2008) Bistetrazoles: nitrogen-rich, high-performing, insensitive energetic compounds. Chem Mater 20:3629–3637
Tao GH, Guo Y, Parrish DA, Shreeve JM (2010) Energetic 1,5-diamino-4H-tetrazolium nitro-substituted azolates. J Mater Chem 20:2999–3005
Gao Y, Gao HX, Twamley B, Shreeve JM (2007) Energetic nitrogen rich salts of N,N-bis[1(2)H-tetrazol-5-yl]amine. Adv Mater 19:2884–2888
Joo YH, Shreeve JM (2010) High-density energetic mono- or bis(oxy)-5-nitroiminotetrazoles. Angew Chem 122:7478–7481
Huynh MH, Hiskey MA, Hartline EL, Montoya DP, Gilardi R (2004) Polyazido high-nitrogen compounds: hydrazoand azo-1,3,5-triazine. Angew Chem Int Ed 43:4924–4928
Klapötke TM, Piercey DG (2011) 1,1′-Azobis(tetrazole): a high energetic nitrogen-rich compound with a N10 chain. Inorg Chem 50:2732–2734
Qi C, Li SH, Li YC, Wang Y, Zhao XX, Pang SP (2012) Synthesis and promising properties of a new family of high-nitrogen compounds: polyazido- and polyamino-substituted N,N′-azo-1,2,4-triazoles. Chem Eur J 18(51):16562–16570
Li YC, Qi C, Li SH, Zhang HH, Sun CH, Yu YZ, Pang SP (2010) 1,1′-Azobis-1,2,3-triazole: a high-nitrogen compound with stable N8 structure and photochromism. J Am Chem Soc 132:12172–12173
Liu W, Li SH, Li YC, Yang YZ, Yu Y, Pang SP (2014) Nitrogen-rich salts based on polyamino substituted N,N′-azo-1,2,4-triazole: a new family of high-performance energetic materials. J Mater Chem A 2:15978–15986
Sivabalan R, Anniyappan M, Pawar SJ, Talawar MB, Gore GM, Venugopalan S, Gandhe BR (2006) Synthesis, characterization and thermolysis studies on triazole and tetrazole based high nitrogen content high energy materials. J Hazard Mater A137:672–680
Li ZM, Zhang JG, Cui Y, Zhang TL, Shu YJ, Sinditskii VP, Serushkin VV, Egorshin VY (2010) A novel nitrogen-rich cadmium coordination compound based on 1,5-diaminotetrazole: synthesis, structure investigation, and thermal properties. J Chem Eng Data 55:3109–3116
Qi C, Li SH, Li YC, Wang Y, Chen XK, Pang SP (2011) A novel stable high-nitrogen energetic material: 4,4′-azobis(1,2,4-triazole). J Mater Chem 21:3221
Zhang QH, Shreeve JM (2013) Growing catenated nitrogen atom chains. Angew Chem Int Ed 52:8792–8794
Yin P, Parrish DA, Shreeve JM (2014) Bis(nitroamino-1,2,4-triazolates): N-bridging strategy toward insensitive energetic materials. Angew Chem Int Ed 53:12889–12892
Politzer P, Lane P, Murray JS (2014) Some interesting aspects of N-oxides. Mol Phys 112:719–725
Fabian J, Lewars E (2004) Azabenzenes (azines)—the nitrogen derivatives of benzene with one to six N atoms: stability, homodesmotic stabilization energy, electron distribution, and magnetic ring current; a computational study. Can J Chem 82:50–69
Politzer P, Lane P, Murray JS (2013) Computational analysis of relative stabilities of polyazine N-oxides. Struct Chem 24:1965–1974
Simon P (2004) Isoconversional methods—fundamentals, meaning and application. J Therm Anal Calorim 76:123–132
Fan MH, Panezai H, Sun JH, Bai SY, Wu X (2014) Thermal and kinetic performance of water desorption for N2 absorption in Li-LSX zeolite. J Phys Chem C 118:23761–23767
Li KY, Huang XY, Fleischmann C, Rein G, Ji J (2014) Pyrolysis of medium-density fiberboard: optimized search for kinetics scheme and parameters via a genetic algorithm driven by Kissinger’s method. Energy Fuel 28:6130–6139
Sbirrazzuoli N, Vincent L, Mija A, Guigo N (2009) Integral, differential and advanced isoconversional methods. Complex mechanisms and isothermal predicted conversion-time curves. Chemom Intell Lab Syst 96:219–226
Slopiecka K, Bartocci P, Fantozzi F (2012) Thermogravimetric analysis and kinetic study of poplar wood pyrolysis. Appl Energy 97:491–497
Zhu P, Sui S, Wang B (2004) A study of pyrolysis and pyrolysis products of flame-retardant cotton fabrics by DSC, TGA, and PY–GC–MS. J Anal Appl Pyrolysis 71:645
Wu BD, Wang SW, Yang L, Zhang TL, Zhang JG, Zhou ZN, Yu KB (2011) Preparation, crystal structures, thermal decomposition and explosive properties of two novel energetic compounds M(IMI)4(N3)2 (M = CuII and NiII, IMI = Imidazole): the new high-nitrogen materials (N > 46%). Eur J Inorg Chem 16:2616–2623
Ma C, Huang J, Ma HX, Xu KZ, Lv XQ, Song JR, Zhao NN, He JY, Zhao YS (2013) Preparation, crystal structure, thermal decomposition, quantum chemical calculations on [K(ZTO)·H2O]∞ and its ligand ZTO. J Mol Struct 1036:521–527
Chauhan NP, Mozafari M, Ameta R, Punjabi PB, Ameta SC (2015) Spectral and thermal characterization of halogen-bonded novel crystalline oligo(p-bromoacetophenone formaldehyde). J Phys Chem B 119:3223–3230
Tan CC, Dalapati GK, Tan HR, Bosman MB, Hui HK, Tripathy S, Chi D (2015) Crystallization of sputter-deposited amorphous (FeSi2)1−xAlx thin films. Cryst Growth Des 15:1692–1696
Vyazovkin S, Burnham AK, Criado JM, Pérez-Maqueda LA, Popescu C, Sbirrazzuoli N (2011) ICTAC Kinetics Committee recommendations for performing kinetic computations on thermal analysis data. Thermochim Acta 520:1–19
Akahira T, Sunose T (1971) Method of determining activation deterioration constant of electrical insulating materials. Res Rep Chiba Inst Technol Sci Technol 16:22–31
Starink MJ (2003) The determination of activation energy from linear heating rate experiments: a comparison of the accuracy of isoconversion methods. Thermochim Acta 404:163–176
Chen HY, Zhang TL, Zhang JG, Qiao XJ, Yu KB (2006) Crystal structure, thermal decomposition mechanism and explosive properties of [Na(H2TNPG)(H2O)]n. J Hazard Mater A129:31–36
Li Y, Wu BD, Zhang TL, Liu ZH, Zhang JG (2010) Preparation, crystal structure, thermal decomposition, and explosive properties of [Cd(en)(N3)2]n. Propellants Explos Pyrotech 35:521–528
Wu BD, Li Y, Wang SW, Zhang TL, Zhang JG, Zhou ZN, Yu KB (2011) Preparation, crystal structure, thermal decomposition, and explosive properties of a novel energetic compound [Zn(N2H4)2(N3)2]n: a new high-nitrogen material (N = 65.60%). Z Anorg Allg Chem 637:450–455
Ma C, Huang J, Zhong YT, Xu KZ, Song JR, Zhang Z (2013) Preparation, structural investigation and thermal decomposition behavior of two high-nitrogen energetic materials: ZTO·2H2O and ZTO(phen)·H2O. Bull Korean Chem Soc 34:2086–2092
Kissinger HE (1957) Reaction kinetics in differential thermal analysis. J Anal Chem 29:1702–1706
Doyle CD (1961) Kinetic analysis of thermogravimetric data. J Appl Polym Sci 15:285–292
Ozawa T (1965) A new method of analyzing thermogravimetric data. Bull Chem Soc Jpn 38:1881–1885
Hu RZ, Gao SL, Zhao FQ, Shi QZ, Zhang TL, Zhang JJ (2008) Thermal analysis kinetics, 2nd edn. Science Press, Beijing
Zhang TL, Hu RZ, Xie Y (1994) The estimation of critical temperatures of thermal explosion for energetic materials using non-isothermal DSC. Thermochim Acta 244:171–176
Zhong YT, Huang J, Song JR, Xu KZ, Zhao D, Wang LQ, Zhang XY (2011) Synthesis, crystal structure and thermal behavior of GZTO·H2O. Chin J Chem 29:1672–1676
CJ, YL and SP conceived and designed the experiments. CJ and SZ performed the experiments. CJ, YL, SZ and TF analyzed the data. CJ completed the manuscript. All authors read and approved the final manuscript.
The authors gratefully acknowledge the financial support from the National Natural Science Foundation of China (21576026 and U153062).
The authors declare that they have no competing interests.
Availability of data and materials
The experiment data supporting the conclusions of this article are included within the article and its additional file.
Ethics approval and consent to participate
All authors approval and consent to the publication.
National Natural Science Foundation of China (21576026 and U153062).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
- Thermal decomposition
- Kinetic study
- Thermogravimetric–differential scanning calorimetry (TG–DSC)
- Pyrolysis–gas chromatography/mass spectrometry (PY–GC/MS)