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:1311.7272 · Physical Review B

Fermionic projected entangled pair states at finite temperature

Piotr Czarnik, Jacek Dziarmaga

Published Jul 30, 2014• 33 citations
arXivHTMLPDF

Abstract

An algorithm for imaginary time evolution of a fermionic projected entangled pair state with ancillas from infinite temperature down to a finite temperature state is presented. As a benchmark application, it is applied to spinless fermions hopping on a square lattice subject to $p$-wave pairing interactions. With a tiny bias it allows one to evolve the system across a high-temperature continuous symmetry-breaking phase transition.

DOI
10.1103/physrevb.90.035144

Related papers

Projected entangled pair states at finite temperature: Imaginary time evolution with ancillas

Piotr Czarnik, Łukasz Cincio, Jacek Dziarmaga

A projected entangled pair state (PEPS) with ancillas is evolved in imaginary time. This tensor network represents a thermal state of a two-dimensional (2D) lattice quantum system. A finite-temperature phase diagram of the 2D quantum Ising model in a transverse field is obtained as a benchmark application.

arXiv:1209.0454•Physical Review B•Dec 5, 2012•100 citations

Projected entangled pair states at finite temperature: Iterative self-consistent bond renormalization for exact imaginary time evolution

Piotr Czarnik, Jacek Dziarmaga

A projected entangled pair state (PEPS) with ancillas can be evolved in imaginary time to obtain thermal states of a strongly correlated quantum system on a two-dimensional lattice. Every application of a Suzuki-Trotter gate multiplies the PEPS bond dimension $D$ by a factor $k$. It has to be renormalized back to the original $D$. In order to preserve the accuracy of the Suzuki-Trotter (ST) decomposition, the renormalization in principle has to take into account full environment made of the new tensors with the bond dimension $k\ifmmode\times\else\texttimes\fi{}D$. Here, we propose a self-consistent renormalization procedure operating with the original bond dimension $D$, but without compromising the accuracy of the ST decomposition. The iterative procedure renormalizes the bond using full environment made of renormalized tensors with the bond dimension $D$. After every renormalization, the new renormalized tensors are used to update the environment, and then the renormalization is repeated again and again until convergence. As a benchmark application, we obtain thermal states of the transverse field quantum Ising model on a square lattice, both infinite and finite, evolving the system across a second-order phase transition at finite temperature.

arXiv:1411.6778•Physical Review B•Jul 9, 2015•37 citations

Time evolution of an infinite projected entangled pair state: An efficient algorithm

Piotr Czarnik, Jacek Dziarmaga, Philippe Corboz

An infinite projected entangled pair state (iPEPS) is a tensor network ansatz to represent a quantum state on an infinite 2D lattice whose accuracy is controlled by the bond dimension $D$. Its real, Lindbladian, or imaginary time evolution can be split into small time steps. Every time step generates a new iPEPS with an enlarged bond dimension ${D}^{\ensuremath{'}}>D$, which is approximated by an iPEPS with the original $D$. In P. Czarnik and J. Dziarmaga, Phys. Rev. B 98, 045110 (2018), an algorithm was introduced to optimize the approximate iPEPS by maximizing directly its fidelity to the one with the enlarged bond dimension ${D}^{\ensuremath{'}}$. In this paper, we implement a more efficient optimization employing a local estimator of the fidelity. For imaginary time evolution of a thermal state's purification, we also consider using unitary disentangling gates acting on ancillas to reduce the required $D$. We test the algorithm simulating Lindbladian evolution and unitary evolution after a sudden quench of transverse field ${h}_{x}$ in the 2D quantum Ising model. Furthermore, we simulate thermal states of this model and estimate the critical temperature with good accuracy: $0.1%$ for ${h}_{x}=2.5$ and $0.5%$ for the more challenging case of ${h}_{x}=2.9$ close to the quantum critical point at ${h}_{x}=3.04438(2)$.

arXiv:1811.05497•Physical review. B./Physical review. B•Jan 8, 2019•125 citations

Simulation of fermionic lattice models in two dimensions with projected entangled-pair states: Next-nearest neighbor Hamiltonians

Philippe Corboz, Jacob Jordan, Guifré Vidal

In a recent contribution [P. Corboz, R. Or\'us, B. Bauer, and G. Vidal, Phys. Rev. B 81, 165104 (2010)] fermionic projected entangled-pair states (PEPSs) were used to approximate the ground state of free and interacting spinless fermion models, as well as the $t\text{\ensuremath{-}}J$ model. This paper revisits these three models in the presence of an additional next-nearest hopping amplitude in the Hamiltonian. First we explain how to account for next-nearest neighbor Hamiltonian terms in the context of fermionic PEPS algorithms based on simulating time evolution. Then we present benchmark calculations for the three models of fermions and compare our results against analytical, mean-field, and variational Monte Carlo results, respectively. Consistent with previous computations restricted to nearest neighbor Hamiltonians, we systematically obtain more accurate (or better converged) results for gapped phases than for gapless ones.

arXiv:1008.3937•Physical Review B•Dec 21, 2010•94 citations

Time evolution of an infinite projected entangled pair state: An algorithm from first principles

Piotr Czarnik, Jacek Dziarmaga

A typical quantum state obeying the area law for entanglement on an infinite two-dimensional (2D) lattice can be represented by a tensor network Ansatz, known as an infinite projected entangled pair state (iPEPS), with a finite bond dimension $D$. Its real-imaginary time evolution can be split into small time steps. An application of a time step generates a new iPEPS with a bond dimension $k$ times the original one. The new iPEPS does not make optimal use of its enlarged bond dimension $kD$; hence, in principle, it can be represented accurately by a more compact Ansatz, preferably with the original $D$. In this work we show how the more compact iPEPS can be optimized variationally to maximize its overlap with the new iPEPS. To compute the overlap we use the corner-transfer-matrix renormalization group. By simulating sudden quench of the transverse field in the 2D quantum Ising model with the proposed algorithm, we provide a proof of principle that real-time evolution can be simulated with iPEPS. A similar proof is provided with the same model for imaginary-time evolution of purification of its thermal states.

arXiv:1804.03872•Physical review. B./Physical review. B•Jul 6, 2018•22 citations

Symmetric projected entangled-pair states analysis of a phase transition in coupled spin- 12 ladders

Juraj Hašík, Glen Bigan Mbeng, Sylvain Capponi +2 more

Infinite projected entangled-pair states (iPEPS) have been introduced to accurately describe many-body wave functions on two-dimensional lattices. In this context, two aspects are crucial: the systematic improvement of the Ansatz by the optimization of its building blocks, i.e., tensors characterized by bond dimension $D$, and the extrapolation scheme to reach the ``thermodynamic'' limit $D\ensuremath{\rightarrow}\ensuremath{\infty}$. Recent advances in variational optimization and scaling based on correlation lengths demonstrated the ability of iPEPS to capture phases with spontaneously broken continuous symmetry such as the antiferromagnetic (N\'eel) one with high fidelity, in addition to valence-bond solids which are already well described by finite-$D$ iPEPS. In contrast, systems in the vicinity of continuous quantum phase transitions still present a challenge for iPEPS, especially when non-Abelian symmetries are involved. Here, we consider the iPEPS Ansatz to describe the continuous transition between the (gapless) antiferromagnet and the (gapped) paramagnet that exists in the $S=1/2$ Heisenberg model on coupled two-leg ladders. In particular, we show how accurate iPEPS results can be obtained down to a narrow interval around criticality and analyze the scaling of the order parameter in the N\'eel phase in a spatially anisotropic situation.

arXiv:2206.06230•Physical review. B./Physical review. B•Sep 30, 2022•3 citations
An algorithm for imaginary time evolution of a fermionic projected entangled pair state with ancillas from infinite temperature down to a finite temperature state is presented. As a benchmark application, it is applied to spinless fermions hopping on a square lattice subject to ppp-wave pairing interactions. With a tiny bias it allows one to evolve the system across a high-temperature continuous symmetry-breaking phase transition.
A projected entangled pair state (PEPS) with ancillas can be evolved in imaginary time to obtain thermal states of a strongly correlated quantum system on a two-dimensional lattice. Every application of a Suzuki-Trotter gate multiplies the PEPS bond dimension DDD by a factor kkk. It has to be renormalized back to the original DDD. In order to preserve the accuracy of the Suzuki-Trotter (ST) decomposition, the renormalization in principle has to take into account full environment made of the new tensors with the bond dimension $k\ifmmode\times\else\texttimes\fi{}D$. Here, we propose a self-consistent renormalization procedure operating with the original bond dimension DDD, but without compromising the accuracy of the ST decomposition. The iterative procedure renormalizes the bond using full environment made of renormalized tensors with the bond dimension DDD. After every renormalization, the new renormalized tensors are used to update the environment, and then the renormalization is repeated again and again until convergence. As a benchmark application, we obtain thermal states of the transverse field quantum Ising model on a square lattice, both infinite and finite, evolving the system across a second-order phase transition at finite temperature.
An infinite projected entangled pair state (iPEPS) is a tensor network ansatz to represent a quantum state on an infinite 2D lattice whose accuracy is controlled by the bond dimension DDD. Its real, Lindbladian, or imaginary time evolution can be split into small time steps. Every time step generates a new iPEPS with an enlarged bond dimension ${D}^{\ensuremath{'}}>D$, which is approximated by an iPEPS with the original DDD. In P. Czarnik and J. Dziarmaga, Phys. Rev. B 98, 045110 (2018), an algorithm was introduced to optimize the approximate iPEPS by maximizing directly its fidelity to the one with the enlarged bond dimension ${D}^{\ensuremath{'}}$. In this paper, we implement a more efficient optimization employing a local estimator of the fidelity. For imaginary time evolution of a thermal state's purification, we also consider using unitary disentangling gates acting on ancillas to reduce the required DDD. We test the algorithm simulating Lindbladian evolution and unitary evolution after a sudden quench of transverse field hx{h}_{x}hx​ in the 2D quantum Ising model. Furthermore, we simulate thermal states of this model and estimate the critical temperature with good accuracy: 0.10.1%0.1 for hx=2.5{h}_{x}=2.5hx​=2.5 and 0.50.5%0.5 for the more challenging case of hx=2.9{h}_{x}=2.9hx​=2.9 close to the quantum critical point at hx=3.04438(2){h}_{x}=3.04438(2)hx​=3.04438(2).
In a recent contribution [P. Corboz, R. Or\'us, B. Bauer, and G. Vidal, Phys. Rev. B 81, 165104 (2010)] fermionic projected entangled-pair states (PEPSs) were used to approximate the ground state of free and interacting spinless fermion models, as well as the $t\text{\ensuremath{-}}J$ model. This paper revisits these three models in the presence of an additional next-nearest hopping amplitude in the Hamiltonian. First we explain how to account for next-nearest neighbor Hamiltonian terms in the context of fermionic PEPS algorithms based on simulating time evolution. Then we present benchmark calculations for the three models of fermions and compare our results against analytical, mean-field, and variational Monte Carlo results, respectively. Consistent with previous computations restricted to nearest neighbor Hamiltonians, we systematically obtain more accurate (or better converged) results for gapped phases than for gapless ones.
A typical quantum state obeying the area law for entanglement on an infinite two-dimensional (2D) lattice can be represented by a tensor network Ansatz, known as an infinite projected entangled pair state (iPEPS), with a finite bond dimension DDD. Its real-imaginary time evolution can be split into small time steps. An application of a time step generates a new iPEPS with a bond dimension kkk times the original one. The new iPEPS does not make optimal use of its enlarged bond dimension kDkDkD; hence, in principle, it can be represented accurately by a more compact Ansatz, preferably with the original DDD. In this work we show how the more compact iPEPS can be optimized variationally to maximize its overlap with the new iPEPS. To compute the overlap we use the corner-transfer-matrix renormalization group. By simulating sudden quench of the transverse field in the 2D quantum Ising model with the proposed algorithm, we provide a proof of principle that real-time evolution can be simulated with iPEPS. A similar proof is provided with the same model for imaginary-time evolution of purification of its thermal states.
Infinite projected entangled-pair states (iPEPS) have been introduced to accurately describe many-body wave functions on two-dimensional lattices. In this context, two aspects are crucial: the systematic improvement of the Ansatz by the optimization of its building blocks, i.e., tensors characterized by bond dimension DDD, and the extrapolation scheme to reach the ``thermodynamic'' limit $D\ensuremath{\rightarrow}\ensuremath{\infty}$. Recent advances in variational optimization and scaling based on correlation lengths demonstrated the ability of iPEPS to capture phases with spontaneously broken continuous symmetry such as the antiferromagnetic (N\'eel) one with high fidelity, in addition to valence-bond solids which are already well described by finite-DDD iPEPS. In contrast, systems in the vicinity of continuous quantum phase transitions still present a challenge for iPEPS, especially when non-Abelian symmetries are involved. Here, we consider the iPEPS Ansatz to describe the continuous transition between the (gapless) antiferromagnet and the (gapped) paramagnet that exists in the S=1/2S=1/2S=1/2 Heisenberg model on coupled two-leg ladders. In particular, we show how accurate iPEPS results can be obtained down to a narrow interval around criticality and analyze the scaling of the order parameter in the N\'eel phase in a spatially anisotropic situation.