Assessing the performance of density functional theory for the dynamic polarizabilities of amino acids: Treatment of correlation and role of exact exchange

The theoretical determination of electric response properties of the biological systems is a field where the application of density functional theory (DFT) appears to be quite promising. In this work, the performance of 41 density functional methods is evaluated in predicting dynamic polarizabilities of an experimental benchmark set of 20 proteinogenic amino acids. The behavior of a large number of density functionals, including various types of the local spin density approximation (LSDA), generalized gradient approximation (GGA), meta‐GGA (m‐GGA), hybrid‐GGA (h‐GGA), hybrid meta‐GGA (hm‐GGA), and range‐separated hybrid‐GGA (rsh‐GGA), has been assessed for the purpose. Analyzing the results of our DFT benchmarking, we found that these computationally economical methods show very diverse predictive capability and a careful selection of DFT functionals is very important in the polarizability calculations. Considering the role of exchange, correlation, dispersion and long‐range corrections, it turned out that in the LSDA class, SVWN3 gives better results than SPL and SVWN5 toward the reference values. Of the GGA methods, OPBE outperforms all other functionals. The M06‐L is the best method of m‐GGA class. The B3LYP and TPSSh are the best functionals of h‐GGA and hm‐GGA lineages, respectively. Finally, CAM‐B3LYP is the best method of rsh‐GGA functionals that predicts the most accurate polarizability for amino acids by a large margin with respect to others. Overall, the best performing functionals turn out to be hm‐GGAs TPSSh, TPSS1KCIS, M05, tau‐HCTHhyb, and h‐GGA B3LYP. Hopefully, the results of this investigation might provide the useful guidance to propose a new exchange‐correlation functional for calculating the optical properties of biomolecular materials. © 2013 Wiley Periodicals, Inc.     

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     

Vibrational frequency scale factors for density functional theory and the polarization consistent basis sets

Calculated harmonic vibrational frequencies systematically deviate from experimental vibrational frequencies. The observed deviation can be corrected by applying a scale factor. Scale factors for: (i) harmonic vibrational frequencies [categorized into low (<1000 cm^?1) and high (>1000 cm^?1)], (ii) vibrational contributions to enthalpy and entropy, and (iii) zero‐point vibrational energies (ZPVEs) have been determined for widely used density functionals in combination with polarization consistent basis sets (pc‐n, n = 0,1,2,3,4). The density functionals include pure functionals (BP86, BPW91, BLYP, HCTH93, PBEPBE), hybrid functionals with Hartree‐Fock exchange (B3LYP, B3P86, B3PW91, PBE1PBE, mPW1K, BH&HLYP), hybrid meta functionals with the kinetic energy density gradient (M05, M06, M05‐2X, M06‐2X), a double hybrid functional with Møller‐Plesset correlation (B2GP‐PLYP), and a dispersion corrected functional (B97‐D). The experimental frequencies for calibration were from 41 organic molecules and the ZPVEs for comparison were from 24 small molecules (diatomics, triatomics). For this family of basis sets, the scale factors for each property are more dependent on the functional selection than on basis set level, and thus allow for a suggested scale factor for each density functional when employing polarization consistent basis sets (pc‐n, n = 1,2,3,4). A separate scale factor is recommended when the un‐polarized basis set, pc‐0, is used in combination with the density functionals. © 2012 Wiley Periodicals, Inc.     

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.     

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.     

Evaluation of exchange‐correlation functionals in comparison to B3LYP for the description of silicon and Cu‐doped silicon clusters

Several developed exchange‐correlation functionals in density functional theory have been systematically applied to describe the geometries and electronic properties of small silicon (Si_n+1, n < 5) and doped silicon (CuSi_n) clusters. The performance of the various approaches is done with their critical comparison with B3LYP and available high level wave function methods. Our calculations indicate that all functional give reasonable results. Further, OLYP/6‐311+G* approach generally agrees with B3LYP results. The good performance of OLYP is of significant interest knowing that the hybrid functionals are computationally more demanding than nonhybrid schemes. So, we recommend OLYP/6‐311+G* approach for studying the doped silicon clusters and understanding the electronic properties of silicon by the presence of doped metal impurities. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Standard grids for high‐precision integration of modern density functionals: SG‐2 and SG‐3

Density‐functional approximations developed in the past decade necessitate the use of quadrature grids that are far more dense than those required to integrate older generations of functionals. This category of difficult‐to‐integrate functionals includes meta‐generalized gradient approximations, which depend on orbital gradients and/or the Laplacian of the density, as well as functionals based on B97 and the popular “Minnesota” class of functionals, each of which contain complicated and/or oscillatory expressions for the exchange inhomogeneity factor. Following a strategy introduced previously by Gill and co‐workers to develop the relatively sparse “SG‐0” and “SG‐1” standard quadrature grids, we introduce two higher‐quality grids that we designate SG‐2 and SG‐3, obtained by systematically “pruning” medium‐ and high‐quality atom‐centered grids. The pruning procedure affords computational speedups approaching a factor of two for hybrid functionals applied to systems of atoms, without significant loss of accuracy. The grid dependence of several popular density functionals is characterized for various properties. © 2017 Wiley Periodicals, Inc.     

Spin‐component‐scaled double hybrids: An extensive search for the best fifth‐rung functionals blending DFT and perturbation theory

Following up on an earlier preliminary communication (Kozuch and Martin, Phys. Chem. Chem. Phys. 2011, 13 , 20104), we report here in detail on an extensive search for the most accurate spin‐component‐scaled double hybrid functionals [of which conventional double hybrids (DHs) are a special case]. Such fifth‐rung functionals approach the performance of composite ab initio methods such as G3 theory at a fraction of their computational cost, and with analytical derivatives available. In this article, we provide a critical analysis of the variables and components that maximize the accuracy of DHs. These include the selection of the exchange and correlation functionals, the coefficients of each component [density functional theory (DFT), exact exchange, and perturbative correlation in both the same spin and opposite spin terms], and the addition of an ad‐hoc dispersion correction; we have termed these parametrizations “DSD‐DFT” (Dispersion corrected, Spin‐component scaled, Double‐hybrid DFT). Somewhat surprisingly, the quality of DSD‐DFT is only mildly dependent on the underlying DFT exchange and correlation components, with even DSD‐LDA yielding respectable performance. Simple, nonempirical GGAs appear to work best, whereas meta‐GGAs offer no advantage (with the notable exception of B95c). The best correlation components appear to be, in that order, B95c, P86, and PBEc, while essentially any good GGA exchange yields nearly identical results. On further validation with a wider variety of thermochemical, weak interaction, kinetic, and spectroscopic benchmarks, we find that the best functionals are, roughly in that order, DSD‐PBEhB95, DSD‐PBEP86, DSD‐PBEPW91, and DSD‐PBEPBE. In addition, DSD‐PBEP86 and DSD‐PBEPBE can be used without source code modifications in a wider variety of electronic structure codes. Sample job decks for several commonly used such codes are supplied as electronic Supporting Information. Copyright © 2013 Wiley Periodicals, Inc.     

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     

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.     

Implementation of nuclear gradients of range‐separated hybrid density functionals and benchmarking on rotational constants for organic molecules

We have implemented the nuclear gradient for several range‐separated hybrid density functionals in the general quantum chemistry code ORCA. To benchmark the performance, we have used a recently published set of back‐corrected gas phase rotational constants, which we extended by three molecules. In our evaluation, CAM‐B3LYP‐D3 and ωB97X‐D3 show great accuracy, and are surpassed by B2PLYP‐D3 only. Lower‐cost alternatives to quadruple‐ζ basis set‐based calculations, among them a smaller basis set and the use of resolution‐of‐the‐identity approaches, are assessed and shown to yield acceptable deviations. In addition, the Hartree‐Fock‐based back‐correction method is compared to a density functional theory alternative, which largely shows consistency between the two. A new, well‐performing, spin‐component scaled MP2 variant is designed and discussed, as well. © 2014 Wiley Periodicals, Inc.     

Analytic derivatives for the XYG3 type of doubly hybrid density functionals: Theory,implementation, and assessment

We present a theoretical development of the equations required to perform an analytic geometry optimization of a molecular system using the XYG3 type of doubly hybrid (xDH) functionals. In contrast to the well‐established B2PLYP type of DH functionals, the energy expressions in the xDH functionals are constructed by using density and orbital information from another standard Kohn–Sham (KS) functional (e.g., B3LYP) for doing the self‐consistent field calculations. Thus, the xDH functionals are nonvariational in both the hybrid density functional part and the second‐order perturbation part, each of which requires formally to solve a coupled‐perturbed KS equation. An implementation is reported here which combines the two parts by defining a total Lagrangian such that only a single set of the Z‐vector equations need to be solved. The computational cost with our implementation is of the same order as those for the conventional Møller–Plesset theory to the second order (MP2) and B2PLYP. Systematic test calculations are provided for covalently bonded molecules as well as compounds involving the intramolecular nonbonded interactions for the main group elements. Satisfactory performance of the xDH functionals demonstrates that the extra computer time on top of the conventional KS procedure is well‐invested, in particular, when the standard KS functionals and MP2 as well, are problematic. © 2013 Wiley Periodicals, Inc.     

Prediction of the reduction potential in transition‐metal containing complexes: How expensive? For what accuracy?

Accurate computationally derived reduction potentials are important for catalyst design. In this contribution, relatively inexpensive density functional theory methods are evaluated for computing reduction potentials of a wide variety of organic, inorganic, and organometallic complexes. Astonishingly, SCRF single points on B3LYP optimized geometries with a reasonably small basis set/ECP combination works quite well‐‐B3LYP with the BS1 [modified‐LANL2DZ basis set/ECP (effective core potential) for metals, LANL2DZ(d,p) basis set/LANL2DZ ECP for heavy nonmetals (Si, P, S, Cl, and Br), and 6‐31G(d') for other elements (H, C, N, O, and F)] and implicit PCM solvation models, SMD (solvation model based on density) or IEFPCM (integral equation formalism polarizable continuum model with Bondi atomic radii and α = 1.1 reaction field correction factor). The IEFPCM‐Bondi‐B3LYP/BS1 methodology was found to be one of the least expensive and most accurate protocols, among six different density functionals tested (BP86, PBEPBE, B3LYP, B3P86, PBE0, and M06) with thirteen different basis sets (Pople split‐valence basis sets, correlation consistent basis sets, or Los Alamos National Laboratory ECP/basis sets) and four solvation models (SMD, IEFPCM, IPCM, and CPCM). The MAD (mean absolute deviation) values of SCRF‐B3LYP/BS1 of 49 studied species were 0.263 V for SMD and 0.233 V for IEFPCM‐Bondi; and the linear correlations had respectable R ² values (R ² = 0.94 for SMD and R ² = 0.93 for IEFPCM‐Bondi). These methodologies demonstrate relatively reliable, convenient, and time‐saving functional/basis set/solvation model combinations in computing the reduction potentials of transition metal complexes with moderate accuracy. © 2017 Wiley Periodicals, Inc.     

Electronic structure and spectroscopic properties of 6‐aminophenanthridine and its derivatives: Insights from density functional theory

6‐Aminophenanthridine (6AP) and its derivatives show important biological activities as antiprion compounds and inhibitors of the protein folding activity of the ribosome. Both of these activities depend on the RNA binding property of these compounds, which has been recently characterized by fluorescence spectroscopy. Hence, fundamental insights into the photophysical properties of 6AP compounds are highly important to understand their biological activities. In this work, we have calculated electronic structures and optical properties of 6AP and its three derivatives 6AP8CF₃, 6AP8Cl, and 6APi by density functional theory (DFT) and time‐dependent density functional theory (TDDFT). Our calculated spectra show a good agreement with the experimental absorption and fluorescence spectra, and thus, provide deep insights into the optical properties of the compounds. Furthermore, comparing the results obtained with four different hybrid functionals, we demonstrate that the accuracy of the functionals varies in the order B3LYP > PBE0 > M062X > M06HF. © 2015 Wiley Periodicals, Inc.     

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.     

Computational Study of Cage Like (ZnO)12 Cluster Using Hybrid and Hybrid Meta Functionals

Density Functional Theory employing hybrid and M06 functionals in combination with three different basis sets is used to calculate the ground state of a cage like (ZnO)₁₂ nanocluster which has been consistently reported as the more stable cluster for its particular size. B3LYP and B3PW91 hybrid functionals combined with 6‐31+G*, Lanl2dz and SDD basis sets are employed to treat the ZnO molecular system. Alternatively, three M06 functionals in combination with three basis sets are employed in the nanostructure calculations. Results obtained by treating ZnO sodalite cage nanocluster with M06 functionals demonstrated comparable quality to results obtained with hybrid functionals. Within this study, efficient theoretical DFT methods with the widely known hybrid and the recently created M06 meta‐hybrid functionals are employed to study nanostructured ZnO. Our resulting parameters provide a fresh approach performance wise on the different theoretical methods to treat transition metal nanostructures, particularly, ZnO nanoclusters geometry and electronic structure.