取代氯苯化合物的电子结构和键离解能的理论研究

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

Theoretical Study of the N-NO2 Bond Dissociation Energies for Energetic Materials with Density Functional Theory

The N-NO2 bond dissociation energies （BDEs） for 7 energetic materials were computed by means of accurate density functional theory （B3LYP, B3PW91 and B3P86） with 6-31G＊＊ and 6-311G＊＊ basis sets. By comparing the computed energies and experimental results, we find that the B3P86/6-311G＊＊ method can give good results of BDE, which has the mean absolute deviation of 1.30kcal/mol. In addition, substituent effects were also taken into account. It is noted that the Hammett constants of substituent groups are related to the BDEs of the N-NO2 bond and the bond dissociation energies of the energetic materials studied decrease when increasing the number of NO2 group.     

Density Functional Calculations of C-NO2 Bond Dissociation Energies for Nitroalkanes Molecules

Bond dissociation energies for the removal of nitrogen dioxide group in some nitroalkane energetic materials have been calculated by using the three hybrid density functional theory （DFT） methods B3LYP, B3PW91 and B3P86 with 6-31g^＊＊ and 6-311g^＊＊ basis sets. The computed BDEs have been compared with the available experimental results. It is found that the B3P86 method with 6-31g^＊＊ and 6-311g^＊＊ basis sets can obtain satisfactory bond dissociation energies （BDEs）, which are in extraordinary agreement with the experimental data. Considering the smaller mean absolute deviation and maximum difference, the reliable B3P86/6-311g^＊＊ method was recommended to compute the BDEs for the removal of nitrogen dioxide group in the nitroalkane energetic materials. Using the method, the BDEs of 8 other nitroalkane energetic materials have been calculated and the maximum difference from experimental value is 1.76 kcal.mo1^-1 （for the BDE of tC4Hg-NOz）, which further proves the reliability of B3P86/6-311g^＊＊ method. In addition, it is noted that the BDEs of C-NO2 bond change slightly for main chain nitroalkane compounds with the maximum difference of only 3.43 kcal mo1^-1.     

Density Functional Calculations of C–NO_2 Bond Dissociation Energies for Nitroalkanes Molecules

Bond dissociation energies for the removal of nitrogen dioxide group in some nit- roalkane energetic materials have been calculated by using the three hybrid density functional theory (DFT) methods B3LYP, B3PW91 and B3P86 with 6-31g** and 6-311g** basis sets. The computed BDEs have been compared with the available experimental results. It is found that the B3P86 method with 6-31g** and 6-311g** basis sets can obtain satisfactory bond dissociation energies (BDEs), which are in extraordinary agreement with the experimental data. Considering the smaller mean absolute deviation and maximum difference, the reliable B3P86/6-311g** method was recommended to compute the BDEs for the removal of nitrogen dioxide group in the nitroalkane energetic materials. Using the method, the BDEs of 8 other nitroalkane energetic materials have been calculated and the maximum difference from experimental value is 1.76 kcal·mol-1 (for the BDE of tC4H9–NO2), which further proves the reliability of B3P86/6-311g** method. In addition, it is noted that the BDEs of C–NO2 bond change slightly for main chain nitroalkane compounds with the maximum difference of only 3.43 kcal mol-1.     

PCM study of the solvent and substituent effects on bond dissociation energies of the C?NO bond

Quantum chemical calculations are used to estimate the equilibrium C? NO bond dissociation energies (BDEs) for eight X? NO molecule (X = CCl₃, C₆F₅, CH₃, CH₃CH₂, iC₃H₇, tC₄H₉, CH₂CHCH₂, and C₆H₅CH₂). These compounds are studied by employing the hybrid density functional theory (B3LYP, B3PW91, B3P86) methods together with 6‐31G** and 6‐311G** basis sets and the complete basis set (CBS‐QB3) method. The obtained results are compared with the available experimental results. It is demonstrated that B3P86/6‐31G** and CBS‐QB3 methods are accurate for computing the reliable BDEs for the X? NO molecule. Considering the inevitably computational cost of CBS‐QB3 method and the reliability of the B3P86 calculations, B3P86 method with 6‐31G** basis set may be more suitable to calculate the BDEs of the C? NO bond. The solvent effects on the BDEs of the C? NO bond are analyzed and it is shown that the C? NO BDEs in a vacuum computed by using B3PW91/6‐311G** method are the closest to the computed values in acetontrile and the average solvent effect is 1.48 kcal/mol. Subsequently, the substituent effects of the BDEs of the C? NO bond are further analyzed and it is found that electron denoting group stabilizes the radical and as a result BDE decreases; whereas electron withdrawing group stabilizes the group state of the molecule and thus increases the BDE from the parent molecule. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

PCM Study of Bond Dissociation Energies of the S–NO Bond: A DFT Study

Quantum chemical calculations are used to estimate the equilibrium S–NO bond dissociation energies (BDEs) for seven S-nitro-N-acetyl-d,l-penicillamine dipeptides (SNAP-two peptides) in acetonitrile solution. These compounds were studied by employing the hybrid density functional theory (B3LYP, B3P86 and B3PW91) methods together with the 6-31G** basis set. The obtained results are compared with available experimental results. It is demonstrated that B3PW91 method is the best method to compute the bond dissociation energies of SNAP-two peptides. The substituent and solvent effects of the S–NO BDEs are further analyzed. The results show that S–NO BDE increases with the increment of isoelectric points of substituted groups. In addition, the S–NO BDE decreases due to the inclusion of solvent effects. Furthermore, SNAP-two peptides and the other NO-donors are compared.     

Investigation of correlation between impact sensitivities and bond dissociation energies in some triazole energetic compounds

In this article, density functional theory has been utilized to study on the correlation between impact sensitivities h _50% and the bond dissociation energies (BDEs) of nine triazole energetic explosives. By employing B3LYP and B3P86 method with the 6-311G** basis set, all the molecules have been fully optimized. The BDEs for removal of the NO₂ group in these compounds have also been calculated at the same level. Computed results show that BDEs calculated by B3LYP method are all less than those by B3P86 method. The relationship between the impact sensitivities and the weakest C–NO₂ bond dissociation energy (BDE) values have been investigated. The results indicate a good linear correlation between the impact sensitivity h _50% and the ratio (BDE/E) of the weakest BDE to the total energy E.     

Bond Dissociation Energies and Electronic Structures in a Series of Peroxy Radicals: A Theoretical Study

Quantum chemical calculations are performed to estimate the bond dissociation energies (BDEs) for 18 peroxy radicals. Since DFT methods are researched to have low basis sets sensitivity, these radicals are studied by utilizing the hybrid density functional theory (DFT) (B3LYP, B3P86, B3PW91 and PBE1PBE) in conjunction with the 6‐311G** basis set and the complete basis set (CBS‐Q) method. On the basis of comparisons of the computational results and the experimental values, we evaluate the effectiveness of above methods. It is demonstrated that CBS‐Q method is the best method for computing the reliable BDEs of C—OO bond, with the average absolute errors of 2.1 kcal/mol. So CBS‐Q method is suitable to predict accurate BDEs of C‐OO bond for peroxy compounds. The computational energy gaps between the HOMO and LUMO of studied compounds are almost identical from the point of view of stability and substantial HOMO‐LUMO gaps for all molecules suggest their electronic stability. In addition, substituent effect on the C—OO BDE of peroxy radicals is analyzed. It is noted that the effects of substitution on the C—OO BDE of peroxy radicals are significant. Our results will shed lights on future theoretical and experimental work.     

PCM study of the solvent and substituent effects on bond dissociation energies of the O?NO Bond—A DFT study

A PCM continuum model, at the B3LYP, B3P86, and B3PW91 three‐parameter hybrid DFT methods with 6‐311G** basis set, is used to study the bond dissociation energies (BDEs) of benzyl nitrites. Compared the computed results with the experimental values, it is noted that B3PW91 functional is the best method to compute the BDEs of benzyl nitrites. The solvent and substituent effects on the BDEs of the O? NO bond are analyzed, and it is shown that the BDE of the O? NO bond decreases with the increment of the Hammett constants of substituent groups on benzene for benzyl nitrites except C₆H₅CH₂O? NO. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

An evaluation of various computational methods for the treatment of organoselenium compounds

A reliable computational method for the prediction of organoselenium geometries and bond dissociation energies (BDEs) has been determined on the basis of the performance of density functional theory (DFT: B3LYP and B3PW91) and ab initio molecular orbital procedures (Hartree-Fock (HF)) in conjunction with various Pople basis sets including (but not limited to) the 6-31G(d), 6-31G(d,p), 6-311G(d), 6-311G(d,p), 6-311G(2df,p), and 6-311G(3df,3pd) sets. Predicted geometries and BDEs are compared with available experimental data and quadratic configuration interaction including single and double substitutions (QCISD) results. The B3PW91/6-311G(2df,p) level of theory is recommended for the prediction of the geometries and energetics of organoselenium compounds.     

N—O键解离焓的密度泛函理论研究

在高精度计算方法G3和G3B3的基础上,比较了密度泛函理论(DFT)十几种方法对N—O键解离焓(BDE)相对于实验值的计算精度,发现用B3P86方法计算15种化合物N—O键的BDE,均方根误差最小,仅为6.36kJ·mol-1,计算值与实验值的线性相关系数为0.991.在此基础上,用该方法分别计算了非芳香化合物及芳香化合物的N—O键BDE.通过自然键轨道分析,发现部分N—O键的BDE与N—O键的键长、原子电荷密度及键级之间存在定量关系.此外,在B3P86方法的基础上预测了几种典型的杂环芳香化合物N—O键BDE值.     

Theoretical studies on four-membered ring compounds with NF2, ONO2, N3, and NO2 groups

Density functional theory (DFT) method has been employed to study the geometric and electronic structures of a series of four-membered ring compounds at the B3LYP/6-311G** and the B3P86/6-311G** levels. In the isodesmic reactions designed for the computation of heats of formation (HOFs), 3,3-dimethyl-oxetane, azetidine, and cyclobutane were chosen as reference compounds. The HOFs for N(3) substituted derivations are larger than those of oxetane compounds with --ONO2 and/or --NF2 substituent groups. The HOFs for oxetane with --ONO2 and/or --NF2 substituent groups are negative, while the HOFs for N3 substituted derivations are positive. For azetidine compounds, the substituent groups within the azetidine ring affect the HOFs, which increase as the difluoroamino group being replaced by the nitro group. The magnitudes of intramolecular group interactions were predicted through the disproportionation energies. The strain energy (SE) for the title compounds has been calculated using homodesmotic reactions. For azetidine compounds, the NF2 group connecting N atom in the ring decrease the SE of title compounds. Thermal stability were evaluated via bond dissociation energies (BDE) at the UB3LYP/6-311G** level. For the oxetane compounds, the O--NO2 bond is easier to break than that of the ring C--C bond. For the azetidine and cyclobutane compounds, the homolyses of C--NX2 and/or N--NX2 (X = O, F) bonds are primary step for bond dissociation. Detonation properties of the title compounds were evaluated by using the Kamlet-Jacobs equation based on the calculated densities and HOFs. It is found that 1,1-dinitro-3,3-bis(difluoroamino)-cyclobutane, with predicted density of ca. 1.9 g/cm(3), detonation velocity (D) over 9 km/s, and detonation pressure (P) of 41 GPa that are lager than those of TNAZ, is expected to be a novel candidate of high energy density materials (HEDMs). The detonation data of nitro-BDFAA and TNCB are also close to the requirements for HEDMs.     

环境污染物中碳氯键离解能的理论研究

利用六种密度泛函理论方法（B3LYP, B3P86, MPW1K, TPSS1KCIS, X3LYP, BMK）对碳氯键离解能进行理论计算,结果发现几种新发展的密度泛函(DFT)方法用于碳氯键离解能的计算比传统的B3LYP有较大的改善,其中对能量估算相对准确的B3P86方法对碳氯键离解能的计算精度最高,对17个分子中碳氯键离解能计算的平均绝对偏差为6.58 kJ/mol。最后运用B3P86方法对一系列环境危害较大,但可通过光化学降解和生物降解的氯代有机物的碳氯键离解能值进行预测,并讨论了影响碳氯键离解能的结构性质关系。     

S=O homolytic bond dissociation enthalpies in sulfoxides

The S=O bond dissociation enthalpies (BDE) were calculated using high-level ab initio methods including G3, G3B3, CBS-Q, CBS-4M, CCSD(T), and MP2. Based on the comparison of these theoretical values and experimental ones, the performances of a number of density functional theory (DFT) methods were then assessed. It was found that the B3P86 method gave the lowest root of mean square error. We therefore used this method to calculate the S=O BDEs of a number of substituted sulfoxides. The electronic effect of the substituents and the remote substituents effect of aryl-substituted sulfoxides on S=O BDE were investigated. In addition, a molecular orbital analysis of typical molecules was conducted in order to investigate the electronic effect on S=O BDEs. We also predicted several S=O BDE values of heteroaromatic substituted sulfoxides using the B3P86 method.     

Theoretical Studies on Energetic Property and Stability of Pyrazine and Pyridine Derivatives

Density functional calculations at the B3LYP level with 6‐311G** and aug‐cc‐pVDZ basis sets were performed to predict the heats of formation (HOFs) for two pyrazine derivatives and eight pyridine derivatives. In the isodesmic reactions designed for the computation of heats of formation (HOFs), pyrazine and pyridine were chosen as reference compounds. The N‐oxidations for the ring nitrogen of pyrazine and pyridine derivatives decrease the HOF values when N‐oxide oxygen is neighboring with amino groups, but increase when it neighbors with nitro groups. Thermal stability was evaluated via bond dissociation energies (BDE) at the UB3LYP/6‐311G** level. As a whole, the homolysis of C–NO₂ bonds is the main step for bond dissociation of the title compounds. The BDE values of title compounds are influenced by intramolecular hydrogen bonds. The hydrogen bond effects associated with the length of the H···O bonds were analyzed by the electron density at the critical points and natural bond orbital.     

Bond dissociation energies for removal of the hydroxyl group in some alcohols from quantum chemical calculations

In this theoretical work, 22 alcohols and their geometric structure properties have been investigated employing quantum chemical methods to calculate the C? OH equilibrium bond distances and bond dissociation energies (BDEs). Since DFT methods have been researched to have low basis sets sensitivity for small and medium molecules in our previous work (Zhao et al., J Mol Struct, 2006, 766, 87), 22 title compounds have been studied by employing the hybrid density functional theory (B3LYP, B3PW91, B3P86, PBE1PBE) in conjunction with the 6‐311G** basis set and the complete basis set (CBS–Q) method. Comparison with the available experimental data shows that CBS–Q and B3P86 methods calculated results agree very well with the experimental values, with the average absolute errors of 1.3 kcal/mol and 3.5 kcal/mol, respectively. So considering the expensive computational time, CBS–Q method can be chosen as a satisfactory method of predicting the accurate BDEs for removal of the OH group in small and medium size alcohols. And B3P86 method may give accurate BDEs for larger alcohols we haven't studied. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Computational DFT studies on a series of toluene derivatives as potential high energy density compounds

Based on the full optimized molecular geometric structures at B3LYP/6-31G**, B3LYP/6-31+G**, B3P86/6-31G**, and B3P86/6-31+G** levels, the densities (ρ), detonation velocities (D), and pressures (P) for a series of toluene derivatives, as well as their thermal stabilities, were investigated to look for high energy density compounds (HEDCs). The heats of formation (HOFs) are also calculated via designed isodesmic reactions. The calculations on the bond dissociation energies (BDEs) indicate that the BDEs of the initial scission step are between 48 and 59 kcal/mol, and pentanitrotoluene is the most reactive compound, while 2,4,6-trinitrotoluene is the least reactive compound for toluene derivatives studied. A good linear relationship between BDE/E and impact sensitivity is also obtained. The condensed phase HOFs are calculated for the title compounds. These results would provide basic information for the further studies of HEDCs. The detonation data of pentanitrotoluene show that it meets the requirement for HEDCs.     

The C-H and alpha(C-X) bond dissociation enthalpies of toluene, C6H5-CH2X (X = F, Cl), and their substituted derivatives: a DFT study

The homolytic C-H bond dissociation enthalpies (BDEs) of toluene and its para- and meta-substituted derivatives have been estimated by using the (RO)B3LYP/6-311++G(2df,2p)//(U)B3LYP/6-311G(d,p) procedure. The performance of two other hybrid functionals of DFT, namely, B3PWP91 and O3LYP, has also been evaluated using the same basis sets and molecules. Our computed results are compared with the available experimental values and are found to be in good agreement. The (RO)B3LYP and (RO)O3LYP procedures are found to produce reliable BDEs for the C-H bonds in toluene and the C-X (X = F, Cl) bond in alpha-substituted toluene (C6H5-CH2X) and their substituted derivatives. The substituent effect on the BDE values has been analyzed in terms of the ground-state effect and the radical effect. The effect of polarization of the C-H bond on the substituent effect is also analyzed. The BDE(C-H) and BDE(C-X) values for alpha-substituted (X = F and Cl) toluenes with a set of para substituents are presented for the first time.     

密度泛函法计算C1-C14正构烃生成焓及C-C键裂解能

采用密度泛函理论(DFT)方法研究了费托石脑油裂解反应中涉及到C₁-C₁₄正构烃和自由基中间体的生成焓及其C-C键解离能(BDE)。 结果表明,在所有评价的密度泛函理论方法（B97-1、BB1K、B1B95、MPWB1K和MPW1B95）中,MPW1B95/6-311G(d,p)方法计算最精确。 以此方法为基准,进一步对高碳烃及其裂解产物的标准生成焓和C-C键解离能进行了预测。 与可得到的实验数据相比,MPW1B95/6-311G(d,p)方法预测的烃和自由基的平均生成焓分别为0.8和2.7 kJ/mol,C-C键解离能的平均绝对误差只有3.1 kJ/mol,表明此方法不仅可准确计算正构烃标准生成焓和C-C键解离能,而且还能正确预测C-C键解离能变化趋势。     

1, 4-Pyrone Effects on O-H Bond Dissociation Energies of Catechols in Flavonoids: A Density Functional Theory Study

In recent years, there has been growing interest in selecting efficient antioxidants with low toxicity to reduce the damage of free radicals. Among these antioxidants, flavonoids have been paid much attention, owing to their excellent antioxidative and pharmacological activities1. Up to now, many efforts have been given to summarize the structure-activity relationships (SAR) for flavonoids. It has been widely accepted that two structural factors are critical for flavonoids to enhance the…