This work presents a new quantitative model to predict the heat of explosion of nitroaromatic compounds using the natural bond orbital (NBO) charge and 15N？NMR chemical shifts of the nitro groups (15 N N i t r o ) as structural parameters. The values of the heat of explosion predicted for 21 nitroaromatic compounds using the model described here were compared with experimental data. The prediction ability of the model was assessed by the leave-one-out cross-validation method. The cross-validation results show that the model is significant and stable and that the predicted accuracy is within 0.146？MJ？kg？1, with an overall root mean squared error of prediction (RMSEP) below 0.183？MJ？kg？1. Strong correlations were observed between the heat of explosion and the charges (R2？=？0.9533) and 15N？NMR chemical shifts (R2？=？0.9531) of the studied compounds. In addition, the dependence of the heat of explosion on the presence of activating or deactivating groups of nitroaromatic explosives was analyzed. All calculations, including optimizations, NBO charges, and 15 N N i t r o NMR chemical shifts analyses, were performed using density functional theory (DFT) and a 6-311+G(2d,p) basis set. Based on these results, this practical quantitative model can be used as a tool in the design and development of highly energetic materials (HEM) based on nitroaromatic compounds. 1. Introduction Aromatic molecules with nitro groups are an important class of highly energetic materials (HEM). One of the most important thermodynamic properties that determine the performance of these explosives and propellants is the heat of explosion (HE). The HE is the quantity of heat liberated when an HEM undergoes detonation as a high explosive or deflagration as a low explosive. The processes of detonation and deflagration occur even in the absence of external oxygen or air because HEM contain oxygen themselves. The HE can be theoretically calculated [1] and experimentally determined [2]. The calculated value is determined by the difference between the energies of formation of the explosive components and the energies of formation of the explosion products. Experimentally, the HE is determined using a calorimetric bomb. The sample quantity is usually chosen to obtain a loading density of 0.1？g/cm3. If a powder will not explode (compounds with heat of explosion lower than 800？cal/g), a “hot” powder with a known heat of explosion is added. Thus, the HE of the sample powder can be calculated from the mixture. The calorimetric values used in this work are based on liquid H2O as a reaction product.
References
[1]
M. H. Keshavarz, “Determining heats of detonation of non-aromatic energetic compounds without considering their heats of formation,” Journal of Hazardous Materials, vol. 142, no. 1-2, pp. 54–57, 2007.
[2]
R.. Meyer, J. K？hler, and A. Homburg, Explosives, Wiley-VHC, Weinheim, Germany, 2007.
[3]
A. Delpuech and J. Cherville, “Relation entre La structure electronique et la sensibilité au choc des explosifs seconaires nitrés,” Propellants Explosives, Pyrotechnics, vol. 4, no. 6, pp. 121–128, 1979.
[4]
C. Cao and S. Gao, “Two dominant factors influencing the impact sensitivities of nitrobenzenes and saturated nitro compounds,” Journal of Physical Chemistry B, vol. 111, no. 43, pp. 12399–12402, 2007.
[5]
B. M. Rice and J. J. Hare, “A quantum mechanical investigation of the relation between impact sensitivity and the charge distribution in energetic molecules,” Journal of Physical Chemistry A, vol. 106, no. 9, pp. 1770–1783, 2002.
[6]
S. Zeman, “New aspects of initiation reactivities of energetic materials demonstrated on nitramines,” Journal of Hazardous Materials, vol. 132, no. 2-3, pp. 155–164, 2006.
[7]
M. A. Hiskey, K. R. Brower, and J. C. Oxley, “Thermal decomposition of nitrate esters,” Journal of Physical Chemistry, vol. 95, no. 10, pp. 3955–3960, 1991.
[8]
C. Zhang, “Review of the establishment of nitro group charge method and its applications,” Journal of Hazardous Materials, vol. 161, no. 1, pp. 21–28, 2009.
[9]
S. P. Hernández-Rivera and R. Infante-Castillo, “A systematic theoretical investigation of the relationship between heats of detonation and NBO charges and 15N？NMR chemical shifts of nitro groups in nitramines and nitro paraffins,” Journal of Molecular Structure, vol. 960, no. 1–3, pp. 57–62, 2010.
[10]
S. P. Hernández-Rivera and R. Infante-Castillo, “Predicting heats of explosion of nitrate esters through their NBO charges and 15N？NMR chemical shifts on the nitro groups,” Computational and Theoretical Chemistry, vol. 963, no. 2-3, pp. 279–283, 2011.
[11]
R. Infante-Castillo and S. P. Hernández-Rivera, “On the choice of optimal protocol for calculation of 13C and 15N？NMR isotropic chemical shifts in nitramine systems,” Journal of Molecular Structure, vol. 940, no. 1–3, pp. 124–128, 2010.
[12]
T. A. Keith and R. F. W. Bader, “Calculation of magnetic response properties using a continuous set of gauge transformations,” Chemical Physics Letters, vol. 210, no. 1–3, pp. 223–231, 1993.
[13]
R. Infante-Castillo and S. P. Hernández-Rivera, “Experimental and theoretical studies of the molecular structure of 1-(2-pyridinylmethyl)-2-methylbenzimidazole,” Journal of Molecular Structure, vol. 968, no. 1–3, pp. 6–12, 2010.
[14]
R. Infante-Castillo and S. P. Hernández-Rivera, “DFT studies of the structure and vibrational and NMR spectra of 1-(2-methylpropenyl)-2-methylbenzimidazole,” Journal of Molecular Structure, vol. 917, no. 2-3, pp. 158–163, 2009.
[15]
B. M. Rice and J. Hare, “Predicting heats of detonation using quantum mechanical calculations,” Thermochimica Acta, vol. 384, no. 1-2, pp. 377–391, 2002.
[16]
M. H. Keshavarz, “Simple procedure for determining heats of detonation,” Thermochimica Acta, vol. 428, no. 1-2, pp. 95–99, 2005.
[17]
M. H. Keshavarz, “A simple way to predict heats of detonation of energetic compounds only from their molecular structures,” Propellants, Explosives, Pyrotechnics, vol. 37, no. 1, pp. 93–99, 2012.
[18]
M. H. Keshavarz, “Quick estimation of heats of detonation of aromatic energetic compounds from structural parameters,” Journal of Hazardous Materials, vol. 143, no. 1-2, pp. 549–554, 2007.
[19]
M. H. Keshavarz, “Estimating heats of detonation and detonation velocities of aromatic energetic compounds,” Propellants, Explosives, Pyrotechnics, vol. 33, no. 6, pp. 448–453, 2008.
[20]
M. H. Keshavarz, “Theoretical prediction of condensed phase heat of formation of nitramines, nitrate esters, nitroaliphatics and related energetic compounds,” Journal of Hazardous Materials, vol. 136, no. 2, pp. 145–150, 2006.
[21]
V. I. Pepekin and S. A. Gubin, “Heat of explosion of commercial and brisant high explosives,” Combustion, Explosion and Shock Waves, vol. 43, no. 2, pp. 212–218, 2007.
[22]
L. M. Hall, L. H. Hall, and L. B. Kier, “QSAR modeling of β-lactam binding to human serum proteins,” Journal of Computer-Aided Molecular Design, vol. 17, no. 2–4, pp. 103–118, 2003.
[23]
M. J. Frisch, G. W. Trucks, H. B. Schlegel, et al., Gaussian 09 Rev. A.1., CT. Gaussian, Wallingford, UK, 2009.
[24]
A. D. Becke, “Density-functional thermochemistry. III. The role of exact exchange,” Journal of Chemical Physics, vol. 98, no. 7, pp. 5648–5652, 1993.
[25]
C. Lee, W. Yang, and R. G. Parr, “Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density,” Physical Review B, vol. 37, no. 2, pp. 785–789, 1988.
[26]
E. D. Glendening, A. E. Reed, J. E. Carpenter, and F. Weinhold, NBO Version 3.1, 2009.
[27]
B. Mennucci and J. Tomasi, “Continuum solvation models: a new approach to the problem of solute's charge distribution and cavity boundaries,” Journal of Chemical Physics, vol. 106, no. 12, pp. 5151–5158, 1997.
[28]
S. Zeman, “Modified Evans-Polanyi-Semenov relationship in the study of chemical micromechanism governing detonation initiation of individual energetic materials,” Thermochimica Acta, vol. 384, no. 1-2, pp. 137–154, 2002.
[29]
S. Zeman, “The relationship between differential thermal analysis data and the detonation characteristics of polynitroaromatic compounds,” Thermochimica Acta, vol. 41, no. 1, pp. 199–212, 1980.
