Predicting Trigger Bonds in Explosive Materials through Wiberg Bond Index Analysis

Understanding the explosive decomposition pathways of high‐energy‐density materials (HEDMs) is important for developing compounds with improved properties. Rapid reaction rates make the detonation mechanisms of HEDMs difficult to understand, so computational tools are used to predict trigger bonds—weak bonds that break, leading to detonation. Wiberg bond indices (WBIs) have been used to compare bond densities in HEDMs to reference molecules to provide a relative scale for the bond strength to predict the activated bonds most likely to break to trigger an explosion. This analysis confirms that X?NO₂ (X=N,C,O) bonds are trigger linkages in common HEDMs such as TNT, RDX and PETN, consistent with previous experimental and theoretical studies. Calculations on a small test set of substituted tetrazoles show that the assignment of the trigger bond depends upon the functionality of the material and that the relative weakening of the bond correlates with experimental impact sensitivities.     

A density functional theory-time-dependent density functional theory investigation of photo-induced hydrogen bond and proton transfer for 2-(3,5-dichloro-2,6-dihydroxy-phenyl)-benzoxazole-6-carboxylicacid

Given the paramount importance of excited-state relaxation in the photochemical process, excited-state hydrogen bonding interactions and excited-state intramolecular proton transfer (ESIPT) are always hot topics. In this work, we theoretically explore the excited-state dynamical behaviors for a novel 2-(3,5-dichloro-2,6-dihydroxy-phenyl)-benzoxazole-6-carboxylicacid (DDPBC) system. As two intramolecular hydrogen bonds (O1 H2⋯N3 and O4 H5⋯O6) exist in the DDPBC structure, we first check if the double proton transfer form cannot be formed in the S1 state. Then, we explore the changes of geometrical parameters involved in hydrogen bonds, based on which we confirm that the dual intramolecular hydrogen bonds are strengthened on photo-excitation. The O1 H2⋯N3 hydrogen bond particularly plays a more important role in excited state. When it comes to the photo-induced excitation, we find charge transfer and electronic density redistribution around O1 H2 and N3 atom moieties. We verify the ESIPT tendency arising from the O1 H2⋯N3 hydrogen bond. In the analysis of the potential energy curves, along with O1 H2⋯N3 and O4 H5⋯O6, we demonstrate that the ESIPT reaction should occur along with O1 H2⋯N3 rather than O4 H5⋯O6. This work not only clarifies the specific ESIPT mechanism for DDPBC system but also paves the way for further novel applications based on DDPBC structure in the future.     

alpha- and beta-FOX-7, polymorphs of a high energy density material, studied by X-ray single crystal and powder investigations in the temperature range from 200 to 423 K

The alpha-beta phase transition in the novel energetic material 1,1-diamino-2,2-dinitroethylene, C2H4N4O4 (FOX-7), has been studied by single-crystal X-ray investigations at five different temperatures over the 200-393 K range. In these investigations, the positions of the hydrogen atoms were experimentally determined without any geometric constraints. In addition, X-ray powder investigations using the Guinier technique have been performed to characterize the beta-phase up to 423 K. The alpha-beta phase transition at 389 K is first order, shows a discontinuous increase of the molar volume and entropy (DeltaV = 1.75 cm3/mol, X-ray investigation; DeltaS = 1.5 cal/K mol, DSC analysis), and can be classified as displacive. The hitherto unknown structure of beta-FOX-7 was solved at 393 K and showed simple structural relations to the alpha-polymorph. The characteristic bonding in wave-shaped layers is now found for beta-FOX-7 (P2(1)2(1)2(1), z = 4, a= 6.9738(7) A, b = 6.635(1) A, c = 11.648(2) A, 393 K), as well as for alpha-FOX-7 (P2(1)/n, z = 4, a = 6.9467(7) A, b = 6.6887(9) A, c = 11.350(1) A, beta = 90.143(13) degrees , 373 K). Interestingly, whereas the intramolecular C-C, C-N, N-O, and N-H bond distances remain nearly unchanged for both polymorphs over the whole temperature range from 200 to 393 K, the two nitro groups deviate strongly from the molecular plane formed by the two carbon and two amino nitrogen atoms. In alpha-FOX-7 at 373 K, the nitro groups are twisted -47 and +6 degrees with respect to the carbon-carbon bond, but in beta-FOX-7 at 393 K, these twist angles are changed to -36 and +20 degrees . Within the layers, the FOX-7 molecules show strong pi-conjugation and extensive intra- and intermolecular hydrogen bonding. In this investigation, we have been able to show that alpha- and beta-FOX-7 build up different nets of intermolecular hydrogen bonds. In alpha-FOX-7, each oxygen atom of the nitro groups is involved in two hydrogen bonds resulting in two intramolecular and six intermolecular hydrogen bonds. But in beta-FOX-7 this coordination changes, and half of the oxygen atoms build up two and the other half build up three hydrogen bonds leading to two intramolecular and eight intermolecular hydrogen bonds. The average intermolecular hydrogen bond distance increases slightly from 2.31 A in alpha-FOX-7 to 2.52 A in beta-FOX-7. The C-NO2 bonds are of particular interest because they are referred to as the detonation trigger. It has been suggested that these bonds could be strengthened by the extensive intermolecular hydrogen bonding within the layers in both polymorphs. Such bond strengthening via cooperative effects was proposed in earlier DFT calculations on FOX-7 and may be one key to understanding its low sensitivity and high activation energy to impact.     

DFT investigation on detonation properties and sensitivities of bridged triazolo[4,5-d]pyridazine based energetic materials

A series of bridged triazolo[4,5-d]pyridazine based energetic materials were optimized at B3LYP/6-311G(d, p) level of density functional theory (DFT), and their detonation properties and sensitivities were calculated. The results show that the  NN bridge/ N₃ group were beneficial to improve values of heats of formation while  NN bridge/ C(NO₂)₃ group can improve detonation properties remarkably. In view of the sensitivities, compound F2 possesses the minimum values of impact sensitivity which reveals that  NHNH bridge/ C(NO₂)₃ group will decrease the stability of the designed compounds. Take both of detonation properties and sensitivities into consideration, compounds C8, E7, E8, F8 were screened as candidates of potential energetic materials since these compounds possess similar detonation properties and sensitivities values to those of RDX. All the calculated results were except to shine lights on the design and synthesis of novel high energy density materials.     

Theoretical studies on polynitrobicyclo[1.1.1]pentanes in search of novel high energy density materials

Bicyclo[1.1.1]pentane is a highly strained hydrocarbon system due to close proximity of nonbonded bridge head carbons. Based on fully optimized molecular geometries at the density functional theory using the B3LYP/6-31G* level, densities, detonation velocities, and pressures for a series of polynitrobicyclo[1.1.1]pentanes, as well as their thermal stabilities were investigated in search for high energy density materials (HEDMs). The designed compounds with more than two nitro groups are characterized by high heat of formation and magnitude correlative with the number and space distance of nitro groups. Density was calculated using the crystal packing calculations and an increase in the number of nitro groups increases the density. The increase in density shows a linear increase in the detonation characteristics. Bond dissociation energy was analyzed to determine thermal stability. Calculations of the bond length and bond dissociation energies of the C-NO₂ bond indicate that this may be the possible trigger bond in the pyrolysis mechanism. 1,2,3-Trinitrobicyclo[1.1.1]pentane (S3), 1,2,3,4-tetranitrobicyclo[1.1.1]pentane (S4), and 1,2,3,4,5-pentanitrobicyclo[1.1.1]pentane (S5) have better energetic characteristics with better stability and insensitivity, and as such may be explored in defense applications as promising candidates of the HEDMs series.     

Theoretical calculation of nitro-1-(2,4,6-trinitrophenyl)-1H-azoles energetic compounds

In order to study the properties of new energetic compounds formed by introducing nitroazoles into 2,4,6-trinitrobezene, the density, heat of formation and detonation properties of 36 nitro-1-(2,4,6-trinitrobenzene)-1H-azoles energetic compounds are studied by density functional theory, and their stability and melting point are predicted. The results show that most of target compounds have good detonation properties and stability. And it is found that nitro-1-(2,4,6-Trinitrophenyl)-1H-pyrrole compounds and nitro-1-(2,4,6-trinitrop-enyl)-1H-Imidazole compounds have good thermal stability, and their weakest bond is C NO₂ bond, the bond dissociation energy of the weakest bond is 222–238 kJ mol⁻¹ and close to 2,4,6-trinitrotoluene (235 kJ mol⁻¹). The weakest bond of the other compounds may be the C NO₂ bond or the N N bond, and the strength of the N N bond is related to the nitro group on azole ring.     

Molecular Dynamics Simulations for Effects of Fluoropolymer Binder Content in CL-20/TNT Based Polymer-Bonded Explosives

In order to better understand the role of binder content, molecular dynamics (MD) simulations were performed to study the interfacial interactions, sensitivity and mechanical properties of 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane/2,4,6-trinitrotoluene (CL-20/TNT) based polymer-bonded explosives (PBXs) with fluorine rubber F₂₃₁₁. The binding energy between CL-20/TNT co-crystal (1 0 0) surface and F₂₃₁₁, pair correlation function, the maximum bond length of the N–NO₂ trigger bond, and the mechanical properties of the PBXs were reported. From the calculated binding energy, it was found that binding energy increases with increasing F₂₃₁₁ content. Additionally, according to the results of pair correlation function, it turns out that H–O hydrogen bonds and H–F hydrogen bonds exist between F₂₃₁₁ molecules and the molecules in CL-20/TNT. The length of trigger bond in CL-20/TNT were adopted as theoretical criterion of sensitivity. The maximum bond length of the N–NO₂ trigger bond decreased very significantly when the F₂₃₁₁ content increased from 0 to 9.2%. This indicated increasing F₂₃₁₁ content can reduce sensitivity and improve thermal stability. However, the maximum bond length of the N–NO₂ trigger bond remained essentially unchanged when the F₂₃₁₁ content was further increased. Additionally, the calculated mechanical data indicated that with the increase in F₂₃₁₁ content, the rigidity of CL-20/TNT based PBXs was decrease, the toughness was improved.     

Computational study of energetic nitrogen‐rich derivatives of 1,1′‐ and 5,5′‐bridged ditetrazoles

Density functional theory method was used to study the heats of formation (HOFs), electronic structure, energetic properties, and thermal stability for a series of bridged ditetrazole derivatives with different linkages and substituent groups. The results show that the ? N₃ group and azo bridge (? N?N? ) play a very important role in increasing the HOF values of the ditetrazole derivatives. The effects of the substituents on the HOMO–LUMO gap are combined with those of the bridge groups. The calculated detonation velocities and detonation pressures indicate that the ? NO₂, ? NF₂, ? N?N? , or ? N(O)?N? group is an effective structural unit for enhancing the detonation performance for the derivatives. An analysis of the bond dissociation energies for several relatively weak bonds suggests that the N? N bond in the ring or outside the ring is the weakest one and the N? N cleavage is possible to happen in thermal decomposition. Overall, the ? CH₂? CH₂? or ? NH? NH? group is an effective bridge for enhancing the thermal stability of the bridged ditetrazoles. Because of their desirable detonation performance and thermal stability, five compounds may be considered as the potential candidates of high‐energy density materials (HEDMs). These results provide basic information for the molecular design of novel HEDMs. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Theoretical study on initial decomposition paths of energetic materials featuring RDX ring

The molecular properties of RDX are affected by the introduction of different functional groups, and the decomposition process of these analogues is studied in this paper. DFT method is used to study the initial decomposition reaction paths of 30 high energy materials based RDX skeleton. In the nitro cleavage reaction, the energy barrier become relatively low by introducing CH(NO₂)₂ or  C(NO₂)₃ groups on the C site of the six membered ring. In the ring opening reaction, the ring opening process is easier to proceed by introducing  NH₂ or  NHNH₂ groups on the C site of the six membered ring.     

Density functional theory studies of effects of boron replacement on the structure and property of RDX and HMX

In this study, based on two model nitramine compounds hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) and octahydro-1,3,5,7-tetranitro-1,3,5, 7-tetrazocine (HMX), two series of new energetic molecules were designed by replacing carbon atoms in the ring with different amounts of boron atoms, their structures and performances were investigated theoretically by the density functional theory method. The results showed that the boron replacement could affect the molecular shape and electronic structure of RDX and HMX greatly, and then would do harm to the main performance like the heat of formation, density, and sensitivity. However, the compound RDX-B2 is an exception; it was formed by replacing two boron atoms into the system of RDX and has the symmetric boat-like structure. Its oxygen balance (4.9%), density (1.91 g/cm³), detonation velocity (8.85 km/s), and detonation pressure (36.9 GPa) are all higher than RDX. Furthermore, RDX-B2 has shorter and stronger N NO₂ bonds than RDX, making it possesses lower sensitivity (45 cm) and better thermal stability (the bond dissociation energy for the N NO₂ bond is 204.7 kJ/mol) than RDX. Besides, RDX-B1 and HMX-B4 also have good overall performance; these three new molecules may be regarded as a new potential candidate for high energy density compounds.     

Computational study of structures and proton transfer in hydrogen‐bonded ammonia complexes using semiempirical valence‐bond approach

In this study, the seGVB method was implemented for the N H bonding system, specifically for hydrogen‐bonded ammonia complexes, and the model well reproduces the MP2 geometries and energetics. A comparison between the ammonia dimer and water dimer is given from the viewpoint of valance‐bond structures in terms of the calculated bond energies and pair–pair interactions. The linear hydrogen bond is found to be stronger than the bent bonds in both cases, with the difference in energy between the linear and cyclic structures being comparable in both cases although the NH bonds are generally weaker. The energy decomposition clearly demonstrates that the changes in electronic energy are quite different in the two cases due to the presence of an additional lone pair on the water molecule, and it is this effect which leads to the net stabilization of the cyclic structure for the ammonia dimer. Proton‐transfer profiles for hydrogen‐bonded ammonia complexes [NH₂ H NH₂]⁻ and [NH₃ H NH₃]⁺ were calculated. The barrier for proton transfer in [NH₃ H NH₃]⁺ is larger than that in [NH₂ H NH₂]⁻, but smaller than that in the protonated water dimer. The different bonding structures substantially affect the barrier to proton transfer, even though they are isoelectronic systems. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 73: 357–367, 1999     

Energy barriers to gas-phase unimolecular decomposition of trinitrotoluenes

Alternative versions of gas-phase unimolecular decomposition of six isomeric trinitrotoluenes, in particular homolytic dissociation of C_arom–NO₂ and C_arom–CH₃ bonds, nitro–nitrite rearrangement, intramolecular hydrogen transfer from the methyl group to nitro group with formation of aci-trinitrotoluenes, and formation of various bicyclic intermediates, have been simulated at the B3LYP/6-31+G(2df,p) level of theory. Except for 3,4,5-trinitrotoluene, the most energetically favorable for all other examined trinitrotoluenes is intramolecular hydrogen transfer. 3,4,5-Trinitrotoluene preferentially decomposes via formation of [6 + 4]-bicyclic intermediates or homolytic dissociation of the C_arom–NO₂ bond.     

Nitrification Progress of Nitrogen-Rich Heterocyclic Energetic Compounds: A Review

As a momentous energetic group, a nitro group widely exists in high-energy-density materials (HEDMs), such as trinitrotoluene (TNT), 1,3,5-triamino-2,4,6-trinitrobenzene (TATB), cyclo-1,3,5-trimethylene-2,4,6-trinitramine (RDX), etc. The nitro group has a significant effect on improving the oxygen balance and detonation performances of energetic materials (EMs). Moreover, the nitro group is a strong electron-withdrawing group, and it can increase the acidity of the acidic hydrogen-containing nitrogen-rich energetic compounds to facilitate the construction of energetic ionic salts. Thus, it is possible to design nitro-nitrogen-rich energetic compounds with adjustable properties. In this paper, the nitration methods of azoles, including imidazole, pyrazole, triazole, tetrazole, and oxadiazole, as well as azines, including pyrazine, pyridazine, triazine, and tetrazine, have been concluded. Furthermore, the prospect of the future development of nitrogen-rich heterocyclic energetic compounds has been stated, so as to provide references for researchers who are engaged in the synthesis of EMs.     

Molecular design of trinitromethyl-substituted nitrogen-rich heterocycle derivatives with good oxygen balance as high-energy density compounds

Density functional theory method was used to study the heats of formation (HOFs), electronic structure, energetic properties, and pyrolysis mechanism of a series of trinitromethyl-substituted heterocycle (including triazole, tetrazole, furazan, tetrazine, and fused heterocycles) derivatives. It is found that the fused ring, tetrazine, and tetrazole are effective structural units for increasing the HOFs of the derivatives. The substitution of the combination of nitro and trinitromethyl is very useful for improving their HOFs. The calculated energetic properties indicate that the combination of the nitro and trinitromethyl is very helpful for improving their detonation properties and oxygen balances (OB). Most of the title compounds have a good OB over zero. The OB of six compounds are very high and over 22. An analysis of the bond dissociation energies for several relatively weak bonds suggests that the N–O bond in the ring is a trigger bond for BIII-1, CI-3, and CI-4, and the ring–NO₂ and (NO₂)₂C–NO₂ bond cleavage is likely to happen in thermal decomposition for the remaining compounds. Considering the detonation performance and thermal stability, seven compounds could be regarded as potential candidates for high-energy compounds. Four compounds may be used as the novel high-energy oxidizers.     

The effect of intramolecular interactions on hydrogen bond acidity

DFT calculations on a range of molecules containing intramolecular hydrogen bonds are reported, with a view to establishing how intramolecular hydrogen bonding affects their intermolecular interactions. It is shown that properties such as the energy of the intramolecular H-bond are unrelated to the ability to form external H-bonds. Conversely, several properties of complexes with a reference base correlate well with an experimental scale of H-bond acidity, and accurate predictive models are determined. A more detailed study, using electrostatic and overlap properties of complexes with a reference base, is used to predict the location, as well as strength, of hydrogen bond acidity. The effects of intramolecular hydrogen bonding on acidity can be seen not just on O-H and N-H, where acidity is greatly reduced, but also on certain C-H groups, which in some cases become the primary source of acidity.     

N‐Diazo‐Bridged Nitroazoles: Catenated Nitrogen‐Atom Chains Compatible with Nitro Functionalities

N‐diazo‐bridged azoles were synthesized based on oxidative coupling of N‐aminoazoles. Incorporation of extended catenated nitrogen‐atom chains with nitro groups led to compounds with favorable functional compatibilities. This combination gives rise to a series of high‐density energetic materials (HEDMs) with high heats of formation, enhanced densities, positive oxygen balances, and good detonation properties while retaining excellent thermal stabilities and relatively low impact sensitivities. Calculated and experimental studies showed the delicate balance between the length of the nitrogen atom chain, energetic performance, and inherent stability, thus, providing a promising strategy for designing advanced energetic materials.     

(E)‐Methyl 2‐anilinomethylene‐3‐oxobutanoate: X‐ray and density functional theory studies

Molecules of the title compound, C₁₂H₁₃NO₃, are not planar and are stabilized by electrostatic interactions, as the dipole moment of the molecule is 3.76 D. They are also stabilized by intramolecular hydrogen bonds of N...O and C...O types, and by a complicated network of weak intermolecular hydrogen bonds of the C...O type. This paper also reports the theoretical investigation of the hydrogen bonding and electronic structure of the title compound using natural bond orbital (NBO) analysis.     

A theoretical charge density study on nitrogen-rich 4,4′,5,5′-tetranitro-2,2′-bi-1H-imidazole (TNBI) energetic molecule

The bond topological and electrostatic properties of nitrogen-rich 4,4′,5,5′-tetranitro-2,2′-bi-1H-imidazole (TNBI) energetic molecule have been calculated from the DFT method with the basis set 6-311G** and the AIM theory. The optimized geometry of this molecule is almost matched with the experimental geometric parameters. The electron density at the bond critical point and the Laplacian of electron density of C–NO₂ bonds are not equal, one of them is much weaker than the other. Similar trend exists in the C–N bonds of the imidazole ring of the molecule. The ratio of the bond dissociation energy (BDE) of the weakest bond to the molecular total energy exhibits nearly a linear correlation with the impact sensitivity; its h ₅₀% value is ~32.01 cm. The electrostatic potential around both the nitro groups are found unequal; the NO₂ group of weakest C–NO₂ bond exhibits an extended electronegative region.     

New theoretically predicted RDX‐ and β‐HMX‐based high‐energy‐density molecules

Theoretically new high‐energy‐density materials (HEDM) in which the hydrogens on RDX and β‐HMX (hexahydro‐1,3,5‐trinitro‐1,3,5‐triazine and octahydro‐1,3,5,7‐tetranitro‐1,3,5,7‐tetrazocine, respectively) were sequentially replaced by (N NO₂)x functional groups were designed and evaluated using density functional theory calculations in combination with the Kamlet–Jacobs equations and an atoms‐in‐molecules (AIM) analysis. Improved detonation properties and reduced sensitivity compared to RDX and β‐HMX were predicted. Interestingly, the RDX and β‐HMX derivatives having one attached N NO₂ group [RDX‐(NNO₂)1 and HMX‐(NNO₂)1] showed excellent detonation properties (detonation velocities: 9.529 and 9.575 km·s⁻¹, and detonation pressures: 40.818 and 41.570 GPa, respectively), which were superior to the parent compounds. Sensitivity estimations obtained by calculating impact sensitivities and HOMO‐LUMO gaps indicated that RDX‐(NNO₂)1 and HMX‐(NNO₂)1 were less stable than RDX and HMX but more stable than any of the other derivatives. This method of sequential NNO₂ group attachment on conventional HEDMs offers a firm basis for further studies on the design of new explosives. Furthermore, the newly found structures may be promising candidates for better HEDMs.     

Polymorphs of gabapentin

Gabapentin [or 1‐(aminomethyl)cyclohexaneacetic acid], C₉H₁₇NO₂, exists as a zwitterion [1‐(ammoniomethyl)cyclohexaneacetate] in the solid state. The crystal structures and bonding networks of two new monoclinic polymorphs (β‐gabapentin and γ‐gabapentin) are studied and compared with a previously reported gabapentin polymorph [α‐gabapentin: Ibers (2001). Acta Cryst. C 57 , 641–643]. All three polymorphs have extensive networks of hydrogen bonds between the NH₃⁺ and COO⁻ groups of neighbouring molecules. In β‐gabapentin, there is an additional weak intramolecular hydrogen bond.