Source author record

Emmanuel Fromager

Emmanuel Fromager appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2022arXiv

Local Potential Functional Embedding Theory: A Self-Consistent Flavor of Density Functional Theory for Lattices without Density Functionals

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.

preprint2022arXiv

Quantum embedding of multi-orbital fragments using the Block-Householder-transformation

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).

preprint2021arXiv

Exact exchange-correlation potentials for calculating the fundamental gap with a fixed number of electrons

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.

preprint2020arXiv

A weight-dependent local correlation density-functional approximation for ensembles

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.

preprint2020arXiv

Individual correlations in ensemble density-functional theory: State-driven/density-driven decompositions without additional Kohn-Sham systems

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.

preprint2020arXiv

N-centered ensemble density-functional theory for open systems

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

preprint2020arXiv

Weight Dependence of Local Exchange-Correlation Functionals in Ensemble Density-Functional Theory: Double Excitations in Two-Electron Systems

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.

preprint2016arXiv

Ghost interaction correction in ensemble density-functional theory for excited states with and without range separation

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.

preprint2015arXiv

Linear interpolation method in ensemble Kohn-Sham and range-separated density-functional approximations for excited states

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.

preprint2014arXiv

Double hybrid density-functional theory using the Coulomb-attenuating method

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.

preprint2014arXiv

On the exact formulation of multi-configuration density-functional theory: electron density versus orbitals occupation

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.

preprint2013arXiv

Alternative separation of exchange and correlation energies in multi-configuration range-separated density-functional theory

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.

preprint2013arXiv

Alternative separation of exchange and correlation energies in range-separated density-functional perturbation theory

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.

preprint2013arXiv

Analysis of double hybrid density-functionals along the adiabatic connection

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.

preprint2013arXiv

Assessment of Charge-Transfer Excitations in Organic Dyes obtained from TD-srDFT Based on Long-Range MP2 and MCSCF Wave Functions

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.

preprint2013arXiv

Generalized adiabatic connection in ensemble density-functional theory for excited states: example of the H2 molecule

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.

preprint2013arXiv

Multi-configuration time-dependent density-functional theory based on range separation

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.

preprint2012arXiv

Rigorous formulation of two-parameter double-hybrid density-functionals

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.