Molecular design of aminopolynitroazole-based high-energy materials

The density functional theory (DFT) was employed to calculate the energetic properties of several aminopolynitroazoles. The calculations were performed to study the effect of amino and nitro substituents on the heats of formation, densities, detonation performances, thermal stabilities, and sensitivity characteristics of azoles. DFT-B3LYP, DFT-B3PW91, and MP2 methods utilizing the basis sets 6-31 G* and 6-311 G (2df, 3p) were adopted to predict HOFs via designed isodesmic reactions. All of the designed aminopolynitroazoles had heats of formation of >220 kJ mol(-1). The crystal densities of the aminopolynitroazoles were predicted with the cvff force field. All of the energetic azoles had densities of >1.83 g/cm(3). The detonation velocities and pressures were evaluated using the Kamlet-Jacobs equations, utilizing the predicted densities and heats of formation. It was found that aminopolynitroazoles have a detonation velocity of about 9.1 km/s and detonation pressure of 36 GPa. The bond dissociation energies for the C-NO(2) and N-NO(2) bonds were analyzed to investigate the stabilities of the designed molecules. The charge on the nitro group was used to assess impact sensitivity in the present study. The results obtained imply that the designed molecules are stable and are expected to be candidates for high-energy materials (HEMs).     

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

     

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.     

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

Theoretical design of energetic nitrogen-rich derivatives of 1,7-diamino-1,7-dinitrimino-2,4,6-trinitro-2,4,6-triazaheptane

The heats of formation (HOFs), energetic properties, and thermal stability of a series of 1,7-diamino-1,7-dinitrimino-2,4,6-trinitro-2,4,6-triazaheptane derivatives with different substituents, different numbers of substituents, and different original chains are found by using the DFT-B3LYP method. The results show that -NO₂ or -NH₂ is an effective substituent for increasing the gas-phase HOFs of the title compounds, especially -NO₂ group. As the numbers of substitutents increase, their HOFs enhance obviously. Increasing the length of original chain is helpful for improving their HOFs. The substitution of -NO₂ is useful for enhancing their detonation performances and the effects of the length of original chains on detonation properties are coupled with those of the substituents. An analysis of the BDE of the weakest bonds indicates that the substitution of the -NH₂ groups and replacing the -NO₂ groups of N-NO₂ by the -NH₂ groups are favorable for improving their thermal stability, while the substitution of -NO₂ and increasing the length of original chain decrease their thermal stability. Considering the detonation performance and thermal stability, seven compounds may be considered as the potential candidates of high energy density compounds.     

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.     

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.     

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.     

Characterization of nitrogen-bridged 1,2,4,5-tetrazine-, furazan-, and 1H-tetrazole-based polyheterocyclic compounds: heats of formation, thermal stability, and detonation properties

The heats of formation (HOFs), thermal stability, and detonation properties for a series of nitrogen-bridged 1,2,4,5-tetrazine-, furazan-, and 1H-tetrazole-based polyheterocyclic compounds (3,6-bis(1H-1,2,3,4-tetrazole-5-ylamino)-1,2,4,5- tetrazine (TST), 3,6-bis(furazan-5-ylamino)-1,2,4,5-tetrazine (FSF), 3,4-bis(1,2,4,5- tetrazine-3-ylamino)-furazan (SFS), 3,4-bis(1H-1,2,3,4-tetrazole-5-ylamino)-furazan (TFT), 1,5-bis(1,2,4,5-tetrazine-3-ylamino)-1H-1,2,3,4-tetrazole (STS), and 1,5-bis(furazan-3-ylamino)-1H-1,2,3,4-tetrazole (FTF) derivatives) were systematically studied by using density functional theory. The results show that the -N(3) or -NHNH(2) group plays a very important role in increasing the HOF values of the derivatives. Among these series, the SFS derivatives have lower energy gaps, while the TFT derivatives have higher ones. Incorporation of the -NH(2) group into the FSF, SFS, STS, or FTF ring is favorable for enhancing its thermal stability, whereas the substitution of the -NHNH(2) group could increase the thermal stability of the TST, SFS, STS, or FTF ring. The calculated detonation properties indicate that the -NO(2) or -NF(2) is very helpful for enhancing the detonation performance for these derivatives. Considering the detonation performance and thermal stability, six derivatives may be regarded as promising candidates of high-energy density materials (HEDMs). These results provide basic information for the molecular design of novel HEDMs.     

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.     

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.     

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.     

New chiral oligopyridines-4,4'-bis(disaccharide)-functionalised 2,2'-bipyridines and 4'-(disaccharide)-functionalised 2,2':6',2'-terpyridines

A new class of oligopyridine ligands bearing disaccharides linked to the 4- and 4'-positions of a 2,2'-bipyridine or the 4'-position of a 2,2':6',2'-terpyridine metal-binding domain are described. Representative ligands with furanosylfuranose and pyranosylpyranose (cellobiose) substituents have been prepared.     

A theoretical study on 1,5-diazido-3-nitrazapentane (DANP) and 1,7-diazido-2,4,6-trinitrazaheptane (DATNH): molecular and crystal structures, thermodynamic and detonation properties, and pyrolysis mechanism

1,5-Diazido-3-nitrazapentane (DANP) and 1,7-diazido-2,4,6-trinitrazaheptane (DATNH) are two energetic plasticizers. To better understand them, a detailed theoretical investigation was carried out using density functional theory and molecular mechanics methods. The crystal structures, spectra, thermodynamic properties, heats of formation, detonation velocity, detonation pressure, specific impulse and thermal stability were estimated. Possible initiation steps of pyrolysis were discussed by considering the bond breaking of N–NO₂, C–N₃, and N–N₂ (via hydrogen transfer) for both compounds and the cyclization of the adjacent nitro and azido groups for DATNH. Results show that the rupture of N–NO₂ and N–N₂ (via hydrogen transfer) may happen simultaneously as the initial step of pyrolysis. Both crystals have P-1 symmetry as was observed experimentally. DANP has higher stability than DATNH, while DATNH has better detonation performance than DANP. In addition, DANP has a lower while DATNH has a higher specific impulse than RDX, which shows their prospects as propellant components.     

Looking for high energy density compounds among polynitraminepurines

A series of purine derivatives with nitramine groups are calculated by using density functional theory (DFT). The molecular theory density, heats of formation, bond dissociation energies and detonation performance are investigated at DFT-B3LYP/6-311G** level. The isodesmic reaction method is employed to calculate the HOFs of the energies obtained from electronic structure calculations. Results show that the position of nitramine groups can influence the values of HOFs. The bond dissociation energies and the impact sensitivity are analyzed to investigate the thermal stability of the purine derivatives. The calculated bond dissociation energies of ring-NHNO₂ and NH-NO₂ bond show that the NH-NO₂ bond should be the trigger bond in pyrolysis processes. The H₅₀ of most compounds are larger than that of CL-20 and RDX.     

A DFT study of tautomers of 3-amino-1-nitroso-4-nitrotriazol-5-one-2-oxide

We report herein the structure and explosive properties of the possible isomers of 3-amino-1-nitroso-4-nitrotriazol-5-one-2-oxide computed from the B3LYP/aug-cc-pVDZ level. The optimized structures, vibrational frequencies and thermodynamic values for triazol-5-one-N-oxides were obtained in the ground state. Several designed compounds have densities varying from 2.103 to 2.177 g/cm³. The detonation properties were evaluated by the Kamlet-Jacob equations based on the predicted density and the calculated heat of explosion. The detonation properties of triazol-5-one-N-oxides (D 9.87 to 10.11 km s^?1 and P 48.95 to 50.61 GPa) appear to be promising compared with those of 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (D 9.20 km s^?1, P 42.0 Gpa) and octanitrocubane (D 9.90 km s^?1, P 48.45 GPa). The substitution of secondary amino hydrogen of the triazole ring by amino group shows better impact sensitivity/or stability however the model compounds seem to be highly sensitive.     

Crystal structure, detonation performance, and thermal stability of a new polynitro cage compound: 2, 4, 6, 8, 10, 12, 13, 14, 15-nonanitro-2, 4, 6, 8, 10, 12, 13, 14, 15-nonaazaheptacyclo [5.5.1.13, 11. 15, 9] pentadecane

A new polynitro cage compound 2, 4, 6, 8, 10, 12, 13, 14, 15-nonanitro-2, 4, 6, 8, 10, 12, 13, 14, 15-nonaazaheptcyclo [5.5.1.1(3,11).1(5,9)] pentadecane (NNNAHP) was designed in the present work. Its molecular structure was optimized at the B3LYP/6-31 G(d,p) level of density functional theory (DFT) and crystal structure was predicted using the Compass and Dreiding force fields and refined by DFT GGA-RPBE method. The obtained crystal structure of NNNAHP belongs to the P-1 space group and the lattice parameters are a = 9.99 ?, b = 10.78 ?, c = 9.99 ?, α = 90.01°, β = 120.01°, γ = 90.00°, and Z = 2, respectively. Based on the optimized crystal structure, the band gap, density of state, thermodynamic properties, infrared spectrum, strain energy, detonation characteristics, and thermal stability were predicted. Calculation results show that NNNAHP has detonation properties close to those of CL-20 and is a high energy density compound with moderate stability.     

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.     

Conjugation in multi-tetrazole derivatives: a new design direction for energetic materials

Multi-tetrazole derivatives with conjugated structures were designed and investigated in this study. Using quantum chemistry methods, the crystal structures, electrostatic potentials (ESPs), multicenter bond orders, HOMO–LUMO energy gaps, and detonation properties of the derivatives were calculated. As expected, these molecules with conjugated structures showed low energies of their crystal structures, molecular layering in their crystals, high average ESPs, high multicenter bond order values, and enhanced detonation properties. The derivative 1,2-di(1H-tetrazol-5-yl)diazene (N2) was predicted to have the best density (1.87 g/cm³), detonation velocity (9006 m/s), and detonation pressure (36.8 GPa) of the designed molecules, while its total crystal energy was low, suggesting that it is relatively stable. Its sensitivity was also low, as the molecular stacking that occurs in its crystal allows external forces to be dissipated into movements of crystal layers. Finally, its multicenter bond order was high, indicating a highly conjugated structure.     

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