反式-1,4,5,8-四硝基-1,4,5,8-四氮杂萘烷异构体热解机理与稳定性的理论研究

运用密度泛函理论和半经验分子轨道方法,对一系列高能杂环硝胺—反式-1,4,5,8-四硝基-1,4,5,8-四氮杂萘烷异构体的热解机理和稳定性进行了系统地计算研究。在B3LYP/6-31G**和PM3水平上,分别计算了标题物的化学键离解能(BDE)和热解反应活化能(Ea),并根据BDE和Ea数值考察了硝胺取代基对化合物稳定性和热解机理的影响;同时,还详细考察了BDE与Ea、化学键重叠布居数、前线轨道能级以及能隙之间的相关性。结果表明,由BDE、Ea和静态电子结构参数推断的标题物热稳定性和热解机理的结论基本是一致的,N-NO2键均裂是标题物的热解引发步骤,间位取代异构体较对位取代异构体稳定,而邻位取代的异构体稳定性最差。     

Molecular dynamics simulations on the structures and properties of ε-CL-20-based PBXs——Primary theoretical studies on HEDM formulation design

Five polymer bonded explosives (PBXs) with the base explosiveε-CL-20 (hexanitrohexaazaisowurtzitane), the most important high energy density compound (HEDC), and five polymer binders (Estane 5703, GAP, HTPB, PEG, and F2314) were constructed. Molecular dynamics (MD) method was employed to investigate their binding energies (Ebind), compatibility, safety, mechanical properties, and energetic properties. The information and rules were reported for choosing better binders and guiding formulation design of high energy density material (HEDM). According to the calculated binding energies, the ordering of compatibility and stability of the five PBXs was predicted as ε-CL-20/PEG ＞ ε-CL-20/ Estane5703 ≈ε-CL-20/GAP ＞ ε-CL-20/HTPB ＞ ε-CL-20/F2314. By pair correlation function g(r) analyses, hydrogen bonds and vdw are found to be the main interactions between the two components. The elasticity and isotropy of PBXs based ε-CL-20 can be obviously improved more than pure ε-CL-20 crystal. It is not by changing the molecular structures of ε-CL-20 for each binder to affect the sensitivity. The safety and energetic properties of these PBXs are mainly influenced by the thermal capability (C°p) and density (ρ) of binders, respectively.     

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 study of energetic nitrogen‐rich derivatives of 1,1′‐ and 5,5′‐bridged ditetrazoles

Density functional theory method was used to study the heats of formation (HOFs), electronic structure, energetic properties, and thermal stability for a series of bridged ditetrazole derivatives with different linkages and substituent groups. The results show that the ? N₃ group and azo bridge (? N?N? ) play a very important role in increasing the HOF values of the ditetrazole derivatives. The effects of the substituents on the HOMO–LUMO gap are combined with those of the bridge groups. The calculated detonation velocities and detonation pressures indicate that the ? NO₂, ? NF₂, ? N?N? , or ? N(O)?N? group 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 suggests that the N? N bond in the ring or outside the ring is the weakest one and the N? N cleavage is possible to happen in thermal decomposition. Overall, the ? CH₂? CH₂? or ? NH? NH? group is an effective bridge for enhancing the thermal stability of the bridged ditetrazoles. Because of their desirable detonation performance and thermal stability, five compounds may be considered as the potential candidates of high‐energy density materials (HEDMs). These results provide basic information for the molecular design of novel HEDMs. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Theoretical studies on new potential high energy density compounds (HEDCs) adamantyl nitrates from gas to solid

A series of adamantyl nitrates have been theoretically studied from gas to solid to search for new po-tential high energy density compounds (HEDCs). The heats of formation (HOFs) for the 26 title com-pounds were calculated by designing isodesmic reactions at the B3LYP/6-31G level. It was found that the HOFs of the 26 isomers with the same number of —ONO2 groups (n) are not correlated well with the corresponding substituted positions. According to the obtained heats of detonation (Q),detonation velocities (D),and detonation pressures (P) using the Kamlet-Jacobs equations,it was found that when n=7～8,the adamantyl nitrates meet the criterion as an HEDC. The calculations on bond dissociation energies of O—N (EO—N) showed that the adamantyl nitrates with gemi —ONO2 always have the worst stability among the isomers,and all the adamantyl nitrates with gemi —ONO2 have similar stability. Due to the complexity of their structures,values of EO—N do not decrease with the increase of the substituent number n obviously,and the stability of adamantyl nitrates is not determined by only one structural parameter. Considering the stability requirement,only 1,2,4,6,8,9,10-adamantyl heptanitrate is recom-mended as a feasible HEDC. Molecular packing searching for 1,2,4,6,8,9,10-adamantyl heptanitrate among 7 most possible space groups (P21/c,P-1,P212121,P21,Pbca,C2/c,and Pna21) using Compass and Dreiding force fields showed that this compound tends to crystallize in P21/c. Ab initio periodic calculations on the electronic structure of the predicted packing showed that the O—NO2 bond is the trigger bond during thermolysis,which agrees with the result derived from the study of dissociation energies of O—N bonds.     

Theoretical predictions on pentaerythritol tetranitrate-based high energy density compounds

A novel family of pentaerythritol tetranitrate (PETN) derivatives based parent PETN skeleton were designed by introducing two energetic groups –NF₂ and –NO₂. Their electronic structure, heats of formation, detonation properties, impact sensitivity, and thermal stability were investigated by using density functional theory. The findings reveal that most of the title compounds have good detonation performance. The –NF₂ group played an important role in improving the densities, heats of detonation, and detonation properties of the designed molecules. The values of h₅₀ for almost all the PETN derivatives are higher than that of octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine. An analysis of bond dissociation energy suggests that the N-NO₂ bond tends to be a trigger bond in thermal decomposition. Taking both detonation properties and thermal stabilities into consideration, the three compounds may be selected as potential high-energy-density compounds.     

Theoretical Studies on the Structures,Thermodynamic Properties,Detonation Performance,and Pyrolysis Mechanisms for Six Dinitrate Esters

Density functional theory (DFT) has been employed to study the geometric and electronic structures of six dinitrate esters including ethylene glycol dinitrate (EGDN), diethylene glycol dinitrate (Di‐EGDN), triethylene glycol dinitrate (Tri‐EGDN), tetraethylene glycol dinitrate (Tetra‐EGDN), pentaethylene glycol dinitrate (Penta‐EGDN) and hexaethylene glycol dinitrate (Hexa‐EGDN) at the B3LYP/6‐31G* level. Their IR spectra were obtained and assigned by vibrational analysis. Based on the frequencies scaled by 0.96 and the principle of statistic thermodynamics, the thermodynamic properties were evaluated, which were linearly related with the number of CH₂CH₂O groups as well as the temperature, obviously showing good group additivity. Detonation performances were evaluated by the Kamlet‐Jacobs equations based on the calculated densities and heats of formation. It was found that density, detonation velocity, detonation pressure decreased with the increase of the number of CH₂CH₂O groups. Thermal stability and the pyrolysis mechanism of the title compounds were investigated by calculating the bond dissociation energies (BDE) at the B3LYP/6‐31G* level. For the nitrate esters, the O‐NO₂ bond is a trigger bond during a thermolysis initiation process.     

New theoretically predicted RDX‐ and β‐HMX‐based high‐energy‐density molecules

Theoretically new high‐energy‐density materials (HEDM) in which the hydrogens on RDX and β‐HMX (hexahydro‐1,3,5‐trinitro‐1,3,5‐triazine and octahydro‐1,3,5,7‐tetranitro‐1,3,5,7‐tetrazocine, respectively) were sequentially replaced by (N NO₂)x functional groups were designed and evaluated using density functional theory calculations in combination with the Kamlet–Jacobs equations and an atoms‐in‐molecules (AIM) analysis. Improved detonation properties and reduced sensitivity compared to RDX and β‐HMX were predicted. Interestingly, the RDX and β‐HMX derivatives having one attached N NO₂ group [RDX‐(NNO₂)1 and HMX‐(NNO₂)1] showed excellent detonation properties (detonation velocities: 9.529 and 9.575 km·s⁻¹, and detonation pressures: 40.818 and 41.570 GPa, respectively), which were superior to the parent compounds. Sensitivity estimations obtained by calculating impact sensitivities and HOMO‐LUMO gaps indicated that RDX‐(NNO₂)1 and HMX‐(NNO₂)1 were less stable than RDX and HMX but more stable than any of the other derivatives. This method of sequential NNO₂ group attachment on conventional HEDMs offers a firm basis for further studies on the design of new explosives. Furthermore, the newly found structures may be promising candidates for better HEDMs.     

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 the structures, thermodynamic properties, detonation properties, and pyrolysis mechanisms of spiro nitramines

Density function theory (DFT) has been employed to study the geometric and electronic structures of a series of spiro nitramines at the B3LYP/6-31G level. The calculated results agree reasonably with available experimental data. Thermodynamic properties derived from the infrared spectra on the basis of statistical thermodynamic principles are linearly correlated with the number of nitramine groups as well as the temperature. Detonation performances were evaluated by the Kamlet-Jacobs equations based on the calculated densities and heats of formation. It is found that some compounds with the predicted densities of ca. 1.9 g/cm3, detonation velocities over 9 km/s, and detonation pressures of about 39 GPa (some even over 40 GPa) may be novel potential candidates of high energy density materials (HEDMs). Thermal stability and the pyrolysis mechanism of the title compounds were investigated by calculating the bond dissociation energies (BDE) at the B3LYP/6-31G level and the activation energies (E(a)) with the selected PM3 semiempirical molecular orbital (MO) based on the unrestricted Hartree-Fock model. The relationships between BDE, E(a), and the electronic structures of the spiro nitramines were discussed in detail. Thermal stabilities and decomposition mechanisms of the title compounds derived from the B3LYP/6-31G BDE and the UHF-PM3 E(a) are basically consistent. Considering the thermal stability, TNSHe (tetranitrotetraazaspirohexane), TNSH (tetranitrotetraazaspiroheptane), and TNSO (tetranitrotetraazaspirooctane) are recommended as the preferred candidates of HEDMs. These results may provide basic information for the molecular design of HEDMs.     

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.     

Molecular design of a new family of bridged bis(multinitro-triazole) with outstanding oxygen balance as high-density energy compounds

A new family of bridged bis(multinitro-triazole) was designed and investigated using the density functional theory method. The density, oxygen balance, heat of formation, detonation performance, and impact sensitivity were calculated systematically. The results show that the multinitromethyl groups play an important role in increasing densities. At the same time, different bridged groups present diverse performances with high density (1.86-1.96 g·cm⁻³), excellent detonation properties (V = 8.72 km·s⁻¹-9.20 km·s⁻¹; P = 34.54 GPa-39.49 GPa), outstanding oxygen balance (0%-11.59%), and acceptably impact sensitivity. Especially, tetrazine (M7)-bridged and diaminofurazan (M9)-bridged groups are very helpful for enhancing their detonation performance (V_(M7) = 9.12 km·s⁻¹, P_(M7) = 38.51 GPa; V_(M9) = 9.20 km·s⁻¹, P_(M9) = 39.49 GPa), respectively, which are better than RDX. They could be seen as the potential candidates of high energy density materials (HEDMs).     

Predicting Trigger Bonds in Explosive Materials through Wiberg Bond Index Analysis

Understanding the explosive decomposition pathways of high‐energy‐density materials (HEDMs) is important for developing compounds with improved properties. Rapid reaction rates make the detonation mechanisms of HEDMs difficult to understand, so computational tools are used to predict trigger bonds—weak bonds that break, leading to detonation. Wiberg bond indices (WBIs) have been used to compare bond densities in HEDMs to reference molecules to provide a relative scale for the bond strength to predict the activated bonds most likely to break to trigger an explosion. This analysis confirms that X?NO₂ (X=N,C,O) bonds are trigger linkages in common HEDMs such as TNT, RDX and PETN, consistent with previous experimental and theoretical studies. Calculations on a small test set of substituted tetrazoles show that the assignment of the trigger bond depends upon the functionality of the material and that the relative weakening of the bond correlates with experimental impact sensitivities.     

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     

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

Structurally uneasy but kinetically stable nitrogens in 1,3-disubstituted cyclotetrazenes: Viable high-energy-density materials

In designing the polynitrogen and nitrogen-rich high-energy-density materials (HEDMs), one continuously challenged issue is to solve the severe conflict between “stability and detonation.” Most known nitrogen-containing and nonsaltlike HEDMs involve the normally bonded skeletal nitrogens. Here, we computationally proposed a nontraditional class of nitrogen-rich compounds with structurally uneasy nitrogens, that is, 1,3-disubstituted cyclotetrazenes N₄R₂ ( I ), which bear the moderate 6π-aromaticity and slight singlet diradical character. The target N₄R₂ ( I ) can preferably balance the “stability and detonation” necessity of HEDM, and can be fabricated into diversified compounds ranging from simple to nano-sized high-energy porous frameworks. The robustness of N₄R₂ ( I ) should much expand our perceiving limitation in the nitrogen-rich HEDM chemistry.     

Theoretical studies on heats of formation,detonation properties,and bond dissociation energies of monofurazan derivatives

The heats of formation (HOFs) for a series of monofurazan derivatives were calculated by using density functional theory. It is found that the ? CN or ? N₃ group plays a very important role in increasing the HOF values of the furazan derivatives. The detonation velocities and detonation pressures of the furazan derivatives are evaluated at two different levels. The results show that the ? NF₂ group is very helpful for enhancing the detonation performance for the furazan derivatives, but the case is quite the contrary for the ? CH₃ group. An analysis of the bond dissociation energies and bond orders for the weakest bonds indicate that the substitutions of ? CN group are favorable and enhances the thermal stability of the furazan derivatives, but the ? NO₂ groups produce opposite effects. These results provide basic information for the molecular design of novel high‐energy density materials. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Theoretical studies on the structures, thermodynamic properties, detonation properties, and pyrolysis mechanisms of four trinitrate esters

Density function theory (DFT) has been employed to study the geometric and electronic structures of four trinitrate ester including nitroglycerin (NG), butanetriol trinitrate (BTTN), trimethanolethane trinitrate (TMETN) and trimethylolpropane trinitrate (TMPTN) at the B3LYP/6-31G^* level. Their IR spectra were obtained and assigned by vibrational analysis. Based on the frequencies scaled by 0.96 and the principle of statistic thermodynamics, the thermodynamic properties were evaluated, which were linearly related with the number of methylene groups as well as the temperature, obviously showing good group additivity. Detonation performances were evaluated by the Kamlet–Jacobs equations based on the calculated densities and heats of formation. It is found that density, detonation velocity, detonation pressure are decrease with the increase of the number methylene groups. Thermal stability and the pyrolysis mechanism of the title compounds were investigated by calculating the bond dissociation energies (BDE) at the B3LYP/6-31G^* level. For the nitrate esters, the ONO₂ bond is a trigger bond during thermolysis initiation process.     

Solvent effects on the geometry,electronic structure,and bonding style of Zn(N5)2: A theoretical study

Nitrogen-rich compounds involving the cyclo-pentazole anion (cyclo-N₅⁻) have attracted extensive attention due to higher energy release and environmental friendliness than traditional high energy density materials (HEDMs). However, the synthesis of stable HEDMs with cyclo-N₅⁻ is still a challenge. In this study, the effect of nine solvents on the geometrical and electronic structures and solvation energies of Zn(N₅)₂, one of the recently synthesized nitrogen-rich compounds, was studied using the density functional theory and the polarized continuum model. The results indicated an increase in the stability of Zn(N₅)₂ in the solution phase compared to the vacuum phase, and the stability of Zn(N₅)₂ increases with increasing dielectric constants. The energy gap of frontier molecular orbitals and the absolute value of total energy in water are the largest, revealing that Zn(N₅)₂ is more stable in water than in other solvents. To understand the stabilization mechanism of Zn(N₅)₂ by water, further studies were performed with the natural bond orbital (NBO) analysis and the quantum theory of atoms in molecules (QTAIM) analysis using the explicit solvent model. The charge transfer and the hydrogen bonds are observed between Zn(N₅)₂ and water, which are beneficial to improvement of the stability of Zn(N₅)₂. This may indicate the solvents that have strong interactions with the cyclo-N₅⁻ candidate may improve the possibility of success of synthesis.