No Arabic abstract
A Neural-Networks-based approach is proposed to construct a new type of exchange-correlation functional for density functional theory. It is applied to improve B3LYP functional by taking into account of high-order contributions to the exchange-correlation functional. The improved B3LYP functional is based on a neural network whose structure and synaptic weights are determined from 116 known experimental atomization energies, ionization potentials, proton affinities or total atomic energies which were used by Becke in his pioneer work on the hybrid functionals [J. Chem. Phys. ${bf 98}$, 5648 (1993)]. It leads to better agreement between the first-principles calculation results and these 116 experimental data. The new B3LYP functional is further tested by applying it to calculate the ionization potentials of 24 molecules of the G2 test set. The 6-311+G(3{it df},2{it p}) basis set is employed in the calculation, and the resulting root-mean-square error is reduced to 2.2 kcal$cdot$mol$^{-1}$ in comparison to 3.6 kcal$cdot$mol$^{-1}$ of conventional B3LYP/6-311+G(3{it df},2{it p}) calculation.
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.
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 train a neural network as the universal exchange-correlation functional of density-functional theory that simultaneously reproduces both the exact exchange-correlation energy and potential. This functional is extremely non-local, but retains the computational scaling of traditional local or semi-local approximations. It therefore holds the promise of solving some of the delocalization problems that plague density-functional theory, while maintaining the computational efficiency that characterizes the Kohn-Sham equations. Furthermore, by using automatic differentiation, a capability present in modern machine-learning frameworks, we impose the exact mathematical relation between the exchange-correlation energy and the potential, leading to a fully consistent method. We demonstrate the feasibility of our approach by looking at one-dimensional systems with two strongly-correlated electrons, where density-functional methods are known to fail, and investigate the behavior and performance of our functional by varying the degree of non-locality.
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$.