Toward a Physically Based Quantitative Modeling of Impact

Feb 14, 2013 - Mariela Inés Profeta , Jorge Marcelo Romero , Nelly Lidia Jorge , André Grand , Alfonso Hernández-Laguna. Journal of Molecular Model...
0 downloads 0 Views 536KB Size
Article pubs.acs.org/JPCA

Toward a Physically Based Quantitative Modeling of Impact Sensitivities Didier Mathieu* CEA, DAM, Le Ripault, F-37260 Monts, France S Supporting Information *

ABSTRACT: Among the subsequent steps leading from impact to explosive decomposition in nitro compounds, the ability of early exothermic reactions to trigger the decomposition of neighboring molecules before the released energy has dissipated away is assumed to be critical. The rate of this process is roughly estimated using as inputs the energy content and the dissociation energy of the weakest X−NO2 bonds. While the sensitivity index thus obtained was previously shown to exhibit striking correlations with gap test pressures, its correlation with drop weight impact test data is poorer. Nevertheless, considering four different subsets of molecules depending on the environment of the most labile nitro groups, straightforward regressions against this sensitivity index yield a practical method to estimate impact sensitivity, whose combination of fair performance and generality is provided by no alternative approach, except purely empirical models based on extensive parametrization. compound and equivalent in practice to OB.8 Over the years, alternative single-variable correlations have been put forward, based on a variety of molecular properties including electronic band gap,17 electronic shakeup promotion energies,18 13C and 15 N NMR chemical shifts,19 molecular electronegativity,20 X− NO2 bond lengths,21 total volumic energy content,22 a parameter derived from the hard and soft acid and bases (HSAB) principle,23 etc. As emphasized by Brill and James for nitroaromatic compounds, such studies based on small data sets are prone to chance correlations and are often misleading.24,25 This problem is even more pregnant for multivariables correlations, especially for those that consider only small data sets. Features related to electron distributions26,27 or molecular electrostatic potentials28−30 are especially popular descriptors. Some authors focus on the electronic structure of covalent X− NO2 bonds (called explosophore bonds).31,32 Some correlations based on constitutive descriptors, in the spirit of group contribution methods, are also available.33,34 Current research increasingly resorts to advanced quantitative structure− property relationships (QSPR) methodology and regression techniques beyond straightforward multilinear regression (MLR), such as partial least-squares (PLS) or artificial neural networks (ANN).35−41 While the use of large data sets should decrease the risk of chance correlations, the fact that input variables are typically selected out of extended pools of descriptors and the flexibility afforded by nonlinear regression techniques does increase this risk. As a result, it is not easy to figure out the true predictive value of such QSPR models.

1. INTRODUCTION The impact sensitivity of energetic materials is most often characterized by the height H50 that a given weight must be dropped onto the sample to trigger an observable decomposition with a 50% probability.1,2 The H50 values measured in this so-called drop weight impact test depend on many factors beyond chemical structure, some of them poorly controlled, such as the detailed microstructure and defects of the sample or the environment of the experiment.3−6 As a result, they exhibit large uncertainties. To some extent, shock sensitivity, as characterized by gap test threshold pressures, provides a more reliable index of safety.7 Both sensitivity tests involve completely different initiation regimes, namely, reaction times in the range of 200−250 μs and pressures between 7−15 kbar for the drop weight impact test, versus corresponding ranges of 0.05−2 μs and 30−200 kbar for the gap test.8 Nevertheless, significant correlations between H50 and gap-test pressures are observed in practice.8 Being more easily amenable to measurements, H50 remains therefore the most popular criterion to characterize mechanical sensitivities of new compounds. A general procedure to evaluate H50 from molecular structure would therefore be of interest not only as a guide to design new explosives9,10 but also to assess reactive hazards associated with organic compounds in the broader context of emerging safety regulations concerning the chemical industry.11,12 For many years, much effort has been devoted to the development of such structure−sensitivity relationships. Early approaches are mainly empirical. Half a century ago, experimental evidence that oxygen rich compounds are especially sensitive stimulated the development of linear correlations between log 10 (H 50 ) and oxygen balance (OB).13−16 In 1990, Storm et al. introduced a sensitivity index defined on the basis of the empirical formula of the © 2013 American Chemical Society

Received: November 27, 2012 Revised: February 14, 2013 Published: February 14, 2013 2253

dx.doi.org/10.1021/jp311677s | J. Phys. Chem. A 2013, 117, 2253−2259

The Journal of Physical Chemistry A

Article

models for H50 have already been reported, they usually fail in demonstrating reasonable predictions for an extended test set of explosives. Not surprisingly, the best correlations with experiment are obtained using state-of-the-art QSPR techniques. The most impressive results to date are probably those reported very recently by Xu et al.41 Using an ANN model with input variables automatically picked out of an extended pool of descriptors, these authors obtain R2 = 0.85 for a training set of 127 explosives and R2 = 0.87 for a test set of 29 explosives. Using a very similar approach, Morill and Byrd obtain R2 = 0.69 for a test set of seven molecules.39 This discrepancy illustrates the difficulty in assessing the performance of such purely empirical models. Notwithstanding the recent ANN of Xu et al.,41 there are only very few models with R2 > 0.80.66,69 These correlations are obtained for small data sets with typically about 10 compounds. Notwithstanding the above-mentioned overfitting issues, the quality of such correlations rapidly decreases as the size/diversity of the data set increases. For instance, a carefully designed model based on charge distributions in molecules yields R2 = 0.54 for the 14 compounds in the test set.26 In fact, notwithstanding the Xu et al. results, the highest R2 values obtained are typically close to 0.70 as extended data sets (with N > 100 compounds) are considered. In this work, we focus on a more physically grounded approach that provides similar performance without resorting to an empirical search for optimal input variables and analytical expression. While not as accurate as the Xu et al. method, it yields more valuable insight into relevant aspects of initiation and the role of molecular parameters.

Considering the limited amount of experimental data at hand, a more physically grounded approach is desirable, even if a fully rigorous and quantitative model will be hard to obtain considering the complex phenomena involved. Indeed, an impact wave propagating into a material interacts with defects such as voids, shear bands, or interfaces through various mechanical processes.42 In order to trigger an explosion, these processes must be able to focus the mechanical energy in localized regions of the material and lead to the formation of hot spots. In the simplest picture, hydrodynamic processes around localized defects are sufficient to generate hot spots where temperature is high enough to initiate thermal bond scission reactions. Assuming such processes to be similar for most organic explosives, one might expect H50 to depend on the decomposition temperature Tdec. In fact, correlations between H50 and Tdec are observed only for the most stable explosives and exhibit significant exceptions.43 Therefore, there is a need to better understand the role of hydrodynamic aspects, microstructure and defects in typical granular explosives,44,45 and/or to identify alternative molecular mechanisms. In the seventies, it was proposed that the decomposition might be a result of shock induced electronic excitations.46,47 Theoretical justifications for such processes were subsequently developed.48−50 In fact, calculated probabilities of mechanically induced electronic transitions are significantly overestimated unless they account for the decoherence of the electronic wavefuntion induced by its interaction with the macroscopic sample.51 Quantitative estimates can only be obtained under the assumption that the electron−phonon coupling time scale is shorter than the reaction time scale.52 This condition implies that the approach is best satisfied for inorganic primary explosives, such as metal azides for which some limited correlations have been reported.17 Other attempts to account for the first steps in the decomposition process include the viscoplastic collapse of the porosity,53 jetting or vaporization of material into voids,54,55 pressure-induced proton tunneling,56 mechanical generation of charge transfer pairs that localize at vacancies,57 or intermolecular metastable trigger reactions leading to early chemical energy release.58 However, these assumed mechanisms did not lead to any general quantitative model that would allow to correlate observed H50 data with microscopic features. Alternatively, it has been assumed that the transfer of the mechanical energy from the phonon bath to intramolecular vibrations is the most determinant process influencing observed sensitivities. This energy transfer depends essentially on the excitation of doorway modes, i.e., of the molecular vibrations of lowest frequency. The corresponding rate can be theoretically estimated using standard approximations.59,60 To some extent, H50 data correlates with values thus calculated61 or with related theoretical quantities.62−65 However, such correlations rely only on small sets of compounds. The alternative assumption that sensitivities depend primarily on the energy required to trigger the first bond scissions yields similar results.66−68 The performance of predictive models for H50 is often characterized in terms of the determination coefficient R2 and root-mean-square deviation (rms) between observed and predicted (rather than f itted) log10(H50) values. Considering H50 directly leads to ill-defined statistics as sensitivity data sets typically exhibit a few insensitive compounds with H50 values in the range 300−500 cm, while most explosives exhibit H50 values N−NO2 and −NH−NO2. Finally, ΔG0diss is assumed to exhibit the same value for all O− NO2 bonds. With such transferability assumptions, only 11 ΔG0diss values are needed for the bonds present in our data set. With the help of the ORCA program,76 these 11 values are computed within the harmonic approximation using the model compounds listed in Table 1. HF/3-21G geometries and

as no thermal equilibrium is assumed. It is expressed here as a constant fraction η of the chemical energy content per atom, an approximation that implies that all compounds considered exhibit similar values of their diffusion coefficient. However, the activation energy E†pr for the propagation increases with the strength of the chemical bonds to be broken for new molecules to decompose, hence on the bond dissociation energy (BDE) of the weakest X−NO2 bonds in unreacted molecules, hereafter denoted Dmin. Since E†pr = Dmin for simple gas phase X−NO2 bond cleavages, a scaling coefficient λ is introduced to account for the role of the environment of the nitro group: E†pr = λDmin. For simplicity, λ is first assumed to be constant for all explosives considered. However, it will prove necessary to introduce different values depending on the nature of the weakest X− NO2 bond, as described in section 3. Under these assumptions, the propagation rate kpr is given by ⎛ E† ⎞ ⎞ ⎛ Dmin pr ⎟ = Zpr exp⎜ − λ ⎟ k pr = Zpr exp⎜⎜ − 0 ⎟ ⎝ η 2Δd H /3N ⎠ ⎝ kBT ⎠

Table 1. Definition of Standard X−NO2 Bond Types and Corresponding ΔG0diss Values

(1)

where kB is the Boltzmann constant, ΔdH0 the energy released by the decomposition of one mole of the compound, N the corresponding number of atoms, and Zpr the prefactor of the process. Since the impact energy that eventually leads to decomposition is proportional to H50, it is natural to assume that log10(H50) linearly correlates with log10(kpr), hence with the shock sensitivity index SI defined as follows SI = NDmin /Δd H 0

(2)

and previously introduced to rationalize gap test data.70,71 This linear relationship exhibits a slope equal to 3λ/2η. Assuming λ close to unity, we can estimate the fraction η of the chemical energy that efficiently contributes to propagate the decomposition. In the present work, the decomposition enthalpy ΔdH0 is obtained as the difference between the formation enthalpy ΔfH0 of the explosive compound and corresponding enthalpies for the decomposition products. The latter are estimated with the help of the H2O−CO2 arbitrary.72 In other words, it is assumed that the main decomposition products are, in order of decreasing formation priority, H2O, CO2, and, of course, N2. For the initial molecule, ΔfH0 is approximated by the value derived from the RM1 Hamiltonian73 as implemented in the MOPAC7 program.74 For the products, standard experimental enthalpies compiled in the NIST database are used.75 For simplicity, Zpr is assumed constant. In that case, E†pr should in principle be replaced in eq 1 by the Gibbs energy of activation ΔG†pr for the propagation of the decomposition. There is no way to rigorously evaluate this quantity. However, in an attempt to capture at least partially some entropic contributions, standard Gibbs free energies of dissociation ΔG0diss are used. The systematic calculation of ΔG0diss for all X− NO2 bonds would be tedious, given the need for many calculations on large radicals. Therefore, the present work resorts to a rough approximation, assuming the transferability of a limited number of ΔG0diss values. For aromatic C−NO2 bonds, ΔG0diss is assumed to depend only on the presence of proton donor groups in adjacent positions on the aromatic ring, as in such cases NO2 oxygens are involved in hydrogen bonding. For nitro groups bonded to sp3 carbons, the energy of the C−NO2 bond is assumed to depend only on the number of nitro groups in geminal positions. For nitramines, two distincts

frequencies are used. The frequencies are scaled using the standard 0.9207 coefficient.77 Total electronic energies are calculated at the B3LYP level using the def2-TZVP basis set,78 the chain of spheres acceleration algorithm,79 and van der Waals corrections.80 It must be emphasized that the present approach relies exclusively on the energy content and bond strengths of isolated molecules. As such, it cannot account for the dependence of the impact sensitivity on such parameters as defects, polymorphism, or the crystal orientation with respect to the shock direction.81 However, this dependence may be addressed with approaches complementary to the present one, such as molecular dynamics simulations based on reactive force fields.82,83

3. RESULTS The X−NO2 dissociation energies obtained in this work are listed in Table 1. As expected, they follow the usual stability order, namely, CA−NO2 > N−NO2 > O−NO2 with CA standing for an aromatic carbon. For nitroalkyls, isolated C− NO2 bonds are somewhat stronger than N−NO2 bonds, while geminal C−NO2 bonds are no stronger than O−NO2 bonds. Despite the fact that X−NO2 bonds in similar chemical 2255

dx.doi.org/10.1021/jp311677s | J. Phys. Chem. A 2013, 117, 2253−2259

The Journal of Physical Chemistry A

Article

this work leave much room for further improvement, e.g., through the use of more accurate BDE data. The poor correlation obtained for nitramines (set S2) stems from the fact that cyclotrimethylene trinitramine (RDX) and cyclotetramethylene tetranitramine (HMX) prove less sensitive than estimated on the basis of present approximate SI values, while the opposite is observed for N,N′-dinitro-methanediamine O2N−NH−CH2−NH−NO2, which is the smallest nitramine in the data set. The fact that these three outliers do not fit in the present correlation does not appear to stem from the use of standard bond energies from Table 1. Indeed, assuming transferable ΔG0diss values prove reasonable for these compounds. For instance, a value of 26.43 kcal/mol is computed for RDX, instead of the assumed value of 28.61 kcal/mol. Similarly, the value of 34.5 kcal/mol calculated for N,N′-dinitro-methanediamine is reasonably close to the assumed value of 32.98 kcal/mol. For NACs (set S1), the correlation is significantly spoiled by the fact that 2,4,6-trinitrobenzene-1,3,5-triol proves significantly more sensitive than expected on the basis of the corresponding sensitivity index. In contrast, other compounds with strong intramolecular hydrogen bonds, including triaminotrinitrobenzene (TATB), fit nicely into the present scheme. Therefore, there is no obvious explanation for the underestimation of the sensitivity of this compound. Finally, the correlation for set S4 is especially sensitive to uncertainties associated with present SI data since all SI values are very similar (in the range 1−1.3) in this data set, while H50 values range from 1 implies that, as the critical step in the propagation of a decomposition reaction to the surrounding material is attained, a regime where chemical energy is released faster than it diffuses away is already reached. Therefore, values of η larger 2256

dx.doi.org/10.1021/jp311677s | J. Phys. Chem. A 2013, 117, 2253−2259

The Journal of Physical Chemistry A

Article

than unity probably arise as a result of present approximations, especially the fact that the prefactor in eq 1 is assumed constant, while ΔG†pr is simply derived from ΔG0diss. Nevertheless, the result that η = 3.4 for NACs is puzzling as long as ΔG†pr ≃ Dmin is assumed. In fact, although this approximation is reasonable for nitramines, nitric esters, and nitroaliphatic compounds, for which the homolysis of X−NO2 bonds is believed to be the dominant process determining the stability of the compound, it might be misleading when applied to NACs for which preliminary rearrangements of ring substituents are often energetically favored.25 In such cases, activation energies significantly lower than Dmin would account for the large value of η presently obtained for NACs. Compared with the Storm et al. sensitivity index defined on the basis of the empirical formula,8 the present one is more physically grounded.70 Figure 2 shows how presently calculated

Table 2. Performance Comparison between Present Linear Regressions and Recent Models Based on MLR, PLS, or ANN Techniques (See Text for Details) training set 2

2

test set 2

model

ref

R (fit)

R (LOO)

rms

R (ext)

rms

MLRa MLRc MLRc PLSa ANNa ANNb ANNc this work

40 41 39 40 40 40 41

0.771 0.768 0.814 0.766 0.816 0.795 0.848 0.710

0.593 0.703 0.795 0.674

0.212 0.194

0.715 0.722 0.691 0.718 0.740 0.756 0.866 0.721

0.251 0.177 0.285 0.250 0.247 0.257 0.130 0.265

0.846 0.658

0.214 0.192 0.203 0.164 0.238

a

Using electrotopological state indices. bUsing the same descriptors as in ref 35. cSelecting variables out of a large pool of descriptors.

combined with automatic variable selection from extended pools of descriptors, as done by Xu et al.41 Unfortunately, such techniques provide good results at the expense of physical insight.

4. CONCLUSIONS AND PERSPECTIVES Although present equations might not be as reliable as neural network models, the present sensitivity index opens new perspectives for the a priori evaluation of drop weight impact test data. A major advantage of the present approach over empirical QSPR methodologies lies in the fact that it can be systematically improved through removing simplifying assumptions. For instance, the present focus on the weakest bond in the compound is clearly a rough approximation. Accounting for all X−NO2 bonds should yield significant improvement. This will be the next step in our present approach. However, in the context of QSPR modeling, the present sensitivity index appears as a promising descriptor in view of its ability to fit a significant body of drop-weight impact test data on the basis of a straightforward and physically motivated partition of the data set. It might be combined with descriptors taking account secondary aspects of impact sensitivity, such as intermolecular interactions, particle size distribution, etc.

Figure 2. Calculated versus measured impact sensitivities for the 156 compounds in the database. Same color code as that in Figure 1. Empty symbols are for compounds in the training set, and filled symbols are for compounds in the test set.



impact sensitivities compare to experiment. The fact that they are not very accurate is to be expected in view of present approximations and keeping in mind experimental uncertainties. However, it is gratifying to observe that the present model does not lead to any significant outlier, suggesting that SI is indeed the primary factor determining observed sensitivities. Table 2 compares the performances of present correlations with the outcome of recent models taking advantage of state-ofthe-art machine learning techniques. R2(fit), R2(ext), and R2(LOO) are the determination coefficients between measured log 10 (H 50 ) values and calculated counterparts derived, respectively, from the training set, the test set, and the leaveone-out (LOO) cross-validation procedure. Values of the rootmean-square error (rms) are provided for both the training and the test sets. For every model, these subsets contain, respectively, 127 and 29 compounds from the present data set, except for the model from ref 39, which used a somewhat larger data set. Although present regressions rely on relatively poor correlations, their predictive value compares well with the best empirical linear models, according to statistics reported in this table for the external test set. In fact, more accurate predictions can only be obtained using neural networks

ASSOCIATED CONTENT

S Supporting Information *

Table with experimental and calculated log10(H50/cm) values for the 156 explosives considered in this work, referred to by their index as defined in ref 43. For every compound, the corresponding molecular subset (S1 to S4) is indicated, as well as a letter T or V depending on whether this compound belongs to the training or the test (validation) set. This material is available free of charge via the Internet at http://pubs.acs.org.



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Notes

The authors declare no competing financial interest.



REFERENCES

(1) Akhavan, J. The Chemistry of Explosives, 2nd ed.; The Royal Society of Chemistry: Cambridge, U.K., 2004. (2) Fried, L. E.; Manaa, M. R.; Pagoria, P. F.; Simpson, R. L. Design and Synthesis of Energetic Materials. Annu. Rev. Mater. Res. 2001, 31, 291−321. 2257

dx.doi.org/10.1021/jp311677s | J. Phys. Chem. A 2013, 117, 2253−2259

The Journal of Physical Chemistry A

Article

(25) Brill, T. B.; James, K. J. Kinetics and Mechanisms of Thermal Decomposition of Nitroaromatic Explosives. Chem. Rev. 1993, 93, 2667−2692. (26) Rice, B. M.; Hare, J. J. A Quantum Mechanical Investigation of the Relation between Impact Sensitivity and the Charge Distribution in Energetic Molecules. J. Phys. Chem. A 2002, 106, 1770−1783. (27) Zhang, C. Review of the Establisment of Nitro Group Charge Method and Its Application. J. Hazard. Mater. 2009, 161, 21−28. (28) Murray, J. S.; Lane, P.; Politzer, P.; Bolduc, P. R. A Relationship between Impact Sensitivity and the Electrostatic Potentials at the Midpoints of C-NO2 Bonds in Nitroaromatics. Chem. Phys. Lett. 1990, 168, 135−139. (29) Murray, J. S.; Lane, P.; Politzer, P. Relationships Between Impact Sensitivities and Molecular-Surface Electrostatic Potentials of Nitroaromatic and Nitroheterocyclic Molecules. Mol. Phys. 1995, 85, 1−8. (30) Murray, J. S.; Concha, M. C.; Politzer, P. Links Between Surface Electrostatic Potentials of Energetic Molecules, Impact Sensitivities and C−NO2/N−NO2 Bond Dissociation Energies. Mol. Phys. 2009, 107, 89−97. (31) Kohno, Y.; Maekawa, K.; Tsuchioka, T.; Hashizume, T.; Imamura, A. A Relationship between the Impact Sensitivity and the Electronic Structures for the Unique N−N Bond in the HMX Polymorphs. Combust. Flame 1994, 96, 343−350. (32) Kohno, Y.; Ueda, K.; Imamura, A. Molecular Dynamics Simulations of Initial Decomposition Process on the Unique N−N Bond in Nitramines in the Crystalline State. J. Phys. Chem. 1996, 100, 4701−4712. (33) Keshavarz, M. H. Simple Relationship for Predicting Impact Sensitivity of Nitroaromatics, Nitramines, and Nitroaliphatics. Propellants, Explos., Pyrotech. 2010, 35, 175−181. (34) Lai, W.-P.; Lian, P.; Wang, B.-Z.; Ge, Z.-X. New Correlations for Predicting Impact Sensitivities of Nitro Energetic Compounds. J. Energ. Mater. 2010, 28, 45−76. (35) Nefati, H.; Cense, J.-M.; Legendre, J.-J. Prediction of the Impact Sensitivity by Neural Networks. J. Chem. Inf. Comput. Sci. 1996, 36, 804−810. (36) Cho, S. G.; No, K. T.; Goh, E. M.; Kim, J. K.; Shin, J. H.; Joo, Y. D.; Seong, S. Optimization of Neural Networks Architecture for Impact Sensitivity of Energetic Molecules. Bull. Korean Chem. Soc. 2005, 26, 399−408. (37) Keshavarz, M. H.; Jaafari, M. Investigation of the Various Structure Parameters for Predicting Impact Sensitivity of Energetic Molecules via Artificial Neural Networks. Propellants, Explos., Pyrotech. 2006, 31, 216−225. (38) Jun, Z.; Xin-Iu, C.; Bi, H.; Xiang-Dong, Y. Neural Networks Study on the Correlation Between Impact Sensitivity and Molecular Structures for Nitramine Explosives. Struct. Chem. 2006, 17, 501−507. (39) Morill, J. A.; Byrd, E. F. C. Development of Quantitative Structure−Property Relationships for Predictive Modeling and Design of Energetic Materials. J. Mol. Graphics Model. 2008, 27, 349−355. (40) Wang, R.; Jiang, J.; Pan, Y.; Cao, H.; Cui, Y. Prediction of Impact Sensitivity of Nitro Energetic Compounds by Neural Network based on Electrotopological-State Indices. J. Hazard. Mater. 2009, 166, 155−186. (41) Xu, J.; Zhu, L.; Fang, D.; Wang, L.; Xiao, S.; Liu, L.; Xu, W. QSPR Studies of Impact Sensitivity of Nitro Energetic Compounds Using Three-Dimensional Descriptors. J. Mol. Graphics Model. 2012, 36, 10−19. (42) Bowden, F. P.; Yoffe, A. D. Initiation and Growth of Explosion in Liquids and Solids; Cambridge University Press: Cambridge, U.K., 1952. (43) Mathieu, D.; Beaucamp, S. Matériaux Energétiques; Encyclopédie ″Techniques de l’Ingénieur″, Traité ″Sciences Fondamentales″; Fascicule AF 6710: Paris, France, 2004. (44) Gonthier, K. A. Predictions for Weak Mechanical Ignition of Strain Hardened Granular Explosive. J. Appl. Phys. 2004, 95, 3482− 3494.

(3) van der Steen, A. C.; Verbeek, H. J.; Meulenbrugge, J. J. In Influence of RDX Crystal Shape on the Shock Sensitivity of PBXs. Proceedings of the 9th International Symposium on Detonation, Portland, Oregon, 1989; pp 83−88. (4) Armstrong, R. W.; Coffey, C. S.; DeVost, V. F.; Elban, W. L. Crystal Size Dependence for Impact in Initiation of Cyclotrimethylenetrinitramine Explosive. J. Appl. Phys. 1990, 68, 979−984. (5) Handley, C. A.; Lambourn, B. D. Predicting the Effect of Porosity on the Shock Sensitivity of Explosives. Shock Compression of Condensed Matter 2009: Proceedings of the American Physical Society Topical Group on Shock Compression of Condensed Matter, Nashville, TN, 28 June−3 July 2009; American Institute of Physics: Melville, New York, 2009; pp 221−224. (6) Bellitto, V. J.; Melnik, M. I. Surface Defects and their Role in the Shock Sensitivity of Cyclotrimethylene-Trinitramine. Appl. Surf. Sci. 2010, 256, 3478−3481. (7) Tan, B.; Long, X.; Peng, R.; Li, H.; Jin, B.; Chu, S. On the Shock Sensitivity of Explosive Compounds With Small-Scale Gap Test. J. Phys. Chem. A 2011, 115, 10610−10616. (8) Storm, C. B.; Stine, J. R.; Kramer, J. F. Sensitivity Relationships in Energetic Materials. In Chemistry and Physics of Energetic Materials; Bulusu, S. N., Ed.; NATO ASI Series C; Kluwer: Dordrecht, The Netherlands, 1990; Vol. 309, pp 605−639. (9) Rice, B. M.; Byrd, E. F. C. Theoretical Chemical Characterization of Energetic Materials. J. Mater. Res. 2006, 21, 2444−2452. (10) Zeman, S. High Energy Density Materials; Structure and Bonding; Springer-Verlag: Heidelberg, Germany, 2007; Vol. 125; pp 195−271. (11) Warhurst, M. A. REACH, a New Approach to Chemicals Regulation in Europe: a Brief History, Key Features, and Expected Outcomes. J. Eur. Environ. Planning Law 2005, 2, 164−172. (12) Silk, J. C. Development of a Globally Harmonized System for Hazard Communication. Int. J. Hyg. Environ. Health 2003, 206, 447− 452. (13) Kamlet, M. J. Third Symposium (International) on Detonation; Office of Naval Research: Washington, DC, 1960; pp 13−34. (14) Kamlet, M. J. Sixth Symposium (International) on Detonation; Office of Naval Research: Washington, DC, 1976; pp 312−322. (15) Kamlet, M. J.; Adolph, H. G. The Relationship of Impact Sensitivity with Structure of Organic High Explosives. II. Polynitroaromatic Explosives. Propellants, Explos., Pyrotech. 1979, 4, 30−34. (16) Kamlet, M. J.; Adolph, H. G. Some Comments Regarding the Sensitivities. Thermal Stabilities, and Explosive Performance Characteristics of Fluorodinitromethyl Compounds 1981, 84−92. (17) Zhu, W.; Xiao, H. First-Principles Band Gap Criterion for Impact Sensitivity of Energetic Crystals: a Review. Struct. Chem. 2010, 21, 657−665. (18) Sharma, J.; Beard, B. C.; Chaykovsky, M. Correlation of Impact Sensitivity With Electronic Levels and Structure of Molecules. J. Phys. Chem. 1991, 95, 1209−1213. (19) Zeman, S. In Energetic Materials; Politzer, P., Murray, J. S., Eds.; Elsevier: Amsterdam, The Netherlands, 2003; Vol. 2; pp 25−52. (20) Mullay, J. A Relationship Between Impact Sensitivity and Molecular Electronegativity. Propellants, Explos., Pyrotech. 1987, 12, 60−63. (21) Politzer, P.; Murray, J. S.; Lane, P.; Sjoberg, P.; Adolph, H. G. Shock-Sensitivity Relationships for Nitramines and Nitroaliphatics. Chem. Phys. Lett. 1991, 181, 78−82. (22) Pepekin, V. I.; Korsunskii, B. L.; Denisaev, A. A. Initiation of Solid Explosives by Mechanical Impact. Combust., Explos. Shock Waves 2008, 44, 586−590. (23) Koch, E. C. Acid-Base Interactions in Energetic Materials: I. The Hard and Soft Acids and Bases (HSAB) Principle − Insights to Reactivity and Sensitivity of Energetic Materials. Propellants, Explos., Pyrotech. 2005, 30, 5−16. (24) Brill, T. B.; James, K. J. Thermal Decomposition of Energetic Materials. 61. Perfidy in the Amino-2,4,6-Trinitribenzene Series of Explosives. J. Phys. Chem. 1993, 97, 8752−8758. 2258

dx.doi.org/10.1021/jp311677s | J. Phys. Chem. A 2013, 117, 2253−2259

The Journal of Physical Chemistry A

Article

(45) Hamate, Y.; Horie, Y. Ignition and Detonation of Solid Explosives: a Micromechanical Burn Model. Shock Waves 2006, 16, 125−147. (46) Delpuech, A.; Cherville, J. Relation between Shock Sensitiveness of Secondary Explosives and their Molecular Electronic Structure. 1. Nitroaromatics and Nitramines. Propellants, Explos., Pyrotech. 1978, 3, 169. (47) Delpuech, A.; Cherville, J. Relation between Shock Sensitiveness of Secondary Explosives and their Molecular Electronic Structure. 2. Nitrate Esters. Propellants, Explos., Pyrotech.. 1979, 4, 121. (48) Mathieu, D.; Martin, P. Molecular Dynamics Simulation of Shockwaves Including some Nonadiabatic Effects. Comput. Mater. Sci. 2000, 17, 347−351. (49) Kuklja, M. M. An Excitonic Mechanism of Detonation Initiation in Explosives. J. Chem. Phys. 2000, 112, 3417−3423. (50) Kuklja, M. M.; Aduev, B. P.; Aluker, E. D.; Krasheninin, V. I.; Krechetov, A. G.; Mitrofanov, A. Y. Role of Electronic Excitations in Explosive Decomposition of Solids. J. Appl. Phys. 2001, 89, 4156− 4166. (51) Mathieu, D.; Martin, P.; La Hargue, J.-P. Simulation of the Electron Dynamics in Shockwaves and Implications for the Sensitivity of Energetic Materials. Phys. Scr., T 2005, 118, 171−173. (52) Reed, E. J. Electron-Ion Coupling in Shocked Energetic Materials. J. Phys. Chem. C 2012, 116, 2205−2211. (53) Belmas, R.; Plotard, J.-P. Physical Origin of Hot Spots in Pressed Explosive Compositions. J. Phys. IV 1995, 5, 61−87. (54) Holian, B. L.; Germann, T. C.; Maillet, J.-B.; White, C. T. Atomistic Mechanism for Hot Spot Initiation. Phys. Rev. Lett. 2002, 89, 285501. (55) Holian, B. L.; Germann, T. C.; Maillet, J.-B.; White, C. T. Atomistic Mechanism for Hot Spot Initiation: Erratum. Phys. Rev. Lett. 2003, 90, 069902. (56) Méreau, R.; Mathieu, D.; Elstner, M.; Frauenheim, T. Theoretical Investigation of the High-Pressure Behavior of Nitric Acid. Phys. Rev. B 2004, 69, 104101. (57) Tsiaouris, D.; Munn, R. W. Energy of Charged States in the RDX Crystal: Trapping of Charge-Transfer Pairs As a Possible Mechanism for Initiating Detonation. J. Chem. Phys. 2005, 122, 184708. (58) Ramaswamy, A. L. Mesoscopic Approach to Energetic Material Sensitivity. J. Energ. Mater. 2006, 24, 35−65. (59) Dlott, D. D.; Fayer, M. D. Shocked Molecular Solids: Vibrational Up Pumping, Defect Hot Spot Formation, and the Onset of Chemistry. J. Chem. Phys. 1990, 92, 3798−3811. (60) Tokmakoff, A.; Fayer, M. D.; Dlott, D. D. Chemical Reaction Initiation and Hot-Spot Formation in Shocked Energetic Molecular Materials. J. Phys. Chem. 1993, 97, 1901−1913. (61) Fried, L. E.; Ruggiero, A. J. Energy Tranfer Rates in Primary, Secondary, and Insensitive Explosives. J. Phys. Chem. 1994, 98, 9786− 9791. (62) McNesby, K. L.; Coffey, C. S. Spectroscopic Determination of Impact Sensitivities of Explosives. J. Phys. Chem. B 1997, 101, 3097− 3104. (63) Ye, S.; Tonokura, K.; Koshi, M. Energy Transfer Rates and Impact Sensitivities of Crystalline Explosives. Combust. Flame 2003, 132, 240−246. (64) Ye, S.; Koshi, M. Theoretical Studies of Energy Transfer Rates of Secondary Explosives. J. Phys. Chem. B 2006, 110, 18515−18520. (65) Su-Hong, G.; Xin-Lu, C.; Li-Sha, W.; Xiang-Dong, Y. Correlation Between Normal Mode Vibrations and Impact Sensitivities of Some Secondary Explosives. J. Mol. Struct. 2007, 809, 55−60. (66) Rice, B. M.; Sahu, S.; Owens, F. J. Density Functional Calculations of Bond Dissociation Energies for NO2 Scission in Some Nitroaromatic Molecules. J. Mol. Struct. 2002, 583, 69−72. (67) Song, X.-S.; Cheng, X.-L.; Yang, X.-D.; He, B. Relationship between the Bond Dissociation Energies and Impact Sensitivities of Some Nitro-Explosives. Propellants, Explos., Pyrotech. 2006, 31, 306− 309.

(68) Jun, Z.; Da-Hai, X.; Xin-Iu, C. Investigation of Correlation between Impact Sensitivities and Bond Dissociation Energies in some Triazole Energetic Compounds. Struct. Chem. 2010, 21, 1235−1240. (69) Anders, G.; Borges, I., Jr. Topological Analysis of the Molecular Charge Density and Impact Sensitivity Models of Energetic Molecules. J. Phys. Chem. A 2011, 115, 9055−9068. (70) Mathieu, D. Theoretical Shock Sensitivity Index for Explosives. J. Phys. Chem. A 2012, 116, 1794−1800. (71) In the original definition of SI provided in ref 70, N was defined as the number of bonds in the molecule. Using the number of atoms clearly makes virtually no difference in practice. However, it is more consistent with the idea that the energy released by decomposing molecules gets spread over all degrees of freedom, whose number is directly proportional to the number of atoms. (72) Kamlet, M. J.; Jacob, S. J. Chemistry of Detonations. I. A Simple Method for Calculating Detonation Properties of C−H−N−O Explosives. J. Chem. Phys. 1968, 48, 23. (73) Rocha, G. B.; Freire, R. O.; Simas, A. M.; Stewart, J. J. P. RM1: a Reparameterization of AM1 for H, C, N, O, P, S, F, Cl, Br, and I. J. Comput. Chem. 2006, 27, 1101−1111. (74) Mopac7. http://sourceforge.net/projects/mopac7 (accessed Feb 8, 2013). (75) NIST Chemistry WebBook. http://webbook.nist.gov/chemistry (accessed Feb 8, 2013). (76) Neese, F. The ORCA Program System. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2012, 2, 73−78. (77) Scott, A. P.; Radom, L. Harmonic Vibrational Frequencies: An Evaluation of Hartree-Fock, Moller-Plesset, Quadratic Configuration Interaction, Density Functional Theory, and Semiempirical Scale Factors. J. Phys. Chem. 1996, 100, 16502−16513. (78) Weigend, F.; Ahlrichs, R. Balanced Basis Sets of Split Valence, Triple Zeta Valence and Quadruple Zeta Valence Quality for H to Rn: Design and Assessment of Accuracy. Phys. Chem. Chem. Phys. 2005, 7, 3297−3305. (79) Neese, F.; Wennmohs, F.; Hansen, A.; Becker, U. Efficient, Approximate and Parallel Hartree−Fock and Hybrid DFT calculations. A Chain of Spheres Algorithm for the Hartree−Fock Exchange. Chem. Phys. 2012, 356, 98−109. (80) Grimme, S. Accurate Description of van der Waals Complexes by Density Functional Theory Including Empirical Corrections. J. Comput. Chem. 2004, 25, 1463−1473. (81) Yoo, C. S.; Holmes, N. C.; Souers, P. C.; Wu, C. J.; Ree, F. H.; Dick, J. J. Anisotropic Shock Sensitivity and Detonation Temperature of Pentaerythritol Tetranitrate Single Crystal. J. Appl. Phys. 2000, 88, 70−75. (82) Zybin, S. V.; Goddard, W. A., III; Xu, P.; van Duin, A. C. T.; Thompson, A. P. Physical Mechanism of Anisotropic Sensitivity in Pentaerythritol Tetranitrate from Compressive-Shear Reaction Dynamics Simulations. Appl. Phys. Lett. 2010, 96, 081918. (83) Zhou, T.-T.; Huang, F.-L. Effects of Defects on Thermal Decomposition of HMX via ReaxFF Molecular Dynamics Simulations. J. Phys. Chem. B 2011, 115, 278−287.

2259

dx.doi.org/10.1021/jp311677s | J. Phys. Chem. A 2013, 117, 2253−2259