含能材料的密度、爆速、爆压和静电感度的理论研究

用密度泛函理论(DFT) B3LYP方法, 在6-31G*基组水平下, 全优化计算了系列硝胺类和硝基芳烃类爆炸物的几何构型, 用Monte-Carlo方法和自编程序, 基于0.001 e•bohr^－3等电子密度面所包围的体积空间求得分子平均摩尔体积(V)和理论密度(ρ). 用Kamlet-Jacobs方程基于理论密度(ρ)和PM3计算生成焓(ΔH_f)估算标题物的爆速(D)和爆压(p), 发现多环硝胺类化合物的爆轰性能优于芳烃硝基类化合物, 故此, 在寻求高能量密度材料(HEDM)时, 我们应特别关注多环硝胺化合物. 与ρ和D文献值比较, 表明本理论计算方法和结果是适用可靠的. 将爆速(D)和爆压(p)计算值与静电感度实验值(E_ES)进行比较和关联, 发现: 若对化合物进行细致分类讨论, 则它们之间存在较好的线性关系. 据此建议, 在含能材料分子设计中, 可通过理论计算爆轰性质(D或p)去预估难以定量求得或尚未合成的含能材料的静电火花感度值(E_ES). 此外, 我们还讨论了取代基对ρ, D和p的影响, 也有助于分子设计.     

Detonation Properties of High Explosives Containing Ammonia Borane

Ammonia borane (AB) is used as a combustion agent to improve the properties of high explosives. The detonation velocity (D_v) and detonation pressure (P) of raw high explosives and of samples containing AB were calculated and compared. The detonation properties, impact sensitivities, thermal sensitivities, and thermal decomposition characteristics of high explosives containing AB were also measured. The results indicated that when the AB content was 20 wt‐%, the optimal detonation velocity and detonation pressure were achieved. Both the detonation velocity and detonation pressure of the high explosives containing AB were clearly increased compared with those of the raw high explosives. Moreover, the detonation velocities of high explosives containing AB were 7078 to 7423 m · s^–1 and their density ranged from 1.570 to 1.589 g · cm^–3. The detonation pressure ranged from 34.5 to 37 GPa and the average heat of detonation was 6688 J · g^–1. Furthermore, the impact and thermal sensitivities were 170 cm and 613 K, respectively, whereas a slight change occurred in the thermal decomposition characteristics. These results suggest that AB can serve as a powerful combustible agent in energetic materials and improve the detonation properties and sensitivities of high explosives.     

Stepwise Degradation of Polyfunctional Compounds

The formal kinetic peculiarities of the decomposition of polyfunctional compounds containing several identical reaction groups were used to determine the stepwise character of secondary reactions and to interpret their chemical mechanism. It was found that the degradation of difluoramines, azides, nitro esters, and aliphatic and aromatic nitro compounds always occurred stepwise via long-lived intermediate products. Only nitramines with closely spaced reaction centers exhibited deviations from this mechanism.     

Molecular Structure and a Density Functional Theoretical Study on 2‐Nitroimino‐5‐Nitro‐Hexahydro‐1,3,5‐Triazine (NNHT)

2‐Nitroimino‐5‐nitro‐hexahydro‐1,3,5‐triazine (NNHT), was synthesized and its structure was determined by single‐crystal X‐ray diffraction. The crystal is monoclinic, space group P2₁/c with crystal parameters of a = 9.4031(13) Å, b = 8.5891(12) Å, c = 9.0200(13) Å, β = 91.213(2)°, V = 728.33(18) ų, Z = 4, F(000) = 392, D_c = 1.734 g/cm³. The experimental geometry of NNHT was input to Gaussian‐03W program and optimized using DFT‐B3LYP/6‐311++G** method. The IR frequencies and NMR chemical shift were carried out and compared well with those of the experimental. The atomic net charges and the population analysis are discussed. The heat of formation (HOF) for NNHT was evaluated by designing an isodesmic reaction. The detonation velocity (D) and detonation pressure (P) were estimated by using the well known Kamlet‐Jacobs equation, based on the theoretical HOF.     

Synthesis of Amino,Azido, Nitro,and Nitrogen‐Rich Azole‐Substituted Derivatives of 1H‐Benzotriazole for High‐Energy Materials Applications

The amino, azido, nitro, and nitrogen‐rich azole substituted derivatives of 1H‐benzotriazole have been synthesized for energetic material applications. The synthesized compounds were fully characterized by ¹H and ¹³C NMR spectroscopy, IR, MS, and elemental analysis. 5‐Chloro‐4‐nitro‐1H‐benzo[1,2,3]triazole ( 2 ) and 5‐azido‐4,6‐dinitro‐1H‐benzo[1,2,3]triazole ( 7 ) crystallize in the Pca2₁ (orthorhombic) and P2₁/c (monoclinic) space group, respectively, as determined by single‐crystal X‐ray diffraction. Their densities are 1.71 and 1.77 g cm^?3, respectively. The calculated densities of the other compounds range between 1.61 and 1.98 g cm^?3. The detonation velocity (D) values calculated for these synthesized compounds range from 5.45 to 8.06 km s^?1, and the detonation pressure (P) ranges from 12.35 to 28 GPa.     

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.     

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.     

多硝基尿酸衍生物含能性质的理论研究

为了寻找兼具优异爆轰性能和良好热力学及动力学稳定性的高能材料, 本文设计了15个硝基尿酸化合物, 运用密度泛函理论, 对其性质进行了研究. 通过半经验的K-J方程和比冲量预测了其爆炸性能, 结果表明, 所设计分子的爆热、 分子密度、 爆炸速率和爆炸压强同硝基取代基数目之间存在较强的线性关系. 三硝基尿酸和四硝基尿酸衍生物的爆炸速率超过了8.0 km/s, 爆炸压强超过了30 GPa, 并且大多数衍生物的比冲量要高于目前经常使用的炸药黑索金. 通过计算N—NO₂键的解离能、 特征落高、 分子的自由空间预判了衍生物的稳定性和撞击感度, 结果显示, 绝大多数分子有大于80 kJ/mol的键解离能. 本文的理论结果可以为实验上设计合成新的高能材料提供一些有用的信息.     

One pot synthesis of <Emphasis Type="Italic">N</Emphasis>-ethylaniline from nitrobenzene and ethanol

A novel method for the one pot synthesis of N-alkyl arylamines from nitro aromatic compounds and alcohols is proposed through the combination of the aqueous-phase reforming of alcohol for hydrogen production, the reduction of nitro aromatic compounds for the synthesis of aromatic amine and the N-alkylation of aromatic amine for the production of N-alkyl arylamine over an identical catalyst under the same conditions of temperature and pressure in a single reactor. In this process, hydrogen generated from the aqueous-phase reforming of alcohols was used in-situ for the hydrogenation of nitro aromatic compounds for aromatic amine synthesis, followed by N-alkylation of aromatic amine with alcohols to form the corresponding N-alkyl arylamines at a low partial pressure of hydrogen. For the system composed of nitrobenzene and ethanol, under the conditions of 413 K and P _N2 = 1 MPa, the conversion degrees of nitrobenzene and aniline were 100%, the selectivity to N-ethylaniline and N, N-diethylaniline were 85.9% and 0%–4%, respectivity, after reaction for 8 h at the volumetric ratio of nitrobenzene:ethanol:water = 10:60:0. The selectivity for N, N-diethylaniline production is much lower than that through the traditional method. In this process, hydrogen and aromatic amines generated from the aqueous-phase reforming of alcohols and hydrogenation of nitro aromatic compounds, respectively, could be promptly removed from the surface of the catalyst due to the occurrence of in-situ hydrogenation and N-alkylation reactions. Thus, this may be a potential approach to increase the selectivity to N-alkyl arylamine. Supported by the Program for New Century Excellent Talents in University (Grant No. NCET-04-0557), and the Specialized Research Fund for the Doctoral Program of High Education (Grant No. SRFDP-20060337001)     

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.     

A New Method for Assessment of Performing Mechanical Works of Energetic Compounds by the Cylinder Test

A new method is introduced for assessment of performing mechanical works of energetic compounds by cylinder wall velocities of CHNOFCl energetic compounds on the basis of the cylinder test. Four suitable decomposition paths are used to evaluate the number of moles of gaseous detonation products per gram of explosive, the average molecular weight of these gases, and the heat of detonation in calories per gram by considering different decomposition products HF, HCl, CO, N₂, H₂O, H₂, and CO₂. For CHNO and fluoro energetic compounds, the predicted cylinder wall velocities of these compounds give more reliable results than one of the best available empirical methods. The predicted root mean square (rms) deviations of cylinder wall velocities of the new model for some chloro explosives at actual radial expansions 0.6 and 1.9 mm are 0.010 and 0.062 km · s^–1, which show high reliability of the new method.     

The New Descriptors Effective on the Detonation Performance of Aluminized Explosives

The relationship between detonation velocity and the elemental composition of components of aluminized explosives are assessed through quantitative structure-property relationship (QSPR). Here, two new reliable, simple models are proposed for estimating aluminized explosives detonation heat and velocity based on molecular structure by applying QSPR. In this methodology it is assumed that these two detonation parameters can be presented as a function of elemental composition, density and several structural parameters. This new correlation of heat detonation has determination coefficient of 0.930, root mean square deviation (RMSD) of 324.4 and average absolute deviation (AAD) of 446kJ · kg^–1 for 36 aluminized explosives with different molecular structures as the training set. The predictive power of this new correlation is checked through a cross validation method. Statistical parameters reveal relatively good result for this correlation. Also, the determination coefficient of detonation velocity for the other new model is 0.960 and it has 151.1 (RMSD) and 107.9 m · s^–1 (AAD) for 42 aluminized explosives with different molecular structures as training set. Reliability and validity of new correlation investigated (Q²_Ext = 0.948, Q²_LOO = 0.938, and Q²_LMO = 0.937). The good ability of this new model for prediction detonation velocity of aluminized explosives are confirmed.     

An empirical method for predicting detonation pressure of CHNOFCl explosives

This paper describes a new method for prediction of the Chapman-Jouguet detonation pressures of CHNOFCl explosives using the heat of detonation, Q_det, the number of moles of gaseous products of detonation per gram of explosive, α, and the average molecular weight of gaseous products, M. The equation has the form: P_CJ=15.88α(MQ_det)^1/2ρ₀²−11.17, where P_CJ is the Chapman-Jouguet detonation pressure and ρ₀ the loading density. Calculated P_CJ by this procedure show good result with respect to measured detonation pressure for any pure or mixture of ideal and some of less ideal CHNOFCl explosives at ρ₀>0.8 g/cm³.     

Electron-impact fragmentation of some secondary aliphatic nitramines. Migration of the nitro group in heterocyclic nitramines

The mass spectra of a series of nitramines, including the explosives RDX and HMX, are reported. The nature of the major fragment ions is established by labeling with ¹⁵N or¹³C and by high resolution mass spectrometry. Electron-impact-induced fragmentation pathways including the migration of the nitro group in the heterocyclic nitramines, RDX and HMX, are postulated.     

Theoretical Study of the Mechanism of the Nitro-Nitrite Rearrangement and Its Role in Gas-Phase Monomolecular Decomposition of C-Nitro Compounds

Semiempirical, ab initio, and density functional theory calculations were used to study the primary act of gas-phase monomolecular decomposition of certain C-nitro compounds and their radical cations, associated with the nitro-nitrile rearrangement. It was shown that the reaction fails to occur with all neutral molecules of aliphatic nitro compounds (except for fluoronitromethane and fluoronitroethene) and has a much lower barrier with the corresponding radical cations. An important role of the nitro-nitrile rearrangement in gas-phase decomposition of aromatic nitro compounds was demonstrated.     

Thermal explosion of nitrophenates of transition metals

The thermal decomposition studies on nitrophenates of copper, nickel and cobalt have been undertaken,-t curves show dehydration of these compounds at lower temperatures whereas dehydration cum decomposition seem to occur at higher temperatures leading to oxidative combustion of aromatic part. NO₂ gas is evolved during decomposition which seems to be responsible for oxidative reactions leading to detonation. The explosion temperature and velocity of detonation have been found to be linearly related with the number of nitro groups. The mechanism of thermal explosion has also been discussed.Thanks are due to Prof. S. Giri for providing laboratory facilities; and Aeronautics R&D Board, Ministry of Defence for financial assistance. Authors are also grateful to Shri I. P. S. Kapoor for helpful discussions.     

Relationship between Electronic Charges at Nitrogen Atoms of Nitro Groups and Thermal Reactivity of Nitramines

Thermal reactivity of eleven nitramines has been examined by means of non-isothermal differential thermal analysis, and the data were analyzed according to the Kissinger method. The reactivity was expressed as the E _a R ⁻¹ slopes of the Kissinger relationship. Electronic charges, q ^N, at nitrogen atoms of the nitramine molecules were calculated by means of ab initio DFT B3LYP/6-31G** method. The relationships were confirmed between the slopes E _a R ⁻¹ and the q values for the nitro groups that are primarily split off. Conclusions are made in relation to the mechanism of initiation of polynitro compounds in general. This revised version was published online in July 2006 with corrections to the Cover Date.     

On the molecular origin of the sensitivity to impact of cyclic nitramines

Nitramine explosives can combine relative insensitivity to initiation and great energy content. In this work, based on a previous approach developed for nitroaromatic explosives, we propose four mathematical models to correlate impact sensitivity, given by the h₅₀ value, to molecular charge properties. Fourteen cyclic nitramines were studied using Density Functional Theory (DFT). Six molecules of the set have measured h₅₀ values, which were used to evaluate the sensitivity models. Converged DFT charge densities of the molecules were partitioned and analyzed according to the distributed multipole analysis (DMA) atom-centered method. The sensitivity models were based on the DMA electric multipole values. The electron withdrawing role of the nitro group and the strong polarization of the charges of the nitrogen atom in the amine group were clearly identified. The influence of the electronic properties on the sensitivity of the explosives was characterized by including in the sensitivity models the charge values of the nitro or the nitramine groups and electron delocalization, the latter quantified by the DMA quadrupole values of the ring atoms. Inclusion of electron delocalization effects can improve the prediction of h₅₀ values for two out of the five strained-ring nitramines in the set. The charge values of the nitramine groups are the most important molecular property affecting the impact sensitivity. The h₅₀ values of eight nitramine explosives of the set not available experimentally were computed.     

35Cl NQR studies of 1‐chloro‐2,4‐dinitrobenzene and 1,2‐dichloro‐3‐nitrobenzene as a function of pressure and temperature

The temperature and pressure dependences of ³⁵Cl nuclear quadrupole resonance (NQR) frequency and spin–lattice relaxation time (T₁) were investigated for 1‐chloro‐2,4‐dinitrobenzene and 1,2‐dichloro‐3‐nitrobenzene. T₁ was measured in the temperature range 77–300 K. Furthermore, the NQR frequency (ν) and T₁ for these compounds were measured as a function of pressure up to 5.1 kbar at 300 K. Relaxation was found to be due to the torsional motion of the molecule and the reorientation motion of the nitro group. By analysing the temperature dependence of T₁, the activation energy for the reorientation motion of the nitro group was obtained. The temperature dependence of the average torsional lifetimes of the molecules and the transition probabilities W₁ and W₂ for the Δm = ±1 and Δm = ±2 transitions, were also obtained. Both compounds showed a non‐linear variation of NQR frequency with pressure. The pressure coefficients were observed to be positive. A thermodynamic analysis of the data was carried out to determine the constant‐volume temperature coefficients of the NQR frequency. The spin–lattice relaxation time T₁ for both the compounds was found to be weakly dependent on pressure, showing that the relaxation is mainly due to the torsional motions. Copyright © 2002 John Wiley & Sons, Ltd.