Development of a ReaxFF Reactive Force Field for Ammonium Nitrate

Jan 30, 2014 - A new chemical approach for preventing the anti-social use of ... Development, applications and challenges of ReaxFF reactive force fie...
1 downloads 0 Views 3MB Size
Article pubs.acs.org/JPCA

Development of a ReaxFF Reactive Force Field for Ammonium Nitrate and Application to Shock Compression and Thermal Decomposition Tzu-Ray Shan,*,† Adri C. T. van Duin,‡ and Aidan P. Thompson† †

Sandia National Laboratories, Albuquerque, New Mexico 87185, United States Department of Mechanical and Nuclear Engineering, Pennsylvania State University, University Park, Pennsylvania 16802, United States



S Supporting Information *

ABSTRACT: We have developed a new ReaxFF reactive force field parametrization for ammonium nitrate. Starting with an existing nitramine/TATB ReaxFF parametrization, we optimized it to reproduce electronic structure calculations for dissociation barriers, heats of formation, and crystal structure properties of ammonium nitrate phases. We have used it to predict the isothermal pressure−volume curve and the unreacted principal Hugoniot states. The predicted isothermal pressure−volume curve for phase IV solid ammonium nitrate agreed with electronic structure calculations and experimental data within 10% error for the considered range of compression. The predicted unreacted principal Hugoniot states were approximately 17% stiffer than experimental measurements. We then simulated thermal decomposition during heating to 2500 K. Thermal decomposition pathways agreed with experimental findings.

1. INTRODUCTION Ammonium nitrate (AN, NH4NO3) is a strong oxidizer. It is commonly used as a high-nitrogen fertilizer, but it becomes an explosive when mixed with a few percent of fuel oil (hydrocarbons), which is widely used in mining, excavation, and other industrial activities.1 AN is also responsible for many large accidents (most recently the 2013 fertilizer factory explosion in West, Texas) as well as notorious terrorist bombings over the past century.2,3 AN forms monovalent polyatomic NH4+−NO3− ion pairs in crystalline solid phases and neutral acid−base pairs in the gas phase.4 AN crystallizes into five well-characterized phases at atmospheric pressure, ranging from the highest temperature phase I to the lowest temperature phase V.5−12 Crystallographic information and stable temperature ranges of these phases are summarized in Table 1. Phase IV of AN crystal is stable at room temperature and atmospheric pressure. Illustrated in Figure 1 is a ball-andstick model of phase IV of AN crystal from experimental crystallographic data.8 A Hugoniot state is a thermodynamic equilibrium state reached behind a shockwave. Given an initial equilibrium state with defined initial pressure, density, and energy, the corresponding Hugoniot is the locus of new thermodynamic states of increasing pressure, density, and energy satisfying the conditions of mass, momentum, and energy conservation.13,14 The principal Hugoniot is the locus of states associated with an initial, uncompressed state that is at ambient conditions (room temperature and atmospheric pressure). The principal © 2014 American Chemical Society

Hugoniot is often expressed in the form of an approximate linear relationship between shock velocity US and particle velocity UP. Several experimental studies have been reported pertaining to the static and dynamic loading of AN. On the basis of measured bulk sound speed in porous AN samples, Dremin et al.15 estimated the Hugoniot for AN at solid density 1.73 g/cm3 to be US = 2.20 + 1.96 × UP, where US and UP are both in units of km/s. Work by Sandstrom et al.16,17 utilized explosively generated plane shock waves, resulting in the estimate US = 3.07 + 1.89 × UP for 0.28 ≤ UP ≤ 0.94 and US = 2.23 + 2.18 × UP for 0.72 ≤ UP ≤ 2.12. Courchinoux and Lalle18 reported US = 1.80 + 1.80 × UP using the impedance mismatch method. Robbins et al.19 estimated the solid AN Hugoniot to be US = 3.10 + 1.60 × UP with pressed AN crystals in a two stage light gas gun-driven plate impact experiment. In all these experiments15−19 no indication of reaction or initiation were found. The substantial variations in these reported Hugoniot lines can be attributed to differences in density and experimental technique. In addition to experimental efforts in investigating the properties of AN, several modeling and simulation studies have been reported in this field. For example, Sorescu et al.20 have performed plane-wave ab initio calculations to investigate the structural and electronic properties of phases II, III, IV, and V Received: August 21, 2013 Revised: January 24, 2014 Published: January 30, 2014 1469

dx.doi.org/10.1021/jp408397n | J. Phys. Chem. A 2014, 118, 1469−1478

The Journal of Physical Chemistry A

Article

Table 1. Crystal Structures of Ammonium Nitrate Phases under Atmospheric Pressurea crystal space group temperature range (K) a

I5

II6

III7

IV8

V9

cubic Pm3m 398−442

tetragonal P421m 357−398

orthorhombic Pnma 305−357

orthorhombic Pmmn 255−305

orthorhombic Pccn 1.0 km/s.

the mean stress from ReaxFF MD simulations as functions of strain. Although deviations up to 15% can be observed for individual stress components, phase IV of the AN crystal exhibits near hydrostatic response under isotropic compression. A similar response is observed for DFT calculations. Therefore, mean stress is used to represent pressure in Figure 4a. It is interesting to note that under isotropic compression, the b-axis exhibits the lowest stress. Hence, if the supercell is compressed under truly hydrostatic conditions, the b-axis will be undergo the largest strain. This agrees with the findings of Chellappa et al.12 Overall, this indicates that ReaxFF with the AN force field correctly describes the high pressure equation of state of AN phase IV. 3.4. Bulk Unreacted Hugoniot States. To investigate shock properties of AN phase IV with ReaxFF using the AN force field, points along the principal Hugoniot are found by isotropically compressing the system to a chosen volume and then equilibrating in a ReaxFF MD simulation using the canonical ensemble (NVT with Nosé−Hoover thermostat) at

simulations with the AN force field compared to three sets of experimental data.15,16,18,19 Assuming a linear relationship between the shock and particle velocities: US = C0 + s × UP, the intercepts C0 and the slopes s are compared in Table 4. Compared to average experimental values of C0 and s of 2.64 and 1.79, respectively, values for C0 and s from ReaxFF are, respectively, 15% and 17% larger. The ReaxFF Hugoniot states for particle velocities below approximately 1.0 km/s fall within experimental variation, while ReaxFF overpredicts the shock velocity at larger shock strengths. This can be attributed to the ReaxFF model having an uncompressed density (1.81 g/cm3) that is 4% to 11% higher than the densities of the material used in the shock experiments (1.63 or 1.73 g/cm3). 1473

dx.doi.org/10.1021/jp408397n | J. Phys. Chem. A 2014, 118, 1469−1478

The Journal of Physical Chemistry A

Article

NH4NO3 → NH3 + HNO3

Table 4. Intercepts and Slopes of Shock Velocity-Particle Velocity (US−UP) Hugoniot Relationships (US = C0 + s × UP) and Uncompressed Densities of AN Structures exptl, Dremin15 exptl, Courchinoux18 Sandstrom17 exptl, Robbins19 ReaxFF AN

C0 (km/s)

s

density (g/cm3)

2.20 1.80 3.07 2.23 3.10 2.64

1.96 1.80 1.89 2.18 1.60 2.10

1.73 1.69 1.63

NH4NO3 →

(R1)

1 N2 + NO + 2H 2O 2

NH4NO3 → N2O + 2H 2O NH4NO3 →

1.72 1.81

Hugoniot states on the pressure−volume (specific) plane from ReaxFF MD simulations compared to available experimental data16,19 are shown in Figure 6. At pressures below 25

NH4NO3 →

(R3)

3 1 N2 + NO2 + 2H 2O 4 2

NH4NO3 → N2 + 2H 2O +

(R2)

1 O2 2

1 1 1 N2 + NO + NO2 + 2H 2O 8 2 4

(R4)

(R5)

(R6)

Reaction R1 is generally accepted as the thermal decomposition initiating reaction, which is an endothermic proton transfer.78 When the molten salt is heated from 470 to 500 K, exothermic decomposition reactions R2 or R3 occurs, depending on the heating rate.87 Reaction R3 is more common and is the basis for commercial production of nitrous oxide (N2O).80 When heated above 500 K, reaction R4 is the dominant decomposition pathway.2,81 Reaction R5 has been suggested as the decomposition pathway in detonation, while reaction R6 has been associated with subsonic explosions.78 The following detailed mechanism has been proposed for reaction R3:80 NH4NO3 → NH3 + HNO3 → NH3 + NO2 + OH (R3a) Figure 6. Bulk unreacted pressure−volume Hugoniot states of ammonium nitrate from experiments16,19 and ReaxFF-AN simulations. Hugoniot states on the P−V plane from ReaxFF MD lie within the range of experimental variation.

OH + NH3 → NH 2 + H 2O

(R3b)

NH 2 + NO2 → H 2NO + NO → N2O + H 2O

(R3c)

For reaction R4, the following mechanism has been reported:2,81,84 NH4NO3 → NH3 + HNO3

GPa, the ReaxFF Hugoniot points in the pressure−volume plane are consistent with the data from the experiments. At the highest pressure (34 GPa, 0.385 cm3/g), the ReaxFF model once again appears to be too stiff. As mentioned above, this is likely due to the uncompressed density being 4% to 11% larger. Another source of deviation is the different ways the cells are compressed: isotropically for ReaxFF MD simulations versus uniaxially for experiments. Isotropic compression significantly limits shear deformation, hence prohibiting the relief of stress via such deformation. A more thorough investigation of bulk unreacted Hugoniot states of AN from ReaxFF MD simulations via uniaxial compression techniques, such as nonequilibrium direct shock simulations and multiscale shock technique (MSST),77 is reserved for a further study. 3.5. Thermal Decomposition of Ammonium Nitrate Phase IV. The decomposition chemistry of AN has been studied and reported quite extensively due to its application as an explosive component; interested readers are referred to review articles by Oommen and Jain78 and Chaturvedi and Dave, 79 as well as several studies of more specific aspects.2,75,80−87 The decomposition chemistry of AN heavily depends upon environmental factors such as temperature, pressure, moisture content, and the presence of additives, as well as upon experimental conditions like sample size, sample purity, and heating rate. No single mechanism can explain all possible aspects of its decomposition chemistry. Reported decomposition pathways are summarized below.

(R4a)

HNO3 + HX → NO2 + H 2O + X (HX = NH4 +, H3O+ , HNO3)

(R4b)

NH3 + 2NO2 → NH3 + NO+NO3− → HNO3 + NH 2NO

(R4c)

NH 2NO → HN2OH → N2 + H 2O

(R4d)

Thermal decomposition simulations of phase IV of AN crystal are performed with ReaxFF using the AN force field. Perfect 4 × 4 × 4 phase IV supercells are heated gradually from 0 K to target temperatures of 500, 1000, 1500, 2000, and 2500 K, in five different simulations. In each case, the heating time is 40 ps, and the volume of the cell is kept constant. Heated samples are maintained at target temperatures for another 1000 ps. Chemical species analyses are based on time-averaged bond order values over 100 fs that are sampled every 1 fs. It is found that NH4NO3 remains stable for target temperatures of 500 and 1000 K even after 1 ns, suggesting that the first step of decomposition (reaction R1), which is the endothermic dissociation to ammonia and nitric acid, does not occur at these conditions. For the target temperature of 1500 K, we observe the endothermic dissociation of NH4NO3 (reaction R1), but no subsequent reactions occur. Stable chemical species 1474

dx.doi.org/10.1021/jp408397n | J. Phys. Chem. A 2014, 118, 1469−1478

The Journal of Physical Chemistry A

Article

found after 1 ns at 1500 K are coexisting NH4+, NO3−, NH3, and HNO3 at roughly equal concentrations. Subsequent reactions involving the dissociation of nitric acid (HNO3) occur for target temperatures of 2000 and 2500 K. Illustrated in Figure 7 are the evolutions of chemical species

Figure 8. Evolution of stress and temperature during the thermal decomposition NPT MD simulation at 2500 K using the ReaxFF AN force field.

4. CONCLUSIONS We present the development of a reactive force field (ReaxFF) description for ammonium nitrate (NH4NO3, AN). The force field was trained by extending the ReaxFF description for nitramines/TATB71,72 to include AN by training the hydrogen shift energy barrier between NH4+ and NO3− against our electronic structure calculations. Lattice parameters for AN phases I, IV, and V are within ±3.0% error compared to electronic structure calculations. The predicted isothermal pressure−volume curve for phase IV AN crystal agrees with our electronic structure calculations and experimental data10−12 within 10% error. Predicted shock velocity−particle velocity (US−UP) and pressure−volume (P−V) Hugoniot states for unreacted phase IV of AN crystal fall within experimental variation15,16,18,19 below particle velocity of 1.0 km/s and pressure of 15 GPa, respectively. Deviations in the higher compression range can be attributed in part to differences between the experimental and simulation starting structures, as well as the use of isotropic compression in the MD simulations. High-temperature thermal decomposition simulations of phase IV show decomposition of molten AN to NH3 and HNO3 as the first step, with subsequent dissociation of HNO 3 , in agreement with experimental findings.2,81,84

Figure 7. Thermal decomposition of ammonium nitrate phase IV from an NPT MD simulation at 2500 K using the ReaxFF AN force field. Reaction R4a: NH4NO3 → NH3 + HNO3. Reaction R4b: HNO3 + HX → NO2 + H2O + X. Reaction R4c + R4d: NH3 + NO2 → N2 + H2O + OH. The results agree with the experimentally observed hightemperature decomposition pathway.2,81,84

found in the systems at 2500 K. The dissociation of NH4NO3 to NH3 and HNO3 (reaction R4a) initiates at approximately 1200 K. The subsequent reaction is between HNO3 and HNO3/NH4+ yielding NO2, H2O, and NO3−/NH3 (reaction R4b, where HX = HNO3 or NH4+), initiating at approximately 2400 K. This explains a faster decrease in HNO3 concentration than NH3. Figure 7 also illustrates the final reaction between NH3 and NO2 forming products H2O and N2 (combining reactions R4c and R4d), initiated after reaching 2500 K. Although not shown in Figure 7, we do find trace amounts of NO+NO3− and NH2NO in the system. Since these two products are metastable, it is expected that their concentrations remain low. Major products found in the system after 1 ns are H2O, N2, NO2, and OH. The presence of NO2 is likely due to reaction R4c not fully completing. The presence of OH is due to high-temperature homolysis of H2O. Plotted in Figure 8 are the stress and temperature profiles during the 2500 K thermal decomposition simulation. Principal stresses as shown in Figure 8 are 5.5 GPa on average. In the MD simulations, the temperature required for the first dissociation step and subsequent chemical reactions are significantly larger than is seen in experiments (1200 K versus 500 K in experiments). In addition to the larger NH4NO3 decomposition barrier (39 kcal/mol versus 30−41 kcal/mol from experiments78), this is also likely due to the relatively small sample size (2 nm × 2 nm × 2 nm), complete absence of initial defects, and fast heating rate, all of which lower the nucleation rate relative to the experimental conditions, which involve large polycrystalline samples with many defects and grain boundaries and much slower heating rates. Nevertheless, the experimentally observed high-temperature decomposition pathway (reactions R4 and R4a−R4d) is reproduced by the ReaxFF AN force field.



ASSOCIATED CONTENT



AUTHOR INFORMATION

S Supporting Information *

Reactive MD force field data. This material is available free of charge via the Internet at http://pubs.acs.org. Corresponding Author

*(T.-R.S.) E-mail: [email protected]. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS T.-R.S. acknowledges helpful discussions with department manager John B. Aidun. T.-R.S. and A.P.T. acknowledge funding support from Department of Energy’s Advanced Simulation and Computing and Sandia National Laboratories’ 1475

dx.doi.org/10.1021/jp408397n | J. Phys. Chem. A 2014, 118, 1469−1478

The Journal of Physical Chemistry A

Article

(22) Manaa, M.; Reed, E.; Fried, L.; Galli, G.; Gygi, F. Early Chemistry in Hot and Dense Nitromethane: Molecular Dynamics Simulations. J. Chem. Phys. 2004, 120, 10146−10153. (23) Reed, E. J.; Manaa, M. R.; Fried, L. E.; Glaesemann, K. R.; Joannopoulos, J. D. A Transient Semimetallic Layer in Detonating Nitromethane. Nat. Phys. 2008, 4, 72−76. (24) Wu, C. J.; Fried, L. E.; Yang, L. H.; Goldman, N.; Bastea, S. Catalytic Behaviour of Dense Hot Water. Nat. Chem. 2009, 1, 57−62. (25) Cawkwell, M. J.; Sewell, T. D.; Zheng, L.; Thompson, D. L. Shock-Induced Shear Bands in an Energetic Molecular Crystal: Application of Shock-Front Absorbing Boundary Conditions to Molecular Dynamics Simulations. Phys. Rev. B 2008, 78, 014107. (26) Eason, R. M.; Sewell, T. D. Shock-Induced Inelastic Deformation in Orientated Crystalline Pentaerythritol Tetranitrate. J. Phys. Chem. C 2012, 116, 2226−2239. (27) van Duin, A.; Dasgupta, S.; Lorant, F.; Goddard, W. ReaxFF: A Reactive Force Field for Hydrocarbons. J. Phys. Chem. A 2001, 105, 9396−9409. (28) Chenoweth, K.; van Duin, A. C. T.; Goddard, I.; William, A. ReaxFF Reactive Force Field for Molecular Dynamics Simulations of Hydrocarbon Oxidation. J. Phys. Chem. A 2008, 112, 1040−1053. (29) Shin, Y. K.; Shan, T.-R.; Liang, T.; Noordhoek, M. J.; Sinnott, S. B.; van Duin, A. C. T.; Phillpot, S. R. Variable Charge Many-Body Interatomic Potentials. MRS Bull. 2012, 37, 504−512. (30) Budzien, J.; Thompson, A. P.; Zybin, S. V. Reactive Molecular Dynamics Simulations of Shock Through a Single Crystal of Pentaerythritol Tetranitrate. J. Phys. Chem. B 2009, 113, 13142− 13151. (31) Zybin, S. V.; Goddard, I.; William, A.; 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. (32) Mattsson, T. R.; Lane, J. M. D.; Cochrane, K. R.; Desjarlais, M. P.; Thompson, A. P.; Pierce, F.; Grest, G. S. First-Principles and Classical Molecular Dynamics Simulation of Shocked Polymers. Phys. Rev. B 2010, 81, 054103. (33) Dlott, D. New Developments in the Physical Chemistry of Shock Compression. Annu. Rev. Phys. Chem. 2011, 62, 575−597. (34) Yan, Q.-L.; Zeman, S. Theoretical Evaluation of Sensitivity and Thermal Stability for High Explosives Based on Quantum Chemistry Methods: A Brief Review. Int. J. Quantum Chem. 2013, 113, 1049− 1061. (35) Ge, N.-N.; Wei, Y.-K.; Ji, G.-F.; Chen, X.-R.; Zhao, F.; Wei, D.Q. Initial Decomposition of the Condensed-Phase β-HMX under Shock Waves: Molecular Dynamics Simulations. J. Phys. Chem. B 2012, 116, 13696−13704. (36) Long, Y.; Liu, Y. G.; Nie, F. D.; Chen, J. Theoretical Study of Impacting and Desensitizing Theoretical Study of Impacting and Desensitizing for HMX-Graphite Mixture Explosive. Shock Waves 2012, 22, 605−614. (37) Zhou, T.; Zybin, S.; Liu, Y.; Huang, F.; Goddard, I.; William, A. Anisotropic Shock Sensitivity for β-Octahydro-1,3,5,7-tetranitro1,3,5,7-tetrazocine Energetic Material Under Compressive-Shear Loading from ReaxFF- Reactive Dynamics Simulations. J. Appl. Phys. 2012, 111, 124904. (38) An, Q.; Liu, Y.; Zybin, S. V.; Kim, H.; Goddard, I.; William, A. Anisotropic Shock Sensitivity of Cyclotrimethylene Trinitramine (RDX) from Compress-and-Shear Reactive Dynamics. J. Phys. Chem. C 2012, 116, 10198−10206. (39) Zhu, W.; Huang, H.; Huang, H.; Xiao, H. Initial Chemical Events in Shocked Cctahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine: A New Initiation Decomposition Mechanism. J. Chem. Phys. 2012, 136, 044516. (40) Reed, E. J. Electron-Ion Coupling in Shocked Energetic Materials. J. Phys. Chem. C 2012, 116, 2205−2211. (41) An, Q.; Zybin, S.; Goddard, I.; William, A.; Jaramillo-Botero, A.; Blanco, M.; Luo, S.-N. Elucidation of the Dynamics for Hot-Spot Initiation at Nonuniform Interfaces of Highly Shocked Materials. Phys. Rev. B 2011, 84, 220101.

Laboratory Directed Research and Development. Sandia National Laboratories is a multiprogram laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-AC04-94AL85000. SAND NO. 2013-1054J



REFERENCES

(1) Cook, M. The Science of Industrial Explosives. IRECO Chemicals: Salt Lake City, UT, 1974 (2) Oxley, J.; Smith, J.; Rogers, E.; Yu, M. Ammonium Nitrate: Thermal Stability and Explosivity Modifiers. Thermochim. Acta 2002, 384, 23−45. (3) Ammonium Nitrate Disasters, 2013. http://en.wikipedia.org/ wiki/Ammonium_nitrate_disasters. (4) Alavi, S.; Thompson, D. L. Theoretical Study of Proton Transfer in Ammonium Nitrate Clusters. J. Chem. Phys. 2002, 117, 2599−2608. (5) Hermann, M.; Engel, W. Phase Transitions and Lattice Dynamics of Ammonium Nitrate. Propellants, Explos., Pyrotech. 1997, 22, 143− 147. (6) Choi, C.; Mapes, J.; Prince, E. The Structure of Ammonium Nitrate (IV). Acta Crystallogr. 1972, B28, 1357−1361. (7) Ahtee, M.; Smolander, K.; Lucas, B.; Hewat, A. The Structure of the Low-Temperature Phase V of Ammonium Nitrate, ND4NO3. Acta Crystallogr. 1983, C39, 651−655. (8) Lucas, B.; Ahtee, M.; Hewat, A. Crystal-Structure of Phase-II Ammonium Nitrate. Acta Crystallogr. 1979, B35, 1038−1041. (9) Ahtee, M.; Kurki-Suonio, K.; Lucas, B.; Hewat, A. Determination of Molecular Orientations in Cubic ND4NO3 by Multipole Analysis. Acta Crystallogr. 1979, A35, 591−597. (10) Davidson, A. J.; Chellappa, R. S.; Dattelbaum, D. M.; Yoo, C.-S. Pressure Induced Isostructural Metastable Phase Transition of Ammonium Nitrate. J. Phys. Chem. A 2011, 115, 11889−11896. (11) Dunuwille, M.; Davidson, A. J.; Chellappa, R. S.; Dattelbaum, D. M.; Yoo, C.-S. Pressure-Induced Isostructural Metastable Phase Transition of Ammonium Nitrate. J. Phys. Chem. A 2012, 116, 7600−7601. (12) Chellappa, R. S.; Dattelbaum, D. M.; Velisavljevic, N.; Sheffield, S. The Phase Diagram of Ammonium Nitrate. J. Chem. Phys. 2012, 137, 064504. (13) Rankine, W. J. M. On the Thermodynamic Theory of Waves of Finite Longitudinal Disturbance. Phil. Trans. R. Soc. Lond. 1870, 160, 277−288. (14) Hugoniot, P. H. Memoire Sur la Propagation du Mouvement Dans les Corps et Plus Spécialement Dans les Gaz Parfaits. J. Ec. Polytech. 1887, 57, 3−97. (15) Dremin, A.; Shvedov, K.; Avdonin, O. Shock Compressibility and Temperature of Certain Explosives in the Porous State. Combust. Explos. Shock Waves 1970, 6, 449−455. (16) Sandstrom, F. W.; Persson, P. A.; Olinger, B. Isothermal and Shock Compression of High Density Ammonium Nitrate and Ammonium Perchlorate. Proc. 10th Inter. De. Symp., Off. Of Naval Res. Report ONR 33395−12 1993, 766−774. (17) Sandstrom, F. W.; Persson, P. A.; Olinger, B. Isothermal and Shock Compression of High Density Ammonium Nitrate and Ammonium Perchlorate. AIP Conf. Proc. 1994, 309, 1409−1412. (18) Courchinoux, R.; Lalle, P. Unreacted Hugoniot of Ammonium Nitrate. AIP Conf. Proc. 1994, 309, 1397−1400. (19) Robbins, D. L.; Sheffield, S. A.; Dattelbaum, D. M.; Velisavljevic, N.; Stahl, D. B. Equation of State of Ammonium Nitrate. AIP Conf. Proc. 2009, 1195, 552−555. (20) Sorescu, D.; Thompson, D. L. Classical and Quantum Mechanical Studies of Crystalline Ammonium Nitrate. J. Phys. Chem. A 2001, 105, 720−733. (21) Velardez, G.; Alavi, S.; Thompson, D. L. Molecular Dynamics Studies of Melting and Solid-State Transitions of Ammonium Nitrate. J. Chem. Phys. 2004, 120, 9151−9159. 1476

dx.doi.org/10.1021/jp408397n | J. Phys. Chem. A 2014, 118, 1469−1478

The Journal of Physical Chemistry A

Article

(42) Chen, X.-F.; Liu, J.-F.; Meng, Z.-H.; Han, K.-L. Thermal Unimolecular Decomposition Mechanism of 2,4,6-Trinitrotoluene: A First-Principles DFT Study. Theor. Chem. Acc. 2010, 127, 327−344. (43) Ojeda, O. U.; Cagin, T. Multiscale Modeling of Crystalline Energetic Materials. Comput. Mater. Continua 2010, 16, 127−173. (44) Cawkwell, M. J.; Ramos, K. J.; Hooks, D. E.; Sewell, T. D. Homogeneous Dislocation Nucleation in Cyclotrimethylene Trinitramine Under Shock Loading. J. Appl. Phys. 2010, 107, 063512. (45) Cui, H.-L.; Ji, G.-F.; Chen, X.-R.; Zhu, W.-H.; Zhao, F.; Wen, Y.; Wei, D.-Q. First-Principles Study of High-Pressure Behavior of Solid Beta-HMX. J. Phys. Chem. A 2010, 114, 1082−1092. (46) Landerville, A. C.; Oleynik, I. I.; White, C. T. Reactive Molecular Dynamics of Hypervelocity Collisions of PETN Molecules. J. Phys. Chem. A 2009, 113, 12094−12104. (47) Bedrov, D.; Hooper, J. B.; Smith, G. D.; Sewell, T. D. ShockInduced Transformations in Crystalline RDX: A Uniaxial ConstantStress Hugoniostat Molecular Dynamics Simulation Study. J. Chem. Phys. 2009, 131, 034712. (48) Zhu, W.; Xiao, J.-J.; Zheng, J.; Zhao, X.-B.; Chen, Z.-E.; Xiao, H.-M. A Theoretical Criterion for Sensitivity of Energetic CompositesMolecular Dynamics Studies on AP/HMX Systems at Various Concentrations and Temperatures. Acta Chim. Sin. 2008, 66, 2592− 2596. (49) Conroy, M. W.; Oleynik, I. I.; Zybin, S. V.; White, C. T. FirstPrinciples Anisotropic Constitutive Relationships in Beta-Cyclotetramethylene Tetranitramine (beta-HMX). J. Appl. Phys. 2008, 104, 053506. (50) Cawkwell, M. J.; Sewell, T. D.; Zheng, L.; Thompson, D. L. Shock-Induced Shear Bands in an Energetic Molecular Crystal: Application of Shock-Front Absorbing Boundary Conditions to Molecular Dynamics Simulations. Phys. Rev. B 2008, 78, 014107. (51) Conroy, M. W.; Oleynik, I. I.; Zybin, S. V.; White, C. T. FirstPrinciples Investigation of Anisotropic Constitutive Relationships in Pentaerythritol Tetranitrate. Phys. Rev. B 2008, 77, 094107. (52) Nomura, K.-I.; Kalia, R. K.; Nakano, A.; Vashishta, P. Reactive Nanojets: Nanostructure-Enhanced Chemical Reactions in a Defected Dnergetic Crystal. Appl. Phys. Lett. 2007, 91, 183109. (53) Vashishta, P.; Kalia, R. K.; Nakano, A.; Homan, B. E.; McNesby, K. L. Multimillion Atom Reactive Simulations of Nanostructured Energetic Materials. J. Propul. Power 2007, 23, 688−692. (54) Plimpton, S. Fast Parallel Algorithms for Short Range Molecular Dynamics. J. Comput. Phys 1995, 117, 1−19. (55) Aktulga, H. M.; Fogarty, J. C.; Pandit, S. A.; Grama, A. Y. Parallel Reactive Molecular Dynamics: Numerical Methods and Algorithmic Techniques. Parallel Comput. 2012, 38, 245−259. (56) Dunning, T. Gaussian Basis Sets for Use in Correlated Molecular Calculations. I. The Atoms Boron through Neon and Hydrogen. J. Chem. Phys. 1989, 90, 1007−1023. (57) Becke, A. Density-Functional Exchange-Energy Approximation with Correct Asymptotic Behavior. Phys. Rev. A 1988, 38, 3098−3100. (58) Lee, C.; Yang, W.; Parr, R. Development of the Colle−Salvetti Correlation-Energy Formula into a Functional of the Electron-Density. Phys. Rev. B 1988, 37, 785−789. (59) Valiev, M.; et al. NWChem: A Comprehensive and Scalable Open-Source Solution for Large Scale Molecular Simulations. Comput. Phys. Commun. 2010, 181, 1477−1489. (60) Kresse, G.; Hafner, J. Ab Initio Molecular Dynamics for Liquid Metals. Phys. Rev. B 1993, 47, 558−561. (61) Kresse, G.; Hafner, J. Ab Initio Molecular Dynamics Simulation of the Liquid-Metal-Amorphous-Semiconductor Transition in Germanium. Phys. Rev. B 1994, 49, 14251−14269. (62) Kresse, G.; Furthmüller, J. Efficient Iterative Schemes for ab Initio Total-Energy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B 1996, 54, 11169−11186. (63) Blochl, P. Projector Augmented-Wave Method. Phys. Rev. B 1994, 50, 17953−17979. (64) Kresse, G.; Joubert, D. From Ultrasoft Pseudopotentials to the Projector Augmented-Wave Method. Phys. Rev. B 1999, 59, 1758− 1775.

(65) Perdew, J.; Chevary, J.; Vosko, S.; Jackson, K.; Pederson, M.; Singh, D.; Fiolhais, C. Atoms, Molecules, Solids, and Surfaces: Applications of the Generalized Gradient Approximation for Exchange and Correlation. Phys. Rev. B 1992, 46, 6671−6687. (66) Perdew, J.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. B 1996, 77, 3865−3868. (67) Monkhorst, H.; Pack, J. Special Points for Brillouin-Zone Integrations. Phys. Rev. B 1976, 13, 5188−5192. (68) Strachan, A.; van Duin, A.; Chakraborty, D.; Dasgupta, S.; Goddard, W. Shock Waves in High-Energy Materials: The Initial Chemical Events in Nitramine RDX. Phys. Rev. Lett. 2003, 91, 098301. (69) van Duin, A.; Zeiri, Y.; Dubnikova, F.; Kosloff, R.; Goddard, W. Atomistic-Scale Simulations of the Initial Chemical Events in the Thermal Initiation of Triacetonetriperoxide. J. Am. Chem. Soc. 2005, 127, 11053−11062. (70) van Duin, A. C. T.; Zybin, S. V.; Chenoweth, K.; Zhang, L.; Han, S.-P.; Strachan, A.; Goddard, I.; William, A. Reactive Force Fields Based on Quantum Mechanics for Applications to Materials at Extreme Conditions. AIP Conf. Proc. 2006, 845, 581−584. (71) Zhang, L.; Zybin, S. V.; van Duin, A. C. T.; Dasgupta, S.; Goddard, I.; William, A.; Kober, E. M. Carbon Cluster Formation during Thermal Decomposition of Octahydro-1,3,5,7-tetranitro1,3,5,7-tetrazocine and 1,3,5-Triamino-2,4,6-trinitrobenzene High Explosives from ReaxFF Reactive Molecular Dynamics Simulations. J. Phys. Chem. A 2009, 113, 10619−10640. (72) Zhang, L.; van Duin, A. C. T.; Zybin, S. V.; Goddard, I.; William, A. Thermal Decomposition of Hydrazines from Reactive Dynamics Using the ReaxFF Reactive Force Field. J. Phys. Chem. B 2009, 113, 10770−10778. (73) Rahaman, O.; van Duin, A. C. T.; Goddard, I.; William, A.; Doren, D. J. Development of a ReaxFF Reactive Force Field for Glycine and Application to Solvent Effect and Tautomerization. J. Phys. Chem. B 2011, 115, 249−261. (74) Liu, L.; Liu, Y.; Zybin, S. V.; Sun, H.; Goddard, I.; William, A. ReaxFF-lg: Correction of the ReaxFF Reactive Force Field for London Dispersion, with Applications to the Equations of State for Energetic Materials. J. Phys. Chem. A 2011, 115, 11016−11022. (75) Cagnina, S.; Rotureau, P.; Fayetb, G.; Adamo, C. The Ammonium Nitrate and Its Mechanism of Decomposition in the Gas Phase: A Theoretical Study and a DFT Benchmark. Phys. Chem. Chem. Phys. 2013, 15, 10849−10858. (76) Shan, T.-R.; Wixom, R. R.; Mattsson, A. E.; Thompson, A. P. Atomistic Simulation of Orientation Dependence in Shock-Induced Initiation of Pentaerythritol Tetranitrate. J. Phys. Chem. B 2013, 117, 928−936. (77) Reed, E.; Fried, L.; Joannopoulos, J. A Method for Tractable Dynamical Studies of Single and Double Shock Compression. Phys. Rev. Lett. 2003, 90, 235503. (78) Oommen, C.; Jain, S. Ammonium Nitrate: A Promising Rocket Propellant Oxidizer. J. Hazard. Mater. 1999, 67, 253−281. (79) Chaturvedi, S.; Dave, P. N. Review on Thermal Decomposition of Ammonium Nitrate. J. Energ. Mater. 2013, 31, 1−26. (80) Park, J.; Lin, M. C. Thermal Decomposition of Gaseous Ammonium Nitrate at Low Pressure: Kinetic Modeling of Product Formation and Heterogeneous Decomposition of Nitric Acid. J. Phys. Chem. A 2009, 113, 13556−13561. (81) Sinditskii, V.; Egorshev, V.; Levshenkov, A.; Serushkin, V. Ammonium Nitrate: Combustion Mechanism and the Role of Additives. Propellants, Explos., Pyrotech. 2005, 30, 269−280. (82) Vyazovkin, S.; Clawson, J.; Wight, C. Thermal Dissociation Kinetics of Solid and Liquid Ammonium Nitrate. Chem. Mater. 2001, 13, 960−966. (83) Rossi, M.; Bottaro, J.; McMillen, D. The Thermal Decomposition of the New Energetic Material Ammonium Dinitramide (NH4 N(NO2 )2 ) in Relaton to Nitramide (NH2NO2) AND NH4NO3. Int. J. Chem. Kinet. 1993, 25, 549−570. (84) Brill, T.; Brush, P.; Patil, D. Chemistry of Ammonium Nitrate and Ammonium Dinitramide Near the Burning Surface Temperature. Combust. Flame 1993, 92, 178−186. 1477

dx.doi.org/10.1021/jp408397n | J. Phys. Chem. A 2014, 118, 1469−1478

The Journal of Physical Chemistry A

Article

(85) Oxley, J.; Kaushik, S.; Gilson, N. Thermal Decomposition of Ammonium Nitrate-Based Composites. Thermochim. Acta 1989, 153, 269−286. (86) Brower, K. R.; Oxley, J. C.; Tewari, M. Evidence for Homolytic Decomposition of Ammonium Nitrate at High Tempearture. J. Phys. Chem. 1989, 93, 4029−4033. (87) Andersen, W.; Bills, K.; Mishuck, E.; Moe, G.; Schultz, R. A Model Describing Combustion of Solid Composite Propellants Containing Ammonium Nitrate. Combust. Flame 1959, 3, 301−317.

1478

dx.doi.org/10.1021/jp408397n | J. Phys. Chem. A 2014, 118, 1469−1478