Sergey V. Zybin

ReaxFF-lg: Correction of the ReaxFF Reactive Force Field for London Dispersion, with Applications to the Equations of State for Energetic Materials

The practical levels of density functional theory (DFT) for solids (LDA, PBE, PW91, B3LYP) are well-known not to account adequately for the London dispersion (van der Waals attraction) so important in molecular solids, leading to equilibrium volumes for molecular crystals ~10-15% too high. The ReaxFF reactive force field is based on fitting such DFT calculations and suffers from the same problem. In the paper we extend ReaxFF by adding a London dispersion term with a form such that it has low gradients (lg) at valence distances leaving the already optimized valence interactions intact but behaves as 1/R(6) for large distances. We derive here these lg corrections to ReaxFF based on the experimental crystal structure data for graphite, polyethylene (PE), carbon dioxide, and nitrogen and for energetic materials: hexahydro-1,3,5-trinitro-1,3,5-s-triazine (RDX), pentaerythritol tetranitrate (PETN), 1,3,5-triamino-2,4,6-trinitrobenzene (TATB), and nitromethane (NM). After this dispersion correction the average error of predicted equilibrium volumes decreases from 18.5 to 4.2% for the above systems. We find that the calculated crystal structures and equation of state with ReaxFF-lg are in good agreement with experimental results. In particular, we examined the phase transition between α-RDX and γ-RDX, finding that ReaxFF-lg leads to excellent agreement for both the pressure and volume of this transition occurring at ~4.8 GPa and ~2.18 g/cm(3) density from ReaxFF-lg vs 3.9 GPa and ~2.21 g/cm(3) from experiment. We expect ReaxFF-lg to improve the descriptions of the phase diagrams for other energetic materials.

Carbon Cluster Formation during Thermal Decomposition of Octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine and 1,3,5-Triamino-2,4,6-trinitrobenzene High Explosives from ReaxFF Reactive Molecular Dynamics Simulations

We report molecular dynamics (MD) simulations using the first-principles-based ReaxFF reactive force field to study the thermal decomposition of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) and octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) at various densities and temperatures. TATB is known to produce a large amount (15-30%) of high-molecular-weight carbon clusters, whereas detonation of nitramines such as HMX and RDX (1,3,5-trinitroperhydro-1,3,5-triazine) generate predominantly low-molecular-weight products. In agreement with experimental observation, these simulations predict that TATB decomposition quickly (by 30 ps) initiates the formation of large carbonaceous clusters (more than 4000 amu, or approximately 15-30% of the total system mass), and HMX decomposition leads almost exclusively to small-molecule products. We find that HMX decomposes readily on this time scale at lower temperatures, for which the decomposition rate of TATB is about an order of magnitude slower. Analyzing the ReaxFF MD results leads to the detailed atomistic structure of this carbon-rich phase of TATB and allows characterization of the kinetics and chemistry related to this phase and their dependence on system density and temperature. The carbon-rich phase formed from TATB contains mainly polyaromatic rings with large oxygen content, leading to graphitic regions. We use these results to describe the initial reaction steps of thermal decomposition of HMX and TATB in terms of the rates for forming primary and secondary products, allowing comparison to experimentally derived models. These studies show that MD using the ReaxFF reactive force field provides detailed atomistic information that explains such macroscopic observations as the dramatic difference in carbon cluster formation between TATB and HMX. This shows that ReaxFF MD captures the fundamental differences in the mechanisms of such systems and illustrates how the ReaxFF may be applied to model complex chemical phenomena in energetic materials. The studies here illustrate this for modestly sized systems and modest periods; however, ReaxFF calculations of reactive processes have already been reported on systems with approximately 10(6) atoms. Thus, with suitable computational facilities, one can study the atomistic level chemical processes in complex systems under extreme conditions.

Physical mechanism of anisotropic sensitivity in pentaerythritol tetranitrate from compressive-shear reaction dynamics simulations

We propose computational protocol (compressive shear reactive dynamics) utilizing the ReaxFF reactive force field to study chemical initiation under combined shear and compressive load. We apply it to predict the anisotropic initiation sensitivity observed experimentally for shocked pentaerythritol tetranitrate single crystals. For crystal directions known to be sensitive we find large stress overshoots and fast temperature increase that result in early bond-breaking processes whereas insensitive directions exhibit small stress overshoot, lower temperature increase, and little bond dissociation. These simulations confirm the model of steric hindrance to shear and capture the thermochemical processes dominating the phenomena of shear-induced chemical initiation.

Reactive Molecular Dynamics Simulations of Shock Through a Single Crystal of Pentaerythritol Tetranitrate

Large-scale molecular dynamics simulations and the reactive force field ReaxFF were used to study shock-induced initiation in crystalline pentaerythritol tetranitrate (PETN). In the calculations, a PETN single crystal was impacted against a wall, driving a shockwave back through the crystal in the [100] direction. Two impact speeds (4 and 3 km/s) were used to compare strong and moderate shock behavior. The primary difference between the two shock strengths is the time required to exhibit the same qualitative behaviors with the lower impact speed lagging behind the faster impact speed. For both systems, the shock velocity exhibits an initial deceleration due to onset of endothermic reactions followed by acceleration due to the onset of exothermic reactions. At long times, the shock velocity reaches a steady value. After the initial deceleration period, peaks are observed in the profiles of the density and axial stress with the strongly shocked system having sharp peaks while the weakly shocked system developed broad peaks due to the slower shock velocity acceleration. The dominant initiation reactions in both systems lead to the formation of NO(2) with lesser quantities of NO(3) and formaldehyde also produced.

Decomposition of Condensed Phase Energetic Materials: Interplay between Uni- and Bimolecular Mechanisms

Activation energy for the decomposition of explosives is a crucial parameter of performance. The dramatic suppression of activation energy in condensed phase decomposition of nitroaromatic explosives has been an unresolved issue for over a decade. We rationalize the reduction in activation energy as a result of a mechanistic change from unimolecular decomposition in the gas phase to a series of radical bimolecular reactions in the condensed phase. This is in contrast to other classes of explosives, such as nitramines and nitrate esters, whose decomposition proceeds via unimolecular reactions both in the gas and in the condensed phase. The thermal decomposition of a model nitroaromatic explosive, 2,4,6-trinitrotoluene (TNT), is presented as a prime example. Electronic structure and reactive molecular dynamics (ReaxFF-lg) calculations enable to directly probe the condensed phase chemistry under extreme conditions of temperature and pressure, identifying the key bimolecular radical reactions responsible for the low activation route. This study elucidates the origin of the difference between the activation energies in the gas phase (~62 kcal/mol) and the condensed phase (~35 kcal/mol) of TNT and identifies the corresponding universal principle. On the basis of these findings, the different reactivities of nitro-based organic explosives are rationalized as an interplay between uni- and bimolecular processes.

Anisotropic shock sensitivity for β-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine energetic material under compressive-shear loading from ReaxFF-ℓg reactive dynamics simulations

We report here the predictions on anisotropy of shock sensitivity and of chemical process initiation in single crystal β-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (β-HMX) using compressive shear reactive dynamics (CS-RD) model with ReaxFF-lg reactive force field. Analysis of resolved shear stress induced by uniaxial compression along three shock directions normal to (110), (011), and (010) planes leads to identify eight slip systems as candidates for shear deformation. For each of the eight slip systems, non-equilibrium reactive dynamics simulations were carried out to determine thermal, mechanical, and chemical responses to shear deformation. Shock direction normal to (010) plane exhibits large shear stress barriers arising from steric hindrance between molecules of adjacent layers leading to local dramatic energy and temperature increases under shear flow that in turn accelerate chemical bond breaking and initial product formation processes, promoting further molecular decomposition and eventually transition to detonation. This suggests that single crystal β-HMX is sensitive to shocks in direction normal to (010) plane. Shock directions normal to (110) and (011) planes reveal significantly less steric hindrance, leading to more modest energy and temperature increases followed by slower chemical reaction initiation. Thus, shock directions normal to (110) and (011) planes are less sensitive than shock direction normal to (010) plane, which agree with interpretations from currently available plate impact experiments on HMX. This validation of CS-RD and ReaxFF for characterizing sensitivity of single crystal energetic materials indicates that these methods can be applied to study sensitivity for more complex polymer bonded explosives and solid composite propellants having complex microstructures, corrugated interfaces, as well as defects.

Density-Dependent Liquid Nitromethane Decomposition: Molecular Dynamics Simulations Based on ReaxFF

The decomposition mechanism of hot liquid nitromethane at various compressions was studied using reactive force field (ReaxFF) molecular dynamics simulations. A competition between two different initial thermal decomposition schemes is observed, depending on compression. At low densities, unimolecular C-N bond cleavage is the dominant route, producing CH(3) and NO(2) fragments. As density and pressure rise approaching the Chapman-Jouget detonation conditions (∼30% compression, >2500 K) the dominant mechanism switches to the formation of the CH(3)NO fragment via H-transfer and/or N-O bond rupture. The change in the decomposition mechanism of hot liquid NM leads to a different kinetic and energetic behavior, as well as products distribution. The calculated density dependence of the enthalpy change correlates with the change in initial decomposition reaction mechanism. It can be used as a convenient and useful global parameter for the detection of reaction dynamics. Atomic averaged local diffusion coefficients are shown to be sensitive to the reactions dynamics, and can be used to distinguish between time periods where chemical reactions occur and diffusion-dominated, nonreactive time periods.

Explanation of the Colossal Detonation Sensitivity of Silicon Pentaerythritol Tetranitrate (Si-PETN) Explosive

DFT calculations have identified the novel rearrangement shown here for decomposition of the Si derivative of the PETN explosive [PentaErythritol TetraNitrate (PETN), C(CH(2)ONO(2))(4)] that explains the very dramatic increase in sensitivity observed experimentally. The critical difference is that Si-PETN allows a favorable five-coordinate transition state in which the new Si-O and C-O bonds form simultaneously, leading to a transition state barrier of 33 kcal/mol (it is 80 kcal/mol for PETN) and much lower than the normal O-NO(2) bond fission observed in other energetic materials (approximately 40 kcal/mol). In addition this new mechanism is very exothermic (45 kcal/mol) leading to a large net energy release at the very early stages of Si-PETN decomposition that facilitates a rapid temperature increase and expansion of the reaction zone.

Thermal Decomposition of Hydrazines from Reactive Dynamics Using the ReaxFF Reactive Force Field

We report reactive dynamics (RD) studies on: the decomposition of bulk hydrazine (N(2)H(4)); the decomposition of bulk monomethyl-hydrazine (CH(3)N(2)H(3)), hereafter referred to simply as methyl-hydrazine; the decomposition of hydrazine in the presence of hydrogen peroxide (H(2)O(2)); and decomposition hydrazine on catalytic surfaces Pt[100] and Pt[111] under various conditions. These studies use the ReaxFF reactive force field to describe the multitude of chemical reactions in these systems for a variety of reaction conditions in order to show that this approach leads to realistic decomposition mechanisms and rates. In particular, we determined how the decomposition of hydrazine is affected by temperature, pressure, and heating rate. We analyzed chemical reaction mechanism of the decomposition of hydrazine at the studied conditions and found that at lower temperatures the initial product from hydrazine decomposition is NH(3), whereas at higher temperatures H(2) and N(2) are the dominant early products. Prominent intermediates observed during these decompositions include N(2)H(3), N(2)H(2,) and NH(2), in agreement with quantum mechanical studies (7.3 ps at 3000 K). As the heating rate is decreased, the onset for hydrazine decomposition shifts to lower temperatures. Using a constant heating rate, we found that higher pressure (increased density) favors formation of NH(3) over N(2) and H(2). In studies of the catalytic decomposition of hydrazine on surfaces Pt[100] and Pt[111], we found that the presence of a Pt-catalyst reduces the initial decomposition temperature of hydrazine by about 50%. We found that the Pt[100]-surface is 20 times more active for hydrazine decomposition than the Pt[111]-surface, in qualitative agreement with experiments. These studies indicate how ReaxFF RD can be useful in understanding the chemical processes involved in bulk and catalytic decomposition and in oxidation of reactive species under various reaction conditions.

Density functional theory calculations of solid nitromethane under hydrostatic and uniaxial compressions with empirical van der Waals correction.

First-principles density functional theory calculations have been performed with and without an empirical van der Waals (vdW) correction to obtain constitutive relationships of solid nitromethane under hydrostatic and uniaxial compressions. The unit-cell parameters at zero pressure and the hydrostatic equation of state at 0 K are in reasonable agreement with experimental data using pure DFT, and the agreement is significantly improved with the inclusion of the vdW dispersion correction. Uniaxial compressions normal to the {100}, {010}, {001}, {110}, {101}, {011}, and {111} planes were performed, and a comparison of the principal stresses, changes in energy, and shear stresses for different compression directions clearly indicate anisotropic behavior of solid nitromethane upon compression. The calculated anisotropic constitutive relationships might help to link the anisotropic shock sensitivity and the underlying atomic-scale properties of solid nitromethane.