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

arXiv:2207.08995 · Physical Review A

Levy-Lieb embedding of density-functional theory and its quantum kernel: Illustration for the Hubbard dimer using near-term quantum algorithms

C. D. Pemmaraju, Amol A. Deshmukh

Published Oct 10, 2022• 2 citations
arXivHTMLPDF

Abstract

The constrained-search formulation of Levy and Lieb provides a concrete mapping from $N$-representable densities to the space of $N$-particle wave functions and explicitly defines the universal functional of density-functional theory. We numerically implement the Levy-Lieb procedure for a paradigmatic lattice system, the Hubbard dimer, using a modified variational quantum eigensolver approach. We demonstrate density variational minimization using the resulting hybrid quantum-classical scheme featuring real-time computation of the Levy-Lieb functional along the search trajectory. We further illustrate a fidelity-based quantum kernel associated with the density to pure-state embedding implied by the Levy-Lieb procedure and employ the kernel for learning observable functionals of the density. We study the kernel's ability to generalize with high accuracy through numerical experiments on the Hubbard dimer.

DOI
10.1103/physreva.106.042807

Related papers

N -representability and stationarity in time-dependent density-functional theory

Morrel H. Cohen, Adam Wasserman

To construct an $N$-representable time-dependent density-functional theory, a generalization to the time domain of the Levy-Lieb (LL) constrained-search algorithm is required. That the action is only stationary in the Dirac-Frenkel variational principle eliminates the possibility of basing the search on the action itself. Instead, we use the norm of the partial functional derivative of the action in the Hilbert space of the wave functions in place of the energy of the LL search. The electron densities entering the formalism are $N$-representable, and the resulting universal action functional has a unique stationary point in the density at that corresponding to the solution of the Schr\"odinger equation. The original Runge-Gross (RG) formulation is subsumed within the current formalism. Concerns in the literature about the meaning of the functional derivatives and the internal consistency of the RG formulation are allayed by clarifying the nature of the functional derivatives entering the formalism.

arXiv:cond-mat/0405532•Physical Review A•Mar 31, 2005•8 citations

Solving the Hubbard model using density matrix embedding theory and the variational quantum eigensolver

Lana Mineh, Ashley Montanaro

Calculating the ground state properties of a Hamiltonian can be mapped to the problem of finding the ground state of a smaller Hamiltonian through the use of embedding methods. These embedding techniques have the ability to drastically reduce the problem size, and hence the number of qubits required when running on a quantum computer. However, the embedding process can produce a relatively complicated Hamiltonian, leading to a more complex quantum algorithm. In this paper we carry out a detailed study into how density matrix embedding theory (DMET) could be implemented on a quantum computer to solve the Hubbard model. We consider the variational quantum eigensolver (VQE) as the solver for the embedded Hamiltonian within the DMET algorithm. We derive the exact form of the embedded Hamiltonian and use it to construct efficient ansatz circuits and measurement schemes. We conduct detailed numerical simulations up to 16 qubits, the largest to date, for a range of Hubbard model parameters and find that the combination of DMET and VQE is effective for reproducing ground state properties of the model.

arXiv:2108.08611•Physical review. B./Physical review. B•Mar 14, 2022•19 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

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

Enhancing density functional theory using the variational quantum eigensolver

Evan Sheridan, Lana Mineh, Raul A. Santos +1 more

Quantum computers open up new avenues for modelling the physical properties of materials and molecules. Density Functional Theory (DFT) is the gold standard classical algorithm for predicting these properties, but relies on approximations of the unknown universal functional, limiting its general applicability for many fundamental and technologically relevant systems. In this work we develop a hybrid quantum/classical algorithm called quantum enhanced DFT (QEDFT) that systematically constructs quantum approximations of the universal functional using data obtained from a quantum computer. We benchmark the QEDFT algorithm on the Fermi-Hubbard model, both numerically and on data from experiments on real quantum hardware. We find that QEDFT surpasses the quality of groundstate results obtained from Hartree-Fock DFT, as well as from direct application of conventional quantum algorithms such as VQE. Furthermore, we demonstrate that QEDFT works even when only noisy, low-depth quantum computation is available, by benchmarking the algorithm on data obtained from Google's quantum computer. We further show how QEDFT also captures quintessential properties of strongly correlated Mott physics for large Fermi-Hubbard systems using functionals generated on much smaller system sizes. Our results indicate that QEDFT can be applied to realistic materials and molecular systems, and has the potential to outperform the direct application of either DFT or VQE alone, without the requirement of large scale or fully fault-tolerant quantum computers.

arXiv:2402.18534•Feb 28, 2024•2 citations

Wave function and external potential from constrained search in density-functional theory

Ji Luo

It is proven that in Levy's constrained search or Lieb's ensemble constrained search, if the minimum is achieved by wave functions that satisfy the steadiness condition for the current density kernel and have continuous second derivatives, the minimizing wave functions are eigenstates of some external potential. If the excited states of this potential can only achieve stationarity and cannot achieve a minimum in the constrained search, the minimizing wave functions are ground states and the corresponding electron density is $\mathit{v}$-representable or ensemble $\mathit{v}$-representable. This is also true for the ensemble system with a varying electron number. In time-dependent density-functional theory, if the stationarity of the action integral is achieved in the constrained search by a wave function that satisfies the continuity equation for the current density kernel and has a continuous first derivative with respect to the time coordinate and continuous second derivatives with respect to the space coordinates, the stationary wave function satisfies the time-dependent Schr\"odinger equation of some time-dependent external potential and the corresponding electron density is time-dependent $\mathit{v}$-representable.

Physical Review A•Jul 16, 2008•0 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
The constrained-search formulation of Levy and Lieb provides a concrete mapping from NNN-representable densities to the space of NNN-particle wave functions and explicitly defines the universal functional of density-functional theory. We numerically implement the Levy-Lieb procedure for a paradigmatic lattice system, the Hubbard dimer, using a modified variational quantum eigensolver approach. We demonstrate density variational minimization using the resulting hybrid quantum-classical scheme featuring real-time computation of the Levy-Lieb functional along the search trajectory. We further illustrate a fidelity-based quantum kernel associated with the density to pure-state embedding implied by the Levy-Lieb procedure and employ the kernel for learning observable functionals of the density. We study the kernel's ability to generalize with high accuracy through numerical experiments on the Hubbard dimer.
To construct an NNN-representable time-dependent density-functional theory, a generalization to the time domain of the Levy-Lieb (LL) constrained-search algorithm is required. That the action is only stationary in the Dirac-Frenkel variational principle eliminates the possibility of basing the search on the action itself. Instead, we use the norm of the partial functional derivative of the action in the Hilbert space of the wave functions in place of the energy of the LL search. The electron densities entering the formalism are NNN-representable, and the resulting universal action functional has a unique stationary point in the density at that corresponding to the solution of the Schr\"odinger equation. The original Runge-Gross (RG) formulation is subsumed within the current formalism. Concerns in the literature about the meaning of the functional derivatives and the internal consistency of the RG formulation are allayed by clarifying the nature of the functional derivatives entering the formalism.
Calculating the ground state properties of a Hamiltonian can be mapped to the problem of finding the ground state of a smaller Hamiltonian through the use of embedding methods. These embedding techniques have the ability to drastically reduce the problem size, and hence the number of qubits required when running on a quantum computer. However, the embedding process can produce a relatively complicated Hamiltonian, leading to a more complex quantum algorithm. In this paper we carry out a detailed study into how density matrix embedding theory (DMET) could be implemented on a quantum computer to solve the Hubbard model. We consider the variational quantum eigensolver (VQE) as the solver for the embedded Hamiltonian within the DMET algorithm. We derive the exact form of the embedded Hamiltonian and use it to construct efficient ansatz circuits and measurement schemes. We conduct detailed numerical simulations up to 16 qubits, the largest to date, for a range of Hubbard model parameters and find that the combination of DMET and VQE is effective for reproducing ground state properties of the model.
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.
Quantum computers open up new avenues for modelling the physical properties of materials and molecules. Density Functional Theory (DFT) is the gold standard classical algorithm for predicting these properties, but relies on approximations of the unknown universal functional, limiting its general applicability for many fundamental and technologically relevant systems. In this work we develop a hybrid quantum/classical algorithm called quantum enhanced DFT (QEDFT) that systematically constructs quantum approximations of the universal functional using data obtained from a quantum computer. We benchmark the QEDFT algorithm on the Fermi-Hubbard model, both numerically and on data from experiments on real quantum hardware. We find that QEDFT surpasses the quality of groundstate results obtained from Hartree-Fock DFT, as well as from direct application of conventional quantum algorithms such as VQE. Furthermore, we demonstrate that QEDFT works even when only noisy, low-depth quantum computation is available, by benchmarking the algorithm on data obtained from Google's quantum computer. We further show how QEDFT also captures quintessential properties of strongly correlated Mott physics for large Fermi-Hubbard systems using functionals generated on much smaller system sizes. Our results indicate that QEDFT can be applied to realistic materials and molecular systems, and has the potential to outperform the direct application of either DFT or VQE alone, without the requirement of large scale or fully fault-tolerant quantum computers.
It is proven that in Levy's constrained search or Lieb's ensemble constrained search, if the minimum is achieved by wave functions that satisfy the steadiness condition for the current density kernel and have continuous second derivatives, the minimizing wave functions are eigenstates of some external potential. If the excited states of this potential can only achieve stationarity and cannot achieve a minimum in the constrained search, the minimizing wave functions are ground states and the corresponding electron density is v\mathit{v}v-representable or ensemble v\mathit{v}v-representable. This is also true for the ensemble system with a varying electron number. In time-dependent density-functional theory, if the stationarity of the action integral is achieved in the constrained search by a wave function that satisfies the continuity equation for the current density kernel and has a continuous first derivative with respect to the time coordinate and continuous second derivatives with respect to the space coordinates, the stationary wave function satisfies the time-dependent Schr\"odinger equation of some time-dependent external potential and the corresponding electron density is time-dependent v\mathit{v}v-representable.
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.