Magnetic exchange couplings evaluated with Rung 3.5 density functionals

Rung 3.5 exchange-correlation functionals are assessed for the calculation of magnetic exchange coupling parameters and atomic spin populations for a variety of inorganic and organic magnetic systems. Density functional theory calculations of exchange couplings sensitively depend on nonlocal contributions to the exchange-correlation functional. Semilocal functionals, Rungs 1-3 on "Jacob's Ladder" of density functional approximations, yield excessively delocalized electrons and overestimated absolute exchange couplings. Fourth-rung hybrid functionals admixing nonlocal exchange improve the results. We show that new "Rung 3.5" functionals give magnetic properties intermediate between semilocal and hybrid functionals, providing additional evidence that these functionals incorporate some desirable aspects of nonlocal exchange. Results for ferromagnetic complexes indicate areas for future improvement.     

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.     

Comparing modern density functionals for conjugated polymer band structures: screened hybrid, Minnesota, and Rung 3.5 approximations

Semiconducting polymers with π-conjugated backbones show promise in fields such as photovoltaics. Practical applications of conjugated polymers require precise control over the polymer's electronic band structure. Several new classes of density functional approximation, including screened hybrids, semilocal Minnesota functionals, and Rung 3.5 functionals, show potential for improved predictions of conjugated polymer band structures. This work compares these methods to standard global hybrid density functionals for bandgaps and band structures of representative conjugated polymers. The new methods exhibit particular promise for modeling three-dimensionally periodic bulk polymers, which can be problematic for global hybrids.     

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.     

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.     

Effective local potentials for orbital-dependent density functionals

Practicality of the Kohn-Sham density functional scheme for orbital-dependent functionals hinges on the availability of an efficient procedure for constructing local exchange-correlation potentials in finite basis sets. We have shown recently that the optimized effective potential (OEP) method, commonly used for this purpose, is not free from difficulties. Here we propose a robust alternative to OEPs, termed effective local potentials (ELPs), based on minimizing the variance of the difference between a given nonlocal potential and its desired local counterpart. The ELP method is applied to the exact-exchange-only problem and shown to be promising for overcoming troubles with OEPs.     

Development of density functionals for thermochemical kinetics

A density functional theory exchange-correlation functional for the exploration of reaction mechanisms is proposed. This functional, denoted BMK (Boese-Martin for Kinetics), has an accuracy in the 2 kcal/mol range for transition state barriers but, unlike previous attempts at such a functional, this improved accuracy does not come at the expense of equilibrium properties. This makes it a general-purpose functional whose domain of applicability has been extended to transition states, rather than a specialized functional for kinetics. The improvement in BMK rests on the inclusion of the kinetic energy density together with a large value of the exact exchange mixing coefficient. For this functional, the kinetic energy density appears to correct "back" the excess exact exchange mixing for ground-state properties, possibly simulating variable exchange.     

Inequalities for fermion density matrices

A partial trace over the occupation numbers of all but k states in the density matrix of an ensemble with an arbitrary number of single-particle states is defined as the (reduced) k-state density matrix. This matrix is used to obtain a complete, practical solution to the problem of determining the representability of the diagonal elements of the one- and two-particle (reduced) density matrices. This solution is expressed as a series of linear inequalities involving the density-matrix elements; the inequalities are identical with those derived previously by Davidson and McCrae by a different method. In addition, our method is used to obtain nonlinear, matrix inequalities on the off-diagonal elements of the density matrices.     

General performance of density functionals

The density functional theory (DFT) foundations date from the 1920s with the work of Thomas and Fermi, but it was after the work of Hohenberg, Kohn, and Sham in the 1960s, and particularly with the appearance of the B3LYP functional in the early 1990s, that the widespread application of DFT has become a reality. DFT is less computationally demanding than other computational methods with a similar accuracy, being able to include electron correlation in the calculations at a fraction of time of post-Hartree-Fock methodologies. In this review we provide a brief outline of the density functional theory and of the historic development of the field, focusing later on the several types of density functionals currently available, and finishing with a detailed analysis of the performance of DFT across a wide range of chemical properties and system types, reviewed from the most recent benchmarking studies, which encompass several well-established density functionals together with the most recent efforts in the field. Globally, an overall picture of the level of performance of the plethora of currently available density functionals for each chemical property is drawn, with particular attention being dedicated to the relative performance of the popular B3LYP density functional.     

Applications of some bounds to density functionals for atoms

In this article, some applications of lower bounds to density-dependent functionals in terms of radial expectation values are performed. These diverse applications include accurate upper bounds to the exact kinetic energy and to conjecture a wide set of relationships among radial and momentum expectation values. Some open questions for the improvement of these results are remarked upon. © 1995 John Wiley & Sons, Inc.     

Scaled density functional theory correlation functionals

We show that a simple one-parameter scaling of the dynamical correlation energy estimated by the density functional theory (DFT) correlation functionals helps increase the overall accuracy for several local and nonlocal functionals. The approach taken here has been described as the "scaled dynamical correlation" (SDC) method [Ramachandran, J. Phys. Chem. A 2006, 110, 396], and its justification is the same as that of the scaled external correlation (SEC) method of Brown and Truhlar. We examine five local and five nonlocal (hybrid) DFT functionals, the latter group including three functionals developed specifically for kinetics by the Truhlar group. The optimum scale factors are obtained by use of a set of 98 data values consisting of molecules, ions, and transition states. The optimum scale factors, found with a linear regression relationship, are found to differ from unity with a high degree of correlation in nearly every case, indicating that the deviation of calculated results from the experimental values are systematic and proportional to the dynamic correlation energy. As a consequence, the SDC scaling of dynamical correlation decreases the mean errors (signed and unsigned) by significant amounts in an overwhelming majority of cases. These results indicate that there are gains to be realized from further parametrization of several popular exchange-correlation functionals.     

Combination of pseudopotentials and density functionals

Semilocal pseudopotentials have been determined for first–row (Li to Ne), second row (Na to Ar), and third-row atoms (K, Ca). Core–valence correlation is included by adjusting the pseudopotentials to experimental energies of ions with a single valence electron. Correlation within the valence shell is taken into account by using the spin–density functional formalism. The approximations involved in this approach are tested for atomic ionization energies as well as binding energies of monohydrides and alkali diatomics, agreement with experiment is usually satisfactory, but in certain applications density functionals should be already included in the fitting of the local part of the pseudopotential. In addition, 3s/3p and 3s/2p basis sets (for first and second row, respectively), designed for use in connection with our pseudopotentials, are given; it is shown that they yield reasonable results for both SCF and correlation energies.     

Approximately N-representable density functional density matrices

For the first time, we obtain practical density matrices approximately N-representable by correlated-determinant wave functions, which are functionals of the electron density and entirely defined by information obtainable from the X-ray coherent diffraction experiment. © 1994 John Wiley & Sons, Inc.     

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.     

The importance of middle-range Hartree-Fock-type exchange for hybrid density functionals

Hybrid functionals are responsible for much of the utility of modern Kohn-Sham density functional theory. When rigorously applied to solid-state metallic and small band gap systems, however, the slow decay of their nonlocal Hartree-Fock-type exchange makes hybrids computationally challenging and introduces unphysical effects. This can be remedied by using a range-separated hybrid which only keeps short-range nonlocal exchange, as in the functional of Heyd et al. [J. Chem. Phys. 118, 8207 (2003)]. On the other hand, many molecular properties require full long-range nonlocal exchange, which can also be included by means of a range-separated hybrid such as the recently introduced LC-omegaPBE functional [O. A. Vydrov and G. E. Scuseria, J. Chem. Phys. 125, 234109 (2006)]. In this paper, we show that a three-range hybrid which mainly includes middle-range Hartree-Fock-type exchange and neglects long- and short-range Hartree-Fock-type exchange yields excellent accuracy for thermochemistry, barrier heights, and band gaps, emphasizing that the middle-range part of the 1/r potential seems crucial to accurately model these properties.     

Chlorine--benzene complexes--the reliability of density functionals for non-covalent radical complexes

The structure of the chlorine atom-benzene complex has been a topic of significant controversy for more than 50 years. We have reexamined the structure of this complex with new density functional methods especially designed for non-covalent complexes, and compared the structures and energetics to those obtained using standard DFT and high accuracy composite methods. We find that the popular B3LYP functional fails to identify stationary points revealed by other functionals, and that the eta(1)-sigma complex appears to be more stable than the eta(1)-pi complex, contrary to other recent work, highlighting the careful selection of methods required in non-covalent radical systems.     

Nonuniversality of commonly used correlation-energy density functionals

The correlation energies of the helium isoelectronic sequence and of Hooke's atom isoelectronic sequence have been evaluated using an assortment of local, gradient, and metagradient density functionals. The results are compared with the exact correlation energies, showing that while several of the more recent density functionals reproduce the exact correlation energies of the helium isoelectronic sequence rather closely, none is satisfactory for Hooke's atom isoelectronic sequence. It is argued that the uniformly acceptable results for the helium sequence can be explained through simple scaling arguments that do not hold for Hooke's atom sequence, so that the latter system provides a more sensitive testing ground for approximate density functionals. This state of affairs calls for further effort towards formulating correlation-energy density functionals that would be truly universal at least for spherically symmetric two-fermion systems.