Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Computation

Levy-Lieb-Based Monte Carlo Study of the Dimensionality Behaviour of the Electronic Kinetic Functional

A. Seshaditya, Luca M. Ghiringhelli, Luigi Delle Site

Published Jun 10, 2017• 1 citation
Publisher

Abstract

We consider a gas of interacting electrons in the limit of nearly uniform density and treat the one dimensional (1D), two dimensional (2D) and three dimensional (3D) cases. We focus on the determination of the correlation part of the kinetic functional by employing a Monte Carlo sampling technique of electrons in space based on an analytic derivation via the Levy-Lieb constrained search principle. Of particular interest is the question of the behaviour of the functional as one passes from 1D to 3D; according to the basic principles of Density Functional Theory (DFT) the form of the universal functional should be independent of the dimensionality. However, in practice the straightforward use of current approximate functionals in different dimensions is problematic. Here, we show that going from the 3D to the 2D case the functional form is consistent (concave function) but in 1D becomes convex; such a drastic difference is peculiar of 1D electron systems as it is for other quantities. Given the interesting behaviour of the functional, this study represents a basic first-principle approach to the problem and suggests further investigations using highly accurate (though expensive) many-electron computational techniques, such as Quantum Monte Carlo.

DOI
10.3390/computation5020030

Related papers

Electronic Energy Functionals: Levy-Lieb principle within the Ground State Path Integral Quantum Monte Carlo

Luigi Delle Site, Luca M. Ghiringhelli, David M. Ceperley

We propose a theoretical/computational protocol based on the use of the Ground State (GS) Path Integral (PI) Quantum Monte Carlo (QMC) for the calculation of the kinetic and Coulomb energy density for a system of $N$ interacting electrons in an external potential. The idea is based on the derivation of the energy densities via the $N-1$-conditional probability density within the framework of the Levy-Lieb constrained search principle. The consequences for the development of energy functionals within the context of Density Functional Theory (DFT) are discussed. We propose also the possibility of going beyond the energy densities and extend this idea to a computational procedure where the $N-1$-conditional probability is an implicit functional of the electron density, independently from the external potential. In principle, such a procedure paves the way for an {\it on-the-fly} determination of the energy functional for any system.

arXiv:1205.6503•May 29, 2012•0 citations

Electronic energy functionals: Levy–Lieb principle within the ground state path integral quantum Monte Carlo

Luigi Delle Site, Luca M. Ghiringhelli, David M. Ceperley

Abstract We propose a theoretical/computational protocol based on the use of the Ground State Path Integral Quantum Monte Carlo for the calculation of the kinetic and Coulomb energy density for a system of N interacting electrons in an external potential. The idea is based on the derivation of the energy densities via the (N − 1)‐conditional probability density within the framework of the Levy–Lieb constrained search principle. The consequences for the development of energy functionals within the context of density functional theory are discussed. We propose also the possibility of going beyond the energy densities and extend this idea to a computational procedure where the (N − 1)‐conditional probability is an implicit functional of the electron density, independently from the external potential. In principle, such a procedure paves the way for an on‐the‐fly determination of the energy functional for any system. © 2012 Wiley Periodicals, Inc.

International Journal of Quantum Chemistry•Aug 30, 2012•10 citations

Interacting electrons, spin statistics, and information theory

Luca M. Ghiringhelli, Ian P. Hamilton, Luigi Delle Site

We consider a nearly (or quasi) uniform gas of interacting electrons for which spin statistics play a crucial role. A previously developed procedure, based on the extension of the Levy-Lieb constrained search principle and Monte Carlo sampling of electron configurations in space, allows us to approximate the form of the kinetic-energy functional. For a spinless electron gas, this procedure led to a correlation term, which had the form of the Shannon entropy, but the resulting kinetic-energy functional does not satisfy the Lieb-Thirring inequality, which is rigorous and one of the most general relations regarding the kinetic energy. In this paper, we show that when the fermionic character of the electrons is included via a statistical spin approach, our procedure leads to correlation terms, which also have the form of the Shannon entropy and the resulting kinetic-energy functional does satisfy the Lieb-Thirring inequality. In this way we further strengthen the connection between Shannon entropy and electron correlation and, more generally, between information theory and quantum mechanics.

The Journal of Chemical Physics•Jan 6, 2010•55 citations

Design of kinetic functionals for many-body electron systems: Combining analytical theory with Monte Carlo sampling of electronic configurations

Luca M. Ghiringhelli, Luigi Delle Site

In a previous work [L. Delle Site, J. Phys. A: Math. Theor. 40, 2787 (2007)], the derivation of an analytic expression for the kinetic functional of a many-body electron system has been proposed. Though analytical, the formula is still nonlocal (multidimensional) and thus not ideal for numerical applications. In this work, by treating the test case of a uniform gas of interacting spinless electrons, we propose a computational protocol which combines the previous analytic results with the Monte Carlo sampling of electronic configurations in space. This, we show, leads to an internally consistent scheme to design well founded local kinetic functionals.

arXiv:0711.1832•Physical Review B•Feb 28, 2008•40 citations

New methods for the ab-initio simulation of correlated systems

Schade, Robert

Strong electronic correlations are at the heart of many interesting phenomena. For the theoretical description of these materials, a proper treatment of the local atomic physics is required. We propose a novel approach combining functionals of the electron density and functionals of the one-particle reduced density matrix to improve this description. The proposed method has a solid foundation in reduced density-matrix functional theory and no double-counting problem arises. It employs a decomposition of the electron-electron interaction in real space. The interaction close to the correlated orbitals, for example, the partially filled 3d-orbitals of transition-metal ions, is described with a density-matrix functional and otherwise with a local or semi-local density functional. We propose to evaluate the density-matrix functional from Levy’s constrained search problem, i.e., via a constrained minimization over an ensemble of many-particle wave functions. In contrast to approximate parametrized density-matrix functionals, this evaluation allows us to systematically improve the functional towards the exact result. In situations where the one-particle basis is too large to evaluate the density-matrix functional from Levy’s constrained minimization problem, we apply a series of approximations that each can be converged to the exact result: the first approximation step is the local approximation of the density-matrix functional proposed by Blöchl, Walther and Pruschke. For the density-matrix functionals within the local interactions, we propose the adaptive cluster approximation (ACA) that systematically truncates non-interacting one-particle states and drastically reduces the computational effort. The resulting density-matrix functional for a local interaction and a small number of non-interacting one-particle states is then evaluated with the constrained minimization problem. We explore different parametrizations of many-particle wave functions in this constrained minimization problem. A parametrization based on a configuration-interaction-like ansatz is shown to converge rapidly if suitable selection criterion for the Slater determinants is chosen. An impurity-bath-separation ansatz is shown to be suitable for single-impurity Anderson models. It is shown that the constrained minimization can be solved for matrix product states with a DMRG-like iterative minimization. Furthermore, we show that Gutzwiller-Jastrow-correlated wave functions can be used with a quantum Monte Carlo procedure as many-particle wave functions. Finally, we formulate an algorithm for the evaluation of the density-matrix functional on near-term quantum computers. Results from the execution of the algorithm on an existing quantum computer with transmon qubits are presented. The proposed approach combining density functionals and density-matrix functionals is implemented in the CP-PAW code based on the projector augmented-wave formalism. We present results for the dissociation curve of the hydrogen molecule as the prototypical case of strong static correlation and show that static correlation is well described with the new approach. We show results for the nonmagnetic state of the transition-metaoxide NiO that is described qualitatively wrong with the DFT+U method. The proposed method properly describes the nonmagnetic state of NiO as an insulator and predicts a qualitatively correct spectral function.

•Jan 1, 2019•1 citation

Density functional theory beyond the linear regime: Validating an adiabatic local density approximation

N. Helbig, Johanna I. Fuks, Michele Casula +4 more

We present a local density approximation (LDA) for one-dimensional (1D) systems interacting via the soft-Coulomb interaction based on quantum Monte Carlo calculations. Results for the ground-state energies and ionization potentials of finite 1D systems show excellent agreement with exact calculations obtained by exploiting the mapping of an $N$-electron system in $d$ dimensions onto a single electron in $N\ifmmode\times\else\texttimes\fi{}d$ dimensions, properly symmetrized by the Young diagrams. We conclude that 1D LDA is of the same quality as its three-dimensional (3D) counterpart, and we infer conclusions about 3D LDA. The linear and nonlinear time-dependent responses of 1D model systems using LDA, exact exchange, and the exact solution are investigated and show very good agreement in both cases, except for the well-known problem of missing double excitations. Consequently, the 3D LDA is expected to be of good quality beyond the linear response. In addition, the 1D LDA should prove useful in modeling the interaction of atoms with strong laser fields, where this specific 1D model is often used.

arXiv:1101.2564•Physical Review A•Mar 14, 2011•83 citations
We consider a gas of interacting electrons in the limit of nearly uniform density and treat the one dimensional (1D), two dimensional (2D) and three dimensional (3D) cases. We focus on the determination of the correlation part of the kinetic functional by employing a Monte Carlo sampling technique of electrons in space based on an analytic derivation via the Levy-Lieb constrained search principle. Of particular interest is the question of the behaviour of the functional as one passes from 1D to 3D; according to the basic principles of Density Functional Theory (DFT) the form of the universal functional should be independent of the dimensionality. However, in practice the straightforward use of current approximate functionals in different dimensions is problematic. Here, we show that going from the 3D to the 2D case the functional form is consistent (concave function) but in 1D becomes convex; such a drastic difference is peculiar of 1D electron systems as it is for other quantities. Given the interesting behaviour of the functional, this study represents a basic first-principle approach to the problem and suggests further investigations using highly accurate (though expensive) many-electron computational techniques, such as Quantum Monte Carlo.
We propose a theoretical/computational protocol based on the use of the Ground State (GS) Path Integral (PI) Quantum Monte Carlo (QMC) for the calculation of the kinetic and Coulomb energy density for a system of NNN interacting electrons in an external potential. The idea is based on the derivation of the energy densities via the N−1N-1N−1-conditional probability density within the framework of the Levy-Lieb constrained search principle. The consequences for the development of energy functionals within the context of Density Functional Theory (DFT) are discussed. We propose also the possibility of going beyond the energy densities and extend this idea to a computational procedure where the N−1N-1N−1-conditional probability is an implicit functional of the electron density, independently from the external potential. In principle, such a procedure paves the way for an {\it on-the-fly} determination of the energy functional for any system.
Abstract We propose a theoretical/computational protocol based on the use of the Ground State Path Integral Quantum Monte Carlo for the calculation of the kinetic and Coulomb energy density for a system of N interacting electrons in an external potential. The idea is based on the derivation of the energy densities via the (N − 1)‐conditional probability density within the framework of the Levy–Lieb constrained search principle. The consequences for the development of energy functionals within the context of density functional theory are discussed. We propose also the possibility of going beyond the energy densities and extend this idea to a computational procedure where the (N − 1)‐conditional probability is an implicit functional of the electron density, independently from the external potential. In principle, such a procedure paves the way for an on‐the‐fly determination of the energy functional for any system. © 2012 Wiley Periodicals, Inc.
We consider a nearly (or quasi) uniform gas of interacting electrons for which spin statistics play a crucial role. A previously developed procedure, based on the extension of the Levy-Lieb constrained search principle and Monte Carlo sampling of electron configurations in space, allows us to approximate the form of the kinetic-energy functional. For a spinless electron gas, this procedure led to a correlation term, which had the form of the Shannon entropy, but the resulting kinetic-energy functional does not satisfy the Lieb-Thirring inequality, which is rigorous and one of the most general relations regarding the kinetic energy. In this paper, we show that when the fermionic character of the electrons is included via a statistical spin approach, our procedure leads to correlation terms, which also have the form of the Shannon entropy and the resulting kinetic-energy functional does satisfy the Lieb-Thirring inequality. In this way we further strengthen the connection between Shannon entropy and electron correlation and, more generally, between information theory and quantum mechanics.
In a previous work [L. Delle Site, J. Phys. A: Math. Theor. 40, 2787 (2007)], the derivation of an analytic expression for the kinetic functional of a many-body electron system has been proposed. Though analytical, the formula is still nonlocal (multidimensional) and thus not ideal for numerical applications. In this work, by treating the test case of a uniform gas of interacting spinless electrons, we propose a computational protocol which combines the previous analytic results with the Monte Carlo sampling of electronic configurations in space. This, we show, leads to an internally consistent scheme to design well founded local kinetic functionals.
Strong electronic correlations are at the heart of many interesting phenomena. For the theoretical description of these materials, a proper treatment of the local atomic physics is required. We propose a novel approach combining functionals of the electron density and functionals of the one-particle reduced density matrix to improve this description. The proposed method has a solid foundation in reduced density-matrix functional theory and no double-counting problem arises. It employs a decomposition of the electron-electron interaction in real space. The interaction close to the correlated orbitals, for example, the partially filled 3d-orbitals of transition-metal ions, is described with a density-matrix functional and otherwise with a local or semi-local density functional. We propose to evaluate the density-matrix functional from Levy’s constrained search problem, i.e., via a constrained minimization over an ensemble of many-particle wave functions. In contrast to approximate parametrized density-matrix functionals, this evaluation allows us to systematically improve the functional towards the exact result. In situations where the one-particle basis is too large to evaluate the density-matrix functional from Levy’s constrained minimization problem, we apply a series of approximations that each can be converged to the exact result: the first approximation step is the local approximation of the density-matrix functional proposed by Blöchl, Walther and Pruschke. For the density-matrix functionals within the local interactions, we propose the adaptive cluster approximation (ACA) that systematically truncates non-interacting one-particle states and drastically reduces the computational effort. The resulting density-matrix functional for a local interaction and a small number of non-interacting one-particle states is then evaluated with the constrained minimization problem. We explore different parametrizations of many-particle wave functions in this constrained minimization problem. A parametrization based on a configuration-interaction-like ansatz is shown to converge rapidly if suitable selection criterion for the Slater determinants is chosen. An impurity-bath-separation ansatz is shown to be suitable for single-impurity Anderson models. It is shown that the constrained minimization can be solved for matrix product states with a DMRG-like iterative minimization. Furthermore, we show that Gutzwiller-Jastrow-correlated wave functions can be used with a quantum Monte Carlo procedure as many-particle wave functions. Finally, we formulate an algorithm for the evaluation of the density-matrix functional on near-term quantum computers. Results from the execution of the algorithm on an existing quantum computer with transmon qubits are presented. The proposed approach combining density functionals and density-matrix functionals is implemented in the CP-PAW code based on the projector augmented-wave formalism. We present results for the dissociation curve of the hydrogen molecule as the prototypical case of strong static correlation and show that static correlation is well described with the new approach. We show results for the nonmagnetic state of the transition-metaoxide NiO that is described qualitatively wrong with the DFT+U method. The proposed method properly describes the nonmagnetic state of NiO as an insulator and predicts a qualitatively correct spectral function.
We present a local density approximation (LDA) for one-dimensional (1D) systems interacting via the soft-Coulomb interaction based on quantum Monte Carlo calculations. Results for the ground-state energies and ionization potentials of finite 1D systems show excellent agreement with exact calculations obtained by exploiting the mapping of an NNN-electron system in ddd dimensions onto a single electron in $N\ifmmode\times\else\texttimes\fi{}d$ dimensions, properly symmetrized by the Young diagrams. We conclude that 1D LDA is of the same quality as its three-dimensional (3D) counterpart, and we infer conclusions about 3D LDA. The linear and nonlinear time-dependent responses of 1D model systems using LDA, exact exchange, and the exact solution are investigated and show very good agreement in both cases, except for the well-known problem of missing double excitations. Consequently, the 3D LDA is expected to be of good quality beyond the linear response. In addition, the 1D LDA should prove useful in modeling the interaction of atoms with strong laser fields, where this specific 1D model is often used.