Computational studies on the crystal structure, thermodynamic properties, detonation performance, and pyrolysis mechanism of 2,4,6,8-tetranitro-1,3,5,7-tetraazacubane as a novel high energy density material

Studies have suggested that octanitrocubane (ONC) is one of the most powerful non-nuclear high energy density material (HEDM) currently known. 2,4,6,8-Tetranitro-1,3,5,7-tetraazacubane (TNTAC) studied in this work may also be a novel HEDM due to its high nitrogen content and crystal density. Density functional theory and molecular mechanics methods have been employed to study the crystal structure, IR spectrum, electronic structure, thermodynamic properties, gas-phase and condensed-phase heat of formation, detonation performance, and pyrolysis mechanism of TNTAC. The TNTAC has a predicted density of about 2.12 g/cm(3), and its detonation velocity (10.42 km/s) and detonation pressure (52.82 GPa) are higher than that of ONC. The crystalline packing is P2(1)2(1)2(1), and the corresponding cell parameters are Z = 4, a = 8.87 ?, b = 8.87 ?, and c = 11.47 ?. Both the density of states of the predicted crystal and the bond dissociation energy of the molecule in gas phase show that the cage C-N bond is the trigger bond during thermolysis. The activation energy of the pyrolysis initiation reaction obtained from the B3LYP/6-311++G(2df,2p) level is 125.98 kJ/mol, which indicates that TNTAC meets the thermal stability request as an exploitable HEDM.     

多硝基四面体烷结构和性能的理论研究

在B3LYP/6-31G**水平下, 对四种四面体烷硝基衍生物进行理论研究. 基于全优化构型, 计算其红外光谱(IR)、热力学性质; 通过设计合理等键反应计算其气相生成热(HOF); 运用Kamlet-Jacobs方程估算其爆速(D)和爆压(p); 通过计算和比较各化合物的两种可能引发键(C—C和C—N)离解能(EC—C和EC—N), 确认该系列化合物的热解引发键和热稳定性. 讨论了各性能参数与其结构参数的关系. 兼顾高能量密度化合物(HEDC)的能量性质和稳定性要求, 最终认为该系列化合物不可作为潜在HEDC.     

Computational characterization of 2, 4, 6, 8, 10, 12, 13-heptaazatetracyclo [5.5.1.03,11.05,9]tridecanes as potential energetic compounds

The characters of high density and high heat of formation of cage molecules have attracted a lot of investigations as potential energetic materials. Several such compounds have been synthesized, e.g., octanitrocubane, hexanitrohexaazaisowurzitane (CL-20), and 4-trinitroethyl-2, 6, 8, 10, 12-pentanitrohexaazaisowurtzitane(TNE-CL-20). In the present study, a new cage compound, namely 2, 4, 6, 8, 10, 12, 13-heptaazatetracyclo [5.5.1.0^3,11.0^5,9] tridecane (HATT), was proposed. Density functional theory has been employed to study the geometric and electronic structures for a series of nitro derivatives of HATT at the B3LYP/6-31G(d,p) level. Thermodynamic properties derived on the basis of statistical thermodynamic principles are linearly correlated with the numbers of nitro group as well as the temperature. Detonation performance was evaluated based on the calculated densities and heats of formation. It is found that some title compounds have high densities of ca. 1.9 g cm^?3, detonation velocities over 9.0 km s^?1, and detonation pressures of about 40.0 GPa and may be novel potential candidates of high energy density compounds (HEDCs). Thermal stability and pyrolysis mechanism of the nitro HATTs were investigated by calculating the bond dissociation energies (BDE). In conjunction with the detonation performance and thermal stability, HATTs with no less than five nitro groups are recommended as the preferred candidates of HEDCs. These results provide basic information for the further studies of cage compounds.     

Computational study on the crystal structure, thermodynamic properties, detonation performance and pyrolysis mechanism of a novel high density cage compound 10-(5-nitrimino-1,2,3,4-tetrazol-1-yl)methyl-2,4,6,8,12-pentanitrohexaazaisowurtzitane

A novel polynitro cage compound 10-(5-nitrimino-1,2,3,4-tetrazol-1-yl)methyl-2,4,6,8,12-pentanitro-hexaazaisowurtzitane, composed of CL-20 and tetrazole framework, has been designed. DFT-B3LYP/6-31G(d) and molecular mechanics methods are employed to calculate its IR spectrum, heat of formation, thermodynamic properties, and crystal structure. Besides, the stability of this compound is evaluated using the bond dissociation energy. The result shows that the initial step of thermal decomposition is the rupture of N–NO₂ bond in the side chain. This compound is most likely to crystallize in the P-1 space group, and corresponding cell parameters are Z = 2, a = 7.65 Å, b = 14.30 Å, c = 10.36 Å, α = 91.53°, β = 50.83°, γ = 89.44°, and ρ = 2.025 g cm^?3. Detonation velocity and detonation pressure of this compound are estimated to be 9.090 km s^?1 and 38.078 GPa using the Kamlet–Jacobs equation, similar to those of CL-20. Considering detonation performance and thermal stability, this compound meets the requirements of exploitable high energy density materials.     

Quantum Chemical Studies on Structure and Detonation Performance of Bis(2,2-dinitropropyl ethylene)formal

Based on the full optimized molecular geometric structures via B3LYP/6-311+G(2d,p) method, a new gem-dinitro energetic plasticizer, bis(2,2-dinitropropyl ethylene)formal was investigated in order to search for high-performance energetic material. IR spectrum, heat of formation, and detonation performances were predicted. 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 the four N-NO₂ BDEs are nearly equal to the values of 164.38 kJ/mol, which shows that the title compound is a stable compound. The detonation velocity and pressure were evaluated by using Kamlet-Jacobs equations basedon the theoretical density and condensed HOF. The crystal structure obtained by molec-ular mechanics belongs to P2₁ space group, with lattice parameters Z=2, a=13.8017 ?, b=13.4072 ?, c=5.5635 ?.     

Studies on CL-20: The most powerful high energy material

CL-20 is an attractive HEM having density (>2 g cm^-3) and velocity of detonation (9400 m s^-1) superior to HMX (1.9 g cm^-3 and 9100 m s^-1). During this study, CL-20 was synthesized to establish viability of efficient synthesis method. The compound synthesized at HEMRL was characterized by FTIR, ¹H NMR and elemental analysis. Thermal studies (dynamic DSC and isothermal TG) were undertaken to determine kinetic parameters and understand the decomposition patterns. An attempt is made to explain the mechanism of decomposition of CL-20 on the basis of the data obtained by the authors and findings of other researchers. The activation energy values obtained during this work by adopting various approaches are close to the values reported for N-NO₂ bond cleavage suggesting that it is global rate determining process rather than the collapse of cage structure. Mass spectra also provides evidences in this regard. Monitoring of decomposition products at high temperature supports these findings and brings out that NO₂ initiates secondary decomposition processes because of entrapment in cage structure. This revised version was published online in August 2006 with corrections to the Cover Date.     

The effect of nitro groups on the structures and energetic properties of the derivatives composed of TATB and cubane

Based on the successful experience of synthesis of the TATB (1, 3, 5-triamino-2, 4, 6-trinitrobenzene) and cubane, we propose to consider their nitro derivatives combined by C–N bond as a series of high energy density compounds. First principles molecular orbital calculations have been used to investigate the structural and energetic properties, including the heat of formation, density, detonation performance, and impact sensitivity. Natural bond orbital analysis was carried out to investigate the influence of substituents on the electron delocalization. The results implied that the inclusion of nitro group will decrease the stability of cage skeleton and weaken the C–NO₂ bond. The calculated heats of formation, density, detonation velocity, and detonation pressure are positive and large. The results revealed that two of five derivatives have the close performance and sensitivity to those of CL-20, indicating that they may be explored as new potential high energy materials. Leave them with the notable value to dig out.     

Design and properties prediction of modified CL-20 energetic derivatives

We designed a series of energetic compounds based on the CL-20 molecular skeleton, and the properties including molecular geometric structures, electronic structures, density, heat of formation, detonation performances, and impact sensitivity were evaluated using density functional theory (DFT). The results indicate that five molecules have higher density values than that of Octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX; 1.91 g/cm³) and A4 has a larger density value (2.07 g/cm³) than that of CL-20 (2.04 g/cm³). In addition, most of the molecules have better detonation performances and stability than those of CL-20, with A4 showing much greater detonation velocity (9.93 km/s) and pressure (47.32 GPa) than those of CL-20 with a h₅₀ value of 14.02 cm. Taking both excellent detonation performance and low sensitivity into consideration, all seven compounds except for A3 and A5 are considered as potential energetic compounds. These theoretically calculated results would be conducive to the design and synthesis of novel nitramine energetic compounds.     

Computational Design of High Energy RDX-Based Derivatives: Property Prediction,Intermolecular Interactions,and Decomposition Mechanisms

A series of new high-energy insensitive compounds were designed based on 1,3,5-trinitro-1,3,5-triazinane (RDX) skeleton through incorporating -N(NO₂)-CH₂-N(NO₂)-, -N(NH₂)-, -N(NO₂)-, and -O- linkages. Then, their electronic structures, heats of formation, detonation properties, and impact sensitivities were analyzed and predicted using DFT. The types of intermolecular interactions between their bimolecular assemble were analyzed. The thermal decomposition of one compound with excellent performance was studied through ab initio molecular dynamics simulations. All the designed compounds exhibit excellent detonation properties superior to 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (CL-20), and lower impact sensitivity than CL-20. Thus, they may be viewed as promising candidates for high energy density compounds. Overall, our design strategy that the construction of bicyclic or cage compounds based on the RDX framework through incorporating the intermolecular linkages is very beneficial for developing novel energetic compounds with excellent detonation performance and low sensitivity.     

Spectroscopic, structural, and conformational properties of (Z)-4,4,4-trifluoro-3-(2-hydroxyethylamino)-1-(2-hydroxyphenyl)-2-buten-1-one, C12H12F3NO3: a trifluoromethyl-substituted β-aminoenone

The (Z)-4,4,4-trifluoro-3-(2-hydroxyethylamino)-1-(2-hydroxyphenyl)-2-buten-1-one (C(12)H(12)F(3)NO(3)) compound was thoroughly studied by IR, Raman, UV-visible, and (13)C and (19)F NMR spectroscopies. The solid-state molecular structure was determined by X-ray diffraction methods. It crystallizes in the P2(1)/c space group with a = 12.1420(4) ?, b = 7.8210(3) ?, c = 13.8970(5) ?, β = 116.162(2)°, and Z = 4 molecules per unit cell. The molecule shows a nearly planar molecular skeleton, favored by intramolecular OH···O and NH···O bonds, which are arranged in the lattice as an OH···O bonded polymer coiled around crystallographic 2-fold screw-axes. The three postulated tautomers were evaluated using quantum chemical calculations. The lowest energy tautomer (I) calculated with density functional theory methods agrees with the observed crystal structure. The structural and conformational properties are discussed considering the effect of the intra- and intermolecular hydrogen bond interactions.     

DFT studies on the four polymorphs of crystalline CL-20 and the influences of hydrostatic pressure on epsilon-CL-20 crystal

Based on density functional theory (DFT), four different methods with the generalized gradient approximation (GGA) have been employed to investigate the structural and electronic properties of the four polymorphs (alpha.H2O, beta, gamma, and epsilon phases) of CL-20, which is a well-known high energy density compound (HEDC). The relaxed crystal structures compare well with experimental data. According to the constitution of the frontier energy bands and the Mulliken population analyses, the N-NO2 bond is predicted to be the trigger bond during thermolysis. The density of states (DOS) of alpha-CL-20.H2O is somewhat different from those of the other three crystals for its inclusion of H2O molecules that contribute the frontier energy bands. The band gaps obtained from the four different methods are consistent with each other. According to the calculated values of band gaps, the sensitivity of the four polymorphs of CL-20 is predicted as epsilon < beta < gamma < alpha.H2O, which agrees well with the experimental result. The effects of hydrostatic compression on the most stable epsilon-CL-20 have also been investigated using the GGA-PBE method in the pressure range of 0-400 GPa. epsilon-CL-20 has anisotropic compressibility at low or high pressure. The band gap is found to decrease with increasing pressure, showing the corresponding sensitivity increase. Based on the changes of the band gap and DOS with pressure, 400 GPa is considered to be the critical pressure for the insulator-metal phase transition.     

几种新型六氮杂异伍兹烷衍生物结构与性能的理论预测—— 高能量密度化合物的寻求

以CN, NC, ONO2, N3, NH2, N2H, NHNH2, N4H和N4H3 9种含氮高能基团为取代基, 分别取代2,4,6,8,10,12-六氮杂异伍兹烷(IW)中亚氨基的6个H原子所形成的9种六氮杂异伍兹烷衍生物作为研究目标分子. 运用密度泛函理论, 在B3LYP/6-31G**水平上求得了它们的分子几何构型、电子结构、解离能(BDE)及IR谱等信息, 并设计等键反应计算了生成热( ). 基于统计热力学原理计算拟合了100～1200 K温度范围内体系的热力学函数, 利用Kamlet-Jacobs方程估算了它们的爆轰性能. 研究结果表明, 9种六氮杂异伍兹烷衍生物存在两种可能的热解引发类型. 在衍生物HNiIW, HBDAIW和HBAIW中, 可能的热解引发键是取代基内部的化学键, 而其余衍生物的热解引发键则可能是骨架N与取代基R之间N—R键. 另外, 硝酸酯基(ONO2)取代所得化合物HNiIW的密度ρ、爆速D及爆压p分别为1.998 g&#8226;cm－3, 9.71 km&#8226;s－1和44.47 GPa, 完全达到高能量密度化合物(HEDC)的基本要求, 且优于已应用的HNIW, 有望成为新型的HEDC.     

定常冲击波作用下六硝基六氮杂异伍兹烷(CL-20)/奥克托今(HMX)含能共晶初始分解机理研究

采用ReaxFF分子动力学方法同时结合多尺度冲击技术(MSST)模拟了4–10 km×s~(-1)定常冲击波加载下含能共晶CL-20/HMX沿不同晶格矢量的初始物理化学响应。获得了系统温度、压力、密度以及粒子速度的时间演化路径,以及初始分解路径,最终稳定反应产物和冲击雨贡纽等。研究结果表明:冲击波入射至含能共晶后,物理上依次经历诱导期、快压缩、慢压缩以及膨胀过程。快压缩和慢压缩过程分别对应反应物的快分解和慢分解。采用指数函数对反应物的衰减曲线进行拟合,并比较了共晶中反应物的衰减速率。整体上,随着冲击波速度的增加,反应物响应的时间逐渐提前,并且,冲击波沿各晶格矢量入射后,共晶中CL-20分子分解的响应时间均早于HMX。CL-20快分解阶段的衰减速率最高,HMX快分解的衰减速率居其次。相对于快分解阶段,慢分解阶段各反应物的衰减速率差异较小。含能共晶的初始反应路径是CL-20聚合形成二聚体,而冲击诱导共晶分解的初始反应路径是CL-20中N-NO_2键断裂形成NO_2。随后产生N_2O,NO,HONO,OH,H等中间小分子。最终稳定产物是N_2,H_2O,CO_2,CO和H_2。晶格矢量b,c方向冲击感度相同,低于晶格矢量a方向的感度。冲击诱导共晶中CL-20和HMX分解的最小冲击波速度(us)分别为6 km×s~(-1)和7 km×s~(-1)。采用冲击雨贡纽关系计算得到沿晶格矢量a,b,c冲击诱导CL-20/HMX共晶起爆的压力分别为16.52 GPa,17.41 GPa和17.41 GPa。爆轰压力范围介于36.75 GPa–47.43 GPa。     

Computational studies on 3,5,7,10,12,14,15,16-octanitro-3,5,7,10,12,14,15,16-octaaza-pentacyclo[7.5.1.12,8.04,13.06,11]hexadecane as potential high-energy-density compound

The B3LYP/6-31G(d) method of density functional theory was used to study molecular geometry, electronic structure, infrared spectrum, and thermodynamic properties. Detonation properties were evaluated using Kamlet–Jacobs equations based on the calculated density and heat of formation. Thermal stability of 3,5,7,10,12,14,15,16-octanitro-3,5,7,10,12,14,15,16-octaaza-pentacyclo[7.5.1.1^2,8.0^4,13.0^6,11]hexadecane (cage-HMX) was investigated by calculating the bond dissociation energy at unrestricted B3LYP/6-31G(d) level. The calculated results show that the first step of pyrolysis is the rupture of the N–NO₂ bond. The crystal structure obtained by molecular mechanics belongs to P2₁ space group, with lattice parameters a = 8.866 Å, b = 11.527 Å, c = 13.011 Å, Z = 4, and ρ = 2.219 g cm^?3. Both the detonation velocity of 9.79 km s^?1 and the detonation pressure of 45.45 GPa are better than those of CL-20. According to the quantitative standard of energetics and stability as a high-energy-density compound, cage-HMX essentially satisfies this requirement. These results provide basic information for molecular design of novel HEDCs.     

新型高能量密度化合物3,6-双(3,5-二硝基-1,2,4-三唑-1)-1,2,4,5-四嗪-1,4-二氧化物的性能预估及合成路线设计

采用C++自编译程序及组合原理,设计并筛选出一种未见报道的新型富氮类高能量密度化合物-3,6-双(3,5.二硝基.1,2,4-三唑.1)-1,2,4,5-四嗪-1,4-二氧化物,用B3LYP法,在6-31G**基组水平上得到该化合物全优化构型;在振动分析的基础上求得体系的振动频率、IR谱;通过键级分析得到热解引发键的键离解能(BDE);采用Monte-Carlo 方法预估了密度;设计等键等电子反应计算了生成焓;运用Kamlet-Jacobs公式预测爆速、爆压和爆热;运用Keshavarz 等推导的预估撞击感度H50的公式预测了撞击感度性能;并利用逆合成分析法设计其合成路线.结果表明:该化合物存在8个强吸收峰,校正后的热解引发键的BDE为264KJ·mol-1,稳定性较优;密度1.955 g·cm-3、生成焓901.72 kJ·mol-1、爆速9191.48 m·s-1、爆压39.32 GPa、爆热6705.15 j·g-1;撞击感度H50为55.85cm,低于黑索金(RDX)和奥克托今(HMx);以上性能均达到了高能量密度化合物的标准,且该化合物设计合成路线步骤较少、原料易得,有望得到广泛应用.     

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.     

A DFT study of cage compounds: 3, 5, 8, 10, 11, 12-hexanitro-3, 5, 8, 10, 11, 12-hexaazatetracyclo [5.5.1.12,6.04,9] dodecane and its derivatives as high energetic materials

A new cage compound, 3, 5, 8, 10, 11, 12-hexanitro-3, 5, 8, 10, 11, 12-hexaazatetracyclo [5.5.1.1^2,6.0^4,9] dodecane (HNHATCD, I) as well as its –ONO₂ (II) and –N₃ (III) derivatives were proposed in the present work. Their molecular structures were optimized at the B3LYP/6-31G(d,p) level of density functional theory. Heat of formation, strain energy, detonation performance, and thermal stability were studied. Results show that the –N₃ group greatly increases the heat of formation, but decreases the strain energy and density, and it is much more helpful for enhancing the detonation energy than the –NO₂ and –ONO₂ groups. An analysis of bond dissociation energies (BDEs) of the weakest bonds implies that the BDE of –N₃ derivatives is the smallest but it is still larger than 120 kJ mol^?1, revealing that these designed compounds have a high thermal stability. Considering the detonation performance and thermal stability, I and II may be potential candidates of high energy density materials.     

Synthesis and Crystal Structure of a Ruthenium(II) Compound of Dinitrogen

The title compound [Ru(bph)(PPh3)2(N2)]·(PF6)·1.5CH2Cl2 (bph = N-(benzoyl)- N?-(picolinylidene)-hydrazine) has been prepared and characterized by X-ray diffraction analysis. It crystallizes in the monoclinic system, space group P21/n with a = 13.5201(6), b = 15.5078(6), c = 25.0656(13) ?, β = 100.257(3)°, C50.50H43Cl3F6N5OP3Ru, Mr = 1150.23, V = 5171.4(4) ?3, Z = 4, F(000) = 2332, Dc = 1.477 g/cm3 and μ(MoKα) = 0.615 mm-1. The structure was refined to R = 0.0524 and wR = 0.1398 for 8060 observed reflections (I > 2σ(I)). The Ru-N2 and N(N bond distances are 1.919(4) and 1.133(5) ?, respectively.     

3,6-二肼基-1,2,4,5-四嗪的晶体结构及理论研究

通过缓慢蒸发溶剂法培养得到3, 6-二肼基-1, 2, 4, 5-四嗪(DHT)的单晶, 用X-射线单晶衍射仪进行了结构测定. 该晶体属于单斜晶系, P2(1)/c空间群, a=4.032 (4) ?, b=5.649 (6) ?, c=12.074 (14) ?. β=99.32°, Z=2, V=271.4(5) ?3. DHT分子中肼基N原子与四嗪环基本位于同一平面呈现轴对称结构, 分子中的大量氢键使之形成箭尾形排列的三维网络结构. 通过实验测得DHT的燃烧热为1787kJ?mol-1, 5s爆发点为454 K. 在DFT-B3LYP/6-311G*水平下对DHT的电子结构和自然键轨道进行了分析. 通过原子化能的方法计算得到DHT的标准生成热为1075 kJ?mol-1, 与实验值接近. 爆轰性能计算表明, DHT在密度为1.64 g?cm-3时, 爆速和爆压分别为9.27 km?s-1和36.02 GPa, 高于TNT和HMX.     

Density functional theory studies of effects of boron replacement on the structure and property of RDX and HMX

In this study, based on two model nitramine compounds hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) and octahydro-1,3,5,7-tetranitro-1,3,5, 7-tetrazocine (HMX), two series of new energetic molecules were designed by replacing carbon atoms in the ring with different amounts of boron atoms, their structures and performances were investigated theoretically by the density functional theory method. The results showed that the boron replacement could affect the molecular shape and electronic structure of RDX and HMX greatly, and then would do harm to the main performance like the heat of formation, density, and sensitivity. However, the compound RDX-B2 is an exception; it was formed by replacing two boron atoms into the system of RDX and has the symmetric boat-like structure. Its oxygen balance (4.9%), density (1.91 g/cm³), detonation velocity (8.85 km/s), and detonation pressure (36.9 GPa) are all higher than RDX. Furthermore, RDX-B2 has shorter and stronger N NO₂ bonds than RDX, making it possesses lower sensitivity (45 cm) and better thermal stability (the bond dissociation energy for the N NO₂ bond is 204.7 kJ/mol) than RDX. Besides, RDX-B1 and HMX-B4 also have good overall performance; these three new molecules may be regarded as a new potential candidate for high energy density compounds.