-
Extended $N$-centered ensemble density functional theory of double electronic excitations
Authors:
Filip Cernatic,
Emmanuel Fromager
Abstract:
A recent work [arXiv:2401.04685] has merged $N$-centered ensembles of neutral and charged electronic ground states with ensembles of neutral ground and excited states, thus providing a general and in-principle exact (so-called extended $N$-centered) ensemble density functional theory of neutral and charged electronic excitations. This formalism made it possible to revisit the concept of density-fu…
▽ More
A recent work [arXiv:2401.04685] has merged $N$-centered ensembles of neutral and charged electronic ground states with ensembles of neutral ground and excited states, thus providing a general and in-principle exact (so-called extended $N$-centered) ensemble density functional theory of neutral and charged electronic excitations. This formalism made it possible to revisit the concept of density-functional derivative discontinuity, in the particular case of single excitations from the highest occupied Kohn-Sham (KS) molecular orbital, without invoking the usual "asymptotic behavior of the density" argument. In this work, we address a broader class of excitations, with a particular focus on double excitations. An exact implementation of the theory is presented for the two-electron Hubbard dimer model. A thorough comparison of the true physical ground- and excited-state electronic structures with that of the fictitious ensemble density-functional KS system is also presented. Depending on the choice of the density-functional ensemble as well as the asymmetry of the dimer and the correlation strength, an inversion of states can be observed. In some other cases, the strong mixture of KS states within the true physical system makes the assignment "single excitation" or "double excitation" irrelevant.
△ Less
Submitted 11 February, 2024;
originally announced February 2024.
-
Neutral electronic excitations and derivative discontinuities: An extended $N$-centered ensemble density functional theory perspective
Authors:
Filip Cernatic,
Pierre-François Loos,
Bruno Senjean,
Emmanuel Fromager
Abstract:
This work merges two different types of many-electron ensembles, namely the Theophilou-Gross-Oliveira-Kohn ensembles of ground and neutrally-excited states, and the more recent $N$-centered ensembles of neutral and charged ground states. On that basis, an in-principle exact and general, so-called extended $N$-centered, ensemble density-functional theory of charged and neutral electronic excitation…
▽ More
This work merges two different types of many-electron ensembles, namely the Theophilou-Gross-Oliveira-Kohn ensembles of ground and neutrally-excited states, and the more recent $N$-centered ensembles of neutral and charged ground states. On that basis, an in-principle exact and general, so-called extended $N$-centered, ensemble density-functional theory of charged and neutral electronic excitations is derived. We revisit in this context the concept of density-functional derivative discontinuity for neutral excitations, without ever invoking nor using the asymptotic behavior of the ensemble electronic density. The present mathematical construction fully relies on the weight dependence of the ensemble Hartree-exchange-correlation density-functional energy, which makes the theory applicable to lattice models and opens new perspectives for the description of gaps in mesoscopic systems.
△ Less
Submitted 27 April, 2024; v1 submitted 9 January, 2024;
originally announced January 2024.
-
Density functional theory beyond the Born-Oppenheimer approximation: Exact mapping onto an electronically non-interacting Kohn-Sham molecule
Authors:
Emmanuel Fromager,
Benjamin Lasorne
Abstract:
This work presents an alternative, general, and in-principle exact extension of electronic Kohn-Sham density functional theory (KS-DFT) to the fully quantum-mechanical molecular problem. Unlike in existing multi-component or exact-factorization-based DFTs of electrons and nuclei, both nuclear and electronic densities are mapped onto a fictitious electronically non-interacting molecule (referred to…
▽ More
This work presents an alternative, general, and in-principle exact extension of electronic Kohn-Sham density functional theory (KS-DFT) to the fully quantum-mechanical molecular problem. Unlike in existing multi-component or exact-factorization-based DFTs of electrons and nuclei, both nuclear and electronic densities are mapped onto a fictitious electronically non-interacting molecule (referred to as KS molecule), where the electrons still interact with the nuclei. Moreover, in the present molecular KS-DFT, no assumption is made about the mathematical form (exactly factorized or not) of the molecular wavefunction. By expanding the KS molecular wavefunction à la Born-Huang, we obtain a self-consistent set of "KS beyond Born-Oppenheimer" electronic equations coupled to nuclear equations that describe nuclei interacting among themselves and with non-interacting electrons. An exact adiabatic connection formula is derived for the Hartree-exchange-correlation energy of the electrons within the molecule and, on that basis, a practical adiabatic density-functional approximation is proposed and discussed.
△ Less
Submitted 2 April, 2024; v1 submitted 22 December, 2023;
originally announced December 2023.
-
Unitary transformations within density matrix embedding approaches: A novel perspective on the self-consistent scheme for electronic structure calculation
Authors:
Quentin Marécat,
Benjamin Lasorne,
Emmanuel Fromager,
Matthieu Saubanère
Abstract:
In this work, we introduce an original self-consistent scheme based on the one-body reduced density matrix ($γ$) formalism. A significant feature of this methodology is the utilization of an optimal unitary transformation of the Hamiltonian, determined through a self-consistently determined, unitary reflection $\mathbf{R}[γ]$. This enables the extraction of all reduced properties of the system fro…
▽ More
In this work, we introduce an original self-consistent scheme based on the one-body reduced density matrix ($γ$) formalism. A significant feature of this methodology is the utilization of an optimal unitary transformation of the Hamiltonian, determined through a self-consistently determined, unitary reflection $\mathbf{R}[γ]$. This enables the extraction of all reduced properties of the system from a smaller, accurately solved embedding cluster, and to systematically reconstruct the reduced density matrix of the system. This process ensures that both extended and embedded systems satisfy the local virial-like relation, providing quantitative insight into the correspondence between the fragment in the extended system and its embedded analogue. The performance and convergence of the method, as well as the N-representability of the resulting correlated density matrix, are evaluated and discussed within the context of the one-dimensional Hubbard model, which provides exact results for a comprehensive comparison.
△ Less
Submitted 25 September, 2023; v1 submitted 13 June, 2023;
originally announced June 2023.
-
A unified density-matrix functional construction of quantum baths in density matrix embedding theory beyond the mean-field approximation
Authors:
Sajanthan Sekaran,
Oussama Bindech,
Emmanuel Fromager
Abstract:
The equivalence in one-electron quantum bath between the practical implementation of density matrix embedding theory (DMET) and the more recent Householder-transformed density matrix functional embedding theory has been shown previously in the standard but special case where the reference full-size (one-electron reduced) density matrix, from which the bath is constructed, is idempotent [J. Chem. P…
▽ More
The equivalence in one-electron quantum bath between the practical implementation of density matrix embedding theory (DMET) and the more recent Householder-transformed density matrix functional embedding theory has been shown previously in the standard but special case where the reference full-size (one-electron reduced) density matrix, from which the bath is constructed, is idempotent [J. Chem. Phys. 157, 214112 (2022)]. We prove mathematically that the equivalence remains valid when the density matrix is not idempotent anymore, thus allowing for the construction of correlated (one-electron) quantum baths. A density-matrix functional exactification of DMET is derived within the present unified quantum embedding formalism. Numerical examples reveal that the embedding cluster can be quite sensitive to the level of density-matrix functional approximation used for computing the reference density matrix.
△ Less
Submitted 27 June, 2023; v1 submitted 28 April, 2023;
originally announced April 2023.
-
Quantum embedding of multi-orbital fragments using the Block-Householder-transformation
Authors:
Saad Yalouz,
Sajanthan Sekaran,
Emmanuel Fromager,
Matthieu Saubanère
Abstract:
Recently, some of the authors introduced the use of the Householder transformation as a simple and intuitive method for the embedding of local molecular fragments (see Sekaran et. al., Phys. Rev. B 104, 035121 (2021), and Sekaran et. al., Computation 10, 45 (2022)). In this work, we present an extension of this approach to the more general case of multi-orbital fragments using the block version of…
▽ More
Recently, some of the authors introduced the use of the Householder transformation as a simple and intuitive method for the embedding of local molecular fragments (see Sekaran et. al., Phys. Rev. B 104, 035121 (2021), and Sekaran et. al., Computation 10, 45 (2022)). In this work, we present an extension of this approach to the more general case of multi-orbital fragments using the block version of the Householder transformation applied to the one-body reduced density matrix, yet unlocking the applicability to general quantum chemistry/condensed-matter physics Hamiltonians.
A step by step construction of the Block-Householder transformation is presented. Both physical and numerical interest of the approach are highlighted. The specific mean-field (non-interacting) case is thoroughly detailed as it is shown that the embedding of a given $N$ spin-orbitals fragment leads to the generation of two separated sub-systems: a $2N$ spin-orbitals "fragment+bath" cluster that exactly contains $N$ electrons, and a remaining cluster's "environment" which is described by so-called core electrons. We illustrate the use of this transformation in different cases of embedding {scheme} for practical applications. We particularly focus on the extension of the previously introduced Local Potential Functional Embedding Theory (LPFET) and Householder-transformed Density Matrix Functional Embedding Theory (Ht-DMFET) to the case of multi-orbital fragments. These calculations are realized on different types of systems such as model Hamiltonians (Hubbard rings) and \textit{ab initio} molecular systems (hydrogen rings).
△ Less
Submitted 5 November, 2022; v1 submitted 21 September, 2022;
originally announced September 2022.
-
Reduced density matrix functional theory from an ab initio seniority-zero wave function: Exact and approximate formulations along adiabatic connection paths
Authors:
Bruno Senjean,
Saad Yalouz,
Naoki Nakatani,
Emmanuel Fromager
Abstract:
Currently, there is a growing interest in the development of a new hierarchy of methods based on the concept of seniority, which has been introduced quite recently in quantum chemistry. Despite the enormous potential of these methods, the accurate description of both dynamical and static correlation effects within a single and in-principle-exact approach remains a challenge. In this work, we propo…
▽ More
Currently, there is a growing interest in the development of a new hierarchy of methods based on the concept of seniority, which has been introduced quite recently in quantum chemistry. Despite the enormous potential of these methods, the accurate description of both dynamical and static correlation effects within a single and in-principle-exact approach remains a challenge. In this work, we propose an alternative formulation of reduced density-matrix functional theory (RDMFT) where the (one-electron reduced) density matrix is mapped onto an ab initio seniority-zero wave function. In this theory, the exact natural orbitals and their occupancies are determined self-consistently from an effective seniority-zero calculation. The latter involves a universal higher-seniority density matrix functional for which an adiabatic connection (AC) formula is derived and implemented under specific constraints that are related to the density matrix. The pronounced curvature of the (constrained) AC integrand, which is numerically observed in prototypical hydrogen chains and the Helium dimer, indicates that a description of higher-seniority correlations within second-order perturbation theory is inadequate in this context. Applying multiple linear interpolations along the AC or connecting second-order perturbation theory to a full-seniority treatment via Padé approximants are better strategies. Such information is expected to serve as a guide in the future design of higher-seniority density-matrix functional approximations.
△ Less
Submitted 12 August, 2022; v1 submitted 1 April, 2022;
originally announced April 2022.
-
Local Potential Functional Embedding Theory: A Self-Consistent Flavor of Density Functional Theory for Lattices without Density Functionals
Authors:
Sajanthan Sekaran,
Matthieu Saubanère,
Emmanuel Fromager
Abstract:
The recently proposed Householder transformed density-matrix functional embedding theory (Ht-DMFET) [Sekaran et al., Phys. Rev. B 104, 035121 (2021)], which is equivalent to (but formally simpler than) density matrix embedding theory (DMET) in the non-interacting case, is revisited from the perspective of density-functional theory (DFT). An in-principle-exact density-functional version of Ht-DMFET…
▽ More
The recently proposed Householder transformed density-matrix functional embedding theory (Ht-DMFET) [Sekaran et al., Phys. Rev. B 104, 035121 (2021)], which is equivalent to (but formally simpler than) density matrix embedding theory (DMET) in the non-interacting case, is revisited from the perspective of density-functional theory (DFT). An in-principle-exact density-functional version of Ht-DMFET is derived for the one-dimensional Hubbard lattice with a single embedded impurity. On the basis of well-identified density-functional approximations, a local potential functional embedding theory (LPFET) is formulated and implemented. Even though LPFET performs better than Ht-DMFET in the low-density regime, in particular when electron correlation is strong, both methods are unable to describe the density-driven Mott-Hubbard transition, as expected. These results combined with our formally exact density-functional embedding theory reveal that a single statically embedded impurity can in principle describe the gap opening, provided that the complementary correlation potential (that describes the interaction of the embedding cluster with its environment, which is simply neglected in both Ht-DMFET and LPFET) exhibits a derivative discontinuity (DD) at half filling. The extension of LPFET to multiple impurities (which would enable to circumvent the modeling of DDs) and its generalization to quantum chemical Hamiltonians are left for future work.
△ Less
Submitted 16 February, 2022;
originally announced February 2022.
-
Ensemble Density Functional Theory of Neutral and Charged Excitations
Authors:
Filip Cernatic,
Bruno Senjean,
Vincent Robert,
Emmanuel Fromager
Abstract:
Recent progress in the field of (time-independent) ensemble density-functional theory (DFT) for excited states are reviewed. Both Gross-Oliveira-Kohn (GOK) and $N$-centered ensemble formalisms, which are mathematically very similar and allow for an in-principle-exact description of neutral and charged electronic excitations, respectively, are discussed. Key exact results like, for example, the equ…
▽ More
Recent progress in the field of (time-independent) ensemble density-functional theory (DFT) for excited states are reviewed. Both Gross-Oliveira-Kohn (GOK) and $N$-centered ensemble formalisms, which are mathematically very similar and allow for an in-principle-exact description of neutral and charged electronic excitations, respectively, are discussed. Key exact results like, for example, the equivalence between the infamous derivative discontinuity problem and the description of weight dependencies in the ensemble exchange-correlation density functional, are highlighted. The variational evaluation of orbital-dependent ensemble Hartree-exchange (Hx) energies is discussed in detail. We show in passing that state-averaging individual exact Hx energies can lead to severe (solvable though) $v$-representability issues. Finally, we explore the possibility to use the concept of density-driven correlation, which has been recently introduced and does not exist in regular ground-state DFT, for improving state-of-the-art correlation density-functional approximations for ensembles. The present review reflects the efforts of a growing community to turn ensemble DFT into a rigorous and reliable low-cost computational method for excited states. We hope that, in the near future, this contribution will stimulate new formal and practical developments in the field.
△ Less
Submitted 1 November, 2021; v1 submitted 10 September, 2021;
originally announced September 2021.
-
Householder transformed density matrix functional embedding theory
Authors:
Sajanthan Sekaran,
Masahisa Tsuchiizu,
Matthieu Saubanère,
Emmanuel Fromager
Abstract:
Quantum embedding based on the (one-electron reduced) density matrix is revisited by means of the unitary Householder transformation. While being exact and equivalent to (but formally simpler than) density matrix embedding theory (DMET) in the non-interacting case, the resulting Householder transformed density matrix functional embedding theory (Ht-DMFET) preserves, by construction, the single-par…
▽ More
Quantum embedding based on the (one-electron reduced) density matrix is revisited by means of the unitary Householder transformation. While being exact and equivalent to (but formally simpler than) density matrix embedding theory (DMET) in the non-interacting case, the resulting Householder transformed density matrix functional embedding theory (Ht-DMFET) preserves, by construction, the single-particle character of the bath when electron correlation is introduced. In Ht-DMFET, the projected "impurity+bath" cluster's Hamiltonian (from which approximate local properties of the interacting lattice can be extracted) becomes an explicit functional of the density matrix. In the spirit of single-impurity DMET, we consider in this work a closed (two-electron) cluster constructed from the full-size non-interacting density matrix. When the (Householder transformed) interaction on the bath site is taken into account, per-site energies obtained for the half-filled one-dimensional Hubbard lattice match almost perfectly the exact Bethe Ansatz results in all correlation regimes. In the strongly correlated regime, the results deteriorate away from half-filling. This can be related to the electron number fluctuations in the (two-site) cluster which are not described neither in Ht-DMFET nor in regular DMET. As expected, the per-site energies dramatically improve when increasing the number of embedded impurities. Formal connections with density/density matrix functional theories have been briefly discussed and should be explored further. Work is currently in progress in this direction.
△ Less
Submitted 8 June, 2021; v1 submitted 6 March, 2021;
originally announced March 2021.
-
Exact exchange-correlation potentials for calculating the fundamental gap with a fixed number of electrons
Authors:
M. J. P. Hodgson,
J. Wetherell,
Emmanuel Fromager
Abstract:
Capturing the discontinuous shift by $Δ$ in the exact exchange-correlation (xc) potential is the standard proposal for calculating the fundamental gap, $E_\mathrm{g}$, from the Kohn-Sham (KS) gap, $\varepsilon_\mathrm{g}$, within KS density functional theory (DFT), as $E_\mathrm{g} = \varepsilon_\mathrm{g} + Δ$, yet this discontinuity is absent from existing approximations. The '$N$-centered' form…
▽ More
Capturing the discontinuous shift by $Δ$ in the exact exchange-correlation (xc) potential is the standard proposal for calculating the fundamental gap, $E_\mathrm{g}$, from the Kohn-Sham (KS) gap, $\varepsilon_\mathrm{g}$, within KS density functional theory (DFT), as $E_\mathrm{g} = \varepsilon_\mathrm{g} + Δ$, yet this discontinuity is absent from existing approximations. The '$N$-centered' formulation of ensemble DFT artificially maintains a total electron number, $N$, in order to yield $E_\mathrm{g}$ not through a discontinuous shift in the xc potential but via the ensemble-weight derivative of the xc energy. Within the $N$-centered approach we calculate exact xc potentials for a one-dimensional finite system and show analytically that $Δ$ can in fact be interpreted as a discontinuous shift in the exact $N$-centered ensemble xc potential, thereby extending to charged excitations an exact property of uncharged excitations. We show that applying the Levy-Zahariev 'shift-in-potential' procedure in this context relocates the discontinuous shift to the unimportant periphery of the system, so that the exact xc potential in effect is free of discontinuities and thus the inability of a local functional to capture discontinuous behavior is inconsequential.
△ Less
Submitted 13 January, 2021; v1 submitted 12 October, 2020;
originally announced October 2020.
-
Weight Dependence of Local Exchange-Correlation Functionals in Ensemble Density-Functional Theory: Double Excitations in Two-Electron Systems
Authors:
Clotilde Marut,
Bruno Senjean,
Emmanuel Fromager,
Pierre-François Loos
Abstract:
Gross--Oliveira--Kohn (GOK) ensemble density-functional theory (GOK-DFT) is a time-\textit{independent} extension of density-functional theory (DFT) which allows to compute excited-state energies via the derivatives of the ensemble energy with respect to the ensemble weights. Contrary to the time-dependent version of DFT (TD-DFT), double excitations can be easily computed within GOK-DFT. However,…
▽ More
Gross--Oliveira--Kohn (GOK) ensemble density-functional theory (GOK-DFT) is a time-\textit{independent} extension of density-functional theory (DFT) which allows to compute excited-state energies via the derivatives of the ensemble energy with respect to the ensemble weights. Contrary to the time-dependent version of DFT (TD-DFT), double excitations can be easily computed within GOK-DFT. However, to take full advantage of this formalism, one must have access to a \textit{weight-dependent} exchange-correlation functional in order to model the infamous ensemble derivative contribution to the excitation energies. In the present article, we discuss the construction of first-rung (i.e., local) weight-dependent exchange-correlation density-functional approximations for two-electron atomic and molecular systems (He and H$_2$) specifically designed for the computation of double excitations within GOK-DFT. In the spirit of optimally-tuned range-separated hybrid functionals, a two-step system-dependent procedure is proposed to obtain accurate energies associated with double excitations.
△ Less
Submitted 16 June, 2020; v1 submitted 13 May, 2020;
originally announced May 2020.
-
A weight-dependent local correlation density-functional approximation for ensembles
Authors:
Pierre-François Loos,
Emmanuel Fromager
Abstract:
We report a local, weight-dependent correlation density-functional approximation that incorporates information about both ground and excited states in the context of density-functional theory for ensembles (eDFT). This density-functional approximation for ensembles is specially designed for the computation of single and double excitations within Gross--Oliveira--Kohn (GOK) DFT (i.e., eDFT for neut…
▽ More
We report a local, weight-dependent correlation density-functional approximation that incorporates information about both ground and excited states in the context of density-functional theory for ensembles (eDFT). This density-functional approximation for ensembles is specially designed for the computation of single and double excitations within Gross--Oliveira--Kohn (GOK) DFT (i.e., eDFT for neutral excitations), and can be seen as a natural extension of the ubiquitous local-density approximation in the context of ensembles. The resulting density-functional approximation, based on both finite and infinite uniform electron gas models, automatically incorporates the infamous derivative discontinuity contributions to the excitation energies through its explicit ensemble weight dependence. Its accuracy is illustrated by computing single and double excitations in one-dimensional many-electron systems in the weak, intermediate and strong correlation regimes. Although the present weight-dependent functional has been specifically designed for one-dimensional systems, the methodology proposed here is general, i.e., directly applicable to the construction of weight-dependent functionals for realistic three-dimensional systems, such as molecules and solids.
△ Less
Submitted 8 May, 2020; v1 submitted 11 March, 2020;
originally announced March 2020.
-
Individual correlations in ensemble density-functional theory: State-driven/density-driven decompositions without additional Kohn-Sham systems
Authors:
Emmanuel Fromager
Abstract:
Gould and Pittalis [Phys. Rev. Lett. 123, 016401 (2019)] recently revealed a density-driven (DD) correlation energy that is specific to many-electron ensembles and must be accounted for by approximations. We derive in this Letter a general and simpler expression in terms of the ensemble weights, the ensemble Kohn-Sham (KS) orbitals, and their linear response to variations in the ensemble weights.…
▽ More
Gould and Pittalis [Phys. Rev. Lett. 123, 016401 (2019)] recently revealed a density-driven (DD) correlation energy that is specific to many-electron ensembles and must be accounted for by approximations. We derive in this Letter a general and simpler expression in terms of the ensemble weights, the ensemble Kohn-Sham (KS) orbitals, and their linear response to variations in the ensemble weights. As no additional state-driven KS systems are needed, its evaluation is greatly simplified. We confirm the importance of DD effects and introduce a direct and promising route to approximations.
△ Less
Submitted 27 May, 2020; v1 submitted 23 January, 2020;
originally announced January 2020.
-
N-centered ensemble density-functional theory for open systems
Authors:
Bruno Senjean,
Emmanuel Fromager
Abstract:
Two (so-called left and right) variants of N-centered ensemble density-functional theory (DFT) [Senjean and Fromager, Phys. Rev. A 98, 022513 (2018)] are presented. Unlike the original formulation of the theory, these variants allow for the description of systems with a fractional electron number. While conventional DFT for open systems uses only the true electron density as basic variable, left/r…
▽ More
Two (so-called left and right) variants of N-centered ensemble density-functional theory (DFT) [Senjean and Fromager, Phys. Rev. A 98, 022513 (2018)] are presented. Unlike the original formulation of the theory, these variants allow for the description of systems with a fractional electron number. While conventional DFT for open systems uses only the true electron density as basic variable, left/right N-centered ensemble DFT relies instead on (i) a fictitious ensemble density that integrates to a central (integral) number N of electrons, and (ii) a grand canonical ensemble weight $α$ which is equal to the deviation of the true electron number from N. Within such a formalism, the infamous derivative discontinuity that appears when crossing an integral number of electrons is described exactly through the dependence in $α$ of the left and right N-centered ensemble Hartree-exchange-correlation density functionals. Incorporating N-centered ensembles into existing density-functional embedding theories is expected to pave the way towards the in-principle-exact description of an open fragment by means of a pure-state N-electron many-body wavefunction. Work is currently in progress in this direction
△ Less
Submitted 1 February, 2020; v1 submitted 15 December, 2019;
originally announced December 2019.
-
Site-occupation--Green's function embedding theory: A density-functional approach to dynamical impurity solvers
Authors:
Laurent Mazouin,
Matthieu Saubanère,
Emmanuel Fromager
Abstract:
A reformulation of site-occupation embedding theory (SOET) in terms of Green's functions is presented. Referred to as site-occupation--Green's function embedding theory (SOGET), this novel extension of density-functional theory for model Hamiltonians shares many features with dynamical mean-field theory (DMFT) but is formally exact (in any dimension). In SOGET, the impurity-interacting correlation…
▽ More
A reformulation of site-occupation embedding theory (SOET) in terms of Green's functions is presented. Referred to as site-occupation--Green's function embedding theory (SOGET), this novel extension of density-functional theory for model Hamiltonians shares many features with dynamical mean-field theory (DMFT) but is formally exact (in any dimension). In SOGET, the impurity-interacting correlation potential becomes a density-functional self-energy which is frequency-dependent and in principle non-local. A simple local density-functional approximation (LDA) combining the Bethe Ansatz (BA) LDA with the self-energy of the two-level Anderson model is constructed and successfully applied to the one-dimensional Hubbard model. Unlike in previous implementations of SOET, no many-body wavefunction is needed, thus reducing drastically the computational cost of the method.
△ Less
Submitted 10 October, 2019; v1 submitted 2 August, 2019;
originally announced August 2019.
-
Ground and excited energy levels can be extracted exactly from a single ensemble density-functional theory calculation
Authors:
Killian Deur,
Emmanuel Fromager
Abstract:
Gross-Oliveira-Kohn density-functional theory (GOK-DFT) for ensembles is the DFT analog of state-averaged wavefunction-based (SA-WF) methods. In GOK-DFT, the state-averaged (so-called ensemble) exchange-correlation (xc) energy is described by a single functional of the density which, for a fixed density, depends on the weights assigned to each state in the ensemble. We show that, if a many-weight-…
▽ More
Gross-Oliveira-Kohn density-functional theory (GOK-DFT) for ensembles is the DFT analog of state-averaged wavefunction-based (SA-WF) methods. In GOK-DFT, the state-averaged (so-called ensemble) exchange-correlation (xc) energy is described by a single functional of the density which, for a fixed density, depends on the weights assigned to each state in the ensemble. We show that, if a many-weight-dependent xc functional is employed, then it becomes possible to extract, in principle exactly, all individual energy levels from a single GOK-DFT calculation, exactly like in a SA-WF calculation. More precisely, starting from the Kohn-Sham energies, a global Levy-Zahariev-type shift as well as a state-specific (ensemble-based) xc derivative correction must be applied in order to reach the energy level of interest. We illustrate with the asymmetric Hubbard dimer the importance and substantial weight dependence of both corrections. A comparison with more standard extraction procedures, which rely on a sequence of ensemble calculations, is made at the ensemble exact exchange level of approximation.
△ Less
Submitted 13 February, 2019; v1 submitted 6 December, 2018;
originally announced December 2018.
-
Multiple impurities and combined local density approximations in Site-Occupation Embedding Theory
Authors:
Bruno Senjean,
Naoki Nakatani,
Masahisa Tsuchiizu,
Emmanuel Fromager
Abstract:
Site-occupation embedding theory (SOET) is an in-principle-exact multi-determinantal extension of density-functional theory for model Hamiltonians. Various extensions of recent developments in SOET [Senjean et al., Phys. Rev. B 97, 235105 (2018)] are explored in this work. An important step forward is the generalization of the theory to multiple impurity sites. We also propose a new single-impurit…
▽ More
Site-occupation embedding theory (SOET) is an in-principle-exact multi-determinantal extension of density-functional theory for model Hamiltonians. Various extensions of recent developments in SOET [Senjean et al., Phys. Rev. B 97, 235105 (2018)] are explored in this work. An important step forward is the generalization of the theory to multiple impurity sites. We also propose a new single-impurity density-functional approximation (DFA) where the density-functional impurity correlation energy of the two-level (2L) Hubbard system is combined with the Bethe ansatz local density approximation (BALDA) to the full correlation energy of the (infinite) Hubbard model. In order to test the new DFAs, the impurity-interacting wavefunction has been computed self-consistently with the density matrix renormalization group method (DMRG). Double occupation and per-site energy expressions have been derived and implemented in the one-dimensional case. A detailed analysis of the results is presented, with a particular focus on the errors induced either by the energy functionals solely or by the self-consistently converged densities. Among all the DFAs (including those previously proposed), the combined 2L-BALDA is the one that performs the best in all correlation and density regimes. Finally, extensions in new directions, like a partition-DFT-type reformulation of SOET, a projection-based SOET approach, or the combination of SOET with Green functions, are briefly discussed as a perspective.
△ Less
Submitted 5 August, 2019; v1 submitted 15 June, 2018;
originally announced June 2018.
-
Unified formulation of fundamental and optical gap problems in density-functional theory for ensembles
Authors:
Bruno Senjean,
Emmanuel Fromager
Abstract:
Solving the fundamental and optical gap problems, which yield information about charged and neutral excitations in electronic systems, is one of the biggest challenge in density-functional theory (DFT). Despite their intrinsic difference, we show that the two problems can be made formally identical by introducing a universal and canonical ensemble weight dependent exchange-correlation (xc) density…
▽ More
Solving the fundamental and optical gap problems, which yield information about charged and neutral excitations in electronic systems, is one of the biggest challenge in density-functional theory (DFT). Despite their intrinsic difference, we show that the two problems can be made formally identical by introducing a universal and canonical ensemble weight dependent exchange-correlation (xc) density functional. The weight dependence of the xc energy turns out to be the key ingredient for describing the infamous derivative discontinuity and represents a new path for its approximation.
△ Less
Submitted 26 August, 2018; v1 submitted 22 May, 2018;
originally announced May 2018.
-
Exploring weight-dependent density-functional approximations for ensembles in the Hubbard dimer
Authors:
Killian Deur,
Laurent Mazouin,
Bruno Senjean,
Emmanuel Fromager
Abstract:
Gross-Oliveira-Kohn density-functional theory (GOK-DFT) is an extension of DFT to excited states where the basic variable is the ensemble density, i.e. the weighted sum of ground- and excited-state densities. The ensemble energy (i.e. the weighted sum of ground- and excited-state energies) can be obtained variationally as a functional of the ensemble density. Like in DFT, the key ingredient to mod…
▽ More
Gross-Oliveira-Kohn density-functional theory (GOK-DFT) is an extension of DFT to excited states where the basic variable is the ensemble density, i.e. the weighted sum of ground- and excited-state densities. The ensemble energy (i.e. the weighted sum of ground- and excited-state energies) can be obtained variationally as a functional of the ensemble density. Like in DFT, the key ingredient to model in GOK-DFT is the exchange-correlation functional. Developing density-functional approximations (DFAs) for ensembles is a complicated task as both density and weight dependencies should in principle be reproduced. In a recent paper [Phys. Rev. B 95, 035120 (2017)], the authors applied exact GOK-DFT to the simple but nontrivial Hubbard dimer in order to investigate (numerically) the importance of weight dependence in the calculation of excitation energies. In this work, we derive analytical DFAs for various density and correlation regimes by means of a Legendre-Fenchel transform formalism. Both functional and density driven errors are evaluated for each DFA. Interestingly, when the ensemble exact-exchange-only functional is used, these errors can be large, in particular if the dimer is symmetric, but they cancel each other so that the excitation energies obtained by linear interpolation are always accurate, even in the strongly correlated regime.
△ Less
Submitted 10 May, 2018; v1 submitted 1 March, 2018;
originally announced March 2018.
-
Site-Occupation Embedding Theory using Bethe Ansatz Local Density Approximations
Authors:
Bruno Senjean,
Naoki Nakatani,
Masahisa Tsuchiizu,
Emmanuel Fromager
Abstract:
Site-occupation embedding theory (SOET) is an alternative formulation of density-functional theory (DFT) for model Hamiltonians where the fully-interacting Hubbard problem is mapped, in principle exactly, onto an impurity-interacting (rather than a non-interacting) one. It provides a rigorous framework for combining wavefunction (or Green function) based methods with DFT. In this work, exact expre…
▽ More
Site-occupation embedding theory (SOET) is an alternative formulation of density-functional theory (DFT) for model Hamiltonians where the fully-interacting Hubbard problem is mapped, in principle exactly, onto an impurity-interacting (rather than a non-interacting) one. It provides a rigorous framework for combining wavefunction (or Green function) based methods with DFT. In this work, exact expressions for the per-site energy and double occupation of the uniform Hubbard model are derived in the context of SOET. As readily seen from these derivations, the so-called bath contribution to the per-site correlation energy is, in addition to the latter, the key density functional quantity to model in SOET. Various approximations based on Bethe ansatz and perturbative solutions to the Hubbard and single impurity Anderson models are constructed and tested on a one-dimensional ring. The self-consistent calculation of the embedded impurity wavefunction has been performed with the density matrix renormalization group method. It has been shown that promising results are obtained in specific regimes of correlation and density. Possible further developments have been proposed in order to provide reliable embedding functionals and potentials.
△ Less
Submitted 10 April, 2018; v1 submitted 9 October, 2017;
originally announced October 2017.
-
Combining extrapolation with ghost interaction correction in range-separated ensemble density functional theory for excited states
Authors:
Md. Mehboob Alam,
Killian Deur,
Stefan Knecht,
Emmanuel Fromager
Abstract:
The extrapolation technique of Savin [J. Chem. Phys. 140, 18A509 (2014)], which was initially applied to range-separated ground-state-density-functional Hamiltonians, is adapted in this work to ghost-interaction-corrected (GIC) range-separated ensemble density-functional theory (eDFT) for excited states. While standard extrapolations rely on energies that decay as $μ^{-2}$ in the large range-separ…
▽ More
The extrapolation technique of Savin [J. Chem. Phys. 140, 18A509 (2014)], which was initially applied to range-separated ground-state-density-functional Hamiltonians, is adapted in this work to ghost-interaction-corrected (GIC) range-separated ensemble density-functional theory (eDFT) for excited states. While standard extrapolations rely on energies that decay as $μ^{-2}$ in the large range-separation-parameter $μ$ limit, we show analytically that (approximate) range-separated GIC ensemble energies converge more rapidly (as $μ^{-3}$) towards their pure wavefunction theory values ($μ\rightarrow+\infty$ limit), thus requiring a different extrapolation correction. The purpose of such a correction is to further improve on the convergence and, consequently, to obtain more accurate excitation energies for a finite (and, in practice, relatively small) $μ$ value. As a proof of concept, we apply the extrapolation method to He and small molecular systems (viz. H$_{2}$, HeH$^{+}$ and LiH), thus considering different types of excitations like Rydberg, charge transfer and double excitations. Potential energy profiles of the first three and four singlet $Σ^+$ excitation energies in HeH$^{+}$ and H$_{2}$, respectively, are studied with a particular focus on avoided crossings for the latter. Finally, the extraction of individual state energies from the ensemble energy is discussed in the context of range-separated eDFT, as a perspective.
△ Less
Submitted 16 October, 2017; v1 submitted 11 August, 2017;
originally announced August 2017.
-
Exact ensemble density functional theory for excited states in a model system: investigating the weight dependence of the correlation energy
Authors:
Killian Deur,
Laurent Mazouin,
Emmanuel Fromager
Abstract:
Ensemble density functional theory (eDFT) is an exact time-independent alternative to time-dependent DFT (TD-DFT) for the calculation of excitation energies. Despite its formal simplicity and advantages in contrast to TD-DFT (multiple excitations, for example, can be easily taken into account in an ensemble), eDFT is not standard which is essentially due to the lack of reliable approximate exchang…
▽ More
Ensemble density functional theory (eDFT) is an exact time-independent alternative to time-dependent DFT (TD-DFT) for the calculation of excitation energies. Despite its formal simplicity and advantages in contrast to TD-DFT (multiple excitations, for example, can be easily taken into account in an ensemble), eDFT is not standard which is essentially due to the lack of reliable approximate exchange-correlation (xc) functionals for ensembles. Following Smith et al. [Phys. Rev. B 93, 245131 (2016)], we propose in this work to construct an exact eDFT for the nontrivial asymmetric Hubbard dimer, thus providing more insight into the weight dependence of the ensemble xc energy in various correlation regimes. For that purpose, an exact analytical expression for the weight-dependent ensemble exchange energy has been derived. The complementary exact ensemble correlation energy has been computed by means of Legendre--Fenchel transforms. Interesting features like discontinuities in the ensemble xc potential in the strongly correlated limit have been rationalized by means of a generalized adiabatic connection formalism. Finally, functional-driven errors induced by ground-state density-functional approximations have been studied. In the strictly symmetric case or in the weakly correlated regime, combining ensemble exact exchange with ground-state correlation functionals gives relatively accurate ensemble energies. However, when approaching the equiensemble in the strongly correlated regime, this approximation leads to highly curved ensemble energies with negative slope which is unphysical. Using both ground-state exchange and correlation functionals gives much better results in that case. In fact, exact ensemble energies are almost recovered in some density domains. The analysis of density-driven errors is left for future work.
△ Less
Submitted 16 December, 2016; v1 submitted 11 October, 2016;
originally announced October 2016.
-
Ghost interaction correction in ensemble density-functional theory for excited states with and without range separation
Authors:
Md. Mehboob Alam,
Stefan Knecht,
Emmanuel Fromager
Abstract:
Ensemble density-functional theory (eDFT) suffers from the so-called "ghost interaction" error when approximate exchange-correlation functionals are used. In this work, we present a rigorous ghost interaction correction (GIC) scheme in the context of range-separated eDFT. The method relies on an exact decomposition of the ensemble short-range exchange-correlation energy into a multideterminantal e…
▽ More
Ensemble density-functional theory (eDFT) suffers from the so-called "ghost interaction" error when approximate exchange-correlation functionals are used. In this work, we present a rigorous ghost interaction correction (GIC) scheme in the context of range-separated eDFT. The method relies on an exact decomposition of the ensemble short-range exchange-correlation energy into a multideterminantal exact exchange term, which involves the long-range interacting ensemble density matrix instead of the Kohn--Sham (KS) one, and a complementary density-functional correlation energy. A generalized adiabatic connection formula is derived for the latter. In order to perform practical calculations, the complementary correlation functional has been simply modeled by its ground-state local density approximation (LDA) while long-range interacting ground- and excited-state wavefunctions have been obtained self-consistently by combining a long-range configuration interaction calculation with a short-range LDA potential. We show that GIC reduces the curvature of approximate range-separated ensemble energies drastically while providing considerably more accurate excitation energies, even for charge-transfer and double excitations. Interestingly, the method performs well also in the context of standard KS-eDFT, which is recovered when the range-separation parameter is set to zero.
△ Less
Submitted 3 June, 2016; v1 submitted 17 March, 2016;
originally announced March 2016.
-
Local density approximation in site-occupation embedding theory
Authors:
Bruno Senjean,
Masahisa Tsuchiizu,
Vincent Robert,
Emmanuel Fromager
Abstract:
Site-occupation embedding theory (SOET) is a density-functional theory (DFT)-based method which aims at modelling strongly correlated electrons. It is in principle exact and applicable to model and quantum chemical Hamiltonians. The theory is presented here for the Hubbard Hamiltonian. In contrast to conventional DFT approaches, the site (or orbital) occupations are deduced in SOET from a partiall…
▽ More
Site-occupation embedding theory (SOET) is a density-functional theory (DFT)-based method which aims at modelling strongly correlated electrons. It is in principle exact and applicable to model and quantum chemical Hamiltonians. The theory is presented here for the Hubbard Hamiltonian. In contrast to conventional DFT approaches, the site (or orbital) occupations are deduced in SOET from a partially-interacting system consisting of one (or more) impurity site(s) and non-interacting bath sites. The correlation energy of the bath is then treated implicitly by means of a site-occupation functional. In this work, we propose a simple impurity-occupation functional approximation based on the two-level (2L) Hubbard model which is referred to as two-level impurity local density approximation (2L-ILDA). Results obtained on a prototypical uniform 8-site Hubbard ring are promising. The extension of the method to larger systems and more sophisticated model Hamiltonians is currently in progress.
△ Less
Submitted 11 April, 2016; v1 submitted 8 February, 2016;
originally announced February 2016.
-
Combining linear interpolation with extrapolation methods in range-separated ensemble density-functional theory
Authors:
Bruno Senjean,
Erik D. Hedegård,
Md. Mehboob Alam,
Stefan Knecht,
Emmanuel Fromager
Abstract:
The combination of a recently proposed linear interpolation method (LIM) [Senjean et al., Phys. Rev. A 92, 012518 (2015)], which enables the calculation of weight-independent excitation energies in range-separated ensemble density-functional approximations, with the extrapolation scheme of Savin [J. Chem. Phys. 140, 18A509 (2014)] is presented in this work. It is shown that LIM excitation energies…
▽ More
The combination of a recently proposed linear interpolation method (LIM) [Senjean et al., Phys. Rev. A 92, 012518 (2015)], which enables the calculation of weight-independent excitation energies in range-separated ensemble density-functional approximations, with the extrapolation scheme of Savin [J. Chem. Phys. 140, 18A509 (2014)] is presented in this work. It is shown that LIM excitation energies vary quadratically with the inverse of the range-separation parameter mu when the latter is large. As a result, the extrapolation scheme, which is usually applied to long-range interacting energies, can be adapted straightforwardly to LIM. This extrapolated LIM (ELIM) has been tested on a small test set consisting of He, Be, H2 and HeH+. Relatively accurate results have been obtained for the first singlet excitation energies with the typical mu=0.4 value. The improvement of LIM after extrapolation is remarkable, in particular for the doubly-excited 2^1Sigma+g state in the stretched H2 molecule. Three-state ensemble calculations in H2 also show that ELIM does not necessarily improves relative excitation energies, even though individual excitation energies are more accurate after extrapolation. Finally, an alternative decomposition of the short-range ensemble exchange-correlation energy is proposed in order to correct for ghost-interaction errors in multideterminant range-separated ensemble density-functional theory calculations. The implementation and calibration of such a scheme is currently in progress.
△ Less
Submitted 6 November, 2015; v1 submitted 26 September, 2015;
originally announced September 2015.
-
Linear interpolation method in ensemble Kohn-Sham and range-separated density-functional approximations for excited states
Authors:
Bruno Senjean,
Stefan Knecht,
Hans Jørgen Aa. Jensen,
Emmanuel Fromager
Abstract:
Gross-Oliveira-Kohn density functional theory (GOK-DFT) for ensembles is in principle very attractive, but has been hard to use in practice. A novel, practical model based on GOK-DFT for the calculation of electronic excitation energies is discussed. The new model relies on two modifications of GOK-DFT: use of range separation and use of the slope of the linearly-interpolated ensemble energy, rath…
▽ More
Gross-Oliveira-Kohn density functional theory (GOK-DFT) for ensembles is in principle very attractive, but has been hard to use in practice. A novel, practical model based on GOK-DFT for the calculation of electronic excitation energies is discussed. The new model relies on two modifications of GOK-DFT: use of range separation and use of the slope of the linearly-interpolated ensemble energy, rather than orbital energies. The range-separated approach is appealing as it enables the rigorous formulation of a multi-determinant state-averaged DFT method. In the exact theory, the short-range density functional, that complements the long-range wavefunction-based ensemble energy contribution, should vary with the ensemble weights even when the density is held fixed. This weight dependence ensures that the range-separated ensemble energy varies linearly with the ensemble weights. When the (weight-independent) ground-state short-range exchange-correlation functional is used in this context, curvature appears thus leading to an approximate weight-dependent excitation energy. In order to obtain unambiguous approximate excitation energies, we propose to interpolate linearly the ensemble energy between equiensembles. It is shown that such a linear interpolation method (LIM) can be rationalized and that it effectively introduces weight dependence effects. As proof of principle, LIM has been applied to He, Be, H$_2$ in both equilibrium and stretched geometries as well as the stretched HeH$^+$ molecule. Very promising results have been obtained for both single (including charge transfer) and double excitations with spin-independent short-range local and semi-local functionals. Even at the Kohn--Sham ensemble DFT level, that is recovered when the range-separation parameter is set to zero, LIM performs better than standard time-dependent DFT.
△ Less
Submitted 12 June, 2015; v1 submitted 24 April, 2015;
originally announced April 2015.
-
On the exact formulation of multi-configuration density-functional theory: electron density versus orbitals occupation
Authors:
Emmanuel Fromager
Abstract:
The exact formulation of multi-configuration density-functional theory (DFT) is discussed in this work. As an alternative to range-separated methods, where electron correlation effects are split in the coordinate space, the combination of Configuration Interaction methods with orbital occupation functionals is explored at the formal level through the separation of correlation effects in the orbita…
▽ More
The exact formulation of multi-configuration density-functional theory (DFT) is discussed in this work. As an alternative to range-separated methods, where electron correlation effects are split in the coordinate space, the combination of Configuration Interaction methods with orbital occupation functionals is explored at the formal level through the separation of correlation effects in the orbital space. When applied to model Hamiltonians, this approach leads to an exact Site-Occupation Embedding Theory (SOET). An adiabatic connection expression is derived for the complementary bath functional and a comparison with Density Matrix Embedding Theory (DMET) is made. Illustrative results are given for the simple two-site Hubbard model. SOET is then applied to a quantum chemical Hamiltonian, thus leading to an exact Complete Active Space Site-Occupation Functional Theory (CASSOFT) where active electrons are correlated explicitly within the CAS and the remaining contributions to the correlation energy are described with an orbital occupation functional. The computational implementation of SOET and CASSOFT as well as the development of approximate functionals are left for future work.
△ Less
Submitted 25 October, 2014; v1 submitted 8 September, 2014;
originally announced September 2014.
-
Double hybrid density-functional theory using the Coulomb-attenuating method
Authors:
Yann Cornaton,
Emmanuel Fromager
Abstract:
A double hybrid approximation using the Coulomb-attenuating method (CAM-DH) is derived within range-separated density-functional perturbation theory, in the spirit of a recent work by Cornaton {\it et al.} [Phys. Rev. A 88, 022516 (2013)]. The energy expression recovered through second order is linear in the parameters $α$ and $β$ that control the Coulomb attenuation. The method has been tested wi…
▽ More
A double hybrid approximation using the Coulomb-attenuating method (CAM-DH) is derived within range-separated density-functional perturbation theory, in the spirit of a recent work by Cornaton {\it et al.} [Phys. Rev. A 88, 022516 (2013)]. The energy expression recovered through second order is linear in the parameters $α$ and $β$ that control the Coulomb attenuation. The method has been tested within the local density approximation on a small test set consisting of rare-gas and alkaline-earth-metal dimers as well as diatomics with single, double and triple bonds. In this context, the semi-empirical $α=0.19$ and $β=0.46$ parameters, that were optimized for the hybrid CAM-B3LYP functional, do not provide accurate interaction and total energies. Using semi-local functionals with density scaling, that was neglected in this work, may lead to different conclusions. Calibration studies on a larger test set would be necessary at this point. This is left for future work. Finally, we propose as a perspective an alternative CAM-DH approach that relies on the perturbation expansion of a partially long-range interacting wavefunction. In this case the energy is not linear anymore in $α$ and $β$. Work is in progress in this direction.
△ Less
Submitted 18 March, 2014; v1 submitted 2 December, 2013;
originally announced December 2013.
-
Assessment of Charge-Transfer Excitations in Organic Dyes obtained from TD-srDFT Based on Long-Range MP2 and MCSCF Wave Functions
Authors:
Erik D. Hedegård,
Frank Heiden,
Stefan Knecht,
Emmanuel Fromager,
Hans Jørgen Aa. Jensen
Abstract:
Charge transfer excitations can be described within TD-DFT, not only by means of long-range corrected exchange functionals but also with a combination of wave function theory and TD-DFT based on range separation. The latter approach enables a rigorous formulation of multi-determinantal TD-DFT schemes where excitation classes, which are absent in conventional TD-DFT spectra (like for example double…
▽ More
Charge transfer excitations can be described within TD-DFT, not only by means of long-range corrected exchange functionals but also with a combination of wave function theory and TD-DFT based on range separation. The latter approach enables a rigorous formulation of multi-determinantal TD-DFT schemes where excitation classes, which are absent in conventional TD-DFT spectra (like for example double excitations), can be addressed. This paper investigates the combination of both the long-range MCSCF and SOPPA ansätze with a short-range DFT (srDFT) description. We find that the combinations of SOPPA or MCSCF with TD-DFT yield better results than could be expected from the pure wave function schemes. For the Time-Dependent MCSCF short-range DFT ansatz (TD-MC-srDFT) excitation energies calculated over a larger benchmark set of molecules with predominantly single reference character yield good agreement with their reference values, and are in general comparable to the long-range corrected functional CAM-B3LYP. The SOPPA-srDFT scheme is tested for a subset of molecules used for benchmarking TD-MC-srDFT and performs slightly better against the reference data for this small subset. Beyond the proof-of-principle calculations comprising the first part of this contribution, we additionally studied the low-lying singlet excited states (S1 and S2) of the retinal chromophore. The chromophore displays multireference character in the ground state and both excited states exhibit considerable double excitation character, which in turn cannot be described within standard TD-DFT, due to the adiabatic approximation. However, a TD-MC-srDFT approach can account for the multireference character, and excitation energies are obtained with accuracy comparable to CASPT2, although using a much smaller active space.
△ Less
Submitted 22 August, 2013;
originally announced August 2013.
-
Generalized adiabatic connection in ensemble density-functional theory for excited states: example of the H2 molecule
Authors:
Odile Franck,
Emmanuel Fromager
Abstract:
A generalized adiabatic connection for ensembles (GACE) is presented. In contrast to the traditional adiabatic connection formulation, both ensemble weights and interaction strength can vary along a GACE path while the ensemble density is held fixed. The theory is presented for non-degenerate two-state ensembles but it can in principle be extended to any ensemble of fractionally occupied excited s…
▽ More
A generalized adiabatic connection for ensembles (GACE) is presented. In contrast to the traditional adiabatic connection formulation, both ensemble weights and interaction strength can vary along a GACE path while the ensemble density is held fixed. The theory is presented for non-degenerate two-state ensembles but it can in principle be extended to any ensemble of fractionally occupied excited states. Within such a formalism an exact expression for the ensemble exchange-correlation density-functional energy, in terms of the conventional ground-state exchange-correlation energy, is obtained by integration over the ensemble weight. Stringent constraints on the functional are thus obtained when expanding the ensemble exchange-correlation energy through second order in the ensemble weight. For illustration purposes, the analytical derivation of the GACE is presented for the H2 model system in a minimal basis, leading thus to a simple density-functional approximation to the ensemble exchange-correlation energy. Encouraging results were obtained with this approximation for the description in a large basis of the first ^1Σ^+_g excitation in H2 upon bond stretching. Finally, a range-dependent GACE has been derived, providing thus a pathway to the development of a rigorous state-average multi-determinant density-functional theory.
△ Less
Submitted 16 October, 2013; v1 submitted 21 August, 2013;
originally announced August 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.
-
Alternative separation of exchange and correlation energies in range-separated density-functional perturbation theory
Authors:
Yann Cornaton,
Alexandrina Stoyanova,
Hans Jørgen Aa. Jensen,
Emmanuel Fromager
Abstract:
An alternative separation of short-range exchange and correlation energies is used in the framework of second-order range-separated density-functional perturbation theory. This alternative separation was initially proposed by Toulouse et al. [Theor. Chem. Acc. 114, 305 (2005)] and relies on a long-range interacting wavefunction instead of the non-interacting Kohn-Sham one. When second-order correc…
▽ More
An alternative separation of short-range exchange and correlation energies is used in the framework of second-order range-separated density-functional perturbation theory. This alternative separation was initially proposed by Toulouse et al. [Theor. Chem. Acc. 114, 305 (2005)] and relies on a long-range interacting wavefunction instead of the non-interacting Kohn-Sham one. When second-order corrections to the density are neglected, the energy expression reduces to a range-separated double-hybrid (RSDH) type of functional, RSDHf, where "f" stands for "full-range integrals" as the regular full-range interaction appears explicitly in the energy expression when expanded in perturbation theory. In contrast to usual RSDH functionals, RSDHf describes the coupling between long- and short-range correlations as an orbital-dependent contribution. Calculations on the first four noble-gas dimers show that this coupling has a significant effect on the potential energy curves in the equilibrium region, improving the accuracy of binding energies and equilibrium bond distances when second-order perturbation theory is appropriate.
△ Less
Submitted 8 August, 2013; v1 submitted 25 June, 2013;
originally announced June 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.
-
Multi-configuration time-dependent density-functional theory based on range separation
Authors:
Emmanuel Fromager,
Stefan Knecht,
Hans Jørgen Aa. Jensen
Abstract:
Multi-configuration range-separated density-functional theory is extended to the time-dependent regime. An exact variational formulation is derived. The approximation, which consists in combining a long-range Multi-Configuration-Self-Consistent Field (MCSCF) treatment with an adiabatic short-range density-functional (DFT) description, is then considered. The resulting time-dependent multi-configur…
▽ More
Multi-configuration range-separated density-functional theory is extended to the time-dependent regime. An exact variational formulation is derived. The approximation, which consists in combining a long-range Multi-Configuration-Self-Consistent Field (MCSCF) treatment with an adiabatic short-range density-functional (DFT) description, is then considered. The resulting time-dependent multi-configuration short-range DFT (TD-MC-srDFT) model is applied to the calculation of singlet excitation energies in H2, Be and ferrocene, considering both short-range local density (srLDA) and generalized gradient (srGGA) approximations. In contrast to regular TD-DFT, TD-MC-srDFT can describe double excitations. As expected, when modeling long-range interactions with the MCSCF model instead of the adiabatic Buijse-Baerends density-matrix functional as recently proposed by Pernal [K. Pernal, J. Chem. Phys. 136, 184105 (2012)], the description of both the 1^1D doubly-excited state in Be and the 1^1Σ^+_u state in the stretched H2 molecule are improved, although the latter is still significantly underestimated. Exploratory TD-MC-srDFT/GGA calculations for ferrocene yield in general excitation energies at least as good as TD-DFT/CAM-B3LYP, and superior to wave-function (TD-MCSCF, symmetry adapted cluster-configuration interaction) and TD-DFT results based on LDA, GGA, and hybrid functionals.
△ Less
Submitted 14 January, 2013; v1 submitted 20 November, 2012;
originally announced November 2012.
-
Rigorous formulation of two-parameter double-hybrid density-functionals
Authors:
Emmanuel Fromager
Abstract:
A two-parameter extension of the density-scaled double hybrid approach of Sharkas et al. [J. Chem. Phys. 134, 064113 (2011)] is presented. It is based on the explicit treatment of a fraction of multideterminantal exact exchange. The connection with conventional double hybrids is made when neglecting density scaling in the correlation functional as well as second-order corrections to the density. I…
▽ More
A two-parameter extension of the density-scaled double hybrid approach of Sharkas et al. [J. Chem. Phys. 134, 064113 (2011)] is presented. It is based on the explicit treatment of a fraction of multideterminantal exact exchange. The connection with conventional double hybrids is made when neglecting density scaling in the correlation functional as well as second-order corrections to the density. In this context, the fraction ac of second-order Møller-Plesset (MP2) correlation energy is not necessarily equal to the square of the fraction ax of Hartree-Fock exchange. More specifically, it is shown that ac \leq ax2, a condition that conventional semi-empirical double hybrids actually fulfill. In addition, a new procedure for calculating the orbitals, which has a better justification than the one routinely used, is proposed. Referred to as λ1 variant, the corresponding double hybrid approximation has been tested on a small set consisting of H2, N2, Be2, Mg2, and Ar2. Three conventional double hybrids (B2-PLYP, B2GP-PLYP, and PBE0-DH) have been considered. Potential curves obtained with λ1- and regular double hybrids can, in some cases, differ significantly. In particular, for the weakly bound dimers, the λ1 variants bind systematically more than the regular ones, which is an improvement in many but not all cases. Including density scaling in the correlation functionals may of course change the results significantly. Moreover, optimized effective potentials based on a partially-interacting system could also be used to generate proper orbitals. Work is currently in progress in those directions.
△ Less
Submitted 11 January, 2012; v1 submitted 7 September, 2011;
originally announced September 2011.
