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.     

DFT study of the structural characteristics of the yttrium(3+) aqua ion

Density functional theory (DFT) using SVWN5, B3LYP, B3P86, O3LYP, B3PW91, B1LYP, B971, MPW1PW91, PBE1PBE, BHandH, and BHandHLYP density functionals was employed to study the structural characteristics of the Y(H₂O) ₈ ³⁺ yttrium aqua ion. The nonlocal hybrid GGA functionals show worse predictive ability in structural calculations of the Y(H₂O) ₈ ³⁺ aqua ion compared to the relatively simple combined functional BHandH and to the simplest SVWN5 functional in LSDA theory.     

Analysis of the performance of DFT functionals in the study of light emission by oxyluciferin analogs

Firefly bioluminescence has attracted much attention from the research community due to its most intriguing multicolor bioluminescence. Computational studies have been fundamental to the clarification of this phenomenon. However, no clear demonstration of the accuracy of the theoretical methods used in this field of research has been provided. In this work, we have used luciferin and dehydroluciferin in water as a model for the calculation of the emission wavelength of oxyluciferin and analogs. This model was used in a density functional theory (DFT) benchmarking in which 24 functionals were tested. A non‐negligible “case to case” dependence was observed. However in general, the better classes of the DFT functionals are generalized gradient approximation (GGA) and meta GGA. The generalized gradient exchange functional TPSSVWN5, the hybrid‐meta GGA functional mPWKCIS, and the hybrid GGA functional MPWLYP1M also present good results. For the contrary, the popular hybrid GGA functionals present a worse performance. © 2012 Wiley Periodicals, Inc.     

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

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     

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.     

Broken symmetry density functional study of a mixed‐valence unsymmetrical dinuclear iron complex

Density functional theory (DFT) calculations with different exchange‐correlation functionals were performed for a mixed valence Fe(II)/Fe(III) binuclear complex with μ‐methoxo and two μ‐carboxylate bridging ligands, (1) with geometry optimizations being performed for all possible spin multiplicities (M_S = 2, 4, 6, 8, and 10). Within the exchange‐correlation functionals studied, only the hybrid GGA functionals B3P and B3LYP and also the pure GGA functional RPBE, predicts the geometry with high spin (S = 9/2) to be more stable than the geometry with low spin state (S = 1/2) by 20 kcal/mol, in agreement with the experimental findings. These functionals also predict the same stability order for the different spin states, being M_S = 10>8>6>2>4. The meta‐GGA functionals TPSS and TPSSh and also the pure GGA functionals BLYP and BP86 predict different stability orders. The computed average EPR g‐tensor, g_av, of 2.03, at the B3LYP level, is in good agreement with the experimental findings. Heisenberg exchange coupling constants, J, were calculated within the broken‐symmetry formalism, at the B3LYP level, showing that the two iron centers are antiferromagnetic coupling, with a very weak coupling constant of about ?7 cm^?1, in good agreement with the experimental value. Additionally, the effect of using different multiplicities of the reference geometries on the computed J value is discussed. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Databases for transition element bonding: metal-metal bond energies and bond lengths and their use to test hybrid, hybrid meta, and meta density functionals and generalized gradient approximations

We propose a data set of bond lengths for 8 selected transition metal dimers (Ag(2), Cr(2), Cu(2), CuAg, Mo(2), Ni(2), V(2), and Zr(2)) and another data set containing their atomization energies and the atomization energy of ZrV, and we use these for testing density functional theory. The molecules chosen for the test sets were selected on the basis of the expected reliability of the data and their ability to constitute a diverse and representative set of transition element bond types while the data sets are kept small enough to allow for efficient testing of a large number of computational methods against a very reliable subset of experimental data. In this paper we test 42 different functionals: 2 local spin density approximation (LSDA) functionals, 12 generalized gradient approximation (GGA) methods, 13 hybrid GGAs, 7 meta GGA methods, and 8 hybrid meta GGAs. We find that GGA density functionals are more accurate for the atomization energies of pure transition metal systems than are their meta, hybrid, or hybrid meta analogues. We find that the errors for atomization energies and bond lengths are not as large if we limit ourselves to dimers with small amounts of multireference character. We also demonstrate the effects of increasing the fraction of Hartree-Fock exchange in multireference systems by computing the potential energy curve for Cr(2) and Mo(2) with several functionals. We also find that BLYP is the most accurate functional for bond energies and is reasonably accurate for bond lengths. The methods that work well for transition metal bonds are found to be quite different from those that work well for organic and other main group chemistry.     

Performance of density functionals for calculation of reductive ring-opening reaction energies of Li<Superscript>+</Superscript>-EC and Li<Superscript>+</Superscript>-VC

Reaction energies were determined for reductive ring-opening reactions of Li⁺-coordinated ethylene carbonate (EC) and vinylene carbonate (VC) by using various density functional theory (DFT) and ab-initio methods applying the basis sets up to Dunnings aug-cc-pVQZ. The methods examined include the local density functional (SVWN), the pure gradient-corrected density functionals (BLYP and BPW91), and the hybrid density functionals (B3LYP, B1LYP, B3PW91, and mPW1PW91). Comparison of the DFT results with ab-initio results indicates that the mPW1PW91 approach introduced by Adamo and Barone, is superior to all the other DFT methods (including B3LYP). The performance of more cost-effective Pople-type basis sets ranging from 6-31G(d,p) to 6-311++G(3df,3pd) was assessed at DFT levels of theory by calibrating them with the aug-cc-pVQZ results     

Parameters for excess electron transfer in DNA. Estimation using unoccupied Kohn-Sham orbitals and TD DFT

We show that the energetics and electronic couplings for excess electron transfer (EET) can be accurately estimated by using unoccupied Kohn-Sham orbitals (UKSO) calculated for neutral pi stacks. To assess the performance of different DFT functionals, we use MS-PT2 results for seven pi stacks of nucleobases as reference data. The DFT calculations are carried out by using the local spin density approximation SVWN, two generalized gradient approximation functionals BP86 and BLYP, and two hybrid functionals B3LYP and BH&HLYP. Best estimations within the UKSO approach are obtained by the B3LYP and SVWN methods. TD DFT calculations provide less accurate values of the EET parameters as compared with the UKSO data. Also, the excess charge distribution in the radical anions is well described by the LUMOs of neutral systems. In contrast, spin-unrestricted DFT calculations of radical anions considerably overestimate delocalization of the excess electron. The excellent results obtained for the ground and excited states of the radical anions (excitation energy, transition dipole moment, electronic coupling, and excess electron distribution) by using UKSO of neutral dimers suggest an efficient strategy to calculate the EET parameters for DNA pi stacks.     

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.     

Scrutinizing “Invisible” astatine: A challenge for modern density functionals

The main‐group 6p elements did not receive much attention in the development of recent density functionals. In many cases it is still difficult to choose among the modern ones a relevant functional for various applications. Here, we illustrate the case of astatine species (At, Z = 85) and we report the first, and quite complete, benchmark study on several properties concerning such species. Insights on geometries, transition energies and thermodynamic properties of a set of 19 astatine species, for which reference experimental or theoretical data has been reported, are obtained with relativistic (two‐component) density functional theory calculations. An extensive set of widely used functionals is employed. The hybrid meta‐generalized gradient approximation (meta‐GGA) PW6B95 functional is overall the best choice. It is worth noting that the range‐separated HSE06 functional as well as the old and very popular B3LYP and PBE0 hybrid‐GGAs appear to perform quite well too. Moreover, we found that astatine chemistry in solution can accurately be predicted using implicit solvent models, provided that specific parameters are used to build At cavities. © 2016 Wiley Periodicals, Inc.     

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.     

On the DFT Ground State of Crystalline Bromine and Iodine

We report on an erroneous ground state within common density functional theory (DFT) methods for the solid elements bromine and iodine. Phonon computations at the GGA level for both molecular crystals yield imaginary vibrational modes, erroneously indicating dynamic instability—that fact alone could easily pass as a computational artefact, but these imaginary modes lead to energetically more favorable and dynamically stable structures, made up of infinite monoatomic chains. In contrast, meta‐GGA and hybrid functionals yield the correct energetic order for bromine, while for iodine, most global hybrids do not improve the GGA result significantly. The qualitatively correct answer, in both cases, is given by the long‐range corrected hybrid LC‐ωPBE, the Minnesota functionals M06L and M06, and by periodic Hartree–Fock and MP2 theory. This poor performance of economic DFT functionals should be kept in mind, for example, during global structure optimizations of systems with significant contributions from halogen bonds.     

Accurate energies calculated by empirical corrections to the local spin density approximation

We tested the efficacy of three empirical correction schemes on atomization energies calculated by the following theories: Kohn–Sham density functional theory (KS-DFT) with local spin density approximation (LSDA), two KS-DFT gradient-corrected methods, one hybrid Hartree–Fock/KS-DFT method similar to B3LYP, and the ab initio extrapolation procedures G1 and G2. Empirical corrections improved the LSDA results greatly, while the other theories were improved slightly or not at all. The best procedure for correcting LSDA atomization energies brings the mean absolute deviation from experiment from 38.3 to 4.0 kcal/mol on a subset of 44 molecules in the G1 dataset that were not used in deriving the empirical parameters. This corrected LSDA is interesting for three reasons: it could be a useful computational tool in some cases, it implies that the LSDA itself gives accurate energies for reactions where atomic coordinations stay unchanged, and it gives insight into the search of better functionals.