No Arabic abstract
Exact density functionals for the exchange and correlation energies are approximated in practical calculations for the ground-state electronic structure of a many-electron system. An important exact constraint for the construction of approximations is to recover the correct non-relativistic large-$Z$ expansions for the corresponding energies of neutral atoms with atomic number $Z$ and electron number $N=Z$, which are correct to leading order ($-0.221 Z^{5/3}$ and $-0.021 Z ln Z$ respectively) even in the lowest-rung or local density approximation. We find that hydrogenic densities lead to $E_x(N,Z) approx -0.354 N^{2/3} Z$ (as known before only for $Z gg N gg 1$) and $E_c approx -0.02 N ln N$. These asymptotic estimates are most correct for atomic ions with large $N$ and $Z gg N$, but we find that they are qualitatively and semi-quantitatively correct even for small $N$ and for $N approx Z$. The large-$N$ asymptotic behavior of the energy is pre-figured in small-$N$ atoms and atomic ions, supporting the argument that widely-predictive approximate density functionals should be designed to recover the correct asymptotics. It is shown that the exact Kohn-Sham correlation energy, when calculated from the pure ground-state wavefunction, should have no contribution proportional to $Z$ in the $Zto infty$ limit for any fixed $N$.
In this chapter, we provide a review of ground-state Kohn-Sham density-functional theory of electronic systems and some of its extensions, we present exact expressions and constraints for the exchange and correlation density functionals, and we discuss the main families of approximations for the exchange-correlation energy: semilocal approximations, single-determinant hybrid approximations, multideterminant hybrid approximations, dispersion-corrected approximations, as well as orbital-dependent exchange-correlation density functionals. The chapter aims at providing both a consistent birds-eye view of the field and a detailed description of some of the most used approximations. It is intended to be readable by chemists/physicists and applied mathematicians.
We present the self-consistent implementation of current-dependent (hybrid) meta generalized gradient approximation (mGGA) density functionals using London atomic orbitals. A previously proposed generalized kinetic energy density is utilized to implement mGGAs in the framework of Kohn--Sham current density-functional theory (KS-CDFT). A unique feature of the non-perturbative implementation of these functionals is the ability to seamlessly explore a wide range of magnetic fields up to 1 a.u. ($sim 235000$T) in strength. CDFT functionals based on the TPSS and B98 forms are investigated and their performance is assessed by comparison with accurate CCSD(T) data. In the weak field regime magnetic properties such as magnetizabilities and NMR shielding constants show modest but systematic improvements over GGA functionals. However, in strong field regime the mGGA based forms lead to a significantly improved description of the recently proposed perpendicular paramagnetic bonding mechanism, comparing well with CCSD(T) data. In contrast to functionals based on the vorticity these forms are found to be numerically stable and their accuracy at high field suggests the extension of mGGAs to CDFT via the generalized kinetic energy density should provide a useful starting point for further development of CDFT approximations.
We review and expand on our work to impose constraints on the effective Kohn Sham (KS) potential of local and semi-local density functional approximations. In this work, we relax a previously imposed positivity constraint, which increased the computational cost and we find that it is safe to do so, except in systems with very few electrons. The constrained minimisation leads invariably to the solution of an optimised effective potential (OEP) equation in order to determine the KS potential. We review briefly our previous work on this problem and demonstrate with numerous examples that despite well-known mathematical issues of the OEP with finite basis sets, our OEP equations are well behaved. We demonstrate that constraining the screening charge of the Hartree, exchange and correlation potential not only corrects its asymptotic behaviour but also allows the exchange and correlation potential to exhibit nonzero derivative discontinuity, a feature of the exact KS potential that is necessary for the accurate prediction of band-gaps in solids but very hard to capture with semi-local approximations.
A single-term density functional model for nondynamic and strong correlation is presented, based on single-determinant Kohn-Sham density functional theory. It is derived from modeling the adiabatic connection and contains only two nonlinear empirical parameters. Preliminary tests show that the model recovers majority of nondynamic correlation during a molecular dissociation and at the same time performs reasonably for atomization energies. It demonstrates the feasibility of developing DFT functionals for nondynamic and strong correlation within the single-determinant KS scheme.
We investigate the performance of the range-separated hybrid (RSH) scheme, which combines long-range Hartree-Fock (HF) and a short-range density-functional approximation (DFA), for calculating photoexcitation/photoionization spectra of the H and He atoms, using a B-spline basis set in order to correctly describe the continuum part of the spectra. The study of these simple systems allows us to quantify the influence on the spectra of the errors coming from the short-range exchange-correlation DFA and from the missing long-range correlation in the RSH scheme. We study the differences between using the long-range HF exchange (nonlocal) potential and the long-range exact exchange (local) potential. Contrary to the former, the latter supports a series of Rydberg states and gives reasonable photoexcitation/photoionization spectra, even without applying linear-response theory. The most accurate spectra are obtained with the linear-response time-dependent range-separated hybrid (TDRSH) scheme. In particular, for the He atom at the optimal value of the range-separation parameter, TDRSH gives slightly more accurate photoexcitation and photoioniza-tion spectra than standard linear-response time-dependent HF. More generally, the present work shows the potential of range-separated density-functional theory for calculating linear and nonlinear optical properties involving continuum states.