Physicochemical analysis of biomass (dried basis).
The efficient and high yielding acid-base and Organosolv methods were studied for cellulose isolation from date palm lignocellulose waste biomass and thereafter analyzed for nonisothermal kinetic and thermodynamic parameter determination using model-free methods. The structural and chemical characterization of the isolated celluloses revealed structures and functional groups characteristics of cellulose. Thermal decomposition analysis revealed one major peak with average mass loss of 72.51 ± 0.7% and 55.82 ± 1.1% for the acid-base and Organosolv method, respectively. This occurred in the temperature region between 250 and 350°C associated with cellulose degradation and contrasted with the three peaks detected in the original biomass. The kinetic and thermodynamic results revealed a strong relationship between the average activation energy and average change in enthalpy with a difference of 5.23 and 147.07 kJmol−1 for Organosolv and acid-base methods, respectively. The Gibbs’s free energy results revealed that Organosolv cellulose pyrolysis would reach equilibrium faster in KAS, Starink and FWO models with average ΔG values of 115.80 ± 36.62, 115.89 ± 36.65, and 119.45 ± 37.98 kJmol−1, respectively. The acid-base method for FWO model gave negative entropy values. The Malek method revealed the acid-base and Organoslv cellulose pyrolysis mechanism as (gα=−ln1−α14) and (gα=−ln1−α13), characterized by random nucleation and growth, respectively.
- lignocellulose biomass
- nonisothermal kinetics
- isolation methods
Lignocellulosic biomass is the most abundant, renewable, and one of the cheapest carbon neutral raw materials in the biosphere that can be used to produce sustainable products such as biofuels, using different technologies . The lignocellulosic biomass consists of mainly cellulose carbohydrate polymer, hemicellulose, and the aromatic component, lignin . Lignocellulose biomass can store up to 47 MJkg−1 more energy than lithium ion batteries (0.8 MJkg−1) . Lignocellulose biomass is considered a potential candidate to sustainable green alternative source of energy and chemicals due to its high energy density, volatile matter content, and global widespread . The release of volatile matter and other contents in biomass has been extensively studied using pyrolysis technology. Pyrolysis involves the conversion of biomass into bio-oil, gases (volatile matter) and biochar, in the absence of oxygen [4, 5]. The technique is robust and essential in providing vital knowledge of kinetics of devolatilization of any biomass prior to further processing via different conversion technologies. In addition, pyrolysis is effective in reducing the bulky biomass into uniform, energy dense, and easily transportable fuel . Despite always being the first stage in most combustion or gasification process, there are no accurate and enough data on the kinetics and reaction mechanisms of different lignocellulosic biomass . This is attributed to the complexity and the varying physico-chemical properties in different lignocellulosic biomass . In addition, there may be many reactions occurring from the extremely complex pyrolysis process of the lignocellulose biomass . Therefore, developing accurate kinetic models to account for all reactions taking place remains a challenge . Isolation of the complex lignocellulosic biomass into individual fractions and characterization of the individual fractions can provide a better understanding of the combustion kinetics and reaction mechanism. Different biomass isolation/extraction techniques such as liquid-liquid, liquid-solid, acid-base, ultrasound, and microwave-assisted extractions, among others, have been reported before . The choice of the method depends on the biomass type and its fraction to be isolated . Hence, each procedure affects the sample’s product yields and physical, chemical, kinetic, and thermodynamic properties differently. Despite the studies on the yields and operating parameters such as solvent and time , less or no information is available regarding the kinetic and thermodynamic parametric studies for the combustion of the isolated lignocellulose fractions to assess the difference in the extraction processes. Different general kinetic models on lignocellulose biomass have been suggested [10, 11].
Date palm waste constitutes about 500,000 metric tons per year from ca. 44 million date palm trees found in the United Arab Emirates where this research was conducted. The aims of this research are to isolate cellulose from date palm lignocellulose complex using low concentration acid-base solutions and Organosolv techniques and to model nonisothermal combustion kinetics using model-free methods and finally to predict the most probable mechanistic reaction mechanism of the isolated celluloses. Using thermal-gravimetric technique at different heating rates, kinetic and thermodynamic parameters were calculated using model-free methods, namely Kissinger-Akahira-Sunose (KAS), Flynn-Wall-Ozawa (FWO), and Starink models. The FWO model-free method compensates the experimental measurement errors. However, the KAS and Starink methods depend on the choice of good constant degree of conversion from the derivative mass loss function to provide precision of the kinetic data . In addition, application of different model-free methods involves wide conversion range that allows for study of change in mechanism during a reaction and reduces mass transfer limitations by using multiple heating rates .
2. Isolation techniques and nonisothermal kinetic studies
The rachis part of adult date palm waste (DPW) (10–15 years old) was supplied by the UAE University farm, Al foah, Al Ain. The samples were ground to 180-micron particle size to reduce the effects of heat and mass transfer limitations. All the solvents (ACS grade) and reagents were supplied by Sigma Aldrich and were used with no further purification. Prior to cellulose isolation, 10 g biomass was valorized with benzene/ethanol (2, 1 v/v) for 48 h using Soxhlet extraction to reduce extractives such as waxes and resins surrounding the lignocellulose complex.
2.1 Acid-base and Organosolv cellulose isolation methods
For acid-base isolation; DPW (5 g) extractive-free sample of particle size 180 μm was placed in a 250-mL beaker and leached with 200 mL of 0.1 M HCl while heating at 100°C for 2 h under stirring at 150 rpm. After vacuum filtration, the cellulose and lignin-rich residue was washed with 20 mL of deionized water to remove any residual hemicellulose and then air dried overnight. The hemicellulose was solubilized by HCl and heating due to its labile nature making it easy to dissolve out of the lignocellulose complex. The cellulose- and lignin-rich residue was further treated with 200 mL of 0.1 M NaOH while heating at 100°C for 2 h under constant stirring at 150 rpm. After subsequent vacuum filtration of the mixture, the cellulose-rich residue was washed with 20 mL of 0.1 M NaOH to remove any residual lignin. The isolated cellulose was air dried under laboratory conditions overnight prior to characterization.
For Organosolv isolation, DPW were isolated using methanol/water solvents as reported in literature, with some modifications . The 6.7 g sample of particle size 180 μm was placed in high pressure/temperature reactor vessel (Parr 4848, U.S.A). A mixture of 84 mL sulfuric acid (0.045 N), 13.4 mL formaldehyde (37 wt.%), and 84 mL methanol was added to the reactor vessel containing the sample. The reactor was sealed and purged with nitrogen gas (6–10 bars), and the reaction was performed for 1 h at 160°C under constant stirring at 700 rpm. The product mixture was vacuum filtered after cooling to room temperature. The cellulose-rich residues were air dried overnight prior to characterization, and the yield was determined by a gravimetric analysis technique. The ultimate analysis was conducted, and the results are recorded in Table 1.
|Proximate and ultimate analyses||Date palm lignocellulose||Organosolv cellulose||Acid-base cellulose|
|Moisture content (wt.%)||6.72 ± 0.4||7.08 ± 0.4||8.72 ± 0.4|
|Volatile matter (wt.%)||78.62 ± 0.04||65.22 ± 0.02||66.92 ± 0.01|
|Ash content (wt.%)||6.12 ± 0.1||7.36 ± 0.04||7.24 ± 0.01|
|Fixed carbon (wt.%)||5.40 ± 0.01||4.40 ± 0.10||4.80 ± 0.14|
|Cellulose yield (wt.%)||—||43.15 ± 2.40||64.15 ± 2.40|
2.2 Fourier transform infrared (FTIR) spectroscopy
The IRTrace-100 FTIR spectrophotometer (Shimadzu, Kyoto, Japan) was used for the FTIR analysis. The extracted celluloses were analyzed to investigate the difference in the functional groups after extraction. The spectral results were recorded within a range of 500–4000 cm−1 wavelength using 4 cm−1 spectral resolution and 34 scans. Figure 1 shows several major absorption bands and the difference between the samples. DPW sample before isolation showed typical lignocellulose strong band absorption bonds. For instance, the bands at 1037 cm−1 indicating C▬O, C〓C and C▬C bond stretching, between 2840 and 2926 cm−1 indicating C▬H stretching and 3200–3474 cm−1 for O▬H stretching, were assigned to cellulose, hemicellulose and lignin, respectively. Similar results were reported for other biomass in the literature . However, there were noticeable changes in the functional groups of celluloses from the same biomass with reduction in the peak intensity, an indication of component(s) removal (plausibly lignin and hemicellulose). For instance, the would-be lignin and hemicellulose band intensity at 1037 cm−1 greatly reduced an indication of component(s) removal. The C▬H bond stretching in the region of 2840 assigned to lignin in DPW was absent in isolated cellulose samples. In addition, the reduced peak intensity between 845 and 1156 cm−1 associated with the C▬O▬C asymmetrical stretching and glycosidic bond, a characteristic of cellulose, was observed. Furthermore, the decrease in OH vibration strength around 3200–3474 cm−1 indicates a reduction in some of the OH-containing compounds which are phenolics from lignin. It is worth to note that both extraction methods showed similar functional groups except that the Organosolv cellulose had C▬H bond assigned to lignin in the region around 2326 and 2363 cm−1 which was absent for the acid-base cellulose samples. The FTIR results showed the effectiveness of the cellulose isolation methods from DPW biomass complex. The samples were further characterized for their morphological differences using SEM imaging technique.
2.3 Scanning electron microscopy (SEM)
The structural morphologies of the isolated cellulose were analyzed using the scanning electron microscope (JEOL Neoscope JCM-5000, Tokyo Japan). The samples were Au/C coated using vacuum spatter while clamped on the sample holder. The images were captured on spot size of 40 using 10 kV. The SEM results in Figure 2 show a difference in the structural morphologies between the cellulose samples from the two methods. Figure 2(a) shows the original DPW with ring-like structures (see the arrow point) plausible to be the cellulose chiral nematic ordering, surrounded by irregular shaped structures which could be assumed to be lignin and hemicellulose. Acid-base cellulose in Figure 2(b) showed porous surface similar to those reported for cellulose from the teak wood . However, Organosolv cellulose (Figure 2(c)) showed an aggregate of cellulose block structure with uneven polished surface.
2.4 Thermogravimetric analysis (TGA)
The combustion characteristics of isolated celluloses were studied using thermogravimetric analysis (TGA). The analysis was done on a TGA (Q500, TA instrument). Samples of 6 mg (±1.0) were first equilibrated at 25°C for 5 min and then heated at specific heating rates of 10, 15, 20, and 25°C/min to 900°C. The process was performed under constant nitrogen environment flowing at 20 mL/min. As the thermal decomposition progressed, the change in weight was recorded continuously as a function of temperature and time. Figure 3 shows the isoconversion versus temperature at different heating rates for the isolated celluloses from DPW. The conversion curves for acid-base (colored) and Organosolv (black) methods below 300 and 340°C, respectively, showed similar thermal decomposition patterns at all heating rates. There was a slight shift toward higher temperature side with increasing heating rates, possibly due to the increasing thermal energy in the system . However, at higher temperatures, the conversion pattern changed for both methods, possibly due to the change in the degradation chemistry of components under pyrolysis. It is worth to note that Organosolv cellulose showed better thermal stability than the acid-base cellulose. Figure 4 shows the differential thermogravimetric (DTG) results against temperature at different heating rates for the DPW and the isolated celluloses. The results showed a typical thermal degradation of lignocellulose biomass. The curves of both samples moved downward as the heating rate increased, owing to the shorter reaction time at increasing temperature, a phenomenon known as thermal hysteresis. However, Figure 4(a) showed a noticeable difference between the peak mass loss patterns compared to Figure 4(b) for the two cellulose methods, which suggests a difference in the degradation chemistry. Table 2 shows the temperature ranges that define the major stages of mass loss in response to increasing temperature for isolated celluloses from both methods, as given in Table 3. Stage I started from minimum temperature, Tmin to T1, the total average celluloses mass loss for acid-base and Organosolv were 9.46 ± 0.1 and 5.28 ± 0.1%, respectively. This was attributed to the inherent moisture and water molecules embedded in the intercellular and intracellular void spaces of the celluloses. Stage II, from T1 to T3 for both methods, there was only one major clear peak (Figure 4(b)) and the average mass loss in this region was 72.51 ± 0.7 and 55.82 ± 1.1%, for acid-base and Organosolv celluloses, respectively. The weight loss in this stage is associated with pyrolysis of mainly cellulose and to a lesser extent hemicellulose . Compared to the three peaks observed in Figure 4(a) for the original DPW, this clearly shows that both methods were effective for cellulose isolation from the complex lignocellulose matrix of DPW. Stage III had total average mass loss of 8.68 ± 1.2 and 33.08 ± 0.8%, for acid-base and Organosolv methods, respectively. This represented combustion of the carbonaceous and some part of char oxidation . In addition, the higher mass loss for Organosolv cellulose was plausibly due to residual lignin. Moreover, the FTIR results showed some lignin functional groups for this method. The last stage was associated with charring process and ash formation. The average total mass loss for acid-base and Organosolv methods in this stage were 9.57 ± 1.3 and 5.20 ± 0.4%. TGA analysis data was used for kinetic modeling using the model-free methods.
|Heating rate (°C min−1)||Tmin (°C)||T1 (°C)||T2 (°C)||T3 (°C)||Tmax (°C)|
|Stages||Temperature||Heating rate (°C min−1)|
|Stage I, WL%||Tmin–T1||10.06||9.17||9.33||9.26|
|Stage II, WL %||T1–T3||71.59||72.76||73.16||72.52|
|Stage III, WL %||T3–Tmax||10.18||8.88||7.30||8.35|
|Final residue at 900-100 °C (%)||7.99||9.19||11.21||9.87|
|Stage I, WL %||Tmin–T1||4.35||5.24||5.72||5.79|
|Stage II, WL %||T1–T3||54.42||56.16||57.16||55.55|
|Stage III, WL %||T3–Tmax||35.57||33.30||32.24||33.69|
|Final residue at 900-100 °C (%)||5.66||5.30||4.88||4.97|
2.5 Nonisothermal kinetic analysis
The TGA data were used to calculate the nonisothermal kinetic and thermodynamic parameters using model-free equations of Flynn-Wall-Ozawa (FWO), Kissinger-Akahila-Sunose (KAS), and Starink reported in the literature .
where is the time to reach a given extent of conversion at temperature . At , the value of is determined from the slope of the plot ), and ln(
where is the initial biomass weight, is the change in weight at a particular time during the experiment, and is the residual weight after time of the experiment.
The choice for these model-free methods is because no previous knowledge about reaction mechanism is required to determine the reaction activation energy . The preexponential factor (A) and thermodynamic parameters [enthalpy (
where is the heating rate, is the activation energy, is the maximum peak temperature, is the boltzman constant, and is the plank constant.
The activation energies for both sample methods were calculated using the three models, namely FWO, KAS, and Starink. These model-free methods avoid the shortcomings during model fitting and kinetic compensation effects. The FWO model-free method compensates the experimental measurement errors. However, the KAS and Starink methods depend on the good constant degree of conversion from the derivative mass loss function to provide precision of the kinetic data . Therefore, application of different model-free methods involves a wide conversion range that allows for the study of change in mechanism during a reaction and reduces mass transfer limitations by using multiple heating rates . Figure 5 shows the relationship of activation energy and enthalpy from the three model-free methods for acid-base and Organosolv celluloses. Results showed little or no difference between
Tables 4–6 show other thermodynamic parameters from the three model-free methods for the acid-base and the Organosolv cellulose samples, respectively. The
|α||Log A (s−1)||ΔG (kJmol−1)||ΔS (Jmol−1)||Log A (s−1)||ΔG (kJmol−1)||ΔS (Jmol−1)|
|α||Log A (s−1)||ΔG (kJmol−1)||ΔS (Jmol−1)||Log A (s−1)||ΔG (kJmol−1)||ΔS (Jmol−1)|
|α||Log A (s−1)||ΔG (kJmol−1)||ΔS (Jmol−1)||Log A (s−1)||ΔG (kJmol−1)||ΔS (Jmol−1)|
2.6 Reaction model determination
Malek method which is the commonly used approach to determine probable reaction mechanism involving heterogeneous reaction was used . The Malek method is described by the following equation.
where is rate of reaction at a given conversion, , and heating rate, β, approximates the temperature integral profile and . The x values used were in a range of 5–20 and the temperature approximation function is defined by the following equation .
The theoretical z(α) plots versus α depend on
The experimental and the fitted z(α) master plots of acid-base cellulose showed a normal distribution behavioral curve trend for all model-free methods at investigated heating rates. However, Organosolv cellulose showed a sigmoid curve skewed more to the left hand side. The correlation coefficient of acid-base method ranged between 0.9789 and 0.9884, while that of Organosolv ranged between 0.9525 and 0.9795, signifying the accuracy in the reported data. It is worth to note that both methods had best fit at 15°C/min. The data were fit with polynomial curves of n = 3 and n = 4 for Organosolv and acid-base celluloses, respectively, implying third and fourth dimension growth as described by general Avrami-Erofeev model of multidimensional nuclei and random growth reaction mechanism (. This type of mechanism is often as a result of hydration, adsorption, dissolution, and defects on the crystallite within particle size of the sample that can cause thermodynamic inhibition leading to varying activation energies . Therefore, random nucleation and growth is the most probable reaction mechanism for the pyrolysis of celluloses isolated from DPW by acid-base and Organosolv methods.
The low cost and high yield acid-base and Organosolv methods were assessed for isolation of cellulose from date palm lignocellulose waste biomass. The structural, chemical, and morphological characterizations of the isolated celluloses were studied. The nonisothermal combustion studies were investigated using three different model-free methods. The reaction mechanism was studied using Malek method.
The SEM images revealed chiral nematic orderings structures distinctive of cellulose. The change in FTIR peak intensity and the difference in the vibrational bond stretching among the isolated celluloses and between original biomass signified component removal from the lignocellulose complex. The TGA results from both methods showed one major decomposition peak assigned to cellulose in contrast to original biomass with three peaks. The results further revealed a possible difference in the degradation chemistry at higher temperature where isoconversion was higher than 0.6.
The FWO model for the acid-base method gave the lowest activation energy (99.77–173.76 kJmol−1) and the Organosolv method by KAS model gave the highest activation energy (419.63–934.49 kJmol−1). There was a strong relationship between activation energy and enthalpy, and the positive enthalpy values confirmed that endothermic reaction took place during the pyrolysis of the cellulose samples. The Gibbs’s free energy, ΔG, results revealed that Organosolv cellulose pyrolysis reaction would easily reach equilibrium, much easier in a trend of KAS > Starink> FWO models. The measure for disorder was less favorable for the acid-base method with negative entropy values in the FWO model-free method.
The reaction mechanism by Malek method was described by Avrami-Erofeev model mechanism ( for the acid-base method and ( for the Organoslv method.
The results of this study confirm the existence of multistep mechanism occurring in solid-state reactions due to variations in activation energy with the heating rates. The study provides important data information and a robust approach to understanding the cellulose pyrolysis structures and mechanisms by different isolation methods across a broad range of temperature and different heating rates.
This work was financially supported by the Emirates Centre for Energy and Environment Research, UAEU (31R107).
Conflict of interest
The authors declare no conflict of interest.
Liu C, Wang H, Karim AM, Sun J, Wang Y. Catalytic fast pyrolysis of lignocellulosic biomass. Chemical Society Reviews. 2014; 43(22):7594-7623
Wu W, Mei Y, Zhang L, Liu R, Cai J. Effective activation energies of lignocellulosic biomass pyrolysis. Energy & Fuels. 2014; 28(6):3916-3923
Liao JC, Mi L, Pontrelli S, Luo S. Fuelling the future: Microbial engineering for the production of sustainable biofuels. Nature Reviews. Microbiology. 2016; 14(5):288
Bridgwater AV. Review of fast pyrolysis of biomass and product upgrading. Biomass and Bioenergy. 2012; 38:68-94
Mohan D, Pittman CU Jr, Steele PH. Pyrolysis of wood/biomass for bio-oil: A critical review. Energy & Fuels. 2006; 20(3):848-889
Huang L, Ding T, Liu R, Cai J. Prediction of concentration profiles and theoretical yields in lignocellulosic biomass pyrolysis. Journal of Thermal Analysis and Calorimetry. 2015; 120(2):1473-1482
Müller-Hagedorn M, Bockhorn H, Krebs L, Müller U. A comparative kinetic study on the pyrolysis of three different wood species. Journal of Analytical and Applied Pyrolysis. 2003; 68:231-249
Ibbett R, Gaddipati S, Davies S, Hill S, Tucker G. The mechanisms of hydrothermal deconstruction of lignocellulose: New insights from thermal–analytical and complementary studies. Bioresource Technology. 2011; 102(19):9272-9278
Segneanu A-E, Cziple F, Vlazan P, Sfirloaga P, Grozescu I, Gherman VD. Biomass extraction methods. In: Biomass Now-Sustainable Growth and Use. IntechOpen; 2013. pp. 390-399. DOI: 10.5772/55338
Koufopanos CA, Papayannakos N, Maschio G, Lucchesi A. Modelling of the pyrolysis of biomass particles. Studies on kinetics, thermal and heat transfer effects. The Canadian journal of chemical engineering. 1991; 69(4):907-915
Wu W, Cai J, Liu R. Isoconversional kinetic analysis of distributed activation energy model processes for pyrolysis of solid fuels. Industrial and Engineering Chemistry Research. 2013; 52(40):14376-14383
Gordina NE, Prokof’ev VY, Hmylova OE, Kul’pina YN. Effect of ultrasound on the thermal behavior of the mixtures for the LTA zeolite synthesis based on metakaolin. Journal of Thermal Analysis and Calorimetry. 2017; 129(3):1415-1427
Guida MY, Bouaik H, ElMouden L, Moubarik A, Aboulkas A, Elharfi K, et al. Utilization of starink approach and avrami theory to evaluate the kinetic parameters of the pyrolysis of olive mill solid waste and olive mill wastewater. Journal of Advanced Chemical Engineering. 2016; 6:1-8
Luo H, Abu-Omar MM. Lignin extraction and catalytic upgrading from genetically modified poplar. Green Chemistry. 2018; 20(3):745-753
Xu F, Yu J, Tesso T, Dowell F, Wang D. Qualitative and quantitative analysis of lignocellulosic biomass using infrared techniques: A mini-review. Applied Energy. 2013; 104:801-809
Rohmawati B, Sya’idah FAN, Alighiri D, Eden WT. Synthesis of bioplastic-based renewable cellulose acetate from teak wood ( Tectona grandis) biowaste using glycerol-chitosan plasticizer. Oriental Journal of Chemistry. 2018; 34(4):1810
Wang G, Li W, Li B, Chen H. TG study on pyrolysis of biomass and its three components under syngas. Fuel. 2008; 87(4–5):552-558
Mohammed MAA, Salmiaton A, Azlina WW, Amran MSM. Gasification of oil palm empty fruit bunches: A characterization and kinetic study. Bioresource Technology. 2012; 110:628-636
Lapuerta MN, Hernández JJ, Rodriguez JN. Kinetics of devolatilisation of forestry wastes from thermogravimetric analysis. Biomass and Bioenergy. 2004; 27(4):385-391
Hu Y, Wang Z, Cheng X, Ma C. Non-isothermal TGA study on the combustion reaction kinetics and mechanism of low-rank coal char. RSC Advances. 2018; 8(41):22909-22916
Shuping Z, Yulong W, Mingde Y, Chun L, Junmao T. Pyrolysis characteristics and kinetics of the marine microalgae Dunaliella tertiolecta using thermogravimetric analyzer. Bioresource Technology. 2010; 101(1):359-365
Kim YS, Kim YS, Kim SH. Investigation of thermodynamic parameters in the thermal decomposition of plastic waste−waste lube oil compounds. Environmental Science & Technology. 2010; 44(13):5313-5317
Vlaev LT, Georgieva VG, Genieva SD. Products and kinetics of non-isothermal decomposition of vanadium (IV) oxide compounds. Journal of Thermal Analysis and Calorimetry. 2007; 88(3):805-812
Coker AK. Ludwig’s Applied Process Design for Chemical and Petrochemical Plants. Burlington, VT, USA: Elsevier, Gulf Professional Publishing; 2014
White JE, Catallo WJ, Legendre BL. Biomass pyrolysis kinetics: A comparative critical review with relevant agricultural residue case studies. Journal of Analytical and Applied Pyrolysis. 2011; 91(1):1-33
Málek J. The kinetic analysis of non-isothermal data. Thermochimica Acta. 1992; 200:257-269
Senum GI, Yang R. Rational approximations of the integral of the Arrhenius function. Journal of Thermal Analysis and Calorimetry. 1977; 11(3):445-447
Stanko M, Stommel M. Kinetic prediction of fast curing polyurethane resins by model-free isoconversional methods. Polymers. 2018; 10(7):698
Burnham AK, Weese RK, Weeks BL. A distributed activation energy model of thermodynamically inhibited nucleation and growth reactions and its application to the β−δ phase transition of HMX. The Journal of Physical Chemistry B. 2004; 108(50):19432-19441