Projected Entangled Pair States Research Articles

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

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.

Chiral topological spin liquids with projected entangled pair states

Topological chiral phases are ubiquitous in the physics of the Fractional Quantum Hall Effect. Non-chiral topological spin liquids are also well known. Here, using the framework of projected entangled pair states (PEPS), we construct a family of chiral spin liquids on the square lattice which are generalized spin-1/2 Resonating Valence Bond (RVB) states obtained from deformed local tensors with $d+i\, d$ symmetry. On a cylinder, we construct four topological sectors with even or odd number of spinons on the boundary and even or odd number of ($\mathbb{Z}_2$) fluxes penetrating the cylinder which, we argue, remain orthogonal in the limit of infinite perimeter. The analysis of the transfer matrix provides evidence of short-range (long-range) triplet (singlet) correlations as for the critical (non-chiral) RVB state. The Entanglement Spectrum exhibits chiral edge modes, which we confront to predictions of Conformal Field Theory, and the corresponding Entanglement Hamiltonian is shown to be long ranged.

Continuum tensor network field states, path integral representations and spatial symmetries

A natural way to generalize tensor network variational classes to quantum field systems is via a continuous tensor contraction. This approach is first illustrated for the class of quantum field states known as continuous matrix-product states (cMPS). As a simple example of the path-integral representation we show that the state of a dynamically evolving quantum field admits a natural representation as a cMPS. A completeness argument is also provided that shows that all states in Fock space admit a cMPS representation when the number of variational parameters tends to infinity. Beyond this, we obtain a well-behaved field limit of projected entangled-pair states (PEPS) in two dimensions that provide an abstract class of quantum field states with natural symmetries. We demonstrate how symmetries of the physical field state are encoded within the dynamics of an auxiliary field system of one dimension less. In particular, the imposition of Euclidean symmetries on the physical system requires that the auxiliary system involved in the class’ definition must be Lorentz-invariant. The physical field states automatically inherit entropy area laws from the PEPS class, and are fully described by the dissipative dynamics of a lower dimensional virtual field system. Our results lie at the intersection many-body physics, quantum field theory and quantum information theory, and facilitate future exchanges of ideas and insights between these disciplines.

Fast convergence of imaginary time evolution tensor network algorithms by recycling the environment

We propose an environment recycling scheme to speed up a class of tensor network algorithms that produce an approximation to the ground state of a local Hamiltonian by simulating an evolution in imaginary time. Specifically, we consider the time-evolving block decimation (TEBD) algorithm applied to infinite systems in 1D and 2D, where the ground state is encoded, respectively, in a matrix product state (MPS) and in a projected entangled-pair state (PEPS). An important ingredient of the TEBD algorithm (and a main computational bottleneck, especially with PEPS in 2D) is the computation of the so-called environment, which is used to determine how to optimally truncate the bond indices of the tensor network so that their dimension is kept constant. In current algorithms, the environment is computed at each step of the imaginary time evolution, to account for the changes that the time evolution introduces in the many-body state represented by the tensor network. Our key insight is that close to convergence, most of the changes in the environment are due to a change in the choice of gauge in the bond indices of the tensor network, and not in the many-body state. Indeed, a consistent choice of gauge in the bond indices confirms that the environment is essentially the same over many time steps and can thus be re-used, leading to very substantial computational savings. We demonstrate the resulting approach in 1D and 2D by computing the ground state of the quantum Ising model in a transverse magnetic field.

Chiral projected entangled-pair state with topological order.

We show that projected entangled-pair states (PEPS) can describe chiral topologically ordered phases. For that, we construct a simple PEPS for spin-1/2 particles in a two-dimensional lattice. We reveal a symmetry in the local projector of the PEPS that gives rise to the global topological character. We also extract characteristic quantities of the edge conformal field theory using the bulk-boundary correspondence.

Numerical detection of symmetry-enriched topological phases with space-group symmetry

Topologically ordered phases of matter, in particular so-called symmetry enriched topological (SET) phases, can exhibit quantum number fractionalization in the presence of global symmetry. In Z_2 topologically ordered states in two dimensions, fundamental translations T_x and T_y acting on anyons can either commute or anticommute. This property, crystal momentum fractionalization, can be seen in a periodicity of the excited-state spectrum in the Brillouin zone. We present a numerical method to detect the presence of this form of symmetry enrichment given a projected entangled pair state (PEPS); we study the minima of spectrum of correlation lengths of the transfer matrix for a cylinder. As a benchmark, we demonstrate our method using a modified toric code model with perturbation. An enhanced periodicity in momentum clearly reveals the anticommutation relation {T_x,T_y}=0$ for the corresponding quasiparticles in the system.

Gauging Quantum States: From Global to Local Symmetries in Many-Body Systems

We present an operational procedure to transform global symmetries into local symmetries at the level of individual quantum states, as opposed to typical gauging prescriptions for Hamiltonians or Lagrangians. We then construct a compatible gauging map for operators, which preserves locality and reproduces the minimal coupling scheme for simple operators. By combining this construction with the formalism of projected entangled-pair states (PEPS), we can show that an injective PEPS for the matter fields is gauged into a G-injective PEPS for the combined gauge-matter system, which potentially has topological order. We derive the corresponding parent Hamiltonian, which is a frustration free gauge theory Hamiltonian closely related to the Kogut-Susskind Hamiltonian at zero coupling constant. We can then introduce gauge dynamics at finite values of the coupling constant by applying a local filtering operation. This scheme results in a low-parameter family of gauge invariant states of which we can accurately probe the phase diagram, as we illustrate by studying a Z2 gauge theory with Higgs matter.

Approximating Gibbs states of local Hamiltonians efficiently with projected entangled pair states

We analyze the error of approximating Gibbs states of local quantum spin Hamiltonians on lattices with Projected Entangled Pair States (PEPS) as a function of the bond dimension ($D$), temperature ($\beta^{-1}$), and system size ($N$). First, we introduce a compression method in which the bond dimension scales as $D=e^{O(\log^2(N/\epsilon))}$ if $\beta<O(\log (N))$. Second, building on the work of Hastings [Phys. Rev. B 73, 085115 (2006)], we derive a polynomial scaling relation, $D=\left(N/\epsilon\right)^{O(\beta)}$. This implies that the manifold of PEPS forms an efficient representation of Gibbs states of local quantum Hamiltonians. From those bounds it also follows that ground states can be approximated with $D=N^{O(\log(N))}$ whenever the density of states only grows polynomially in the system size. All results hold for any spatial dimension of the lattice.

Striped critical spin liquid in a spin-orbital entangled RVB state in a projected entangled-pair state representation

We introduce a spin-orbital entangled (SOE) resonating valence bond (RVB) state on a square lattice of spins $\frac{1}{2}$ and orbitals represented by pseudospins $\frac{1}{2}$. Like the standard RVB state, it is a superposition of nearest-neighbor hard-core coverings of the lattice by spin singlets, but adjacent singlets are favored to have perpendicular orientations and, more importantly, an orientation of each singlet is entangled with orbitals state on its two lattice sites. The SOE-RVB state can be represented by a projected entangled pair state (PEPS) with a bond dimension $D=4$. This representation helps to reveal that the state is a superposition of striped coverings conserving a topological quantum number. The stripes are a critical quantum spin liquid. We propose a spin-orbital Hamiltonian supporting a SOE-RVB ground state.

Perturbative 2-body parent Hamiltonians for projected entangled pair states

We construct parent Hamiltonians involving only local 2-body interactions for a broad class of projected entangled pair states (PEPS). Making use of perturbation gadget techniques, we define a perturbative Hamiltonian acting on the virtual PEPS space with a finite order low energy effective Hamiltonian that is a gapped, frustration-free parent Hamiltonian for an encoded version of a desired PEPS. For topologically ordered PEPS, the ground space of the low energy effective Hamiltonian is shown to be in the same phase as the desired state to all orders of perturbation theory. An encoded parent Hamiltonian for the double semion string net ground state is explicitly constructed as a concrete example.

Advances on tensor network theory: symmetries, fermions, entanglement, and holography

This is a short review on selected theory developments on Tensor Network (TN) states for strongly correlated systems. Specifically, we briefly review the effect of symmetries in TN states, fermionic TNs, the calculation of entanglement Hamiltonians from Projected Entangled Pair States (PEPS), and the relation between the Multi-scale Entanglement Renormalization Ansatz (MERA) and the AdS/CFT or gauge/gravity duality. We stress the role played by entanglement in the emergence of several physical properties and objects through the TN language. Some recent results along these lines are also discussed.

Entanglement Hamiltonian of the quantum Néel state

2D projected entangled pair states (PEPS) provide a unique framework giving access to detailed entanglement features of correlated (spin or electronic) systems. For a bi-partitioned quantum system, it has been argued that the entanglement spectrum (ES) is in a one-to-one correspondence with the physical edge spectrum on the cut and that the structure of the corresponding entanglement Hamiltonian (EH) reflects closely bulk properties (finite correlation length, criticality, topological order, etc). However, entanglement properties of systems with spontaneously broken continuous symmetry are still not fully understood. The spin-1/2 square lattice Heisenberg antiferromagnet provides a simple example showing spontaneous breaking of SU(2) symmetry down to U(1). The ground state can be viewed as a ‘quantum Néel state’ where the classical (Néel) staggered magnetization is reduced by quantum fluctuations. Here I consider the (critical) resonating valence bond (RVB) state doped with spinons to describe such a state; this enables the use of the associated PEPS representation (with virtual bond dimension D = 3) to compute the EH and the ES for a partition of an (infinite) cylinder. In particular, I find that the EH is (almost exactly) a chain of a dilute mixture of heavy (↓ spins) and light (↑ spins) hardcore bosons, where light particles are subject to long-range hoppings. The corresponding ES shows drastic differences with the typical ES obtained previously for ground states with restored SU(2)-symmetry (on finite systems).

Twisted injectivity in projected entangled pair states and the classification of quantum phases

We introduce a class of projected entangled pair states (PEPS) which is based on a group symmetry twisted by a 3-cocycle of the group. This twisted symmetry is expressed as a matrix product operator (MPO) with bond dimension greater than 1 and acts on the virtual boundary of a PEPS tensor. We show that it gives rise to a new standard form for PEPS from which we construct a family of local Hamiltonians which are gapped, frustration-free and include fixed points of the renormalization group flow. Based on this insight, we advance the classification of 2D gapped quantum spin systems by showing how this new standard form for PEPS determines the emergent topological order of these local Hamiltonians. Specifically, we identify their universality class as Dijkgraaf–Witten topological quantum field theory (TQFT).

Symmetries and boundary theories for chiral projected entangled pair states

We investigate the topological character of lattice chiral Gaussian fermionic states in two dimensions possessing the simplest descriptions in terms of projected entangled pair states (PEPS). They are ground states of two different kinds of Hamiltonians. The first one, ${\mathcal{H}}_{\mathrm{ff}}$, is local, frustration free, and gapless. It can be interpreted as describing a quantum phase transition between different topological phases. The second one, ${\mathcal{H}}_{\mathrm{fb}}$, is gapped, and has hopping terms scaling as $1/{r}^{3}$ with the distance $r$. The gap is robust against local perturbations, which allows us to define a Chern number for the PEPS. As for (nonchiral) topological PEPS, the nontrivial topological properties can be traced down to the existence of a symmetry in the virtual modes that are used to build the state. Based on that symmetry, we construct stringlike operators acting on the virtual modes that can be continuously deformed without changing the state. On the torus, the symmetry implies that the ground state space of the local parent Hamiltonian is twofold degenerate. By adding a string wrapping around the torus, one can change one of the ground states into the other. We use the special properties of PEPS to build the boundary theory and show how the symmetry results in the appearance of chiral modes, and a universal correction to the area law for the zero R\'enyi entropy.

Semionic resonating valence-bond states

The nature of the kagome Heisenberg antiferromagnet (HAFM) is under ongoing debate. While recent evidence points towards a Z_2 topological spin liquid, the exact nature of the topological phase is still unclear. In this paper, we introduce semionic Resonating Valence Bond (RVB) states, this is, Resonating Valence Bond states which are in the Z_2 ordered double-semion phase, and study them using Projected Entangled Pair States (PEPS). We investigate their physics and study their suitability as an ansatz for the HAFM, as compared to a conventional RVB state which is in the Toric Code Z_2 topological phase. In particular, we find that a suitably optimized "semionic simplex RVB" outperforms the equally optimized conventional "simplex RVB" state, and that the entanglement spectrum (ES) of the semionic RVB behaves very differently from the ES of the conventional RVB, which suggests to use the ES to discriminate the two phases. Finally, we also discuss the possible relevance of space group symmetry breaking in valence bond wavefunctions with double-semion topological order.

Algorithms for finite projected entangled pair states

Projected Entangled Pair States (PEPS) are a promising ansatz for the study of strongly correlated quantum many-body systems in two dimensions. But due to their high computational cost, developing and improving PEPS algorithms is necessary to make the ansatz widely usable in practice. Here we analyze several algorithmic aspects of the method. On the one hand, we quantify the connection between the correlation length of the PEPS and the accuracy of its approximate contraction, and discuss how purifications can be used in the latter. On the other, we present algorithmic improvements for the update of the tensor that introduce drastic gains in the numerical conditioning and the efficiency of the algorithms. Finally, the state-of-the-art general PEPS code is benchmarked with the Heisenberg and quantum Ising models on lattices of up to $21 \times 21$ sites.

Uncontrolled disorder effects in fabricating photonic quantum simulators on a kagome geometry: A projected-entangled-pair-state versus exact-diagonalization analysis

We propose a flexible numerical framework for extracting the energy spectra and photon transfer dynamics of a unit kagome cell with disordered cavity-cavity couplings under realistic experimental conditions. A projected-entangled pair state (PEPS) ansatz to the many-photon wavefunction allows to gain a detailed understanding of the effects of undesirable disorder in fabricating well-controlled and scalable photonic quantum simulators. The correlation functions associated with the propagation of two-photon excitations reveal intriguing interference patterns peculiar to the kagome geometry and promise at the same time a highly tunable quantum interferometry device with a signature for the formation of resonant or Fabry-Pe\'rot-like transmission of photons. Our results justify the use of the proposed PEPS technique for addressing the role of disorder in such quantum simulators in the microwave regime and promises a sophisticated numerical machinery for yet further explorations of the scalability of the resulting kagome arrays. The introduced methodology and the physical results may also pave the way for unraveling exotic phases of correlated light on a kagome geometry.

Topology and criticality in the resonating Affleck-Kennedy-Lieb-Tasaki loop spin liquid states

We exploit a natural projected entangled-pair state (PEPS) representation for the resonating Affleck-Kennedy-Lieb-Tasaki loop (RAL) state. By taking advantage of PEPS-based analytical and numerical methods, we characterize the RAL states on various two-dimensional lattices. On square and honeycomb lattices, these states are critical since the dimer-dimer correlations decay as a power law. On the kagome lattice, the RAL state has exponentially decaying correlation functions, supporting the scenario of a gapped spin liquid. We provide further evidence that the RAL state on the kagome lattice is a ${\mathbb{Z}}_{2}$ spin liquid, by identifying the four topological sectors and computing the topological entropy. Furthermore, we construct a one-parameter family of PEPS states interpolating between the RAL state and a short-range resonating valence bond state and find a critical point, consistent with the fact that the two states belong to two different phases. We also perform a variational study of the spin-1 kagome Heisenberg model using this one-parameter PEPS.

Unifying projected entangled pair state contractions

The approximate contraction of a tensor network of projected entangled pair states (PEPS) is a fundamental ingredient of any PEPS algorithm, required for the optimization of the tensors in ground state search or time evolution, as well as for the evaluation of expectation values. An exact contraction is in general impossible, and the choice of the approximating procedure determines the efficiency and accuracy of the algorithm. We analyze different previous proposals for this approximation, and show that they can be understood via the form of their environment, i.e. the operator that results from contracting part of the network. This provides physical insight into the limitation of various approaches, and allows us to introduce a new strategy, based on the idea of clusters, that unifies previous methods. The resulting contraction algorithm interpolates naturally between the cheapest and most imprecise and the most costly and most precise method. We benchmark the different algorithms with finite PEPS, and show how the cluster strategy can be used for both the tensor optimization and the calculation of expectation values. Additionally, we discuss its applicability to the parallelization of PEPS and to infinite systems.

Tensor Renormalization of Quantum Many-Body Systems Using Projected Entangled Simplex States

We propose a new class of tensor-network states, which we name projected entangled simplex states (PESS), for studying the ground-state properties of quantum lattice models. These states extend the pair-correlation basis of projected entangled pair states (PEPS) to a simplex. PESS are an exact representation of the simplex solid states and provide an efficient trial wave function that satisfies the area law of entanglement entropy. We introduce a simple update method for evaluating the PESS wave function based on imaginary-time evolution and the higher-order singular-value decomposition of tensors. By applying this method to the spin-1/2 antiferromagnetic Heisenberg model on the kagome lattice, we obtain accurate and systematic results for the ground-state energy, which approach the lowest upper bounds yet estimated for this quantity.