Molecular dynamics calculation on structures,stabilities, mechanical properties,and energy density of CL-20/FOX-7 cocrystal explosives

In this article, different CL-20/FOX-7 cocrystal models were established by the substitution method based on the molar ratios of CL-20:FOX-7. The structures and comprehensive properties, including mechanical properties, stabilities, and energy density, of different cocrystal models were obtained and compared with each other. The main aim was to estimate the influence of molar ratios on properties of cocrystal explosives. The molecular dynamics (MD) simulation results show that the cocrystal model with molar ratio 1:1 has the best mechanical properties and highest binding energy, so the CL-20/FOX-7 cocrystal model is more likely to form in 1:1 M ratio. The detonation parameters show that the cocrystal explosive exhibited preferable energy density and excellent detonation performance. In a word, the 1:1 cocrystal model has the best comprehensive properties, is very promising, and worth more theoretical investigations and experimental tests. This paper gives some original theories to better understand the cocrystal mechanism and provides some helpful guidance and useful instructions to help design CL-20 cocrystal explosives.     

Theoretical investigation of the structures and properties of CL-20/DNB cocrystal and associated PBXs by molecular dynamics simulation

In this work, a CL-20/DNB cocrystal explosive model was established and six different kinds of fluoropolymers, i.e., PVDF, PCTFE, F₂₃₁₁, F₂₃₁₂, F₂₃₁₃ and F₂₃₁₄ were added into the (1 0 0), (0 1 0), (0 0 1) crystal orientations to obtain the corresponding polymer bonded explosives (PBXs). The influence of fluoropolymers on PBX properties (energetic property, stability and mechanical properties) was investigated and evaluated using molecular dynamics (MD) methods. The results reveal a decrease in engineering moduli, an increase in Cauchy pressure (i.e., rigidity and stiffness is lessened), and an increase in plastic properties and ductility, thus indicating that the fluoropolymers have a beneficial influence on the mechanical properties of PBXs. Of all the PBXs models tested, the mechanical properties of CL-20/DNB/F₂₃₁₁ were the best. Binding energies show that CL-20/DNB/F₂₃₁₁ has the highest intermolecular interaction energy and best compatibility and stability. Therefore, F₂₃₁₁ is the most suitable fluoropolymer for PBXs. The mechanical properties and binding energies of the three crystal orientations vary in the order (0 1 0)?>?(0 0 1)?>?(1 0 0), i.e., the mechanical properties of the (0 1 0) crystal orientation are best, and this is the most stable crystal orientation. Detonation performance results show that the density and detonation parameters of PBXs are lower than those of the pure CL-20 and CL-20/DNB cocrystal explosive. The power and energetic performance of PBXs are thus weakened; however, these PBXs still have excellent detonation performance and are very promising. The results and conclusions provide some helpful guidance and novel instructions for the design and manufacture of PBXs.     

Theoretical investigation of the effects of the molar ratio and solvent on the formation of the pyrazole–nitroamine cocrystal explosive 3,4-dinitropyrazole (DNP)/2,4,6,8,10,12-hexanitrohexaazaisowurtzitane (CL-20)

The effects of the molar ratio, temperature, and solvent on the formation of the cocrystal explosive DNP/CL-20 were investigated using molecular dynamics (MD) simulation. The cocrystal structure was predicted through Monte Carlo (MC) simulation and using first-principles methods. The results showed that the DNP/CL-20 cocrystal might be more stable in the molar ratio 1:1 near to 318 K, and the most probable cocrystal crystallizes in the triclinic crystal system with the space group P\( \overline{1} \). Cocrystallization was more likely to occur in methanol and ethanol at 308 K as a result of solvent effects. The optimized structure and the reduced density gradient (RDG) of the DNP/CL-20 complex confirmed that the main driving forces for cocrystallization were a series of hydrogen bonds and van der Waals forces. Analyses of the trigger bonds, the charges on the nitro groups, the electrostatic surface potential (ESP), and the free space per molecule in the cocrystal lattice were carried out to further explore their influences on the sensitivity of CL-20. The results indicated that the DNP/CL-20 complex tended to be more stable and insensitive than pure CL-20. Moreover, an investigation of the detonation performance of the DNP/CL-20 cocrystal indicated that it possesses high power.

Open image in new window

Graphical abstract DNP/CL-20 cocrystal models with different molar ratios were investigated at different temperatures using molecular dynamics (MD) simulation methods. Binding energies and mechanical properties were probed to determine the stability and performance of each cocrystal model. Solvated DNP/CL-20 models were established by adding solvent molecules to the cocrystal surface. The binding energies of the models in various solvents were calculated in order to identify the most suitable solvent and temperature for preparing the cocrystal explosive DNP/CL-20

     

Theoretical insights into the structures and mechanical properties of HMX/NQ cocrystal explosives and their complexes,and the influence of molecular ratios on their bonding energies

Molecular dynamics (MD) methods were employed to study the binding energies and mechanical properties of selected crystal planes of 1,3,5,7-tetranitro-1,3,5,7-tetrazacyclooctane (HMX)/nitroguanidine (NQ) cocrystals at different molecular molar ratios. The densities and detonation velocities of the cocrystals at different molar ratios were estimated. The intermolecular interaction and bond dissociation energy (BDE) of the N–NO₂ bond in the HMX:NQ (1:1) complex were calculated using the B3LYP, MP2(full) and M06-2X methods with the 6-311++G(d,p) and 6-311++G(2df,2p) basis sets. The results indicated that the HMX/NQ cocrystal prefers cocrystalizing in a 1:1 molar ratio, and the cocrystallization is dominated by the (0 2 0) and (1 0 0) facets. The K, G, and E values of the ratio of 1:1 are smaller than those of the other ratios, and the 1:1 cocrystal has the best ductility. The N–NO₂ bond becomes stronger upon the formation of the intermolecular H-bonding interaction and the sensitivity of HMX decreases in the cocrystal. This sensitivity change in the HMX/NQ cocrystal originates not only from the formation of the intermolecular interaction but also from the increment of the BDE of N–NO₂ bond in comparison with isolated HMX. The HMX/NQ (1:1) cocrystal exhibits good detonation performance. Reduced density gradient (RDG) reveals the nature of cocrystallization. Analysis of the surface electrostatic potential further confirmed that the sensitivity decreases in complex (or cocrystal) in comparison with that in isolated HMX.

Open image in new window

Graphical Abstract Binding energies and mechanical properties of HMX/NQ cocrystals in different molecular molar ratios were studied using molecular dynamics methods. The origin of the sensitivity change in the HMX/NQ cocrystal originates from formation of intermolecular interactions and the bond dissociation energy increment of the N–NO₂ bond

     

Theoretical insights into the stabilities,detonation performance,and electrostatic potentials of cocrystals containing α- or β-HMX and TATB,FOX-7, NTO,or DMF in various molar ratios

A molecular dynamics method was employed to study the binding energies associated with the cocrystallization (at selected crystal planes) of either 1,3,5-triamino-2,4,6-trinitro-benzene (TATB), 1,1-diamino-2,2-dinitroethylene, 3-nitro-1,2,4-triazol-5-one (TATB, FOX-7, and NTO, respectively, all of which are explosives), or N,N-dimethylformamide (DMF, a nonenergetic solvent) in various molar ratios with 1,3,5,7-tetranitro-1,3,5,7-tetrazacyclooctane in its α and β conformations (α-HMX and β-HMX, respectively). The results showed that the cocrystals with low molar ratios (2:1, 1:1, 1:2, and 1:3) were the most stable. The binding energies of HMX/NTO and HMX/DMF were larger than those of HMX/TATB and HMX/FOX-7. According to the calculated stabilities, HMX prefers to adopt its α form in HMX/TATB and its β form in HMX/NTO, whereas the two forms coexist in HMX/FOX-7. For HMX/TATB, HMX/NTO, and α-HMX/FOX-7, increasing the proportion of the cocrystal component with the highest detonation heat (HMX in the first two cases, FOX-7 in the latter) increases the detonation heat, velocity, and pressure of the cocrystal. However, increasing the proportion of the component with the highest detonation heat in β-HMX/FOX-7 and γ-CL-20/FOX-7 increases the detonation heat of the cocrystal but decreases its detonation velocity. An investigation of the surface electrostatic potential revealed how the sensitivity changes upon cocrystal formation.

Open image in new window

Graphical Abstract Surface electrostatic potential of HMX/TATB

     

Easy methods to study the smart energetic TNT/CL-20 co-crystal

2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (CL-20) is a high-energy nitramine explosive with high mechanical sensitivity. 2,4,6-trinitrotoluene (TNT) is insensitive but by no means a high performance explosive. To reveal the significant importance and smart-material functionality of the energetic-energetic co-crystals, the stability, mechanical and explosive properties TNT/CL-20 co-crystal, TNT crystal and CL-20 crystal were studied. Non-hydrogen bonded non-covalent interactions govern the structures of energetic-energetic co-crystals. However, it is very difficult to accurately calculate the non-covalent intermolecular interaction energies. In this paper, the local conformation and the intricate non-covalent interactions were effectively mapped and analyzed from the electron density (ρ) and its derivatives. The results show that the two components TNT and CL-20 are connected mainly by nitro–aromatic interactions, and nitro–nitro interactions. The steric interactions in TNT/CL-20 could not be confronted with the attractive interactions. Moreover, the scatter graph of TNT crystal reveals the reason why TNT is brittle. The detailed electrostatic potential analysis predicted that the detonation velocities (D) and impact sensitivity for the compounds both increase in the sequence of CL-20 > TNT/CL-20 co-crystal > TNT. Additionally, TNT/CL-20 co-crystal has better malleability than its pure components. This demonstrates the capacity and the feasibility of realizing explosive smart materials by co-crystallization, even if strong hydrogen bonding schemes are generally lacking in energetic materials.

Growth morphology of CL-20/HMX cocrystal explosive: insights from solvent behavior under different temperatures

A 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (CL-20) /1,3,5,7-tetranitro-1,3,5,7-tetrazacyclooctane (HMX)-isopropanol (IPA) interfacial model was constructed to investigate the effect of temperature on cocrystal morphology. A constant volume and temperature molecular dynamics (NVT-MD) simulation was performed on the interfacial model at various temperatures (295–355 K, 20 K intervals). The surface electrostatic potential (ESP) of the CL-20/HMX cocrystal structure and IPA molecule were studied by the B3LYP method at 6–311++G (d, p) level. The surface energies, polarities, adsorption energy, mass density distribution, radial distribution function (RDF), mean square displacement (MSD) and relative changes of attachment energy were analyzed. The results show that polarities of (1 0 0) and (0 1 1) cocrystal surfaces may be more negative and affected by IPA solvent. The adsorption energy per area indicates that growth of the (1 0–2) face in IPA conditions may be more limited, while the (1 0 0) face tends to grow more freely. MSD and diffusion coefficient (D) analyses demonstrated that IPA molecules gather more easily on the cocrystal surface at lower temperatures, and hence have a larger effect on the growth of cocrystal faces. RDF analysis shows that, with the increasing of temperature, the strength of hydrogen bond interactions between cocrystal and solvent becomes stronger, being highest at 335 K for the (1 0 0) and (0 1 1) interfacial models. Results of relative changes of modified attachment energy show that (1 0 0) and (0 1 1) faces tends to be larger than other faces. Moreover, the predicted morphologies at 295 and 355 K are consistent with experimental values, proving that the CL-20/HMX-IPA interfacial model is a reasonable one for this study.

Open image in new window

Graphical Abstract Construction of 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (CL-20) /1,3,5,7-tetranitro-1,3,5,7-tetrazacyclooctane (HMX)-isopropanol (IPA) interfacial model, analysis, and morphology prediction of cocrystal.

     

Theoretical studies on the structures and detonation properties of nitramine explosives containing benzene ring

The nitramine compounds containing benzene ring were optimized to obtain their molecular geometries and electronic structures at DFT-B3LYP/6-31+G(d) level. The theoretical molecular density (ρ), heat of formation (HOF), energy gap (ΔE(LUMO-HOMO)), charge on the nitro group (-Q(NO2)), detonation velocity (D) and detonation pressure (P), estimated using Kamlet-Jacobs equations, showed that the detonation properties of these compounds were excellent. It is found that there are good linear relationships between density, heat of formation, detonation velocity, detonation pressure and the number of nitro group. The simulation results reveal that molecule G performs similarly to famous explosive HMX, and molecule H outperforms HMX. According to the quantitative standard of energetics as an HEDC (high energy density compound), molecule H essentially satisfies this requirement. These results provide basic information for molecular design of novel high energetic density compounds.     

Theoretical studies on the thermodynamic properties, densities, detonation properties, and pyrolysis mechanisms of trinitromethyl-substituted aminotetrazole compounds

Trinitromethyl-substituted aminotetrazoles with –NH₂, –NO₂, –N₃, and –NHC(NO₂)₃ groups were investigated at the B3LYP/6-31G(d) level of density functional theory. Their sublimation enthalpies, thermodynamic properties, and heats of formation were calculated. The thermodynamic properties of these compounds increase with temperature as well as with the number of nitro groups attached to the tetrazole ring. In addition, the detonation velocities and detonation pressures of these compounds were successfully predicted using the Kamlet–Jacobs equations. It was found that these compounds exhibit good detonation properties, and that compound G (D = 9.2 km/s, P = 38.8 GPa) has the most powerful detonation properties, which are similar to those of the well-known explosive HMX (1,3,5,7-tetranitro-1,3,5,7-tetrazocine). Finally, the electronic structures and bond dissociation energies of these compounds were calculated. The BDEs of their C–NO₂ bonds were found to range from 101.9 to 125.8 kJ/mol^-1. All of these results should provide useful fundamental information for the design of novel HEDMs.     

Molecular design of new nitramine explosives: 1,3,5,7-tetranitro-8-(nitromethyl)-4-imidazolino[4,5-b]4-imidazolino-[4,5-e] pyridine and its N-oxide

Two new nitramine compounds containing pyridine, 1,3,5,7-tetranitro-8-(nitromethyl) -4-imidazolino[4,5-b]4-imidazolino-[4,5-e]pyridine and its N-oxide 1,3,5,7-tetranitro-8- (nitromethyl)-4-imidazolino[4,5-b]4-imidazolino-[4,5-e]pyridine-4-ol were proposed. Density functional theory (DFT) has been employed to study the molecular geometries, electronic structures, infrared spectra, and thermodynamic properties at the B3LYP/6-31G* level. Their detonation performances evaluated using the Kamlet-Jacobs equations with the calculated densities and heats of formation are superior to those of HMX. The predicted densities of them were ca. 2 g*cm^-3, detonation velocities were over 9 km*s^-1, and detonation pressures were about 40 GPa, showing that they may be potential candidates of high energy density materials (HEDMs). The natural bond orbital analysis indicated that N-NO₂ bond is the trigger bond during thermolysis process. The stability of the title compounds is slightly lower than that of 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12- hexaazaisowurtzitane (CL-20). The results of this study may provide basic information for the molecular design of new HEDMs.     

Density functional calculations for a high energy density compound of formula C6H6?n (NO2) n

A series of polynitroprismanes, C(6)H(6-n )(NO(2))(n) (n?=?1-6) intended for use as high energy density compounds (HEDCs) were designed computationally. Their electronic structures, heats of formation, interactions between nitro groups, specific enthalpies of combustion, bond dissociation energies, and explosive performances (detonation velocities and detonation pressures) were calculated using density functional theory (DFT) with the 6-311 G** basis set. The results showed that all of the polynitroprismanes had high positive heats of formation that increased with the number of substitutions for the prismane derivatives, while the specific enthalpy of combustion decreased as the number of nitro groups increased. In addition, the range of enthalpy of combustion reducing is getting smaller. Interactions between ortho (vicinal) groups deviate from the group additivity rule and decrease as the number of nitro groups increases. In terms of thermodynamic stability, all of the polynitroprismanes had higher bond dissociation energies (BDEs) than RDX and HMX. Detonation velocities and detonation pressures were estimated using modified Kamlet-Jacobs equations based on the heat of detonation (Q) and the theoretical density of the molecule (ρ). It was found that ρ, D, and P are strongly linearly related to the number of nitro groups. Taking both their energetic properties and thermal stabilities into account, pentanitroprismane and hexanitroprismane are potential candidate HEDCs.     

Dinitroamino benzene derivatives: a class new potential high energy density compounds

Dinitroamino benzene derivatives are designed and studied in detail with quantum chemistry method. The molecular theory density, heats of formation, bond dissociation energies, impact sensitive and detonation performance are investigated at DFT-B3LYP/6-311G** level. The results of detonation performance indicated most of the compounds have better detonation velocity and pressure than RDX and HMX. The N-N bond can be regard as the trigger bond in explosive reaction, and the bond dissociation energies of trigger bond are almost not affected by the position and number of substituent group. The impact sensitive are calculated by two different theory methods. It is found that the compounds, which can become candidates of high energy materials, have smaller H₅₀ values than RDX and HMX. It is hoped that this work can provide some basis information for further theory and experiment studies of benzene derivatives.     

Theoretical investigation on the heats of formation and detonation performance in polydinitroaminocubanes

A series of polydinitroaminocubanes have been designed computationally. We calculated the heats of formation, the detonation velocity (D) and detonation pressure (P) of the title compounds by density function theory (DFT) with 6-311?G** basis set. The relationship between the heats of formation and the molecular structures is discussed. The result shows that all cubane derivatives have high and positive heats of formation, which increase with increasing number of dinitroamino groups. The detonation performances of the title compound were estimated by Kamlet-Jacobs equation, and the result indicated that most cubane derivatives have good detonation performance over RDX (hexahydro-1,3,5-trinitro-1,3,5-trizine) and HMX (1,3,5,7-tetranitro-1,3,5,7-tetraazacyclooctane). In addition, we also found that the heat of detonation (Q) is another very important impact in increasing detonation performance except density. The relative stabilities of the title compound are discussed in the terms of the calculated heats of formation, and the energy gaps between the frontier orbitals. The results have not only shown that these compounds may be used as high energy density compounds (HEDCs), but also provide some useful information for further investigation.     

Is 1-nitro-1-triazene a high energy density material?

An azo bridge (–N?=?N–) can not only desensitize explosives but also dramatically increase their heats of formation and explosive properties. Amino and nitro are two important high energy density functional groups. Here, we present calculations on 1-nitro-1-triazene (NH₂–N?=?N–NO₂). Thermal stability and detonation parameters were predicted theoretically at CCSD(T)/6-311G* level, based on the geometries optimized at MP2/6-311G* level. It was found that the p?→?π conjugation interaction and the intramolecular hydrogen bonding that exist in the system together increase the thermal stability of the molecule. Moreover, the detonation parameters were evaluated to be better than those of the famous HMX and RDX. Finally, the compound was demonstrated to be a high energy density material.     

Theoretical studies on the heats of formation, densities, and detonation properties of substituted s-tetrazine compounds

Substituted s-tetrazine compounds were designed and investigated in order to find comprehensive relationships between the structures and performances of high-nitrogen energetic compounds. Density functional theory (DFT) was used to predict the optimized geometries, electronic structures, heats of formation and densities, and the detonation properties were evaluated by using the VLW equation of state (EOS). Calculation results show that there are good linear relationships between heats of formation, densities, detonation properties and the number of N atom in all designed high-nitrogen compounds. Furthermore, several designed high-nitrogen compounds show good detonation velocities and pressures compared with octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocane (HMX), making them potential candidates for high-energy-density materials (HEDM).     

A first principle study of the structural,electronic, and temperature-dependent thermodynamic properties of graphene/MoS2 heterostructure

Journal of Molecular Modeling - Based on the unique advantages of terahertz (THz) spectrum on the detection of energetic cocrystals, the low-temperature dependent THz spectra of CL-20/TNT cocrystal...     

Quantum chemical studies on the aminopolynitropyrazoles

We have explored the geometric and electronic structures, band gap, thermodynamic properties, density, detonation velocity and detonation pressure of aminopolynitropyrazoles using the density functional theory (DFT) at the B3LYP/aug-cc-pVDZ level. The calculated detonation velocity and detonation pressure, stability and sensitivity of model compounds appear to be promising compared to the known explosives 3,4-dinitro-1 H-pyrazole (3,4-DNP), 3,5-dinitro-1 H-pyrazole (3,5-DNP), hexahydro-1,3,5-trinitro-1,3,5-triazinane (RDX) and octahydro-1,3,5,7-tetranitro-l,3,5,7-tetraazocane (HMX). The position of NH₂ group in the polynitropyrazoles presumably determines the structure, stability, sensitivity, density, detonation velocity and detonation pressure.     

Metal-doped graphene layers composed with boron nitride–graphene as an insulator: a nano-capacitor

A new family of energetic azacubane N-oxides were designed by introducing N-oxides into azacubanes and investigated by using density functional theory. Introducing the N-oxides into the azacubanes could improve their detonation performance significantly due to the increase of the OB and ρ but would also increase the sensitivity to some extent. These effects would be further enhanced as the numbers of N-oxides increase. Among all the designed azacubane N-oxides, D6-4 (1,3,5,7-tetraazacubane-1,3,5,7-tetraoxides) has higher detonation performance than one famous high explosive HMX (1,3,5,7-tetranitro-1,3,5,7-tetrazocane) and lower sensitivity than one very insensitive explosive TNT (1-methyl-2,4,6-trinitrobenzene), suggesting that its overall performance is outstanding and may be considered as the potential candidate of insensitive high explosives. The internal small cage C-N skeleton of D6-4 is surrounded by the external big cage hydrogen bonds and this special double cage structure may be an important reason why it has low sensitivity.     

Chemotaxis-mediated biodegradation of cyclic nitramine explosives RDX, HMX, and CL-20 by Clostridium sp. EDB2

Cyclic nitramine explosives, RDX, HMX, and CL-20 are hydrophobic pollutants with very little aqueous solubility. In sediment and soil environments, they are often attached to solid surfaces and/or trapped in pores and distribute heterogeneously in aqueous environments. For efficient bioremediation of these explosives, the microorganism(s) must access them by chemotaxis ability. In the present study, we isolated an obligate anaerobic bacterium Clostridium sp. strain EDB2 from a marine sediment. Strain EDB2, motile with numerous peritrichous flagella, demonstrated chemotactic response towards RDX, HMX, CL-20, and NO(2)(-). The three explosives were biotransformed by strain EDB2 via N-denitration with concomitant release of NO(2)(-). Biotransformation rates of RDX, HMX, and CL-20 by the resting cells of strain EDB2 were 1.8+/-0.2, 1.1+/-0.1, and 2.6+/-0.2nmol h(-1)mgwet biomass(-1) (mean+/-SD; n=3), respectively. We found that commonly seen RDX metabolites such as TNX, methylenedinitramine, and 4-nitro-2,4-diazabutanal neither produced NO(2)(-) during reaction with strain EDB2 nor they elicited chemotaxis response in strain EDB2. The above data suggested that NO(2)(-) released from explosives during their biotransformation might have elicited chemotaxis response in the bacterium. Biodegradation and chemotactic ability of strain EDB2 renders it useful in accelerating the bioremediation of explosives under in situ conditions.     

Comparative theoretical studies of differently bridged nitramino-substituted ditetrazole 2-<Emphasis Type="Italic">N</Emphasis>-oxides with high detonation performance and an oxygen balance of around zero

In this work, six (A–F) nitramino (–NHNO₂)-substituted ditetrazole 2-N-oxides with different bridging groups (–CH₂–, –CH₂–CH₂–, –NH–, –N=N–, and –NH–NH–) were designed. The six compounds were based on the parent compound tetrazole 2-N-oxide, which possesses a high oxygen balance and high density. The structure, heat of formation, density, detonation properties (detonation velocity D and detonation pressure P), and the sensitivity of each compound was investigated systematically via density functional theory, by studying the electrostatic potential, and using molecular mechanics. The results showed that compounds A–F all have outstanding energetic properties (D: 9.1–10.0 km/s; P: 38.0–46.7 GPa) and acceptable sensitivities (h ₅₀: 28–37 cm). The bridging group present was found to greatly affect the detonation performance of each ditetrazole 2-N-oxide, and the compound with the –NH–NH– bridging group yielded the best results. Indeed, this compound (F) was calculated to have comparable sensitivity to the famous and widely used high explosive 1,3,5,7-tetranitro-1,3,5,7-tetrazocane (HMX), but with values of D and P that were about 8.7% and 19.4% higher than those for HMX, respectively. The present study shows that tetrazole 2-N-oxide is a useful parent compound which could potentially be used in the design of new and improved high-energy compounds to replace existing energetic compounds such as HMX.