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:2106.13583 · Physical review. B./Physical review. B

Simulating both parity sectors of the Hubbard model with tensor networks

Manuel Schneider, Johann Ostmeyer, Karl Jansen, Thomas Luu, Carsten Urbach

Published Oct 11, 2021• 15 citations
arXivHTMLPDF

Abstract

Tensor networks are a powerful tool to simulate a variety of different physical models, including those that suffer from the sign problem in Monte Carlo simulations. The Hubbard model on the honeycomb lattice with nonzero chemical potential is one such problem. Our method is based on projected entangled pair states using imaginary-time evolution. We demonstrate that it provides accurate estimators for the ground state of the model, including cases where Monte Carlo simulations fail miserably. In particular, it shows near to optimal, that is linear, scaling in lattice size. We also present an approach to directly simulate the subspace with an odd number of fermions. It allows to independently determine the ground state in both sectors. Without a chemical potential, this corresponds to half-filling and the lowest-energy state with one additional electron or hole. We identify several stability issues, such as degenerate ground states and large single-particle gaps, and provide possible fixes.

DOI
10.1103/physrevb.104.155118

Related papers

The Hubbard model with fermionic tensor networks

Manuel Schneider, Johann Ostmeyer, Karl Jansen +2 more

Many electromagnetic properties of graphene can be described by the Hubbard model on a honeycomb lattice. However, this system suffers strongly from the sign problem if a chemical potential is included. Tensor network methods are not affected by this problem. We use the imaginary time evolution of a fermionic projected entangled pair state, which allows to simulate both parity sectors independently. Incorporating the fermionic nature on the level of the tensor network allows to fix the particle number to be either even or odd. This way we can access the states at half filling and with one additional electron. We calculate the energy and other observables of both states, which was not possible before with Monte Carlo methods.

arXiv:2110.15340•International Symposium on Lattice Field Theory (LATTICE)•May 16, 2022•3 citations

Accurate Simulation of the Hubbard Model with Finite Fermionic Projected Entangled Pair States

Wen-Yuan Liu, Huanchen Zhai, Ruojing Peng +2 more

We demonstrate the use of finite-size fermionic projected entangled pair states, in conjunction with variational Monte Carlo, to perform accurate simulations of the ground state of the 2D Hubbard model. Using bond dimensions of up to D=28, we show that we can surpass state-of-the-art density matrix renormalization group energies that use up to m=32 000 SU(2) multiplets on eight-leg ladders. We further apply our methodology to 10×16, 12×16, and 16×16 lattices at 1/8 hole doping and observe the dimensional crossover between stripe orientations. Our Letter shows the power of finite-size fermionic tensor networks to resolve the physics of the 2D Hubbard model and related problems.

arXiv:2502.13454•Physical Review Letters•Jun 3, 2025•22 citations

Stochastic and Tensor Network simulations of the Hubbard Model

Johann Ostmeyer

The Hubbard model is an important tool to understand the electrical properties of various materials. More specifically, on the honeycomb lattice it is used to describe graphene predicting a quantum phase transition from a semimetal to a Mott insulating state. In this work two different numerical techniques are presented that have been employed for simulations of the Hubbard model: The Hybrid Monte Carlo algorithm on the one hand allowed us to simulate unprecedentedly large lattices, whereas Tensor Networks can be used to completely avoid the sign problem. Respective strengths and weaknesses of the methods are discussed.

arXiv:2210.06874•International Symposium on Lattice Field Theory (LATTICE)•Nov 7, 2022•1 citation

Thermal tensor network simulations of lattice fermions with fixed filling

Qiaoyi Li, Dai-Wei Qu, Bin-Bin Chen +2 more

Numerical simulations of strongly correlated fermions at finite temperature are essential for studying high-temperature superconductivity and other quantum many-body phenomena. The recently developed tangent-space tensor renormalization group (tanTRG) provides an efficient and accurate framework by representing thermal density operators as matrix product operators. However, the particle number generally varies during the cooling process. The conventional strategy of fine-tuning chemical potentials to reach a target filling is computationally demanding. Here we propose a fixed-$N$ tanTRG algorithm that stabilizes the average particle number by adaptively tuning the chemical potential within the imaginary-time evolution. We benchmark its accuracy on exactly solvable free fermions, and further apply it to the square-lattice Hubbard model. For hole-doped cases, we study the temperature evolution of charge and spin correlations, identifying several characteristic temperature scales for stripe formation. Our results establish fixed-$N$ tanTRG as an efficient and reliable tool for finite-temperature studies of correlated fermion systems.

arXiv:2511.07303•Physical review. B./Physical review. B•Feb 2, 2026•2 citations

Projection of infinite- U Hubbard model and algebraic sign structure

Yunqing Ouyang, Xiao Yan Xu

We propose a projection approach to perform quantum Monte Carlo (QMC) simulation on the infinite-$U$ Hubbard model at some integer fillings where either it is sign problem free or surprisingly has an algebraic sign structure---a power law dependence of average sign on system size. We demonstrate our scheme on the infinite-$U\phantom{\rule{4pt}{0ex}}SU(2N)$ fermionic Hubbard model on both a square and honeycomb lattice at half filling, where it is sign problem free, and suggest possible correlated ground states. The method can be generalized to study certain extended Hubbard models applying to cluster Mott insulators or two-dimensional Moir\'e systems; among one of them at certain non-half-integer filling, the sign has an algebraic behavior such that it can be numerically solved within a polynomial time. Further, our projection scheme can also be generalized to implement the Gutzwiller projection to spin basis such that $SU(2N)$ quantum spin models and Kondo lattice models may be studied in the framework of fermionic QMC simulations.

arXiv:2108.04832•Physical review. B./Physical review. B•Dec 14, 2021•15 citations

Finite-temperature tensor network study of the Hubbard model on an infinite square lattice

Aritra Sinha, Marek M. Rams, Piotr Czarnik +1 more

The Hubbard model is a longstanding problem in the theory of strongly correlated electrons and a very active one in the experiments with ultracold fermionic atoms. Motivated by current and prospective quantum simulations, we apply a two-dimensional tensor network---an infinite projected entangled pair state---evolved in imaginary time by the neighborhood tensor update algorithm working directly in the thermodynamic limit. With $\mathrm{U}(1)\ifmmode\times\else\texttimes\fi{}\mathrm{U}(1)$ symmetry and the bond dimensions up to 29, we generate thermal states down to the temperature of 0.17 times the hopping rate. We obtain results for spin and charge correlators, unaffected by boundary effects. The spin correlators---measurable in prospective ultracold atoms experiments attempting to approach the thermodynamic limit---provide evidence of disruption of the antiferromagnetic background with mobile holes in a slightly doped Hubbard model. The charge correlators reveal the presence of hole-doublon pairs near half filling and signatures of hole-hole repulsion on doping. We also obtain specific heat in the slightly doped regime.

arXiv:2209.00985•Physical review. B./Physical review. B•Nov 2, 2022•24 citations
Tensor networks are a powerful tool to simulate a variety of different physical models, including those that suffer from the sign problem in Monte Carlo simulations. The Hubbard model on the honeycomb lattice with nonzero chemical potential is one such problem. Our method is based on projected entangled pair states using imaginary-time evolution. We demonstrate that it provides accurate estimators for the ground state of the model, including cases where Monte Carlo simulations fail miserably. In particular, it shows near to optimal, that is linear, scaling in lattice size. We also present an approach to directly simulate the subspace with an odd number of fermions. It allows to independently determine the ground state in both sectors. Without a chemical potential, this corresponds to half-filling and the lowest-energy state with one additional electron or hole. We identify several stability issues, such as degenerate ground states and large single-particle gaps, and provide possible fixes.
Many electromagnetic properties of graphene can be described by the Hubbard model on a honeycomb lattice. However, this system suffers strongly from the sign problem if a chemical potential is included. Tensor network methods are not affected by this problem. We use the imaginary time evolution of a fermionic projected entangled pair state, which allows to simulate both parity sectors independently. Incorporating the fermionic nature on the level of the tensor network allows to fix the particle number to be either even or odd. This way we can access the states at half filling and with one additional electron. We calculate the energy and other observables of both states, which was not possible before with Monte Carlo methods.
We demonstrate the use of finite-size fermionic projected entangled pair states, in conjunction with variational Monte Carlo, to perform accurate simulations of the ground state of the 2D Hubbard model. Using bond dimensions of up to D=28, we show that we can surpass state-of-the-art density matrix renormalization group energies that use up to m=32 000 SU(2) multiplets on eight-leg ladders. We further apply our methodology to 10×16, 12×16, and 16×16 lattices at 1/8 hole doping and observe the dimensional crossover between stripe orientations. Our Letter shows the power of finite-size fermionic tensor networks to resolve the physics of the 2D Hubbard model and related problems.
The Hubbard model is an important tool to understand the electrical properties of various materials. More specifically, on the honeycomb lattice it is used to describe graphene predicting a quantum phase transition from a semimetal to a Mott insulating state. In this work two different numerical techniques are presented that have been employed for simulations of the Hubbard model: The Hybrid Monte Carlo algorithm on the one hand allowed us to simulate unprecedentedly large lattices, whereas Tensor Networks can be used to completely avoid the sign problem. Respective strengths and weaknesses of the methods are discussed.
Numerical simulations of strongly correlated fermions at finite temperature are essential for studying high-temperature superconductivity and other quantum many-body phenomena. The recently developed tangent-space tensor renormalization group (tanTRG) provides an efficient and accurate framework by representing thermal density operators as matrix product operators. However, the particle number generally varies during the cooling process. The conventional strategy of fine-tuning chemical potentials to reach a target filling is computationally demanding. Here we propose a fixed-NNN tanTRG algorithm that stabilizes the average particle number by adaptively tuning the chemical potential within the imaginary-time evolution. We benchmark its accuracy on exactly solvable free fermions, and further apply it to the square-lattice Hubbard model. For hole-doped cases, we study the temperature evolution of charge and spin correlations, identifying several characteristic temperature scales for stripe formation. Our results establish fixed-NNN tanTRG as an efficient and reliable tool for finite-temperature studies of correlated fermion systems.
We propose a projection approach to perform quantum Monte Carlo (QMC) simulation on the infinite-UUU Hubbard model at some integer fillings where either it is sign problem free or surprisingly has an algebraic sign structure---a power law dependence of average sign on system size. We demonstrate our scheme on the infinite-USU(2N)U\phantom{\rule{4pt}{0ex}}SU(2N)USU(2N) fermionic Hubbard model on both a square and honeycomb lattice at half filling, where it is sign problem free, and suggest possible correlated ground states. The method can be generalized to study certain extended Hubbard models applying to cluster Mott insulators or two-dimensional Moir\'e systems; among one of them at certain non-half-integer filling, the sign has an algebraic behavior such that it can be numerically solved within a polynomial time. Further, our projection scheme can also be generalized to implement the Gutzwiller projection to spin basis such that SU(2N)SU(2N)SU(2N) quantum spin models and Kondo lattice models may be studied in the framework of fermionic QMC simulations.
The Hubbard model is a longstanding problem in the theory of strongly correlated electrons and a very active one in the experiments with ultracold fermionic atoms. Motivated by current and prospective quantum simulations, we apply a two-dimensional tensor network---an infinite projected entangled pair state---evolved in imaginary time by the neighborhood tensor update algorithm working directly in the thermodynamic limit. With $\mathrm{U}(1)\ifmmode\times\else\texttimes\fi{}\mathrm{U}(1)$ symmetry and the bond dimensions up to 29, we generate thermal states down to the temperature of 0.17 times the hopping rate. We obtain results for spin and charge correlators, unaffected by boundary effects. The spin correlators---measurable in prospective ultracold atoms experiments attempting to approach the thermodynamic limit---provide evidence of disruption of the antiferromagnetic background with mobile holes in a slightly doped Hubbard model. The charge correlators reveal the presence of hole-doublon pairs near half filling and signatures of hole-hole repulsion on doping. We also obtain specific heat in the slightly doped regime.