-
Capturing the electron-electron cusp with the coupling-constant averaged exchange-correlation hole: A case study for Hooke's atoms
Authors:
Lin Hou,
Tom J. P. Irons,
Yanyong Wang,
James W. Furness,
Andrew M. Teale,
Jianwei Sun
Abstract:
In density functional theory the exchange-correlation (XC) energy functional can be defined exactly through the coupling-constant ($λ$) averaged XC hole $\bar{n}_\text{xc}(\mathbf{r},\mathbf{r}')$, representing the probability depletion of finding an electron at $\mathbf{r}$ due to an electron at $\mathbf{r}$. Accurate knowledge of $\bar{n}_\text{xc}(\mathbf{r},\mathbf{r}')$ has been crucial for d…
▽ More
In density functional theory the exchange-correlation (XC) energy functional can be defined exactly through the coupling-constant ($λ$) averaged XC hole $\bar{n}_\text{xc}(\mathbf{r},\mathbf{r}')$, representing the probability depletion of finding an electron at $\mathbf{r}$ due to an electron at $\mathbf{r}$. Accurate knowledge of $\bar{n}_\text{xc}(\mathbf{r},\mathbf{r}')$ has been crucial for developing various XC energy density functional approximations and understanding their performance for real molecules and materials. However, there are very few systems for which accurate XC holes have been calculated, since this requires evaluating the one- and two-particle reduced density matrices for a reference wave function over a range of $λ$ whilst the electron density remains fixed at the physical ($λ=1$) density. Although the coupled-cluster singles and doubles (CCSD) method can yield exact results for a two-electron system in the complete basis set limit, it cannot capture the electron-electron cusp with commonly used finite basis sets. In this study, focusing on the Hooke's atom as a two-electron model system for which certain analytic solutions are known, we examine the effect of this cusp error on the XC hole calculated using CCSD. The Lieb functional is calculated at a range of coupling constants to determine the $λ$-integrated XC hole. Our results indicate that, for the Hooke's atoms, the error introduced by the description of the electron-electron cusp using Gaussian basis sets at the CCSD level is negligible compared to the basis set incompleteness error. The system-, angle- and coupling-constant-averaged XC hole is calculated using the same approach and provides a benchmark against which the Perdew-Burke-Ernzerhof (PBE) and local density approximation (LDA) XC hole models are assessed.
△ Less
Submitted 30 March, 2023;
originally announced March 2023.
-
Exact constraints and appropriate norms in machine learned exchange-correlation functionals
Authors:
Kanun Pokharel,
James W. Furness,
Yi Yao,
Volker Blum,
Tom J. P. Irons,
Andrew M. Teale,
Jianwei Sun
Abstract:
Machine learning techniques have received growing attention as an alternative strategy for developing general-purpose density functional approximations, augmenting the historically successful approach of human designed functionals derived to obey mathematical constraints known for the exact exchange-correlation functional. More recently efforts have been made to reconcile the two techniques, integ…
▽ More
Machine learning techniques have received growing attention as an alternative strategy for developing general-purpose density functional approximations, augmenting the historically successful approach of human designed functionals derived to obey mathematical constraints known for the exact exchange-correlation functional. More recently efforts have been made to reconcile the two techniques, integrating machine learning and exact-constraint satisfaction. We continue this integrated approach, designing a deep neural network that exploits the exact constraint and appropriate norm philosophy to deorbitalize the strongly constrained and appropriately normed SCAN functional. The deep neural network is trained to replicate the SCAN functional from only electron density and local derivative information, avoiding use of the orbital dependent kinetic energy density. The performance and transferability of the machine learned functional are demonstrated for molecular and periodic systems.
△ Less
Submitted 27 May, 2022;
originally announced May 2022.
-
Lieb variation principle in density-functional theory
Authors:
Trygve Helgaker,
Andrew M. Teale
Abstract:
Lieb's convex formulation of density-functional theory is presented in a pedagogical manner, emphasizing its connection to Hohenberg-Kohn theory and to Levy's constrained-search theory. The Hohenberg-Kohn and Lieb variation principles are discussed, highlighting the dual relationship between the ground-state energy and the universal density functional. Applications of the Lieb variation principle…
▽ More
Lieb's convex formulation of density-functional theory is presented in a pedagogical manner, emphasizing its connection to Hohenberg-Kohn theory and to Levy's constrained-search theory. The Hohenberg-Kohn and Lieb variation principles are discussed, highlighting the dual relationship between the ground-state energy and the universal density functional. Applications of the Lieb variation principle are reviewed, demonstrating how it may be utilized to calculate the Kohn-Sham potential of atoms and molecules, to study the exchange-correlation functional and the adiabatic connection by high-precision many-body methods, and to calculate the exchange-correlation hole and energy densities of atoms and molecules.
△ Less
Submitted 26 April, 2022;
originally announced April 2022.
-
Robust All-Electron Optimization in Orbital-Free Density Functional Theory Using the Trust Region Image Method
Authors:
Matthew S. Ryley,
Michael Withnall,
Tom J. P. Irons,
Trygve Helgaker,
Andrew M. Teale
Abstract:
We present a Gaussian-basis implementation of orbital-free density-functional theory (OF-DFT) in which the trust-region image method (TRIM) is used for optimization. This second-order optimization scheme has been constructed to provide benchmark all-electron results with very tight convergence of the particle number constraint, associated chemical potential and electron density. It is demonstrated…
▽ More
We present a Gaussian-basis implementation of orbital-free density-functional theory (OF-DFT) in which the trust-region image method (TRIM) is used for optimization. This second-order optimization scheme has been constructed to provide benchmark all-electron results with very tight convergence of the particle number constraint, associated chemical potential and electron density. It is demonstrated that, by preserving the saddle-point nature of the optimization and simultaneously optimizing the density and chemical potential, an order of magnitude reduction in the number of iterations required for convergence is obtained. The approach is compared and contrasted with a new implementation of the nested optimization scheme put forward by Chan, Cohen and Handy. Our implementation allows for semi-local kinetic-energy (and exchange-correlation) functionals to be handled self-consistently in all-electron calculations. The all-electron Gaussian-basis setting for these calculations will enable direct comparison with a wide range of standard high-accuracy quantum-chemical methods as well as with Kohn-Sham density-functional theory. We expect that the present implementation will provide a useful tool for analysing the performance of approximate kinetic-energy functionals in finite systems.
△ Less
Submitted 22 December, 2020;
originally announced December 2020.
-
Excitation energies from G{ö}rling-Levy perturbation theory along the range-separated adiabatic connection
Authors:
Elisa Rebolini,
Andrew M. Teale,
Trygve Helgaker,
Andreas Savin,
Julien Toulouse
Abstract:
A G{ö}rling-Levy (GL)-based perturbation theory along the range-separated adiabatic connection is assessed for the calculation of electronic excitation energies. In comparison with the Rayleigh-Schr{ö}dinger (RS)-based perturbation theory introduced in a previous work [E. Rebolini, J. Toulouse, A. M. Teale, T. Helgaker, A. Savin, Mol. Phys. 113, 1740 (2015)], this GL-based perturbation theory keep…
▽ More
A G{ö}rling-Levy (GL)-based perturbation theory along the range-separated adiabatic connection is assessed for the calculation of electronic excitation energies. In comparison with the Rayleigh-Schr{ö}dinger (RS)-based perturbation theory introduced in a previous work [E. Rebolini, J. Toulouse, A. M. Teale, T. Helgaker, A. Savin, Mol. Phys. 113, 1740 (2015)], this GL-based perturbation theory keeps the ground-state density constant at each order and thus gives the correct ionization energy at each order. Excitation energies up to first order in the perturbation have been calculated numerically for the helium and beryllium atoms and the hydrogen molecule without introducing any density-functional approximations. In comparison with the RS-based perturbation theory, the present GL-based perturbation theory gives much more accurate excitation energies for Rydberg states but similar excitation energies for valence states.
△ Less
Submitted 28 December, 2017; v1 submitted 23 November, 2017;
originally announced November 2017.
-
Uniform magnetic fields in density-functional theory
Authors:
Erik I. Tellgren,
Andre Laestadius,
Trygve Helgaker,
Simen Kvaal,
Andrew M. Teale
Abstract:
We construct a density-functional formalism adapted to uniform external magnetic fields that is intermediate between conventional Density Functional Theory and Current-Density Functional Theory (CDFT). In the intermediate theory, which we term LDFT, the basic variables are the density, the canonical momentum, and the paramagnetic contribution to the magnetic moment. Both a constrained-search formu…
▽ More
We construct a density-functional formalism adapted to uniform external magnetic fields that is intermediate between conventional Density Functional Theory and Current-Density Functional Theory (CDFT). In the intermediate theory, which we term LDFT, the basic variables are the density, the canonical momentum, and the paramagnetic contribution to the magnetic moment. Both a constrained-search formulation and a convex formulation in terms of Legendre--Fenchel transformations are constructed. Many theoretical issues in CDFT find simplified analogues in LDFT. We prove results concerning $N$-representability, Hohenberg--Kohn-like mappings, existence of minimizers in the constrained-search expression, and a restricted analogue to gauge invariance. The issue of additivity of the energy over non-interacting subsystems, which is qualitatively different in LDFT and CDFT, is also discussed.
△ Less
Submitted 29 September, 2017;
originally announced September 2017.
-
Interpolated energy densities, correlation indicators and lower bounds from approximations to the strong coupling limit of DFT
Authors:
Stefan Vuckovic,
Lucas O. Wagner,
Andrew M. Teale,
Paola Gori-Giorgi
Abstract:
We investigate the construction of approximated exchange-correlation functionals by interpolating locally along the adiabatic connection between the weak- and the strong-coupling regimes, focussing on the effect of using approximate functionals for the strong-coupling energy densities. The gauge problem is avoided by dealing with quantities that are all locally defined in the same way. Using exact…
▽ More
We investigate the construction of approximated exchange-correlation functionals by interpolating locally along the adiabatic connection between the weak- and the strong-coupling regimes, focussing on the effect of using approximate functionals for the strong-coupling energy densities. The gauge problem is avoided by dealing with quantities that are all locally defined in the same way. Using exact ingredients at weak coupling we are able to isolate the error coming from the approximations at strong coupling only. We find that the nonlocal radius model, which retains some of the non-locality of the exact strong-coupling regime, yields very satisfactory results. We also use interpolation models and quantities from the weak- and strong-coupling regimes to define a correlation-type indicator and a lower bound to the exact exchange-correlation energy. Open problems, related to the nature of the local and global slope of the adiabatic connection at weak coupling, are also discussed.
△ Less
Submitted 4 March, 2017;
originally announced March 2017.
-
Exchange-correlation functionals via local interpolation along the adiabatic connection
Authors:
Stefan Vuckovic,
Tom Irons,
Andreas Savin,
Andrew M. Teale,
Paola Gori-Giorgi
Abstract:
The construction of density-functional approximations is explored by modeling the adiabatic connection em locally, using energy densities defined in terms of the electrostatic potential of the exchange-correlation hole. These local models are more amenable to the construction of size-consistent approximations than their global counterparts. In this work we use accurate input local ingredients to a…
▽ More
The construction of density-functional approximations is explored by modeling the adiabatic connection em locally, using energy densities defined in terms of the electrostatic potential of the exchange-correlation hole. These local models are more amenable to the construction of size-consistent approximations than their global counterparts. In this work we use accurate input local ingredients to assess the accuracy of range of local interpolation models against accurate exchange-correlation energy densities. The importance of the strictly-correlated electrons (SCE) functional describing the strong coupling limit is emphasized, enabling the corresponding interpolated functionals to treat strong correlation effects. In addition to exploring the performance of such models numerically for the helium and beryllium isoelectronic series and the dissociation of the hydrogen molecule, an approximate analytic model is presented for the initial slope of the local adiabatic connection. Comparisons are made with approaches based on global models and prospects for future approximations based on the local adiabatic connection are discussed.
△ Less
Submitted 2 May, 2016; v1 submitted 17 February, 2016;
originally announced February 2016.
-
Current Density-Functional Theory using meta-Generalized Gradient Exchange--Correlation Functionals
Authors:
James W. Furness,
Joachim Verbeke,
Erik I. Tellgren,
Stella Stopkowicz,
Ulf Ekström,
Trygve Helgaker,
Andrew M. Teale
Abstract:
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 thes…
▽ More
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.
△ Less
Submitted 29 May, 2015;
originally announced May 2015.
-
Calculating excitation energies by extrapolation along adiabatic connections
Authors:
Elisa Rebolini,
Julien Toulouse,
Andrew M. Teale,
Trygve Helgaker,
Andreas Savin
Abstract:
In this paper, an alternative method to range-separated linear-response time-dependent density-functional theory and perturbation theory is proposed to improve the estimation of the energies of a physical system from the energies of a partially interacting system. Starting from the analysis of the Taylor expansion of the energies of the partially interacting system around the physical system, we…
▽ More
In this paper, an alternative method to range-separated linear-response time-dependent density-functional theory and perturbation theory is proposed to improve the estimation of the energies of a physical system from the energies of a partially interacting system. Starting from the analysis of the Taylor expansion of the energies of the partially interacting system around the physical system, we use an extrapolation scheme to improve the estimation of the energies of the physical system at an intermediate point of the range-separated or linear adiabatic connection where either the electron--electron interaction is scaled or only the long-range part of the Coulomb interaction is included. The extrapolation scheme is first applied to the range-separated energies of the helium and beryllium atoms and of the hydrogen molecule at its equilibrium and stretched geometries. It improves significantly the convergence rate of the energies toward their exact limit with respect to the range-separation parameter. The range-separated extrapolation scheme is compared with a similar approach for the linear adiabatic connection, highlighting the relative strengths and weaknesses of each approach.
△ Less
Submitted 19 March, 2015;
originally announced March 2015.
-
Excitation energies along a range-separated adiabatic connection
Authors:
Elisa Rebolini,
Julien Toulouse,
Andrew M. Teale,
Trygve Helgaker,
Andreas Savin
Abstract:
We present a study of the variation of total energies and excitationenergies along a range-separated adiabatic connection. This connectionlinks the non-interacting Kohn-Sham electronic system to the physicalinteracting system by progressively switching on theelectron-electron interactions whilst simultaneously adjusting aone-electron effective potential so as to keep the ground-statedensity consta…
▽ More
We present a study of the variation of total energies and excitationenergies along a range-separated adiabatic connection. This connectionlinks the non-interacting Kohn-Sham electronic system to the physicalinteracting system by progressively switching on theelectron-electron interactions whilst simultaneously adjusting aone-electron effective potential so as to keep the ground-statedensity constant. The interactions are introduced in arange-dependent manner, first introducing predominantly long-range,and then all-range, interactions as the physical system is approached,as opposed to the conventional adiabatic connection where theinteractions are introduced by globally scaling the standard Coulomb interaction.Reference data are reported for the He and Be atoms and the H2molecule, obtained by calculating the short-range effective potentialat the full configuration-interaction level using Lieb'sLegendre-transform approach. As the strength of the electron-electroninteractions increases, the excitation energies, calculated for thepartially interacting systems along the adiabatic connection, offerincreasingly accurate approximations to the exact excitation energies.Importantly, the excitation energies calculated at an intermediatepoint of the adiabatic connection are much better approximations tothe exact excitation energies than are the corresponding Kohn-Shamexcitation energies. This is particularly evident in situationsinvolving strong static correlation effects and states with multipleexcitation character, such as the dissociating H2 molecule. Theseresults highlight the utility of long-range interacting referencesystems as a starting point for the calculation of excitation energiesand are of interest for developing and analyzing practical approximaterange-separated density-functional methodologies.
△ Less
Submitted 12 December, 2014;
originally announced December 2014.
-
Excited states from range-separated density-functional perturbation theory
Authors:
Elisa Rebolini,
Julien Toulouse,
Andrew M. Teale,
Trygve Helgaker,
Andreas Savin
Abstract:
We explore the possibility of calculating electronic excited states by using perturbation theory along a range-separated adiabatic connection. Starting from the energies of a partially interacting Hamiltonian, a first-order correction is defined with two variants of perturbation theory: a straight-forward perturbation theory, and an extension of the G{ö}rling--Levy one that has the advantage of ke…
▽ More
We explore the possibility of calculating electronic excited states by using perturbation theory along a range-separated adiabatic connection. Starting from the energies of a partially interacting Hamiltonian, a first-order correction is defined with two variants of perturbation theory: a straight-forward perturbation theory, and an extension of the G{ö}rling--Levy one that has the advantage of keeping the ground-state density constant at each order in the perturbation. Only the first, simpler, variant is tested here on the helium and beryllium atoms and on the dihydrogene molecule. The first-order correction within this perturbation theory improves significantly the total ground-and excited-state energies of the different systems. However, the excitation energies are mostly deterio-rated with respect to the zeroth-order ones, which may be explained by the fact that the ionization energy is no longer correct for all interaction strengths. The second variant of the perturbation theory should improve these results but has not been tested yet along the range-separated adiabatic connection.
△ Less
Submitted 12 December, 2014;
originally announced December 2014.
-
Differentiable but exact formulation of density-functional theory
Authors:
Simen Kvaal,
Ulf Ekström,
Andrew M. Teale,
Trygve Helgaker
Abstract:
The universal density functional $F$ of density-functional theory is a complicated and ill-behaved function of the density-in particular, $F$ is not differentiable, making many formal manipulations more complicated. Whilst $F$ has been well characterized in terms of convex analysis as forming a conjugate pair $(E,F)$ with the ground-state energy $E$ via the Hohenberg-Kohn and Lieb variation princi…
▽ More
The universal density functional $F$ of density-functional theory is a complicated and ill-behaved function of the density-in particular, $F$ is not differentiable, making many formal manipulations more complicated. Whilst $F$ has been well characterized in terms of convex analysis as forming a conjugate pair $(E,F)$ with the ground-state energy $E$ via the Hohenberg-Kohn and Lieb variation principles, $F$ is nondifferentiable and subdifferentiable only on a small (but dense) set of its domain. In this article, we apply a tool from convex analysis, Moreau-Yosida regularization, to construct, for any $ε>0$, pairs of conjugate functionals $({}^ε\!E,{}^ε\!F)$ that converge to $(E,F)$ pointwise everywhere as $ε\rightarrow 0^+$, and such that ${}^ε\!F$ is (Fréchet) differentiable. For technical reasons, we limit our attention to molecular electronic systems in a finite but large box. It is noteworthy that no information is lost in the Moreau-Yosida regularization: the physical ground-state energy $E(v)$ is exactly recoverable from the regularized ground-state energy ${}^ε\!E(v)$ in a simple way. All concepts and results pertaining to the original $(E,F)$ pair have direct counterparts in results for $({}^ε\! E, {}^ε\!F)$. The Moreau-Yosida regularization therefore allows for an exact, differentiable formulation of density-functional theory. In particular, taking advantage of the differentiability of ${}^ε\!F$, a rigorous formulation of Kohn-Sham theory is presented that does not suffer from the noninteracting representability problem in standard Kohn-Sham theory.
△ Less
Submitted 12 February, 2014; v1 submitted 13 December, 2013;
originally announced December 2013.
-
A Density Difference Based Analysis of Orbital--Dependent Exchange--Correlation Functionals
Authors:
I. Grabowski,
A. M. Teale,
E. Fabiano,
S. Smiga,
A. Buksztel,
F. Della Sala
Abstract:
We present a density difference based analysis for a range of orbital--dependent Kohn--Sham functionals. Results for atoms, some members of the neon isoelectronic series and small molecules are reported and compared with ab initio wave-function calculations. Particular attention is paid to the quality of approximations to the exchange--only optimized effective potential (OEP) approach: we consider…
▽ More
We present a density difference based analysis for a range of orbital--dependent Kohn--Sham functionals. Results for atoms, some members of the neon isoelectronic series and small molecules are reported and compared with ab initio wave-function calculations. Particular attention is paid to the quality of approximations to the exchange--only optimized effective potential (OEP) approach: we consider both the Localized Hartree Fock as well as the Krieger-Li-Iafrate methods. Analysis of density differences at the exchange--only level reveals the impact the approximations have on the resulting electronic densities. These differences are further quantified in terms of the ground state energies, frontier orbital energy differences and highest occupied orbital energies obtained. At the correlated level an OEP approach based on a perturbative second--order correlation energy expression is shown to deliver results comparable with those from traditional wave function approaches, making it suitable for use as a benchmark against which to compare standard density--functional approximations.
△ Less
Submitted 6 November, 2013;
originally announced November 2013.
-
Alternative separation of exchange and correlation energies in multi-configuration range-separated density-functional theory
Authors:
Alexandrina Stoyanova,
Andrew M. Teale,
Julien Toulouse,
Trygve Helgaker,
Emmanuel Fromager
Abstract:
The alternative separation of exchange and correlation energies proposed by Toulouse et al. [Theor. Chem. Acc. 114, 305 (2005)] is explored in the context of multi-configuration range-separated density-functional theory. The new decomposition of the short-range exchange-correlation energy relies on the auxiliary long-range interacting wavefunction rather than the Kohn-Sham (KS) determinant. The ad…
▽ More
The alternative separation of exchange and correlation energies proposed by Toulouse et al. [Theor. Chem. Acc. 114, 305 (2005)] is explored in the context of multi-configuration range-separated density-functional theory. The new decomposition of the short-range exchange-correlation energy relies on the auxiliary long-range interacting wavefunction rather than the Kohn-Sham (KS) determinant. The advantage, relative to the traditional KS decomposition, is that the wavefunction part of the energy is now computed with the regular (fully-interacting) Hamiltonian. One potential drawback is that, because of double counting, the wavefunction used to compute the energy cannot be obtained by minimizing the energy expression with respect to the wavefunction parameters. The problem is overcome by using short-range optimized effective potentials (OEPs). The resulting combination of OEP techniques with wavefunction theory has been investigated in this work, at the Hartree-Fock (HF) and multi-configuration self-consistent-field (MCSCF) levels. In the HF case, an analytical expression for the energy gradient has been derived and implemented. Calculations have been performed within the short-range local density approximation on H2, N2, Li2 and H2O. Significant improvements in binding energies are obtained with the new decomposition of the short-range energy. The importance of optimizing the short-range OEP at the MCSCF level when static correlation becomes significant has also been demonstrated for H2, using a finite-difference gradient. The implementation of the analytical gradient for MCSCF wavefunctions is currently in progress.
△ Less
Submitted 11 September, 2013; v1 submitted 7 August, 2013;
originally announced August 2013.
-
Analysis of double hybrid density-functionals along the adiabatic connection
Authors:
Yann Cornaton,
Odile Franck,
Andrew M. Teale,
Emmanuel Fromager
Abstract:
We present a graphical analysis of the adiabatic connections underlying double-hybrid density-functional methods that employ second-order perturbation theory. Approximate adiabatic connection formulae relevant to the construction of these functionals are derived and compared directly with those calculated using accurate ab initio methods. The discontinuous nature of the approximate adiabatic integ…
▽ More
We present a graphical analysis of the adiabatic connections underlying double-hybrid density-functional methods that employ second-order perturbation theory. Approximate adiabatic connection formulae relevant to the construction of these functionals are derived and compared directly with those calculated using accurate ab initio methods. The discontinuous nature of the approximate adiabatic integrands is emphasized, the discontinuities occurring at interaction strengths which mark the transitions between regions that are: (i) described predominantly by second- order perturbation theory (ii) described by a mixture of density-functional and second-order perturbation theory contributions and (iii) described purely by density-functional theory. Numerical examples are presented for a selection of small molecular systems and van der Waals dimers. The impacts of commonly used approximations in each of the three sections of the adiabatic connection are discussed along with possible routes for the development of improved double-hybrid methodologies.
△ Less
Submitted 26 February, 2013;
originally announced February 2013.
-
The choice of basic variables in current-density functional theory
Authors:
Erik I. Tellgren,
Simen Kvaal,
Espen Sagvolden,
Ulf Ekström,
Andrew M. Teale,
Trygve Helgaker
Abstract:
The selection of basic variables in current-density functional theory and formal properties of the resulting formulations are critically examined. Focus is placed on the extent to which the Hohenberg--Kohn theorem, constrained-search approach and Lieb's formulation (in terms of convex and concave conjugation) of standard density-functional theory can be generalized to provide foundations for curre…
▽ More
The selection of basic variables in current-density functional theory and formal properties of the resulting formulations are critically examined. Focus is placed on the extent to which the Hohenberg--Kohn theorem, constrained-search approach and Lieb's formulation (in terms of convex and concave conjugation) of standard density-functional theory can be generalized to provide foundations for current-density functional theory. For the well-known case with the gauge-dependent paramagnetic current density as a basic variable, we find that the resulting total energy functional is not concave. It is shown that a simple redefinition of the scalar potential restores concavity and enables the application of convex analysis and convex/concave conjugation. As a result, the solution sets arising in potential-optimization problems can be given a simple characterization. We also review attempts to establish theories with the physical current density as a basic variable. Despite the appealing physical motivation behind this choice of basic variables, we find that the mathematical foundations of the theories proposed to date are unsatisfactory. Moreover, the analogy to standard density-functional theory is substantially weaker as neither the constrained-search approach nor the convex analysis framework carry over to a theory making use of the physical current density.
△ Less
Submitted 8 October, 2012;
originally announced October 2012.