Molecular Theory of Detonation Initiation: Insight from First Principles Modeling of the Decomposition Mechanisms of Organic Nitro Energetic Materials

Photodissociation Dynamics of the Highly Stable ortho-Nitroaniline Cation

0rtho-Nitroaniline (ONA) is a model for the insensitive high explosive 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) that shares strong hydrogen bonding character between adjacent nitro and amino groups. This work reports femtosecond time-resolved mass spectrometry (FTRMS) measurements and theoretical calculations that explain the high stability of the ONA cation compared with related nitroaromatic molecules. Ab initio calculations found that the lowest-lying electronic excited state of the ONA cation, D1, lies more than 2 eV above the ground state, and the energetic barriers to rearrangement and dissociation reactions exceed this D1 energy. These theoretical results were confirmed by FTRMS pump-probe measurements showing that (1) fragment ions represented less than 30% of the total ion yield when a 1014 W cm-2, 1300 nm, 20 fs pump pulse was used to ionize ONA; and (2) 3.1 eV (400 nm) photons were required to induce dissociation of the ONA cation. Stronger coupling between the ground D0 and excited D4 states of the ONA cation at the geometry of neutral ONA resulted in a transient enhancement of fragment ion yields at <300 fs pump-probe delay times, prior to relaxation of the ONA cation to its optimal geometry.

The effect of {O,N}=X⋯M={Ti,Zr,Hf} interactions on the sensitivity of CNO2 trigger bonds in FOX-7: Approach based on the QTAIM/EDA-NOCV analysis

The local chemical reactivity of FOX-7 (1,1-diamino-2,2-nitroethylene, also known as DADNE from DiAminoDiNitroEthylene) was elucidated through a quantitative study of the electrostatic potential on the molecular surface, topological analysis based on Bader's theory, and the EDA-NOCV method. Unlike (O2N)2CC(NH2)H2N⋯Cp2MCH3+ complexes, which exhibit both σ-donor and π-acceptor features, the situation is different concerning the (H2N)2CC(NO2)(O)NO⋯Cp2MCH3+ complexes, where both charge transfers correspond to the σ-donation. The two charge transfers reinforce each other, resulting in increased stability for (H2N)2CC(NO2)(O)NO⋯Cp2MCH3+. This seems to strengthen the (H2N)2CC(NO2)(O)NO⋯M={Ti,Zr,Hf} bond, which may explain the high stability of (H2N)2CC(NO2)(O)NO⋯Cp2MCH3+ compared to (O2N)2CC(NH2)-H2N⋯Cp2MCH3+. Results from topological analysis revealed that the decreased sensitivity to decomposition of CNO2 bonds depends on the chemical nature of the interacting metal, and the best achievements are obtained for the Hf-based complex. Our results demonstrate that the interaction of M={Ti,Zr,Hf} with CNO2 is more favourable than that with CNH2, this specific action on the trigger bond may support the use of Metallocene Methyl Cations (MMC) as possible neutralisers.

Investigation of the decomposition mechanism of MTNP melt-cast explosive at different temperatures and pressures through ReaxFF/lg molecular dynamics simulations

CONTEXT

Thermal decomposition of 1-methyl-3,4,5-trinitropyrazole (MTNP), a melt-cast explosive, was investigated at different temperatures (2500, 2750, 3000, 3250, and 3500 K) and pressures (3000 K/0.5 GPa, 3000 K/1 GPa) using the ReaxFF/lg force field. The study aimed to analyze the changes in reactant quantities, initial reaction pathways, and final product yields. The results demonstrated that complete decomposition of MTNP molecules occurred within a timeframe of 200 ps, with shorter decomposition times observed as the temperature increased. The high-temperature thermal decomposition of MTNP was found to follow two primary reaction pathways. Reaction 1 involved denitration, while reaction 2 proceeded with nitro group isomerization. DFT calculations indicated that nitro group isomerization was the most favorable reaction. During the initial stages, higher quantities of NO2, NO, and N2 were observed compared to other species. This can be attributed to the relatively higher nitrogen and oxygen content in the MTNP structure. Among the five reaction temperatures, it was observed that the quantities of small molecules followed the order of NO2 > NO > N2 > CO. Moreover, with increasing temperature, the quantities of all four small molecules increased, indicating that higher temperatures promoted the progression of the reactions. However, as the pressure increased, there was a trend of initially increasing and then decreasing to zero for the quantities of NO2 and NO. This suggests that high temperature accelerated the high-temperature thermal decomposition of NO2 and NO, leading to a significant increase in the content of N2.

METHODS

A 3 × 5 × 5 supercell model of MTNP was constructed in Materials Studio, consisting of 75 unit cells and 300 MTNP molecules. The model was then subjected to a 20 ps geometric optimization using the Polak-Ribiere version of the conjugate gradient (CG) algorithm in the large-scale atomic/molecular massively parallel simulator (LAMMPS) under the isothermal-isobaric (NPT) ensemble at 1 atm pressure and 300 K temperature. Following the optimization, molecular dynamics simulations were performed on the model at five temperatures (2500, 2750, 3000, 3200, and 3500 K) under 1 atm using the NPT ensemble for a total duration of 1 ns. During the simulations, atomic trajectories, as well as information on atomic and molecular species, were output every 500 steps. Subsequently, a custom script was utilized to analyze the thermal decomposition pathways and products. A time step of 0.1 fs was employed for the calculations, and periodic boundary conditions were applied to eliminate boundary effects.

Thermochemistry, Tautomerism, and Thermal Stability of 5,7-Dinitrobenzotriazoles

Nitro derivatives of benzotriazoles are safe energetic materials with remarkable thermal stability. In the present study, we report on the kinetics and mechanism of thermal decomposition for 5,7-dinitrobenzotriazole (DBT) and 4-amino-5,7-dinitrobenzotriazole (ADBT). The pressure differential scanning calorimetry was employed to study the decomposition kinetics of DBT experimentally because the measurements under atmospheric pressure are disturbed by competing evaporation. The thermolysis of DBT in the melt is described by a kinetic scheme with two global reactions. The first stage is a strong autocatalytic process that includes the first-order reaction (Ea1I = 173.9 ± 0.9 kJ mol−1, log(A1I/s−1) = 12.82 ± 0.09) and the catalytic reaction of the second order with Ea2I = 136.5 ± 0.8 kJ mol−1, log(A2I/s−1) = 11.04 ± 0.07. The experimental study was complemented by predictive quantum chemical calculations (DLPNO-CCSD(T)). The calculations reveal that the 1H tautomer is the most energetically preferable form for both DBT and ADBT. Theory suggests the same decomposition mechanisms for DBT and ADBT, with the most favorable channels being nitro-nitrite isomerization and C–NO2 bond cleavage. The former channel has lower activation barriers (267 and 276 kJ mol−1 for DBT and ADBT, respectively) and dominates at lower temperatures. At the same time, due to the higher preexponential factor, the radical bond cleavage, with reaction enthalpies of 298 and 320 kJ mol−1, dominates in the experimental temperature range for both DBT and ADBT. In line with the theoretical predictions of C–NO2 bond energies, ADBT is more thermally stable than DBT. We also determined a reliable and mutually consistent set of thermochemical values for DBT and ADBT by combining the theoretically calculated (W1-F12 multilevel procedure) gas-phase enthalpies of formation and experimentally measured sublimation enthalpies.

Theoretical Kinetic Studies on Thermal Decomposition of Glycerol Trinitrate and Trimethylolethane Trinitrate in the Gas and Liquid Phases

Glycerol trinitrate (NG) and trimethylolethane trinitrate (TMETN), as typical nitrate esters, are important energetic plasticizers in solid propellants. With the aid of high-precision quantum chemical calculations, the Rice-Ramsperger-Kassel-Marcus (RRKM)/master equation theory and the transition state theory have been employed to investigate the decomposition kinetics of NG and TMETN in the gas phase (over the temperature range of 300-1000 K and pressure range of 0.01-100 atm) and liquid phase (using water as the solvent). The continuum solvation model based on solute electron density (SMD) was used to describe the solvent effect. The thermal decomposition mechanism is closely relevant to the combustion properties of energetic materials. The results show that the RO-NO₂ dissociation channel overwhelmingly favors other reaction pathways, including HONO elimination for the decomposition of NG and TMETN in both the gas phase and liquid phase. At 500 K and 1 atm, the rate coefficient of gas phase decomposition of TMETN is 5 times higher than that of NG. Nevertheless, the liquid phase decomposition of TMETN is a factor of 5835 slower than that of NG at 500 K. The solvation effect caused by vapor pressure and solubility can be used to justify such contradictions. Our calculations provide detailed mechanistic evidence for the initial kinetics of nitrate ester decomposition in both the gas phase and liquid phase, which is particularly valuable for understanding the multiphase decomposition behavior and building detailed kinetic models for nitrate ester.

Control of Explosive Chemical Reactions by Optical Excitations: Defect-Induced Decomposition of Trinitrotoluene at Metal Oxide Surfaces

Interfaces formed by high energy density materials and metal oxides present intriguing new opportunities for a large set of novel applications that depend on the control of the energy release and initiation of explosive chemical reactions. We studied the role of structural defects at a MgO surface in the modification of electronic and optical properties of the energetic material TNT (2-methyl-1,3,5-trinitrobenzene, also known as trinitrotoluene, C₇H₅N₃O₆) deposited at the surface. Using density functional theory (DFT)-based solid-state periodic calculations with hybrid density functionals, we show how the control of chemical explosive reactions can be achieved by tuning the electronic structure of energetic compound at an interface with oxides. The presence of defects at the oxide surface, such as steps, kinks, corners, and oxygen vacancies, significantly affects interfacial properties and modifies electronic spectra and charge transfer dynamics between the oxide surface and adsorbed energetic material. As a result, the electronic and optical properties of trinitrotoluene, mixed with an inorganic material (thus forming a composite), can be manipulated with high precision by interactions between TNT and the inorganic material at composite interfaces, namely, by charge transfer and band alignment. Also, the electron charge transfer between TNT and MgO surface reduces the decomposition barriers of the energetic material. In particular, it is shown that surface structural defects are critically important in the photodecomposition processes. These results open new possibilities for the rather precise control over the decomposition initiation mechanisms in energetic materials by optical excitations.

Understanding Explosive Sensitivity with Effective Trigger Linkage Kinetics

We present a simple linear model for ranking the drop weight impact sensitivity of organic explosives that is based explicitly on chemical kinetics. The model is parameterized to specific heats of explosion, Q, and Arrhenius kinetics for the onset of chemical reactions that are obtained from gas-phase Born-Oppenheimer molecular dynamics simulations for a chemically diverse set of 24 molecules. Reactive molecular dynamics simulations sample all possible decomposition pathways of the molecules with the appropriate probabilities to provide an effective reaction barrier. In addition, the calculations of effective trigger linkage kinetics can be accomplished without prior physical intuition of the most likely decomposition pathways. We found that the specific heat of explosion tends to reduce the effective barrier for decomposition in accordance with the Bell-Evans-Polanyi principle, which accounts naturally for the well-known correlations between explosive performance and sensitivity. Our model indicates that sensitive explosives derive their properties from a combination of weak trigger linkages that react at relatively low temperatures and large specific heats of explosion that further reduce the effective activation energy.

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.

Anisotropic Reaction Properties for Different HMX/HTPB Composites: A Theoretical Study of Shock Decomposition

Plastic-bonded explosives (PBXs) consisting of explosive grains and a polymer binder are commonly synthesized to improve mechanical properties and reduce sensitivity, but their intrinsic chemical behaviors while subjected to stress are not sufficiently understood yet. Here, we construct three composites of β-HMX bonded with the HTPB binder to investigate the reaction characteristics under shock loading using the quantum-based molecular dynamics method. Six typical interactions between HMX and HTPB molecules are detected when the system is subjected to pressure. Although the initial electron structure is modified by the impurity states from HTPB, the metallization process for HMX does not significantly change. The shock decompositions of HMX/HTPB along the (100) and (010) surface are initiated by molecular ring dissociation and hydrogen transfer. The initial oxidations of C and H within HTPB possess advantages. As for the (001) surface, the dissociation is started with alkyl dehydrogenation oxidation, and a stronger hydrogen transfer from HTPB to HMX is detected during the following process. Furthermore, considerable fragment aggregation is observed, which mainly derives from the formation of new C-C and C-N bonds under high pressure. The effect of cluster evolution on the progression of the following reaction is further studied by analyzing the bonded structure and displacement rate.

Photolysis of the Insensitive Explosive 1,3,5-Triamino-2,4,6-trinitrobenzene (TATB)

The explosive 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) is of particular interest due to its extreme insensitivity to impact, shock and heat, while providing a good detonation velocity. To determine its fate under environmental conditions, TATB powder was irradiated with simulated sunlight and, in water, under UV light at 254 nm. The hydrolysis of particles submerged in neutral and alkaline solutions was also examined. We found that, by changing experimental conditions (e.g., light source, and mass and physical state of TATB), the intermediates and final products were slightly different. Mono-benzofurazan was the major transformation product in both irradiation systems. Two minor transformation products, the aci-nitro form of TATB and 3,5-diamino-2,4,6-trinitrophenol, were detected under solar light, while 1,3,5-triamino-2-nitroso-4,6-dinitrobenzene, 1,3,5-triamino-2,4-dinitrobenzene and mono-benzofuroxan were produced under UV light. The product identified as 3,5-diamino-2,4,6-trinitrophenol was identical to the one formed in the dark under alkaline conditions (pH 13) and in water incubated at either 50 °C or aged at ambient conditions. Interestingly, when only a few milligrams of TATB were irradiated with simulated sunlight, the aci-isomer and mono-benzofurazan derivative were detected; however, the hydrolysis product 3,5-diamino-2,4,6-trinitrophenol formed only much later in the absence of light. This suggests that the water released from TATB to form mono-benzofurazan was trapped in the interstitial space between the TATB layers and slowly hydrolyzed the relatively stable aci-nitro intermediate to 3,5-diamino-2,4,6-trinitrophenol. This environmentally relevant discovery provides data on the fate of TATB in surface environments exposed to sunlight, which can transform the insoluble substrate into more soluble and corrosive derivatives, such as 3,5-diamino-2,4,6-trinitrophenol, and that some hydrolytic transformation can continue even without light.

Molecular Dynamics Simulations of the Thermal Decomposition of 3,4-Bis(3-nitrofurazan-4-yl)furoxan

When stimulated, for example, by a high temperature, the physical and chemical properties of energetic materials (EMs) may change, and, in turn, their overall performance is affected. Therefore, thermal stability is crucial for EMs, especially the thermal dynamic behavior. In the past decade, significant efforts have been made to study the thermal dynamic behavior of 3,4-bis(3-nitrofurazan-4-yl)furoxan (DNTF), one of the new high-energy-density materials (HEDMs). However, the thermal decomposition mechanism of DNTF is still not specific or comprehensive. In this study, the self-consistent-charge density-functional tight-binding method was combined with molecular dynamics (MD) simulations to reveal the differences in the thermal decomposition of DNTF under four heating conditions. The O-N (O) bond would fracture first during DNTF initial thermal decomposition at medium and low temperatures, thus triggering the cracking of the whole structure. At 2000 and 2500 K, NO₂ loss on outer ring I is the fastest initial thermal decomposition pathway, and it determines that the decomposition mechanism is different from that of a medium-low temperature. NO₂ is found to be the most active intermediate product; large molecular fragments, such as C₂N₂O, are found for the first time. Hopefully, these results could provide some insights into the decomposition mechanism of new HEDMs.

Recruiting Perovskites to Degrade Toxic Trinitrotoluene

Everybody knows TNT, the most widely used explosive material and a universal measure of the destructiveness of explosions. A long history of use and extensive manufacture of toxic TNT leads to the accumulation of these materials in soil and groundwater, which is a significant concern for environmental safety and sustainability. Reliable and cost-efficient technologies for removing or detoxifying TNT from the environment are lacking. Despite the extreme urgency, this remains an outstanding challenge that often goes unnoticed. We report here that highly controlled energy release from explosive molecules can be accomplished rather easily by preparing TNT-perovskite mixtures with a tailored perovskite surface morphology at ambient conditions. These results offer new insight into understanding the sensitivity of high explosives to detonation initiation and enable many novel applications, such as new concepts in harvesting and converting chemical energy, the design of new, improved energetics with tunable characteristics, the development of powerful fuels and miniaturized detonators, and new ways for eliminating toxins from land and water.

Predicting the impact sensitivity of a polymorphic high explosive: the curious case of FOX-7

The impact sensitivity (IS) of FOX-7 polymorphs is predicted by phonon up-pumping to decrease as layers of FOX-7 molecules flatten. Experimental validation proved anomalous owing to a phase transition during testing, raising questions regarding impact sensitivity measurement and highlighting the need for models to predict IS of polymorphic energetic materials.

Thermal Decomposition of Gas-Phase Bis(1,2,4-oxadiazole)bis(methylene) Dinitrate (BODN): A CCSD(T)-F12/DFT-Based Study of Reaction Pathways

Electronic structure methods based on density functional theory and coupled-cluster theory were employed to characterize elementary steps for the gas-phase thermal decomposition of bis(1,2,4-oxadiazole)bis(methylene) dinitrate (BODN). As typically found for nitrate ester-functionalized compounds, NO₂ and HONO eliminations were the most energetically favorable unimolecular paths for the parent molecule's decomposition. From there, sequences of unimolecular reactions for daughters of the initiation steps were postulated and characterized. For intermediates found to have barriers to unimolecular decomposition that would make their rate at the temperatures and time scales of interest negligible, their decomposition via H-atom abstraction and radical-addition reactions was characterized. Creating a comprehensive network that can be employed to develop a detailed finite-rate chemical kinetics mechanism for simulating BODN's decomposition, the results provide a basis for modeling BODN's combustion, as well as its response to thermal loads germane to its aging, storage, and handling.

Predicting the impact sensitivities of energetic materials through zone-center phonon up-pumping

The development of new energetic materials (EMs) is accompanied by significant hazards, prompting interest in their computational design. Before reliable in silico design strategies can be realized, however, approaches to understand and predict EM response to mechanical impact must be developed. We present here a fully ab initio model based on phonon up-pumping that successfully ranks the relative impact sensitivity of a series of organic EMs. The methodology depends only on the crystallographic unit cell and Brillouin zone center vibrational frequencies. We, therefore, expect this approach to become an integral tool in the large-scale screening of potential EMs.

Recent Advances in Synthesis and Properties of Nitrated-Pyrazoles Based Energetic Compounds

Nitrated-pyrazole-based energetic compounds have attracted wide publicity in the field of energetic materials (EMs) due to their high heat of formation, high density, tailored thermal stability, and detonation performance. Many nitrated-pyrazole-based energetic compounds have been developed to meet the increasing demands of high power, low sensitivity, and eco-friendly environment, and they have good applications in explosives, propellants, and pyrotechnics. Continuous and growing efforts have been committed to promote the rapid development of nitrated-pyrazole-based EMs in the last decade, especially through large amounts of Chinese research. Some of the ultimate aims of nitrated-pyrazole-based materials are to develop potential candidates of castable explosives, explore novel insensitive high energy materials, search for low cost synthesis strategies, high efficiency, and green environmental protection, and further widen the applications of EMs. This review article aims to present the recent processes in the synthesis and physical and explosive performances of the nitrated-pyrazole-based Ems, including monopyrazoles with nitro, bispyrazoles with nitro, nitropyrazolo[4,3-c]pyrazoles, and their derivatives, and to comb the development trend of these compounds. This review intends to prompt fresh concepts for designing prominent high-performance nitropyrazole-based EMs.

Theoretical Investigation of Energetic Salts with Pentazolate Anion

Energetic salts based on pentazolate anion (cyclo-N₅⁻) have attracted much attention due to their high nitrogen contents. However, it is an enormous challenge to efficiently screen out an appropriate cation that can match well with cyclo-N₅⁻. The vertical electron affinity (VEA) of the cations and vertical ionization potential (VIP) of the anions for 135 energetic salts and some cyclo-N₅⁻ salts were calculated by the density functional theory (DFT). The magnitudes of VEA and VIP, and their matchability were analyzed. The results based on the calculations at the B3LYP/6-311++G(d,p) and B3LYP/aug-cc-pVTZ levels indicate that there is an excellent compatibility between cyclo-N₅⁻ and cation when the difference between the VEA of cation and the VIP of cyclo-N₅⁻ anion is −2.8 to −1.0 eV. The densities of the salts were predicted by the DFT method. Relationship between the calculated density and the experimental density was established as ρ_Expt = 1.111ρ_cal − 0.06067 with a correlation coefficient of 0.905. This regression equation could be in turn used to calibrate the calculated density of the cyclo-N₅⁻ energetic salts accurately. This work provides a favorable way to explore the energetic salts with excellent performance based on cyclo-N₅⁻.

First-principle calculations of electronic, vibrational, and thermodynamic properties of 1,3-diamino-2,4,6-trinitrobenzene

Energy-containing materials have aroused people's widespread concern because of its admirable performance in recent years. In this paper, the electronic structure, vibrational, and thermodynamic properties of 1,3-diamino-2,4,6-trinitrobenzene (DATB) are systematically investigated by adopting the first-principle calculations. We find that lattice parameters are in excellent agreement with the previous calculated and experimental values. The vibration spectra are described in detail and the peaks in the Raman and infrared spectra are assigned to different vibration modes. Phonon dispersion curves indicate that the DATB is dynamically stable. According to the vibrational properties, the thermodynamic functions such as enthalpy (H), constant volume heat capacity (C_V), Helmholtz free energy (F), Debye temperature (Θ), and entropy (S) are analyzed. No corresponding experimental values have been found so far, and therefore, knowledge of these properties will provide a reference and guidance for the follow-up research.

Dissociative adsorption modes of TATB on the Al (111) surface: a DFT investigation

Herein, the adsorption modes and electronic structures of TATB/Al (111) systems were investigated using the density functional theory (DFT) approach. We found that chemical adsorption led to the decomposition of the TATB molecule on the Al surface by four adsorption modes. All the adsorption configurations were accompanied by fractures of the N-O bonds in the nitro groups. In addition, there was a hydrogen atom transfer for 5P. For parallel and vertical adsorptions, the TATB molecules favored planar or quasi-planar structures. The order of total energy with BSSE correction matches well with the order of adsorption energy. The absolute values of energy and adsorption energy of 6P and 6V are highest in the parallel and vertical adsorption systems, respectively. Electrons are transferred from the Al (111) surface to the TATB molecule; this results in the activation of TATB on the Al (111) surface and obvious augmentation of the PDOS (partial density of states) peaks of the N and O atoms. From the Al (111) surface to the TATB molecule, the transfer of the electrons of 4P (14.00e) and 6V (9.04e) is largest for the parallel and vertical adsorptions, respectively.

Reaction kinetic properties of 1,3,5-triamino-2,4,6-trinitrobenzene: a DFTB study of thermal decomposition

C–NH₂bond breakage and N–N bond formation are the rate-controlling steps for TATB thermal decomposition.

Vibrationally induced metallisation of the energetic azide α-NaN3

As initiation of an energetic material requires rupture of a covalent bond, and therefore population of antibonding electronic states, consideration of the electronic band gap has dominated initiation mechanisms for solid state materials. Most prominent are models based on metallisation, where static mechanical perturbation leads to closing of the electronic band gap. This work explores an alternative mechanism for the dynamic metallisation of a model energetic material, where vibrational excitation resulting from mechanical impact is found to induce transient metallisation of α-NaN3. The normal coordinates associated with bending the azido anion close the electronic band gap, facilitating the formation of highly reactive species important for initiation of energetic materials. The DFT simulated vibrational spectrum of α-NaN3 exhibits excellent reproduction of the experimental low-temperature inelastic neutron scattering spectrum (INS).

Predicting the Initial Thermal Decomposition Path of Nitrobenzene Caused by Mode Vibration at Moderate-Low Temperatures: Temperature-Dependent Anti-Stokes Raman Spectra Experiments and First-Principals Calculations

The lack of understanding of the initial decomposition micromechanism of energetic materials subjected to external stimulation has hindered its safe storage, usage, and development. The initial thermal decomposition path of nitrobenzene triggered by molecular thermal motion is investigated using temperature-dependent anti-Stokes Raman spectra experiments and first-principles calculations to clarify the initial thermal decomposition micromechanism. The experiment shows that the symmetric nitro stretching, antisymmetric nitro stretching, and phenyl ring stretching vibration modes are active as increasing temperature below 500 K. The DFT method is used to examine the effects of the three mode vibrations on the initial decomposition of nitrobenzene by relaxed scan for each relevant change in bond lengths and bond angles to obtain the optimal reaction channel leading to initial thermal decomposition of nitrobenzene. The results demonstrate that the initial thermal decomposition is the isomerization of nitrobenzene to phenyl nitrite. The optimal reaction channel leading to the initial isomerization is the increase or decrease of angle O-N-C from the antisymmetric nitro stretching vibration, which causes the torsion of nitro group and the subsequent oxygen atom attacking carbon atom. The scanning energy barrier related to angle O-N-C is about 62.1 kcal/mol, which is very consistent with the calculated activation barrier of isomerization of nitrobenzene. This proves the reliability of our conclusions.

Trigger bond analysis of nitroaromatic energetic materials using wiberg bond indices

The identification of trigger bonds, bonds that break to initiate explosive decomposition, using computational methods could help direct the development of novel, "green" and efficient high energy density materials (HEDMs). Comparing bond densities in energetic materials to reference molecules using Wiberg bond indices (WBIs) provides a relative scale for bond activation (%ΔWBIs) to assign trigger bonds in a set of 63 nitroaromatic conventional energetic molecules. Intramolecular hydrogen bonding interactions enhance contributions of resonance structures that strengthen, or deactivate, the CNO₂ trigger bonds and reduce the sensitivity of nitroaniline-based HEDMs. In contrast, unidirectional hydrogen bonding in nitrophenols strengthens the bond to the hydrogen bond acceptor, but the phenol lone pairs repel and activate an adjacent nitro group. Steric effects, electron withdrawing groups and greater nitro dihedral angles also activate the CNO₂ trigger bonds. %ΔWBIs indicate that nitro groups within an energetic molecule are not all necessarily equally activated to contribute to initiation. %ΔWBIs generally correlate well with impact sensitivity, especially for HEDMs with intramolecular hydrogen bonding, and are a better measure of trigger bond strength than bond dissociation energies (BDEs). However, the method is less effective for HEDMs with significant secondary effects in the solid state. Assignment of trigger bonds using %ΔWBIs could contribute to understanding the effect of intramolecular interactions on energetic properties. © 2018 Wiley Periodicals, Inc.

Bimolecular Reactions between Dimethylnitramine and Its Radical Decomposition Products

Bimolecular reactions between intact nitramines and their radical decomposition products can accelerate thermal decomposition, yet the detailed mechanisms of such reactions are not well understood. We have used density functional theory at the M06/6-311++G(3df,3pd) level to locate transition structures and compute 0 K activation barriers for various gas-phase reactions that may contribute to radical-assisted decomposition of dimethylnitramine (DMNA, (CH₃)₂NNO₂). Our calculations indicate that H abstraction from DMNA is the lowest-barrier mechanism for most radicals and a subsequent N-N β-scission in the alkyl radical 3 leads to an imine intermediate and NO₂. H abstraction is thus responsible for conversion of most radicals to NO₂. Also, among the nine radicals considered, NO is found to be least reactive and its reactions with DMNA yield dimethylnitrosoamine (DMNSA, (CH₃)₂NNO), a known product of DMNA decomposition.