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:1912.11907 · Physical Review Letters

Variational Approach for Many-Body Systems at Finite Temperature

Tao Shi, Eugene Demler, J. I. Cirac

Published Oct 28, 2020• 23 citations
arXivHTMLPDF

Abstract

We introduce an equation for density matrices that ensures a monotonic decrease of the free energy and reaches a fixed point at the Gibbs thermal. We build a variational approach for many-body systems that can be applied to a broad class of states, including all bosonic and fermionic Gaussian, as well as their generalizations obtained by unitary transformations, such as polaron transformations in electron-phonon systems. We apply it to the Holstein model on 20×20 and 50×50 square lattices, and predict phase separation between the superconducting and charge-density wave phases in the strong interaction regime.

DOI
10.1103/physrevlett.125.180602

Related papers

Variational study of fermionic and bosonic systems with non-Gaussian states: Theory and applications

Shi, T., Eugene Demler, J. I. Cirac

We present a new variational method for investigating the ground state and out of equilibrium dynamics of quantum many-body bosonic and fermionic systems. Our approach is based on constructing variational wavefunctions which extend Gaussian states by including generalized canonical transformations between the fields. The key advantage of such states compared to simple Gaussian states is presence of non-factorizable correlations and the possibility of describing states with strong entanglement between particles. In contrast to the commonly used canonical transformations, such as the polaron or Lang-Firsov transformations, we allow parameters of the transformations to be time dependent, which extends their regions of applicability. We derive equations of motion for the parameters characterizing the states both in real and imaginary time using the differential structure of the variational manifold. The ground state can be found by following the imaginary time evolution until it converges to a steady state. Collective excitations in the system can be obtained by linearizing the real-time equations of motion in the vicinity of the imaginary time steady-state solution. Our formalism allows us not only to determine the energy spectrum of quasiparticles and their lifetime, but to obtain the complete spectral functions and to explore far out of equilibrium dynamics such as coherent evolution following a quantum quench. We illustrate and benchmark this framework with several examples: a single polaron in the Holstein and Su-Schrieer-Heeger models, non-equilibrium dynamics in the spin-boson and Kondo models, the superconducting to charge density wave phase transitions in the Holstein model.

•Mar 1, 2018•106 citations

Holstein polaron

J. Bonča, S. A. Trugman, I. Batistić

We describe a variational method to solve the Holstein model for an electron coupled to dynamical, quantum phonons on an infinite lattice. The variational space can be systematically expanded to achieve high accuracy with modest computational resources (12-digit accuracy for the one-dimensional polaron energy at intermediate coupling). We compute ground-state and low-lying excited-state properties of the model at continuous values of the wave vector k in essentially all parameter regimes. Our results for the polaron energy band, effective mass, and correlation functions compare favorably with those of other numerical techniques, including the density-matrix renormalization-group technique, the global-local method, and the exact diagonalization technique. We find a phase transition for the first excited state between a bound and unbound system of a polaron and an additional phonon excitation. The phase transition is also treated in strong-coupling perturbation theory.

arXiv:cond-mat/9812252•Physical review. B, Condensed matter•Jul 15, 1999•281 citations

Non-Gaussian variational wave functions for interacting bosons on a lattice

Tian Qian, José J. Fernández-Melgarejo, David Zueco +1 more

A variational method for studying the ground state of strongly interacting quantum many-body bosonic systems is presented. Our approach constructs a class of extensive variational non-Gaussian wave functions which extend Gaussian states by means of nonlinear canonical transformations (NLCTs) on the fields of the theory under consideration. We illustrate this method with the one-dimensional Bose-Hubbard model for which the proposal presented here provides a family of approximate ground states at arbitrarily large values of the interaction strength. We find that, for different values of the interaction, the non-Gaussian NLCT-trial states sensibly improve the ground-state energy estimation when the system is in the Mott phase.

arXiv:2211.04320•Physical review. B./Physical review. B•Jan 17, 2023•1 citation

Fermionic neural Gibbs states

Jannes Nys, Juan Carrasquilla

We introduce fermionic neural Gibbs states (fNGS), a variational framework for modeling finite-temperature properties of strongly interacting fermions. fNGS starts from a reference mean-field thermofield-double state and uses neural-network transformations together with imaginary-time evolution to systematically build strong correlations. Applied to the doped Fermi-Hubbard model, a minimal lattice model capturing essential features of strong electronic correlations, fNGS accurately reproduces thermal energies over a broad range of temperatures, interaction strengths, even at large dopings, for system sizes beyond the reach of exact methods. These results demonstrate a scalable route to studying finite-temperature properties of strongly correlated fermionic systems beyond one dimension with neural-network representations of quantum states.

quant-pharXiv:2512.04663•Dec 5, 2025•0 citations

Reduction of the N -Particle Variational Problem

Claude Garrod, Jerome K. Percus

A variational method is presented which is applicable to N-particle boson or fermion systems with two-body interactions. For these systems the energy may be expressed in terms of the two-particle density matrix: Γ(1, 2 | 1′, 2′)=(Ψ |a2+a1+a1′a2′| Ψ). In order to have the variational equation: δE/δΓ = 0 yield the correct ground-state density matrix one must restrict Γ to the set of density matrices which are actually derivable from N-particle boson (or fermion) systems. Subsidiary conditions are presented which are necessary and sufficient to insure that Γ is so derivable. These conditions are of a form which render them unsuited for practical application. However the following necessary (but not sufficient) conditions are shown by some applications to yield good results: It is proven that if Γ(1, 2 | 1′, 2′) and γ(1 | 1′) are the two-particle and one-particle density matrices of an N-particle system [normalized by tr Γ = N(N − 1) and trγ = N] then the associated operator: G(1,2 | 1′,2′)=δ(1−1′)γ(2 | 2′)+σ Γ(1′,2 | 1,2′)−γ(2 | 1)γ(1′ | 2′) is a nonnegative operator. [Here σ is + 1 or − 1 for bosons or fermions respectively.]

Journal of Mathematical Physics•Dec 1, 1964•477 citations

Superconductivity, charge density waves, and bipolarons in the Holstein model

Benjamin Nosarzewski, Edwin W. Huang, Philip M. Dee +5 more

The electron-phonon $(e\text{\ensuremath{-}}\mathrm{ph})$ interaction remains of great interest in condensed matter physics and plays a vital role in realizing superconductors, charge density waves (CDW), and polarons. We study the two-dimensional Holstein model for $e\text{\ensuremath{-}}\mathrm{ph}$ coupling using determinant quantum Monte Carlo across a wide range of its phase diagram as a function of temperature, electron density, dimensionless $e\text{\ensuremath{-}}\mathrm{ph}$ coupling strength, and the adiabatic ratio of the phonon frequency to the Fermi energy. We describe the behavior of the CDW correlations, the competition between superconducting and CDW orders and polaron formation, the optimal conditions for superconductivity, and the transition from the weak-coupling regime to the strong-coupling regime. Superconductivity is optimized at intermediate $e\text{\ensuremath{-}}\mathrm{ph}$ coupling strength and intermediate electron density, and the superconducting correlations increase monotonically with phonon frequency. The global maximum for superconductivity in the Holstein model occurs at large phonon frequency, the limit where an attractive Hubbard model effectively describes the physics.

arXiv:2103.09242•Physical review. B./Physical review. B•Jun 24, 2021•37 citations
We introduce an equation for density matrices that ensures a monotonic decrease of the free energy and reaches a fixed point at the Gibbs thermal. We build a variational approach for many-body systems that can be applied to a broad class of states, including all bosonic and fermionic Gaussian, as well as their generalizations obtained by unitary transformations, such as polaron transformations in electron-phonon systems. We apply it to the Holstein model on 20×20 and 50×50 square lattices, and predict phase separation between the superconducting and charge-density wave phases in the strong interaction regime.
We present a new variational method for investigating the ground state and out of equilibrium dynamics of quantum many-body bosonic and fermionic systems. Our approach is based on constructing variational wavefunctions which extend Gaussian states by including generalized canonical transformations between the fields. The key advantage of such states compared to simple Gaussian states is presence of non-factorizable correlations and the possibility of describing states with strong entanglement between particles. In contrast to the commonly used canonical transformations, such as the polaron or Lang-Firsov transformations, we allow parameters of the transformations to be time dependent, which extends their regions of applicability. We derive equations of motion for the parameters characterizing the states both in real and imaginary time using the differential structure of the variational manifold. The ground state can be found by following the imaginary time evolution until it converges to a steady state. Collective excitations in the system can be obtained by linearizing the real-time equations of motion in the vicinity of the imaginary time steady-state solution. Our formalism allows us not only to determine the energy spectrum of quasiparticles and their lifetime, but to obtain the complete spectral functions and to explore far out of equilibrium dynamics such as coherent evolution following a quantum quench. We illustrate and benchmark this framework with several examples: a single polaron in the Holstein and Su-Schrieer-Heeger models, non-equilibrium dynamics in the spin-boson and Kondo models, the superconducting to charge density wave phase transitions in the Holstein model.
We describe a variational method to solve the Holstein model for an electron coupled to dynamical, quantum phonons on an infinite lattice. The variational space can be systematically expanded to achieve high accuracy with modest computational resources (12-digit accuracy for the one-dimensional polaron energy at intermediate coupling). We compute ground-state and low-lying excited-state properties of the model at continuous values of the wave vector k in essentially all parameter regimes. Our results for the polaron energy band, effective mass, and correlation functions compare favorably with those of other numerical techniques, including the density-matrix renormalization-group technique, the global-local method, and the exact diagonalization technique. We find a phase transition for the first excited state between a bound and unbound system of a polaron and an additional phonon excitation. The phase transition is also treated in strong-coupling perturbation theory.
A variational method for studying the ground state of strongly interacting quantum many-body bosonic systems is presented. Our approach constructs a class of extensive variational non-Gaussian wave functions which extend Gaussian states by means of nonlinear canonical transformations (NLCTs) on the fields of the theory under consideration. We illustrate this method with the one-dimensional Bose-Hubbard model for which the proposal presented here provides a family of approximate ground states at arbitrarily large values of the interaction strength. We find that, for different values of the interaction, the non-Gaussian NLCT-trial states sensibly improve the ground-state energy estimation when the system is in the Mott phase.
We introduce fermionic neural Gibbs states (fNGS), a variational framework for modeling finite-temperature properties of strongly interacting fermions. fNGS starts from a reference mean-field thermofield-double state and uses neural-network transformations together with imaginary-time evolution to systematically build strong correlations. Applied to the doped Fermi-Hubbard model, a minimal lattice model capturing essential features of strong electronic correlations, fNGS accurately reproduces thermal energies over a broad range of temperatures, interaction strengths, even at large dopings, for system sizes beyond the reach of exact methods. These results demonstrate a scalable route to studying finite-temperature properties of strongly correlated fermionic systems beyond one dimension with neural-network representations of quantum states.
A variational method is presented which is applicable to N-particle boson or fermion systems with two-body interactions. For these systems the energy may be expressed in terms of the two-particle density matrix: Γ(1, 2 | 1′, 2′)=(Ψ |a2+a1+a1′a2′| Ψ). In order to have the variational equation: δE/δΓ = 0 yield the correct ground-state density matrix one must restrict Γ to the set of density matrices which are actually derivable from N-particle boson (or fermion) systems. Subsidiary conditions are presented which are necessary and sufficient to insure that Γ is so derivable. These conditions are of a form which render them unsuited for practical application. However the following necessary (but not sufficient) conditions are shown by some applications to yield good results: It is proven that if Γ(1, 2 | 1′, 2′) and γ(1 | 1′) are the two-particle and one-particle density matrices of an N-particle system [normalized by tr Γ = N(N − 1) and trγ = N] then the associated operator: G(1,2 | 1′,2′)=δ(1−1′)γ(2 | 2′)+σ Γ(1′,2 | 1,2′)−γ(2 | 1)γ(1′ | 2′) is a nonnegative operator. [Here σ is + 1 or − 1 for bosons or fermions respectively.]
The electron-phonon $(e\text{\ensuremath{-}}\mathrm{ph})$ interaction remains of great interest in condensed matter physics and plays a vital role in realizing superconductors, charge density waves (CDW), and polarons. We study the two-dimensional Holstein model for $e\text{\ensuremath{-}}\mathrm{ph}$ coupling using determinant quantum Monte Carlo across a wide range of its phase diagram as a function of temperature, electron density, dimensionless $e\text{\ensuremath{-}}\mathrm{ph}$ coupling strength, and the adiabatic ratio of the phonon frequency to the Fermi energy. We describe the behavior of the CDW correlations, the competition between superconducting and CDW orders and polaron formation, the optimal conditions for superconductivity, and the transition from the weak-coupling regime to the strong-coupling regime. Superconductivity is optimized at intermediate $e\text{\ensuremath{-}}\mathrm{ph}$ coupling strength and intermediate electron density, and the superconducting correlations increase monotonically with phonon frequency. The global maximum for superconductivity in the Holstein model occurs at large phonon frequency, the limit where an attractive Hubbard model effectively describes the physics.