Density functionals for inorganometallic and organometallic chemistry

We present a database of 21 bond dissociation energies for breaking metal-ligand bonds. The molecules in the metal-ligand bond energy database are AgH, CoH, CoO+, CoOH+, CrCH3+, CuOH2+, FeH, Fe(CO)5, FeO, FeS, LiCl, LiO, MgO, MnCH3NiCH2+, Ni(CO)4, RhC, VCO+, VO, and VS. We have also created databases of metal-ligand bond lengths and atomic ionization potentials. The molecules used for bond lengths are AgH, BeO, CoH, CoO+, FeH, FeO, FeS, LiCl, LiO, MgO, RhC, VO, and VS and the ionization potentials are for the following atoms: C, Co, Cr, Cu, Ni, O, and V. The data were chosen based on their diversity and expected reliability, and they are used along with three previously developed databases (transition metal dimer bond energies and bond lengths and main-group molecular atomization energies) for assessing the accuracy of several kinds of density functionals. In particular, we report tests for 42 previously defined functionals: 2 local spin density approximation (LSDA) functionals, 14 generalized gradient approximation (GGA) methods, 13 hybrid GGA methods, 7 meta GGA methods, and 8 hybrid meta GGA methods. In addition to these functionals, we also examine the effectiveness of scaling the correlation energy by testing 13 functionals with scaled or no gradient-corrected correlation energy, and we find that functionals of this kind are more accurate for metal-metal and metal-ligand bonds than any of the functionals already in the literature. We also present a readjusted GGA and a hybrid GGA with parameters adjusted for metals. When we consider these 57 functionals for metal-ligand and metal-metal bond energies simultaneously with main-group atomization energies, atomic ionization potentials, and bond lengths we find that the most accurate functional is G96LYP, followed closely by MPWLYP1M (new in this article), XLYP, BLYP, and MOHLYP (also new in this article). Four of these five functionals have no Hartree-Fock exchange, and the other has only 5%. As a byproduct of this work we introduce a convenient diagnostic, called the B1 diagnostic, for ascertaining the multireference character in a bond.     

Energy landscapes of nucleophilic substitution reactions: a comparison of density functional theory and coupled cluster methods

We have carried out a detailed evaluation of the performance of all classes of density functional theory (DFT) for describing the potential energy surface (PES) of a wide range of nucleophilic substitution (SN2) reactions involving, amongst others, nucleophilic attack at carbon, nitrogen, silicon, and sulfur. In particular, we investigate the ability of the local density approximation (LDA), generalized gradient approximation (GGA), meta-GGA as well as hybrid DFT to reproduce high-level coupled cluster (CCSD(T)) benchmarks that are close to the basis set limit. The most accurate GGA, meta-GGA, and hybrid functionals yield mean absolute deviations of about 2 kcal/mol relative to the coupled cluster data, for reactant complexation, central barriers, overall barriers as well as reaction energies. For the three nonlocal DFT classes, the best functionals are found to be OPBE (GGA), OLAP3 (meta-GGA), and mPBE0KCIS (hybrid DFT). The popular B3LYP functional is not bad but performs significantly worse than the best GGA functionals. Furthermore, we have compared the geometries from several density functionals with the reference CCSD(T) data. The same GGA functionals that perform best for the energies (OPBE, OLYP), also perform best for the geometries with average absolute deviations in bond lengths of 0.06 A and 0.6 degrees, even better than the best meta-GGA and hybrid functionals. In view of the reduced computational effort of GGAs with respect to meta-GGAs and hybrid functionals, let alone coupled cluster, we recommend the use of accurate GGAs such as OPBE or OLYP for the study of SN2 reactions.     

Is combining meta‐GGA correlation functionals with the OPTX exchange functional useful?

The recent generalized gradient approximation (GGA) density functional OCS1 of Handy and Cohen is implemented in the deMon code and tested on a carefully selected set of problems. OCS1 is found to be accurate for molecular atomization energies, transition metal–ligand bonds, and systems with intramolecular hydrogen bonds. However, OCS1 encounters problems for systems with intermolecular hydrogen bonds. It also tends to elongate bond lengths systematically, and sometimes significantly. The OPTX exchange is combined with three meta‐GGA correlation functionals, Lap3, τ1, and τ2, the latter reported for the first time. The new meta‐GGA scheme OPTX exchange plus τ2 correlation called Oτ2 yields improved molecular geometries, NMR shielding constants, and an improved barrier height for the H+H₂ reaction. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

How reliable is DFT in predicting relative energies of polycyclic aromatic hydrocarbon isomers? comparison of functionals from different rungs of jacob's ladder

Density functional theory (DFT) is the only quantum‐chemical avenue for calculating thermochemical/kinetic properties of large polycyclic aromatic hydrocarbons (PAHs) such as graphene nanoflakes. Using CCSD(T)/CBS PAH isomerization energies, we find that all generalized gradient approximation (GGA) and meta GGA DFT functionals have severe difficulties in describing isomerization energies in PAHs. The poor performance of these functionals is demonstrated by the following root‐mean‐square deviations (RMSDs) obtained for a database of C₁₄H₁₀ and C₁₈H₁₂ isomerization energies. The RMSDs for the GGAs range between 6.0 (BP86‐D3) and 23.0 (SOGGA11) and for the meta GGAs they range between 3.5 (MN12‐L) and 11.9 (τ‐HCTH) kJ mol⁻¹. These functionals (including the dispersion‐corrected methods) systematically and significantly underestimate the isomerization energies. A consequence of this behavior is that they all predict that chrysene (rather than triphenylene) is the most stable C₁₈H₁₂ isomer. A general improvement in performance is observed along the rungs of Jacob's Ladder; however, only a handful of functionals from rung four give good performance for PAH isomerization energies. These include functionals with high percentages (40–50%) of exact Hartree–Fock exchange such as the hybrid GGA SOGGA11‐X (RMSD = 1.7 kJ mol⁻¹) and the hybrid‐meta GGA BMK (RMSD = 1.3 kJ mol⁻¹). Alternatively, the inclusion of lower percentages (20–30%) of exact exchange in conjunction with an empirical dispersion correction results in good performance. For example, the hybrid GGA PBE0‐D3 attains an RMSD of 1.5 kJ mol⁻¹, and the hybrid‐meta GGAs PW6B95‐D3 and B1B95‐D3 result in RMSDs below 1 kJ mol⁻¹. © 2016 Wiley Periodicals, Inc.     

Construction of a generalized gradient approximation by restoring the density-gradient expansion and enforcing a tight Lieb-Oxford bound

Recently, a generalized gradient approximation (GGA) to the density functional, called PBEsol, was optimized (one parameter) against the jellium-surface exchange-correlation energies, and this, in conjunction with changing another parameter to restore the first-principles gradient expansion for exchange, was sufficient to yield accurate lattice constants of solids. Here, we construct a new GGA that has no empirical parameters, that satisfies one more exact constraint than PBEsol, and that performs 20% better for the lattice constants of 18 previously studied solids, although it does not improve on PBEsol for molecular atomization energies (a property that neither functional was designed for). The new GGA is exact through second order, and it is called the second-order generalized gradient approximation (SOGGA). The SOGGA functional also differs from other GGAs in that it enforces a tighter Lieb-Oxford bound. SOGGA and other functionals are compared to a diverse set of lattice constants, bond distances, and energetic quantities for solids and molecules (this includes the first test of the M06-L meta-GGA for solid-state properties). We find that classifying density functionals in terms of the magnitude mu of the second-order coefficient of the density gradient expansion of the exchange functional not only correlates their behavior for predicting lattice constants of solids versus their behavior for predicting small-molecule atomization energies, as pointed out by Perdew and co-workers [Phys. Rev. Lett. 100, 134606 (2008); Perdew ibid. 80, 891 (1998)], but also correlates their behavior for cohesive energies of solids, reaction barriers heights, and nonhydrogenic bond distances in small molecules.     

The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals

We present two new hybrid meta exchange- correlation functionals, called M06 and M06-2X. The M06 functional is parametrized including both transition metals and nonmetals, whereas the M06-2X functional is a high-nonlocality functional with double the amount of nonlocal exchange (2X), and it is parametrized only for nonmetals.The functionals, along with the previously published M06-L local functional and the M06-HF full-Hartree–Fock functionals, constitute the M06 suite of complementary functionals. We assess these four functionals by comparing their performance to that of 12 other functionals and Hartree–Fock theory for 403 energetic data in 29 diverse databases, including ten databases for thermochemistry, four databases for kinetics, eight databases for noncovalent interactions, three databases for transition metal bonding, one database for metal atom excitation energies, and three databases for molecular excitation energies. We also illustrate the performance of these 17 methods for three databases containing 40 bond lengths and for databases containing 38 vibrational frequencies and 15 vibrational zero point energies. We recommend the M06-2X functional for applications involving main-group thermochemistry, kinetics, noncovalent interactions, and electronic excitation energies to valence and Rydberg states. We recommend the M06 functional for application in organometallic and inorganometallic chemistry and for noncovalent interactions. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Contribution to the Mark S. Gordon 65th Birthday Festschrift Issue.     

Ground- and excited-state diatomic bond lengths, vibrational levels, and potential-energy curves from conventional and localized Hartree-Fock-based density-functional theory

Ground- and excited-state diatomic bond lengths, vibrational levels, and potential-energy curves are determined using conventional and localized Hartree-Fock (LHF)-based density-functional theory. Exchange only and hybrid functionals (with various fractions of exchange) are considered, together with a standard generalized gradient approximation (GGA). Ground-state bond lengths and vibrational wave numbers are relatively insensitive to whether orbital exchange is treated using the conventional or LHF approach. Excited-state calculations are much more sensitive. For a standard fraction of orbital exchange, N2 and CO vertical excitation energies at experimental bond lengths are accurately described by both conventional and LHF-based approaches, providing an asymptotic correction is present. Excited-state bond lengths and vibrational levels are more accurate with the conventional approach. The best quality, however, is obtained with an asymptotically corrected GGA functional. For the ground and lowest four singlet excited states, the GGA mean absolute errors in bond lengths are 0.006 A (0.5%) and 0.011 A (0.8%) for N2 and CO, respectively. Mean absolute errors in fundamental vibrational wavenumbers are 49 cm(-1) (2.7%) and 68 cm(-1) (5.0%), respectively. The GGA potential-energy curves are compared with near-exact Rydberg-Klein-Rees curves. Agreement is very good for the ground and first excited state, but deteriorates for the higher states.     

A new database and benchmark of the bond energies of noble-gas-containing molecules

We have developed a new database of structures and bond energies of 59 noble-gas-containing molecules. The structures were calculated by CCSD(T)/aug-cc-pVTZ methods and the bond energies were obtained using the CCSD(T)/complete basis set method. Many wavefunction-based and density functional theory methods have been benchmarked against the 59 accurate bond energies. Our results show that the MPW1B95, B2GP-PLYP, and DSD-BLYP functionals with the aug-cc-pVTZ basis set excel in predicting the bond energies of noble-gas molecules with mean unsigned errors (MUEs) of 2.0 to 2.1 kcal/mol. When combinations of Dunning's basis sets are used, the MPW1B95, B2GP-PLYP, DSD-BLYP, and BMK functionals give significantly lower MUEs of 1.6 to 1.9 kcal/mol. Doubly hybrid methods using B2GP-PLYP and DSD-BLYP functionals and MP2 calculation also provide satisfactory accuracy with MUEs of 1.4 to 1.5 kcal/mol. If the Ng bond energies and the total atomization energies of a group of 109 main-group molecules are considered at the same time, the MPW1B95/aug-cc-pVTZ single-level method (MUE = 2.7 kcal/mol) and the B2GP-PLYP and DSD-PLYP functionals with combinations of basis sets or using the doubly hybrid method (MUEs = 1.9-2.2 kcal/mol) give the overall best result.     

Double-hybrid density functional theory for excited electronic states of molecules

Double-hybrid density functionals are based on a mixing of standard generalized gradient approximations (GGAs) for exchange and correlation with Hartree-Fock (HF) exchange and a perturbative second-order correlation part (PT2) that is obtained from the Kohn-Sham (GGA) orbitals and eigenvalues. This virtual orbital-dependent functional (dubbed B2PLYP) contains only two empirical parameters that describe the mixture of HF and GGA exchange (ax) and of the PT2 and GGA correlation (ac), respectively. Extensive testing has recently demonstrated the outstanding accuracy of this approach for various ground state problems in general chemistry applications. The method is extended here without any further empirical adjustments to electronically excited states in the framework of time-dependent density functional theory (TD-DFT) or the closely related Tamm-Dancoff approximation (TDA-DFT). In complete analogy to the ground state treatment, a scaled second-order perturbation correction to configuration interaction with singles (CIS(D)) wave functions developed some years ago by Head-Gordon et al. [Chem. Phys. Lett. 219, 21 (1994)] is computed on the basis of density functional data and added to the TD(A)-DFTGGA excitation energy. The method is implemented by applying the resolution of the identity approximation and the efficiency of the code is discussed. Extensive tests for a wide variety of molecules and excited states (of singlet, triplet, and doublet multiplicities) including electronic spectra are presented. In general, rather accurate excitation energies (deviations from reference data typically <0.2 eV) are obtained that are mostly better than those from standard functionals. Still, systematic errors are obtained for Rydberg (too low on average by about 0.3 eV) and charge-transfer transitions but due to the relatively large ax parameter (0.53), B2PLYP outperforms most other functionals in this respect. Compared to conventional HF-based CIS(D), the method is more robust in electronically complex situations due to the implicit account of static correlation effects by the GGA parts. The (D) correction often works in the right direction and compensates for the overestimation of the transition energy at the TD level due to the elevated fraction of HF exchange in the hybrid GGA part. Finally, the limitations of the method are discussed for challenging systems such as transition metal complexes, cyanine dyes, and multireference cases.     

Simulating periodic trends in the structure and catalytic activity of coinage metal nanoribbons

We present a systematic density functional theory (DFT) study of the structure and catalytic activity of group 10 (Ni, Pd, Pt) and group 11 (Cu, Ag, Au) coinage metal nanoribbons. These infinite, periodic, quasi‐one‐dimensional structures are conceptually important as intermediates between small metal clusters and close‐packed metal surfaces, and have been shown experimentally to be practical catalysts. We find that nanoribbons have significantly higher predicted H₂ dissociation activity than close‐packed metal surfaces consistent with their lower coordination numbers. Computed periodic trends are reasonable, with late transition states and low barriers for H₂ dissociation over late group 10 nanoribbons, suggesting their promise as practical catalysts. These trends are consistent with the isolated nanoribbons' computed molecular electrostatic potentials. Calculations also predict nearly linear Brønsted–Evans–Polanyi relationships between the nanoribbons' H₂ dissociation energies and dissociation barriers. We also test new meta‐generalized gradient approximation (GGA) and hybrid DFT approximations for H₂ dissociation over these nanoribbons. These new functionals increase the (generally underestimated) dissociation barriers predicted by standard GGAs, motivating their continued application in surface chemistry. © 2015 Wiley Periodicals, Inc.     

The performance of semilocal and hybrid density functionals in 3d transition-metal chemistry

We investigate the performance of contemporary semilocal and hybrid density functionals for bond energetics, structures, dipole moments, and harmonic frequencies of 3d transition-metal (TM) compounds by comparison with gas-phase experiments. Special attention is given to the nonempirical metageneralized gradient approximation (meta-GGA) of Tao, Perdew, Staroverov, and Scuseria (TPSS) [Phys. Rev. Lett. 91, 146401 (2003)], which has been implemented in TURBOMOLE for the present work. Trends and error patterns for classes of homologous compounds are analyzed, including dimers, monohydrides, mononitrides, monoxides, monofluorides, polyatomic oxides and halogenides, carbonyls, and complexes with organic pi ligands such as benzene and cyclopentadienyl. Weakly bound systems such as Ca(2), Mn(2), and Zn(2) are discussed. We propose a reference set of reaction energies for benchmark purposes. Our all-electron results with quadruple zeta valence basis sets validate semilocal density-functional theory as the workhorse of computational TM chemistry. Typical errors in bond energies are substantially larger than in (organic) main group chemistry, however. The Becke-Perdew'86 [Phys. Rev. A 38, 3098 (1988); Phys. Rev. B 33, 8822 (1986)] GGA and the TPSS meta-GGA have the best price/performance ratio, while the TPSS hybrid functional achieves a slightly lower mean absolute error in bond energies. The popular Becke three-parameter hybrid B3LYP underbinds significantly and tends to overestimate bond distances; we give a possible explanation for this. We further show that hybrid mixing does not reduce the width of the error distribution on our reference set. The error of a functional for the s-d transfer energy of a TM atom does not predict its error for TM bond energies and bond lengths. For semilocal functionals, self-interaction error in one- and three-electron bonds appears to be a major source of error in TM reaction energies. Nevertheless, TPSS predicts the correct ground-state symmetry in the vast majority of cases and rarely fails qualitatively. This further confirms TPSS as a general purpose functional that works throughout the periodic table. We also give workstation timing comparisons for the 645-atom protein crambin.     

Improved density functionals for water

The accuracy of existing density functional methods for describing the noncovalent interaction energies in small water clusters is investigated by testing 25 density functionals against a data set of 28 water dimers and 8 water trimers whose structures are taken from the literature and from simulations. The most accurate functionals are found to be PW6B95 with a mean unsigned error of 0.13 kcal/mol and MPWB1K and B98 with mean unsigned errors of 0.15 kcal/mol; the best functional with no Hartree-Fock exchange is mPWLYP, which is a GGA with a mean unsigned error of 0.28 kcal/mol. In comparison, the most popular GGA functionals, PBE and BLYP, have mean unsigned errors of 0.52 and 1.03 kcal/mol, respectively. Since GGAs are very cost efficient for both condensed-phase simulations and electronic structure calculations on large systems, we optimized four new GGAs for water. The best of these, PBE1W and MPWLYP1W, have mean unsigned errors of 0.12 and 0.17 kcal/mol, respectively. These new functionals are well suited for use in condensed-phase simulations of water and ice.     

Highly accurate first-principles benchmark data sets for the parametrization and validation of density functional and other approximate methods. Derivation of a robust, generally applicable, double-hybrid functional for thermochemistry and thermochemical kinetics

We present a number of near-exact, nonrelativistic, Born-Oppenheimer reference data sets for the parametrization of more approximate methods (such as DFT functionals). The data were obtained by means of the W4 ab initio computational thermochemistry protocol, which has a 95% confidence interval well below 1 kJ/mol. Our data sets include W4-08, which are total atomization energies of over 100 small molecules that cover varying degrees of nondynamical correlations, and DBH24-W4, which are W4 theory values for Truhlar's set of 24 representative barrier heights. The usual procedure of comparing calculated DFT values with experimental atomization energies is hampered by comparatively large experimental uncertainties in many experimental values and compounds errors due to deficiencies in the DFT functional with those resulting from neglect of relativity and finite nuclear mass. Comparison with accurate, explicitly nonrelativistic, ab initio data avoids these issues. We then proceed to explore the performance of B2x-PLYP-type double hybrid functionals for atomization energies and barrier heights. We find that the optimum hybrids for hydrogen-transfer reactions, heavy-atoms transfers, nucleophilic substitutions, and unimolecular and recombination reactions are quite different from one another: out of these subsets, the heavy-atom transfer reactions are by far the most sensitive to the percentages of Hartree-Fock-type exchange y and MP2-type correlation x in an (x, y) double hybrid. The (42,72) hybrid B2K-PLYP, as reported in a preliminary communication, represents the best compromise between thermochemistry and hydrogen-transfer barriers, while also yielding excellent performance for nucleophilic substitutions. By optimizing for best overall performance on both thermochemistry and the DBH24-W4 data set, however, we find a new (36,65) hybrid which we term B2GP-PLYP. At a slight expense in performance for hydrogen-transfer barrier heights and nucleophilic substitutions, we obtain substantially better performance for the other reaction types. Although both B2K-PLYP and B2GP-PLYP are capable of 2 kcal/mol quality thermochemistry, B2GP-PLYP appears to be the more robust toward nondynamical correlation and strongly polar character. We additionally find that double-hybrid functionals display excellent performance for such problems as hydrogen bonding, prototype late transition metal reactions, pericyclic reactions, prototype cumulene-polyacetylene system, and weak interactions.     

Assessment of new meta and hybrid meta density functionals for predicting the geometry and binding energy of a challenging system: the dimer of H2S and benzene

Noncovalent interactions of a hydrogen bond donor with an aromatic pi system present a challenge for density functional theory, and most density functionals do not perform well for this kind of interaction. Here we test seven recent density functionals from our research group, along with the popular B3LYP functional, for the dimer of H 2S with benzene. The functionals considered include the four new meta and hybrid meta density functionals of the M06 suite, three slightly older hybrid meta functionals, and the B3LYP hybrid functional, and they were tested for their abilities to predict the dissociation energies of three conformations of the H 2S-benzene dimer and to reproduce the key geometric parameters of the equilibrium conformation of this dimer. All of the functionals tested except B3LYP correctly predict which of the three conformations of the dimer is the most stable. The functionals that are best able to reproduce the geometry of the equilibrium conformation of the dimer with a polarized triple-zeta basis set are M06-L, PWB6K, and MPWB1K, each having a mean unsigned relative error across the two experimentally verifiable geometric parameters of only 8%. The success of M06-L is very encouraging because it is a local functional, which reduces the cost for large simulations. The M05-2X functional yields the most accurate binding energy of a conformation of the dimer for which a binding energy calculated at the CCSD(T) level of theory is available; M05-2X gives a binding energy for the system with a difference of merely 0.02 kcal/mol from that obtained by the CCSD(T) calculation. The M06 functional performs well in both categories by yielding a good representation of the geometry of the equilibrium structure and by giving a binding energy that is only 0.19 kcal/mol different from that calculated by CCSD(T). We conclude that the new generation of density functionals should be useful for a variety of problems in biochemistry and materials where aromatic functional groups can serve as hydrogen bond acceptors.     

Role of nonlocal exchange in molecular crystals: The case of two proton‐ordered phases of ice

We present a periodic density functional theory investigation of twoproton‐ordered phases of ice. Their equilibrium lattice parameters,relative stabilities, formation energies, and densities of states havebeen evaluated. Nine exchange‐correlation functionals, representativeof the generalized gradient approximation (GGA), global hybrids,range‐separated hybrids, meta‐GGA, and hybrid meta‐GGA families havebeen taken into account, considering two oxygen basis sets. Althoughthe hydrogen‐bond network of ice is well reproduced at the B3LYP,M06‐L, or LC‐ wPBE levels, formation energies are only correctlyevaluated with the two former functionals. Band gaps on the other handare only quantitatively reproduced at the B3LYP level. These resultsindicate that this last functional, a de facto reference formolecular calculations, gives in average the most accurate results forthe considered ice properties. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Comparison of DFT methods for molecular orbital eigenvalue calculations

We report how closely the Kohn-Sham highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) eigenvalues of 11 density functional theory (DFT) functionals, respectively, correspond to the negative ionization potentials (-IPs) and electron affinities (EAs) of a test set of molecules. We also report how accurately the HOMO-LUMO gaps of these methods predict the lowest excitation energies using both time-independent and time-dependent DFT (TD-DFT). The 11 DFT functionals include the local spin density approximation (LSDA), five generalized gradient approximation (GGA) functionals, three hybrid GGA functionals, one hybrid functional, and one hybrid meta GGA functional. We find that the HOMO eigenvalues predicted by KMLYP, BH&HLYP, B3LYP, PW91, PBE, and BLYP predict the -IPs with average absolute errors of 0.73, 1.48, 3.10, 4.27, 4.33, and 4.41 eV, respectively. The LUMOs of all functionals fail to accurately predict the EAs. Although the GGA functionals inaccurately predict both the HOMO and LUMO eigenvalues, they predict the HOMO-LUMO gap relatively accurately (approximately 0.73 eV). On the other hand, the LUMO eigenvalues of the hybrid functionals fail to predict the EA to the extent that they include HF exchange, although increasing HF exchange improves the correspondence between the HOMO eigenvalue and -IP so that the HOMO-LUMO gaps are inaccurately predicted by hybrid DFT functionals. We find that TD-DFT with all functionals accurately predicts the HOMO-LUMO gaps. A linear correlation between the calculated HOMO eigenvalue and the experimental -IP and calculated HOMO-LUMO gap and experimental lowest excitation energy enables us to derive a simple correction formula.