适用于TATB,RDX,HMX含能材料的全原子力场的建立与验证

报道一个适用于三种常见的含能材料分子三硝基三氨基苯(TATB),环三亚甲基三硝胺(RDX),环四亚甲基四硝胺(HMX)的全原子力场.力场采用广泛使用的力场函数形式,其中键参数通过拟合量子化学密度泛函计算的数据获得,电荷参数和范德华参数通过拟合相应的分子晶体的物理性质(密度和升华焓)优化得到.通过计算分子和分子晶体的性质显示该力场可以用来准确地预测分子结构、分子振动频率和分子晶体的晶胞参数、密度和升华焓.进一步的验证显示该力场可用来较为准确地预测分子晶体的状态方程和机械模量.     

Theoretical study of local structure for Ni2+ ions at tetragonal sites in K2ZnF4:Ni2+ system

A theoretical method for studying the local lattice structure of Ni2+ ions in (NiF6)(4-) coordination complex is presented. Using the ligand-field model, the formulas relating the microscopic spin Hamiltonian parameters with the crystal structure parameters are derived. Based on the theoretical formulas, the 45 x 45 complete energy matrices for d8 (d2) configuration ions in a tetragonal ligand-field are constructed. By diagonalizing the complete energy matrices, the local distortion structure parameters (R perpendicular and R || ) of Ni2+ ions in K2ZnF4:Ni2+ system have been investigated. The theoretical results are accorded well with the experimental values. Moreover, to understand the detailed physical and chemical properties of the fluoroperovskite crystals, the theoretical values of the g factor of K2ZnF4:Ni2+ system at 78 and 290 K are reported first.     

β-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晶体热膨胀各向异性的本质即特定的分子堆积模式.     

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.     

High-pressure studies of 1,3,5,7-cyclooctatetraene: experiment and theory

High-pressure studies of 1,3,5,7-cyclooctatetraene have been performed by using Raman spectroscopy up to 16 GPa and compared with complementary density functional calculations. Angular-dispersive X-ray diffraction studies were also performed in the solid state at elevated pressure. The lattice constants of solid 1,3,5,7-cyclooctatetraene obtained from the X-ray diffraction pattern taken at 3.8 GPa and room temperature are in good agreement with theoretical results. At least two phase transitions were observed during pressure increase followed by the loss of long-range crystallographic order, which was also associated with a strong pressure-induced luminescence that allowed estimation of band gap alterations with pressure.     

A rigorous comparison of theoretical autocorrelation functions with dielectric relaxation parameters for the alpha process of polymers

Recently, tables of parameters used to represent experimental dielectric relaxation data as well as autocorrelation functions have become available. The experimental and autocorrelation function data were represented with the Havriliak–Negami function using rigorous statistical techniques. These tables include not only parameters and their temperature dependencies, but also the confidence intervals for all of the parameters. The important parameters for this work are the two shape parameters, α and β, which represent the width and skewness of the relaxation process, respectively. A comparison is made between the experimental data and several autocorrelation functions by calculating the minimum distance (in units of standard deviations) between the experimental values of αβ parameters for a specific polymer and the αβ parameters corresponding to the autocorrelation functions reported in these tables. Quantities derived from these minimum distances (a distribution function and an error function) are reported for each of the autocorrelation functions. These results are discussed in terms of the basic assumptions of the mode coupling theory. © 1997 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys 35: 1887–1897, 1997     

Separation of 1,3,5,7‐tetranitro‐1,3,5,7‐tetraazacyclooctane and 1,3,5‐trinitro‐1,3,5‐ triazacyclohexane by molecularly imprinted solid‐phase extraction

Synthesis of 1,3,5,7‐tetranitro‐1,3,5,7‐tetraazacyclooctane and 1,3,5‐trinitro‐1,3,5‐triazacyclohexane by the Bachmann process leads to a mixture of both. The separation of 1,3,5,7‐tetranitro‐1,3,5,7‐tetraazacyclooctane and 1,3,5‐trinitro‐1,3,5‐triazacyclohexane from their mixture is difficult because the sizes and physical properties of these homologous compounds are similar. For this purpose, seven molecularly imprinted polymers have been synthesized for each explosive, and a selective solid‐phase extraction procedure has been developed. A molecularly imprinted polymer, synthesized with 1,3,5,7‐tetranitro‐1,3,5,7‐tetraazacyclooctane as the template, methacrylic acid as the monomer and trimethylolpropane trimethacrylate as the cross‐linking agent in a molar ratio of 1:8:8 showed the best separation capability. A packed cartridge containing this polymer can be reused for 23 solid‐phase extraction cycles without repacking, and the total separation capability toward 1,3,5,7‐tetranitro‐1,3,5,7‐tetraazacyclooctane reached 6.81 mg per gram of polymer. 1,3,5‐Trinitro‐1,3,5‐triazacyclohexane was not detected in the separated 1,3,5,7‐tetranitro‐1,3,5,7‐tetraazacyclooctane by high‐performance liquid chromatography and vice versa. This newly developed method had the advantages of high recovery (100%) and purity, environmental friendliness, and room temperature operability. This study showed that some molecularly imprinted polymers that cannot absorb target analytes well in the solvent in which the polymers were polymerized might have high‐binding capacity for the analytes and show imprinting effects in other solvents.     

Ab initio kinetics of gas phase decomposition reactions

The thermal and kinetic aspects of gas phase decomposition reactions can be extremely complex due to a large number of parameters, a variety of possible intermediates, and an overlap in thermal decomposition traces. The experimental determination of the activation energies is particularly difficult when several possible reaction pathways coexist in the thermal decomposition. Ab initio calculations intended to provide an interpretation of the experiment are often of little help if they produce only the activation barriers and ignore the kinetics of the decomposition process. To overcome this ambiguity, a theoretical study of a complete picture of gas phase thermo-decomposition, including reaction energies, activation barriers, and reaction rates, is illustrated with the example of the β-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) molecule by means of quantum-chemical calculations. We study three types of major decomposition reactions characteristic of nitramines: the HONO elimination, the NONO rearrangement, and the N-NO(2) homolysis. The reaction rates were determined using the conventional transition state theory for the HONO and NONO decompositions and the variational transition state theory for the N-NO(2) homolysis. Our calculations show that the HMX decomposition process is more complex than it was previously believed to be and is defined by a combination of reactions at any given temperature. At all temperatures, the direct N-NO(2) homolysis prevails with the activation barrier at 38.1 kcal/mol. The nitro-nitrite isomerization and the HONO elimination, with the activation barriers at 46.3 and 39.4 kcal/mol, respectively, are slow reactions at all temperatures. The obtained conclusions provide a consistent interpretation for the reported experimental data.     

Evolution of Structural, Electronic and Optical Properties of Monoclinic ZrO2 under High Pressure： A First Principles Study

The structural, electronic and optical properties of the monoclinic ZrO2 were studied by ab initio calculations based on the density functional theory and pseudopotential method. The calculated lattice parameters and band gap are in agreement with the experimental and other theoretical values. The evolution of lattice parameters and electronic properties were illustrated under high pressure. Meanwhile, the optical properties, such as adsorption coefficients, imaginary part of dielectric function, and energy loss function, were investigated under both ambient and high pressures.     

Anisotropy of the elastic properties of crystalline cellulose Iβ from first principles density functional theory with Van der Waals interactions

In spite of the significant potential of cellulose nanocrystals as functional nanoparticles for numerous applications, a fundamental understanding of the mechanical properties of defect-free, crystalline cellulose is still lacking. In this paper, the elasticity matrix for cellulose I_β with hydrogen bonding network A was calculated using ab initio density functional theory with a semi-empirical correction for van der Waals interactions. The computed Young’s modulus is found to be 206 GPa along [001] (c-axis), 98 GPa along [010] (b-axis), and 19 GPa along [100] (a-axis). Full compliance matrices are reported for 1.0, 1.5 and 2.0 % applied strains Color contour surfaces that show variations of the Young’s modulus and average Poisson’s ratio with crystallographic direction revealed the extreme anisotropies of these important mechanical properties. The sensitivity of the elastic parameters to misalignments in the crystal were examined with 2D polar plots within selected planes containing specific bonding characteristics; these are used to explain the substantial variability in the reported experimental Young’s moduli values. Results for the lattice directions [001], [010] and [100] are within the range of reported experimental and other numerical values.     

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%以上,可直接作为结构鉴定的标准样品.     

Photochemical preparation of 1,3,5,7-tetracyanoadamantane and its conversion to 1,3,5,7-tetrakis(aminomethyl)adamantane

New adamantane derivatives 1 and 2 that bear functionalized one-carbon extensions at all four bridgehead positions have been prepared. Radical nucleophilic substitution (S(RN)1) reaction of 1,3,5,7-tetrabromoadamantane with cyanide produces 1,3,5,7-tetracyanoadamantane (1), which was reduced with borane reagents to 1,3,5,7-tetrakis(aminomethyl)adamantane (2). Improvements in the preparation of 1,3,5,7-tetrahaloadamantanes (halogen = Br, Cl, I) are also reported. [structure--see text]     

碳氢化合物燃烧中间体Lennard-Jones系数的计算

在已有的基团贡献法公式的基础上,提出了一种新的基团贡献法公式,并通过拟合250种化合物(包括185种稳定化合物临界性质的实验值和65种自由基临界性质的计算值)的临界性质得到了40种基团的贡献值,并用于预测未知化合物的临界性质.选取了训练集以外的、有临界性质实验值的30种化合物作为独立测试集,用于验证所建模型对临界性质的预测能力,T_C和P_C平均绝对偏差分别为8.52%和16.83%.结果表明,预测结果和实验值相吻合,该模型可以用于大分子化合物及自由基的临界性质预测.根据临界性质与Lennard-Jones(L-J)系数的经验关系式,预测了碳氢化合物燃烧中间体的L-J系数,得到独立测试集46种碳氢化合物的L-J系数,与文献值接近,T_C和P_C的平均绝对偏差分别为9.88%和9.96%.比较了训练集中烷烃自由基·C_6H_(13)、烯烃自由基·C_5H_9和炔烃自由基·C_5H_7同分异构体的L-J系数,同时,将己烷自由基·C_6H_(13)与相似的邻近烷烃C_6H_(14)的L-J系数进行比较,发现同分异构体之间或相似化合物之间L-J系数有较大偏差.此外,对缺少L-J系数的114种常见碳氢化合物自由基进行了预测.这对于碳氢化合物的燃烧模拟及基元反应中压强相关的速率常数计算有重要意义.     

Molecular structure and spectroscopic properties of polyaromatic heterocycles by first principle calculations: spectroscopic shifts with the adsorption of thiophene on phyllosilicate surface

The molecular and electronic structure and spectroscopic properties of thiophene, benzothiophene, and dibenzothiophene molecules have been studied theoretically, using different levels of quantum-mechanical calculations based on Hartree–Fock and density functional theory (DFT) approximations. In all calculation levels used, the molecular structures, dipole moments, thermodynamic properties, and vibration normal modes agree well with the available experimental data and predict the non-available values. The calculated frequencies agree with experimental values within scale factors in the range of 0.94–1.0 for thiophene. Our calculations have refined the mode assignments previously reported. The adsorption of thiophene on a pyrophyllite surface along the 001 plane was also investigated by means of a quantum-mechanical method based on the DFT approximation in a periodical crystal lattice model. The adsorbed thiophene adopts a parallel orientation with respect to the phyllosilicate surface. The shifts observed experimentally in the IR bands of thiophene upon adsorption on pyrophyllite are well reproduced in the thiophene–phyllosilicate complex, finding explanations of this phenomenon and new effects with our calculations.     

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

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

Structure and bonding in beta-HMX-characterization of a trans-annular N...N interaction

Chemical bonding in the beta-phase of the 1,3,5,7-tetranitro-1,3,5,7-tetraazacyclooctane (HMX) crystal based on the experimental electron density obtained from X-ray diffraction data at 20 K, and solid state theoretical calculations, has been analyzed in terms of the quantum theory of atoms in molecules. Features of the intra- and intermolecular bond critical points and the oxygen atom lone-pair locations are discussed. An unusual N...N bonding interaction across the 8-membered ring has been discovered and characterized. Hydrogen bonding, O...O and O...C intermolecular interactions are reported. Atomic charges and features of the electrostatic potential are discussed.     

Controlled strain on a double-templated textured polymer film: a new approach to patterned surfaces with bravais lattices and chains

A double template-assisted fabrication method for making surface patterns with tunable lattice geometries on a polymer surface is reported. This technique is based on a locally nonuniform strain produced in a double-templated polymer film that has a strong modulation in thickness. It can produce all 2D primitive Bravais lattices as well as chains on the surface of a polymer. The lattice parameters are controllable with nanoprecision by varying the direction and amount of the applied strain.     

Force constants and many-body electrostatic interactions in two-dimensional colloidal crystal

A model of a two-dimensional colloidal crystal with a hexagonal lattice, the electrostatic interactions in which are described by the nonlinear Poisson-Boltzmann equation, is considered. The calculation procedure for force constants of this crystal is treated in detail. Properties of system symmetry, which make it possible to significantly decrease the volume of calculations and to classify force constants, are analyzed. Numerical data for force constants of a crystal as functions of lattice parameters at different particle sizes are reported. A method that allows us to disclose the presence of many-body interactions in a system by the behavior of force constants at some interval of the values of lattice parameters is proposed. The application of this method to the system under consideration demonstrated that electrostatic interparticle interactions in the system cannot be reduced to simply a pair interaction of any kind; the introduction of many-body potentials is required for the adequate representation of the elastic properties of a crystal.     

Dendritic PtCo alloy nanoparticles as high performance oxygen reduction catalysts

We studied the physical properties and the concentration profile of benzene+water+caprolactam mixtures near the fluid-fluid interface using self-consistent field (SCF) theory. This yields the interfacial tension which plays an important role in describing the stability of transient liquid droplets of one phase in the other. The studies were performed at a fixed temperature of 313K. Flory-Huggins binary interaction parameters and the compound lattice segment numbers are input parameters for the applied SCF theory. These parameters were derived from activity coefficient relations, which are used to describe experimental liquid-liquid and vapor-liquid phase equilibrium measurements. Using first principles, the benzene-water interface was studied and the resulting interfacial tension was found to be in agreement with experimental values. This study illustrates that caprolactam accumulates at the benzene-water interface, acting as a weak surfactant. The interfacial tension is also demonstrated to be affected by the caprolactam concentration and the SCF results are in fair agreement with the experimental observations.     

Application of global kinetic models to HMX beta-delta transition and cookoff processes

The reduction of the number of reactions in kinetic models for both the HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine) beta-delta phase transition and thermal cookoff provides an attractive alternative to traditional multi-stage kinetic models due to reduced calibration effort requirements. In this study, we use the LLNL code ALE3D to provide calibrated kinetic parameters for a two-reaction bidirectional beta-delta HMX phase transition model based on Sandia instrumented thermal ignition (SITI) and scaled thermal explosion (STEX) temperature history curves, and a Prout-Tompkins cookoff model based on one-dimensional time to explosion (ODTX) data. Results show that the two-reaction bidirectional beta-delta transition model presented here agrees as well with STEX and SITI temperature history curves as a reversible four-reaction Arrhenius model yet requires an order of magnitude less computational effort. In addition, a single-reaction Prout-Tompkins model calibrated to ODTX data provides better agreement with ODTX data than a traditional multistep Arrhenius model and can contain up to 90% fewer chemistry-limited time steps for low-temperature ODTX simulations. Manual calibration methods for the Prout-Tompkins kinetics provide much better agreement with ODTX experimental data than parameters derived from differential scanning calorimetry (DSC) measurements at atmospheric pressure. The predicted surface temperature at explosion for STEX cookoff simulations is a weak function of the cookoff model used, and a reduction of up to 15% of chemistry-limited time steps can be achieved by neglecting the beta-delta transition for this type of simulation. Finally, the inclusion of the beta-delta transition model in the overall kinetics model can affect the predicted time to explosion by 1% for the traditional multistep Arrhenius approach, and up to 11% using a Prout-Tompkins cookoff model.