Chinese Journal of Chemical Physics  2016, Vol. 29 Issue (3): 349-357

#### The article information

Mei Zheng, Xiao-hong Li, Hong-ling Cui, Rui-zhou Zhang

Comparative Theoretical Studies on Several Energetic Substituted Dioxin-imidazole Derivatives

Chinese Journal of Chemical Physics, 2016, 29(3): 349-357

http://dx.doi.org/10.1063/1674-0068/29/cjcp1506118

### Article history

Accepted on August 06, 2015
Comparative Theoretical Studies on Several Energetic Substituted Dioxin-imidazole Derivatives
Mei Zheng, Xiao-hong Li , Hong-ling Cui, Rui-zhou Zhang
Dated: Received on June 06, 2015; Accepted on August 06, 2015
College of Physics and Engineering, Henan University of Science and Technology, Luoyang 471003, China
Author: Xiao-hong Li, E-mail:lorna639@126.com
Abstract: The molecular structures, infrared spectra, heats of formation (HOFs), detonation properties, chemical and thermal stabilities of several tetrahydro-[1,4]dioxino[2,3-d:5,6-d'] diimidazole derivatives with different substituents were studied using DFT-B3LYP method. The properties of the compounds with different groups such as -NO2, -NH2, -N3, and -ONO2 were further compared. The -NO2 and -ONO2 groups are effective substituents for increasing the densities of the compounds, while the substitution of -N3 group can produce the largest HOF. The compound with -NO2 group has the same detonation properties as 1,3,5,7-tetranitro-1,3,5,7-tetraazacyclooctane, while the compound with -ONO2 group has lower detonation properties than those of hexahydro-1,3,5-trinitro-1,3,5-triazine. The nature bond orbital analysis reveals that the relatively weak bonds in the molecules are the bonds between substituent groups and the molecular skeletons as well as C-O bonds in the dioxin rings. The electron withdrawing groups (-NO2, -N3, and -ONO2) have inductive effects on the linkages between the groups and molecular skeletons. In addition, researches show that the electronegativities of the groups are related with the stabilities of the compounds. Considering detonation performance and thermal stability, the 1,5-dinitro-2,6-bis(trinitromethyl)-3a,4a,7a,8a-tetrahydro-[1,4]dioxino-[2,3-d:5,6-d'] diimidazole satisfies the requirements of high energy density materials.
Key words: Density functional theory    Condensed density    Heats of formation    Detonation property    Thermal stability
Ⅰ. INTRODUCTION

The use of imidazole derivatives as high energy density materials (HEDMs) to satisfy both civil and military applications was extensive during the past decades [1, 2, 3, 4]. Imidazole rings are relatively stable and rich in C and N atoms,so they are considered as HEDMs when linked with functional groups [5]. Various nitroimidazole derivatives were synthesized and have been investigated experimentally and theoretically. Results show that many of nitroimidazole derivatives meet excellent properties as HEDMs,such as cis-syn-cis-2,6-dioxo-1,3,4,5,7,8-hexanitrodecahydro-1H,5H-diimidazo[4,5-b:4',5'-e] pyr-azine (CL-20). It is polycyclic and is rich in nitro groups,which results in high density (2.07 g/cm3),and good detonation properties [5, 6, 7].

The synthesis of a nitrogen-rich imidazole poly-tri-heterocyclic compound,1,5-dinitro-2,6-bis(trinitromethyl)-3a,4a,7a,8a-tetrahydro-[1, 4]dioxino-[2,3-d:5,6-d'] diimidazole (DNTNDI) was reported by Wu et al. [8]. DNTNDI is constituted by a 1,4-dioxin ring connecting two methyl-imidazole rings,and all the H atoms on the two methyl-imidazole rings are substituted by nitro groups except the H atoms on the dioxin ring. Experiments show that DNTNDI has an excellent density of 1.929 g/cm3 and its oxygen balance is zero,which imply that the moles of gas production and the heat of detonation may reach maximum.

Generally,nitrate ester group (-ONO2) is an organic nitric acid group that can contain enormous explosive force and is often used to modify the mechanical properties [9]. The introduction of an amino group (-NH2) is one of the simplest means to improve the thermal stability of an energetic material [10]. In addition,azido group (-N3) has a high positive heat of formation and can produce more gases in the combustion products,thereby increasing the work capacity of propellant [11]. Considering the importance of the three groups,we wonder what changes the properties of the new energetic materials if the -NO2 groups are substituted by -NH2,-ONO2,and -N3 groups,respectively.

Heat of formation (HOF) is crucial thermodynamic quantity and is frequently taken to be indicative of the "energy content" of a HEDM [12, 13]. Detonation velocity (D) and pressure (P) are the most important targets of scaling the detonation characteristics of energetic materials. The bond dissociation energy (BDE) of the weakest bond for an explosive molecule is a key factor in many decomposition processes [14],and has been expected to play an important role in the initiation of detonation. Lots of calculations indicate that C-NO2 and N-NO2 bonds are the weakest bonds in energetic ring molecules and are the first rupture bond during decomposition process. The BDEs are also related to the stabilities and impact/shock sensitivities of HEDMs [15, 16, 17, 18]. Therefore,it is necessary to understand the properties and the internal mechanism of DNTNDI by using the theoretical investigations.

In this work,the molecular structures,infrared spectroscopy,HOF,thermal stabilities,and detonation properties are calculated for DNTNDI and the derivatives with -NH2,-ONO2,and -N3 groups by using density functional theory B3LYP method. The related properties of these compounds are predicted and their molecular structures are listed in Fig. 1.

 FIG. 1 Molecular frameworks of DNTNDI and derivatives 1, 2, 3.
Ⅱ. COMPUTATIONAL METHOD AND DETAILS

In this work,the theoretical data for DNTNDI and its derivatives compounds 1,2,3 (Fig. 1) were calculated by using density functional theory (DFT) B3LYP method with 6-311++G** basis set. The DFT method,particularly the B3LYP [19, 20, 21, 22] method was testified as a credible approach which not only produces reliable geometries,energies and electron population but also is efficient [23]. All the calculations are performed by Gaussian 03 program package [24].

Infrared Radiation (IR) spectrum is a measurable characteristic of a material which can effectively identify groups and/or ingredients. And it is directly related to the thermodynamic properties of compounds. Thus,the IR spectra for DNTNDI and compounds 1,2,3 were calculated based on optimized molecular structures. Besides,the previous investigations revealed that the calculated vibrational frequencies conform to the experimental frequencies with an appropriate scaling factor [25, 26, 27, 28].

The isodesmic reactions were designed to evaluate HOFs of DNTNDI and compounds 1,2,3. Isodesmic reaction method has been widely and successfully used in predicting HOFs [29, 30, 31, 32]. The isodesmic reactions remain the number of each kind of formal bond during being performed. For the title compounds,the HOFs were derived from the isodesmic reactions shown in Fig. 2.

 FIG. 2 Isodemic reactions of compound 1, 2, 3.

To calculate the detonation properties of DNTNDI and compounds 1,2,3,solid-phase HOFs (△Hf,solid) are required. The △Hf,solid can be evaluated by the employment of gas-phase HOFs (△Hf,gas) and heats of sublimation (△Hsub),as Eq(1):

 $\Delta {{H}_{\text{f},\text{solid}}}=\Delta {{H}_{\text{f},\text{gas}}}-\Delta {{H}_{\text{sub}}}$ (1)

Politzer et al. [33, 34] pointed out that the heats of sublimation of energetic materials containing CHNO atoms are related to their molecular surface areas and electrostatic interaction index. They presented the following relationship:

 $\Delta H_\mathrm{sub} = aA^2 + b(\nu \sigma _\mathrm{tot}^2)^{1/2} + c$ (2)

where A is the 0.001 electrons/bohr3 isosurface area of electronic density of the molecule,ν describes the degree of balance between positive and negative potential on the isosurface,and σtot2 is a measure of variability of the electrostatic potential on the molecular surface. The coefficients a,b and c were determined by Rice et al.,a=2.670×10-4 kcal/mol,b=1.650 kcal/mol,and c=2.966 kcal/mol [35]. Furthermore,Politzer et al. [33, 34, 35] found that the solid density of energetic materials containing CHNO atoms can be corrected by applying parameter νσtot2 for the theoretical density:

 $\rho = \beta _1 \frac{M}{V} + \beta _2 \nu \sigma _\mathrm{tot}^2 + \beta _3$ (3)

where M is the molecular weight in g/mol,V is the volume inside a contour of 0.001 electrons/bohr3 electron density for the molecule. The volume was computed using a Monte Carlo integration that set 2000 point/bohr3 for adequate accuracy. The coefficients β1,β2,and β3 are 0.9183,0.0028,and 0.0443,respectively.

The detonation velocity D and pressure P can be estimated by the semi-empirical Kamlet-Jacobs equations [36]:

 $D=1.01{{(N{{{\bar{M}}}^{1/2}}{{Q}^{1/2}})}^{1/2}}(1+1.30\rho )$ (4)
 $P=1.558{{\rho }^{2}}N{{{\bar{M}}}^{1/2}}{{Q}^{1/2}}$ (5)

where N is the moles of gas phase products from per gram of the explosive, $\overline{M}$ is the average molecular mass of these gases,Q is the standard enthalpy change of the detonation from per gram explosive,ρ is the loaded density of explosives.

The frontier molecular orbitals (FMO) and the average binding energies (ABEs) for title compounds were calculated using B3LYP/6-311++G** method in order to evaluate the chemical stabilities. ABE is defined as the quotient that the molecular binding energy divides the number of atoms in a molecule. The ABE can be calculated using the equation [37]:

 $\mathrm{ABE} = \frac{1}{n}{\sum\limits_{i = 1}^n {E(\mathrm{atom}_i)} - E(\mathrm{molecule})}$ (6)

where E(atomi) is the energy of the ith atom,E(molecule) is the total energy of the molecule,n is the number of atoms in the molecule.

In order to understand the stability of each bond in the title compounds,the natural bond orbital (NBO) [38] calculation was performed by using B3LYP/6-311++G** method. Based on the NBO analysis,pyrolysis mechanism and thermal stability of the title compounds were investigated via BDE. The BDE is the difference between the energies of parent molecule and radicals which are from the fission of the molecule [39, 40, 41]. The BDE can be calculated by using the following equation [42]:

 $\mathrm{BDE}(\mathrm{R}_1 \mathrm{R}_2) = E(\mathrm{R}_1) + E(\mathrm{R}_2) - E(\mathrm{R}_1 \mathrm{R}_2)$ (7)

where R1 and R2 are the radicals produced by the breaking of R1-R2 bond of R1R2 compound.

Ⅲ. RESULTS AND DISCUSSION A. Molecular geometry

The molecular geometries of the title compounds were optimized by B3LYP/6-311++G** method. For DNTNDI and compounds 2,3,the imidazole rings and the oxygen bridges are generally planar,while the oxygen atoms on dioxin ring are slightly deviated from the plane. Compound 1 looks like a curving boat with imidazole rings and oxygen bridges being not on a plane.

Selected bond lengths of DNTNDI and compounds 1,2,3 are listed in Table Ⅰ. For comparison,the bond lengths of the normal single bonds (experimental and calculated bond length) are also shown in Table Ⅰ. Obviously,the calculated bond lengths are very close to the experimental normal single bond lengths. This shows that our calculational method is reliable.

TABLE Ⅰ Selected bond lengths d (Å ) in DNTNDI, and compounds 1, 2, 3.

It is noted that the C-NO2 bond lengths (including C7-N12,13,14 and C31-N34,35,36 single bonds) and N-NO2 bond lengths in DNTNDI are much longer than normal single C-N bond and N-N bond,respectively. Besides,the N-NH2 bond lengths in compound 1 and the N-ONO2 bond lengths in compound 3 are shorter than the bond lengths of N-N and N-O normal single bonds,respectively. The bond lengths of the C-NH2 bonds in compound 1 and the C-ONO2 bonds in compound 3 are slightly shorter than those of C-N and C-O normal bonds,respectively.

B. Infrared spectra and vibrations

The infrared spectra of DNTNDI and compounds 1,2,3 were predicted using B3LYP/6-311++G** method. Figure 3 presents the calculated infrared spectra.

 FIG. 3 The calculated infrared spectra of the title compounds. (a) DNTNDI, (b) compound 1, (c) compound 2, and (d) compound 3.

The infrared vibration modes of these compounds correspond to the peaks. The asymmetric stretching vibrations of O-C-O and C-C-O in the dioxin ring of the title compounds are associated with the peaks of 1156,1129,1164,and 1156 cm-1 for DNINDI,and compounds 1,2,3,respectively. There is no split for these peaks,which implies that the vibration coupling effects of the molecular skeletons are not obvious.

Besides,other evident peaks are all associated with the -NO2,-NH2,-N3,and -ONO2 groups. It is worth noting that,most of the strongest peaks of these compounds correspond to the asymmetric stretching vibration of the groups (-NO2,-N3,and -ONO2),except compound 1. Compound 1 has the strongest peak at 877 cm-1,which is caused by the wagging vibration of -NH2 groups.

The investigations reveal that most of the strong peaks are associated with C-O single bonds in dioxin rings or the bonds between the substituent groups and the imidazole rings. Thereby,the vibrations of these parts in the molecules are easier to be enhanced when the compounds are heated.

C. Heats of formation

HOF is the change of enthalpy from the formation of 1 mole of the compound from its constituent elements,with all substances in their standard states at 1 atmosphere. So in general,the more HOF a compound contains,the more energy is released during the thermal decomposition. From Hess’s law,HOF can be used to estimate the amount of energy release in a detonation. Thereby,HOF is a very important factor which determines the detonation properties.

Isodesmic reaction is used to predict the HOFs of DNTNDI,and compounds 1,2,3. For this method,precise HOF data of reference compounds led to accurate prediction [44]. The energies and experimental HOFs of the reference compounds are listed in Table Ⅱ.

TABLE Ⅱ Calculated total energies (E0), zero-point energies (ZPE), thermal corrections (HT), and the heats of formation (HOF) of the reference compounds.

Table Ⅲ summarizes the total energies,ZPE,thermal corrections,HOFs (in gas and solid phase) of compounds DNTNDI,1,2,3. The results reveal that the DNTNDI and compound 2 exhibit positive HOFs and compound 2 with -N3 groups has the largest HOF among all these compounds. The energy of groups can be reflected by the energy of the corresponding methyl compounds. The comparison of experimental HOFs [46] of CH3N3 (298.3 kJ/mol),CH3NO2 (-80.8 kJ/mol),CH3NH2 (-22.5 kJ/mol),and CH3ONO2 (-122 kJ/mol) shows that the -N3 group has the largest energy among the four groups. Similarly,the -ONO2 group has the smallest energy because of the smallest HOF of CH3ONO2.

TABLE Ⅲ Calculated total energies (E0), zero-point energies (ZPE), thermal corrections (HT), and HOFs of DNTNDI and compounds 1, 2, 3.
D. Detonation propertie

D and P are key indicators to describe detonation properties for explosives. To estimate the D and P for explosive only containing CHON elements,$K-J$ equations can be applied [36]. Obviously,the crystal density ρ is the most crucial parameter which directly affects D and P. The crystal density of 1.883 g/cm3 for DNTNDI was obtained and is slightly smaller than experimental density 1.929 g/cm3 [8]. The reason may be that theoretical calculation only involves the single molecule,while the experiment involves the more molecules and considers the interactions among them.

In order to investigate the effect of the different groups on the properties,the molecular surface area A,the electrostatic potential parameters νσtot2,the oxygen balance OB,the condensed density ρ,and the predicted detonation velocity D and pressure P for the title compounds are summarized and listed in Table Ⅳ. The densities and detonation properties of hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) and 1,3,5,7-tetranitro-1,3,5,7-tetraazacyclooctane (HMX) are also listed in Table Ⅳ for comparison.

TABLE Ⅳ Predicted crystal densities ρ, detonation properties D, the electrostatic potential parameters νσ2tot , the oxygen balance OB, the molecular surface area A, and pressure P of the title compounds.

The predicted D and P of DNTNDI are 9.15 km/s and 38.14 GPa,respectively,which are close to those of HMX and all larger than those of RDX. The oxygen balance of DNTNDI is zero,which implies that the compound can combust completely,produce the maximum gas products and release the maximum heat. The D (8.15 km/s) and P (30.30 GPa) of compound 3 are smaller than those of RDX. Nevertheless,for compounds 1 and 2,the detonation properties are unsatisfactory. Compounds 1 and 2 have less oxygen and the C atoms cannot be combusted completely,so both the gas phase products and the heat release are reduced. It is observed that imidazole ring doesn't contain O atoms,so the groups containing oxygen atoms,such as -NO2 and -ONO2,are necessary to be involved in order to increase the oxygen balance for imidazole derivatives. In addition,it is noted from Table Ⅳ that ρ and △Hf,gas both play the important role in determining D and P.

E. Frontier molecular orbital and average binding energy analysis

FMO are the electron distributions of the highest occupied energy level (HOMO) and the lowest un-occupied energy level (LUMO),which are associated with the electric and optical properties,as well as the chemical reactivities [51]. The FMO of the title compounds were computed using B3LYP/6-311++G** method and are shown in Fig. 4.

 FIG. 4 The calculated FMO and corresponding energy levels of the title compounds based on B3LYP/6-311++G** level.

The FMOs of all these compounds except compound 1 are localized on the external functional groups,such as -NO2,-N3,and -ONO2 groups. Besides,the FMO of compound 1 are localized on the C-N bonds in the imidazole rings. No apparent delocalization effects of electron arise for all the FMOs. The values of the energy gaps between the HOMO and LUMO are 4.49,5.74,4.94,and 4.23 eV for DNINDI and compound 1,2,3,respectively.

The ABEs for DNINDI and compound 1,2,3 were also calculated and are shown in Fig. 5. The electronegativities of -NO2,-NH2,-N3,and -ONO2 groups are from Ref.[52]. The ABEs,HOMO-LUMO gap,and the reciprocal of the electronegativities for corresponding groups are shown in Fig. 4 for comparison. From Fig. 5,the ABEs,HOMO-LUMO gaps,and the reciprocal of the electronegativities for the title compounds follow the similar trends. Thus,for the dioxin-imidazole derivatives,the electronegativities of the substituents directly affect the FMO gaps and the ABEs. A larger electronegativity for the substituent leads to a lower HOMO-LUMO gap and a smaller ABE. The FMO gaps and ABEs reveal the higher chemical stability for compound 1 and the larger chemical reactivity for compound 3.

 FIG. 5 The calculated ABEs and HOMO-LUMO gaps for DNTNDI and compounds 1, 2, 3, as well as reciprocal of the electronegativities for the corresponding substituents 52.
F. Thermal stability

BDE is often used to describe the strength of bonds,discovers the weakest bond and investigates the pyrolysis mechanism and thermal stability. The NBO analyses of the title compounds were performed to select the weakest bonds. The NBO analysis not only provides the charge transfer interactions between donor NBO and acceptor NBO,but also exhibits the WBI,which gives an indication to the stability of a bond. A smaller WBI value indicates a weaker connection for bonds of same type,while a bigger WBI value implies a stronger linkage.

The bonds of the same type are selected and their WBIs are listed in Table Ⅴ for comparison. In the Table Ⅴ,the bonds with the smallest WBI value are considered as the weakest bonds (in bold) for the bonds of same type in each compound and their BDEs are calculated and listed in Table Ⅵ.

TABLE Ⅴ The Wiberg bond index of the title compounds from NBO analysis.
TABLE Ⅵ The calculated BDEs of the selected bonds of the title compounds.

The electronegativities [52] of -NO2,-N3,and -NH2 groups are 3.30,3.04,and 2.42,respectively. Obviously,there are significant inductive effects between the groups and other parts of the molecules. The electron withdrawing groups attract electrons,so the electron populations of their neighbor bonds C-R are reduced. Thus,if a bond is closer to electron withdrawing group or the group is more electronegative,then the bond is weaker. According to the ascending order of electronegativities,the ascending order for BDE of the C-N bonds should be DNTNDI,compound 2,and compound 1,as Table Ⅵ shows.

For DNTNDI,compounds 1,and 2,the N-N bonds have the smallest BDE values,but compound 3 locates its weakest bond at the C-O bond in the dioxin ring. The breaking of these bonds may be the first step of thermal decomposition of the compounds. Moreover,the -NH2 groups make the compounds the most stable,while the -NO2 groups make the compounds the most thermodynamic reactive. Based on data in Table Ⅵ,DNTNDI is the least stable compound,while compound 1 is the most stable compound.

Ⅳ. CONCLUSION

In this work,we investigated the molecular structures,infrared spectroscopy,heats of formation,detonation properties,chemical and thermal stabilities for a series of tetrahydro-[1, 4]dioxino[2,3-d:5,6-d'] diimidazole derivatives with four different substituents using hybrid DFT method. Results reveal that the compounds substituted with -NO2 and -ONO2 groups have larger densities than those with -NH2 and -N3 groups. The -N3 group is a better substituent than the others for increasing the HOF. Most of the strong infrared peaks are associated with the deformation of the substituent groups and the dioxin rings,which indicates that these parts are relatively more reactive areas in the molecules.

The D and P values for the substitution with -NO2 group are closer to those of HMX,the D and P values of the derivative with -ONO2 group are lower than those of RDX. The investigation for calculated ABEs,HOMO-LUMO gaps and the electronegativities of the substituent groups exhibits that the electronegativity of substituent group directly affects the chemical stability of the compound. For the dioxin-imidazole derivatives,a larger electronegativity for the substituent leads to a lower HOMO-LUMO gap and a smaller ABE. The calculations for BDEs indicate that the N-N bonds are the weakest bonds for DNTNDI and compounds 1,2. The C-O bonds in the dioxin rings have the smallest BDE for the compound 3. Considering the detonation properties and thermal stabilities,the DNTNDI can be candidate for high energy density materials.

Ⅴ. ACKNOWLEDGMENTS

This work was supported the National Natural Science Foundation of China (No.U1304111),and the Program for Science & Technology Innovation Talents in Universities of Henan Province (No.14HASTIT039) and the Innovation Team of Henan University of Science and Technology (No.2015XTD001).

References