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）,可作为热和撞击感度相对大小的理论判据.     

Relation between the N-NO2 bond length and stability of the secondary nitramines

The fact of the constancy of activation entropy of N-NO₂ bond homolysis in a series of secondary nitramines was utilized for correction of the experimental values of activation energy E of this process proceeding from the reliable data for the rate constants of the nitramines decomposition in solutions. When comparing the refined values of E (kJ mol^?1) with the N-N bond length d _N-N (Å) the following correlations were obtained: for cyclic and framework nitramines E = 663 ? 356d _N-N, and for the aromatic nitramines E = 1810 ? 1227d _N-N. A linear relationship between E and d is observed in the series of similar compounds. It depends on the electronic and steric effects of substituents.     

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配方设计显然具有指导作用.     

HMX晶体和HMX/F2311 PBXs力学性能的MD模拟研究

用分子动力学方法, 在295 K NVT系综和COMPASS力场下, 对环四甲撑四硝胺(HMX)晶体和F₂₃₁₁沿HMX (001), (010)和(100)晶面所构成PBXs模型的力学性能进行模拟研究. 结果表明, 加入F₂₃₁₁降低了HMX的刚性, 增强了它的延展性. 在HMX (100)面上添加F₂₃₁₁对提高体系的延展性较显著. 为考察温度对力学性能的影响及其机理, 在245～445 K范围完成对HMX (100)/F₂₃₁₁ PBX的MD模拟. 力学分析表明, 随温度增加HMX (100)/F₂₃₁₁的延展性呈抛物线变化规律, 归因于F₂₃₁₁分子链的运动及其构象随温度的变化.     

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.     

Molecular dynamics simulations on the structures and properties of ε-CL-20-based PBXs——Primary theoretical studies on HEDM formulation design

Five polymer bonded explosives (PBXs) with the base explosiveε-CL-20 (hexanitrohexaazaisowurtzitane), the most important high energy density compound (HEDC), and five polymer binders (Estane 5703, GAP, HTPB, PEG, and F2314) were constructed. Molecular dynamics (MD) method was employed to investigate their binding energies (Ebind), compatibility, safety, mechanical properties, and energetic properties. The information and rules were reported for choosing better binders and guiding formulation design of high energy density material (HEDM). According to the calculated binding energies, the ordering of compatibility and stability of the five PBXs was predicted as ε-CL-20/PEG ＞ ε-CL-20/ Estane5703 ≈ε-CL-20/GAP ＞ ε-CL-20/HTPB ＞ ε-CL-20/F2314. By pair correlation function g(r) analyses, hydrogen bonds and vdw are found to be the main interactions between the two components. The elasticity and isotropy of PBXs based ε-CL-20 can be obviously improved more than pure ε-CL-20 crystal. It is not by changing the molecular structures of ε-CL-20 for each binder to affect the sensitivity. The safety and energetic properties of these PBXs are mainly influenced by the thermal capability (C°p) and density (ρ) of binders, respectively.     

Predicting Trigger Bonds in Explosive Materials through Wiberg Bond Index Analysis

Understanding the explosive decomposition pathways of high‐energy‐density materials (HEDMs) is important for developing compounds with improved properties. Rapid reaction rates make the detonation mechanisms of HEDMs difficult to understand, so computational tools are used to predict trigger bonds—weak bonds that break, leading to detonation. Wiberg bond indices (WBIs) have been used to compare bond densities in HEDMs to reference molecules to provide a relative scale for the bond strength to predict the activated bonds most likely to break to trigger an explosion. This analysis confirms that X?NO₂ (X=N,C,O) bonds are trigger linkages in common HEDMs such as TNT, RDX and PETN, consistent with previous experimental and theoretical studies. Calculations on a small test set of substituted tetrazoles show that the assignment of the trigger bond depends upon the functionality of the material and that the relative weakening of the bond correlates with experimental impact sensitivities.     

Thermal behavior and decomposition kinetics of Formex-bonded explosives containing different cyclic nitramines

Thermal behavior and decomposition kinetics of Formex-bonded PBXs based on some attractive cyclic nitramines, such as 1,3,5-trinitro-1,3,5-triazinane (RDX) and 1,3,5,7-tetranitro-1,3,5,7-tetrazocane (HMX). Actually, cis-1,3,4,6-tetranitrooctahy droimidazo-[4,5-d]imidazole (BCHMX) and 2,4,6,8,10,12-hexanitro-2,4,6,8,10, 12-hexaazaisowurtzitane (CL-20), was investigated by means of nonisothermal thermogravimetry (TG) and differential scanning calorimetry (DSC). It was found that the mass loss rate of PBXs involved in this research depends greatly on heating rate and the residue of the decomposition of these PBXs decreases with the heating rate. The onset of the exotherms was noticed at 215.4, 278.7, 231.2 and 233.7 °C with the peak maximum at 235.1, 279.0, 231.2 and 233.7 °C for RDX-Formex, HMX-Formex, CL-20-Formex, and BCHMX-Formex, respectively. Their corresponding exothermic changes were 1788, 1237, 691, and 1583 J g^?1. It was also observed that the dependence on the heating rate for onset temperatures of HMX- and BCHMX-based PBXs was almost the same due to their similar molecular structure. In addition, based on nonisothermal TG data, the kinetic parameters for thermal decomposition of these PBXs were calculated by isoconversional methods. It was shown that the Formex base has great effects on the activation energy distribution of nitramines. It was further found that the kinetic compensation effects occurred during the thermal decomposition of nitramine-based PBXs, and they almost have the same compensation effects due to similar decomposition mechanism.     

Thermodynamique de composes azotes. VII. Etude thermochimique du pyrazole et de l'imidazole

Enthalpies of sublimation for pyrazole and imidazole have been obtained by calorimetry at 298.15K. The ΔH⁰_sub (298.15 K) values for these two compounds are, respectively, 69.16 ± 0.32 and 74.50 ± 0.40 kJ mole^?1. From literature data obtained by combustion calorimetry for ΔH⁰_f (c, 298.15 K), the enthalpies of formation of these compounds in the gaseous state (pyrazole: 185.1 ± 2.3 kJ mole^?, imidazole: 133.0 ± 1.7 kJ mole^?1) have been derived. Several energy values related to the molecular structure of these two compounds (as resonance energy, enthalpy of isomerization, …) have been determined. The study of pyrazole has enabled us to contribute to the evaluation of some characteristics of the NN bond.     

Notes on the use of the vacuum stability test in the study of initiation reactivity of attractive cyclic nitramines in the C4 matrix

The zero-order reaction rates (specific rate constants) for isothermal decomposition at 120 °C of plastic bonded explosives (PBXs) were measured by means of the Czech vacuum stability test, STABIL. The PBXs are based on 1,3,5-trinitro-1,3,5-triazinane (RDX), β-1,3,5,7-tetranitro-1,3,5,7-tetrazocane (HMX), cis-1,3,4,6-tetranitro-octahydroimidazo-[4,5-d] imidazole (BCHMX), and ε-2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (ε-HNIW, ε-CL-20) with 9 wt.% of the C4 type matrix, i.e., a matrix containing 25 wt.% of polyisobutylene, 59 wt.% of dioctyl sebacate and 16 wt.% of the oil HM46. Dependencies were found between the specific rate constants mentioned and the detonation velocities of the PBXs, and consequently between these constants and the impact of pure explosive fillers, i.e., RDX, β-HMX, ε-HNIW, RS-ε-HNIW, and BCHMX and, at the same time, the corresponding PBXs. The results obtained are compared with those from a recent analogous study of PBXs using an SBR (Formex P1) binder which could increase the PBXs’ reactivity in comparison with C4 matrix. These results also help to dispel a widely held view about HNIW being a relatively sensitive explosive.     

Large-Scale Molecular Dynamics Simulations of Energetic Ni Nanocluster Impact onto the Surface

On the atomic scale, Molecular Dynamics (MD) Simulation of Nano Ni cluster impact on Ni (100) substrate surface have been carried out for energies of E _a = 1–5 eV/atom and total energy of E _T = 195 eV (the total energy of cluster is E _T = nE _a, n is the number of cluster atoms) to understand quantitatively the interaction mechanisms between the cluster atoms and the substrate atoms. The many-body Embedded Atom Method (EAM) was used in this simulation. We investigated the maximum substrate temperature T _max and the time t _max within which this temperature is reached as a function of cluster sizes and the total energy E _T. The temperature T _max is linearly proportional to total cluster energy. For the constant energy per atom and for the cluster size increase, the correlated collisions rapidly transfers energy to the substrate, and the time t _max approached a constant value. For constant total energy the temperature T _max and the time t _max versus different cluster sizes was studied. We showed that the cluster implantation and sputtering atoms from the surface are affected by the cluster size and total kinetic energy of the clusters. Finally time dependence of the number N _dis of disordered atoms in the substrate was observed.     

Spectroscopic and electrochemical properties of cyclopalladated complexes derived from 2,3-diphenylpyrazine

A comparative study of complexes [Pd(dphpz)(N∧N)]PF₆ [dphpz^? is the deprotonated form of 2,3-diphenylpyrazine; (N∧N) is ethylenediamine (En), 2,2′-bipyridine (bpy), o-phenanthroline (phen), dipyrido[a,c]phenazine (dppz), 6,7-dicyanopyrido[f,h]quinoxaline (dicnq)] was made, using ¹H NMR, electronic absorption, and emission spectroscopy, and also cyclic voltammetry. Steric interaction of the dphpz^? phenyl rings leads to significant proton shielding in the carbanionic moiety of the cyclometallated ligand. Introduction of heterocyclic diimines instead of ethylenediamine decreases the desheilding of the dphpz^? protons adjacent to the coordination center. Irrespective of the nature of the N∧N ligands, the cyclopalladated complexes are characterized by specific parameters of photo-and electrostimulated electron transfer processes involving the Pd(dphpz) orbitals, namely, by the long-wave absorption band with λ_max 395±6 nm and ε (2.2±1.2) × 10³ 1 mol^?1 cm^?1, the vibrationally structured low-temperature (77 K) luminescence resulting from the spinforbidden optical transfer from the excited to the ground state of the complex (energy E ⁰⁰ 19.27±0.07 kK, lifetime τ 160±30 μs), and the one-electron electroreduction wave with E ₁?(2.0±0.1)V. For the [Pd(dphpz)·(N∧N)]⁺ complexes containing diazine derivatives of phenanthroline (dppz, dicnq), the degradation of the photoexcitation energy from two electronically excited states can occur as isolated process with successive transfer of electrons to the π orbitals localized on the remote moieties: [Pd(dphpz)] and diazine fragments of the N∧N ligands.     

DFT quantum-chemical calculations of nitrogen oxide chemisorption and reactivity on the Cu(100) surface

A DFT quantum-chemical study of NO adsorption and reactivity on the Cu₂₀ and Cu₁₆ metal clusters showed that only the molecular form of NO is stabilized on the copper surface. The heat of monomolecular adsorption was calculated to be ΔH _m = ?49.9 kJ/mol, while dissociative adsorption of NO is energetically unfavorable, ΔH _d = + 15.7 kJ/mol, and dissociation demands a very high activation energy, E _a = + 125.4 kJ/mol. Because of the absence of NO dissociation on the copper surface, the formation mechanism of the reduction products, N₂ and N₂O, is debatable since the surface reaction ultimately leads to N-O bond cleavage. As the reaction occurs with a very low activation energy, E _a = 7.3 kJ/mol, interpretation of the NO direct reduction mechanism is both an important and intriguing problem because the binding energy in the NO molecule is high (630 kJ/mol) and the experimental studies revealed only physically adsorbed forms on the copper surface. It was found that the formation mechanism of the N₂ and N₂O reduction products involves formation (on the copper surface) of the (O_adN-NO_ad) dimer intermediate that is chemisorbed via the oxygen atoms and characterized by a stable N-N bond (r _N-N ～1.3 Å). The N-N binding between the adsorbed NO molecules occurs through electron-accepting interaction between the oxygen atoms in NO and the metal atoms on the “defective” copper surface. The electronic structure of the (O_adN-NO_ad) surface dimer is characterized by excess electron density (ON-NO)^δ? and high reactivity in N-O_ad bond dissociation. The calculated activation energy of the destruction of the chemisorbed intermediate (O_adN-NO_ad) is very low (E _a = 5–10 kJ/mol), which shows that it is kinetically unstable against the instantaneous release of the N₂ and N₂O reduction products into the gas phase and cannot be identified by modern experimental methods of metal surface studies. At the same time, on the MgO surface and in the individual (Ph₃P)₂Pt(O₂N₂) complex, a stable (O_adN-NO_ad) dimer was revealed experimentally.     

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.     

Note on the use of the vacuum stability test in the study of initiation reactivity of attractive cyclic nitramines in Formex P1 matrix

The zero-order reaction rates (specific rate constants) of isothermal decomposition at 120 °C of plastic bonded explosives (PBXs) were measured by means of the Czech vacuum stability test, STABIL. The PBXs are based on 1,3,5-trinitro-1,3,5-triazinane (RDX), 1,3,5,7-tetranitro-1,3,5,7-tetrazocane (HMX), cis-1,3,4,6-tetranitro-octahydroimidazo-[4,5-d]imidazole (BCHMX), and ε 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (ε-HNIW, ε-CL-20) with 13 wt% of the Formex P1 type matrix, i.e., a matrix of the explosive with pentaerythritol tetranitrate (PETN) bound by 13 wt% of a mixture of 25 wt% of styrene–butadiene rubber and 75 wt% of an oily material. Dependencies were found between the specific rate constants mentioned and the detonation velocities of PBXs, and consequently between these constants and the impact and electric spark sensitivities of pure explosive fillers, i.e., RDX, HMX, HNIW, BCHMX, and PETN. It is stated that the higher impact or electric spark sensitivity of their pure explosive fillers corresponds to the higher thermal reactivity of the given PBXs.     

A new approach to angular momentum constraints in bimolecular reactions

The well-known analysis of angular momentum constraints for diatomic systems in terms of a bond potential containing a J-dependent centrifugal part is here extended to larger systems. The new procedure is based on the identification of the minimal kinetic energy associated with a given spatial configuration and J value. Alternatively, for a given total energy E and spatial configuration the maximum J value can be identified. Constant J_max curves have been obtained by a Monte Carlo sampling scheme for a model of the KNaCl system at several different total internal energies. Particular attention is given to the estimation of J-dependent reaction thresholds.     

Molecular dynamics simulation of local motion of polystyrene chain end—comparison with the fluorescence depolarization study

Molecular dynamics (MD) simulation of the local motion of a polystyrene (PS) chain with anthryl group at the chain end surrounded by benzene molecules was performed and the results were compared with those obtained experimentally by the fluorescence depolarization method. The molecular weight dependence of the relaxation time of the probe obtained by the MD simulation was qualitatively in agreement with the results obtained by the fluorescence depolarization method. We also estimated the molecular weight dependence of the relaxation time for the end-to-end vector. Below the degree of polymerization (DP)≤3, the mean relaxation time T_m for the end-to-end vector was similar to that for the vector corresponding to the transition moment of the probe. With the increase of DP, the T_m for the probe tended to reach an asymptotic value unlike that for the end-to-end vector, which monotonically increased with DP. This indicates that the entire motion of a polymer coil contributes to the local motion to a lesser extent as the molecular weight increases. The MD simulations using artificial restraints showed that the rotational relaxation of the probe at the chain end for a dynamically stiff PS chain is realized by the cooperative rotation of the main chain bonds. The internal modes which takes place below 5 monomer units mainly led to the rotational relaxation of the probe at the PS chain end. Finally, the change of T_m with the position along the PS main chain was examined.     

The regression of isothermal thermogravimetric data to evaluate degradation Ea values of polymers: A comparison with literature methods and an evaluation of lifetime predictions reliability. Part II

The kinetics of the isothermal degradation in static air atmosphere of four well known polymers, polyethylene (PE), polystyrene (PS), polycarbonate (PC) and poly(methyl methacrylate) (PMMA) was studied by both a long-term (more than three years) experiment at relatively low temperature (423 K) and a set of short-term experiments at higher temperatures. The activation energy (E_a) values of degradation were determined by both the MacCallum and Wilkinson literature methods, and were compared with those obtained through a new very simple method we set up, based on the direct regression of TG mass loss data. About two years ago we published the results concerning PE and PS because their mass losses during long-term experiments were sufficiently high. The long-term degradation experiments were continued until now and in this second part we report the results concerning PC and PMMA. The degradation E_a values calculated from short-term experimental data through the three different methods were in good agreement with each other for both PC and PMMA, thus confirming the general applicability of our simple method for the determination of E_a. The experimental data at lower temperature of PC were not in agreement with those at higher temperatures, thus confirming the low reliability of the kinetic parameters (and then of lifetime predictions) at low temperature determined by experiments at higher temperatures. Partially disagreeing results were obtained for PMMA, which were discussed and interpreted.     

Energy transfer in reactive and nonreactive collisions of Na with Na2

Stimulated by the experimental finding of vibrationally and rotationally cold dimers in supersonic nozzle molecular beams of sodium, we have studied energy transfer in collisions of Na with Na₂ over a wide range of initial relative translation energies E and impact parameters b by a classical mechanical trajectory method. The vibrational and rotational energies were initialized using Boltzmann distributions characterized by temperatures T_vib = 150 K, T_rot = 50 K. We find that for large values of E the energy transfer in reactive collisions increases with b while it decreases with b for the nonreactive collisions. For low values of E, energy transfer is a decreasing function of b for both reactive and nonreactive encounters. Both the reactive and nonreactive mechanisms are very efficient in effecting transfer, between 40–70% of the initial relative translational energy is converted into internal energy of the diatom, leading to the conclusion that the reverse collisions would result in the rapid relaxation observed in experiment.     

Molecular dynamics simulation and density functional theory insight into the cocrystal explosive of hexaazaisowurtzitane/nitroguanidine

Theoretical methods involving molecular dynamics (MD) simulation and density functional theory were performed to investigate the different molecular ratios, mechanical Properties, structure, trigger bond, and intermolecular interaction of hexaazaisowurtzitane (CL‐20)/nitroguanidine (NQ) cocrystal explosive. Results of MD simulation show that CL‐20 and NQ packed in ratios of 1:1 present the larger binding energy and better mechanical properties than any other molecular ratios, which indicates 1:1 cocrystal can form the stable crystal structure. Shorter length and larger dissociation energy of trigger bond in composite structure than in isolated CL‐20 component suggests that the cocrystal may exhibit less sensitive than CL‐20. Analyses of atoms in molecules, reduced density gradient, and natural bond orbital confirm that intermolecular interactions are mainly derived from a series of weak hydrogen bond and strong vdW forces, involving of NH···O, CH···O, CH···N, O···N, and O···O. Additionally, composite structures of 2 and 3 bringing us more attractive performance will act as a key role in constructing of CL‐20/NQ cocrystal explosive. © 2015 Wiley Periodicals, Inc.