Basis Set Requirements of σ-Functionals for Gaussian- and Slater-Type Basis Functions and Comparison with Range-Separated Hybrid and Double Hybrid Functionals.

σ-Functionals belong to the class of Kohn-Sham (KS) correlation functionals based on the adiabatic-connection fluctuation-dissipation theorem and are technically closely related to the random phase approximation (RPA). They have the same computational demand as the latter, with the computational effort of an energy evaluation for both methods being lower than that of a preceding hybrid DFT calculation for typical systems but yield much higher accuracy, reaching chemical accuracy of 1 kcal/mol for quantities such as reactions and transition energies in main group chemistry. In previous work on σ-functionals, rather large Gaussian basis sets have been used. Here, we investigate the actual basis set requirements of σ-functionals and present three setups that employ smaller Gaussian basis sets ranging from quadruple-ζ (QZ) to triple-ζ (TZ) quality and represent a good compromise between accuracy and computational efficiency. Furthermore, we introduce an implementation of σ-functionals based on Slater-type basis sets and present two setups of QZ and TZ quality for this implementation. We test the accuracy of these setups on a large database of various physical properties and types of reactions, as well as equilibrium geometries and vibrational frequencies. As expected, the accuracy of σ-functional calculations becomes somewhat lower with a decreasing basis set size. However, for all setups considered here, calculations with σ-functionals are clearly more accurate than those within the RPA and even more so than those of the conventional KS methods. For the smallest setup using Gaussian-type basis functions and Slater-type basis functions, we introduce a reparametrization that reduces the loss in accuracy due to the basis set error to some extent. A comparison with the range-separated hybrid ωB97X-V and the double hybrid DSD-BLYP-D3 shows that σ functionals outperform in accuracy both of these accurate and, for their class, representative functionals.

Avoiding spin contamination and spatial symmetry breaking by exact-exchange-only optimized-effective-potential methods within the symmetrized Kohn-Sham framework.

For open-shell atoms and molecules, Kohn-Sham (KS) methods typically resort to spin-polarized approaches that exhibit spin-contamination and often break spatial symmetries. As a result, the KS Hamiltonian operator and the KS orbitals do not exhibit the space and spin symmetry of the physical electron system. The KS formalism can be symmetrized in a rigorous way only in real space, only in spin space, or both in real and spin space. Within such symmetrized KS frameworks, we present exact-exchange-only optimized-effective-potential (OEP) methods that are free of spin contamination and/or spatial symmetry breaking. The effect of symmetrizations on the total energy and its parts and on the exchange potential is analyzed. The presented exact-exchange-only OEP methods may serve as a starting point for high-level symmetrized KS methods based, e.g., on the adiabatic-connection fluctuation-dissipation theorem.

Scaled σ-functionals for the Kohn-Sham correlation energy with scaling functions from the homogeneous electron gas.

The recently introduced σ-functionals constitute a new type of functionals for the Kohn-Sham (KS) correlation energy. σ-Functionals are based on the adiabatic-connection fluctuation-dissipation theorem, are computationally closely related to the well-known direct random phase approximation (dRPA), and are formally rooted in many-body perturbation theory along the adiabatic connection. In σ-functionals, the function of the eigenvalues σ of the Kohn-Sham response matrix that enters the coupling constant and frequency integration in the dRPA is replaced by another function optimized with the help of reference sets of atomization, reaction, transition state, and non-covalent interaction energies. σ-Functionals are highly accurate and yield chemical accuracy of 1 kcal/mol in reaction or transition state energies, in main group chemistry. A shortcoming of σ-functionals is their inability to accurately describe processes involving a change of the electron number, such as ionizations or electron attachments. This problem is attributed to unphysical self-interactions caused by the neglect of the exchange kernel in the dRPA and σ-functionals. Here, we tackle this problem by introducing a frequency- and σ-dependent scaling of the eigenvalues σ of the KS response function that models the effect of the exchange kernel. The scaling factors are determined with the help of the homogeneous electron gas. The resulting scaled σ-functionals retain the accuracy of their unscaled parent functionals but in addition yield very accurate ionization potentials and electron affinities. Moreover, atomization and total energies are found to be exceptionally accurate. Scaled σ-functionals are computationally highly efficient like their unscaled counterparts.

Numerically stable inversion approach to construct Kohn-Sham potentials for given electron densities within a Gaussian basis set framework.

We present a Kohn-Sham (KS) inversion approach to construct KS exchange-correlation potentials corresponding to given electron densities. This method is based on an iterative procedure using linear response to update potentials. All involved quantities, i.e., orbitals, potentials, and response functions, are represented by Gaussian basis functions. In contrast to previous KS inversion methods relying on Gaussian basis sets, the method presented here is numerically stable even for standard basis sets from basis set libraries due to a preprocessing of the auxiliary basis used to represent an exchange-correlation charge density that generates the exchange-correlation potential. The new KS inversion method is applied to reference densities of various atoms and molecules obtained by full configuration interaction or CCSD(T) (coupled cluster singles doubles perturbative triples). The considered examples encompass cases known to be difficult, such as stretched hydrogen or lithium hydride molecules or the beryllium isoelectronic series. For the stretched hydrogen molecule, potentials of benchmark quality are obtained by employing large basis sets. For the carbon monoxide molecule, we show that the correlation potential from the random phase approximation (RPA) is in excellent qualitative and quantitative agreement with the correlation potential from the KS inversion of a CCSD(T) reference density. This indicates that RPA correlation potentials, in contrast to those from semi-local density-functionals, resemble the exact correlation potential. Besides providing exchange-correlation potentials for benchmark purposes, the proposed KS inversion method may be used in density-partition-based quantum embedding and in subsystem density-functional methods because it combines numerical stability with computational efficiency.

Chemical accuracy with σ-functionals for the Kohn-Sham correlation energy optimized for different input orbitals and eigenvalues.

Recently, a new type of orbital-dependent functional for the Kohn-Sham (KS) correlation energy, σ-functionals, was introduced. Technically, σ-functionals are closely related to the well-known direct random phase approximation (dRPA). Within the dRPA, a function of the eigenvalues σ of the frequency-dependent KS response function is integrated over purely imaginary frequencies. In σ-functionals, this function is replaced by one that is optimized with respect to reference sets of atomization, reaction, transition state, and non-covalent interaction energies. The previously introduced σ-functional uses input orbitals and eigenvalues from KS calculations with the generalized gradient approximation (GGA) exchange-correlation functional of Perdew, Burke, and Ernzerhof (PBE). Here, σ-functionals using input orbitals and eigenvalues from the meta-GGA TPSS and the hybrid-functionals PBE0 and B3LYP are presented and tested. The number of reference sets taken into account in the optimization of the σ-functionals is larger than in the first PBE based σ-functional and includes sets with 3d-transition metal compounds. Therefore, also a reparameterized PBE based σ-functional is introduced. The σ-functionals based on PBE0 and B3LYP orbitals and eigenvalues reach chemical accuracy for main group chemistry. For the 10 966 reactions from the highly accurate W4-11RE reference set, the B3LYP based σ-functional exhibits a mean average deviation of 1.03 kcal/mol compared to 1.08 kcal/mol for the coupled cluster singles doubles perturbative triples method if the same valence quadruple zeta basis set is used. For 3d-transition metal chemistry, accuracies of about 2 kcal/mol are reached. The computational effort for the post-self-consistent evaluation of the σ-functional is lower than that of a preceding PBE0 or B3LYP calculation for typical systems.

Numerically stable optimized effective potential method with standard Gaussian basis sets.

We present a numerically stable optimized effective potential (OEP) method based on Gaussian basis sets. The key point of the approach is a sequence of preprocessing steps of the auxiliary basis set used to represent exchange or correlation potentials, the Kohn-Sham (KS) response function, and the right-hand side of the OEP equation in conjunction with a representation of exchange or correlation potentials via exchange or correlation charge densities whose electrostatic potentials generate the potentials. Due to the preprocessing, standard Gaussian basis sets from basis set libraries can be used in OEP calculations. As examples, we present numerical stable computational setups based on aux-cc-pwCVXZ basis sets with X = T, Q, 5 for the orbitals and aux-cc-pVDZ/mp2fit and aux-cc-pVTZ/mp2fit auxiliary basis sets and use them to calculate KS exchange potentials with the exact exchange-only KS method for various atoms and molecules. The resulting exchange potentials not only are numerically stable and physically reasonable but also show convergence with increasing quality of the orbital basis sets. The effect of incorporating exact conditions that the KS exchange potential has to obey is discussed. Moreover, it is briefly demonstrated that the presented approach not only works for KS exchange potentials but equally well for correlation potentials within the direct random phase approximation. Besides for OEP methods, the introduced preprocessing of auxiliary basis sets should also be beneficial in procedures to calculate back effective KS potentials from given electron densities.

Toward chemical accuracy at low computational cost: Density-functional theory with σ-functionals for the correlation energy.

We introduce new functionals for the Kohn-Sham correlation energy that are based on the adiabatic-connection fluctuation-dissipation (ACFD) theorem and are named σ-functionals. Like in the well-established direct random phase approximation (dRPA), σ-functionals require as input exclusively eigenvalues σ of the frequency-dependent KS response function. In the new functionals, functions of σ replace the σ-dependent dRPA expression in the coupling-constant and frequency integrations contained in the ACFD theorem. We optimize σ-functionals with the help of reference sets for atomization, reaction, transition state, and non-covalent interaction energies. The optimized functionals are to be used in a post-self-consistent way using orbitals and eigenvalues from conventional Kohn-Sham calculations employing the exchange-correlation functional of Perdew, Burke, and Ernzerhof. The accuracy of the presented approach is much higher than that of dRPA methods and is comparable to that of high-level wave function methods. Reaction and transition state energies from σ-functionals exhibit accuracies close to 1 kcal/mol and thus approach chemical accuracy. For the 10 966 reactions of the W4-11RE reference set, the mean absolute deviation is 1.25 kcal/mol compared to 3.21 kcal/mol in the dRPA case. Non-covalent binding energies are accurate to a few tenths of a kcal/mol. The presented approach is highly efficient, and the post-self-consistent calculation of the total energy requires less computational time than a density-functional calculation with a hybrid functional and thus can be easily carried out routinely. σ-Functionals can be implemented in any existing dRPA code with negligible programming effort.

Lieb-Oxford bound and pair correlation functions for density-functional methods based on the adiabatic-connection fluctuation-dissipation theorem.

Compliance with the Lieb-Oxford bound for the indirect Coulomb energy and for the exchange-correlation energy is investigated for a number of density-functional methods based on the adiabatic-connection fluctuation-dissipation (ACFD) theorem to treat correlation. Furthermore, the correlation contribution to the pair density resulting from these methods is compared with highly accurate reference values for the helium atom and for the hydrogen molecule at several bond distances. For molecules, the Lieb-Oxford bound is obeyed by all considered methods. For the homogeneous electron gas, it is violated by all methods for low electron densities. The simplest considered ACFD method, the direct random phase approximation (dRPA), violates the Lieb-Oxford bound much earlier than more advanced ACFD methods that, in addition to the simple Hartree kernel, take into account the exchange kernel and an approximate correlation kernel in the calculation of the correlation energy. While the dRPA yields quite poor correlation contributions to the pair density, those from more advanced ACFD methods are physically reasonable but still leave room for improvements, particularly in the case of the stretched hydrogen molecule.

Topological Phase Transitions in Zinc-Blende Semimetals Driven Exclusively by Electronic Temperature.

We show that electronic phase transitions in zinc-blende semimetals with quadratic band touching (QBT) at the center of the Brillouin zone, like GaBi, InBi, or HgTe, can occur exclusively due to a change of the electronic temperature without the need to involve structural transformations or electron-phonon coupling. The commonly used Kohn-Sham density-functional methods based on local and semilocal density functionals employing the local density approximation (LDA) or generalized gradient approximations (GGAs), however, are not capable of describing such phenomena because they lack an intrinsic temperature dependence and account for temperature only via the occupation of bands, which essentially leads only to a shift of the Fermi level without changing the shape or topology of bands. Kohn-Sham methods using the exact temperature-dependent exchange potential, not to be confused with the Hartree-Fock exchange potential, on the other hand, describe such phase transitions. A simple modeling of correlation effects can be achieved by screening of the exchange. In the considered zinc-blende compounds the QBT is unstable at low temperatures and a transition to electronic states without QBT takes place. In the case of HgTe and GaBi Weyl points of type I and type II, respectively, emerge during the transitions. This demonstrates that Kohn-Sham methods can describe such topological phase transitions provided they are based on functionals more accurate than those within the LDA or GGA. Moreover, the electronic temperature is identified as a handle to tune topological materials.

Understanding band gaps of solids in generalized Kohn-Sham theory.

The fundamental energy gap of a periodic solid distinguishes insulators from metals and characterizes low-energy single-electron excitations. However, the gap in the band structure of the exact multiplicative Kohn-Sham (KS) potential substantially underestimates the fundamental gap, a major limitation of KS density-functional theory. Here, we give a simple proof of a theorem: In generalized KS theory (GKS), the band gap of an extended system equals the fundamental gap for the approximate functional if the GKS potential operator is continuous and the density change is delocalized when an electron or hole is added. Our theorem explains how GKS band gaps from metageneralized gradient approximations (meta-GGAs) and hybrid functionals can be more realistic than those from GGAs or even from the exact KS potential. The theorem also follows from earlier work. The band edges in the GKS one-electron spectrum are also related to measurable energies. A linear chain of hydrogen molecules, solid aluminum arsenide, and solid argon provide numerical illustrations.

The nature of one-dimensional carbon: polyynic versus cumulenic.

A question of both fundamental as well as practical importance is the nature of one-dimensional carbon, in particular whether a one-dimensional carbon allotrope is polyynic or cumulenic, that is, whether bond-length alternation occurs or not. By combining the concept of aromaticity and antiaromaticity with the rule of Peierls distortion, the occurrence and magnitude of bond-length alternation in carbon chains with periodic boundary conditions and corresponding carbon rings as a function of the chain or ring length can be explained. The electronic properties of one-dimensional carbon depend crucially on the bond-length alternation. Whereas it is generally accepted that carbon chains in the limit of infinite length have a polyynic structure at the minimum of the potential energy surface with bond-length alternation, we show here that zero-point vibrations lead to an effective equalization of all carbon-carbon bond lengths and thus to a cumulenic structure.