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

Accurate simulation for finite projected entangled pair states in two dimensions

Wen-Yuan Liu, Yi-Zhen Huang, Shou-Shu Gong, Zheng‐Cheng Gu

Published Jun 24, 2021• 35 citations
arXivHTMLPDF

Abstract

Based on the scheme of variational Monte Carlo sampling, we develop an accurate and efficient two-dimensional tensor-network algorithm to simulate quantum lattice models. We find that Monte Carlo sampling shows huge advantages in dealing with finite projected entangled pair states, which allows significantly enlarged system size and improves the accuracy of tensor network simulation. We demonstrate our method on the square-lattice antiferromagnetic Heisenberg model up to $32\ifmmode\times\else\texttimes\fi{}32$ sites, as well as a highly frustrated ${J}_{1}\ensuremath{-}{J}_{2}$ model up to $24\ifmmode\times\else\texttimes\fi{}24$ sites. The results, including ground state energy and spin correlations, are in excellent agreement with those of the available quantum Monte Carlo or density matrix renormalization group methods. Therefore, our method substantially advances the calculation of 2D tensor networks for finite systems, and potentially opens a new door towards resolving many challenging strongly correlated quantum many-body problems.

DOI
10.1103/physrevb.103.235155

Related papers

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

Green's function Monte Carlo combined with projected entangled pair state approach to the frustrated J 1 − J 2 Heisenberg model

He-Yu Lin, Yibin Guo, Rong-Qiang He +2 more

The tensor network algorithm, a family of prevalent numerical methods for quantum many-body problems, aptly captures the entanglement properties intrinsic to quantum systems, enabling precise representation of quantum states. However, its computational cost is notably high, particularly in calculating physical observables like correlation functions. To surmount the computational challenge and enhance efficiency, we propose integrating the Green's function Monte Carlo (GFMC) method with the projected entangled pair state (PEPS) ansatz. This approach combines the high-efficiency characteristics of Monte Carlo with the sign-free nature of tensor network states and proves effective in addressing the computational bottleneck. To showcase its prowess, we apply this hybrid approach to investigate the antiferromagnetic ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ Heisenberg model on the square lattice, a model notorious for its sign problem in quantum Monte Carlo simulations. Our results reveal a substantial improvement in the accuracy of ground-state energy when utilizing a preliminary PEPS as the guiding wave function for GFMC. By calculating the structure factor and spin-spin correlation functions, we further characterize the phase diagram, identifying a possible columnar valence-bond state phase within the intermediate parameter range of $0.52<{J}_{2}/{J}_{1}<0.58$. This comprehensive study underscores the efficacy of our combined approach, demonstrating its ability to accurately simulate frustrated quantum spin systems while ensuring computational efficiency.

arXiv:2406.12207•Physical review. B./Physical review. B•Jun 17, 2024•6 citations

Green's function Monte Carlo combined with projected entangled pair state approach to the frustrated J1-J2 Heisenberg model

Lin, He-Yu, Yibin Guo, Rong-Qiang He +2 more

The tensor network algorithm, a family of prevalent numerical methods for quantum many-body problems, aptly captures the entanglement properties intrinsic to quantum systems, enabling precise representation of quantum states. However, its computational cost is notably high, particularly in calculating physical observables like correlation functions. To surmount the computational challenge and enhance efficiency, we propose integrating the Green's function Monte Carlo (GFMC) method with the projected entangled pair state (PEPS) ansatz. This approach combines the high-efficiency characteristics of Monte Carlo with the sign-free nature of tensor network states and proves effective in addressing the computational bottleneck. To showcase its prowess, we apply this hybrid approach to investigate the antiferromagnetic J1-J2 Heisenberg model on the square lattice, a model notorious for its sign problem in quantum Monte Carlo simulations. Our results reveal a substantial improvement in the accuracy of ground-state energy when utilizing a preliminary PEPS as the guiding wave function for GFMC. By calculating the structure factor and spin-spin correlation functions, we further characterize the phase diagram, identifying a possible columnar valence-bond state phase within the intermediate parameter range of 0.52<J2/J1<0.58. This comprehensive study underscores the efficacy of our combined approach, demonstrating its ability to accurately simulate frustrated quantum spin systems while ensuring computational efficiency.

•Jan 1, 2024•1 citation

Scalable tensor network algorithm for thermal quantum many-body systems in two dimension

Meng Zhang, Hao Zhang, Chao Wang +1 more

Simulating strongly-correlated quantum many-body systems at finite temperatures is a significant challenge in computational physics. In this work, we present a scalable finite-temperature tensor network algorithm for two-dimensional quantum many-body systems. We employ the (fermionic) projected entangled pair state (PEPS) to represent the vectorization of the quantum thermal state and utilize a stochastic reconfiguration method to cool down the quantum states from infinite temperature. We validate our method by benchmarking it against the 2D antiferromagnetic Heisenberg model, the $J_1$-$J_2$ model, and the Fermi-Hubbard model, comparing physical properties such as internal energy, specific heat, and magnetic susceptibility with results obtained from stochastic series expansion (SSE), exact diagonalization, and determinant quantum Monte Carlo (DQMC).

arXiv:2409.05285•Sep 9, 2024•0 citations

Direct sampling of projected entangled-pair states

Tom Vieijra, Jutho Haegeman, Frank Verstraete +1 more

Variational Monte Carlo studies employing projected entangled-pair states (PEPS) have recently shown that they can provide answers to long-standing questions such as the nature of the phases in the two-dimensional ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ model. The sampling in these Monte Carlo algorithms is typically performed with Markov chain Monte Carlo algorithms employing local update rules, which often suffer from long autocorrelation times and interdependent samples. We propose a sampling algorithm that generates independent samples from a PEPS, bypassing all problems related to finite autocorrelation times. This algorithm is a generalization of an existing direct sampling algorithm for unitary tensor networks. We introduce an auxiliary probability distribution from which independent samples can be drawn, and combine it with importance sampling in order to evaluate expectation values accurately. We benchmark our algorithm on the classical Ising model and on variational optimization of two-dimensional quantum spin models.

arXiv:2109.07356•Physical review. B./Physical review. B•Dec 21, 2021•37 citations

Scalable tensor network algorithm for thermal quantum many-body systems in two dimensions

Meng Zhang, Hao Zhang, Chao Wang +1 more

Simulating strongly correlated quantum many-body systems at finite temperatures is a significant challenge in computational physics. In this work, we present a scalable finite-temperature tensor network algorithm for two-dimensional quantum many-body systems. We employ the (fermionic) projected entangled pair state to represent the vectorization of the quantum thermal state, and we utilize a stochastic reconfiguration method to cool down the quantum states from infinite temperature. We validate our method by benchmarking it against the two-dimensional antiferromagnetic Heisenberg model, the ${J}_{1}\ensuremath{-}{J}_{2}$ model, and the Fermi-Hubbard model, comparing physical properties such as internal energy, specific heat, and magnetic susceptibility with results obtained from stochastic series expansion, exact diagonalization, and determinant quantum Monte Carlo.

Physical review. B./Physical review. B•Feb 19, 2025•2 citations
Based on the scheme of variational Monte Carlo sampling, we develop an accurate and efficient two-dimensional tensor-network algorithm to simulate quantum lattice models. We find that Monte Carlo sampling shows huge advantages in dealing with finite projected entangled pair states, which allows significantly enlarged system size and improves the accuracy of tensor network simulation. We demonstrate our method on the square-lattice antiferromagnetic Heisenberg model up to 32\ifmmode×\else\texttimes\fi3232\ifmmode\times\else\texttimes\fi{}3232\ifmmode×\else\texttimes\fi32 sites, as well as a highly frustrated J1\ensuremath−J2{J}_{1}\ensuremath{-}{J}_{2}J1​\ensuremath−J2​ model up to 24\ifmmode×\else\texttimes\fi2424\ifmmode\times\else\texttimes\fi{}2424\ifmmode×\else\texttimes\fi24 sites. The results, including ground state energy and spin correlations, are in excellent agreement with those of the available quantum Monte Carlo or density matrix renormalization group methods. Therefore, our method substantially advances the calculation of 2D tensor networks for finite systems, and potentially opens a new door towards resolving many challenging strongly correlated quantum many-body problems.
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 tensor network algorithm, a family of prevalent numerical methods for quantum many-body problems, aptly captures the entanglement properties intrinsic to quantum systems, enabling precise representation of quantum states. However, its computational cost is notably high, particularly in calculating physical observables like correlation functions. To surmount the computational challenge and enhance efficiency, we propose integrating the Green's function Monte Carlo (GFMC) method with the projected entangled pair state (PEPS) ansatz. This approach combines the high-efficiency characteristics of Monte Carlo with the sign-free nature of tensor network states and proves effective in addressing the computational bottleneck. To showcase its prowess, we apply this hybrid approach to investigate the antiferromagnetic ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ Heisenberg model on the square lattice, a model notorious for its sign problem in quantum Monte Carlo simulations. Our results reveal a substantial improvement in the accuracy of ground-state energy when utilizing a preliminary PEPS as the guiding wave function for GFMC. By calculating the structure factor and spin-spin correlation functions, we further characterize the phase diagram, identifying a possible columnar valence-bond state phase within the intermediate parameter range of 0.52<J2/J1<0.580.52<{J}_{2}/{J}_{1}<0.580.52<J2​/J1​<0.58. This comprehensive study underscores the efficacy of our combined approach, demonstrating its ability to accurately simulate frustrated quantum spin systems while ensuring computational efficiency.
The tensor network algorithm, a family of prevalent numerical methods for quantum many-body problems, aptly captures the entanglement properties intrinsic to quantum systems, enabling precise representation of quantum states. However, its computational cost is notably high, particularly in calculating physical observables like correlation functions. To surmount the computational challenge and enhance efficiency, we propose integrating the Green's function Monte Carlo (GFMC) method with the projected entangled pair state (PEPS) ansatz. This approach combines the high-efficiency characteristics of Monte Carlo with the sign-free nature of tensor network states and proves effective in addressing the computational bottleneck. To showcase its prowess, we apply this hybrid approach to investigate the antiferromagnetic J1-J2 Heisenberg model on the square lattice, a model notorious for its sign problem in quantum Monte Carlo simulations. Our results reveal a substantial improvement in the accuracy of ground-state energy when utilizing a preliminary PEPS as the guiding wave function for GFMC. By calculating the structure factor and spin-spin correlation functions, we further characterize the phase diagram, identifying a possible columnar valence-bond state phase within the intermediate parameter range of 0.52<J2/J1<0.58. This comprehensive study underscores the efficacy of our combined approach, demonstrating its ability to accurately simulate frustrated quantum spin systems while ensuring computational efficiency.
Simulating strongly-correlated quantum many-body systems at finite temperatures is a significant challenge in computational physics. In this work, we present a scalable finite-temperature tensor network algorithm for two-dimensional quantum many-body systems. We employ the (fermionic) projected entangled pair state (PEPS) to represent the vectorization of the quantum thermal state and utilize a stochastic reconfiguration method to cool down the quantum states from infinite temperature. We validate our method by benchmarking it against the 2D antiferromagnetic Heisenberg model, the J1J_1J1​-J2J_2J2​ model, and the Fermi-Hubbard model, comparing physical properties such as internal energy, specific heat, and magnetic susceptibility with results obtained from stochastic series expansion (SSE), exact diagonalization, and determinant quantum Monte Carlo (DQMC).
Variational Monte Carlo studies employing projected entangled-pair states (PEPS) have recently shown that they can provide answers to long-standing questions such as the nature of the phases in the two-dimensional ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ model. The sampling in these Monte Carlo algorithms is typically performed with Markov chain Monte Carlo algorithms employing local update rules, which often suffer from long autocorrelation times and interdependent samples. We propose a sampling algorithm that generates independent samples from a PEPS, bypassing all problems related to finite autocorrelation times. This algorithm is a generalization of an existing direct sampling algorithm for unitary tensor networks. We introduce an auxiliary probability distribution from which independent samples can be drawn, and combine it with importance sampling in order to evaluate expectation values accurately. We benchmark our algorithm on the classical Ising model and on variational optimization of two-dimensional quantum spin models.
Simulating strongly correlated quantum many-body systems at finite temperatures is a significant challenge in computational physics. In this work, we present a scalable finite-temperature tensor network algorithm for two-dimensional quantum many-body systems. We employ the (fermionic) projected entangled pair state to represent the vectorization of the quantum thermal state, and we utilize a stochastic reconfiguration method to cool down the quantum states from infinite temperature. We validate our method by benchmarking it against the two-dimensional antiferromagnetic Heisenberg model, the ${J}_{1}\ensuremath{-}{J}_{2}$ model, and the Fermi-Hubbard model, comparing physical properties such as internal energy, specific heat, and magnetic susceptibility with results obtained from stochastic series expansion, exact diagonalization, and determinant quantum Monte Carlo.