Theoretical studies on four-membered ring compounds with NF2, ONO2, N3, and NO2 groups

Density functional theory (DFT) method has been employed to study the geometric and electronic structures of a series of four-membered ring compounds at the B3LYP/6-311G** and the B3P86/6-311G** levels. In the isodesmic reactions designed for the computation of heats of formation (HOFs), 3,3-dimethyl-oxetane, azetidine, and cyclobutane were chosen as reference compounds. The HOFs for N(3) substituted derivations are larger than those of oxetane compounds with --ONO2 and/or --NF2 substituent groups. The HOFs for oxetane with --ONO2 and/or --NF2 substituent groups are negative, while the HOFs for N3 substituted derivations are positive. For azetidine compounds, the substituent groups within the azetidine ring affect the HOFs, which increase as the difluoroamino group being replaced by the nitro group. The magnitudes of intramolecular group interactions were predicted through the disproportionation energies. The strain energy (SE) for the title compounds has been calculated using homodesmotic reactions. For azetidine compounds, the NF2 group connecting N atom in the ring decrease the SE of title compounds. Thermal stability were evaluated via bond dissociation energies (BDE) at the UB3LYP/6-311G** level. For the oxetane compounds, the O--NO2 bond is easier to break than that of the ring C--C bond. For the azetidine and cyclobutane compounds, the homolyses of C--NX2 and/or N--NX2 (X = O, F) bonds are primary step for bond dissociation. Detonation properties of the title compounds were evaluated by using the Kamlet-Jacobs equation based on the calculated densities and HOFs. It is found that 1,1-dinitro-3,3-bis(difluoroamino)-cyclobutane, with predicted density of ca. 1.9 g/cm(3), detonation velocity (D) over 9 km/s, and detonation pressure (P) of 41 GPa that are lager than those of TNAZ, is expected to be a novel candidate of high energy density materials (HEDMs). The detonation data of nitro-BDFAA and TNCB are also close to the requirements for HEDMs.     

Theoretical Studies on Energetic Property and Stability of Pyrazine and Pyridine Derivatives

Density functional calculations at the B3LYP level with 6‐311G** and aug‐cc‐pVDZ basis sets were performed to predict the heats of formation (HOFs) for two pyrazine derivatives and eight pyridine derivatives. In the isodesmic reactions designed for the computation of heats of formation (HOFs), pyrazine and pyridine were chosen as reference compounds. The N‐oxidations for the ring nitrogen of pyrazine and pyridine derivatives decrease the HOF values when N‐oxide oxygen is neighboring with amino groups, but increase when it neighbors with nitro groups. Thermal stability was evaluated via bond dissociation energies (BDE) at the UB3LYP/6‐311G** level. As a whole, the homolysis of C–NO₂ bonds is the main step for bond dissociation of the title compounds. The BDE values of title compounds are influenced by intramolecular hydrogen bonds. The hydrogen bond effects associated with the length of the H···O bonds were analyzed by the electron density at the critical points and natural bond orbital.     

Theoretical Study of Heterocyclic Nitro Compounds for Use as Novel Explosives and Propellants

The heats of formation (HOFs) of heterocyclic nitro compounds were obtained by using a density functional theory B3LYP method with 6‐31G* and 6‐311+G** basis sets. The isodesmic reactions designed for the evaluation of HOFs keep most of the basic ring structures of the title compounds and thus ensure the credibility of the results. The values of HOFs are 567.90, 874.29 and 975.83 kJ/mol at the B3LYP/6‐31G* level for hexanitrohexazaadamantane ( A ), nonanitrononaza‐tetracyclo[7.3.1.1^3,7.1^5,11] pentadecane ( B ) and tetranitrotetrazacubane ( C ) respectively. The predicted detonation velocities of the title compounds are larger than, and detonation pressures are much larger than that of the widely used 1,3,5,7‐tetranitro‐1,3,5,7‐tetraazacyclooctane (HMX). The dissociation energy for the weakest C‐N bonds in the cage skeleton of the title compounds are 137‐144 kJ/mol at the B3LYP/6‐31G* level.     

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.     

Theoretical studies of 1,2,4,5-tetrazine-based energetic nitrogen-rich compounds

Density functional theory calculations were performed to find the relationships between the structures and performance of a series of 1,2,4,5-tetrazine-based energetic derivatives. The isodesmic reaction method was employed to estimate the heats of formation (HOFs). The result shows that the azo or azoxy group is one of the most energetic functional groups known and its substitution can drastically increase HOFs of a molecule. The detonation properties were also evaluated by the Kamlet–Jacobs equations based on the theoretical densities and HOFs. Results show that NO₂ group is an effective substituent for enhancing the detonation performance. There exist better correlations between OB and detonation velocities and OB and detonation pressures. The energy gaps between the HOMO and LUMO of the studied compounds are also investigated, and from the data we estimated the relative thermal stability ordering of the title compounds.     

Density functional calculation on a high energy density compound having the formula C2OH4?n (NO2) n

A series of ethylene oxide derivations, C₂OH_4?Cn (NO2) _n (x?=?1?C4), has been designed computationally. We calculated the heats of formation (HOFs), bond dissociation energy (BDE), and explosive performances (detonation velocity and detonation pressure) of the title compounds by using density functional theory with 6-311G** basis set. The results show that most of ethylene oxide derivations have positive HOFs values except I. All the calculated BDE and the bond dissociation energies without zero-point energy corrections (BDE⁰) are larger than 200?kJ?mol^?1, which indicate that all the ethylene oxide derivations have good thermal stabilities. The explosive performances of most of ethylene oxide derivations would rank up with cyclotrimethylenetrinitramine (RDX). The results have not only shown that these compounds may be used as high energy density compounds, but also provide some useful information for further syntheses.     

Molecular design of prismane-based potential energetic materials with high detonation performance and low impact sensitivity

To develop new energetic materials, the eleven nitroester substitution derivatives of prismane were investigated at the B3LYP/6-311G** level of density functional theory (DFT). The gas phase heats of formation were calculated by isodesmic reactions and the solid-state heats of formation were obtained by the Politzer approach using the heats of sublimation for the designed compounds. The detonation velocities and pressures of all molecules were calculated by Kamlet–Jacobs equations based on molecular density and heat of detonation. The results show that the nitroester group in prismane is helpful for enhancing molecular detonation properties and power index. Among all molecules, 1,2,3,4-tetrnitroesterprismane has excellent detonation properties (detonation pressure = 40.05 GPa, detonation velocity = 9.28 km/s) and large power index value. The molecular stabilities were evaluated by calculating bond dissociation energies and characteristic heights (H₅₀). The results indicate that the bond dissociation energies of all molecules are above 80 kJ/mol, and all molecules have a larger H₅₀ value than hexanitrohexaazaisowurtzitane (CL-20, 12 cm). The obtained structure–property relationships may provide basic information for the molecular design of novel high-energy materials.     

Theoretical investigation on the structures,densities, and detonation properties of polynitrotetraazaoctahydroanthracenes

The polynitrotetraazaoctahydroanthracenes were optimized to obtain their molecular geometries and electronic structures at density functional theory–B3LYP/6‐31+G(d) level. Detonation velocities (D) and detonation pressures (P) were estimated for this nitramine compounds using Kamlet‐Jacobs equations, based on the theoretical densities (ρ) and heats of formation. It is found that there are good linear relationships between volume, density, detonation velocity, detonation pressure and the number of nitro group. Thermal stability of the compounds was investigated by calculating the bond dissociation energies and energy gap (ΔE_LUMO–HOMO). The simulation results reveal that molecule H performs similarly to famous explosive RDX. These results provide basic information for molecular design of novel high energetic density compounds. © 2011 Wiley Periodicals, Inc.     

Computational study on energetic properties of nitro derivatives of furan substituted azoles

Density functional theory has been used to investigate geometries, heats of formation (HOFs), C-NO₂ bond dissociation energies (BDEs), and relative energetic properties of nitro derivatives of azole substituted furan. HOFs for a series of molecules were calculated by using density functional theory (DFT) and Møller–Plesset (MP2) methods. The density is predicted using crystal packing calculations; all the designed compounds show density above 1.71 g/cm³. The calculated detonation velocities and detonation pressures indicate that the nitro group is very helpful for enhancing the detonation performance for the designed compounds. Thermal stabilities have been evaluated from the bond dissociation energies. Charge on the nitro group was used to assess the impact sensitivity in this study. According to the results of the calculations, tri- and tetra-nitro substituted derivatives reveal high performance with better thermal stability.     

Theoretical investigation of a high density cage compound 1,3,5,7,9,11-hexanitrotetradecahydro-1H-1,3,4,5,7,7b,9,11,12a,12b1,12b2,13-dodecaaza-4,8,12-(epimethanetriyl)cyclohepta[l]cyclopenta[def] phenanthrene

The B3LYP/6-31G** method was used to investigate IR and Raman spectra, heat of formation, and thermodynamic properties of a new designed polynitro cage compound 1,3,5,7,9,11-hexanitrotetradecahydro-1H-1,3,4,5,7,7b,9,11,12a,12b¹,12b²,13-dodecaaza-4,8,12-(epimethanetriyl)cyclohepta[l]cyclopenta[def]phenanthrene. The detonation and pressure were evaluated using the Kamlet–Jacobs equations based on the theoretical density and HOFs. The bond dissociation energies and bond orders for the weakest bonds were analyzed to investigate the thermal stability of the title compound. The results show that N₈–NO₂ bond is predicted to be the trigger bond during pyrolysis. There exists an essentially linear relationship between the WBIs of N–NO₂ bonds and the charges – $ Q_{{{\text{NO}}_{ 2} }} $ on the nitro groups. The crystal structure obtained by molecular mechanics belongs to the P2₁ space group, with lattice parameters Z = 2, a = 11.4658 Å, b = 15.2442 Å, c = 10.2451 Å, ρ = 2.07 g cm^?3. The designed compound has high thermal stability and good detonation properties and is a promising high-energy density compound.     

高能量密度材料3,3′-偶氮-1,2,4,5-四嗪衍生物的分子设计

运用密度泛函理论(DFT)方法,计算系列3,3′-偶氮-1,2,4,5-四嗪衍生物的生成热.结果显示:—N3取代基在增加3,3′-偶氮-1,2,4,5-四嗪衍生物的生成热方面起了非常重要的作用.通过分析标题化合物的最弱键离解能发现:—NH2或—N3取代基非常有利于增加衍生物的热稳定性.计算的爆速(D)和爆压(p)数值表明:—NO2或—NF2取代基有利于提高3,3′-偶氮-1,2,4,5-四嗪衍生物的爆轰性能.综合爆轰性能和热稳定性的计算结果,3种3,3′-偶氮-1,2,4,5-四嗪衍生物可以作为潜在的品优高能量密度材料(HEDM)候选物.     

Theoretic design of 1,2,3,4-tetrazine-1,3-dioxide-based high-energy density compounds with oxygen balance close to zero

Density functional theory method was used to study the heats of formation (HOFs), electronic structure, energetic properties, and thermal stability for a series of 1,2,3,4-tetrazine-1,3-dioxide derivatives with different substituents and bridge groups. It is found that the groups –NO₂, –C(NO₂)₃, and –N=N– play a very important role in increasing the HOFs of the derivatives. The effects of the substituents on the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) energy levels and HOMO–LUMO gaps are coupled to those of different substituents and bridges. The calculated detonation velocities and pressures indicate that the group –NO₂, –NF₂, –ONO₂, –C(NO₂)₃, or –NH– 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 indicates that incorporating the groups –NO₂, –NF₂, –ONO₂, –C(NO₂)₃, and –N=N– into parent ring decreases their thermal stability. Considering the detonation performance and thermal stability, 18 compounds may be considered as the target compounds holding the greatest potential for synthesis and use as high-energy density compounds. Among them, the oxygen balances of four compounds are equal to zero. These results provide basic information for the molecular design of the novel high-energy compounds.     

DFT Studies on Detonation Properties,Pyrolysis Mechanism,and Stability of the Nitro and Hydroxyl Derivatives of Benzene

Ninety‐one nitro and hydroxyl derivatives of benzene were studied at the B3LYP/6‐31G?? level of density functional theory. Detonation properties were calculated using the Kamlet‐Jacobs equation. Three candidates (pentanitrophenol, pentanitrobenzene, and hexanitrobenzene) were recommended as potential high energy density compounds for their perfect detonation performances and reasonable stability. The pyrolysis mechanism was studied by analyzing the bond dissociation energy (BDE) and the activation energy (E_a) of hydrogen transfer (H–T) reaction for those with adjacent nitro and hydroxyl groups. The results show that E_a is much lower than BDEs of all bonds, so when there are adjacent nitro and hydroxyl groups in a molecule, the stability of the compound will decrease and the pyrolysis will be initiated by the H–T process. Otherwise, the pyrolysis will start from the breaking of the weakest C–NO₂ bond, and only under such condition, the Mulliken population or BDE of the C–NO₂ bond can be used to assess the relative stability of the compound.     

Computational investigation of the heat of formation, detonation properties of furazan-based energetic materials

The heats of formation (HOFs) for a series of furazan-based energetic materials were calculated by density functional theory. The isodesmic reaction method was employed to estimate the HOFs. The result shows that the introductions of azo and azoxy groups can increase the HOF, but the introduction of azo group can increase the more HOF, when compared with azoxy group. The detonation velocities and detonation pressures of the furazan-based energetic materials are further evaluated at B3LYP/6-31G* level. Dioxoazotetrafurazan and azoxytetrafurazan may be regarded as the potential candidates of high-energy density materials because of good detonation performance. In addition, there are good linear correlations between OB and detonation velocities, and OB and detonation pressures. The energy gaps between the HOMO and LUMO of the studied compounds are also investigated. These results provide basic information for the molecular design of novel high-energy density materials.     

Computational studies on tetrazole derivatives as potential high energy materials

The tetrazole is an important functionality of the most of energetic materials due to 80% nitrogen content, stability, and high enthalpy of formation. The present structure–property relationship study focuses on the optimized geometries of tetrazole derivatives obtained from density functional theory (DFT) calculations at B3LYP/6-31G* levels. The heat of formation (HOF) of tetrazole derivatives have been calculated by designing the appropriate isodesmic reactions. The increase in nitro groups on azole rings shows the remarkable increase in HOF. Density has been predicted by using CVFF force field. Increase in the nitro group increases the density. Detonation properties of the designed compounds were evaluated by using the Kamlet–Jacobs equation based on predicted densities and HOFs. Designed tetrazole derivatives show detonation velocity (D) over 8 km/s and detonation pressure (P) of about 32 GPa. Thermal stability was evaluated via bond dissociation energies (BDE) of the weakest C–NO₂ bond at B3LYP/6-31G* level. Charge on the nitro group has been used to assess the sensitivity correlation. Overall, the study implies that designed compounds of this series are found to be stable and expected to be the novel candidates of high energy materials (HEMs).     

A theoretical investigation on the structures,densities, detonation properties,and pyrolysis mechanism of the nitro derivatives of phenols

The nitro derivatives of phenols are optimized to obtain their molecular geometries and electronic structures at the DFT‐B3LYP/6‐31G* level. Detonation properties are evaluated using the modified Kamlet–Jacobs equations based on the calculated densities and heats of formation. It is found that there are good linear relationships between density, detonation velocity, detonation pressure, and the number of nitro and hydroxy groups. Thermal stability and pyrolysis mechanism of the title compounds are investigated by calculating the bond dissociation energies (BDEs) at the unrestricted B3LYP/6‐31G* level. The activation energies of H‐transfer reaction is smaller than the BDEs of all bonds and this illustrates that the pyrolysis of the title compounds may be started from breaking O? H bond followed by the isomerization reaction of H transfer. Moreover, the C? NO₂ bond with the smaller bond overlap population and the smaller BDE will also overlap may be before homolysis. According to the quantitative standard of energetics and stability as a high‐energy density compound, pentanitrophenol essentially satisfies this requirement. In addition, we have discussed the effect of the nitro and hydroxy groups on the static electronic structural parameters and the kinetic parameter. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Computational DFT studies on a series of toluene derivatives as potential high energy density compounds

Based on the full optimized molecular geometric structures at B3LYP/6-31G**, B3LYP/6-31+G**, B3P86/6-31G**, and B3P86/6-31+G** levels, the densities (ρ), detonation velocities (D), and pressures (P) for a series of toluene derivatives, as well as their thermal stabilities, were investigated to look for high energy density compounds (HEDCs). The heats of formation (HOFs) are also calculated via designed isodesmic reactions. The calculations on the bond dissociation energies (BDEs) indicate that the BDEs of the initial scission step are between 48 and 59 kcal/mol, and pentanitrotoluene is the most reactive compound, while 2,4,6-trinitrotoluene is the least reactive compound for toluene derivatives studied. A good linear relationship between BDE/E and impact sensitivity is also obtained. The condensed phase HOFs are calculated for the title compounds. These results would provide basic information for the further studies of HEDCs. The detonation data of pentanitrotoluene show that it meets the requirement for HEDCs.     

Theoretical investigation on density,detonation properties,and pyrolysis mechanism of nitro derivatives of benzene and aminobenzenes

Nitro derivatives of benzene and aminobenzenes are optimized at the DFT‐B3LYP/6‐31G* level. The heat of formation (ΔH_f) and crystal theoretical density (ρ) are estimated to evaluate the detonation properties using the modified Kamlet–Jacobs equations. Thermal stability and the pyrolysis mechanism of the title compounds are investigated by calculating the bond dissociation energies (BDE) at the unrestricted B3LYP/6‐31G* level. The kinetic parameter and the static electronic structural parameters can be used to predict the stability and the relative magnitude of the impact sensitivity of homologues. According to the quantitative standard of the energy and the stability as an HEDC, the title compounds having more than four nitro groups satisfy this requirement. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Nitro Derivatives of Parabanic Acid as Potential Explosives: A Theoretical Study

This work deals with certain parabanic acid (PA) derivatives because they possess great calculated density (>1.8 g · cm^–3) and high content of nitrogen (26 %). Computed ballistic properties of eight different parabanic acid derivatives are presented. The structures were optimized at the B3LYP/6‐31G(d, p) level. The calculated data for PA are found to be compatible with the experimental X‐ray data. The detonation performance analyses were done using empirical Kamlet‐Jacobs equations. Additionally, detonation products were assigned and power index were calculated. All the compounds considered are powerful candidates for high energy materials.     

Theoretical studies of furoxan-based energetic nitrogen-rich compounds

Density functional theory calculations were performed to study the effects of different substituents and nitrogen-containing heterocyclic rings on the heats of formation (HOFs), energetic properties, and thermal stability for a series of furoxan derivatives. The isodesmic reaction method was employed to calculate the HOFs of the derivatives using total energies obtained from electronic structure calculations. The detonation velocities and pressures were evaluated by using the semiempirical Kamlet–Jacobs equations, based on the theoretical densities and HOFs. The bond dissociation energies and bond orders for the weakest bonds were analyzed to investigate the thermal stability of the furoxan derivatives. These results provide basic information for the molecular design of novel high energy density materials.