硝酸酯热分解的机理分析

使用加速量热仪研究硝酸正丙酯（NPN）、硝酸异丙酯（IPN）和硝酸异辛酯（EHN）的热分解过程。分析了在绝热条件下热分解反应动力学, 根据实验数据计算出表观活化能、指前因子和反应热等参数。根据NPN、IPN和EHN热分解的起始温度和反应热数据,给出了三种硝酸酯在75℃时的反应风险指数,分析得到三种物质的热失稳风险度。分别在B3LYP/6-31+G(2df,2p)和B3P86/6-31+G(2df,2p)的理论水平下,计算得到NPN、IPN和EHN的O-NO2键离解能(BDE)。由NPN、IPN和EHN的O-NO2键离解能在很大程度上符合由加速量热仪测试得到的活化能,推知三种硝酸酯的热分解反应只是单分子O-NO2键的均裂反应。     

两种太根(Tri-EGDN和Tetra-EGDN)的几何结构和键裂解能的理论研究

利用密度泛函理论(DFT)中的B3LYP、B3PW91、B3P86方法研究了二缩三乙二醇二硝酸酯(Tri-EGDN)和三缩四乙二醇二硝酸酯(Tetra-EGDN)的平衡分子构型、前线轨道能(EHOMO、ELUMO)、能差(ΔE=ELUMO-EHOMO)和O-NO2键的裂解能(BDE),并分析了两种硝酸酯分子和对应自由基的热稳定性.Tri-EGDN和Tetra-EGDN中的O-NO2键的裂解能用B3P86方法计算最接近真实值.由O-NO2键的裂解能很好地符合文献中对应物质的表观活化能,推知两种物质的热分解反应只是单分子O-NO2键的均裂反应.     

氢键和分子-离子相互作用对奥克托今(HMX)引发键影响的理论研究

为了寻找能够降低炸药感度的方法,借助B3LYP和MP2(full)理论,在B3LYP/6-311+ +G(2df,2p)、B3LYP/aug-cc-pVTZ和MP2(full)/6-311++G(2df,2p)三种水平下,分别研究了与HMX形成的12种复合物（六种氢键复合物和六种分子-离子复合物）,并对引发键N-NO2键长和强度的变化及硝基电荷的变化等进行了详细的考量,最后借助分子中原子(AIM)理论揭示了结构和能量变化的本质.结果表明,在形成复合物后,引发键键长变短,离解能增大,硝基电荷增多,引发键增强,感度降低,同时引发键离解能的变化(ΔBDEs)与这两种相互作用能呈良好的线性关系(R2=0.9984).     

硝酸丙酯键离解能和热解机理的密度泛函理论研究

 采用B3LYP方法、在6-31G*基组水平上,优化了硝酸丙酯单体的平衡几何构型,计算了分子中各键的离解能。结果表明该分子中最弱的键为O—NO₂键,次弱键为C—O键。同时研究了热解机理,探索性的预测了断裂O—NO₂键及C—O键后的产物,并采用从头算法、半经验方法和密度泛函理论分别计算了硝酸丙酯的生成热,由半经验方法中的PM3得到的数值和实验结果符合较好。     

硝酸溴与三重态氧原子反应机理的密度泛函研究

用密度泛函理论方法 (B3LYP) ,在 6 311+G(d ,p)水平上对硝酸溴与三重态氧原子的反应进行了研究 ,计算了反应中各驻点物种的平衡构型、振动频率、总能量和零点能 (ZPE) .对计算得到各可能反应途径的过渡态经内禀反应坐标分析加以了证实 ,对反应途径中的键长和能量的变化作了IRC解析 .在B3LYP优化的基础上利用了耦合簇理论方法 (CCSD(T) )在 6 311+G(d ,p)水平上对各驻点物种的单点能进行了修正 .研究表明 ,存在三种可能的反应途径 ,其产物分别为 :cis BrONO和 3 O2 、trans BrONO和 3 O2 以及BrOO和NO2 .其中第三个通道由于活化能垒较低 ,是主要反应 .     

氢键和分子-离子相互作用对奥克托今(HMX)引发键影响的理论研究

为了寻找能够降低炸药感度的方法,借助B3LYP和MP2(full)理论,在B3LYP/6-311++G(2df,2p)、B3LYP/aug-cc-p VTZ和MP2(full)/6-311++G(2df,2p)三种水平下,分别研究了与HMX形成的12种复合物(六种氢键复合物和六种分子-离子复合物),并对引发键N-NO2键长和强度的变化及硝基电荷的变化等进行了详细的考量,最后借助分子中原子(AIM)理论揭示了结构和能量变化的本质.结果表明,在形成复合物后,引发键键长变短,离解能增大,硝基电荷增多,引发键增强,感度降低,同时引发键离解能的变化(ΔBDEs)与这两种相互作用能呈良好的线性关系(R2=0.9984).     

多硝基苯酸酯炸药分子键离解能与撞击感度关系研究

用密度泛函方法,在B3LYP-6-31G*、B3P86/6-31G*和B3LYP-6-311G*三种理论水平对四个含硝基烷基的苯酸酯炸药分子进行了几何结构全优化、能量和频率计算。并对这些炸药分子苯环上的C-NO2和烷基上的C-NO2键离解能分别进行了三种理论水平的计算。结果表明,这类分子中的最弱键是烷基上的C-NO2键。进一步分析实验撞击感度与分子中最弱键离解能量的关系,结果表明,分子最弱键离解能与分子总能量的比值BDE/E和实验撞击感度h50%之间存在一个几乎线性的关联关系。BDE/E是表征炸药撞击感度的一个实用的合理指针。     

木质素模型化合物热解机理的理论研究

采用密度泛函理论方法B3LYP/6-31G (d,p),对β-O-4型木质素二聚体模化物的热解反应机理进行了量子化学理论研究。提出了三种可能的热解反应途径：Cβ-O键均裂的后续反应、Cα-Cβ键均裂的后续反应以及协同反应。计算了各热解反应途径的标准动力学参数,分析了各种主要热解产物的形成演化机理。计算结果表明,β-O-4型模化物中Cβ-O的键离解能最低,其次是Cα-Cβ的。Cβ-O键均裂的后续反应和协同反应路径（3）是主要的反应路径,而Cα-Cβ键均裂的后续和协同反应路径（1）和（2）是热解过程中主要的竞争反应路径。     

BCl和BCl2(2A1)的结构与势能函数

采用多种方法,配有多种基组对BCl和BCl2分子的基态结构进行优化计算,优选出B3P86/6-311++G(3df,3pd)方法对BCl分子进行计算得到基态为 、键长 =1.7159nm,谐振频率为 ;优选出QCISD/6-31G(d,p)方法对BCl2分子进行计算得到基态为 ,平衡核间距RBCl=0.17284nm、键角 =125.3466o、离解能 =8.0592eV,并计算出了谐振频率和力常数．在此基础上,运用多体展式理论方法,推导出BCl2分子的解析势能函数,其等值势能面图准确呈现出BCl2分子的结构特征及能量变化曲线．由此讨论了Cl+BCl和B+ClCl分子反应的势能面特征．可用于研究该分子的微观反应动力学特性．     

AlC,SiC基态分子结构与分析势能函数的量子力学计算

用密度泛函理论的B3LYP方法和二次组态相互作用(QCISD(T))方法,选择6-31G(d,p)、6-311 G(2df,2pd)、6-311 G(3df,3pd)、cc-PVTZ、AUG-cc-PVTZ基组,优化计算了AlC和SiC分子基态的能量,平衡结构,谐振频率.根据原子分子反应静力学原理,导出了AlC和SiC分子基态的合理离解极限和离解能.通过优化计算结果和实验数据的对比,选择QCISD(T)/6-311 G(3df,3pd)方法对AlC和SiC分子基态的势能面进行了单点能扫描.采用最小二乘法拟合得到了AlC和SiC分子基态的Murell-Sor-bie势能函数.同时计算了光谱参数(Be,eα,ωe,ωeχe)和力常数(f2,f3,f4),并与实验结果进行比较.结果表明,计算结果与实验数据吻合的较好.     

Theoretical Study of the C-Cl Bond Dissociation Enthalpy and Electronic Structure of Substituted Chlorobenzene Compounds

利用密度泛函(DFT)三种交换/相关函数(B3LYP, B3PW91,B3P86)结合6-31G**和6-311G**基组,计算了13个取代氯苯化合物的键离解能. 结果表明B3P86/6-311G**方法是计算取代氯苯化合物键离解能的可信方法,研究发现C-Cl键的键离解能与所使用的基组和计算方法密切相关,取代基对C-Cl键的键离解能的影响不明显. 研究了目标化合物的前线轨道能级差,并对取代氯苯化合物的热稳定性做了评估.     

Ignition and combustion of isopropyl nitrate

The ignition delay times for mixtures of isopropyl nitrate (IPN) with air and argon are measured in a rapid-injection reactor at a pressure of 1 atm and in a shock tube at 2–3 atm. It is shown that the ignition delay time τ of mixtures in which heat is largely released due to oxidation by the oxygen contained in the IPN molecule is determined by the unimolecular decomposition of IPN over the entire temperature range covered (500–730 K). For mixtures in which heat is mainly produced by oxidation reactions involving air oxygen, the ignition delay time at high temperatures is controlled by secondary reactions of oxidation of the hydrocarbon moiety of the IPN molecule, leading to an increase in τ by more than an order of magnitude. Liquid IPN burns in a nitrogen atmosphere only at pressures above 40 atm, at a linear rate of ~4 mm/s. The measured flame temperatures are in close agreement with the respective values calculated using a thermodynamic code.     

Quantum chemical calculations of bond dissociation energies for COOH scission and electronic structure in some acids

Quantum chemical calculations are performed to investigate the equilibrium C-COOH bond distances and the bond dissociation energies(BDEs) for 15 acids.These compounds are studied by utilizing the hybrid density functional theory(DFT)(B3LYP,B3PW91,B3P86,PBE1PBE) and the complete basis set(CBS-Q) method in conjunction with the 6311G** basis as DFT methods have been found to have low basis sets sensitivity for small and medium molecules in our previous work.Comparisons between the computational results and the experimental values reveal that CBS-Q method,which can produce reasonable BDEs for some systems in our previous work,seems unable to predict accurate BDEs here.However,the B3P86 calculated results accord very well with the experimental values,within an average absolute error of 2.3 kcal/mol.Thus,B3P86 method is suitable for computing the reliable BDEs of C-COOH bond for carboxylic acid compounds.In addition,the energy gaps between the highest occupied molecular orbital(HOMO) and the lowest unoccupied molecular orbital(LUMO) of studied compounds are estimated,based on which the relative thermal stabilities of the studied acids are also discussed.     

A combined quantum mechanical and experimental approach towards chiral diketopiperazine hydroperoxides

In order to improve hydroperoxide formation from heterocyclic compounds relating to the formation rate and to allow a suitable choice of starting materials for autoxidation, theoretical studies on a set of different amino acid‐derived diketopiperazines and pyrazinoquinazolines were carried out. To estimate their reactivity towards hydroperoxide formation, bond dissociation enthalpies (BDEs) of tertiary α‐C? H bonds as well as reaction enthalpies to the corresponding hydroperoxides were calculated at the B3LYP/TZVP and RMP2/aug‐cc‐pVTZ level of theory. The Evans–Polanyi relation was then used to correlate substrate reactivity with calculated BDEs. Thermal and zero point vibrational energy (ZPE) corrections were determined in the classical harmonic oscillator‐rigid rotor‐particle in a box model. While for the investigated set of diketopiperazines BDEs of 318.8–327.0 kJ mol^?1 were found, BDEs for pyrazinoquinazolines spread between 248.4 and 368.4 kJ mol^?1 at the B3LYP/TZVP level of theory. A selected subset of heterocycles was converted to the corresponding hydroperoxides and the diketopiperazines were obtained in up to 39% yield after 5–7 days, whereas the pyrazinoquinazoline hydroperoxides were isolated in up to 67% yield after 24 h. Thus, replacing an amido moiety in an N‐aryl‐imino moiety when using pyrazinoquinazolines instead of diketopiperazines leads indeed to an improved captodative stabilization of the radical intermediate. Furthermore the theoretical calculations allowed a distinctive forecast of the preferred regioisomeric hydroperoxide. Copyright © 2010 John Wiley & Sons, Ltd.     

组氨酸电离能与红外光谱的密度泛涵理论计算研究

采用混合密度泛涵理论中的B3LYP方法,结合4种基组6-31G(d),6-31G(df, p),6-31+G(d)和6-311+G(2d, 2p),系统计算了光合反应中心叶绿素的配位体-组氨酸在空气、四碌化碳、四氢呋喃、水和蛋白质模拟环境中的几何结构、电离能、红外光谱及同位素标记谱。计算及分析结果表明：组氨酸分子的几何参数在不同计算基组和介质中略有不同,且C₂—N₃,N₃—C₄的键长在空气中最大;同一介质中,增大计算基组和采用扩散函数,均使计算的单点势能和振动频率降低,电离能增加,对应光谱强度增高;而同一计算方法下,介质的电介常数越高,分子的单点势能越低,电离能越小,对应的振动频率减小强度增加;另外,电离能和主要特征峰位及其15N和13C标记谱的计算结果与文献中的实验结果相吻合。所有计算均显示,高基组和施加扩散函数的计算结果与实验更接近。该研究为深入探索叶绿素与组氨酸配位后在光合反应中心的功能与振动光谱特性提供理论参考依据。     

Heterocyclic analogs of phenol as novel potential antioxidants

Five density functional theory (DFT) methods including B3LYP, B3PW91, MPW1K, MPWB, TPSS1KCIS have been evaluated by comparing with the experimental O? H bond dissociation enthalpies (BDEs) of substituted phenols. B3PW91 is found to be the best method, for which the calculation error was 3.62 kJ/mol. Subsequently, the BDEs (O? H) of hydroxyl groups on five‐ and six‐membered heteroatomic aromatic rings have been calculated using the (RO)B3PW91/6‐311++G(2df,2p)//(U)B3LYP/6‐311g(d,p) procedure. In addition, the ionization energy (IE) and proton affinity [PA(O^?)] of these compounds have also been examined. On the basis of our theoretical study, a series of imidazolols, thiazolols, and oxazolols were studied to assess their antioxidant activities. It was found that 5‐oxazolol could be a promising novel antioxidant precursor. Copyright © 2009 John Wiley & Sons, Ltd.     

Theoretical studies on a series of nitroaliphatic energetic compounds

Density functional theory calculations at the B3LYP/6-311G＊＊ level are performed to study the geometric and elec- tronic structures of a series of nitroaliphatic compounds. The heats of formation （HOF） are predicted through the designed isodesmic reactions. Thermal stabilities are evaluated via the homolytic bond dissociation energies （BDEs）. Further, the correlation is developed between impact sensitivity h50% and the ratio （BDE/E） of the weakest BDE to the total energy E containing zero point energy correction. In addition, the relative stability of the title compounds is evaluated based on the analysis of calculated Mulliken population and the energy gaps between the frontier orbitals. The calculated BDEs, HOFs, and energy gaps consistently indicate that compound 1,1,1,6,6,6-hexanitro-3-hexyne is the most unstable and the compound 3,3,4,4,-tetranitro-hexane is the most stable. These results provide basic information for the molecular design of novel high energetic density materials.     

The evaluation of bond dissociation energies for NO2 scission in nitro compounds using density functional and complete basis set methods

By using the density functional theory （B3LYP） and four highly accurate complete basis set （CBS-Q, CBS-QB3, CBS-Lq, and CBS-4M）ab initio methods, the X（C, N, O）-NO2 bond dissociation energies （BDEs） for CH3NO2, C2H3NO2, C2H5NO2, HONO2, CH3ONO2, C2H5ONO2, NH2NO2 （CH3）2NNO2 are computed. By comparing the computed BDEs and experimental results, it is found that the B3LYP method is unable to predict satisfactorily the results of bond dissociation energy （BDE）; however, all four CBS models are generally able to give reliable predication of the X（C, N, O）-NO2 BDEs for these nitro compounds. Moreover, the CBS-4M calculation is the least computationally demanding among the four CBS methods considered, Therefore, we recommend CBS-4M method as a reliable method of computing the BDEs for this nitro compound system.     

Interfacial reaction of silver ultra-thin film deposited on interpenetrating polymer network substrate by liquor-phase reduction

The interfacial reaction, metal transformations, and nonmetal bond types of silver ultra-thin film deposited on polyurethane (PU) based interpenetrating polymer networks (IPN) substrate by the liquor-phase reduction at room temperatures were studied by atomic force microscope (AFM), X-ray photoelectron spectroscopy (XPS) and attenuated total reflection Fourier transform infrared spectroscopy (ATR-FTIR). The IPN substrate was prepared by dip-pulling precursors onto a silicon wafer or a glass plate, followed by solidification at room temperature. The interpenetrate structures of IPN with two crosslinked networks restricted the aggregation of silver during the reduction and deposition. The devised -OH terminal group in PU simplified the determination of reactive site in IPN and reinforced the adhesion between IPN and silver through interfacial reaction. The XPS results, which matched well with the ATR-FTIR results, verified the chemical reactive site of PU in IPN with silver in the oxide state.     

双原子分子内氛重组能的Ab Initio-Morse函数研究

基于重组织现象的自交换模型,提出一种内氛重组能精确的Morse函数标度方案,并对气相电子转移过程中双原子分子体系进行了计算,给出了反应物的结构重组织指标。在HFSCF水平上利用四种基组(6-31G^*、6-31GE、DZ和DZP(d))对标度中所含参数进行了计算,并与实验结果进行比较.结果表明,ab initio计算与实验标度值吻合较好,误差在±1.skcal/mol以内。该方法不仅改进了经典模型而且避免了获得有效的实验数据的困难,给出了一种精确地计算气相电子转移过程中双原子分子内氛重组能的ab initio算法。