TATB基PBX结合能的分子动力学模拟

用分子动力学（MD）方法, 模拟计算了四种氟聚合物(聚偏二氟乙烯（PVDF）、聚三氟氯乙烯（PCTFE）、氟橡胶（F₂₃₁₁)、氟树脂(F₂₃₁₄))与TATB(1,3,5- 三氨基- 2,4,6- 三硝基苯)晶体的相互作用. 结果发现, 四种氟聚物与TATB的结合能大小排序为PVDF>F₂₃₁₁>F₂₃₁₄>PCTFE, 各氟聚物在TATB不同晶面上的结合能大小排序为(001)>(010)>(100), 结合能主要由分子间氢键决定.     

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时, 材料的刚性减弱, 柔性增强.     

碳F4F6多面体遵循独立四元环原则

为了研究由四元环和六元环构成的碳多面体(F₄F₆多面体)的结构和稳定性之间的关系, 本文采用密度泛函方法对所有C₈~C₆₀之间的F₄F₆多面体进行了系统的计算研究. 结果表明, 能量最低的异构体都满足独立四元环原则, 能量较低的满足四元环比邻惩罚原则. 这两条原则与经典富勒烯所遵循的独立五元环原则和五元环比邻惩罚原则具有同等地位, 能够使研究者仅仅从形貌上就可以对碳F₄F₆多面体的稳定性进行简单高效的判定. 结构分析表明, 四元环之间共用的顶点的锥化角大于其他顶点的锥化角并决定了相应分子的稳定性.     

一个新颖的三维异核锑(III)-镨(III)配合物[Sb2-μ4-(EDTA)2Pr(H2O)5]NO3•4H2O的合成、晶体结构及热分解研究

合成了锑-镨与乙二胺四乙酸形成的新颖三维配合物[Sb₂-μ₄-(EDTA)₂Pr(H₂O)₅]NO₃•4H₂O, 用元素分析、红外光谱、热分析及X射线单晶衍射法等进行了组成和结构表征. 结果表明该配合物属正交晶系, 空间群Pnn2; 晶胞参数: a＝1.07031(2) nm, b＝2.30805(4) nm, c＝0.72343(2) nm, V＝1.78711(7) nm³, Z＝2, D_c＝2.202 g/cm³, F(000)＝1164, μ＝2.955 mm^－1, GOF＝1.000, 最终偏离因子R₁＝0.0203, wR₂＝0.0545 [I＞2σ(I)]. 在标题化合物中, 每个镨(III)离子的配位数为9, 与五个水分子中的五个氧原子和四个羧基氧原子配位, 形成三帽三角棱柱空间配位多面体. 锑(III)与EDTA离子中的四个氧原子和两个氮原子配位, 在赤道平面上有一孤对电子. 同时讨论了配合物的热分解过程.     

HMX/F2311 NPT高聚物粘合炸药力学性能的NPT系综分子动力学模拟研究

本文用分子动力学模拟方法研究了以β型奥克托金(β-HMX)含能材料为基,以氟聚物F2311为粘合剂的高聚物粘合炸药(PBX)的力学性能。β-HMX晶体和以其为基的高聚物粘合炸药的弹性常数均由静态弹性常数分析法计算求得，而工程模量和泊松比则由Reuss平均法导出。根据柯西压和本体模量与剪切模量的比值，发现通过加入少量该种聚合物可有效地提高HMX晶体的延展性。     

TATB/氟聚物PBX沿不同晶面力学性能的分子动力学模拟

用分子动力学(MD)方法, 模拟计算了著名钝感炸药TATB(1, 3, 5-三氨基-2, 4, 6-三硝基苯)与四种氟聚合物[聚偏二氟乙烯(PVDF)、聚三氟氯乙烯(PCTFE)、氟橡胶(F₂₃₁₁)、氟树脂(F₂₃₁₄)]构成的高聚物粘结炸药(PBX)的力学性能. 结果表明, 在TATB中添加少量氟聚物能有效改善其力学性能; 沿TATB不同晶面与氟聚物“粘结”, 构成PBX的力学性能有所不同, 改善力学性能的整体效应为(010)≈(100)>(001).     

不同温度下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种晶体的弹性力学性能, 发现拉伸、 体积和剪切模量均随温度的升高而递减, 与实验结果一致.     

5-[N-乙酸根(4-吡啶基)]四唑为桥联配体的铜(Ⅱ)、镉(Ⅱ)配聚物

以5-[N-乙酸根(4-吡啶基)]四唑为桥联配体，分别与Cu(ClO)₂·6H₂O和CdCl₂的水溶液反应，获得配聚物[Cu(a4-ptz)₂·2H₂O]_n (1)和[{Cd(a4-ptz)·(H₂O)₂Cl}·H₂O]_n (2)，2个化合物通过元素分析、IR、TGA等表征，并测定了它的晶体结构。结果表明：化合物1是一个二维四方格子结构的配位聚合体，化合物2含有四元环(Cd₂Cl₂)的一维Z-型聚合链。化合物2在280 ℃以下热稳定性好，并且不溶于一般溶剂，因此可成为潜在的荧光材料。     

高氯酸稀土Sm(Ⅲ)开环冠醚(EO4)配合物的晶体结构

The crystal of title compounds belongs to orthorhombic system with space group P2₁2₁2₁ (No.19). The cell parmeters are a=1.432 16(3) nm, b=1.457 20(3) nm， c=3.005 85(6) nm, and V=6 273.0(2) nm³, Z=8, μ(Mo Kα)=22.21 cm^-1， D_c=1.802 g·cm^-3, F(000)=3 416.00, Finally R₁=0.033 and wR₂=0.073 (I>2.00σ(I)), S=1.023. Every center atom Sm(Ⅲ) is surrounded by the two aether chains, the title compound has ten-coordinate structure and the coordination polyhedron shows a Bicapped Square Antiprism. In the crystal, the title compound has the flack of double moleculars and non-symmetry. The anion ClO₄^- does not take part in the coordination. CCDC: 238709.     

配合物Zn(Met)SO4•H2O(s)的低温热容和标准摩尔生成焓

利用精密自动绝热热量计直接测定了配合物Zn(Met)SO₄•H₂O(s) 在78～370 K温区的摩尔热容. 通过热容曲线的解析得到该配合物的起始脱水温度为T₀＝329.50 K. 将该温区的摩尔热容实验值用最小二乘法拟合得到摩尔热容 (C_p,m)对温度(T)的多项式方程, 并且在此基础上计算出了它的舒平热容值和各种热力学函数值. 依据Hess定律, 通过设计热化学循环, 选择体积为100 cm³、浓度为2 mol•L^－1的盐酸作为量热溶剂, 利用等温环境溶解-反应热量计, 测定和推算出该配合物的标准摩尔生成焓为Δ_fH_m⁰＝－(2069.30±0.74) kJ•mol^－1.     

HMX和HMX/HTPB PBX的晶体缺陷理论研究

建立空位和掺杂点缺陷模型, 用分子动力学(MD)方法, 研究晶体缺陷对β-环四亚甲基硝胺(HMX)和β-HMX/HTPB(端羟基聚丁二烯)高聚物粘结炸药(PBX)的力学性能和爆炸性能的影响. 结果表明, 相对于HMX“完美”晶体(1)考察缺陷晶体(2和3), 以及相对于HMX完美晶体基PBX(1)考察缺陷PBX 2和PBX 3, 均发现弹性系数和(拉伸、体积、剪切)模量下降, 导致体系刚性减弱, 延展性和韧性增强. 这与在基炸药HMX晶体(1, 2和3)中分别加入HTPB高聚物粘结剂形成PBX 1, PBX 2和PBX 3呈现类似的相应的变化趋势和效果. 此外, 研究表明, 爆炸性质也依赖于体系的组成和结构. 因加入的是低能高聚物, 故PBX(1), PBX(2)和PBX(3)的爆热、爆速和爆压均比相应的基炸药(1, 2和3)低, 即晶体(1)＞PBX(1), 晶体(2)＞PBX(2), 晶体(3)＞PBX(3). PBX(1), PBX(2), PBX(3)与对应基炸药(1, 2, 3)的爆速和爆压取相同变化次序, 亦即PBX(1)＞PBX(2)＞PBX(3)对应于晶体(1)＞晶体(2)＞晶体(3). 这些计算结果和规律对PBX配方设计显然具有指导作用.     

TATB/氟聚物PBX力学性能、结合能和爆炸性能的模拟计算

Molecular dynamics （MD） method was used to simulate 1,3,5-triamino-2,4,6-trinitrobenzene （TATB） coated with fluorine containing polymers. The mechanical properties and binding energies of PBXs were obtained. It was found that when the number of chain monomers of fluorine containing polymers was the same, the elasticity of TATB/F2314 was increased more greatly than others and the binding energy of TATB/F2311 was the largest among four PBXs. Detonation heat and velocity of such four PBXs were calculated according to theoretical and empirical formulas. The results show that the order of detonation heat is TATB〉TATB/PVDF〉TATB/F2311〉TATB/ F2314 〉 TATB/PCTFE while the order of detonation velocity is TATB/PVDF 〈 TATB/F2311 〈 TATB/F2314 〈 TATB/PCTFE 〈TATB.     

Anisotropy of Thermal-expansion for β-Octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine: Quantum Chemistry Calculation and Molecular Dynamics Simulation

Molecular dynamics simulations on octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) at 303-383 K and atmospheric pressure are carried out under NPT ensemble and COMPASS force field, the equilibrium structures at elevated temperatures were obtained and showed that the stacking style of molecules don't change. The coefficient of thermal expansion (CTE) values were calculated by linear fitting method. The results show that the CTE values are close to the experimental results and show anisotropy. The total energies of HMX cells with separately increasing expansion rates (100%-105%) along each crystallographic axis was calculated by periodic density functional theory method, the results of the energy change rates are anisotropic, and the correlation equations of energy change-CTE values are established. Thus the hypostasis of the anisotropy of HMX crystal's thermal expansion, the determinate molecular packing style, is elucidated.     

β-1,3,5,7-四硝基-1,3,5,7-四氮杂环辛烷热膨胀各向异性的量子化学计算与分子动力学模拟

在303～383 K和NPT系综和COMPASS力场下对β-1,3,5,7-四硝基-1,3,5,7-四氮杂环辛烷（HMX）超晶胞初始结构的分子动力学模拟,得到常压下各温度的晶体平衡构型并发现分子的堆积方式不变;通过线性拟合求算出线膨胀系数与实验值相近,体现出明显的各向异性. 采用密度泛函理论方法对沿各晶轴方向膨胀率变化（100%～105%）的HMX单胞模型进行了总能计算,得到的能量变化率体现各向异性并与热膨胀系数值关联,建立了关联方程. 由此阐释了HMX晶体热膨胀各向异性的本质即特定的分子堆积模式.     

Molecular Dynamics Simulations for Effects of Fluoropolymer Binder Content in CL-20/TNT Based Polymer-Bonded Explosives

In order to better understand the role of binder content, molecular dynamics (MD) simulations were performed to study the interfacial interactions, sensitivity and mechanical properties of 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane/2,4,6-trinitrotoluene (CL-20/TNT) based polymer-bonded explosives (PBXs) with fluorine rubber F₂₃₁₁. The binding energy between CL-20/TNT co-crystal (1 0 0) surface and F₂₃₁₁, pair correlation function, the maximum bond length of the N–NO₂ trigger bond, and the mechanical properties of the PBXs were reported. From the calculated binding energy, it was found that binding energy increases with increasing F₂₃₁₁ content. Additionally, according to the results of pair correlation function, it turns out that H–O hydrogen bonds and H–F hydrogen bonds exist between F₂₃₁₁ molecules and the molecules in CL-20/TNT. The length of trigger bond in CL-20/TNT were adopted as theoretical criterion of sensitivity. The maximum bond length of the N–NO₂ trigger bond decreased very significantly when the F₂₃₁₁ content increased from 0 to 9.2%. This indicated increasing F₂₃₁₁ content can reduce sensitivity and improve thermal stability. However, the maximum bond length of the N–NO₂ trigger bond remained essentially unchanged when the F₂₃₁₁ content was further increased. Additionally, the calculated mechanical data indicated that with the increase in F₂₃₁₁ content, the rigidity of CL-20/TNT based PBXs was decrease, the toughness was improved.     

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.     

Molecular dynamics simulations of trans-1,4,5,8-tetranitro-1,4,5,8-tetraazadecalin-based polymer-bonded explosives

Molecular dynamics has been applied to investigate the low-sensitivity explosive TNAD (trans-1,4,5,8-tetranitro-1,4,5,8-tetraazadecalin)-based polymer-bonded explosives (PBXs) with four typical fluorine polymers, PVDF (polyvinylidenedifluoride), PCTFE (polychlorotrifluoroethylene), F(2311) (fluorine rubber), and F(2314) (fluorine resin). The elastic constants, mechanical properties (tensile modulus, bulk modulus, shear modulus, and Poission ratio), binding energies, and detonation performances are first reported for the TNAD-based PBXs. The results show that the mechanical properties of TNAD can be effectively improved by the addition of small amounts of fluorine polymers, and the overall effect of fluorine polymers on the mechanical properties of the PBXs along three crystalline surfaces is (001) > (010) > (100). On each crystal surface, improvement in the ductibility made by the fluorine polymers changes approximately in the sequence of PVDF > F(2311) > F(2314) > PCTFE. The binding energies between different TNAD crystalline surfaces and different polymer binders with the same chain segment or mass fraction both decrease in the order of (010) > (100) > (001). The binding properties of the polymers with the same chain segment on each crystal surface of TNAD increase as PVDF < F(2311) < F(2314) < PCTFE, while those of different polymers in the same content decrease in the sequence of PVDF > F(2311) > F(2314) > PCTFE. The detonation performances of the PBXs decrease in comparison with the pure crystal, but they are superior to those of TNT.     

高能混合物的感度理论判别—— 不同配比和不同温度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)最大键长随温度升高而递增, 与感度随温度升高而增大的实验事实相符. 为此我们建议, 把高能复合材料中易爆燃组分引发键的最大键长作为其热和撞击等感度的理论判据, 借以阐明、比较或预示它们的相对安全性.     

Molecular dynamics simulations for pure epsilon-CL-20 and epsilon-CL-20-based PBXs

Molecular dynamics has been employed to simulate the well-known high energy density compound epsilon-CL-20 (hexanitrohexaazaisowurtzitane) crystal and 12 epsilon-CL-20-based PBXs (polymer bonded explosives) with four kinds of typical fluorine polymers, i.e., polyvinylidenedifluoride, polychlorotrifluoroethylene, fluorine rubber (F(2311)), and fluorine resin (F(2314)) individually. The elastic coefficients, isotropic mechanical properties (tensile moduli, bulk moduli, shear moduli, and poission's ratios), and bonding energy are first reported for epsilon-CL-20 crystal and epsilon-CL-20-based polymer bonded explosives (PBXs). The mechanical properties of epsilon-CL-20 can be effectively improved by blending with a small amount of fluorine polymers, and the whole effect of the adding fluorine polymers to improve mechanical properties of PBXs along the three crystalline surfaces of epsilon-CL-20 is found to be (100) approximately (001) > (010). The interaction between each of the crystalline surfaces and each of the fluorine polymers is different, and the ordering of binding energy for the three surfaces is (001) > (100) > (010); F(2314) always has the strongest binding ability with the three different surfaces. F(2314) can best improve the ductibility and tenacity of PBX when it is positioned on epsilon-CL-20 (001) crystal surface. The calculations on detonation performances for pure epsilon-CL-20 crystal and the four epsilon-CL-20-based PBXs show that adding a small amount of fluorine polymer into pure epsilon-CL-20 will lower detonation performance, but each detonation parameter of the obtained PBXs is still excellent.