NO_2,OH和OH~-对环四甲撑四硝胺初始热解的影响

在密度泛函理论的(DFT)B3LYP/6-31g(d)水平上,优化得到了环四甲撑四硝胺(β-HMX)及其与高氯酸铵(AP)裂解产物NO2、OH及OH-分别形成复合物的各种稳定构型,计算了β-HMX及各复合物中最弱的N—NO2键解离能.结果发现:β-HMX与NO2、OH结合后构型变化不是很大,但对称性降低;β-HMX与OH-结合后,HMX构型发生较大变化,原有的对称性明显遭到破坏.计算表明:NO2易与HMX骨架环上亚甲基(—CH2—)中的H作用,"置换"出H而引发HMX的热解,从而改变了HMX的初始分解通道;OH对HMX的N—NO2键解离影响不大,而OH-与β-HMX结合后其N—NO2键解离能比β-HMX降低近200kJ.mol-1,表明OH-对其裂解有明显的促进作用.NO2、OH-的存在可使HMX的分解温度大大降低.     

H+、NH+4对HMX的N—NO2键解离能的影响

采用密度泛函理论的B3P86/6-31G**方法, 优化了β-HMX及其与H⁺、NH^+₄分别形成的复合物的稳定结构, 计算了β-HMX以及复合物中最弱的N—NO₂键解离能. 结果发现, HMX与H⁺、NH^+₄形成复合物后, 使HMX的构型产生较大变化; 与H⁺结合后, HMX的一个N—NO₂键显著伸长, 键级变小; 但与NH^+₄形成复合物后, HMX中键级最小的N—NO₂键长变化不大. 键解离能计算表明, 同β-HMX相比, 与H⁺形成的两种复合物中N—NO₂键解离能分别降低了近20 和82 kJ·mol^-1, 而HMX与NH^+₄形成的复合物中N—NO₂键解离能仅降低了约8 kJ·mol^-1, 表明H⁺对β-HMX的N—NO₂键的初始热裂解反应有促进作用, 而NH^+₄影响不明显.     

CH3SH…HOO开壳型氢键复合物的结构及其电子密度拓扑性质

在DFT-B3LYP/6-311++G**水平下求得CH3SH…HOO复合物势能面上的稳定构型. 计算结果表明, 在HOO以其O8—H7作为质子供体与CH3SH分子中的S5原子为质子受体形成的氢键复合物1和2中, O8—H7明显被“拉长”, 且其伸缩振动频率发生显著的红移, 红移值分别为330.1和320.4 cm-1; 在CH3SH分子以其S5—H6作为质子供体与HOO的端基O9原子为质子受体形成的氢键复合物3和4中, 也存在类似的情况, 但S5—H6伸缩振动频率红移不大. 经MP2/6-311++G**水平计算的4种复合物含BSSE校正的相互作用能分别为-20.81, -20.10, -4.46和-4.52 kJ/mol. 自然键轨道理论(NBO)分析表明, 在CH3SH…HOO复合物1和2中, 引起H7—O8键长增加的因素包括两种电荷转移, 即孤对电子n1(S5)→σ*(H7—O8)和孤对电子n2(S5)→σ*(H7—O8), 其中后者为主要作用. 在复合物3和4中也有相似的电荷转移情况, 但轨道间的相互作用要弱一些. AIM理论分析结果表明, 4个复合物中的S5…H7间和O9…H6间都存在键鞍点, 且其Laplacian量▽2ρ(r)都是很小的正值, 说明这种相互作用介于共价键和离子键之间, 偏静电作用为主.     

Pd催化甲醇裂解制氢的反应机理

基于密度泛函理论(DFT), 研究了甲醇在Pd(111)面上首先发生O—H键断裂的反应历程(CH3OH(s)→CH3O(s)+H(s)→CH2O(s)+2H(s)→CHO(s)+3H(s)→CO(s)+4H(s)). 优化了裂解过程中各反应物、中间体、过渡态和产物的几何构型, 获得了反应路径上各物种的吸附能及各基元反应的活化能数据. 另外, 对甲醇发生C—O键断裂生成CH3(s)和OH(s)的分解过程也进行了模拟计算. 计算结果表明, O—H键的断裂(活化能为103.1 kJ·mol-1)比C—O键的断裂(活化能为249.3 kJ·mol-1)更容易; 甲醇在Pd(111)面上裂解的主要反应历程是: 甲醇首先发生O—H键的断裂, 生成甲氧基中间体(CH3O(s)), 然后甲氧基中间体再逐步脱氢生成CO(s)和H(s). 甲醇发生O—H键断裂的活化能为103.1 kJ·mol-1, 甲氧基上脱氢的活化能为106.7 kJ·mol-1, 两者均有可能是整个裂解反应的速控步骤.     

1,3,4,5,7,8-六硝基八氢化二咪唑[4,5-b∶4',5'-e]吡嗪-2,6-(1H,3H)-N,N'-二亚硝胺分子结构与性能的量子化学研究

设计了一种新型高能量密度化合物1,3,4,5,7,8-六硝基八氢化二咪唑[4,5-b∶4',5'-e]吡嗪-2,6-(1H,3H)-N,N'-二亚硝胺(ONIP).运用密度泛函理论(DFT),在B3PW91/6-31G++(d,p)水平下进行优化并计算出了ONIP的一些重要性质.通过键级的分析,母环的五元环侧链处N—NO2键为分解引发键,其解离能为107.8 kJ/mol;该化合物理论密度为2.00 g/cm3,生成热为1693.71 kJ/mol,爆速为10.21 km/s,爆压为49.17 GPa,表明爆轰性能优异;其撞击感度为33 cm,优于黑索金(RDX)、奥克托金(HMX)和六硝基六氮杂异伍兹烷(CL-20);能级差为3.67 eV,表明分子稳定性较高.给出了2条合成路线,均具有步骤少且原料易得的优点.     

奥克托金(HMX)的T-Jump/FTIR快速热裂解研究

采用T-Jump/FTIR快速热裂解原位红外光谱联用技术研究了奥克托金(HMX)在0.1,0.2,0.3和0.4MPa的Ar气条件下,以1000℃·s-1的升温速率快速升温至设定的反应温度,用快速扫描傅立叶变换红外光谱跟踪分析分解产物的种类和相对摩尔浓度的变化,研究了温度及压力对初始检测产物的影响.结果表明,HMX在快速热裂解5s过程中红外所检测到的主要气相产物为CO,CO2,NO,NO2,N2O,HCONH2,CH2O,H2O,HNCO及HCN,并给出了这些产物相对摩尔浓度随时间变化的曲线.根据气体产物相对摩尔浓度的比率N2O/HCN,研究了压力和反应温度对HMX的快速热裂解过程及机理的影响,认为在低温HMX分解的C—N键断裂在两竞争反应中占优,通过压力的变化证明了气相产物之间存在二次反应.     

N-H…O红移氢键和蓝移氢键的理论研究

用理论方法研究了复合物HCHO…HNO, HCOOH…HNO, HCHO…NH₃, HCOOH…NH₃, HCHO…NH₂F和HCOOH…NH₂F分子间氢键. 在MP2/6-31+G(d,p), MP2/6-311++G(d,p), B3LYP/ 6-31+G(d,p)和B3LYP/6-311++G(d,p)水平上, 利用标准方法和均衡校正方法对6种复合物进行了几何优化和振动频率计算. 计算结果表明: 在复合物HCHO…HNO和HCOOH…HNO中, HNO中N—H键强烈收缩, 存在显著的N—H…O蓝移氢键. 在复合物HCHO…NH₃, HCOOH…NH₃, HCHO…NH₂F和HCOOH…NH₂F中, NH3和NH₂F中N—H键伸长, 存在N—H…O红移氢键. 利用自然键轨道(NBO)分析表明, 电子供体轨道和电子受体轨道之间相互作用的稳定化能、电子密度重排、轨道再杂化和结构重组是决定氢键红移和蓝移的主要因素. 其中, 轨道间稳定化能属于键伸长效应, 电子密度重排、轨道再杂化和结构重组属于键收缩效应. 在复合物HCHO…HNO和HCOOH…HNO中, 由于键收缩效应处于优势地位导致N—H…O蓝移氢键存在. 在复合物HCHO…NH₃, HCOOH…NH₃, HCHO…NH₂F和HCOOH…NH₂F中, 键伸长效应居于主导地位,N—H…O红移氢键出现.     

高压下β-HMX热分解机理的ReaxFF反应分子动力学模拟

采用ReaxFF反应分子动力学方法研究了不同压缩态β-HMX晶体(ρ=1.89、2.11、2.22、2.46、2.80、3.20 g·cm^-3)在T=2500 K时的热分解机理, 分析了压力对初级和次级化学反应速率的影响、高压与低压下初始分解机理的区别以及造成反应机理发生变化的原因. 发现HMX的初始分解机理与压力(或密度)相关. 低压下(ρ<2.80 g·cm^-3)以分子内反应为主, 即N－NO₂键的断裂、HONO的生成以及分子主环的断裂(C－N键的断裂). 高压下(ρ≥2.80 g·cm^-3)分子内反应被显著地抑制, 而分子间反应得到促进, 生成了较多的O₂、HO等小分子和大分子团簇. 初始分解机理随压力的变化导致不同密度下的反应速率和势能也有所不同. 本文在原子水平对高压下HMX反应机理的深入研究对于认识含能材料在极端条件下的起爆、化学反应的发展以及爆轰等具有重要意义.     

高能混合物的感度理论判别—— 不同配比和不同温度AP/HMX的MD研究

为了寻求高能复合材料的感度理论判据, 对高氯酸铵(AP)和HMX(环四甲撑四硝胺)所构成的不同质量比的二元混合体系, 用分子动力学(MD)方法和修正的PCFF力场, 作正则系综(NVT)下的周期性模拟计算, 求得其结合能和HMX热解引发键(N—NO2)的平均键长和最大键长. 结果表明, 结合能随质量比不同呈复杂变化趋势; HMX(N—NO2)引发键的最大键长随体系中HMX配比增加先增后减, 而当AP/HMX为1∶1时其值最大, 恰与此配比下感度最大的实验事实相一致. 选择质量比为1∶1的AP/HMX作不同温度下的NVT-MD模拟, 发现引发键(N—NO2)最大键长随温度升高而递增, 与感度随温度升高而增大的实验事实相符. 为此我们建议, 把高能复合材料中易爆燃组分引发键的最大键长作为其热和撞击等感度的理论判据, 借以阐明、比较或预示它们的相对安全性.     

吡啶-BH~3相互作用复合物的理论研究

对吡啶-BH~3复合物分别用MP2/6-31+G^*和B3LYP/6-31+G^*进行理论计算以预测该复合物的构型及解离能，得到四种构型，在MP2优化构型基础上作CCSD/6-31+G^*单点能量计算以验证MP2与B3LYP结果的可靠性，然后用B3LYP作振动频率分析，计算了各构型的垂直电离势，最后用更大基组作单点能量计算和自然键轨道(NBO)分析。结果表明，N-B直接相连的构型最稳定，其解离能为141.50kJ/mol,MP2和B3LYP对N-H接近的构型结果相关较大，另外两种构型稳定性介于二者之间，解离能分别为15.18kJ/mol,14.06kJ/mol(MP2/6-31+G^*)。     

环四甲撑四硝胺(HMX)结构和性质的DFT研究

用密度泛函理论(DFT)B3LYP方法,取6-31G基组,求得环四甲撑四硝胺分子的几何构型、电子结构、 IR谱和298～1200 K的热力学性质.全优化几何构型和电子结构均具有Ci对称性.在相邻原子之间以N-NO2键的Mulliken集居数最小,表明其间电子分布较少,预示其为热解和起爆的引发键.IR谱与实验结果良好相符.     

Identification of octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) pyrolysis products by simultaneous thermogravimetric modulated beam mass spectrometry and time-of-flight velocity-spectra measurements

The pyrolysis products formed during the isothermal decomposition of HMX at 211°C are H₂O, HCN, CO, CH₂O, NO, N₂O, methylformamide, C₂H₆N₂O, octahydro-1-nitroso-3,5,7-trinitro-1,3,5,7-tetrazocine, and a nonvolatile residue. The temporal behaviors of these products during the decomposition are presented. The method for using time-of-flight (TOF) velocity spectra to assist mass-spectrometry measurements in identifying the different gaseous products formed from the pyrolysis of a material by determining the approximate molecular weights of the different gaseous products contributing to the different m/z values in the mass spectrum of the mixture is described. The ion fragmentation of HMX as a function of electron energy shows complete fragmentation of the HMX molecular ion for electron energies ≥ 12.4 eV. No fragments from the pyrolysis of HMX other than those mentioned above are observed.     

Decomposition of nitramine energetic materials in excited electronic states: RDX and HMX

Ultraviolet excitation (8-ns duration) is employed to study the decomposition of RDX (1,3,5-trinitro-1,3,5-triazacyclohexane) and HMX (1,3,5,7-tetranitro-1,3,5,7-tetrazacyclooctane) from their first excited electronic states. Isolated RDX and HMX are generated in the gas phase utilizing a combination of matrix-assisted laser desorption and supersonic jet expansion techniques. The NO molecule is observed as one of the initial dissociation products by both time-of-flight mass spectroscopy and laser-induced fluorescence spectroscopy. Four different vibronic transitions of NO are observed: A (2)Sigma(v(') = 0)<--X (2)Pi(v(") = 0,1,2,3). Simulations of the NO rovibronic intensities for the A<--X transitions show that dissociated NO from RDX and HMX is rotationally cold (approximately 20 K) and vibrationally hot (approximately 1800 K). Another potential initial product of RDX and HMX excited state dissociation could be OH, generated along with NO, perhaps from a HONO intermediate species. The OH radical is not observed in fluorescence even though its transition intensity is calculated to be 1.5 times that found for NO per radical generated. The HONO intermediate is thereby found not to be an important pathway for the excited electronic state decomposition of these cyclic nitramines.     

First-principles study of the four polymorphs of crystalline octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine

The electronic structure and vibrational properties of the four polymorphs of crystalline octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) have been studied using density functional theory within the local density approximation. The results show that the states of N in the ring make more important contributions to the valence bands than these of C and N of NO2 and so N in the ring acts as an active center. From the low frequency to high-frequency region, the molecular motions of the vibrational frequencies for the four HMX polymorphs present unique features. It is also noted that there is a relationship between the band gap and impact sensitivity for the four HMX polymorphs. From the cell bond order per unit volume, we may infer the variation order of crystal bonding for the four polymorphs and so predict their impact sensitivity order as follows: beta-HMX < gamma-HMX < alpha-HMX < delta-HMX, which is in agreement with their experimental order.     

Liquid chromatography/electrospray ionization tandem mass spectrometry analysis of octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX)

A quantitative liquid chromatography/electrospray ionization tandem mass spectrometry (LC/ESI-MS/MS) method was developed for the analysis of the explosive, octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX). In negative ionization mode, HMX forms an acetate adduct ion [M + CH(3)COO](-), m/z 355, in the presence of a small amount of acetic acid in the mobile phase. The ESI collision-induced dissociation (CID) spectrum of m/z 355 was acquired and the transitions m/z 355 --> 147 and m/z 355 --> 174 were chosen for the determination of HMX in samples. Using this quantification technique, the method detection limit was 1.57 microg/L and good linearity was achieved in the range 5-500 microg/L. This method will help to unambiguously analyze environmentally relevant concentrations of HMX.     

Intermolecular interactions,thermodynamic properties,crystal structure,and detonation performance of HMX/NTO cocrystal explosive

Previous studies have shown that the design of cocrystal explosives is one of the most promising approaches to decrease the sensitivity and maintain the detonation performance of existing explosives. As is well‐known, octahydro‐1,3,5,7‐tetranitro‐1,3,5,7‐tetrazocine (HMX) is a high energy density material (HEDM). But the application of HMX is limited, due to its high sensitivity. Thus, an insensitive explosive 5‐nitro‐1,2,4‐triazol‐3‐one (NTO) is proposed as a cocrystal former (CCF) to cocrystallize with HMX in the present work. The binding energies, heat of formations (HOFs), thermodynamic properties, atoms in molecules, and natural bond orbital analysis of four HMX/NTO complexes have been calculated using density functional theory methods, including meta‐hybrid functional (M062X) and dispersion‐corrected density functionals (B97D, ωB97XD). In addition, the crystal structure of HMX/NTO cocrystal has been investigated using Monte Carlo simulation and first principles methods. The HMX/NTO cocrystal is most likely to crystallize in triclinic crystal system with P1 space group, and corresponding cell parameters are Z = 2, a = 9.06 Å, b = 8.19 Å, c = 10.27 Å, α = 81.94^°, β = 98.42^°, γ = 82.03^°, and ρ = 1.92 g/cm³. The detonation velocity and detonation pressure of HMX/NTO cocrystal are 8.73 km/s and 35.14 GPa, respectively, a little lower than those of HMX. Finally, bond dissociation energies (BDEs) of the weakest trigger bond in HMX/NTO complexes are calculated. The results show that HMX/NTO complexes are thermally stable and meet the thermal requirement of HEDMs (BDE > 120 kJ/mol). © 2012 Wiley Periodicals, Inc.     

Initial chemical events in shocked octahydro-1,3,5,7-tetranitro-1,3,5,7- tetrazocine: a new initiation decomposition mechanism

We have performed ab initio molecular dynamics simulations in conjunction with the multiscale shock technique to study the initial chemical processes of octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) under shock wave loading. The results show that the initial decomposition of shocked HMX is triggered by the N-O bond breaking and the ring opening. This will initiate many decomposition reactions and lead to the production of many small radicals at a moment. As the shock compression continues, these small radicals recombine to produce many large radicals and further form ring-shaped radicals. Then, these radicals begin to further decompose. It is also found that the system transiently produces a large number of metallic states under the shock compression. Our simulations thus suggest a new mechanism for the initial chemical processes of shocked HMX and provide fundamental insight into the initial mechanism at the atomistic level, which is of important implication for understanding and development of energetic materials.