Assessment of density functional theory methods for the computation of heats of formation and ionization potentials of systems containing third row transition metals

The performance of several different density functional theory (DFT) methods, including GGA, hybrid-GGA, meta-GGA, and hybrid-meta-GGA methods, have been assessed in terms of their ability to accurately compute both heats of formation and ionization potentials of systems containing third row transition metals. Two different basis sets were used in this study: 6-31G** and TZVP. It is found that the triple-zeta quality TZVP basis set generally produces the best results for both heats of formation and ionization potentials. One important observation made in this study is that the inclusion of exact exchange terms in DFT methods generally results in more consistently accurate results for both heats of formation and ionization potentials of transition metal systems. In general, DFT methods do not yield good ionization potential results for systems containing titanium or zinc. For heats of formation, it is found that the hybrid-meta-GGA functional, TPSS1KCIS, yields the best overall results when combined with the TZVP basis set, while PBE1PBE (hybrid-GGA) gives the best results for the 6-31G** basis. The hybrid-GGA functional, B3LYP, is found to produce the lowest overall errors for ionization potentials when combined with both 6-31G** and TZVP.     

An extended hybrid density functional (X3LYP) with improved descriptions of nonbond interactions and thermodynamic properties of molecular systems

We derive here the form for the exact exchange energy density for a density that decays with Gaussian-type behavior at long range. This functional is intermediate between the B88 and the PW91 exchange functionals. Using this modified functional to match the form expected for Gaussian densities, we propose the X3LYP extended functional. We find that X3LYP significantly outperforms Becke three parameter Lee-Yang-Parr (B3LYP) for describing van der Waals and hydrogen bond interactions, while performing slightly better than B3LYP for predicting heats of formation, ionization potentials, electron affinities, proton affinities, and total atomic energies as validated with the extended G2 set of atoms and molecules. Thus X3LYP greatly enlarges the field of applications for density functional theory. In particular the success of X3LYP in describing the water dimer (with R(e) and D(e) within the error bars of the most accurate determinations) makes it an excellent candidate for predicting accurate ligand-protein and ligand-DNA interactions.     

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.     

Electronic structures of 4d transition metal monoxides by density functional theory

Bond distances, dissociation energies, ionization potentials and electron affinities of 4d transition metal monoxides from YO to CdO and their positive and negative ions were studied by use of density functional methods B3LYP, BLYP, B3PW91, BPW91, B3P86, BP86, SVWN, MPW1PW91 and PBE1PBE. It was found that calculated properties are highly dependent on the functionals employed, especially for dissociation energy. For most neutral species, pure density functionals BLYP, BPW91 and BP86 have good performance in predicting dissociation energy than hybrid density functionals B3LYP, B3PW91 and B3P86. In addition, BLYP gives the largest bond distance compared with other density functional methods, while SVWN gives shortest bond distance, largest dissociation energy and electron affinity. For the ground state, the spin multiplicity of the charged species can be obtained by ± 1 of their corresponding neutral species.     

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     

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.     

Doubly hybrid density functional xDH-PBE0 from a parameter-free global hybrid model PBE0

Following the XYG3 model which uses orbitals and density from B3LYP, an empirical doubly hybrid (DH) functional is developed by using inputs from PBE0. This new functional, named xDH-PBE0, has been tested on a number of different molecular properties, including atomization energies, bond dissociation enthalpies, reaction barrier heights, and nonbonded interactions. From the results obtained, xDH-PBE0 not only displays a significant improvement with respect to the parent PBE0, but also shows a performance that is comparable to XYG3. Arguably, while PBE0 is a parameter-free global hybrid (GH) functional, the B3LYP GH functional contains eight fit parameters. From a more general point of view, the present work points out that reliable and general-purpose DHs can be obtained with a limited number of fit parameters.     

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     

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

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

Density functional study of structural and electronic properties of AlnN (1 ≤ n ≤ 12) clusters

Low‐lying equilibrium geometric structures of Al_nN (n = 1–12) clusters obtained by an all‐electron linear combination of atomic orbital approach, within spin‐polarized density functional theory, are reported. The binding energy, dissociation energy, and stability of these clusters are studied within the local spin density approximation (LSDA) and the three‐parameter hybrid generalized gradient approximation (GGA) due to Becke–Lee–Yang–Parr (B3LYP). Ionization potentials, electron affinities, hardness, and static dipole polarizabilities are calculated for the ground‐state structures within the GGA. It is observed that symmetric structures with the nitrogen atom occupying the internal position are lowest‐energy geometries. Generalized gradient approximation extends bond lengths as compared with the LSDA lengths. The odd–even oscillations in the dissociation energy, the second differences in energy, the highest occupied molecular orbital–lowest unoccupied molecular orbital (HOMO–LUMO) gaps, the ionization potential, the electron affinity, and the hardness are more pronounced within the GGA. The stability analysis based on the energies clearly shows the Al₇N cluster to be endowed with special stability. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

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     

Electronic structures of 3d-metal mononitrides

Bond distances, vibrational frequencies, electron affinities, ionization potentials, and dissociation energies of the title molecules in neutral, positively, and negatively charged ions were studied by use of density functional methods B3LYP, BLYP, BHLYP, BPW91, and B3PW91. The calculated results are compared with experiments and previous theoretical studies. It was found that the calculated properties are highly dependent on the functionals employed, in particular for the dissociation energy and vibrational frequency. For neutral species, pure density functional methods BLYP and BPW91 have relatively good performance in reproducing the experimental bond distance and vibrational frequency. For cations, hybrid exchange functional methods B3LYP and B3PW91 are good in predicting the dissociation energy. For both neutral and charged species, BHLYP tends to give smaller dissociation energy.     

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

Electronic Structures of 5d Transition Metal Monoxides by Density Functional Theory

Bond distances, vibrational frequencies, electron affinities, ionization potentials, and dissociation energies of the diatomic 5d transition metal (except La) monoxides and their positively and negatively charged ions were studied by use of density functional methods B3LYP, BLYP, B3PW91, BPW91, B3P86, BP86, MPW1PW91, PBE1PBE, and SVWN. Our calculation shows that for each individual species, the calculated properties are quite sensitive to the method used. Compared with hybrid density functional method B3PW91 (B3P86), pure density functional method BPW91 (BP86) gives longer bond distance (lower vibrational frequency) from HfO to PtO for neutral species, HfO⁺ to IrO⁺ for cationic species, and HfO⁻ to AuO⁻ for anionic species. While for B3LYP and BLYP, the trend was observed for cationic species from HfO⁺ to IrO⁺ and anionic species from HfO⁻ to AuO⁻ (except TaO⁻), but not for neutrals. Pure density function methods BLYP, BPW91, and BP86 give larger dissociation energy compared with hybrid density functional methods B3LYP, B3PW91, and B3P86. SVWN in most cases gives the smallest bond distance, while BLYP gives the largest value. MPW1PW91 and PBE1PBE show the same performance in predicting the spectroscopic constants. In addition, useful empirical criteria that one has obtained the ground states of a species and its ions are the spin multiplicities of a neutral and its single charged ions which differs by ±1.     

Performance of density functionals for first row transition metal systems

This article investigates the performance of five commonly used density functionals, B3LYP, BP86, PBE0, PBE, and BLYP, for studying diatomic molecules consisting of a first row transition metal bonded to H, F, Cl, Br, N, C, O, or S. Results have been compared with experiment wherever possible. Open-shell configurations are found more often in the order PBE0>B3LYP>PBE approximately BP86>BLYP. However, on average, 58 of 63 spins are correctly predicted by any functional, with only small differences. BP86 and PBE are slightly better for obtaining geometries, with errors of only 0.020 A. Hybrid functionals tend to overestimate bond lengths by a few picometers and underestimate bond strengths by favoring open shells. Nonhybrid functionals usually overestimate bond energies. All functionals exhibit similar errors in bond energies, between 42 and 53 kJmol. Late transition metals are found to be better modeled by hybrid functionals, whereas nonhybrid functionals tend to have less of a preference. There are systematic errors in predicting certain properties that could be remedied. BLYP performs the best for ionization potentials studied here, PBE0 the worst. In other cases, errors are similar. Finally, there is a clear tendency for hybrid functionals to give larger dipole moments than nonhybrid functionals. These observations may be helpful in choosing and improving existing functionals for tasks involving transition metals, and for designing new, improved functionals.