苯和苯胺类硝基衍生物热解机理和撞击感度的理论研究

在DFT-B3LYP/6-31G*水平获得苯和苯胺类硝基衍生物的全优化分子几何和电子结构.通过非限制性(U)B3LYP/6-31G*计算求得标题物各化学键离解能(BDE).用UHF-PM3 MO方法求得引发键C-NO2键均裂反应的活化能(Ea).以静态指标(键集居数、前线轨道能级差和硝基上净电荷)和动态理论指标(BDE和Ea)阐明了热解引发机理,关联了实验撞击感度.运用SPSS程序关联静态和动态理论指标,表明它们均可以用作预示标题物的热解引发机理和撞击感度.     

反式-1,4,5,8-四硝基-1,4,5,8-四氮杂萘烷异构体热解机理与稳定性的理论研究

运用密度泛函理论和半经验分子轨道方法,对一系列高能杂环硝胺—反式-1,4,5,8-四硝基-1,4,5,8-四氮杂萘烷异构体的热解机理和稳定性进行了系统地计算研究。在B3LYP/6-31G**和PM3水平上,分别计算了标题物的化学键离解能(BDE)和热解反应活化能(Ea),并根据BDE和Ea数值考察了硝胺取代基对化合物稳定性和热解机理的影响;同时,还详细考察了BDE与Ea、化学键重叠布居数、前线轨道能级以及能隙之间的相关性。结果表明,由BDE、Ea和静态电子结构参数推断的标题物热稳定性和热解机理的结论基本是一致的,N-NO2键均裂是标题物的热解引发步骤,间位取代异构体较对位取代异构体稳定,而邻位取代的异构体稳定性最差。     

新型高能量密度化合物3,6-双(3,5-二硝基-1,2,4-三唑-1)-1,2,4,5-四嗪-1,4-二氧化物的性能预估及合成路线设计

采用C++自编译程序及组合原理,设计并筛选出一种未见报道的新型富氮类高能量密度化合物-3,6-双(3,5.二硝基.1,2,4-三唑.1)-1,2,4,5-四嗪-1,4-二氧化物,用B3LYP法,在6-31G**基组水平上得到该化合物全优化构型;在振动分析的基础上求得体系的振动频率、IR谱;通过键级分析得到热解引发键的键离解能(BDE);采用Monte-Carlo 方法预估了密度;设计等键等电子反应计算了生成焓;运用Kamlet-Jacobs公式预测爆速、爆压和爆热;运用Keshavarz 等推导的预估撞击感度H50的公式预测了撞击感度性能;并利用逆合成分析法设计其合成路线.结果表明:该化合物存在8个强吸收峰,校正后的热解引发键的BDE为264KJ·mol-1,稳定性较优;密度1.955 g·cm-3、生成焓901.72 kJ·mol-1、爆速9191.48 m·s-1、爆压39.32 GPa、爆热6705.15 j·g-1;撞击感度H50为55.85cm,低于黑索金(RDX)和奥克托今(HMx);以上性能均达到了高能量密度化合物的标准,且该化合物设计合成路线步骤较少、原料易得,有望得到广泛应用.     

哒嗪裂解反应机理的密度泛函理论研究

利用密度泛函理论(DFT)方法, 以哒嗪为煤的模型化合物, 在(U)B3LYP/6-31G(d,p)水平下计算了标题化合物的键布居数和键裂解能, 并对其热解机理进行了探讨. 在(U)B3LYP/6-311++G(d,p)水平下计算了热解过程中各物种的单点能, 并对总能量进行了校正. 结果表明, 标题物N-N键的布居数和裂解能均最小, 且各键裂解能随温度变化不大, 热解最终产物为HCN和乙炔, NH₃可能是HCN经过二次转换得到的. 采用过渡态理论计算了300~1900 K温度范围内热解的速率常数, 求得lnk与1/T的线性关系.     

双环-HMX结构和性质的理论研究

在DFT-B3LYP/6-311G*水平上, 计算研究了高能化合物四硝基四氮杂双环辛烷(双环-HMX) α和β两种异构体的结构和性质. 比较分子对称性、分子内氢键和环张力等几何参数以及分子总能量和前线轨道能级等电子结构参数, 发现α比β稳定. 分子中N—N键较长, N—N键集居数较小, 预示该键为热解和起爆的引发键. 基于简谐振动分析求得IR谱频率和强度. 运用统计热力学方法求得200～1000 K温度的热力学性质. 以非限制性半经验PM3方法探讨其热解机理, 求得各反应通道的过渡态和活化能, 发现热解始于侧链N—NO₂键的均裂. 还从理论上预测了该化合物的密度、爆速和爆压, 有助于寻求高能量密度材料(HEDM).     

1,3,5,7-四硝基金刚烷结构和性能的理论研究

在DFT-B3LYP/6-31G*水平下求得1,3,5,7-四硝基金刚烷的全优化分子几何和电子结构. 经简谐振动分析求得其IR谱并作归属. 由统计热力学求得其不同温度下的热力学性质. 以非限制性半经验MO方法探讨其热解机理, 求得各反应通道的过渡态和活化能, 发现热解始于侧链C—NO₂键的均裂. 还基于理论计算密度和生成热, 以Kamlet-Jacobs方程估算其爆速和爆压.     

环三甲撑三硝胺(RDX)结构和性质的DFT研究

环三甲撑三硝胺(RDX)是综合性能好、使用极广的高能炸药。本文用密度泛函理论(DFT)B3LYP方法,取6-31G*基组,求得其分子几何构型、电子结构、IR谱和热力学性质。全优化几何构型和电子结构均具有C3V对称性。在相邻原子之间以NNO2键的Mulliken集居数最小,表明其间电子分布较少,预示其为热解和起爆的引发键。IR谱与实验结果良好相符。计算所得298-1200K的热力学性质,对进一步深入研究RDX的反应和性质有助。     

Li(NTO)(H2O)2的热分解行为及其结构与性质的关系研究

摘要在水中合成了3-硝基-1,2,4．三唑-5-酮(NTO)的锂盐Li(NTO)(H2O)2,利用DSC,TG／DTG和IR方法研究了其金属配合物的热分解机理,并用Kissinger法、Ozawa法、积分法和微分法对标题配合物进行了非等温动力学研究,得到了热分解反应的动力学参数,确定了热分解第一阶段的动力学方程及配合物的热爆炸临界温度(Tb)为289．33℃．采取RHF／6-31G,DFT-RB3LYP／6-31G方法对标题化合物进行了几何全优化,并对其成键情况、电荷分布和化合物的稳定性进行了分析．     

苯的硝基和叠氮基衍生物的理论研究

在密度泛函理论B3LYP/6-31G*水平下优化了91个苯的硝基(NO₂)和叠氮基(N₃)衍生物的分子几何构型, 预测了它们的密度和生成热, 采用Kamlet-Jacobs方法计算了爆速和爆压, 筛选得到11种爆轰性能较好的高能量密度化合物(HEDC), 计算了它们的多个可能的热解引发键的键离解能(BDE)以及按“氧化呋咱机理”分解时的活化能(E_a). 结果表明, 当分子中有NO₂与N₃相邻时, 分解按“氧化呋咱机理”进行, 分解反应的E_a均大于100 kJ/mol|分子中没有NO₂和N₃相邻时, 热解始于C-NO₂或C-N₃均裂, 裂解的BDE都大于200 kJ/mol. 只含NO₂或N₃的7个物质的稳定性好于同时含NO₂和N₃的物质, 而只含N₃的物质的稳定性又好于只含NO₂的物质, 五叠氮苯和六叠氮苯具有很出色的爆轰性能和稳定性. 无论是能量还是稳定性方面, 筛选得到的11种物质基本符合HEDC的要求.     

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

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

Theoretical investigation on density,detonation properties,and pyrolysis mechanism of nitro derivatives of benzene and aminobenzenes

Nitro derivatives of benzene and aminobenzenes are optimized at the DFT‐B3LYP/6‐31G* level. The heat of formation (ΔH_f) and crystal theoretical density (ρ) are estimated to evaluate the detonation properties using the modified Kamlet–Jacobs equations. Thermal stability and the pyrolysis mechanism of the title compounds are investigated by calculating the bond dissociation energies (BDE) at the unrestricted B3LYP/6‐31G* level. The kinetic parameter and the static electronic structural parameters can be used to predict the stability and the relative magnitude of the impact sensitivity of homologues. According to the quantitative standard of the energy and the stability as an HEDC, the title compounds having more than four nitro groups satisfy this requirement. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

A theoretical investigation on the structures,densities, detonation properties,and pyrolysis mechanism of the nitro derivatives of phenols

The nitro derivatives of phenols are optimized to obtain their molecular geometries and electronic structures at the DFT‐B3LYP/6‐31G* level. Detonation properties are evaluated using the modified Kamlet–Jacobs equations based on the calculated densities and heats of formation. It is found that there are good linear relationships between density, detonation velocity, detonation pressure, and the number of nitro and hydroxy groups. Thermal stability and pyrolysis mechanism of the title compounds are investigated by calculating the bond dissociation energies (BDEs) at the unrestricted B3LYP/6‐31G* level. The activation energies of H‐transfer reaction is smaller than the BDEs of all bonds and this illustrates that the pyrolysis of the title compounds may be started from breaking O? H bond followed by the isomerization reaction of H transfer. Moreover, the C? NO₂ bond with the smaller bond overlap population and the smaller BDE will also overlap may be before homolysis. According to the quantitative standard of energetics and stability as a high‐energy density compound, pentanitrophenol essentially satisfies this requirement. In addition, we have discussed the effect of the nitro and hydroxy groups on the static electronic structural parameters and the kinetic parameter. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Theoretical Studies on the Reaction Mechanism of 1-Chloroethane with Hydroxyl Radical

The reaction mechanism of 1-chloroethane with hydroxyl radical has been inves- tigated by using density functional theory (DFT) B3LYP/6-31G (d, p) method. All bond dissociation enthalpies were computed at the same theoretical level. It was found that hydrogen abstraction pathway is the most favorable. There are two hydrogen abstraction pathways with activation barriers of 0.630 and 4.988 kJ/mol, respectively, while chlorine abstraction pathway was not found. It was observed that activation energies have a more reasonable correlation with the reaction enthalpy changes (△Hr) than with bond dissociation enthalpies (BDE).     

Kinetics and mechanism of 2-pyridylacetic acid pyrolysis in the gas phase: A joint experimental and theoretical study

A combination of the experimental and theoretical study was carried out on the reaction mechanism associated with the pyrolysis of 2-pyridylacetic acid in the gas phase. Methylpyridine and carbon dioxide were analyzed as the products, using a static system over the pressure range of 18–55 torr and the temperature of 541.2–583.4 K. The experimental kinetic data show that the pyrolysis process is homogeneous, unimolecular and proceeds through a concerted mechanism. Theoretical studies at the B3LYP level using the 6-31G^* basis set confirmed an asynchronous concerted mechanism for the reaction. Computed kinetic and activation parameters are in good agreement with the experimental one.     

Theoretical studies on the structures, thermodynamic properties, detonation properties, and pyrolysis mechanisms of spiro nitramines

Density function theory (DFT) has been employed to study the geometric and electronic structures of a series of spiro nitramines at the B3LYP/6-31G level. The calculated results agree reasonably with available experimental data. Thermodynamic properties derived from the infrared spectra on the basis of statistical thermodynamic principles are linearly correlated with the number of nitramine groups as well as the temperature. Detonation performances were evaluated by the Kamlet-Jacobs equations based on the calculated densities and heats of formation. It is found that some compounds with the predicted densities of ca. 1.9 g/cm3, detonation velocities over 9 km/s, and detonation pressures of about 39 GPa (some even over 40 GPa) may be novel potential candidates of high energy density materials (HEDMs). Thermal stability and the pyrolysis mechanism of the title compounds were investigated by calculating the bond dissociation energies (BDE) at the B3LYP/6-31G level and the activation energies (E(a)) with the selected PM3 semiempirical molecular orbital (MO) based on the unrestricted Hartree-Fock model. The relationships between BDE, E(a), and the electronic structures of the spiro nitramines were discussed in detail. Thermal stabilities and decomposition mechanisms of the title compounds derived from the B3LYP/6-31G BDE and the UHF-PM3 E(a) are basically consistent. Considering the thermal stability, TNSHe (tetranitrotetraazaspirohexane), TNSH (tetranitrotetraazaspiroheptane), and TNSO (tetranitrotetraazaspirooctane) are recommended as the preferred candidates of HEDMs. These results may provide basic information for the molecular design of HEDMs.     

基于平均影响值的反向传播神经网络方法用于提高密度泛函理论计算Y—NO键均裂能精度

将基于平均影响值(Mean impact value,MIV)的反向传播神经网络(Back propagation neural netowrk,BPNN)(MIV-BPNN)方法用于提高密度泛函理论(Density functional theory,DFT)计算Y—NO(Y=N,S,O及C)键均裂能的精度.利用量子化学计算和MIV-BPNN联合方法计算92个含Y—NO键的有机分子体系的均裂能.结果表明,相对于单一的密度泛函理论B3LYP/6-31G(d)方法,利用全参数下的BPNN方法计算92个有机分子均裂能的均方根误差从22.25 kJ/mol减少到1.84 kJ/mol,而MIV-BPNN方法使均方根误差减少到1.36 kJ/mol,可见B3LYP/6-31G(d)和MIV-BPNN联合方法可以提高均裂能的量子化学计算精度,并可预测实验上无法获取的均裂能值.     

Water effect on the o-h dissociation enthalpy of para-substituted phenols: a DFT study

The effect of water on the O-H bond dissociation enthalpy (BDE) of para-substituted phenols has been investigated by means of DFT calculations. It is shown that the experimental BDE values are fairly well-reproduced by simple B3LYP/6-31G* calculations carried out on the phenol/phenoxyl-water complexes taking into account only hydrogen-bonding (HB) interactions of water molecules with molecular sites (HB model). On the contrary, the BDE values computed with the polarizable continuum model (PCM/B3LYP/6-31G*)8 are overestimated by about 3-4 kcal/mol. Discrepancy between theory and experiment increases using the PCM method in addition to the HB model. Calculations show that, in general, the HB interaction with water molecules decreases the BDE of phenols bearing electron-releasing groups while increasing the BDE of phenols bearing electron-withdrawing substituents. This opposite effect is explained by considering the resonance structures with charge separation both in phenols and in phenoxyl radicals. With electron donors, the phenoxyl radical is preferentially stabilized by the HB acceptor interaction with two water molecules, while with electron acceptors the phenol is preferentially stabilized by the HB donor interaction with one water molecule.     

Which DFT functional performs well in the calculation of methylcobalamin? Comparison of the B3LYP and BP86 functionals and evaluation of the impact of empirical dispersion correction

Controversy remains regarding the suitable density functionals for the calculation of vitamin B(12) systems that contain cobalt. To identify the optimum functionals, geometry optimization calculations were performed on a full-size model of methylcobalamin (MeCbl) using the B3LYP, B3LYP-D, BP86, and BP86-D methods in conjunction with the 6-31G* basis set. Single-point energy evaluations were also performed with the 6-311+G(2d,p) basis set. Consistent with previous studies, the BP86-optimized geometry showed fairly good agreement with the experimental geometry. Various factors that may influence the homolytic bond dissociation energy (BDE) of the Co-C bond of MeCbl were systematically evaluated with these methods. Our analysis demonstrated that dispersion was the largest correction term that influenced the magnitude of BDE. Previous studies have shown that B3LYP significantly underestimates BDE, whereas BP86 gives BDE values that are fairly close to the experimental values (36-37 kcal/mol). The same trend in the relative magnitudes of the BDEs was observed in the present calculations. However, BP86 underestimated the BDE for a full model of MeCbl. When the amount of Hartree-Fock exchange in the B3LYP functional was reduced to 15% and the dispersion correction was made (i.e., B3LYP*-D), the calculated BDE was in good accord with experimental values. B3P86-D also performed well. A detailed analysis was undertaken to determine which atoms in cobalamin have large dispersion interactions with a methyl fragment of MeCbl.     

Theoretical Studies on the Structures,Thermodynamic Properties,Detonation Performance,and Pyrolysis Mechanisms for Six Dinitrate Esters

Density functional theory (DFT) has been employed to study the geometric and electronic structures of six dinitrate esters including ethylene glycol dinitrate (EGDN), diethylene glycol dinitrate (Di‐EGDN), triethylene glycol dinitrate (Tri‐EGDN), tetraethylene glycol dinitrate (Tetra‐EGDN), pentaethylene glycol dinitrate (Penta‐EGDN) and hexaethylene glycol dinitrate (Hexa‐EGDN) at the B3LYP/6‐31G* level. Their IR spectra were obtained and assigned by vibrational analysis. Based on the frequencies scaled by 0.96 and the principle of statistic thermodynamics, the thermodynamic properties were evaluated, which were linearly related with the number of CH₂CH₂O groups as well as the temperature, obviously showing good group additivity. Detonation performances were evaluated by the Kamlet‐Jacobs equations based on the calculated densities and heats of formation. It was found that density, detonation velocity, detonation pressure decreased with the increase of the number of CH₂CH₂O groups. Thermal stability and the pyrolysis mechanism of the title compounds were investigated by calculating the bond dissociation energies (BDE) at the B3LYP/6‐31G* level. For the nitrate esters, the O‐NO₂ bond is a trigger bond during a thermolysis initiation process.