# Kinetic Analysis of the Thermal Decomposition of Polymer-Bonded Explosive Based on PETN: Model-Fitting Method and Isoconversional Method

Authors: Trung Toan Nguyen; Duc Nhan Phan; Van Thom Do; Hoang Nam Nguyen
Journal: Advances in Materials Science and Engineering
 (2020)

---

## Abstract

This work investigates kinetics and thermal decomposition behaviors of pentaerythritol tetranitrate (PETN) and two polymer-bonded explosive (PBX) samples created from PETN (named as PBX-PN-85 and PBX-PP-85) using the vacuum stability test (VST) and thermogravimetry (TG/DTG) techniques. Both model-free (isoconversional) and model-fitting methods were applied to determine the kinetic parameters of the thermal decomposition. It was found that kinetic parameters obtained by the modified Kissinger–Akahira–Sunose method (using non-isothermal TG/DTG data) were close to those obtained by the isoconversional and model-fitting methods that use isothermal VST data. The activation energy values of thermal decomposition reactions were 125.6–137.1, 137.3–144.9, and 143.9–152.4 kJ·mol−1 for PBX-PN-85, PETN, and PBX-PP-85, respectively. The results demonstrate the negative effect of the nitrocellulose-based binder in reducing the thermal stability of single PETN, while the polystyrene-based binder seemingly shows no adverse influence on the thermal decomposition of PETN in our presented PBX compositions.


---

## Body

## 1. Introduction

Pentaerythritol tetranitrate (also called PETN), which is a popular nitrate explosive, is extensively used in civilian and military applications due to its high energy characteristics [1, 2]. The main drawback of PETN is that it is highly sensitive to impact and friction pulses, and thus there may be explosion risks in the processes of production, storage, and especially in the transportation of PETN. Therefore, PETN is often used in the form of a polymer-bonded explosive (also called PBX) [3–5] to overcome this shortcoming. The PBX sample based on PETN is a composite explosive in which the PETN particles are covered and bonded by the polymeric binder. The polymeric binder in PBX composition can reduce the sensitivity of explosive to mechanical pulses and make the charge compression easier and safer [6–9].Like all other high-energy materials, the PBX composition is also an unstable thermodynamic system [10]. Compared to the single explosive such as PETN, the thermal decomposition of PBX is generally more complex because of the incompatible reactions between explosives and the polymeric binders [11–14], and the safety and performance of PBX sample may change during aging. Therefore, the behaviors and kinetics of thermal decomposition of PBX samples should be evaluated and clarified to assess the safety and the thermal hazard potential during storage. The decomposition kinetics and thermal behaviors of several composite explosives based on PETN have been carefully examined, and it was shown that the presence of the additive materials exerts powerful impacts on the thermal decomposition of these explosives [15–19].In this study, we investigate the thermal decomposition kinetics of two PETN-based PBX samples named as PBX-PN-85 and PBX-PP-85. These results were compared with those of a single PETN to evaluate the influence of binders on the thermal decomposition behaviors of explosives. Thermal analyses were conducted by thermogravimetry (TG/DTG) and vacuum stability test (VST), where the VST method (developed by STABIL, Czech Republic) is commonly used to examine the chemical stability, compatibility, and shelf-life of energetic materials in recent years [10, 14, 20]. Based on TG/DTG curves and VST results, we determined kinetic parameters by using isoconversional and model-fitting methods according to the ICTAC kinetics committee recommendations [21]. It is noteworthy that while model-fitting and isoconversional methods are both applicable to isothermal data, the model-fitting method is generally not applicable to non-isothermal techniques since the kinetic parameters obtained for non-isothermal data utilizing model-fitting technique are highly variable [22].


## 2. Materials and Methods

### 2.1. Materials

We imported Class-1 PETN with a melting temperature of over 139.0°C from South Korea. Good quality nitrocellulose (NC), that has a nitrogen content of 12.20%, was supplied from the Vietnam factory. Polystyrene (PS), which has an average molecular weight of 80,000u, was produced in our laboratory. The widely used plasticizer, dioctyl phthalate (DOP), was utilized to make the binder from the PS and NC (the ratios of DOP/NC and DOP/PS were 3/1 and 2/1, respectively). Besides, some suitable solvents such as toluene and ethyl acetate are used in this study.


### 2.2. Experimental Techniques and Methods

#### 2.2.1. Preparation of PBX Samples

PS and NC were plasticized by DOP (the ratios of DOP/NC and DOP/PS are 3/1 and 2/1, respectively), then dissolved in the suitable solvents (toluene for PS and ethyl acetate for NC), maintaining the polymer/solvent ratio of 1/15 (w/v), and kept for approximately 5 hours to obtain a homogeneous solution. The sample preparation method was carried out by mixing PETN crystals with the binders based on PS and NC. The mixing process took place on the heating device at 70°C for 30 minutes, and the PETN crystals were covered by the binder layers [6, 23]. After that, the mixtures were vacuum-dried at 60 to 70°C for 5 hours to remove the solvents. Table 1 shows the compositions of two proposed PBX samples.Table 1
The composition of two PBX samples based on PETN.
CompositionContent of materials (wt%)PETNDOPPSNCPBX-PP85.0010.005.00—PBX-PN85.0011.25—3.75


#### 2.2.2. Experimental Techniques

The thermal decomposition behaviors of PBX samples were determined by employing thermogravimetry (TG/DTG) and vacuum stability test (VST). The experimental conditions are as follows.TG/DTG analysis was conducted using a NETZSCH STA 409 PC/PG (NETZSCH-Gerätebau GmbH, Selb, Germany). PBX samples (around 5.0 mg) were put in an aluminum oxide crucible and heated at various heating rates of 3, 5, 7, and 10 K·min−1. The sample vial was heated from 30 to 350°C. Experimental heating processes were carried on under a dynamic nitrogen atmosphere with a flow rate of 30 mL·min−1.The VST test was carried out in a STABIL apparatus (OZM Research, Pardubice, Czech Republic). During this test, PBX samples were heated isothermally at different temperatures of 135, 140, 145, and 150°C under vacuum pressure (at most 0.672 kPa). A pressure sensor and a computer are used to record the relationship of the released gas volume versus heating time. The sample mass was 20.0 mg, and we performed the test from 5 to 10 days (depending on the temperature test), until the end of the thermal decomposition reaction.


#### 2.2.3. Kinetics Analysis

The kinetics of thermal decomposition reactions in solids can be expressed as the following equation:(1)dφdt=kT⋅cφ,where t is time (s); φ is the extent of conversion; T is the temperature (K); and kT is a function of conversion representing the reaction model. For the TG tests, the extent of conversion φ is defined as follows:(2)φt=mi−mtmi−mf,where mi, mt, and mf are the initial sample mass, sample mass at any time t, and sample mass at the end of the decomposition, respectively. For the VST test, the value of φ is expressed as [13, 24](3)φt=pt−p0pmax−p0,where pmax is the pressure at the end of decomposition, pt is the pressure at any time t, and p0 is initial pressure. Following the Arrhenius equation, kT is computed by(4)kT=A⋅e−E/RT,where R, T, E, and A are the universal gas constant, absolute temperature, activation energy, and pre-exponential factor, respectively.(1) Model-Fitting Method. Model-fitting methods employed different forms of reaction model for finding the referred decomposition kinetic parameters; some of these reaction models are illustrated in Table 2.Table 2
Different reaction models generally applied to describe the thermal decomposition in solids [21, 22, 24, 25].
Model no.Reaction modelcφiφ1Avarami-Erofe’ev (A2)21−φ−ln1−φ1/2[−ln1−φ]1/22Avarami-Erofe’ev (A3)31−φ−ln1−φ2/3[−ln1−φ]1/33Avarami-Erofe’ev (A4)41−φ−ln1−φ3/4[−ln1−φ]1/44Power law (P2)2φ(1/2)φ(1/2)5Power law (P3)3φ(2/3)φ(1/3)6Power law (P4)4φ(3/4)φ(1/4)71-D diffusion (D1)1/(2φ)φ282-D diffusion (D2)−ln1−φ−1[1−φln1−φ] + φ93-D diffusion-Jander (D3)[31−φ2/3]/21−1−φ1/3[1−1−φ1/3]210First-order (F1)1−φ−ln1−φ11Second-order (F2)1−φ21−φ−1−112Third-order (F3)1−φ31−φ−2−1/213Contracting area (R2)21−φ1/21−1−φ1/214Contracting volume (R3)31−φ2/31−1−φ1/3For isothermal conditions, equation (1) can be used in the integral form as follows [22, 24]:(5)iφ=∫0φcφ−1dφ=kTt,where iφ is the integrated form of the reaction model cφ. Replacing each reaction model from Table 1 into equation (5), the rate constant can be computed from the slope of the straight line when plotting iφ against t. Then, we can base on iφ and cφ formulations in the reaction model to determine kT. For each selected reaction model, kT is determined at different temperatures, and the kinetic parameters can be calculated by the Arrhenius equation:(6)lnkT=lnA−ERT.(2) Isoconversional (Model-Free) Method. The reaction model (as in equation (1)) in the isoconversional method is considered as being independent of temperature. Under the isothermal condition of the VST test, combining equations (4) and (5), we obtain [13, 22, 24](7)−lntφ,j=lnAiφ−EφRTj,where the activation energy (Eφ) is estimated from the slope of the straight line when plotting (−lntφ,j) against (1/Tj). Thus, the dependencies of activation energy (Eφ) on the conversion of the thermal decomposition (φ) could be obtained.At a linear heating rateβ=dT/dt and under non-isothermal condition, equation (1) can be described as the following:(8)dφdt=βdφdT=Ae−E/RTcφ,and Kissinger, Akahira, and Sunose introduced the model-free method using the following linear equation [22, 26]:(9)lnβkTφ,k2=lnARiφEφ−EφRTφ.For each extent of conversionφ, the values of E and A can be found from the slope and the intercept when plotting lnβ/T2 against 1/T. The function iφ can be determined according to the methodology reported in the literature [27, 28].


## 2.1. Materials

We imported Class-1 PETN with a melting temperature of over 139.0°C from South Korea. Good quality nitrocellulose (NC), that has a nitrogen content of 12.20%, was supplied from the Vietnam factory. Polystyrene (PS), which has an average molecular weight of 80,000u, was produced in our laboratory. The widely used plasticizer, dioctyl phthalate (DOP), was utilized to make the binder from the PS and NC (the ratios of DOP/NC and DOP/PS were 3/1 and 2/1, respectively). Besides, some suitable solvents such as toluene and ethyl acetate are used in this study.


## 2.2. Experimental Techniques and Methods

### 2.2.1. Preparation of PBX Samples

PS and NC were plasticized by DOP (the ratios of DOP/NC and DOP/PS are 3/1 and 2/1, respectively), then dissolved in the suitable solvents (toluene for PS and ethyl acetate for NC), maintaining the polymer/solvent ratio of 1/15 (w/v), and kept for approximately 5 hours to obtain a homogeneous solution. The sample preparation method was carried out by mixing PETN crystals with the binders based on PS and NC. The mixing process took place on the heating device at 70°C for 30 minutes, and the PETN crystals were covered by the binder layers [6, 23]. After that, the mixtures were vacuum-dried at 60 to 70°C for 5 hours to remove the solvents. Table 1 shows the compositions of two proposed PBX samples.Table 1
The composition of two PBX samples based on PETN.
CompositionContent of materials (wt%)PETNDOPPSNCPBX-PP85.0010.005.00—PBX-PN85.0011.25—3.75


### 2.2.2. Experimental Techniques

The thermal decomposition behaviors of PBX samples were determined by employing thermogravimetry (TG/DTG) and vacuum stability test (VST). The experimental conditions are as follows.TG/DTG analysis was conducted using a NETZSCH STA 409 PC/PG (NETZSCH-Gerätebau GmbH, Selb, Germany). PBX samples (around 5.0 mg) were put in an aluminum oxide crucible and heated at various heating rates of 3, 5, 7, and 10 K·min−1. The sample vial was heated from 30 to 350°C. Experimental heating processes were carried on under a dynamic nitrogen atmosphere with a flow rate of 30 mL·min−1.The VST test was carried out in a STABIL apparatus (OZM Research, Pardubice, Czech Republic). During this test, PBX samples were heated isothermally at different temperatures of 135, 140, 145, and 150°C under vacuum pressure (at most 0.672 kPa). A pressure sensor and a computer are used to record the relationship of the released gas volume versus heating time. The sample mass was 20.0 mg, and we performed the test from 5 to 10 days (depending on the temperature test), until the end of the thermal decomposition reaction.


### 2.2.3. Kinetics Analysis

The kinetics of thermal decomposition reactions in solids can be expressed as the following equation:(1)dφdt=kT⋅cφ,where t is time (s); φ is the extent of conversion; T is the temperature (K); and kT is a function of conversion representing the reaction model. For the TG tests, the extent of conversion φ is defined as follows:(2)φt=mi−mtmi−mf,where mi, mt, and mf are the initial sample mass, sample mass at any time t, and sample mass at the end of the decomposition, respectively. For the VST test, the value of φ is expressed as [13, 24](3)φt=pt−p0pmax−p0,where pmax is the pressure at the end of decomposition, pt is the pressure at any time t, and p0 is initial pressure. Following the Arrhenius equation, kT is computed by(4)kT=A⋅e−E/RT,where R, T, E, and A are the universal gas constant, absolute temperature, activation energy, and pre-exponential factor, respectively.(1) Model-Fitting Method. Model-fitting methods employed different forms of reaction model for finding the referred decomposition kinetic parameters; some of these reaction models are illustrated in Table 2.Table 2
Different reaction models generally applied to describe the thermal decomposition in solids [21, 22, 24, 25].
Model no.Reaction modelcφiφ1Avarami-Erofe’ev (A2)21−φ−ln1−φ1/2[−ln1−φ]1/22Avarami-Erofe’ev (A3)31−φ−ln1−φ2/3[−ln1−φ]1/33Avarami-Erofe’ev (A4)41−φ−ln1−φ3/4[−ln1−φ]1/44Power law (P2)2φ(1/2)φ(1/2)5Power law (P3)3φ(2/3)φ(1/3)6Power law (P4)4φ(3/4)φ(1/4)71-D diffusion (D1)1/(2φ)φ282-D diffusion (D2)−ln1−φ−1[1−φln1−φ] + φ93-D diffusion-Jander (D3)[31−φ2/3]/21−1−φ1/3[1−1−φ1/3]210First-order (F1)1−φ−ln1−φ11Second-order (F2)1−φ21−φ−1−112Third-order (F3)1−φ31−φ−2−1/213Contracting area (R2)21−φ1/21−1−φ1/214Contracting volume (R3)31−φ2/31−1−φ1/3For isothermal conditions, equation (1) can be used in the integral form as follows [22, 24]:(5)iφ=∫0φcφ−1dφ=kTt,where iφ is the integrated form of the reaction model cφ. Replacing each reaction model from Table 1 into equation (5), the rate constant can be computed from the slope of the straight line when plotting iφ against t. Then, we can base on iφ and cφ formulations in the reaction model to determine kT. For each selected reaction model, kT is determined at different temperatures, and the kinetic parameters can be calculated by the Arrhenius equation:(6)lnkT=lnA−ERT.(2) Isoconversional (Model-Free) Method. The reaction model (as in equation (1)) in the isoconversional method is considered as being independent of temperature. Under the isothermal condition of the VST test, combining equations (4) and (5), we obtain [13, 22, 24](7)−lntφ,j=lnAiφ−EφRTj,where the activation energy (Eφ) is estimated from the slope of the straight line when plotting (−lntφ,j) against (1/Tj). Thus, the dependencies of activation energy (Eφ) on the conversion of the thermal decomposition (φ) could be obtained.At a linear heating rateβ=dT/dt and under non-isothermal condition, equation (1) can be described as the following:(8)dφdt=βdφdT=Ae−E/RTcφ,and Kissinger, Akahira, and Sunose introduced the model-free method using the following linear equation [22, 26]:(9)lnβkTφ,k2=lnARiφEφ−EφRTφ.For each extent of conversionφ, the values of E and A can be found from the slope and the intercept when plotting lnβ/T2 against 1/T. The function iφ can be determined according to the methodology reported in the literature [27, 28].


## 2.2.1. Preparation of PBX Samples

PS and NC were plasticized by DOP (the ratios of DOP/NC and DOP/PS are 3/1 and 2/1, respectively), then dissolved in the suitable solvents (toluene for PS and ethyl acetate for NC), maintaining the polymer/solvent ratio of 1/15 (w/v), and kept for approximately 5 hours to obtain a homogeneous solution. The sample preparation method was carried out by mixing PETN crystals with the binders based on PS and NC. The mixing process took place on the heating device at 70°C for 30 minutes, and the PETN crystals were covered by the binder layers [6, 23]. After that, the mixtures were vacuum-dried at 60 to 70°C for 5 hours to remove the solvents. Table 1 shows the compositions of two proposed PBX samples.Table 1
The composition of two PBX samples based on PETN.
CompositionContent of materials (wt%)PETNDOPPSNCPBX-PP85.0010.005.00—PBX-PN85.0011.25—3.75


## 2.2.2. Experimental Techniques

The thermal decomposition behaviors of PBX samples were determined by employing thermogravimetry (TG/DTG) and vacuum stability test (VST). The experimental conditions are as follows.TG/DTG analysis was conducted using a NETZSCH STA 409 PC/PG (NETZSCH-Gerätebau GmbH, Selb, Germany). PBX samples (around 5.0 mg) were put in an aluminum oxide crucible and heated at various heating rates of 3, 5, 7, and 10 K·min−1. The sample vial was heated from 30 to 350°C. Experimental heating processes were carried on under a dynamic nitrogen atmosphere with a flow rate of 30 mL·min−1.The VST test was carried out in a STABIL apparatus (OZM Research, Pardubice, Czech Republic). During this test, PBX samples were heated isothermally at different temperatures of 135, 140, 145, and 150°C under vacuum pressure (at most 0.672 kPa). A pressure sensor and a computer are used to record the relationship of the released gas volume versus heating time. The sample mass was 20.0 mg, and we performed the test from 5 to 10 days (depending on the temperature test), until the end of the thermal decomposition reaction.


## 2.2.3. Kinetics Analysis

The kinetics of thermal decomposition reactions in solids can be expressed as the following equation:(1)dφdt=kT⋅cφ,where t is time (s); φ is the extent of conversion; T is the temperature (K); and kT is a function of conversion representing the reaction model. For the TG tests, the extent of conversion φ is defined as follows:(2)φt=mi−mtmi−mf,where mi, mt, and mf are the initial sample mass, sample mass at any time t, and sample mass at the end of the decomposition, respectively. For the VST test, the value of φ is expressed as [13, 24](3)φt=pt−p0pmax−p0,where pmax is the pressure at the end of decomposition, pt is the pressure at any time t, and p0 is initial pressure. Following the Arrhenius equation, kT is computed by(4)kT=A⋅e−E/RT,where R, T, E, and A are the universal gas constant, absolute temperature, activation energy, and pre-exponential factor, respectively.(1) Model-Fitting Method. Model-fitting methods employed different forms of reaction model for finding the referred decomposition kinetic parameters; some of these reaction models are illustrated in Table 2.Table 2
Different reaction models generally applied to describe the thermal decomposition in solids [21, 22, 24, 25].
Model no.Reaction modelcφiφ1Avarami-Erofe’ev (A2)21−φ−ln1−φ1/2[−ln1−φ]1/22Avarami-Erofe’ev (A3)31−φ−ln1−φ2/3[−ln1−φ]1/33Avarami-Erofe’ev (A4)41−φ−ln1−φ3/4[−ln1−φ]1/44Power law (P2)2φ(1/2)φ(1/2)5Power law (P3)3φ(2/3)φ(1/3)6Power law (P4)4φ(3/4)φ(1/4)71-D diffusion (D1)1/(2φ)φ282-D diffusion (D2)−ln1−φ−1[1−φln1−φ] + φ93-D diffusion-Jander (D3)[31−φ2/3]/21−1−φ1/3[1−1−φ1/3]210First-order (F1)1−φ−ln1−φ11Second-order (F2)1−φ21−φ−1−112Third-order (F3)1−φ31−φ−2−1/213Contracting area (R2)21−φ1/21−1−φ1/214Contracting volume (R3)31−φ2/31−1−φ1/3For isothermal conditions, equation (1) can be used in the integral form as follows [22, 24]:(5)iφ=∫0φcφ−1dφ=kTt,where iφ is the integrated form of the reaction model cφ. Replacing each reaction model from Table 1 into equation (5), the rate constant can be computed from the slope of the straight line when plotting iφ against t. Then, we can base on iφ and cφ formulations in the reaction model to determine kT. For each selected reaction model, kT is determined at different temperatures, and the kinetic parameters can be calculated by the Arrhenius equation:(6)lnkT=lnA−ERT.(2) Isoconversional (Model-Free) Method. The reaction model (as in equation (1)) in the isoconversional method is considered as being independent of temperature. Under the isothermal condition of the VST test, combining equations (4) and (5), we obtain [13, 22, 24](7)−lntφ,j=lnAiφ−EφRTj,where the activation energy (Eφ) is estimated from the slope of the straight line when plotting (−lntφ,j) against (1/Tj). Thus, the dependencies of activation energy (Eφ) on the conversion of the thermal decomposition (φ) could be obtained.At a linear heating rateβ=dT/dt and under non-isothermal condition, equation (1) can be described as the following:(8)dφdt=βdφdT=Ae−E/RTcφ,and Kissinger, Akahira, and Sunose introduced the model-free method using the following linear equation [22, 26]:(9)lnβkTφ,k2=lnARiφEφ−EφRTφ.For each extent of conversionφ, the values of E and A can be found from the slope and the intercept when plotting lnβ/T2 against 1/T. The function iφ can be determined according to the methodology reported in the literature [27, 28].


## 3. Results and Discussions

### 3.1. The TG/DTG Results

TG/DTG curves of PETN, PBX-PN-85, and PBX-PP-85 at different heating rates (i.e., 3, 5, 7, and 10 K·min−1) were recorded and are shown in Figure 1.Figure 1
TG/DTG curves of PETN, PBX-PN-85, and PBX-PP-85 at several heating rates.
(a)(b)(c)For comparison, the characteristic parameters of TG/DTG curves of PETN, PBX-PN-85, and PBX-PP-85 are summarized in Table3.Table 3
Kinetic parameters from TG/DTG data of PBX-PN-85, PBX-PP-85, and PETN.
Materialβ (K·min−1)TG/DTG curveTot (°C)Ti (°C)Tp (°C)Mass loss (%)Residue (%)PETN3.0165.5130.2186.098.31.75.0169.6132.0191.197.22.87.0172.7133.2196.096.33.710.0178.9135.6200.896.04.0PBX-PN-853.0162.3127.1183.995.34.75.0166.8129.8188.294.45.67.0169.0130.9194.594.06.010.0175.3133.5198.893.16.9PBX-PP-853.0168.6136.2189.493.76.35.0172.5138.0194.293.16.97.0177.0139.9197.491.58.510.0183.2142.5203.691.98.1Tot—onset temperature of decomposition exotherm peak, Ti—the initial temperature of the mass loss, and Tp—the peak temperature of the DTG curve.It has been shown that only a single decomposition step is related to the thermal decomposition of a single PETN (i.e., one-step reaction model) in all samples. The peak temperatures in DTG curves for the thermal decomposition of PETN, PBX-PN-85, and PBX-PP-85 were observed in a range of 186–201, 183–199, and 189–204°C, respectively. Generally, the decomposition residue of all samples increases with the increasing heating rate, and the residue mass of PBX-PN-85 and PBX-PP-85 may be less than that of PETN because of the presence of inert ingredient (i.e., DOP and PS) in their formulation.As seen, the thermal decomposition of PETN occurred at a higher temperature than PBX-PN-85 but lower than that of PBX-PP-85. DTG curves in Figure1 also demonstrated the role of the binder on the thermal decomposition behaviors of PETN. Thermal stability and decomposition temperature of single PETN are clearly affected by the binders. Specifically, the energetic binder (i.e., NC-based binder) reduces the decomposition temperature of single PETN, while the nonenergetic binder (i.e., PS-based binder) hardly affects the thermal stability of PETN. Interestingly, the PBX-PN-85 sample also shows an early decomposition. This may be caused by the exothermic decomposition effect of NC and the release of the gaseous products in the early decomposition of PETN in the binder matrix [13]. Contrarily, the decomposition temperature of PBX-PP-85 has a slight increase because this PBX sample has less amount of PETN (compared to single PETN) in its composition [8].According to equation (2), the corresponding curve of the conversion φ versus temperature T was obtained and is presented in Figure 2. It can be observed that all of the φ−T curves have the form of the sigmoid function (i.e., the “S” shape curve), which is often obtained in the thermal analysis of high energetic material [29].Figure 2
φ versus T curves of PETN, PBX-PN-85, and PBX-PP-85 at several heating rates.
(a)(b)(c)As recommended by the ICTAC kinetics committee [21, 30], the modified Kissinger–Akahira–Sunose (KAS) technique was applied for investigating the relationship between the activation energy and the conversion fraction. The fitting lines at each conversion fraction are presented in Figure 3 by plotting lnβ/T2 against 1/T.Figure 3
Isoconversional plots at selected conversion degree of PETN, PBX-PN-85, and PBX-PP-85 according to the modified KAS method.
(a)(b)(c)The activation energyEφ and the pre-exponential factor A of the thermal decomposition of all samples were evaluated based on the slopes and intercepts of the plots shown in Figure 3 and are expressed in Table 4. The dependence of activation energy on the conversion by the modified KAS method is shown in Figure 4.Table 4
Kinetic parameters of PETN and PBXs by the modified KAS method.
φPETNPBX-PN-85PBX-PP-85Eφ (kJ⋅mol−1)lnA, (min−1)R2Eφ (kJ⋅mol−1)lnA, (min−1)R2Eφ (kJ⋅mol−1)lnA, (min−1)R20.1133.824.70.9955124.521.30.9835139.327.80.98800.2138.326.30.9915136.325.40.9905150.432.50.99500.3141.327.20.9850136.524.90.9945152.633.20.99150.4145.228.40.9884138.626.80.9964152.132.40.98350.5145.328.60.9935137.924.80.9924155.534.10.99800.6147.129.90.9828136.624.60.9920153.232.20.99090.7145.828.90.9845135.824.70.9960148.531.30.99450.8141.427.50.9899137.926.40.9949147.629.80.99520.9141.827.70.9854127.222.60.9935146.727.80.9960Mean144.928.6137.125.2152.432.6Figure 4
The dependence of activation energy on the conversion by the modified KAS method.The mean activation energy (Eφ) values (calculated according to the modified KAS method with the conversion range of 0.3–0.7 [13, 16]) of single PETN, PBX-PN-85, and PBX-PP-85 are 144.9, 137.1, and 152.4 kJ·mol−1, respectively.To verify these results obtained by the KAS method, the kinetic parameters were recalculated using NETZSCH Thermokinetic software according to the ASTM E698 [31]. The Eφ values of single PETN, PBX-PN-85, and PBX-PP-85 are 138.2, 129.7, and 150.9 kJ·mol−1, respectively. The results were close to those obtained from the KAS method with the maximum difference of approximately 5%, which indicates the accuracy and reliability of measurements.Compared to single PETN, the introduction of NC in PBX-PN-85 formulation leads to the decrease of its activation energy because the thermal stability of the nitrate group in NC is less than that in PETN [32, 33]. For PBX-PP-85 formulation, besides the reducing mechanical sensitivity of the explosive, an inert binder system (i.e., DOP/PS) acts as a stabilizer by trapping the active thermal decomposition product [29, 34], resulting in the higher activation energy and thus higher stability versus single PETN.


### 3.2. The VST Results

Results of TG/DTG were verified by VST studies, which were conducted at several isothermal temperatures (e.g., 125, 130, 135, and 140°C). From the relationship of the volume of released gas versus heating time in VST tests, the conversion fraction values (φ) for each time (t) were computed according to equation (3), and the curve of the conversion fraction (φ) versus heating time (t) is plotted for all samples in Figure 5.Figure 5
Conversion fractions (φ) versus heating time t for PBX-PN-85, PBX-PP-85, and PETN.
(a)(b)(c)It is obvious from Figure5 that the obtained φ−t curves of PBX-PN-85, PBX-PP-85, and PETN obey the logarithmic trend, and the decomposition rate of all samples increased when the temperature was increased. In the case of PBX-PN-85 and PETN, it is apparent that the conversion grew rapidly during the first 2000 minutes and the conversion rate is low in the next 8000 minutes until the samples were completely decomposed. Similarly, the conversion rate of PBX-PP-85 can also be derived from Figure 5. In the case of PBX-PP-85, the conversion also increased rapidly in the first 2000 minutes at 135 and 140°C, while the conversion took longer time when tested at 125 and 130°C. By comparing the results of all samples, it was found that PBX-PP-85 sample took a longer time for completing the conversion than PBX-PN-85 and PETN samples. For example, at an isothermal temperature 125°C, the time for completing the conversion of PBX-PN-85, PETN, and PBX-PP-85 was 9990, 10485, and 10970 min, respectively. The differences in the thermal stability of three samples were assessed by the activation energy employing the model-free (isoconversional) and model-fitting methods.


#### 3.2.1. Isothermal Isoconversional Method from VST Data

Next, we utilize VST data to find kinetics parameters of PBX samples. By drawing−lnt versus 1/T curve at each temperature according to equation (7) for each conversion φ, the kinetic thermal decomposition parameters of three samples were calculated and are expressed in Table 5. The dependence of activation energy on the conversion is presented in Figure 6.Table 5
Kinetic parameter of PETN and PBXs by the isothermal isoconversional method.
φPETNPBX-PN-85PBX-PP-85Eφ (kJ⋅mol−1)lnA, (min−1)R2Eφ (kJ⋅mol−1)lnA, (min−1)R2Eφ (kJ⋅mol−1)lnA, (min−1)R20.1130.229.40.9835123.924.10.9966135.631.10.97670.2137.431.80.9830127.226.30.9921138.432.50.97730.3137.932.00.9962128.627.30.9938145.634.90.98970.4139.232.50.9949132.328.90.9939147.335.70.98900.5142.633.70.9889132.329.20.9924147.936.20.98830.6143.634.00.9930135.930.50.9967147.536.40.98850.7136.832.00.9854137.531.10.9858145.536.10.99000.8136.631.90.9868128.728.60.9885145.336.50.98570.9131.230.20.9831121.426.60.9941142.036.20.9879Mean137.331.9129.728.0143.935.0Figure 6
The dependence of activation energy on the conversion using the VST technique.TheseEφ values are in the same range of the Eφ values calculated from TG/DTG data using the modified KAS method, thus demonstrating the reliability of the VST test in evaluating the Eφ values of studied samples by the isothermal isoconversional method. The Eφ values of PETN, PBX-PN-85, and PBX-PP-85 were 137.3, 129.7, and 143.9 kJ·mol−1, respectively. Besides that, the lnA values of PETN, PBX-PN-85, and PBX-PP-85 were 31.9, 28.0, and 35.0 min−1, respectively.


#### 3.2.2. Isothermal Model-Fitting Method from VST Data

The model-fitting method includes two fits:the first fit for determination of the model rate constant kT that best fits the VST data and the second fit for determination of several kinetic parameters such as pre-exponential factor A from the Arrhenius equation and the activation energy Eφ [13, 22]. Fourteen reaction models from Table 1 were selected for fitting the conversion of the thermal decomposition. Table 6 shows the results for the pre-exponential factor, the activation energy, and squares of the correlation coefficients obtained by plotting the integrated form iφ versus heating time t at each isothermal temperature.Table 6
Kinetic parameters of PETN and PBXs calculated from the model-fitting method.
ModelPETNPBX-PN-85PBX-PP-85Eφ (kJ mol−1)lnA, (min−1)R2Eφ (kJ mol−1)lnA, (min−1)R2Eφ (kJ mol−1)lnA, (min−1)R21137.432.40.9818122.027.60.9828144.434.60.98962138.532.50.9815121.627.20.9830144.534.40.98963138.332.30.9827121.927.10.9821144.534.00.98924139.133.60.9900123.528.70.9846146.435.80.99085138.533.00.9890122.828.20.9833146.235.40.99076134.131.40.9845122.927.90.9838146.135.10.99047139.433.00.9838123.928.30.9853146.935.30.99068139.132.70.9935123.928.10.9844146.935.10.99089139.432.30.9862123.728.50.9862147.635.80.990510139.833.30.9754119.727.20.9855146.735.20.988611140.632.60.9860125.628.00.9918148.534.90.991312139.734.40.9853124.929.80.9853148.236.60.991413141.336.20.9846124.831.10.9846150.636.90.994014142.638.20.9834129.634.20.9834152.440.00.9903From Table6, it was observed that the kinetic parameters (specifically, pre-exponential factor and activation energy), which are derived from the VST data utilizing the model-fitting method, are essentially independent of the reaction model used. It means that there is an insignificant difference between the kinetic parameters (i.e., Eφ and A) calculated using different model reactions. For selecting the best model fitting, the correlation coefficient R2 was used as a parameter. It was found that the best model that describes the thermal decomposition of PETN is contracting volume (R3) and the best model for PBX-PN-85 is 3-D diffusion-Jander (D3) while the best model for PBX-PP-85 is second-order (F2).On the other hand, we found that the activation energy values calculated by the isothermal isoconversional method are approximate to those derived from the isothermal model-fitting method using the VST technique. To make a comparison, kinetic parameters obtained by different methods are presented in Table7.Table 7
Kinetic parameters of PETN and PBXs obtained by several methods.
SampleTG/DTG using KAS methodVacuum stability test (VST)Isoconversional methodModel-fitting methodEφ (kJ·mol−1)lnA (min−1)Eφ (kJ·mol−1)lnA (min−1)Eφ (kJ·mol−1)lnA (min−1)PETN144.928.6137.331.9139.132.7PBX-PN-85137.125.2129.728.0125.628.0PBX-PP-85152.432.6143.935.0150.636.9The values in Table7 demonstrate the role of the binder on the initiation activation energy of PETN. The energetic binder based on NC reduces the activation energy (i.e., PBX-PN-85 has lower thermal stability than that of PETN), while the introduction of the nonenergetic material based on PS increases the activation energy of the explosive.The ranges ofEφ values obtained by different methods were 125.6–137.1 kJ·mol−1 for PBX-PN-85, 137.3–144.9 kJ·mol−1 for PETN, and 143.9–152.4 kJ·mol−1 for PBX-PP-85. These ranges are narrow, proving the equivalence of these methods with each other. It is also noteworthy that Eφ of PETN obtained in this work is in good agreement with that reported by Pouretedal et al. [15] (135.1–136.9 kJ·mol−1) and by Künzel et al. [16] (137.4 ± 3.4 kJ·mol−1).According to Brill et al. [35], the Eφ and lnA values—determined from isothermal techniques in the same conditions—usually compensate each other. They found that all of Eφ−lnA pairs that lie in the regression line represent legitimate measurement and reliable results. We have plotted all Eφ values of PETN, PBX-PN-85, and PBX-PP-85 versus corresponding lnA values (presented in Table 7) in Figure 7.Figure 7
The linear dependence ofEφ with lnA of all samples obtained from the VST test.
(a)(b)It can be observed the linear relationship ofEφ with lnA values, which confirm the accuracy of the kinetic parameters obtained by the isoconversional method and model-fitting method using the VST test.


## 3.1. The TG/DTG Results

TG/DTG curves of PETN, PBX-PN-85, and PBX-PP-85 at different heating rates (i.e., 3, 5, 7, and 10 K·min−1) were recorded and are shown in Figure 1.Figure 1
TG/DTG curves of PETN, PBX-PN-85, and PBX-PP-85 at several heating rates.
(a)(b)(c)For comparison, the characteristic parameters of TG/DTG curves of PETN, PBX-PN-85, and PBX-PP-85 are summarized in Table3.Table 3
Kinetic parameters from TG/DTG data of PBX-PN-85, PBX-PP-85, and PETN.
Materialβ (K·min−1)TG/DTG curveTot (°C)Ti (°C)Tp (°C)Mass loss (%)Residue (%)PETN3.0165.5130.2186.098.31.75.0169.6132.0191.197.22.87.0172.7133.2196.096.33.710.0178.9135.6200.896.04.0PBX-PN-853.0162.3127.1183.995.34.75.0166.8129.8188.294.45.67.0169.0130.9194.594.06.010.0175.3133.5198.893.16.9PBX-PP-853.0168.6136.2189.493.76.35.0172.5138.0194.293.16.97.0177.0139.9197.491.58.510.0183.2142.5203.691.98.1Tot—onset temperature of decomposition exotherm peak, Ti—the initial temperature of the mass loss, and Tp—the peak temperature of the DTG curve.It has been shown that only a single decomposition step is related to the thermal decomposition of a single PETN (i.e., one-step reaction model) in all samples. The peak temperatures in DTG curves for the thermal decomposition of PETN, PBX-PN-85, and PBX-PP-85 were observed in a range of 186–201, 183–199, and 189–204°C, respectively. Generally, the decomposition residue of all samples increases with the increasing heating rate, and the residue mass of PBX-PN-85 and PBX-PP-85 may be less than that of PETN because of the presence of inert ingredient (i.e., DOP and PS) in their formulation.As seen, the thermal decomposition of PETN occurred at a higher temperature than PBX-PN-85 but lower than that of PBX-PP-85. DTG curves in Figure1 also demonstrated the role of the binder on the thermal decomposition behaviors of PETN. Thermal stability and decomposition temperature of single PETN are clearly affected by the binders. Specifically, the energetic binder (i.e., NC-based binder) reduces the decomposition temperature of single PETN, while the nonenergetic binder (i.e., PS-based binder) hardly affects the thermal stability of PETN. Interestingly, the PBX-PN-85 sample also shows an early decomposition. This may be caused by the exothermic decomposition effect of NC and the release of the gaseous products in the early decomposition of PETN in the binder matrix [13]. Contrarily, the decomposition temperature of PBX-PP-85 has a slight increase because this PBX sample has less amount of PETN (compared to single PETN) in its composition [8].According to equation (2), the corresponding curve of the conversion φ versus temperature T was obtained and is presented in Figure 2. It can be observed that all of the φ−T curves have the form of the sigmoid function (i.e., the “S” shape curve), which is often obtained in the thermal analysis of high energetic material [29].Figure 2
φ versus T curves of PETN, PBX-PN-85, and PBX-PP-85 at several heating rates.
(a)(b)(c)As recommended by the ICTAC kinetics committee [21, 30], the modified Kissinger–Akahira–Sunose (KAS) technique was applied for investigating the relationship between the activation energy and the conversion fraction. The fitting lines at each conversion fraction are presented in Figure 3 by plotting lnβ/T2 against 1/T.Figure 3
Isoconversional plots at selected conversion degree of PETN, PBX-PN-85, and PBX-PP-85 according to the modified KAS method.
(a)(b)(c)The activation energyEφ and the pre-exponential factor A of the thermal decomposition of all samples were evaluated based on the slopes and intercepts of the plots shown in Figure 3 and are expressed in Table 4. The dependence of activation energy on the conversion by the modified KAS method is shown in Figure 4.Table 4
Kinetic parameters of PETN and PBXs by the modified KAS method.
φPETNPBX-PN-85PBX-PP-85Eφ (kJ⋅mol−1)lnA, (min−1)R2Eφ (kJ⋅mol−1)lnA, (min−1)R2Eφ (kJ⋅mol−1)lnA, (min−1)R20.1133.824.70.9955124.521.30.9835139.327.80.98800.2138.326.30.9915136.325.40.9905150.432.50.99500.3141.327.20.9850136.524.90.9945152.633.20.99150.4145.228.40.9884138.626.80.9964152.132.40.98350.5145.328.60.9935137.924.80.9924155.534.10.99800.6147.129.90.9828136.624.60.9920153.232.20.99090.7145.828.90.9845135.824.70.9960148.531.30.99450.8141.427.50.9899137.926.40.9949147.629.80.99520.9141.827.70.9854127.222.60.9935146.727.80.9960Mean144.928.6137.125.2152.432.6Figure 4
The dependence of activation energy on the conversion by the modified KAS method.The mean activation energy (Eφ) values (calculated according to the modified KAS method with the conversion range of 0.3–0.7 [13, 16]) of single PETN, PBX-PN-85, and PBX-PP-85 are 144.9, 137.1, and 152.4 kJ·mol−1, respectively.To verify these results obtained by the KAS method, the kinetic parameters were recalculated using NETZSCH Thermokinetic software according to the ASTM E698 [31]. The Eφ values of single PETN, PBX-PN-85, and PBX-PP-85 are 138.2, 129.7, and 150.9 kJ·mol−1, respectively. The results were close to those obtained from the KAS method with the maximum difference of approximately 5%, which indicates the accuracy and reliability of measurements.Compared to single PETN, the introduction of NC in PBX-PN-85 formulation leads to the decrease of its activation energy because the thermal stability of the nitrate group in NC is less than that in PETN [32, 33]. For PBX-PP-85 formulation, besides the reducing mechanical sensitivity of the explosive, an inert binder system (i.e., DOP/PS) acts as a stabilizer by trapping the active thermal decomposition product [29, 34], resulting in the higher activation energy and thus higher stability versus single PETN.


## 3.2. The VST Results

Results of TG/DTG were verified by VST studies, which were conducted at several isothermal temperatures (e.g., 125, 130, 135, and 140°C). From the relationship of the volume of released gas versus heating time in VST tests, the conversion fraction values (φ) for each time (t) were computed according to equation (3), and the curve of the conversion fraction (φ) versus heating time (t) is plotted for all samples in Figure 5.Figure 5
Conversion fractions (φ) versus heating time t for PBX-PN-85, PBX-PP-85, and PETN.
(a)(b)(c)It is obvious from Figure5 that the obtained φ−t curves of PBX-PN-85, PBX-PP-85, and PETN obey the logarithmic trend, and the decomposition rate of all samples increased when the temperature was increased. In the case of PBX-PN-85 and PETN, it is apparent that the conversion grew rapidly during the first 2000 minutes and the conversion rate is low in the next 8000 minutes until the samples were completely decomposed. Similarly, the conversion rate of PBX-PP-85 can also be derived from Figure 5. In the case of PBX-PP-85, the conversion also increased rapidly in the first 2000 minutes at 135 and 140°C, while the conversion took longer time when tested at 125 and 130°C. By comparing the results of all samples, it was found that PBX-PP-85 sample took a longer time for completing the conversion than PBX-PN-85 and PETN samples. For example, at an isothermal temperature 125°C, the time for completing the conversion of PBX-PN-85, PETN, and PBX-PP-85 was 9990, 10485, and 10970 min, respectively. The differences in the thermal stability of three samples were assessed by the activation energy employing the model-free (isoconversional) and model-fitting methods.


### 3.2.1. Isothermal Isoconversional Method from VST Data

Next, we utilize VST data to find kinetics parameters of PBX samples. By drawing−lnt versus 1/T curve at each temperature according to equation (7) for each conversion φ, the kinetic thermal decomposition parameters of three samples were calculated and are expressed in Table 5. The dependence of activation energy on the conversion is presented in Figure 6.Table 5
Kinetic parameter of PETN and PBXs by the isothermal isoconversional method.
φPETNPBX-PN-85PBX-PP-85Eφ (kJ⋅mol−1)lnA, (min−1)R2Eφ (kJ⋅mol−1)lnA, (min−1)R2Eφ (kJ⋅mol−1)lnA, (min−1)R20.1130.229.40.9835123.924.10.9966135.631.10.97670.2137.431.80.9830127.226.30.9921138.432.50.97730.3137.932.00.9962128.627.30.9938145.634.90.98970.4139.232.50.9949132.328.90.9939147.335.70.98900.5142.633.70.9889132.329.20.9924147.936.20.98830.6143.634.00.9930135.930.50.9967147.536.40.98850.7136.832.00.9854137.531.10.9858145.536.10.99000.8136.631.90.9868128.728.60.9885145.336.50.98570.9131.230.20.9831121.426.60.9941142.036.20.9879Mean137.331.9129.728.0143.935.0Figure 6
The dependence of activation energy on the conversion using the VST technique.TheseEφ values are in the same range of the Eφ values calculated from TG/DTG data using the modified KAS method, thus demonstrating the reliability of the VST test in evaluating the Eφ values of studied samples by the isothermal isoconversional method. The Eφ values of PETN, PBX-PN-85, and PBX-PP-85 were 137.3, 129.7, and 143.9 kJ·mol−1, respectively. Besides that, the lnA values of PETN, PBX-PN-85, and PBX-PP-85 were 31.9, 28.0, and 35.0 min−1, respectively.


### 3.2.2. Isothermal Model-Fitting Method from VST Data

The model-fitting method includes two fits:the first fit for determination of the model rate constant kT that best fits the VST data and the second fit for determination of several kinetic parameters such as pre-exponential factor A from the Arrhenius equation and the activation energy Eφ [13, 22]. Fourteen reaction models from Table 1 were selected for fitting the conversion of the thermal decomposition. Table 6 shows the results for the pre-exponential factor, the activation energy, and squares of the correlation coefficients obtained by plotting the integrated form iφ versus heating time t at each isothermal temperature.Table 6
Kinetic parameters of PETN and PBXs calculated from the model-fitting method.
ModelPETNPBX-PN-85PBX-PP-85Eφ (kJ mol−1)lnA, (min−1)R2Eφ (kJ mol−1)lnA, (min−1)R2Eφ (kJ mol−1)lnA, (min−1)R21137.432.40.9818122.027.60.9828144.434.60.98962138.532.50.9815121.627.20.9830144.534.40.98963138.332.30.9827121.927.10.9821144.534.00.98924139.133.60.9900123.528.70.9846146.435.80.99085138.533.00.9890122.828.20.9833146.235.40.99076134.131.40.9845122.927.90.9838146.135.10.99047139.433.00.9838123.928.30.9853146.935.30.99068139.132.70.9935123.928.10.9844146.935.10.99089139.432.30.9862123.728.50.9862147.635.80.990510139.833.30.9754119.727.20.9855146.735.20.988611140.632.60.9860125.628.00.9918148.534.90.991312139.734.40.9853124.929.80.9853148.236.60.991413141.336.20.9846124.831.10.9846150.636.90.994014142.638.20.9834129.634.20.9834152.440.00.9903From Table6, it was observed that the kinetic parameters (specifically, pre-exponential factor and activation energy), which are derived from the VST data utilizing the model-fitting method, are essentially independent of the reaction model used. It means that there is an insignificant difference between the kinetic parameters (i.e., Eφ and A) calculated using different model reactions. For selecting the best model fitting, the correlation coefficient R2 was used as a parameter. It was found that the best model that describes the thermal decomposition of PETN is contracting volume (R3) and the best model for PBX-PN-85 is 3-D diffusion-Jander (D3) while the best model for PBX-PP-85 is second-order (F2).On the other hand, we found that the activation energy values calculated by the isothermal isoconversional method are approximate to those derived from the isothermal model-fitting method using the VST technique. To make a comparison, kinetic parameters obtained by different methods are presented in Table7.Table 7
Kinetic parameters of PETN and PBXs obtained by several methods.
SampleTG/DTG using KAS methodVacuum stability test (VST)Isoconversional methodModel-fitting methodEφ (kJ·mol−1)lnA (min−1)Eφ (kJ·mol−1)lnA (min−1)Eφ (kJ·mol−1)lnA (min−1)PETN144.928.6137.331.9139.132.7PBX-PN-85137.125.2129.728.0125.628.0PBX-PP-85152.432.6143.935.0150.636.9The values in Table7 demonstrate the role of the binder on the initiation activation energy of PETN. The energetic binder based on NC reduces the activation energy (i.e., PBX-PN-85 has lower thermal stability than that of PETN), while the introduction of the nonenergetic material based on PS increases the activation energy of the explosive.The ranges ofEφ values obtained by different methods were 125.6–137.1 kJ·mol−1 for PBX-PN-85, 137.3–144.9 kJ·mol−1 for PETN, and 143.9–152.4 kJ·mol−1 for PBX-PP-85. These ranges are narrow, proving the equivalence of these methods with each other. It is also noteworthy that Eφ of PETN obtained in this work is in good agreement with that reported by Pouretedal et al. [15] (135.1–136.9 kJ·mol−1) and by Künzel et al. [16] (137.4 ± 3.4 kJ·mol−1).According to Brill et al. [35], the Eφ and lnA values—determined from isothermal techniques in the same conditions—usually compensate each other. They found that all of Eφ−lnA pairs that lie in the regression line represent legitimate measurement and reliable results. We have plotted all Eφ values of PETN, PBX-PN-85, and PBX-PP-85 versus corresponding lnA values (presented in Table 7) in Figure 7.Figure 7
The linear dependence ofEφ with lnA of all samples obtained from the VST test.
(a)(b)It can be observed the linear relationship ofEφ with lnA values, which confirm the accuracy of the kinetic parameters obtained by the isoconversional method and model-fitting method using the VST test.


## 3.2.1. Isothermal Isoconversional Method from VST Data

Next, we utilize VST data to find kinetics parameters of PBX samples. By drawing−lnt versus 1/T curve at each temperature according to equation (7) for each conversion φ, the kinetic thermal decomposition parameters of three samples were calculated and are expressed in Table 5. The dependence of activation energy on the conversion is presented in Figure 6.Table 5
Kinetic parameter of PETN and PBXs by the isothermal isoconversional method.
φPETNPBX-PN-85PBX-PP-85Eφ (kJ⋅mol−1)lnA, (min−1)R2Eφ (kJ⋅mol−1)lnA, (min−1)R2Eφ (kJ⋅mol−1)lnA, (min−1)R20.1130.229.40.9835123.924.10.9966135.631.10.97670.2137.431.80.9830127.226.30.9921138.432.50.97730.3137.932.00.9962128.627.30.9938145.634.90.98970.4139.232.50.9949132.328.90.9939147.335.70.98900.5142.633.70.9889132.329.20.9924147.936.20.98830.6143.634.00.9930135.930.50.9967147.536.40.98850.7136.832.00.9854137.531.10.9858145.536.10.99000.8136.631.90.9868128.728.60.9885145.336.50.98570.9131.230.20.9831121.426.60.9941142.036.20.9879Mean137.331.9129.728.0143.935.0Figure 6
The dependence of activation energy on the conversion using the VST technique.TheseEφ values are in the same range of the Eφ values calculated from TG/DTG data using the modified KAS method, thus demonstrating the reliability of the VST test in evaluating the Eφ values of studied samples by the isothermal isoconversional method. The Eφ values of PETN, PBX-PN-85, and PBX-PP-85 were 137.3, 129.7, and 143.9 kJ·mol−1, respectively. Besides that, the lnA values of PETN, PBX-PN-85, and PBX-PP-85 were 31.9, 28.0, and 35.0 min−1, respectively.


## 3.2.2. Isothermal Model-Fitting Method from VST Data

The model-fitting method includes two fits:the first fit for determination of the model rate constant kT that best fits the VST data and the second fit for determination of several kinetic parameters such as pre-exponential factor A from the Arrhenius equation and the activation energy Eφ [13, 22]. Fourteen reaction models from Table 1 were selected for fitting the conversion of the thermal decomposition. Table 6 shows the results for the pre-exponential factor, the activation energy, and squares of the correlation coefficients obtained by plotting the integrated form iφ versus heating time t at each isothermal temperature.Table 6
Kinetic parameters of PETN and PBXs calculated from the model-fitting method.
ModelPETNPBX-PN-85PBX-PP-85Eφ (kJ mol−1)lnA, (min−1)R2Eφ (kJ mol−1)lnA, (min−1)R2Eφ (kJ mol−1)lnA, (min−1)R21137.432.40.9818122.027.60.9828144.434.60.98962138.532.50.9815121.627.20.9830144.534.40.98963138.332.30.9827121.927.10.9821144.534.00.98924139.133.60.9900123.528.70.9846146.435.80.99085138.533.00.9890122.828.20.9833146.235.40.99076134.131.40.9845122.927.90.9838146.135.10.99047139.433.00.9838123.928.30.9853146.935.30.99068139.132.70.9935123.928.10.9844146.935.10.99089139.432.30.9862123.728.50.9862147.635.80.990510139.833.30.9754119.727.20.9855146.735.20.988611140.632.60.9860125.628.00.9918148.534.90.991312139.734.40.9853124.929.80.9853148.236.60.991413141.336.20.9846124.831.10.9846150.636.90.994014142.638.20.9834129.634.20.9834152.440.00.9903From Table6, it was observed that the kinetic parameters (specifically, pre-exponential factor and activation energy), which are derived from the VST data utilizing the model-fitting method, are essentially independent of the reaction model used. It means that there is an insignificant difference between the kinetic parameters (i.e., Eφ and A) calculated using different model reactions. For selecting the best model fitting, the correlation coefficient R2 was used as a parameter. It was found that the best model that describes the thermal decomposition of PETN is contracting volume (R3) and the best model for PBX-PN-85 is 3-D diffusion-Jander (D3) while the best model for PBX-PP-85 is second-order (F2).On the other hand, we found that the activation energy values calculated by the isothermal isoconversional method are approximate to those derived from the isothermal model-fitting method using the VST technique. To make a comparison, kinetic parameters obtained by different methods are presented in Table7.Table 7
Kinetic parameters of PETN and PBXs obtained by several methods.
SampleTG/DTG using KAS methodVacuum stability test (VST)Isoconversional methodModel-fitting methodEφ (kJ·mol−1)lnA (min−1)Eφ (kJ·mol−1)lnA (min−1)Eφ (kJ·mol−1)lnA (min−1)PETN144.928.6137.331.9139.132.7PBX-PN-85137.125.2129.728.0125.628.0PBX-PP-85152.432.6143.935.0150.636.9The values in Table7 demonstrate the role of the binder on the initiation activation energy of PETN. The energetic binder based on NC reduces the activation energy (i.e., PBX-PN-85 has lower thermal stability than that of PETN), while the introduction of the nonenergetic material based on PS increases the activation energy of the explosive.The ranges ofEφ values obtained by different methods were 125.6–137.1 kJ·mol−1 for PBX-PN-85, 137.3–144.9 kJ·mol−1 for PETN, and 143.9–152.4 kJ·mol−1 for PBX-PP-85. These ranges are narrow, proving the equivalence of these methods with each other. It is also noteworthy that Eφ of PETN obtained in this work is in good agreement with that reported by Pouretedal et al. [15] (135.1–136.9 kJ·mol−1) and by Künzel et al. [16] (137.4 ± 3.4 kJ·mol−1).According to Brill et al. [35], the Eφ and lnA values—determined from isothermal techniques in the same conditions—usually compensate each other. They found that all of Eφ−lnA pairs that lie in the regression line represent legitimate measurement and reliable results. We have plotted all Eφ values of PETN, PBX-PN-85, and PBX-PP-85 versus corresponding lnA values (presented in Table 7) in Figure 7.Figure 7
The linear dependence ofEφ with lnA of all samples obtained from the VST test.
(a)(b)It can be observed the linear relationship ofEφ with lnA values, which confirm the accuracy of the kinetic parameters obtained by the isoconversional method and model-fitting method using the VST test.


## 4. Conclusions

The thermal decomposition behavior and kinetics of PETN, PBX-PN-85 (85% PETN, 15% NC-based binder), and PBX-PP-85 (85% PETN, 15% PS-based binder) were examined by TG/DTG and VST techniques using model-fitting and isoconversional methods. The activation energy values computed by the non-isothermal technique (modified KAS method) were close to those calculated by model-fitting and isothermal isoconversional methods, ranging from 125.6–137.1 kJ·mol−1 for PBX-PN-85, 137.3–144.9 kJ·mol−1 for PETN, and 143.9–152.4 kJ·mol−1 for PBX-PP-85. The results indicate that the NC-based binder accelerates the thermal decomposition of PBX, thus reducing the activation energy of single PETN. One more valuable conclusion of this study is the influence of the PS-based binder. Specifically, the presence of the PS-based binder causes no adverse effect on the thermal stability of the main explosive in PBX composition.



---

*Source: 9260818-2020-07-11.xml*