X-H (X = C, N, O, Si, P, S) 键离解能计算的密度泛函方法研究

使用了不同密度泛函方法计算X-H (X = C, N, O, Si, P, S) 键离解能，并分析不同密度泛函方法的计算精度。研究发现大多数密度泛函方法包括B3LYP, B3P86, B3PW91, G96LYP, PBE1PBE，和BH&HLYP都明显低估键离解能13-25 kJ/mol。该现象与是否使用无限基组无关，因为即使使用无限基组键离解能仍然被低估。因此密度泛函方法不适合用于键离解能的估算。其中B3P86方法的偏差最小。进一步分析表明，使用限制性开壳层计算并无任何优势，在大多数情况下非限制性开壳层计算实际上比限制性开壳层计算要好。最后，我们发现了密度泛函方法对键离解能的低估是系统的，因此建议利用校准后的UDFT/6-311++G(d, p)方法计算化学键离解能。     

碳氟键均裂解离能的密度泛函理论研究

运用6种密度泛函方法(B3LYP, B3P86, B3PW91, PBE1PBE, MPW1B95, MPW1K)对15个含氟有机化合物的碳氟键均裂解离能进行理论计算, 得到的理论值与实验值比较, 发现B3P86方法用于碳氟键均裂解离能的计算相对可靠. 使用验证后的理论方法对含氟杂环有机化合物和卤氟烃中的碳氟键均裂解离能进行了预测和分析, 并进一步讨论了α-取代基效应以及Hammett型取代基效应对碳氟键均裂解离能的影响.     

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

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

ClO/ClO~-耦合体系结构性质及反应通道研究

利用从头算法和密度泛函理论对ClO/ClO~-体系进行了研究。首先应用密度泛 函理论的六种方法(B3LYP,BLYP,B3P86,BP86,BHLYP,LSDA)和从头算理论的CCSD方法 在6-3+G~*,6-311+G~*及6-311+G(3df)基组水平上对单体进行了优化,选出最适合 该体系的方法和基组B3LYP/6-311+G(3df);然后在B3LYP/6-311+G(3df)水平上计算 了沿各种反应通道的离解能,并且进行了校正。找出了存在的中间体及中间体异构 化的过渡态,并进行了IRC路径解析。对各稳定体进行了频率分析和成键分析。结 果表明,单体ClO和单体ClO~-结合为稳定的中间体后,其离解方式主要是向着 ClOO+Cl~-和ClOO~-+Cl两种方式进行,两种离解方式的离解能分别为-33.39和82. 88kJ/mol,并且前者是一个离解性电子转移过程,经历一个电子转移过渡态。     

精确计算化学键解离能的ONIOM-G3B3方法及其在抗氧化剂研究中的应用

通过对179个有机化合物C—H, N—H, O—H键离解能的理论计算, 系统评估了高精度组合从头算方法(ONIOM-G3B3)和密度泛函理论方法(B3LYP)在预测键离解能上的可靠性. 研究发现ONIOM-G3B3方法可以准确预测各类有机化合物的键解离能, 精度达到5.9 kJ/mol. 运用ONIOM-G3B3方法成功预测了两类重要的天然抗氧化剂维生素E族和茶多酚族的键解离能, 并进一步讨论了抗氧化活性、自由基清除机理及其构效关系.     

精确计算化学键解离能的ONIOM-G3B3方法及其在抗氧化剂研究中的应用

通过对179个有机化合物C—H,N—H,O—H键离解能的理论计算,系统评估了高精度组合从头算方法(ONIOM-G3B3)和密度泛函理论方法(B3LYP)在预测键离解能上的可靠性.研究发现ONIOM-G3B3方法可以准确预测各类有机化合物的键解离能,精度达到5.9 kJ/mol.运用ONIOM-G3B3方法成功预测了两类重要的天然抗氧化剂维生素E族和茶多酚族的键解离能,并进一步讨论了抗氧化活性、自由基清除机理及其构效关系.     

木质素模化物键离解能的理论研究

采用密度泛函理论B3P86方法,在6-31G(d,p)基组水平上,对木质素结构中的6种连接方式(β-O-4、α-O-4、4-O-5、β-1、α-1、5-5)的63个木质素模化物的醚键(C-O)和C-C键的键离解能E_B进行了理论计算研究。分析了不同取代基对键离解能的影响以及键长与键离解能的相关性。计算结果表明,C-O键的键离解能通常比C-C键的小,在各种醚键中C_α-O键的平均键离解能最小,为182.7 kJ/mol;其次是β-O-4连接中的C_β-O键,苯环和烷烃基上的取代基对醚键的键离解能有较强的弱化作用,C-O键的键长和键离解能的相关性较差。与C-O键相比,C-C键的键离解能受苯环上取代基的影响很小,而烷烃基上的取代基对C-C键的键离解能有较大的影响,C-C键的键离解能和键长之间存在较强的线性关系,C-C键的键长越长,其键离解能越小。     

以双亚甲胺席夫碱为配体的双核吡唑铜配合物的磁性理论研究

采用对称性破损态方法结合密度泛函理论,以双亚甲胺席夫碱为配体的双核吡唑铜配合物[Cu_2L(PZ)](L=双甲亚胺三阴离子,由1-苯基-3-甲基-4-甲酰吡唑和1,3-二氨-2-丙醇衍生而成;PZ=吡唑阴离子)作为研究对象,通过与实验数据相比较,讨论了不同密度泛函方法与基组对金属铜配合物磁交换耦合常数的影响.结果表明,4种混合密度泛函方法(B3PW91,B3LYP,B3P86和PBE)及3种单一密度泛函方法(BPW91,BLYP和BP86)的计算结果都能与实验值符号一致,且B3P86方法所得到的结果与实验值最为接近,而单一密度泛函的计算结果误差较大,与实验值吻合程度不好.同时采用B3P86方法计算所得交换耦合常数Jab对基组的依赖性较大.研究表明,2个Cu(Ⅱ)离子之间弱的反铁磁相互作用主要源于单占据分子轨道SOMOs小的能量劈裂.     

铝原子团簇Al5、Al5-和Al5+稳定结构的密度泛函理论研究

采用密度泛函理论的四种方法:杂化密度泛函B3LYP与B3PW91、Perdew-Wang91交换与相关泛函WP91PW91、局域自旋密度近似SVWN,研究了A15、Al5-和Al5+团簇的多种可能结构,找到了它们稳定的结构与自旋态,与已有的理论结果作了比较,并计算了Al5-的绝热与垂直电子离解能、Al5的绝热与垂直电离势,同有关的实验数据比较,符合较好.同时对四种密度泛函方法的计算结果作了一些比较与讨论.     

Au(II)化合物[Au(CH_2)_2PH_2]_2X_2(X=F,Cl,Br,I)的量子化学理论研究

采用从头计算HF,MP2方法和密度泛函理论,对Au(II)系列化合物[Au(CH2)2PH2]2X2(X=F,Cl,Br,I)的几何结构、电子结构和振动频率进行了研究.研究表明Au的5d和6s电子参与Au—Au以及Au—X之间的成键.Au—Au,Au—X键强烈的电子相关作用使HF方法不适于该体系的研究,BP86和B3LYP两种泛函给出较大的Au—Au和Au—X键长,而MP2方法和局域的密度泛函方法则给出了合理的结构参数.局域密度泛函方法计算得到的Au—Au键和Au—X键振动频率也与实验数据符合较好.还运用含时密度泛函理论计算了[Au(CH2)2PH2]2X2的电子激发能,对分子在紫外-可见光谱范围内的电子跃迁进行了分析,考察了卤素配体对激发能的影响,并结合分子轨道能级的变化对此给予了解释.     

DFT study of the C-Cl bond dissociation enthalpies and electronic structure of substituted chlorobenzene compounds

Quantum chemical calculations were used to estimate the bond dissociation energies (BDEs) for 13 substituted chlorobenzene compounds. These compounds were studied by employing the hybrid density functional theory methods (B3LYP, B3PW1, B3P86) with 6-31G** and 6-311G** basis sets. It was demonstrated that B3P86/6-311G** method is the best method for computing the reliable BDEs for substituted chlorobenzene compounds which contain the C-Cl bond. It was found that the C-Cl BDE depends strongly on a computational method and basis set used. Substitution effect on the C-Cl BDE of substituted chlorobenzene compounds is further discussed. It is shown that the effects of substitution on the C-Cl BDE of substituted chlorobenzene compounds are very insignificant. Frontier orbital energy gap of studied compounds was also investigated. From the data on frontier orbital energies gap, we estimated the relative thermal stability of substituted chlorobenzene compounds.     

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.     

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     

Theoretical studies of bond dissociation energies for removal of the nitrile group in some nitrile compounds

The CCN bond distances and bond dissociation energies (BDEs) are estimated by utilizing quantum chemical calculations for 16 nitrile compounds. Since DFT methods have been researched to have low basis sets sensitivity for small and medium molecules in our earlier work [Jun Zhao, Xinlu Cheng, Xiangdong. Yang, J. Mol. Struct. (Theochem) 766 (2006) 87] 16 nitrile compounds are studied by employing the hybrid density functional theory (B3LYP, B3PW91, B3P86) and the complete basis set (CBS-Q) method in conjunction with the 6-311G^** basis set. The obtained results are compared with the available experimental data. It is demonstrated that CBS-Q method, which can produce reasonable BDEs for some systems, seems unable to predict accurate BDEs here. While, the B3P86 calculated results agree very well with the experimental values. So B3P86 method is suitable for computing the reliable BDEs of CCN bond for nitrile compounds.     

Assessment of experimental bond dissociation energies using composite ab initio methods and evaluation of the performances of density functional methods in the calculation of bond dissociation energies

Composite ab initio CBS-Q and G3 methods were used to calculate the bond dissociation energies (BDEs) of over 200 compounds listed in CRC Handbook of Chemistry and Physics (2002 ed.). It was found that these two methods agree with each other excellently in the calculation of BDEs, and they can predict BDEs within 10 kJ/mol of the experimental values. Using these two methods, it was found that among the examined compounds 161 experimental BDEs are valid because the standard deviation between the experimental and theoretical values for them is only 8.6 kJ/mol. Nevertheless, 40 BDEs listed in the Handbook may be highly inaccurate as the experimental and theoretical values for them differ by over 20 kJ/mol. Furthermore, 11 BDEs listed in the Handbook may be seriously flawed as the experimental and theoretical values for them differ by over 40 kJ/mol. Using the 161 cautiously validated experimental BDEs, we then assessed the performances of the standard density functional (DFT) methods including B3LYP, B3P86, B3PW91, and BH&HLYP in the calculation of BDEs. It was found that the BH&HLYP method performed poorly for the BDE calculations. B3LYP, B3P86, and B3PW91, however, performed reasonably well for the calculation of BDEs with standard deviations of about 12.1-18.0 kJ/mol. Nonetheless, all the DFT methods underestimated the BDEs by 4-17 kJ/mol in average. Sometimes, the underestimation by the DFT methods could be as high as 40-60 kJ/mol. Therefore, the DFT methods were more reliable for relative BDE calculations than for absolute BDE calculations. Finally, it was observed that the basis set effects on the BDEs calculated by the DFT methods were usually small except for the heteroatom-hydrogen BDEs.     

Theoretical calculation of bond dissociation energies and heats of formation for nitromethane and polynitromethanes with density functional theory

The C? NO₂ bond dissociation energies (BDEs) and the heats of formation (HOFs) of nitromethane and polynitromethanes (dinitromethane, trinitromethane, and tetranitromethane) system in gas phase at 298.15 K were calculated theoretically. Density functional theory (DFT) B3LYP, B3P86, B3PW91, and PBE0 methods in combination with different basis sets were employed. It was found that the C? NO₂ bond BDEs can be improved from B3LYP to B3PW91 to B3P86 or PBE0 functional. Levels of theory employing B3P86 and PBE0 functionals were found to be sufficiently reliable without the presence of diffusion functions. As the number of NO₂ groups on the same C atom increases, the PBE0 functional performs better than the B3P86 functional. Regarding the calculated HOFs, all four functionals can yield satisfactory results with deviations of <2 kcal mol^?1 from experimental ones for CH₂(NO₂)₂ and CH(NO₂)₃, when the diffusion functions are not augmented. For the C(NO₂)₄ molecule, the large basis sets augmented with polarization functions and diffusion functions are required to yield a good result. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

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.     

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.