Theoretical studies on nitrogen rich energetic azoles

Different nitro azole isomers based on five membered heterocyclics were designed and investigated using computational techniques in order to find out the comprehensive relationships between structure and performances of these high nitrogen compounds. Electronic structure of the molecules have been calculated using density functional theory (DFT) and the heat of formation has been calculated using the isodesmic reaction approach at B3LYP/6-31G* level. All designed compounds show high positive heat of formation due to the high nitrogen content and energetic nitro groups. The crystal densities of these energetic azoles have been predicted with different force fields. All the energetic azoles show densities higher than 1.87 g/cm³. Detonation properties of energetic azoles are evaluated by using Kamlet-Jacobs equation based on the calculated densities and heat of formations. It is found that energetic azoles show detonation velocity about 9.0 km/s, and detonation pressure of 40GPa. Stability of the designed compounds has been predicted by evaluating the bond dissociation energy of the weakest C-NO₂ bond. The aromaticity using nucleus independent chemical shift (NICS) is also explored to predict the stability via delocalization of the π-electrons. Charge on the nitro group is used to assess the impact sensitivity in the present study. Overall, the study implies that all energetic azoles are found to be stable and expected to be the novel candidates of high energy density materials (HEDMs).     

Computational studies on polynitropurines as potential high energy density materials

As part of a search for high energy density materials (HEDMs), a series of purine derivatives with nitro groups were designed computationally. The relationship between the structures and the performances of these polynitropurines was studied. Density functional theory (DFT) at the B3LYP/6-311G** level was employed to evaluate the heats of formation (HOFs) of the polynitropurines by designing an isodesmic reaction method. Results indicated that the HOFs were influenced by the number and positions of substituent groups. Detonation properties were evaluated using the Kamlet–Jacobs equations, based on the theoretical densities and heats of formation of the polynitropurines. The relative stabilities of the polynitropurines were studied via the pyrolysis mechanism and the UB3LYP/6-311G** method. Homolysis of the ring–NO₂ bond is predicted to be the initial step in the thermal decomposition of these purine derivatives. Considering their detonation properties and relative stabilities, the tetranitropurine (D1) derivatives may be regarded as potential candidates for practical HEDCs. These results may provide useful information for further investigations.     

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.     

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 of -NH2 and -NO2 substituted dipyridines

In this work, the experimental synthesized bipyridines 3,3'-Dinitro-2,2'-bipyridine (DNBPy), 3,3'-Dinitro-2,2'-bipyridine-1,1'-dioxide (DNBPyO), and (3-Nitro-2-pyridyl)(5-nitro-2-pyridyl) amine (NPyA), and a set of designed dipyridines that have similar frameworks but different linkages and substituents with NPyA were studied theoretically at the B3LYP/6-31G* level of density functional theory. The gas-phase heats of formation were predicted based on the isodesmic reactions and the condensed-phase heats of formation and heats of sublimation were estimated in the framework of the Politzer approach. The crystal densities have been computed from molecular packing. Results show that this method gives a good estimation of density in comparison with the available experimental data for DNBPy, DNBPyO, and NPyA. The predicted detonation velocities and pressures indicate that the performance of dipyridines linked with -O-, -NH-, or -CH(2)- bridges have not been improved compared with that of the directly linked dipyridines, but all derivatives have better detonation properties than DNBPy, DNBPyO, and NPyA because of the presence of more nitro groups. An analysis of the bond dissociation energies (BDEs) or the impact sensitivity (h (50)) suggests that introduction of different bridges but not substituents has little influence on thermal stability. The calculated h (50) may be more reliable than BDE for predicting stability. Four bridged bipyridines have quite good detonation performance and low sensitivity.     

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).     

Comparative theoretical studies of energetic pyrazole-pyridine derivatives

The pyrazole-pyridine derivatives were optimized to obtain their molecular geometries and electronic structures at the DFT-B3LYP/6-31G(d,p) and DFT-B3P86/6-31G(d,p) levels. Molecular mechanics (MM) calculations were performed for the title compounds. Heats of formation (HOFs) were predicted through designed isodesmic reactions. Detonation performance was evaluated by using the Kamlet-Jacobs equations based on the calculated densities and heats of formation. The thermal stability of the title compounds was investigated via the bond dissociation energies (BDEs). The simulation results reveal that the compound with one pyrazole ring that is fully nitro-substituted performs similarly to the famous explosive HMX, and the compound with two pyrazole rings that are fully nitro-substituted outperforms HMX. According to the quantitative standard of energetics and stability as high energy density materials (HEDMs), the compound with two pyrazole rings that are fully nitro-substituted essentially satisfies this requirement.     

Computational study of energetic nitrogen-rich derivatives of 1,4-bis(1-azo-2,4-dinitrobenzene)-iminotetrazole

The heats of formation (HOFs), electronic structure, energetic properties, and thermal stabilities for a series of 1,4-bis(1-azo-2,4-dinitrobenzene)-iminotetrazole derivatives with different substituents and substitution positions and numbers of nitrogen atoms in the nitrobenzene rings were studied using the DFT-B3LYP method. All the substituted compounds have higher HOFs than their parent compounds. As the number of nitrogen atoms in the nitrobenzene ring increases, the HOFs of the derivatives with the same substituent rise gradually. Replacing carbon atoms in the nitrobenzene with nitrogen atoms to form N–N bonds is very helpful in improving their HOFs. Most of the substituted compounds have higher HOMO–LUMO gaps than the corresponding unsubstituted compounds. Substitution of the –NO₂, –NF₂, or –ONO₂ group and an increase in the number of nitrogen atoms in the nitrobenzene rings are useful for enhancing their detonation performance. The substituents’ substitution is not favorable for improving thermal stability. Considering detonation performance and thermal stability, five compounds may be considered potential candidates for high energy density compounds (HEDCs).     

Design and selection of nitrogen-rich bridged di-1,3,5-triazine derivatives with high energy and reduced sensitivity

The heats of formation (HOFs), electronic structures, energetic properties, and thermal stabilities of a series of energetic bridged di-1,3,5-triazine derivatives with different substituents and linkages were studied using density functional theory. It was found that the groups -N(3) and -N=N- are effective structural units for improving the HOF values of the di-1,3,5-triazine derivatives. The effects of the substituents on the HOMO-LUMO gap combine with those of the bridge groups. The calculated detonation velocities and detonation pressures indicate that substituting the -ONO(2), -NF(2), or -N=N- group is very useful for enhancing the detonation performance of these derivatives. Analysis of the bond dissociation energies for several relatively weak bonds suggests that most of the derivatives have good thermal stability. On the whole, the -NH(2), -N(3), -NH-, and -CH=CH- groups are effective structural units for increasing the thermal stabilities of the derivatives. Based on detonation performance and thermal stability, nine of the compounds can be considered potential candidates for high energy density materials with reduced sensitivity.     

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.     

Computational design and structure–property relationship studies on heptazines

This study aimed to design novel nitrogen-rich heptazine derivatives as high energy density materials (HEDM) by exploiting systematic structure–property relationships. Molecular structures with diverse energetic substituents at varying positions in the basic heptazine ring were designed. Density functional techniques were used for prediction of gas phase heat of formation by employing an isodesmic approach, while crystal density was assessed by packing calculations. The results reveal that nitro derivatives of heptazine possess a high heat of formation and further enhancement was achieved by the substitution of nitro heterocycles. The crystal packing density of the designed compounds varied from 1.8 to 2 g cm⁻³, and hence, of all the designed molecules, nitro derivatives of heptazine exhibit better energetic performance characteristics in terms of detonation velocity and pressure. The calculated band gap of the designed molecules was analyzed to establish sensitivity correlations, and the results reveal that, in general, amino derivatives possess better insensitivity characteristics. The overall performance of the designed compounds was moderate, and such compounds may find potential applications in gas generators and smoke-free pyrotechnic fuels as they are rich in nitrogen content.     

Density functional study on the derivatives of purine

The derivatives of purine are designed through substituting the hydrogen atoms on it for nitro and amino functional groups. Geometries and frequency are analyzed at the B3LYP/6-31 G** level of density functional theory(DFT). Heats of formation (HOF), bond dissociation energy(BDE) and detonation parameters (detonation velocity and detonation pressure) are obtained in detail at the same level. It is found that the BDE values of all derivatives are over 120KJ·mol(-1), and have high positive heats of formation. These derivatives possess excellent detonation properties, for B1, B2, and C, the detonation velocity are 9.58, 9.57,and 9.90 km·s(-1), and the detonation pressure are 43.40,46.05, and 46.37 Gpa, respectively, the detonation performances are better than cyclotrimethylenetrinitramine (RDX)and cyclotetramethylenetetranitramine (HMX). Hence, the derivations of purine may be promising well-behaved high energy density materials.     

Looking for high energy density compounds among polynitraminecubanes

Based on fully optimized geometric structures at DFT-B3LYP/6-311G** level, we calculated electronic structures, heats of formation, strain energies, bond dissociation energies and detonation performance (detonation velocity and detonation pressure) for a series of polynitraminecubanes. Our results have shown that energy gaps of cubane derivatives are much higher than that of triaminotrinitrobenzene (TATB), which means that cubane derivatives may be more sensitive than TATB. Polynitraminecubanes have high and positive heats of formation, and a good linear relationship between heats of formation and nitramine group numbers was presented. As the number of nitramine groups in the molecule increases, the enthalpies of combustion values are increasingly negative, but the specific enthalpy of combustion values decreases. It is found that all cubane derivatives have high strain energies, which are affected by the number and position of nitramine group. The calculated bond dissociation energies of C-NHNO₂ and C-C bond show that the C-C bond should be the trigger bond in the pyrolysis process. It is found that detonation velocity (D), detonation pressure (P) and molecule density (ρ) have good linear relationship with substituented group numbers. Heptanitraminecubane and octanitraminecubane have good detonation performance over 1,3,5,7-tetranitro-1,3,5,7-tetraazacyclooctane (HMX), and they can be regarded as potential candidates of high energy density compounds (HEDCs). The results have not only shown that these compounds may be used as HEDCs, but also provide some useful information for further investigation.     

Theoretical design of novel energetic salts derived from bicyclo-HMX

We designed three novel cage energetic anions by introducing ionic bridges containing NΘ, N(OΘ) and N(NΘNO₂) into cis-2,4,6,8-tetranitro-1H,5H-2,4,6,8- tetraazabicyclo[3.3.0] octane (bicyclo-HMX or BCMHX). The properties of 21 energetic salts, based on cage anions and ammonium-based cations, were studied by density functional theory (DFT) and volume-based thermodynamics (VBT) calculations. Compared to the parent nonionic BCHMX, most title salts have lower predicted impact sensitivities, higher predicted densities, larger predicted heats of formation (HOFs) and better predicted detonation properties. In particular, 11 energetic salts not only exhibit excellent predicted energetic properties, superior to 2,4,6,8,10,12-hexanitrohexaazaisowurtzitane (CL-20), but also have lower predicted sensitivity than CL-20. The best salt had a predicted detonation velocity of 10.06 km s^?1, a predicted detonation pressure of 48.54 GPa and a predicted sensitivity (h₅₀) of 23.99 cm. By introducing ionic bridges into highly nitrated rings, or modifying the original bridge with ionic bridges, some highly nitrated cage compounds with both excellent performance and low sensitivity can be developed strategically.

Open image in new window

Graphical abstract Heats of detonation, detonation velocities, and detonation pressures of salts derived from bicyclo-HMX

     

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 studies of energetic nitrogen-rich ionic salts composed of substituted 5-nitroiminotetrazolate anions and various cations

We have performed density functional theory and volume-based thermodynamics calculations to study the effects of different combinations of energetic anions and cations on the crystal densities, heats of formation, energetic properties, and thermodynamics of formation for a series of 5-nitroiminotetrazolate-based ionic salts. The results show that the substitution of the -NO₂, -NF₂, -N₃, or -C(NO₂)₃ group is helpful for increasing the densities of the salts. Incorporating every substituent (-NH₂, -NO₂, -NF₂, -N₃, or -C(NO₂)₃) into the salt is favorable for improving its HOF and detonation performance. Incorporating the cation B1, B2, B10, or B11 into the salts is helpful for improving its detonation properties. Increasing negative charge for the 5-nitroiminotetrazolate-based salts is unfavorable for enhancing the density and detonation performance, but is helpful for increasing the HOFs. Many salts present comparable detonation performance with commonly used explosives RDX or HMX. Among them, 21 salts have near or better properties than HMX. The thermodynamics of formation of the salts show that the majority of the 5-nitroiminotetrazolate salts with the cation B1, B3, B9, B10, B12 could be synthesized by the proposed reactions.     

Theoretical study on novel nitrogen-rich energetic compounds of bis(amino)-azobis(azoles) with tetrazene unit

Six stereoisomers of 5,5'-bis(amino)-1,1'-azobis(tetrazoles) and 30 other structures, including all possible bis(amino)-azobis(azoles) with an N-N=N-N unit, were designed. The molecular geometries were fully optimized at the DFT-B3LYP level with the 6-31++g (d, p) basis set. From the absence of any imaginary frequency in the infrared vibration frequency spectrum, it is predicted that all these studied structures may exist in stable forms. The results of the total energies of the stereoisomers of 5,5'-bis(amino)-1,1'-azobis(tetrazoles) indicate that the two symmetric trans-form structures are more likely to exist than the other four. The pyrolysis process, chemical stability and molecular electrostatic potential were studied via the investigation of their electronic structure. Heats of formation (HOFs) were calculated using the atomization energy method based on the results of the harmonic vibration frequencies, and a linear relationship was found between the HOF and nitrogen chain or nitrogen content. Densities of the title compounds were predicted with the Monte Carlo method. Finally, according to the results of the calculated HOFs and densities, the explosive parameters of these compounds were calculated using the Kamlet-Jacobs formula. 5,5'-Bis(amino)-1,1'-azobis(tetrazoles) and its isomer 5,5'-bis(amino)-2,2'-azobis(tetrazoles) may have potential for use as energetic compounds.     

The search for new powerful energetic transition metal complexes based on 3,3′-dinitro-5,5′-bis-1,2,4-triazole-1,1′-diolate anion: a DFT study

In this study, employing a new high oxygen balance energetic 3,3′-dinitro-5,5′-bis-1,2,4-triazole-1,1′-diolate anion (DNBTDO) as the bidentate ligand, NH₃ and NH₂NO₂ as short energetic ligands, and Cu/Ni as the metal atoms, two series of novel energetic metal complexes were computationally designed. Their structures and properties were studied by density functional theory, electrostatic potential data, and molecular mechanics methods. The results showed that the designed metal complexes have high detonation performance and acceptable sensitivity: Cu/Ni(DNBTDO)(NH₂NO₂)₂ (A3/B3) have better detonation properties and lower sensitivity than the most powerful CHNO explosive hexanitrohexaazaisowurtzitane, Cu/Ni(DNBTDO)(NH₃)(NH₂NO₂) (A2/B2) have comparable energetic performance and sensitivity with 1,3,5,7-tetranitro-1,3,5,7-tetrazocane, Ni(DNBTDO)(NH₃)₂ (B1) has comparative energy level and sensitivity with 1,3,5-trinitro-1,3,5-triazinane. These five energetic metal complexes may be attractive to energetic materials researchers. Besides, both the energetic ligands and metal atoms could have a great influence on the structures, heats of formation, detonation properties, and stability of energetic metal complexes, and the effects are coupled with each other. This study may be helpful in the search for and development of new improved energetic materials.     

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.     

Theoretical study on benzoheterocycle based energetic materials,effect of heterocyclic-fused,conjugation, hydrogen bond,and substitutional group on the detonation performance

In this paper, four series of benzoheterocycle based energetic materials (EMs) have been designed to plan out a strategy to improve the density and safety of EMs, such as combining the insensitive group with aminobenzene ring and the high energetic nitramine explosives, benzo-heterocycle mother ring, designing multi-nitrogen heterocycles with a conjugated system containing N-N and C-N high energy bonds, and hydrogen bonding. Their optimized structure and detonation properties were first calculated and discussed using DFT methods. After calculation, these designed explosives all showed good detonation from 7352 m/s to 8788 m/s. Among them, the compounds with six nitro groups, 1c, 2c, 3c, and 4c, exhibit better performance and rather poor impact sensitivity. However, we found that the compounds with five nitro groups and one amino group have a limited performance reduction and a rapid stability improvement. These four compounds, 1b, 2b, 3b, and 4b, have good detonation performance and better stability. Moreover, the synthesis routes for these four compounds were also designed. The precursor 4–0 and mononitro product 4–1 were successfully synthesized. Their 1H NMR, single crystal, and elemental analysis were also done to verify the structures.