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

Physical review. B./Physical review. B

Improved summations of n -point correlation functions of projected entangled-pair states

Boris Ponsioen, Juraj Hašík, Philippe Corboz

Published Nov 7, 2023• 11 citations
Publisher

Abstract

Numerical treatment of two-dimensional strongly-correlated systems is both extremely challenging and of fundamental importance. Infinite projected entangled-pair states (PEPS), a class of tensor networks, have demonstrated cutting-edge performance for ground-state calculations, working directly in the thermodynamic limit. Furthermore, in recent years the application of PEPS has been extended to also low-lying excited states, using an ansatz that targets quasiparticle states above the ground state with high accuracy. A major technical challenge for those simulations is the accurate evaluation of summations of two- and three-point correlation functions with reasonable computational cost. In this paper, we show how a reformulation of $n$-point functions in the context of PEPS leads to extra contributions to the results that prove to play an important role. Benchmarks for the frustrated ${J}_{1}\ensuremath{-}{J}_{2}$ Heisenberg model illustrate the improved precision, efficiency, and stability of the simulations compared to previous approaches. Leveraging automatic differentiation to generate the most tedious and error-prone parts of the computation, the straightforward implementation presented here is a step towards broader adoption of the PEPS excitation ansatz in future applications.

DOI
10.1103/physrevb.108.195111

Related papers

Improved summations of nnn-point correlation functions of projected entangled-pair states

Boris Ponsioen, Juraj Hašík, Philippe Corboz

Numerical treatment of two dimensional strongly-correlated systems is both extremely challenging and of fundamental importance. Infinite projected entangled-pair states (PEPS), a class of tensor networks, have demonstrated cutting-edge performance for ground state calculations, working directly in the thermodynamic limit. Furthermore, in recent years the application of PEPS has been extended to also low-lying excited states, using an ansatz that targets quasiparticle states above the ground state with high accuracy. A major technical challenge for those simulations is the accurate evaluation of summations of two- and three-point correlation functions with reasonable computational cost. In this work, we show how a reformulation of $n$-point functions in the context of PEPS leads to extra contributions to the results that prove to play an important role. Benchmarks for the frustrated $J_1-J_2$ Heisenberg model illustrate the improved precision, efficiency and stability of the simulations compared to previous approaches. Leveraging automatic differentiation to generate the most tedious and error-prone parts of the computation, the straightforward implementation presented here is a step towards broader adoption of the PEPS excitation ansatz in future applications.

arXiv:2306.13327•Jun 23, 2023•0 citations

Automatic differentiation applied to excitations with projected entangled pair states

Boris Ponsioen, Fakher F. Assaad, Philippe Corboz

The excitation ansatz for tensor networks is a powerful tool for simulating the low-lying quasiparticle excitations above ground states of strongly correlated quantum many-body systems. Recently, the two-dimensional tensor network class of infinite projected entangled-pair states gained new ground state optimization methods based on automatic differentiation, which are at the same time highly accurate and simple to implement. Naturally, the question arises whether these new ideas can also be used to optimize the excitation ansatz, which has recently been implemented in two dimensions as well. In this paper, we describe a straightforward way to reimplement the framework for excitations using automatic differentiation, and demonstrate its performance for the Hubbard model at half filling.

arXiv:2107.03399•SciPost Physics•Jan 6, 2022•48 citations

Quantum critical phase with infinite projected entangled paired states

Didier Poilblanc, Matthieu Mambrini

Whether two-dimensional tensor networks, such as the so-called projected entangled pair states (PEPS), have the potential to describe a genuine quantum critical point or quantum critical phases separating two ordered phases is of much interest for the understanding of frustrated magnets and strongly correlated systems. Here, the authors identify a specific family of SU(2)-invariant PEPS, which provides excellent variational energies for the ${J}_{1}-{J}_{2}$ frustrated Heisenberg model at large frustration. Using a full-update infinite-PEPS scheme directly in the thermodynamic limit, the authors provide evidence of diverging correlation lengths in both the singlet and triplet channels. It is argued that such PEPS give a qualitative description of the critical point or critical phase of the ${J}_{1}-{J}_{2}$ model.

arXiv:1702.05950•Physical review. B./Physical review. B•Jul 12, 2017•81 citations

Excitations with projected entangled pair states using the corner transfer matrix method

Boris Ponsioen, Philippe Corboz

We present an extension of a framework for simulating single quasiparticle or collective excitations on top of strongly correlated quantum many-body ground states using infinite projected entangled pair states, a tensor network ansatz for two-dimensional wave functions in the thermodynamic limit. Our approach performs a systematic summation of locally perturbed states in order to obtain excited eigenstates localized in momentum space, using the corner transfer matrix method, and generalizes the framework to arbitrary unit cell sizes, the implementation of global Abelian symmetries and fermionic systems. Results for several test cases are presented, including the transverse Ising model, the spin-$1/2$ Heisenberg model, and a free fermionic model, to demonstrate the capability of the method to accurately capture dispersions. We also provide insight into the nature of excitations at the $k=(\ensuremath{\pi},0)$ point of the Heisenberg model.

arXiv:2001.02645•Physical review. B./Physical review. B•May 5, 2020•54 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

Simulation of strongly correlated fermions in two spatial dimensions with fermionic projected entangled-pair states

Philippe Corboz, Román Orús, Bela Bauer +1 more

We explain how to implement, in the context of projected entangled-pair states (PEPSs), the general procedure of fermionization of a tensor network introduced in P. Corboz and G. Vidal, Phys. Rev. B 80, 165129 (2009). The resulting fermionic PEPS, similar to previous proposals, can be used to study the ground state of interacting fermions on a two-dimensional lattice. As in the bosonic case, the cost of simulations depends on the amount of entanglement in the ground state and not directly on the strength of interactions. The present formulation of fermionic PEPS leads to a straightforward numerical implementation that allowed us to recycle much of the code for bosonic PEPS. We demonstrate that fermionic PEPS are a useful variational ansatz for interacting fermion systems by computing approximations to the ground state of several models on an infinite lattice. For a model of interacting spinless fermions, ground state energies lower than Hartree-Fock results are obtained, shifting the boundary between the metal and charge-density wave phases. For the $t\text{\ensuremath{-}}J$ model, energies comparable with those of a specialized Gutzwiller-projected ansatz are also obtained.

arXiv:0912.0646•Physical Review B•Apr 8, 2010•323 citations
Numerical treatment of two-dimensional strongly-correlated systems is both extremely challenging and of fundamental importance. Infinite projected entangled-pair states (PEPS), a class of tensor networks, have demonstrated cutting-edge performance for ground-state calculations, working directly in the thermodynamic limit. Furthermore, in recent years the application of PEPS has been extended to also low-lying excited states, using an ansatz that targets quasiparticle states above the ground state with high accuracy. A major technical challenge for those simulations is the accurate evaluation of summations of two- and three-point correlation functions with reasonable computational cost. In this paper, we show how a reformulation of nnn-point functions in the context of PEPS leads to extra contributions to the results that prove to play an important role. Benchmarks for the frustrated J1\ensuremath−J2{J}_{1}\ensuremath{-}{J}_{2}J1​\ensuremath−J2​ Heisenberg model illustrate the improved precision, efficiency, and stability of the simulations compared to previous approaches. Leveraging automatic differentiation to generate the most tedious and error-prone parts of the computation, the straightforward implementation presented here is a step towards broader adoption of the PEPS excitation ansatz in future applications.
Numerical treatment of two dimensional strongly-correlated systems is both extremely challenging and of fundamental importance. Infinite projected entangled-pair states (PEPS), a class of tensor networks, have demonstrated cutting-edge performance for ground state calculations, working directly in the thermodynamic limit. Furthermore, in recent years the application of PEPS has been extended to also low-lying excited states, using an ansatz that targets quasiparticle states above the ground state with high accuracy. A major technical challenge for those simulations is the accurate evaluation of summations of two- and three-point correlation functions with reasonable computational cost. In this work, we show how a reformulation of nnn-point functions in the context of PEPS leads to extra contributions to the results that prove to play an important role. Benchmarks for the frustrated J1−J2J_1-J_2J1​−J2​ Heisenberg model illustrate the improved precision, efficiency and stability of the simulations compared to previous approaches. Leveraging automatic differentiation to generate the most tedious and error-prone parts of the computation, the straightforward implementation presented here is a step towards broader adoption of the PEPS excitation ansatz in future applications.
The excitation ansatz for tensor networks is a powerful tool for simulating the low-lying quasiparticle excitations above ground states of strongly correlated quantum many-body systems. Recently, the two-dimensional tensor network class of infinite projected entangled-pair states gained new ground state optimization methods based on automatic differentiation, which are at the same time highly accurate and simple to implement. Naturally, the question arises whether these new ideas can also be used to optimize the excitation ansatz, which has recently been implemented in two dimensions as well. In this paper, we describe a straightforward way to reimplement the framework for excitations using automatic differentiation, and demonstrate its performance for the Hubbard model at half filling.
Whether two-dimensional tensor networks, such as the so-called projected entangled pair states (PEPS), have the potential to describe a genuine quantum critical point or quantum critical phases separating two ordered phases is of much interest for the understanding of frustrated magnets and strongly correlated systems. Here, the authors identify a specific family of SU(2)-invariant PEPS, which provides excellent variational energies for the J1−J2{J}_{1}-{J}_{2}J1​−J2​ frustrated Heisenberg model at large frustration. Using a full-update infinite-PEPS scheme directly in the thermodynamic limit, the authors provide evidence of diverging correlation lengths in both the singlet and triplet channels. It is argued that such PEPS give a qualitative description of the critical point or critical phase of the J1−J2{J}_{1}-{J}_{2}J1​−J2​ model.
We present an extension of a framework for simulating single quasiparticle or collective excitations on top of strongly correlated quantum many-body ground states using infinite projected entangled pair states, a tensor network ansatz for two-dimensional wave functions in the thermodynamic limit. Our approach performs a systematic summation of locally perturbed states in order to obtain excited eigenstates localized in momentum space, using the corner transfer matrix method, and generalizes the framework to arbitrary unit cell sizes, the implementation of global Abelian symmetries and fermionic systems. Results for several test cases are presented, including the transverse Ising model, the spin-1/21/21/2 Heisenberg model, and a free fermionic model, to demonstrate the capability of the method to accurately capture dispersions. We also provide insight into the nature of excitations at the $k=(\ensuremath{\pi},0)$ point of the Heisenberg model.
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.
We explain how to implement, in the context of projected entangled-pair states (PEPSs), the general procedure of fermionization of a tensor network introduced in P. Corboz and G. Vidal, Phys. Rev. B 80, 165129 (2009). The resulting fermionic PEPS, similar to previous proposals, can be used to study the ground state of interacting fermions on a two-dimensional lattice. As in the bosonic case, the cost of simulations depends on the amount of entanglement in the ground state and not directly on the strength of interactions. The present formulation of fermionic PEPS leads to a straightforward numerical implementation that allowed us to recycle much of the code for bosonic PEPS. We demonstrate that fermionic PEPS are a useful variational ansatz for interacting fermion systems by computing approximations to the ground state of several models on an infinite lattice. For a model of interacting spinless fermions, ground state energies lower than Hartree-Fock results are obtained, shifting the boundary between the metal and charge-density wave phases. For the $t\text{\ensuremath{-}}J$ model, energies comparable with those of a specialized Gutzwiller-projected ansatz are also obtained.