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

Infinite projected entangled pair states algorithm improved: Fast full update and gauge fixing

Ho N. Phien, Johann A. Bengua, Hoang Duong Tuan, Philippe Corboz, Román Orús

Published Jul 24, 2015• 226 citations
arXivHTMLPDF

Abstract

The infinite projected entangled pair states (iPEPS) algorithm [J. Jordan et al., Phys. Rev. Lett. 101, 250602 (2008)] has become a useful tool in the calculation of ground-state properties of two-dimensional quantum lattice systems in the thermodynamic limit. Despite its many successful implementations, the method has some limitations in its present formulation which hinder its application to some highly entangled systems. The purpose of this paper is to unravel some of these issues, in turn enhancing the stability and efficiency of iPEPS methods. For this, we first introduce the fast full update scheme, where effective environment and iPEPS tensors are both simultaneously updated (or evolved) throughout time. As we shall show, this implies two crucial advantages: (i) dramatic computational savings and (ii) improved overall stability. In addition, we extend the application of the local gauge fixing, successfully implemented for finite-size PEPS [M. Lubasch et al., Phys. Rev. B 90, 064425 (2014)], to the iPEPS algorithm. We see that the gauge fixing not only further improves the stability of the method but also accelerates the convergence of the alternating least-squares sweeping in the (either ``full'' or ``fast full'') tensor update scheme. The improvement in terms of computational cost and stability of the resulting ``improved'' iPEPS algorithm is benchmarked by studying the ground-state properties of the quantum Heisenberg and transverse-field Ising models on an infinite square lattice.

DOI
10.1103/physrevb.92.035142

Related papers

Optimization of infinite projected entangled pair states: The role of multiplets and their breaking

Juraj Hašík, Federico Becca

The infinite projected entangled pair states (iPEPS) technique [J. Jordan et al., Phys. Rev. Lett. 101, 250602 (2008)] has been widely used in the recent years to assess the properties of two-dimensional quantum systems, working directly in the thermodynamic limit. This formalism, which is based upon a tensor-network representation of the ground-state wave function, has several appealing features, e.g., encoding the so-called area law of entanglement entropy by construction; still, the method presents critical issues when dealing with the optimization of tensors in order to find the best possible approximation to the exact ground state of a given Hamiltonian. Here, we discuss the obstacles that arise in the optimization by imaginary-time evolution within the so-called simple and full updates and connect them to the emergence of a sharp multiplet structure in the ``virtual'' indices of tensors. In this case, a generic choice of the bond dimension $D$ is not compatible with the multiplets and leads to a symmetry breaking (e.g., generating a finite magnetic order). In addition, varying the initial guess, different final states may be reached, with very large deviations in the magnetization value. In order to exemplify this behavior, we show the results of the $S=1/2$ Heisenberg model on an array of coupled ladders, for which a vanishing magnetization below the critical interladder coupling is recovered only for selected values of $D$, while a blind optimization with a generic $D$ gives rise to a finite magnetization down to the limit of decoupled ladders.

arXiv:1905.02164•Physical review. B./Physical review. B•Aug 21, 2019•6 citations

Towards next-generation methods to optimize two-dimensional tensor networks: Algorithmic differentiation and applications to quantum magnets

Juraj Hašík

Infinite Projected entangled-pair state (iPEPS) is a two-dimensional tensor network (TN) ansatz targeting the ground states of lattice models directly in the thermodynamic limit. It is the most successful among the true two-dimensional TN methods, challenging the more common density matrix renormalization group computations conducted on a finite-width cylinders. iPEPS, being a variational method does not suffer from the infamous sign problem, and hence could be readily applied to both fermionic problems off the half-filling and frustrated spin system. Using the imaginary-time optimization, the so called Full Update, led to the celebrated results of iPEPS such as the analysis of stripes in doped t-J and Hubbard models or the insight into the magnetization plateaus in Shastry-Sutherland model. However, the physics of the continuous phase transitions from the symmetry-broken Neel phase into the valence-bond solids or spin-liquids in frustrated magnets is still eluding iPEPS. The main difficulty lies in the optimization, since Full Update is just not precise enough. In the pursuit of the best variational iPEPS states, the gradient methods were put forward in 2016 by Corboz and Vanderstraeten. However, due to their complex nature, they have not been widely adopted especially for iPEPS ansatze with larger unit cells. We will introduce an alternative approach first applied in the context of iPEPS by Liao et al., based on the well established method of Algorithmic differentiation (AD) which allows us to evaluate gradients for iPEPS in a robust and conceptually simple way. Finally, we will present the results of our investigation into the physics of frustrated magnets putting AD to work.

•Oct 25, 2019•3 citations

Spectral functions with infinite projected entangled-pair states

Juan Diego Arias Espinoza, Philippe Corboz

Infinite projected entangled-pair states (iPEPS) provide a powerful tool to study two-dimensional strongly correlated systems directly in the thermodynamic limit. In this work, we extend the iPEPS toolbox by a method to efficiently evaluate nonequal time two-point correlators, enabling the computation of spectral functions. It is based on an iPEPS ansatz of the ground state in a large unit cell, with an operator applied in the center of the cell, which is evolved in real-time using the fast-full update method. At every time step, the two-point correlators within a cell are computed based on the corner transfer matrix renormalization group method. Benchmark results for the two-dimensional transverse field Ising model show that the main features of the dynamical structure factor can already be reproduced at relatively small bond dimensions and unit cell sizes. The results for the magnon dispersion are found to be in good agreement with previous data obtained with the iPEPS excitation ansatz.

arXiv:2405.10628•Physical review. B./Physical review. B•Sep 26, 2024•9 citations

Variational optimization with infinite projected entangled-pair states

Philippe Corboz

An infinite projected entangled-pair state (iPEPS) is a powerful variational tensor network ansatz for two-dimensional ground states in the thermodynamic limit, and can be seen as a natural generalization of a matrix-product state to two dimensions. One of the main challenges in iPEPS simulations is the optimization of the tensors, i.e., finding the optimal variational parameters, in order to have the best representation of the ground state of a given Hamiltonian. The author presents a variational optimization scheme, in which the energy is minimized in an iterative way by sweeping over all the tensors in the ansatz, in a similar spirit as done in the density-matrix renormalization group method. Benchmark results for challenging problems are presented that show that the variational scheme yields considerably more accurate results than the previously best imaginary-time evolution algorithm, with a similar computational cost and with a faster convergence towards the ground state.

arXiv:1605.03006•Physical review. B./Physical review. B•Jul 14, 2016•234 citations

Time evolution of an infinite projected entangled pair state: A gradient tensor update in the tangent space

Jacek Dziarmaga

Time evolution of an infinite two-dimensional (2D) many body quantum lattice system can be described by the Suzuki-Trotter decomposition applied to the infinite projected entangled pair state (iPEPS). Each Trotter gate increases the bond dimension of the tensor network, $D$, that has to be truncated back in a way that minimizes a suitable error measure. This paper goes beyond simplified error measures--like the one used in the full update, the simple update, and their intermediate neighborhood tensor update--and directly maximizes an overlap between the exact iPEPS with the increased bond dimension and the new iPEPS with the truncated one. The optimization is performed in a tangent space of the iPEPS variational manifold. This gradient tensor update is benchmarked by a simulation of a sudden quench of a transverse field in the 2D quantum Ising model and the quantum Kibble-Zurek mechanism in the same 2D system.

arXiv:2205.11067•Physical review. B./Physical review. B•Jul 11, 2022•19 citations

Infinite projected entangled-pair state algorithm for ruby and triangle-honeycomb lattices

Saeed S. Jahromi, Román Orús, Mehdi Kargarian +1 more

The infinite projected entangled-pair state (iPEPS) algorithm is one of the most efficient techniques for studying the ground-state properties of two-dimensional quantum lattice Hamiltonians in the thermodynamic limit. Here, we show how the algorithm can be adapted to explore nearest-neighbor local Hamiltonians on the ruby and triangle-honeycomb lattices, using the corner transfer matrix (CTM) renormalization group for 2D tensor network contraction. Additionally, we show how the CTM method can be used to calculate the ground-state fidelity per lattice site and the boundary density operator and entanglement entropy (EE) on an infinite cylinder. As a benchmark, we apply the iPEPS method to the ruby model with anisotropic interactions and explore the ground-state properties of the system. We further extract the phase diagram of the model in different regimes of the couplings by measuring two-point correlators, ground-state fidelity, and EE on an infinite cylinder. Our phase diagram is in agreement with previous studies of the model by exact diagonalization.

arXiv:1711.04798•Physical review. B./Physical review. B•Mar 29, 2018•32 citations
The infinite projected entangled pair states (iPEPS) algorithm [J. Jordan et al., Phys. Rev. Lett. 101, 250602 (2008)] has become a useful tool in the calculation of ground-state properties of two-dimensional quantum lattice systems in the thermodynamic limit. Despite its many successful implementations, the method has some limitations in its present formulation which hinder its application to some highly entangled systems. The purpose of this paper is to unravel some of these issues, in turn enhancing the stability and efficiency of iPEPS methods. For this, we first introduce the fast full update scheme, where effective environment and iPEPS tensors are both simultaneously updated (or evolved) throughout time. As we shall show, this implies two crucial advantages: (i) dramatic computational savings and (ii) improved overall stability. In addition, we extend the application of the local gauge fixing, successfully implemented for finite-size PEPS [M. Lubasch et al., Phys. Rev. B 90, 064425 (2014)], to the iPEPS algorithm. We see that the gauge fixing not only further improves the stability of the method but also accelerates the convergence of the alternating least-squares sweeping in the (either full'' or fast full'') tensor update scheme. The improvement in terms of computational cost and stability of the resulting ``improved'' iPEPS algorithm is benchmarked by studying the ground-state properties of the quantum Heisenberg and transverse-field Ising models on an infinite square lattice.
The infinite projected entangled pair states (iPEPS) technique [J. Jordan et al., Phys. Rev. Lett. 101, 250602 (2008)] has been widely used in the recent years to assess the properties of two-dimensional quantum systems, working directly in the thermodynamic limit. This formalism, which is based upon a tensor-network representation of the ground-state wave function, has several appealing features, e.g., encoding the so-called area law of entanglement entropy by construction; still, the method presents critical issues when dealing with the optimization of tensors in order to find the best possible approximation to the exact ground state of a given Hamiltonian. Here, we discuss the obstacles that arise in the optimization by imaginary-time evolution within the so-called simple and full updates and connect them to the emergence of a sharp multiplet structure in the ``virtual'' indices of tensors. In this case, a generic choice of the bond dimension DDD is not compatible with the multiplets and leads to a symmetry breaking (e.g., generating a finite magnetic order). In addition, varying the initial guess, different final states may be reached, with very large deviations in the magnetization value. In order to exemplify this behavior, we show the results of the S=1/2S=1/2S=1/2 Heisenberg model on an array of coupled ladders, for which a vanishing magnetization below the critical interladder coupling is recovered only for selected values of DDD, while a blind optimization with a generic DDD gives rise to a finite magnetization down to the limit of decoupled ladders.
Infinite Projected entangled-pair state (iPEPS) is a two-dimensional tensor network (TN) ansatz targeting the ground states of lattice models directly in the thermodynamic limit. It is the most successful among the true two-dimensional TN methods, challenging the more common density matrix renormalization group computations conducted on a finite-width cylinders. iPEPS, being a variational method does not suffer from the infamous sign problem, and hence could be readily applied to both fermionic problems off the half-filling and frustrated spin system. Using the imaginary-time optimization, the so called Full Update, led to the celebrated results of iPEPS such as the analysis of stripes in doped t-J and Hubbard models or the insight into the magnetization plateaus in Shastry-Sutherland model. However, the physics of the continuous phase transitions from the symmetry-broken Neel phase into the valence-bond solids or spin-liquids in frustrated magnets is still eluding iPEPS. The main difficulty lies in the optimization, since Full Update is just not precise enough. In the pursuit of the best variational iPEPS states, the gradient methods were put forward in 2016 by Corboz and Vanderstraeten. However, due to their complex nature, they have not been widely adopted especially for iPEPS ansatze with larger unit cells. We will introduce an alternative approach first applied in the context of iPEPS by Liao et al., based on the well established method of Algorithmic differentiation (AD) which allows us to evaluate gradients for iPEPS in a robust and conceptually simple way. Finally, we will present the results of our investigation into the physics of frustrated magnets putting AD to work.
Infinite projected entangled-pair states (iPEPS) provide a powerful tool to study two-dimensional strongly correlated systems directly in the thermodynamic limit. In this work, we extend the iPEPS toolbox by a method to efficiently evaluate nonequal time two-point correlators, enabling the computation of spectral functions. It is based on an iPEPS ansatz of the ground state in a large unit cell, with an operator applied in the center of the cell, which is evolved in real-time using the fast-full update method. At every time step, the two-point correlators within a cell are computed based on the corner transfer matrix renormalization group method. Benchmark results for the two-dimensional transverse field Ising model show that the main features of the dynamical structure factor can already be reproduced at relatively small bond dimensions and unit cell sizes. The results for the magnon dispersion are found to be in good agreement with previous data obtained with the iPEPS excitation ansatz.
An infinite projected entangled-pair state (iPEPS) is a powerful variational tensor network ansatz for two-dimensional ground states in the thermodynamic limit, and can be seen as a natural generalization of a matrix-product state to two dimensions. One of the main challenges in iPEPS simulations is the optimization of the tensors, i.e., finding the optimal variational parameters, in order to have the best representation of the ground state of a given Hamiltonian. The author presents a variational optimization scheme, in which the energy is minimized in an iterative way by sweeping over all the tensors in the ansatz, in a similar spirit as done in the density-matrix renormalization group method. Benchmark results for challenging problems are presented that show that the variational scheme yields considerably more accurate results than the previously best imaginary-time evolution algorithm, with a similar computational cost and with a faster convergence towards the ground state.
Time evolution of an infinite two-dimensional (2D) many body quantum lattice system can be described by the Suzuki-Trotter decomposition applied to the infinite projected entangled pair state (iPEPS). Each Trotter gate increases the bond dimension of the tensor network, DDD, that has to be truncated back in a way that minimizes a suitable error measure. This paper goes beyond simplified error measures--like the one used in the full update, the simple update, and their intermediate neighborhood tensor update--and directly maximizes an overlap between the exact iPEPS with the increased bond dimension and the new iPEPS with the truncated one. The optimization is performed in a tangent space of the iPEPS variational manifold. This gradient tensor update is benchmarked by a simulation of a sudden quench of a transverse field in the 2D quantum Ising model and the quantum Kibble-Zurek mechanism in the same 2D system.
The infinite projected entangled-pair state (iPEPS) algorithm is one of the most efficient techniques for studying the ground-state properties of two-dimensional quantum lattice Hamiltonians in the thermodynamic limit. Here, we show how the algorithm can be adapted to explore nearest-neighbor local Hamiltonians on the ruby and triangle-honeycomb lattices, using the corner transfer matrix (CTM) renormalization group for 2D tensor network contraction. Additionally, we show how the CTM method can be used to calculate the ground-state fidelity per lattice site and the boundary density operator and entanglement entropy (EE) on an infinite cylinder. As a benchmark, we apply the iPEPS method to the ruby model with anisotropic interactions and explore the ground-state properties of the system. We further extract the phase diagram of the model in different regimes of the couplings by measuring two-point correlators, ground-state fidelity, and EE on an infinite cylinder. Our phase diagram is in agreement with previous studies of the model by exact diagonalization.