Charge Distributions of Nitro Groups Within Organic Explosive Crystals: Effects on Sensitivity and Modeling

Theoretical Investigation into Polymorphic Transformation between β-HMX and δ-HMX by Finite Temperature String

Polymorphic transformation is important in chemical industries, in particular, in those involving explosive molecular crystals. However, due to simulating challenges in the rare event method and collective variables, understanding the transformation mechanism of molecular crystals with a complex structure at the molecular level is poor. In this work, with the constructed order parameters (OPs) and K-means clustering algorithm, the potential of mean force (PMF) along the minimum free-energy path connecting β-HMX and δ-HMX was calculated by the finite temperature string method in the collective variables (SMCV), the free-energy profile and nucleation kinetics were obtained by Markovian milestoning with Voronoi tessellations, and the temperature effect on nucleation was also clarified. The barriers of transformation were affected by the finite-size effects. The configuration with the lower potential barrier in the PMF corresponded to the critical nucleus. The time and free-energy barrier of the polymorphic transformation were reduced as the temperature increased, which was explained by the pre-exponential factor and nucleation rate. Thus, the polymorphic transformation of HMX could be controlled by the temperatures, as is consistent with previous experimental results. Finally, the HMX polymorph dependency of the impact sensitivity was discussed. This work provides an effective way to reveal the polymorphic transformation of the molecular crystal with a cyclic molecular structure, and further to prepare the desired explosive by controlling the transformation temperature.

Study on the regulation of ε-CL-20 by an external electric field

CONTEXT

To explore the influence of the external electric field (EEF) on ε-CL-20. The molecular structure, frontier molecular orbitals (FMOs), global reactivity parameters (GRP), surface electrostatic potential, nitro charge, and UV-Vis spectra of ε-CL-20 under EEF were studied using density functional theory (DFT). The calculation results show that the electric field applied along N16-N24 has a significant effect on the structure of ε-CL-20. With an increase in the positive EEF, the bond length of the initiating bond decreases, and the bond order and bond dissociation energy increase, which increases the thermal stability of ε-CL-20 to a certain extent. In addition, with an increase in the positive EEF intensity, the LUMO migrates from both sides of the positive electric field to one side of the nitro group, and the HOMO migrates from the skeleton to the nitro group. It is worth noting that in the negative EEF, when the electric field strength changed from 0 to 0.016 a.u., the negative charge of the total nitro group gradually decreased. When the electric field strength becomes 0.02 a.u., the negative charge of the total nitro group suddenly increases, and ε-CL-20 is significantly polarized. When the electric field strength is sufficiently strong, the occupied and unoccupied orbitals of the ε-CL-20 molecule change, resulting in a change in the energy level difference between the occupied and unoccupied orbitals, which further excites the corresponding excited state, resulting in a new UV-Vis absorption peak.

METHODS

Based on the density functional theory (DFT), the structural optimization and energy calculation were carried out by using B3LYP/6-311 + G(d, p) and B3LYP/def2-TZVPP methods, respectively. After optimization convergence, vibration analysis was performed without imaginary frequencies to obtain stable configurations. Then the molecular structure, frontier molecular orbitals (FMOs), global reactivity parameters (GRP), surface electrostatic potential, nitro charge, and UV-Vis spectra were analyzed.

High-pressure structural studies and pressure-induced sensitisation of 3,4,5-trinitro-1H-pyrazole

Herein we report the first high-pressure study of the energetic material 3,4,5-trinitro-1H-pyrazole (3,4,5-TNP) using neutron powder diffraction and single-crystal X-ray diffraction. A new high-pressure phase, termed Form II, was first identified through a substantial change in the neutron powder diffraction patterns recorded over the range 4.6-5.3 GPa, and was characterised further by compression of a single crystal to 5.3 GPa in a diamond-anvil cell using X-ray diffraction. 3,4,5-TNP was found to be sensitive to initiation under pressure, as demonstrated by its unexpected and violent decomposition at elevated pressures in successive powder diffraction experiments. Initiation coincided with the sluggish phase transition from Form I to Form II. Using a vibrational up-pumping model, its increased sensitivity under pressure can be explained by pressure-induced mode hardening. These findings have potential implications for the safe handling of 3,4,5-TNP, on the basis that shock- or pressure-loading may lead to significantly increased sensitivity to initiation.

Unimolecular Decomposition Reactions of Picric Acid and Its Methylated Derivatives─A DFT Study

To handle energetic materials safely, it is important to have knowledge about their sensitivity. Density functional theory (DFT) has proven a valuable tool in the study of energetic materials, and in the current work, DFT is employed to study the thermal unimolecular decomposition of 2,4,6-trinitrophenol (picric acid, PA), 3-methyl-2,4,6-trinitrophenol (methyl picric acid, mPA), and 3,5-dimethyl-2,4,6-trinitrophenol (dimethyl picric acid, dmPA). These compounds have similar molecular structures, but according to the literature, mPA is far less sensitive to impact than the other two compounds. Three pathways believed important for the initiation reactions are investigated at 0 and 298.15 K. We compare the computed energetics of the reaction pathways with the objective of rationalizing the unexpected sensitivity behavior. Our results reveal a few if any significant differences in the energetics of the three molecules, and thus do not reflect the sensitivity deviations observed in experiments. These findings point toward the potential importance of crystal structure, crystal morphology, bimolecular reactions, or combinations thereof on the impact sensitivity of nitroaromatics.

Analysis of Ignition Sites for the Explosives 3,3'-Diamino-4,4'-azoxyfurazan (DAAF) and 1,3,5,7-Tetranitro-1,3,5,7-tetrazoctane (HMX) Using Crush Gun Impact Testing

The handling safety characteristics of energetic materials must be measured in order to ensure the safe transport and use of explosives. Drop-weight impact sensitivity measurements are one of the first standardized tests performed for energetics. They utilize a small amount of the explosive sample and a standard weight, which is dropped on the material from various heights to determine its sensitivity. While multiple laboratories have used the impact sensitivity test as an initial screening tool for explosive sensitivity for the past 60 years, variability exists due to the use of different instruments, different methods to determine the initiation, and the scatter commonly associated with less-sensitive explosives. For example, standard explosives such as 1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (HMX) initiate reliably and consistently on the drop-weight impact test, whereas insensitive explosives such as 3,3'-diamino-4,4'-azoxyfurazan (DAAF) exhibit variability in sound levels and the expended material. Herein we investigate the impact sensitivity of DAAF and HMX along with a more detailed investigation of ignition sites using a novel "crush gun" apparatus: a pneumatically powered drop-weight tower with advanced diagnostics, including high-speed visual and infrared cameras. Using this crush gun assembly, the ignition sites in HMX and DAAF were analyzed with respect to the effects of particle size and the presence of a source of grit. The formation of ignition sites was observed in both explosives; however, only HMX showed ignition sites that propagated to a deflagration at lower firing speeds. Finally, the presence of grit particles was shown to increase the occurrence of ignition sites in DAAF at lower firing speeds, though propagation to a full reaction was not observed on the time scale of the test. These results enable a better understanding of how ignition and propagation occurs during the impact testing of DAAF.

A non-empirical intermolecular force-field for trinitrobenzene and its application in crystal structure prediction

An anisotropic atom-atom distributed intermolecular force-field (DIFF) for rigid trinitrobenzene (TNB) is developed using distributed multipole moments, dipolar polarizabilities, and dispersion coefficients derived from the charge density of the isolated molecule. The short-range parameters of the force-field are fitted to first- and second-order symmetry-adapted perturbation theory dimer interaction energy calculations using the distributed density-overlap model to guide the parameterization of the short-range anisotropy. The second-order calculations are used for fitting the damping coefficients of the long-range dispersion and polarization and also for relaxing the isotropic short-range coefficients in the final model, DIFF-srL2(rel). We assess the accuracy of the unrelaxed model, DIFF-srL2(norel), and its equivalent without short-range anisotropy, DIFF-srL0(norel), as these models are easier to derive. The model potentials are contrasted with empirical models for the repulsion-dispersion fitted to organic crystal structures with multipoles of iterated stockholder atoms (ISAs), FIT(ISA,L4), and with Gaussian Distributed Analysis (GDMA) multipoles, FIT(GDMA,L4), commonly used in modeling organic crystals. The potentials are tested for their ability to model the solid state of TNB. The non-empirical models provide more reasonable relative lattice energies of the three polymorphs of TNB and propose more sensible hypothetical structures than the empirical force-field (FIT). The DIFF-srL2(rel) model successfully has the most stable structure as one of the many structures that match the coordination sphere of form III. The neglect of the conformational flexibility of the nitro-groups is a significant approximation. This methodology provides a step toward force-fields capable of representing all phases of a molecule in molecular dynamics simulations.

First-Principles Many-Body Nonadditive Polarization Energies from Monomer and Dimer Calculations Only: A Case Study on Water

The many-body polarization energy is the major source of nonadditivity in strongly polar systems such as water. This nonadditivity is often considerable and must be included, if only in an average manner, to correctly describe the physical properties of the system. Models for the polarization energy are usually parametrized using experimental data, or theoretical estimates of the many-body effects. Here we show how many-body polarization models can be developed for water complexes using data for the monomer and dimer only using ideas recently developed in the field of intermolecular perturbation theory and state-of-the-art approaches for calculating distributed molecular properties based on the iterated stockholder atoms (ISA) algorithm. We show how these models can be calculated, and we validate their accuracy in describing the many-body nonadditive energies of a range of water clusters. We further investigate their sensitivity to the details of the polarization damping models used. We show how our very best polarization models yield many-body energies that agree with those computed with coupled-cluster methods, but at a fraction of the computational cost.