RDX基PBX的模型、结构、能量及其与感度关系的分子动力学研究

以RDX（环三亚甲基三硝胺）为基、PS（聚苯乙烯）为粘结剂构成PBX（高聚物粘结炸药）的MD（分子动力学）模拟初始模型.比较分别以1根46链节和2根23链节PS置于RDX（001）晶面上的两种（PBX1和PBX2）模型下的MD模拟结果,发现二者的结构、相互作用能和力学性能均很接近.取PBX2进行5种温度（195,245,295,345和395 K）下的NPT系综、MD模拟系统研究,发现随温度依次升高,各体系中RDX引发键N NO2键的最大键长（Lmax）递增,N–N键连的N与N之间的双原子作用能（EN-N）和内聚能密度（CED）递减,与感度随温度升高而增大的实验事实相一致.综合已有工作,对高能复合材料（如PBX和固体推进剂等）的感度理论研究,建议关注其中易爆燃组分在外界刺激下的结构和能量变化,其引发键Lmax和作为引发键强度度量的双原子作用能（如EN-N）,可作为热和撞击感度相对大小的理论判据.     

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的分解温度大大降低.     

不同温度下HMX和RDX晶体的感度判别和力学性能预估分子动力学比较研究

为探讨和比较炸药黑索金(RDX)和奥克托今(HMX)晶体的结构、 能量和力学性质随温度的递变规律, 在COMPASS力场和NPT系综下, 对其合适的等原子数超晶胞模型分别进行5个温度(195, 245, 295, 345和395 K)下的分子动力学(MD)周期性模拟研究. 结果表明, 随着温度的升高, RDX和HMX晶体引发键N-NO₂的最大键长(L_max)的逐渐增大以及引发键连双原子作用能(E_N-N)和内聚能密度(CED)的逐渐减小均与感度随温度升高而增大的实验结果一致; 且在各温度下, RDX晶体的L_max均大于HMX晶体的L_max, 与HMX相比, RDX的E_N-N和CED均较小, 上述结果与RDX比HMX感度大的实验结果相符. 由此表明, 在一定条件下, L_max, E_N-N和CED可用于高能物质的热和撞击感度的相对大小的判别. 基于MD模拟原子运动轨迹, 用静态法求得2种晶体的弹性力学性能, 发现拉伸、 体积和剪切模量均随温度的升高而递减, 与实验结果一致.     

外电场对不同组分比例的HMX/FOX-7共晶炸药感度影响的理论研究

借助分子动力学方法对不同电场不同组分比例的HMX/FOX-7共晶炸药晶体表面成键能和力学性质进行了研究。在不同外电场作用下,利用B3LYP和MP2(full)方法在6-311++G(d,p)和6-311++G(2df,2p)基组水平上,计算了HMX/FOX-7(1∶1)复合物中HMX的撞击感度h_(50)及其N-NO_2键离解能。结果表明:1∶1比例的共晶炸药最稳定、力学性质较好。共晶主要在HMX(0 2 0)和(1 0 0)晶面上形成。共晶复合物的形成使HMX引发键N-NO_2增强,感度降低。正方向外电场增大了晶面成键能和HMX中N-NO_2键离解能,升高了h_(50)值,降低了HMX的撞击感度;而负方向外电场对其影响正好相反。正方向外电场使G和E值增大,炸药延展性降低;负方向外电场效应正好相反。     

NO2, 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降低近200 kJ·mol-1, 表明OH-对其裂解有明显的促进作用. NO2、OH-的存在可使HMX的分解温度大大降低.     

HMX/TATB复合材料弹性性能的MD模拟

用分子动力学(MD)方法COMPASS力场, 分别在正则系综(NVT)和等温等压系综(NPT)下, 模拟计算了著名常用高能炸药HMX(环四甲撑四硝胺)与著名钝感炸药TATB (1,3,5-三氨基-2,4,6三硝基苯)所构成的混合体系在室温时的弹性性能和结合能. 结果表明, 在NVT和NPT两种系综下模拟所得结果呈平行一致的趋势; 与纯HMX相比, HMX/TATB复合材料的拉伸模量、体模量和剪切模量均有所下降; 在NVT系综下, 还完成了HMX/TATB混合体系的不同温度的MD模拟. 发现当温度在245～345 K范围时, 体系的刚性和弹性变化很小; 但当温度达到395 K时, 材料的刚性减弱, 柔性增强.     

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条合成路线,均具有步骤少且原料易得的优点.     

胞嘧啶…NO复合物结构与性质的理论研究

用密度泛函理论在BL3YP/6-311＋G*基组水平上对胞嘧啶…NO复合物体系进行了理论计算, 发现了6个能量极小的复合物. 其结合方式是NO的N或O原子与胞嘧啶的N—H键形成氢键, 最稳定的复合物的结合能为－9.65 kJ/mol. NO的N原子与胞嘧啶的结合具有更强的优势, N结合的复合物中NO的键长缩短, 而O结合的复合物中NO键长伸长. 同时, 对复合物的振动分析发现, 在胞嘧啶中所有的与NO结合的N—H键的伸缩频率下降, 而所有氨基的面内弯曲振动频率是增加的.     

高能体系分子间相互作用研究: 含NNO~2和NH~2混合物

以abinitioHF/6-31G^*计算求得NH~3+NH~2NO~2的两种优化构型,经MP4电子相关能校正和Boys-Bernardi方案校正基组叠加误差求得精确的分子间相互作用能。还用PM3方法计算研究TATB(均三氨基三硝基苯)分别与HMX(奥克托金)和RDX(黑索金)的混合体系,经色散能校正电子相关近似地求得分子间相互作用能。结果表明,NH~3与NH~2NO~2之间的最大结合能为-38.32kJ/mol;分子间相互作用增强了N-NO~2键强度;TATB与HMX,RDX的结合能远大于石墨与HMX或RDX的结合能,表明TATB对HMX和RDX的润湿和钝感作用较石墨更强。     

TATB及其衍生物的撞击感度与C-NO2键级的关系

对TATB及其系列衍生物进行了半经验分子轨道PM3和AM1计算，求得全优化几何构型和电子结构．结果发现各化合物分子中C-NO2键的Wiberg键级(WCN)均小于其他C—N或C—C键的键级．实验撞击感度h50%与WCN计算结果之间分别呈三阶成二阶多项式关系，相关系数r≥0．92．当WCN＜0．90时，该值对撞击感度影响甚微，表明属高感度炸药；而当WCN＞0．90时，撞击感度与分子中最小的G—NO2键级WCN呈显的线性关系，WCN增大，撞击感度线性递减，符合判别感度相对大小的“最小键级原理”．     

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.     

Molecular dynamics study on the relationships of modeling, structural and energy properties with sensitivity for RDX-based PBXs

In this paper,a primary model is established for MD(molecular dynamics) simulation for the PBXs(polymer-bonded explosives) with RDX(cyclotrimethylene trinitramine) as base explosive and PS as polymer binder.A series of results from the MD simulation are compared between two PBX models,which are represented by PBX1 and PBX2,respectively,including one PS molecular chain having 46 repeating units and two PS molecular chains with each having 23 repeating units.It has been found that their structural,interaction energy and mechanical properties are basically consistent between the two models.A systematic MD study for the PBX2 is performed under NPT conditions at five different temperatures,i.e.,195 K,245 K,295 K,345 K,and 395 K.We have found that with the temperature increase,the maximum bond length(L max) of RDX N N trigger bond increases,and the interaction energy(E N-N) between two N atoms of the N-N trigger bond and the cohesive energy density(CED) decrease.These phenomena agree with the experimental fact that the PBX becomes more sensitive as the temperature increases.Therefore,we propose to use the maximum bond length L max of the trigger bond of the easily decomposed and exploded component and the interaction energy E N-N of the two relevant atoms as theoretical criteria to judge or predict the relative degree of heat and impact sensitivity for the energetic composites such as PBXs and solid propellants.     

Ab initio molecular dynamics simulations study on initial decompositions of β‐HMX at high temperature coupled with high pressures

We performed ab initio molecular dynamics simulations to investigate initial decomposition mechanisms and subsequent chemical processes of β‐HMX (cyclotetramethylene tetranitramine) (octahydro‐1,3,5,7‐tetranitro‐1,3,5,7‐tetrazocine) crystals at high temperature coupled with high pressures. It was found that the initial decomposition step is the simultaneous C–H and N–NO₂ bond cleavage at 3,500 K. When the pressure (1–10 GPa) is applied, the first reaction steps are primarily the C–N and C–H bond fission at 3,500 K. The C–H bond cleavage is a triggering decomposition step of the HMX crystals at 3,500 K coupled with 16 GPa. This indicates that the C–H bonds are much easier to be broken and the hydrogen radicals are much more active. The applied pressures (1–10 GPa) accelerate the decompositions of HMX at 3,500 K. The decomposition pathways and time evolution of the main chemical species demonstrate that the temperature is the foremost factor that affects the decomposition at high pressures (1–10 GPa). However, the decomposition of HMX is dependent on both the temperature (3,500 K) and the pressure (16 GPa). This work will enrich the knowledge of the decompositions of condensed energetic materials under extreme conditions.     

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

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

1,3,5,7-四乙酰基-1,3,5,7-四氮杂环辛烷合成反应中间体的制备分离

1,3,5,7-四乙酰基-1,3,5,7-四氮杂环辛烷(TAT)是高能炸药奥克托今(HMX)重要的合成反应前体之一,1,5-二乙酰基-1,3,5,7-四氮杂环辛烷(DAPT)醋酐酰解法和乙腈甲醛小分子缩合法均能制备得到TAT.采用半制备高效液相色谱制备分离了上述两种制备方法反应原液中的六种中间体,分别是1-乙酰氧甲基-3,5,7-三乙酰基-1,3,5,7-四氮杂环辛烷、1,3,5-三乙酰基-1,3,5,7-四氮杂环辛烷、1-乙酰氧甲基-3,5,7-三乙酰基-1,3,5,7-四氮杂环辛烷、N,N’-亚甲基二乙酰胺、三乙酰胺二亚甲基三胺和四乙酰胺三亚甲基四胺,为反应机理研究提供了有力的证据.DAPT醋酐酰解法反应原液样品的最佳制备色谱条件为:正相硅胶柱(20 mm×250 mm,10～20μm),流动相:V(乙腈)∶V(甲醇)=95∶5,流速为10 mL·min-1,检测波长为215 nm,进样量为1 mL,样品浓度30 mg/mL.乙腈甲醛小分子缩合法反应原液样品的最佳制备色谱条件为:正相硅胶柱(20 mm×250 mm,10～20μm),流动相:V(乙腈)∶V(水)=90∶10,流速为10 mL·min-1,检测波长为215 nm,进样量为1 mL,样品浓度30 mg/mL.制备得到的中间体经HPLC分析知纯度均在97%以上,可直接作为结构鉴定的标准样品.     

Theoretical studies on high energetic density nitramine explosives containing pyridine

The insensitive property of explosives containing pyridine is combined with the high energy of nitramine explosives,and the concept of new nitramine explosives containing pyridine is proposed,into which nitramine group with N N bonds is introduced as much as possible.Based on molecular structures of nitramine compounds containing pyridine,density functional theory(DFT) calculation method was applied to study designed molecules at B3LYP/6-31+G(d) level.The geometric and electronic structures,density,heats of formation(HOF),detonation performance and bond dissociation energies(BDE) were investigated and comparable to 1,3,5-trinitro-1,3,5-triazinane(RDX) and 1,3,5,7-tetranitro-1,3,5,7-tetrazocane(HMX).The simulation results reveal that molecules B and D perform similarly to traditionally used RDX.Molecule E outperform RDX,with performance that approach that of HMX and may be considered as potential candidate of high energy density compound(HEDC).These results provide basic information for molecular design of novel high energetic density compounds.     

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.