Rung 3.5 density functionals: Another step on Jacob's ladder

Applications of density functional theory (DFT) to computational chemistry and solid‐state physics rely on a “Jacob's Ladder” of progressively more complicated approximations to the many‐body exchange‐correlation (XC) density functional. Accurate, computationally tractable DFT calculations on large and periodic systems remain challenging for existing XC functionals. Simple XC functionals on the three lowest rungs of Jacob's Ladder are insufficiently accurate for many properties, while fourth‐rung hybrid functionals incorporating nonlocal information can be prohibitively expensive. This perspective presents our work toward a compromise, a new class of “Rung 3.5” functionals that incorporate a linear dependence on the nonlocal one‐particle density matrix. This work reviews these functionals' formal underpinning, numerical performance, and prospects for modeling solids and surfaces. © 2012 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.     

About the difference between density functionals defined by energy criterion and density functionals defined by density criterion: Exchange functionals

The difference between density functionals defined by energy criterion and density functionals defined by density criterion is studied for the exchange functional. It is shown that Slater potentials are exact exchange potentials in the sense that they yield the Hartree–Fock electron density if all operators are given by local expressions. © 2016 Wiley Periodicals, Inc.     

Performance of range‐separated hybrid exchange–correlation functionals for the calculation of magnetic exchange coupling constants of organic diradicals

The prediction of magnetic behavior is important for the design of new magnetic materials. Kohn–Sham density functional theory is popular for this purpose, although one should be careful about choosing the right exchange–correlation functional. Here, we perform a statistical analysis to test different range‐separated hybrid density functionals for the calculation of magnetic exchange coupling constants J of fourteen organic diradicals. Our analysis suggests that in absolute terms the MN12SX functional performs best among the series of twelve functionals studied here (including the popular B3LYP), followed by N12SX functionals along with Scuseria's HSE series of functionals. LC‐ PBE was found to be the least accurate, which is in contrast with its good performance for calculating J for transition metal complexes. The HSE family of functionals and B3LYP are the only functionals to reproduce the qualitative trends of the coupling constants correctly for the ferromagnetically coupled diradicals under study. © 2017 Wiley Periodicals, Inc.     

Description of core excitations by time-dependent density functional theory with local density approximation, generalized gradient approximation, meta-generalized gradient approximation, and hybrid functionals

Time-dependent density functional theory (TDDFT) is employed to investigate exchange-correlation-functional dependence of the vertical core-excitation energies of several molecules including H, C, N, O, and F atoms. For the local density approximation (LDA), generalized gradient approximation (GGA), and meta-GGA, the calculated X1s-->pi* excitation energies (X = C, N, O, and F) are severely underestimated by more than 13 eV. On the other hand, time-dependent Hartree-Fock (TDHF) overestimates the excitation energies by more than 6 eV. The hybrid functionals perform better than pure TDDFT because HF exchange remedies the underestimation of pure TDDFT. Among these hybrid functionals, the Becke-Half-and-Half-Lee-Yang-Parr (BHHLYP) functional including 50% HF exchange provides the smallest error for core excitations. We have also discovered the systematic trend that the deviations of TDHF and TDDFT with the LDA, GGA, and meta-GGA functionals show a strong atom-dependence. Namely, their deviations become larger for heavier atoms, while the hybrid functionals are significantly less atom-dependent.     

libKEDF: An accelerated library of kinetic energy density functionals

Kinetic energy density functionals (KEDFs) approximate the kinetic energy of a system of electrons directly from its electron density. They are used in electronic structure methods that lack direct access to orbitals, for example, orbital‐free density functional theory (OFDFT) and certain embedding schemes. In this contribution, we introduce libKEDF, an accelerated library of modern KEDF implementations that emphasizes nonlocal KEDFs. We discuss implementation details and assess the performance of the KEDF implementations for large numbers of atoms. We show that using libKEDF, a single computing node or (GPU) accelerator can provide easy computational access to mesoscale chemical and materials science phenomena using OFDFT algorithms. © 2017 Wiley Periodicals, Inc.     

Molecular structure and vibrational spectra of indole and 5-aminoindole by density functional theory and ab initio Hartree–Fock calculations

The molecular geometry and vibrational frequencies of indole and 5-aminoindole in the ground state have been calculated by using the Hartree–Fock and density functional method with 6-311++G(d,p) basis set. The optimized geometrical parameters obtained by DFT calculations are in good agreement with the experimental values. Comparison of the observed fundamental vibrational frequencies of indole and 5-aminoindole with the calculated results by density functional and Hartree–Fock methods indicates that B3LYP is superior to the scaled Hartree–Fock approach for molecular vibrational problems. The theoretical spectrograms for FT-IR spectrum of 5-aminoindole have been constructed.     

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.     

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     

Conclusive evidence on the insensitivity of additive rules to the combinational details of exchange and correlation functional in hybrid DFT methods

This account is a quest to evaluate the sensitivity of additive rules among total electronic energies of some triple molecular sets to the employed computational level. To achieve the goal, seven randomly generated hybrid functionals have been constructed and employed for calculation on a test set of 29 molecules containing second and third row elements. The computational results based on these nonoptimized DFT methods in conjunction with those of standard B3LYP method clearly demonstrates that in contrast to atomization or isomerization energies, additive rules among total electronic energies are insensitive to reliability or precision of computational method. This insensitivity to details of exchange‐correlation functional indicates that additive rules can be generally well considered using simple computational methods. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Highly efficient implementation of pseudospectral time‐dependent density‐functional theory for the calculation of excitation energies of large molecules

We have developed and implemented pseudospectral time‐dependent density‐functional theory (TDDFT) in the quantum mechanics package Jaguar to calculate restricted singlet and restricted triplet, as well as unrestricted excitation energies with either full linear response (FLR) or the Tamm–Dancoff approximation (TDA) with the pseudospectral length scales, pseudospectral atomic corrections, and pseudospectral multigrid strategy included in the implementations to improve the chemical accuracy and to speed the pseudospectral calculations. The calculations based on pseudospectral time‐dependent density‐functional theory with full linear response (PS‐FLR‐TDDFT) and within the Tamm–Dancoff approximation (PS‐TDA‐TDDFT) for G2 set molecules using B3LYP/6‐31G*^* show mean and maximum absolute deviations of 0.0015 eV and 0.0081 eV, 0.0007 eV and 0.0064 eV, 0.0004 eV and 0.0022 eV for restricted singlet excitation energies, restricted triplet excitation energies, and unrestricted excitation energies, respectively; compared with the results calculated from the conventional spectral method. The application of PS‐FLR‐TDDFT to OLED molecules and organic dyes, as well as the comparisons for results calculated from PS‐FLR‐TDDFT and best estimations demonstrate that the accuracy of both PS‐FLR‐TDDFT and PS‐TDA‐TDDFT. Calculations for a set of medium‐sized molecules, including C_n fullerenes and nanotubes, using the B3LYP functional and 6‐31G^** basis set show PS‐TDA‐TDDFT provides 19‐ to 34‐fold speedups for C_n fullerenes with 450–1470 basis functions, 11‐ to 32‐fold speedups for nanotubes with 660–3180 basis functions, and 9‐ to 16‐fold speedups for organic molecules with 540–1340 basis functions compared to fully analytic calculations without sacrificing chemical accuracy. The calculations on a set of larger molecules, including the antibiotic drug Ramoplanin, the 46‐residue crambin protein, fullerenes up to C₅₄₀ and nanotubes up to 14×(6,6), using the B3LYP functional and 6‐31G^** basis set with up to 8100 basis functions show that PS‐FLR‐TDDFT CPU time scales as N^2.05 with the number of basis functions. © 2016 Wiley Periodicals, Inc.     

Toward assessment of density functionals for vibronic coupling in two‐photon absorption: A case study of 4‐nitroaniline

In this study, we predict vibronic two‐photon absorption (TPA) spectra for 4‐nitroaniline in vacuo. The simulations are performed using density functional theory and the approximate second‐order coupled‐cluster singles and doubles model CC2. Thereby we also demonstrate the possibility of simulations of vibronic TPA spectra with ab initio wavefunction methods that include electron correlation for medium‐sized systems. A special focus is put on the geometric derivatives of the second‐order transition moment and the dipole moment difference between the charge‐transfer excited state and the ground state. The results of CC2 calculations bring new insight into the vibronic coupling mechanism in TPA spectra of 4‐nitroniline and demonstrate that the mixed term is quite large and that it also exhibits a negative interference with the Franck‐Condon contribution. © 2015 Wiley Periodicals, Inc.     

Structure and conformational stability of CH2CHCH2X (X=F, Cl and Br) molecules—post Hartree–Fock and density functional theory methods

The molecular structure and conformational stability of CH₂CHCH₂X (X=F, Cl and Br) molecules were studied using ab initio and density functional theory (DFT) methods. The molecular geometries of 3-fluoropropene were optimized employing BLYP and B3LYP levels of theory of DFT method implementing 6-311+G(d,p) basis set. The MP2/6-31G^*, BLYP and B3LYP levels of theory of ab initio and DFT methods were used to optimize the 3-chloropropene and 3-bromopropene molecules. The structural and physical parameters of the molecules are discussed with the available experimental values. The rotational potential energy surface of the above molecules were obtained at MP2/6-31G^* and B3LYP/6-311+G(d,p) levels of theory. The Fourier decomposition of the rotational potentials were analyzed. The HF/6-31G^* and MP2/6-31G^* levels of theory have predicted the cis conformer as the minimum energy structure for 3-fluoropropene, which is in agreement with the experimental values, whereas the BLYP/6-311+G(d,p) and B3LYP/6-311+G(d,p) levels of theory reverses the order of conformation. The ΔE values calculated for 3-chloropropene at MP2/6-31G^*, BLYP/6-311+G(d,p) and B3LYP/6-311+G(d,p) levels of theory show that the gauche form is more stable than the cis form, which is in agreement with the experimental value. The same levels of theory have also predicted that the gauche form is stable than cis for 3-bromopropene molecule. The maximum hardness principle has been able to predict the stable conformer of 3-fluoropropene at HF/6-31G^* level of theory, but the same level of theory reverses the conformational stability of 3-chloropropene and 3-bromopropene molecules and MP2/6-31G^* level of theory predicted the stable conformer correctly.     

Some applications of local density functional theory to the calculation of reaction energetics

Summary We report the results of a local density functional investigation of the energetics of some isomerization reactions, involving the conversions of several unsaturated systems to highly strained molecules related to triprismane and tetrahedrane. The program DMol was used at the DNP level to compute the activation barriers and total energy changes associated with these processes. We also show, for more than 70 first- and second-row atoms and molecules, that the errors (non-local corrections) in their energies correlate very well with the number of electrons, within isonuclear series. This should provide a useful empirical means for improving dissociation energies obtained within the local approximation.     

Performance of DFT hybrid functionals in the theoretical treatment of H‐bonds: Analysis term‐by‐term

The performance of DFT to reproduce noncovalent interactions like H‐bonds was assessed. Three DFT hybrid functionals (B3LYP, BHandHLYP from Jaguar5.5, and BHandHLYP from G03) were used to calculate the interaction energies and H‐bond distances of several host–guest complexes theoretically designed. Two reference systems (whose experimental data were available) were also calculated at the same levels of theory. In all the cases B3LYP and BHandHLYP from G03 gave rather the same results of interaction energy and distances, whereas the functional BHandHLYP from Jaguar5.5 overestimated the interaction energies and underestimated the H‐bond distances. Considering the construction of each functional, the terms responsible for such differences are the gradient correction to the exchange functional (Becke88) and the correlation functional (VWN, LYP) and not the Hartree–Fock contribution. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Link between the kinetic‐ and exchange‐energy functionals in the generalized gradient approximation

An approximate kinetic‐energy functional of the generalized gradient approximation form was derived following the “conjointness conjecture” of Lee, Lee, and Parr. The functional shares the analytical form of its gradient dependency with the exchange‐energy functionals of Becke and Perdew, Burke, and Ernzerhof. The two free parameters of this functional were determined using the exact values of the kinetic energy of He and Xe atoms. A set of 12 closed‐shell atoms was used to test the accuracy of the proposed functional and more than 30 others taken from the literature. It is shown that the conjointness conjecture leads to a very good class of kinetic‐energy functionals. Moreover, the functional developed in this work is shown to be one of the most accurate despite its simple analytical form. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2002     

Attractive electron–electron interactions within robust local fitting approximations

An analysis of Dunlap's robust fitting approach reveals that the resulting two‐electron integral matrix is not manifestly positive semidefinite when local fitting domains or non‐Coulomb fitting metrics are used. We present a highly local approximate method for evaluating four‐center two‐electron integrals based on the resolution‐of‐the‐identity (RI) approximation and apply it to the construction of the Coulomb and exchange contributions to the Fock matrix. In this pair‐atomic resolution‐of‐the‐identity (PARI) approach, atomic‐orbital (AO) products are expanded in auxiliary functions centered on the two atoms associated with each product. Numerical tests indicate that in 1% or less of all Hartree–Fock and Kohn–Sham calculations, the indefinite integral matrix causes nonconvergence in the self‐consistent‐field iterations. In these cases, the two‐electron contribution to the total energy becomes negative, meaning that the electronic interaction is effectively attractive, and the total energy is dramatically lower than that obtained with exact integrals. In the vast majority of our test cases, however, the indefiniteness does not interfere with convergence. The total energy accuracy is comparable to that of the standard Coulomb‐metric RI method. The speed‐up compared with conventional algorithms is similar to the RI method for Coulomb contributions; exchange contributions are accelerated by a factor of up to eight with a triple‐zeta quality basis set. A positive semidefinite integral matrix is recovered within PARI by introducing local auxiliary basis functions spanning the full AO product space, as may be achieved by using Cholesky‐decomposition techniques. Local completion, however, slows down the algorithm to a level comparable with or below conventional calculations. © 2013 Wiley Periodicals, Inc.     

Carbondioxide rare‐gas systems: Sensitivity of basis sets and double‐hybrid density functionals

This study emphasizes on the performance of six newly developed double‐hybrid density functionals (DHDF) in explaining the potential energy curves of different carbondioxide rare‐gas systems. The basis set sensitivity has also been explored with the use of three basis sets. Our results suggest that for lighter He/Ne‐CO₂ complexes, proper choice of DHDF and basis set lead to results those matches exactly with earlier calculations and also with the experiment. On the other hand, for heavier Ar/Kr‐CO₂ complexes although the equilibrium separation distance matches exactly with earlier observations, the interaction energy values lie far apart. The overall investigation emphasizes on the fact that one has to tune the methods and basis sets properly to achieve good and satisfactory results. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Order of stabilities in water nanoclusters: Insight from some recent double‐hybrid functionals

In the present work, the applicability of some of the recently proposed and modern double‐hybrid (DH) models and other density functional theory (DFT) approximations has been analyzed for a difficult test, the order of stability in low‐energy isomers of water nanoclusters. In particular, we aim to systematically investigate for these functionals the role played by several factors such as dispersion correction, integrand functions upon which the DHs are based, and different spin scaling for the perturbative term in DH calculations of the relative energies for various isomers of water nanoclusters (H₂O)₂₀. From the obtained results, the superior performance of DHs with respect to the functionals from previous rungs is confirmed. It is shown that the dispersion corrected DHs perform better than noncorrected counterparts. Plus, the DH models based on cubic integrand (CI) and quadratic integrand (QI) functions are nearly equivalent in performance. We also find that using only contributions of electron pairs with opposite spin for the perturbative correlation part through scaled opposite spin scheme does not represent a significant improvement on accuracy of the results. Putting all the results together, the dispersion corrected parameterized DHs and parameter‐free DH models involving CI and QI functions outperform other approximations for relative energies of water 20‐mers. Altogether, predicting the correct order of the stability in water nanoclusters may be considered as another Achilles' heel in DFT calculations, although more analyses in this context are still needed. © 2016 Wiley Periodicals, Inc.     

Fitting elephants in the density functionals zoo: Statistical criteria for the evaluation of density functional theory methods as a suitable replacement for counting parameters

Counting parameters has become customary in the density functional theory community as a way to infer the transferability of popular approximations to the exchange‐correlation functionals. Recent work in data science, however, has demonstrated that the number of parameters of a fitted model is not related to the complexity of the model itself, nor to its eventual overfitting. Using similar arguments, here, we show that it is possible to represent every modern exchange‐correlation functional approximations using just one single parameter. This procedure proves the futility of the number of parameters as a measure of transferability. To counteract this shortcoming, we introduce and analyze the performance of three statistical criteria for the evaluation of the transferability of exchange‐correlation functionals. The three criteria are called Akaike information criterion, Vapnik‐Chervonenkis criterion, and cross‐validation criterion and are used in a preliminary assessment to rank 60 exchange‐correlation functional approximations using the ASCDB database of chemical data.