Accurate heats of formation of polycyclic saturated hydrocarbons predicted by using the XYG3 type of doubly hybrid functionals

Polycyclic saturated hydrocarbons (PSHs) are attractive candidates as hydrocarbon propellants. To assess their potential values, one of the key factors is to determine their energy contents, such as to calculate their heats of formation (HOF). In this work, we have calculated HOFs for a set of 36 PSHs including exo-Tricyclo[5.2.1.0^(2,6)] decane, the principal component of the high-energy density hydrocarbon fuel commonly identified as JP-10. The results from B3LYP, B3LYP-D3BJ, M06-2X, B2PLYP, B2PLYP-D3BJ, and the XYG3 type of doubly hybrid (xDH) functionals are presented. It is demonstrated here that the xDH functionals yield accurate HOFs in good agreement with those from experiments or the G4 theory. In particular, XYGJ-OS, a low scaling xDH functional, is shown to hold the promise for accurate prediction of HOFs for PSHs of larger sizes. © 2018 Wiley Periodicals, Inc.     

Accurate potential energy surfaces for hydrogen abstraction reactions: A benchmark study on the XYG3 doubly hybrid density functional

The potential energy surface (PES) for the H + CH₄ system has been constructed with the recently developed XYG3 doubly hybrid functional, while those with the standard B3LYP hybrid functional, and the Møller–Plesset perturbation theory up to the second order (MP2) are also presented for comparison. Quantum dynamics studies demonstrated that satisfactory results on the reaction probabilities and the rate coefficients can be obtained on top of the XYG3‐PES, as compared to the results based on the highly accurate, yet expensive, CCSD(T)‐PES (Li et al., J. Chem. Phys. 2015, 142, 204302). Further investigation suggested that the XYG3 functional is useful in providing accurate rate coefficients for some larger systems involving H atom abstractions. © 2017 Wiley Periodicals, Inc.     

新一代密度泛函方法XYG3

大量的应用研究测试表明,以B3LYP为代表的传统密度泛函方法在反应能垒以及非键相互作用等重要性质的预测上存在困难,并且预测精度随研究体系的增大而不断变差。开发越来越精确的交换相关泛函是现代密度泛函理论发展的主线。近年来,引入未占轨道信息的新一代双杂化密度泛函方法的研究受到越来越多的关注。本篇综述回顾了这一领域已取得的一些进展;推导了双杂化泛函方法在Kohn-Sham密度泛函理论框架下的理论基础;依据各自特点,将现已提出的双杂化泛函划分为三种类别。本文着重测评了以XYG3为代表的一类双杂化泛函的表现。最后,就双杂化泛函的未来发展方向,提出了一些设想与建议。     

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.     

A computational study of bond dissociation enthalpies in haloethenes

Carbon–hydrogen bond dissociation enthalpies (BDEs) were computed for all haloethenes, C₂H_4−nX_n (n=0–3, X=F, Cl, Br, I), at the B3LYP/6-311+G(3df,2p) level using isodesmic reactions. It was found that C–H bond strengths in the monohaloethenes varied substantially, by as much as 18 kJ mol⁻¹, dependent upon the bond's stereochemical position relative to the halogen. BDEs in the dihaloethanes varied in the order CX₂CH–H>(E)-CHXCX–H>(Z)-CHXCX–H. Trends in the computed bond enthalpies were discussed and explained on the basis of relative steric repulsions and hyperconjugative delocalization interactions, as determined from Natural Bond Orbital analysis.     

Agreement between experiment and hybrid DFT calculations for O--H bond dissociation enthalpies in manganese complexes

Information on the accuracy of DFT functionals for redox reactions in transition metal systems is rather limited. To analyze the performance of some popular functionals for redox reactions in manganese systems, calculated O--H bond dissociation enthalpies for Mn-ligands in six different complexes are compared to experimental results. In this benchmark, B3LYP performs well with a mean absolute error of 3.0 kcal/mol. B98 gives similar results to B3LYP (error of 3.8 kcal/mol). B3LYP* gives lower O--H bond strengths than B3LYP and has a mean error of 5.0 kcal/mol. Compared to B98 and B3LYP, B3LYP* has an error trend for the manganese ligands that is more similar to the error for a free water molecule. The nonhybrid functional BLYP consistently and significantly underestimates the O--H bond strengths by approximately 20 kcal/mol. HCTH407 has a rather large mean error of 9.4 kcal/mol and shows no consistent trend. The results support the use of hybrid functionals and the present computational method for large model systems containing manganese. An example is the oxygen evolving complex in photosystem II where hybrid functionals predict the appearance of a Mn(IV)-oxyl radical before the O--O bond formation step.     

Substituent effects on O–H and S–H bond dissociation enthalpies of disubstituted phenols and thiophenols

The O–H and S–H homolytic bond dissociation enthalpies of a set of disubstituted phenols and thiophenols (NH₂, OH, CH₃, Cl, CF₃, and NO₂) have been computed by a density functional theory procedure with the 6‐311++G(d,p) basis set. A very good agreement between our results and available experimental ones is observed. The effect of substituents on structure, charges and BDEs are investigated and their correlation with Hammett parameters is studied. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

精确计算化学键解离能的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族和茶多酚族的键解离能, 并进一步讨论了抗氧化活性、自由基清除机理及其构效关系.     

Calculations of ionization energies and electron affinities for atoms and molecules: A comparative study with different methods

In the present work, we examined the performance of 36 density functionals, including the newly developed doubly hybrid density functional XYG3 (Y. Zhang, X. Xu, and W. A. Goddard III, Proc. Natl. Acad. Sci, USA, 2009, 106, 4963), to calculate ionization energies (IEs) and electron affinities (EAs). We used the well-established G2-1 set as reference, which contains 14 atoms and 24 molecules for IE, along with 7 atoms and 18 molecules for EA. XYG3 leads to mean absolute deviations (MADs) of 0.057 and 0.080 eV for IEs and EAs, respectively, using the basis set of 6-311 + G (3df,2p). In comparison with some other functionals, MADs for IEs are 0.109 (B2PLYP), 0.119 (M06-2X), 0.159 (X3LYP), 0.161 (PBE), 0.162 (B3LYP), 0.165 (PBE0), 0.173 (TPSS), 0.200 (BLYP), and 0.215 eV (LC-BLYP). MADs for EAs are 0.090 (X3LYP), 0.090 (B2PLYP), 0.102 (PBE), 0.103 (M06-2X), 0.104 (TPSS), 0.105 (BLYP), 0.106 (B3LYP), 0.126 (LC-BLYP), and 0.128 eV (PBE0).     

Accurate prediction of heats of formation by a combined method of B3LYP and neural network correction

Recently, we proposed the X1 method which combines the B3LYP/6‐311+G(3df,2p)//B3LYP/6‐311+G(d,p) method with a neural network correction for an accurate yet efficient prediction of heats of formation (Wu and Xu, J Chem Phys 2007, 127, 214105). In this contribution, we discuss in detail how to set up the X1 neural network. We give examples, showing how to apply the X1 method and how the applicability of X1 can be extended. The overall mean absolute deviation of the X1 method from experiment for the 488 heats of formation is 1.52 kcal/mol compared with 9.44 kcal/mol for the original B3LYP results. © 2008 Wiley Periodicals, Inc. J Comput Chem 2009     

The applicability of AM1 and PM3 semi-empirical methods for the study of N-H bond dissociation enthalpies and ionisation potentials of amine type antioxidants

A study of 40 para-substituted anilines is presented. These serve as model structures of amine type antioxidants. Molecules and their radical structures were studied using the AM1 and PM3 quantum chemical methods in order to calculate the N-H bond dissociation enthalpies (BDEs) and ionisation potentials (IPs) which are among the most important characteristics of antioxidants. Calculated BDEs were compared with available experimental values and the results of DFT calculations to ascertain the suitability of AM1 and PM3 methods for amine BDEs calculation. The results show that both methods reproduce experimental BDEs and DFT data satisfactorily. Comparison with experimental data shows that AM1 and PM3 methods overestimate the IP values. The results also indicate that dependences of N-H bond BDEs and IPs on Hammett constants of the substituents are linear.     

Simultaneous determination of enthalpies of dilution and dissociation of protocatechuic acid in aqueous salt solutions by calorimetric measurements

Protocatechuic acid, a major metabolite of antioxidant polyphenols found in green tea, has been shown to prevent carcinogenesis or antitumor growth in vitro and in vivo studies. The dilution enthalpies Δ_dilH_m, dissociation enthalpies Δ_disH_m and thermodynamic dissociation constants K_a for protocatechuic acid in aqueous NaCl or KCl solutions were simultaneously determined by mixing-flow microcalorimetry at T = 298.15 K. In order to verify the reliability of the fitted dissociation parameters, the values of dissociation enthalpies and thermodynamic dissociation constants were also determined by isothermal titration calorimetry and electrical conductivity methods. The K_a values obtained through the proposed method agree well with those reported in literatures and obtained by other techniques. Enthalpic interaction coefficients (h₂, h₃ and h₄) were computed according to the McMillan–Mayer model. The trends of h₂ and Δ_disH_m for protocatechuic acid with increasing salt molality in both the salt solutions were obtained. The different influence of KCl and NaCl on the values of h₂ and Δ_disH_m were discussed in terms of (solute + solute) and (solute + solvent) interactions. The results showed that it is possible to perform in a single calorimetric experiment the simultaneous determination of dilution enthalpies, dissociation enthalpies and constants in a given solvent.     

Substituent effects on O--H bond dissociation enthalpies and ionization potentials of catechols: a DFT study and its implications in the rational design of phenolic antioxidants and elucidation of structure-activity relationships for flavonoid antioxidants

Density functional theory (DFT) on B3LYP/6-31G(d,p) level was employed to investigate the substituent effects on O--H bond dissociation enthalpies (BDEs) and ionization potentials (IPs) of catechols. It was revealed that the ortho hydroxyl of catechol was effective for the reduction of the O--H BDE; however, the group had little influence on the IP. The para substituent effects upon O--H BDEs and IPs for catechols were roughly the same as those for monophenols, and this gave the catechol moiety more potential than monophenol to be used as a lead compound in rational design of phenolic antioxidants. In addition, the 1,4-pyrone effects on O--H BDEs of catecholic rings A or B of flavonoids were also investigated. Although 1,4-pyrone extended the conjugation system of flavonoids, it was not beneficial to reduce the O--H BDE as a result of its electron-withdrawing property. Thus, 1,4-pyrone was unlikely to be favorable to enhance the H-abstraction activity of flavonoids.     

Interaction of two functional groups through the benzene ring: theory and experiment

Energies of 132 benzene para bis-derivatives calculated within the framework of the density functional theory at the level B3LYP/6-311+G(d,p)//B3LYP/6-311+G(d,p) were used for correlations of two types. Correlation with the experimental enthalpies of formation clearly revealed that the published experimental data are generally not dependable and may be loaded with errors of more than 10 kJ mol(-1). On the other hand, the calculated relative energies are biased so that the interaction of the two substituents is systematically overestimated. This shortcoming was insignificant for our correlations of the second type, in which the interaction of substituents expressed in terms of isodesmic reactions was analyzed depending on the effects of inductive and resonance. The results depended strongly on the character of substituents. When one substituent is an electron donor and the other is an acceptor, the inductive-resonance model works and the classical resonance picture is adequate. With two acceptor substituents, this model is still acceptable with lower precision (as crossed conjugation), but with two donors it fails completely and may be acceptable only for a much restricted subclass of strong donors. Many correlations described in the literature must be viewed with great caution when they are based only on a relatively small number of data, in which substituents of different types are not represented in a comparable number.     

常见三氟甲基源释放三氟甲基自由基能力的理论计算研究

三氟甲基自由基源（·CF₃）是发展自由基三氟甲基化反应的基础.采用密度泛函理论（M06-2X）,系统研究35个常见三氟甲基源释放CF₃自由基的能力（TR^·DA：trifluoromethyl radical donor ability）及可能的途径.计算结果表明,35个三氟甲基源通过化学键均裂释放三氟甲基自由基所需能量跨度为-21.5至95.2 kcal·mol^-1.单电子转移、卤键和硫键给体（halogen/chalcogen-bond donor）可促进CF₃自由基释放.相关研究结果将有助于理解和设计新自由基三氟甲基化试剂和反应.     

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     

Basis set dependence of phosphate frequencies in density functional theory calculations

The addition of extravalence, polarization and diffuse functions, were studied in order to conclude how they affect the P? O stretching frequencies of several biological relevant phosphate molecules. The results show that the polarization and the diffuse functions have opposite effects on the frequencies: the polarization functions downshift while the diffuse functions upshift the frequencies. The effect of the valence functions was more difficult to interpret. The effect of the conductor‐like screening model (CPCM)‐continuum model was also studied. The results show that the CPCM‐continuum model has a substantial effect on the frequencies for these small molecules. The continuum model's efficiency is mainly due to its effect on the geometries and not on the frequencies. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012