Accurate bond dissociation enthalpies by using doubly hybrid XYG3 functional

In this work, we examine the performance of XYG3, a newly developed doubly hybrid density functional (Zhang, Xu, and Goddard III, Proc Natl Acad Sci USA 2009, 106, 4963), to calculate covalent bond dissociation enthalpy (BDE). We use 5 atoms, 32 molecular radicals, and 116 closed-shell molecules to set up 142 bond dissociation reactions. For the total of 148 heats of formation (HOFs) and 142 BDEs, XYG3 leads to mean absolute deviations (MADs) of 1.45 and 1.87 kcal/mol, respectively. In comparison with some other functionals, MADs for HOFs are 2.31 (M06-2X), 2.98 (B2PLYP-D), 3.04 (BMK), 3.96 (B3LYP), 4.47 (B2PLYP), 5.42 (B2GP-PLYP), 6.46 (PBE0), and 29.93 kcal/mol (B3P86), and the corresponding errors for BDEs are 2.06 (M06-2X), 2.25 (BMK), 2.51 (B2PLYP-D), 2.89 (B2GP-PLYP), 3.30 (B3P86), 3.44 (B2PLYP), 3.87 (PBE0), and 6.14 kcal/mol (B3LYP).     

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).     

XYG3型双杂化密度泛函方法新进展：从能量到能量的解析梯度

本文回顾了现代密度泛函理论的基础,着重评述了XYG3型双杂化（XYG3typeofdoublyhybrid,xDH）泛函的最新进展,解析能量梯度的实现．XYG3是首个依照绝热途径理论建立的双杂化泛函,在具体实现上具有独特的构架．该类型泛函利用常用泛函（如B3LYP或PBE0等）作母泛函来进行自洽计算,以期获得更好的密度和轨道,然后将所得到的轨道和密度信息带入到xDH泛函中以得到最终能量．由于自洽泛函和最终能量泛函不同,因而在计算解析能量梯度时需要求解耦合微扰Kohn—Sham方程．在此基础上,还评述了xDH泛函在能量,尤其是构型优化方面的具体表现．测试的构型集包括以共价键键合的分子和非键相互作用体系的平衡结构,以及反应过渡态结构．结果表明,xDH双杂化泛函总体上给出了比母泛函更好的能量和几何构型．     

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.     

Performance of recently developed kinetic energy density functionals for the calculation of hydrogen binding strengths and hydrogen-bonded structures

 The accuracy of predicted hydrogen binding energies and equilibrium structures for a benchmark set of molecules is compared for some recently developed density functionals, Becke's three parameter hybrid method with the Lee, Yang, and Parr (LYP) correlation functional (B3LYP), Becke's half and half functional combined with the LYP correlation functional (BHLYP), Perdew, Burke and Ernzerhof functional (PBE), Van Voorhis, Scuseria exchange correlation functional (VSXC), the hybrid Perdew, Burke and Ernzerhof functional (PBE1PBE), and meta-generalized gradient approximation (meta-GGA). Overall, the hybrid functionals which contain a portion of Hartree–Fock exchange (B3LYP, BHLYP, and PBE1PBE) yield the most accurate results. The kinetic-energy-density-dependent functionals, VSXC and meta-GGA, are significantly less accurate. Received: 10 December 1999 / Accepted: 5 March 2000 / Published online: 21 June 2000     

Computational study of static first hyperpolarizability of donor-acceptor substituted (E)-benzaldehyde phenylhydrazone

Various hybrid functionals (B3LYP, B97-2, PBE0, BMK, BH&HLYP, CAM-B3LYP, and LC-ωPBE) implemented in density functional theory were applied to give estimate of static first hyperpolarizabilty (β(0)) of (E)-benzaldehyde phenylhydrazone designated as (E)-BPH. Against those of MP2 computations as a function of the underlying density functional, good agreement was obtained with the BH&HLYP and CAM-B3LYP functionals. The LC-ωPBE functional and the B3LYP, PBE0, B97-2, and BMK functionals underestimated and overestimated β(0), respectively. The basis set effect on the calculated β(0) was also investigated. It turned out that the 6-311+G(2d,p) basis set provided excellent converged value of β(0). On the basis of the calculated results, we investigated the substituent effect on β(0) of donor-acceptor (D-A) substituted (E)-BPH systematically by using the BH&HLYP and CAM-B3LYP computations with the 6-311+G(2d,p) basis set. We proposed a Zwitterion structure to explain the calculated trend in the substituent effect and the enhanced hyperpolarizability of type II compounds (A-(E)-BPH-D) than type I compounds (D-(E)-BPH-A). Natural bonding orbital analysis carried out at BH&HLYP/6-311+G(2d,p)//B3LYP/6-31G(2df,p) level of theory substantiated the claim.     

Performance of six functionals (LDA, PBE, PBESOL, B3LYP, PBE0, and WC1LYP) in the simulation of vibrational and dielectric properties of crystalline compounds. The case of forsterite Mg2SiO4

The performance of six different density functionals (LDA, PBE, PBESOL, B3LYP, PBE0, and WC1LYP) in describing the infrared spectrum of forsterite, a crystalline periodic system with orthorhombic unit cell (28 atoms in the primitive cell, Pbmn space group), is investigated by using the periodic ab initio CRYSTAL09 code and an all‐electron Gaussian‐type basis set. The transverse optical (TO) branches of the 35 IR active modes are evaluated at the equilibrium geometry together with the oscillator strengths and the high‐frequency dielectric tensor ?_∞. These quantities are essential to compute the dielectric function ?(ν), and then the reflectance spectrum R(ν), which is compared with experiment. It turns out that hybrid functionals perform better than LDA and GGA, in general; that B3LYP overperforms WC1LYP and, in turn, PBE0; that PBESOL is better than PBE; that LDA is the worst performing functional among the six under study. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Search for suitable approximation methods for fullerene structure and relative stability studies: case study with C50

Local density approximation (LDA), several popular general gradient approximation (GGA), hybrid module based density functional theoretical methods: SVWN, BLYP, PBE, HCTH, B3LYP, PBE1PBE, B1LYP, and BHandHLYP, and some nonstandard hybrid methods are applied in geometry prediction for C60 and C70. HCTH with 3-21G basis set is found to be one of the best methods for fullerene structural prediction. In the predictions of relative stability of C50 isomers, PM3 is an efficient method in the first step for sorting out the most stable isomers. HCTH with 3-21G predicts very good geometries for C50, similar to the performance of B3LYP6-31G(d). The gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital from the predictions of all the density functional theory methods has the following descending order: E(gap)(half-and-half hybrid)>E(gap)(B3LYP)>E(gap)(HCTH)(GGA)>E(gap)(SVWN)(LDA).     

Performance of conventional and range-separated hybrid density functionals in calculations of electronic circular dichroism spectra of transition metal complexes

A number of density functionals, including 'pure' (nonhybrid) functionals, global hybrids, and range-separated hybrids, were used to calculate the electronic circular dichroism (CD) spectra of 10 tris-bidentate transition metal complexes. The results are compared to one another and to experimental CD spectra, in an effort to illustrate the shortcomings of particular approximations in time-dependent density functional theory (TDDFT). The use of an origin invariant formalism to calculate magnetic transition dipole moments with the help of gauge-including atomic orbitals (GIAOs) is also investigated. With valence basis sets of moderate flexibility, good agreement between GIAO results and rotatory strengths calculated from the dipole-velocity representation is obtained for selected test cases. Empirically broadened vertical CD spectra calculated with the global hybrid functionals B3LYP and PBE0 are found to agree overall the best with experimental CD spectra.     

Optical rotation calculated with time-dependent density functional theory: the OR45 benchmark

Time-dependent density functional theory (TDDFT) computations are performed for 42 organic molecules and three transition metal complexes, with experimental molar optical rotations ranging from 2 to 2 × 10(4) deg cm(2) dmol(-1). The performances of the global hybrid functionals B3LYP, PBE0, and BHLYP, and of the range-separated functionals CAM-B3LYP and LC-PBE0 (the latter being fully long-range corrected), are investigated. The performance of different basis sets is studied. When compared to liquid-phase experimental data, the range-separated functionals do, on average, not perform better than B3LYP and PBE0. Median relative deviations between calculations and experiment range from 25 to 29%. A basis set recently proposed for optical rotation calculations (LPol-ds) on average does not give improved results compared to aug-cc-pVDZ in TDDFT calculations with B3LYP. Individual cases are discussed in some detail, among them norbornenone for which the LC-PBE0 functional produced an optical rotation that is close to available data from coupled-cluster calculations, but significantly smaller in magnitude than the liquid-phase experimental value. Range-separated functionals and BHLYP perform well for helicenes and helicene derivatives. Metal complexes pose a challenge to first-principles calculations of optical rotation.     

Is the spin‐orbit coupling important in the prediction of the 51V hyperfine coupling constants of VIVO2+ species? ORCA versus Gaussian performance and biological applications

Density functional theory calculations of the (51)V hyperfine coupling (HFC) tensor A, have been completed for eighteen V(IV)O(2+) complexes with different donor set, electric charge and coordination geometry. A tensor was calculated with ORCA software with several functionals and basis sets taking into account the spin-orbit coupling contribution. The results were compared with those obtained with Gaussian 03 software using the half-and-half functional BHandHLYP and 6-311g(d,p) basis set. The order of accuracy of the functionals in the prediction of A(iso), A(z) and dipolar term A(z,anis) is BHandHLYP > PBE0 > B3PW > TPSSh > B3LYP > BP86 > VWN5 (for A(iso)), BHandHLYP > PBE0 > B3PW > TPSSh > B3LYP > BP86 > VWN5 (for A(z)), B3LYP > PBE0 ～ B3PW ～ BHandHLYP > TPSSh > BP86 ～ VWN5 (for A(z,anis)). The good agreement in the prediction of A(z) with BHandHLYP is due to a compensation between the overestimation of A(iso) and underestimation of A(z,anis) (A(z) = A(iso) + A(z,anis)), whereas among the hybrid functionals PBE0 performs better than the other ones. BHandHLYP functional and Gaussian software are recommended when the V(IV)O(2+) species contains only V-O and/or V-N bonds, whereas PBE0 functional and ORCA software for V(IV)O(2+) complexes with one or more V-S bonds. Finally, the application of these methods to the coordination environment of V(IV)O(2+) ion in V-proteins, like vanadyl-substituted insulin, carbonic anhydrase, collagen and S-adenosylmethionine synthetase, was discussed.     

Semiempirical hybrid functional with improved performance in an extensive chemical assessment

It is demonstrated that there is still scope for improvement in the quality of conventional, semiempirical hybrid exchange-correlation functionals in density-functional theory. A new functional, denoted B97-3, is determined from a fit to eight chemical properties (316 data points). For a series of 25 chemical assessments (850 data points) including 17 assessments and 10 chemical properties absent from the fitting data, B97-3 provides the lowest or joint-lowest mean absolute error on 15 occasions, compared to 6, 5, and 4 occasions for B3LYP, PBE0, and B97-2, respectively [A. D. Becke, J. Chem. Phys. 98, 5648 (1993); M. Ernzerhof and G. E. Scuseria, J. Chem. Phys. 110, 5029 (1999); C. Adamo and V. Barone, J. Chem. Phys. 110, 6158 (1999); P. J. Wilson, T. J. Bradley, and D. J. Tozer, J. Chem. Phys. 115, 9233 (2001)]. Mean absolute errors from B97-3 are, on average, 21%, 18%, and 12% smaller than from these three functionals. The most notable improvements are obtained for classical reaction barriers, where the error reductions are 60%, 54%, and 27%.     

Systematic investigation on the geometric dependence of the calculated nuclear magnetic shielding constants

We perform a systematic examination on the dependence of the calculated nuclear magnetic shielding constants on the chosen geometry for a selective set of density functional methods of B3LYP, PBE0, and OPBE. We find that the OPBE exchange-correlation functional performs remarkably well when either the optimized geometries or the experimental geometries are used. The popular B3LYP and PBE0 functionals have a clear tendency of deshielding, giving shieldings that are usually too low and shifts that are usually too high, at the experimental geometries. Combined with the Hartree-Fock geometries, however, much improved magnetic constants are obtained for B3LYP and PBE0, due to the compensation effect from the systematic underestimation of bond lengths by the Hartree-Fock method.     

一些密度泛函方法预测电子亲和势

选取了杂化泛函B3LYP, B3PW91, O3LYP, PBE0, 以及与之相对应的GGA泛函BLYP, BPW91, OLYP和PBE, 还选取了能更好地兼顾强相互作用和弱相互作用的X3LYP泛函和在预测NMR的化学位移有较好表现的OPBE泛函, 以及两种meta-GGA泛函VSXC和TPSS, 共12种泛函, 详细地考察了这些泛函在预测EA方面的准确性.     

XYG3 and XYGJ-OS performances for noncovalent binding energies relevant to biomolecular structures

XYG3 and XYGJ-OS doubly hybrid density functionals (DHDFs) are applied to calculate intermolecular interaction energies of biological relevance. The results are compared with those of the other types of DHDFs (MC3BB and B2PLYP), as well as their dispersion corrected variants (MC3BB-D B2PLYP-D, and B2PLYP-D3), in addition to those obtained from B3LYP, B3LYP-D and B3LYP-D3. The reference data are taken from the S22, S22x5 and JSCH-2005 benchmark sets consisting mainly of DNA base pairs and amino acid pairs. The overall good agreement with the reference values of the extrapolated CCSD(T) complete basis set limit suggests that the XYG3 and XYGJ-OS functionals are valuable tools for applications in biochemistry.