Source author record

Frank Verstraete

Frank Verstraete appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

quant-ph cond-mat.str-el cond-mat.stat-mech hep-lat hep-th math-ph math.MP cond-mat.other Computational Complexity physics.chem-ph Computer Vision cond-mat.dis-nn Machine Learning math.QA math.ST physics.comp-ph Statistics Theory

Catalog footprint

What is connected

88works

17topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2023arXiv

Contrasting pseudo-criticality in the classical two-dimensional Heisenberg and $\mathrm{RP}^2$ models: zero-temperature phase transition versus finite-temperature crossover

Tensor-network methods are used to perform a comparative study of the two-dimensional classical Heisenberg and $\mathrm{RP}^2$ models. We demonstrate that uniform matrix product states (MPS) with explicit $\mathrm{SO}(3)$ symmetry can probe correlation lengths up to $\mathcal{O}(10^3)$ sites accurately, and we study the scaling of entanglement entropy and universal features of MPS entanglement spectra. For the Heisenberg model, we find no signs of a finite-temperature phase transition, supporting the scenario of asymptotic freedom. For the $\mathrm{RP}^2$ model we observe an abrupt onset of scaling behaviour, consistent with hints of a finite-temperature phase transition reported in previous studies. A careful analysis of the softening of the correlation length divergence, the scaling of the entanglement entropy and the MPS entanglement spectra shows that our results are inconsistent with true criticality, but are rather in agreement with the scenario of a crossover to a pseudo-critical region which exhibits strong signatures of nematic quasi-long-range order at length scales below the true correlation length. Our results reveal a fundamental difference in scaling behaviour between the Heisenberg and $\mathrm{RP}^2$ models: Whereas the emergence of scaling in the former shifts to zero temperature if the bond dimension is increased, it occurs at a finite bond-dimension independent crossover temperature in the latter.

preprint2022arXiv

Direct sampling of projected entangled-pair states

Variational Monte Carlo studies employing projected entangled-pair states (PEPS) have recently shown that they can provide answers on long-standing questions such as the nature of the phases in the two-dimensional $J_1 - 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.

preprint2022arXiv

Efficient MPS methods for extracting spectral information on rings and cylinders

Based on the MPS formalism, we introduce an ansatz for capturing excited states in finite systems with open boundary conditions, providing a very efficient method for computing, e.g., the spectral gap of quantum spin chains. This method can be straightforwardly implemented on top of an existing DMRG or MPS ground-state code. Although this approach is built on open-boundary MPS, we also apply it to systems with periodic boundary conditions. Despite the explicit breaking of translation symmetry by the MPS representation, we show that momentum emerges as a good quantum number, and can be exploited for labeling excitations on top of MPS ground states. We apply our method to the critical Ising chain on a ring and the classical Potts model on a cylinder. Finally, we apply the same idea to compute excitation spectra for 2-D quantum systems on infinite cylinders. Again, despite the explicit breaking of translation symmetry in the periodic direction, we recover momentum as a good quantum number for labeling excitations. We apply this method to the 2-D transverse-field Ising model and the half-filled Hubbard model; for the latter, we obtain accurate results for, e.g., the hole dispersion for cylinder circumferences up to eight sites.

preprint2022arXiv

Generative modeling with projected entangled-pair states

We argue and demonstrate that projected entangled-pair states (PEPS) outperform matrix product states significantly for the task of generative modeling of datasets with an intrinsic two-dimensional structure such as images. Our approach builds on a recently introduced algorithm for sampling PEPS, which allows for the efficient optimization and sampling of the distributions.

preprint2022arXiv

Matrix quantum groups as matrix product operator representations of Lie groups

We demonstrate that the matrix quantum group $SL_q(2)$ gives rise to nontrivial matrix product operator representations of the Lie group $SL(2)$, providing an explicit characterization of the nontrivial global $SU(2)$ symmetry of the XXZ model with periodic boundary conditions. The matrix product operators are non-injective and their set is closed under multiplication. This allows to calculate the fusion tensors acting on the virtual or quantum degrees of freedom and to obtain the recoupling coefficients, which satisfy a type of pentagon relation. We argue that the combination of this data with the well known $q$-deformed Clebsch-Gordan coefficients and 6j-symbols is consistent with a description of this quantum group in terms of bimodule categories.

preprint2022arXiv

Spinon confinement and deconfinement in spin-1 chains

Motivated by the rich phase diagrams recently unveiled in frustrated spin-1 chains with competing next-nearest neighbour and four-spin interactions, we investigate the nature of the elementary excitations of spin-1 chains in the vicinity of phase transitions using the matrix-product state (MPS) representation of elementary excitations. First, we show that both spinons and magnons naturally arise in SU(2) invariant spin chains when describing ground states and elementary excitations using MPS. Then we investigate the nature of the elementary excitations across the first-oder transition between the Haldane phase and the topologically trivial next-nearest neighbor Haldane phase. We show explicitly that spinons deconfine at the transition, and we calculate the dispersion of these deconfined spinons with MPS. We also show that, immediately away from the transition line, spinons confine, forming dispersive spinon/anti-spinon bound states on both sides of the transition. Finally, we show that deconfined spinons also appear at the transition between the Haldane phase and the spontaneously dimerized phase when this transition is first order.

preprint2022arXiv

The Haldane gap in the SU(3) [3 0 0] Heisenberg chain

We calculate the Haldane gap of the $\mathrm{SU}(3)$ spin $[3~0~0]$ Heisenberg model using variational uniform fully symmetric $\mathrm{SU}(3)$ matrix product states, and find that the minimal gap $Δ/J = 0.0263 $ is obtained in the $[2~1~0]$ sector at momentum $2π/3$. We also discuss the symmetry protected topological order of the ground state, and determine the full dispersion relation of the elementary excitations and the correlation lengths of the system.

preprint2022arXiv

Variational methods for contracting projected entangled-pair states

The norms or expectation values of infinite projected entangled-pair states (PEPS) cannot be computed exactly, and approximation algorithms have to be applied. In the last years, many efficient algorithms have been devised -- the corner transfer matrix renormalization group (CTMRG) and variational uniform matrix product state (VUMPS) algorithm are the most common -- but it remains unclear whether they always lead to the same results. In this paper, we identify a subclass of PEPS for which we can reformulate the contraction as a variational problem that is algorithm independent. We use this variational feature to assess and compare the accuracy of CTMRG and VUMPS contractions. Moreover, we devise a new variational contraction scheme, which we can extend to compute general N-point correlation functions.

preprint2021arXiv

Matrix product operator symmetries and intertwiners in string-nets with domain walls

We provide a description of virtual non-local matrix product operator (MPO) symmetries in projected entangled pair state (PEPS) representations of string-net models. Given such a PEPS representation, we show that the consistency conditions of its MPO symmetries amount to a set of six coupled equations that can be identified with the pentagon equations of a bimodule category. This allows us to classify all equivalent PEPS representations and build MPO intertwiners between them, synthesising and generalising the wide variety of tensor network representations of topological phases. Furthermore, we use this generalisation to build explicit PEPS realisations of domain walls between different topological phases as constructed by Kitaev and Kong [Commun. Math. Phys. 313 (2012) 351-373]. While the prevailing abstract categorical approach is sufficient to describe the structure of topological phases, explicit tensor network representations are required to simulate these systems on a computer, such as needed for calculating thresholds of quantum error-correcting codes based on string-nets with boundaries. Finally, we show that all these string-net PEPS representations can be understood as specific instances of Turaev-Viro state-sum models of topological field theory on three-manifolds with a physical boundary, thereby putting these tensor network constructions on a mathematically rigorous footing.

preprint2021arXiv

Tangent-space methods for truncating uniform MPS

A central primitive in quantum tensor network simulations is the problem of approximating a matrix product state with one of a lower bond dimension. This problem forms the central bottleneck in algorithms for time evolution and for contracting projected entangled pair states. We formulate a tangent-space based variational algorithm to achieve this for uniform (infinite) matrix product states. The algorithm exhibits a favourable scaling of the computational cost, and we demonstrate its usefulness by several examples involving the multiplication of a matrix product state with a matrix product operator.

preprint2021arXiv

Variational optimization of continuous matrix product states

Just as matrix product states represent ground states of one-dimensional quantum spin systems faithfully, continuous matrix product states (cMPS) provide faithful representations of the vacuum of interacting field theories in one spatial dimension. Unlike the quantum spin case however, for which the density matrix renormalization group and related matrix product state algorithms provide robust algorithms for optimizing the variational states, the optimization of cMPS for systems with inhomogeneous external potentials has been problematic. We resolve this problem by constructing a piecewise linear parameterization of the underlying matrix-valued functions, which enables the calculation of the exact reduced density matrices everywhere in the system by high-order Taylor expansions. This turns the variational cMPS problem into a variational algorithm from which both the energy and its backwards derivative can be calculated exactly and at a cost that scales as the cube of the bond dimension. We illustrate this by finding ground states of interacting bosons in external potentials, and by calculating boundary or Casimir energy corrections of continuous many-body systems with open boundary conditions.

preprint2020arXiv

Real-time scattering of interacting quasiparticles in quantum spin chains

We develop a method based on tensor networks to create localized single particle excitations on top of strongly-correlated quantum spin chains. In analogy to the problem of creating localized Wannier modes, this is achieved by optimizing the gauge freedom of momentum excitations on top of matrix product states. The corresponding wavepackets propagate almost dispersionless. The time-dependent variational principle is used to scatter two such wavepackets, and we extract the phase shift from the collision data. We also study reflection and transmission coefficients of a wavepacket impinging on an impurity.

preprint2020arXiv

Restricted Boltzmann Machines for Quantum States with Nonabelian or Anyonic Symmetries

Although artificial neural networks have recently been proven to provide a promising new framework for constructing quantum many-body wave functions, the parameterization of a quantum wavefunction with nonabelian symmetries in terms of a Boltzmann machine inherently leads to biased results due to the basis dependence. We demonstrate that this problem can be overcome by sampling in the basis of irreducible representations instead of spins, for which the corresponding ansatz respects the nonabelian symmetries of the system. We apply our methodology to find the ground states of the one-dimensional antiferromagnetic Heisenberg (AFH) model with spin-half and spin-1 degrees of freedom, and obtain a substantially higher accuracy than when using the $s_z$-basis as input to the neural network. The proposed ansatz can target excited states, which is illustrated by calculating the energy gap of the AFH model. We also generalize the framework to the case of anyonic spin chains.

preprint2020arXiv

Solving frustrated Ising models using tensor networks

Motivated by the recent success of tensor networks to calculate the residual entropy of spin ice and kagome Ising models, we develop a general framework to study frustrated Ising models in terms of infinite tensor networks %, i.e. tensor networks that can be contracted using standard algorithms for infinite systems. This is achieved by reformulating the problem as local rules for configurations on overlapping clusters chosen in such a way that they relieve the frustration, i.e. that the energy can be minimized independently on each cluster. We show that optimizing the choice of clusters, including the weight on shared bonds, is crucial for the contractibility of the tensor networks, and we derive some basic rules and a linear program to implement them. We illustrate the power of the method by computing the residual entropy of a frustrated Ising spin system on the kagome lattice with next-next-nearest neighbour interactions, vastly outperforming Monte Carlo methods in speed and accuracy. The extension to finite-temperature is briefly discussed.

preprint2020arXiv

Uncertainty and Trade-offs in Quantum Multiparameter Estimation

Uncertainty relations in quantum mechanics express bounds on our ability to simultaneously obtain knowledge about expectation values of non-commuting observables of a quantum system. They quantify trade-offs in accuracy between complementary pieces of information about the system. In Quantum multiparameter estimation, such trade-offs occur for the precision achievable for different parameters characterizing a density matrix: an uncertainty relation emerges between the achievable variances of the different estimators. This is in contrast to classical multiparameter estimation, where simultaneous optimal precision is attainable in the asymptotic limit. We study trade-off relations that follow from known tight bounds in quantum multiparameter estimation. We compute trade-off curves and surfaces from Cramér--Rao type bounds which provide a compelling graphical representation of the information encoded in such bounds, and argue that bounds on simultaneously achievable precision in quantum multiparameter estimation should be regarded as measurement uncertainty relations. From the state-dependent bounds on the expected cost in parameter estimation, we derive a state independent uncertainty relation between the parameters of a qubit system.

preprint2019arXiv

Approaching the Kosterlitz-Thouless transition for the classical XY model with tensor networks

We apply variational tensor-network methods for simulating the Kosterlitz-Thouless phase transition in the classical two-dimensional XY model. In particular, using uniform matrix product states (MPS) with non-abelian O(2) symmetry, we compute the universal drop in the spin stiffness at the critical point. In the critical low-temperature regime, we focus on the MPS entanglement spectrum to characterize the Luttinger-liquid phase. In the high-temperature phase, we confirm the exponential divergence of the correlation length and estimate the critical temperature with high precision. Our MPS approach can be used to study generic two-dimensional phase transitions with continuous symmetries.

preprint2019arXiv

Galois conjugated tensor fusion categories and non-unitary CFT

We provide a generalisation of the matrix product operator (MPO) formalism for string-net projected entangled pair states (PEPS) to include non-unitary solutions of the pentagon equation. These states provide the explicit lattice realisation of the Galois conjugated counterparts of (2+1) dimensional TQFTs, based on tensor fusion categories. Although the parent Hamiltonians of these renormalisation group fixed point states are non-Hermitian, many of the topological properties of the states still hold, as a result of the pentagon equation. We show by example that the the topological sectors of the Yang-Lee theory (the non-unitary counterpart of the Fibonacci fusion category) can be constructed even in the absence of closure under Hermitian conjugation of the basis elements of the Ocneanu tube algebra. We argue that this can be generalised to the non-unitary solutions of all $SU(2)$ level $k$ models. The topological sector construction is demonstrated by applying the concept of strange correlators to the Yang-Lee model, giving rise to a non-unitary version of the classical hard hexagon model in the Yang-Lee universality class and obtaining all generalised twisted boundary conditions on a finite cylinder of the Yang-Lee edge singularity.

preprint2019arXiv

Symmetric cluster expansions with tensor networks

Cluster expansions for the exponential of local operators are constructed using tensor networks. In contrast to other approaches, the cluster expansion does not break any spatial or internal symmetries and exhibits a very favourable prefactor to the error scaling versus bond dimension. This is illustrated by time evolving a matrix product state using very large time steps, and by constructing a novel robust algorithm for finding ground states of 2-dimensional Hamiltonians using projected entangled pair states as fixed points of 2-dimensional transfer matrices.

preprint2016arXiv

A tensor product state approach to spin-1/2 square $J_1$-$J_2$ antiferromagnetic Heisenberg model: evidence for deconfined quantum criticality

The ground state phase of spin-1/2 $J_1$-$J_2$ antiferromagnetic Heisenberg model on square lattice around the maximally frustrated regime ($J_2\sim 0.5J_1$) has been debated for decades. Here we study this model using the cluster update algorithm for tensor product states (TPSs). The ground state energies at finite sizes and in the thermodynamic limit (with finite size scaling) are in good agreement with exact diagonalization study. Through finite size scaling of the spin correlation function, we find the critical point $J_2^{c_1}=0.572(5)J_1$ and critical exponents $ν=0.50(8)$, $η_s=0.28(6)$. In the range of $0.572 < J_2/J_1 \leqslant 0.6 $ we find a paramagnetic ground state with exponentially decaying spin-spin correlation. Up to $24\times 24$ system size, we observe power law decaying dimer-dimer and plaquette-plaquette correlations with an anomalous plaquette scaling exponent $η_p=0.24(1)$ and an anomalous columnar scaling exponent $η_c=0.28(1)$ at $J_2/J_1=0.6$. These results are consistent with a potential gapless $U(1)$ spin liquid phase. However, since the $U(1)$ spin liquid is unstable due to the instanton effect, a VBS order with very small amplitude might develop in the thermodynamic limit. Thus, our numerical results strongly indicate a deconfined quantum critical point (DQCP) at $J_2^{c_1}$. Remarkably, all the observed critical exponents are consistent with the $J-Q$ model.

preprint2016arXiv

Boundary-field-driven control of discontinuous phase transitions on hyperbolic lattices

The multistate Potts models on two-dimensional hyperbolic lattices are studied with respect to various boundary effects. The free energy is numerically calculated by Corner Transfer Matrix Renormalization Group method. We analyze phase transitions of the Potts models in the thermodynamic limit with respect to contracted boundary layers. A false phase transition is present even if a couple of the boundary layers are contracted. Its significance weakens, as the number of the contracted boundary layers increases, until the correct phase transition (deep inside the bulk) prevails over the false one. For this purpose we derive a thermodynamic quantity, the so-called bulk excess free energy, which depends on the contracted boundary layers and memorizes additional boundary effects. In particular, the magnetic field is imposed on the outermost boundary layer. While the boundary magnetic field does not affect the second-order phase transition in the bulk if suppressing all the boundary effects on the hyperbolic lattices, the first-order (discontinuous) phase transition is significantly sensitive to the boundary magnetic field. Contrary to the phase transition on the Euclidean lattices, the discontinuous phase transition on the hyperbolic lattices can be continuously controlled (within a certain temperature coexistence region) by varying the boundary magnetic field.

preprint2016arXiv

Confinement and string breaking for QED$_2$ in the Hamiltonian picture

The formalism of matrix product states is used to perform a numerical study of 1+1 dimensional QED -- also known as the (massive) Schwinger model -- in the presence of an external static `quark' and `antiquark'. We obtain a detailed picture of the transition from the confining state at short interquark distances to the broken-string `hadronized' state at large distances and this for a wide range of couplings, recovering the predicted behavior both in the weak and strong coupling limit of the continuum theory. In addition to the relevant local observables like charge and electric field, we compute the (bipartite) entanglement entropy and show that subtraction of its vacuum value results in a UV-finite quantity. We find that both string formation and string breaking leave a clear imprint on the resulting entropy profile. Finally, we also study the case of fractional probe charges, simulating for the first time the phenomenon of partial string breaking.

preprint2016arXiv

Gradient methods for variational optimization of projected entangled-pair states

We present a conjugate-gradient method for the ground-state optimization of projected entangled-pair states (PEPS) in the thermodynamic limit, as a direct implementation of the variational principle within the PEPS manifold. Our optimization is based on an efficient and accurate evaluation of the gradient of the global energy functional by using effective corner environments, and is robust with respect to the initial starting points. It has the additional advantage that physical and virtual symmetries can be straightforwardly implemented. We provide the tools to compute static structure factors directly in momentum space, as well as the variance of the Hamiltonian. We benchmark our method on Ising and Heisenberg models, and show a significant improvement on the energies and order parameters as compared to algorithms based on imaginary-time evolution.

preprint2016arXiv

Hamiltonian simulation of the Schwinger model at finite temperature

Using Matrix Product Operators (MPO) the Schwinger model is simulated in thermal equilibrium. The variational manifold of gauge invariant MPO is constructed to represent Gibbs states. As a first application the chiral condensate in thermal equilibrium is computed and agreement with earlier studies is found. Furthermore, as a new application the Schwinger model is probed with a fractional charged static quark-antiquark pair separated infinitely far from each other. A critical temperature beyond which the string tension is exponentially suppressed is found, which is in qualitative agreement with analytical studies in the strong coupling limit. Finally, the CT symmetry breaking is investigated and our results strongly suggest that the symmetry is restored at any nonzero temperature.

preprint2016arXiv

Matrix product operators for symmetry-protected topological phases: Gauging and edge theories

Projected entangled pair states (PEPS) provide a natural ansatz for the ground states of gapped, local Hamiltonians in which global characteristics of a quantum state are encoded in properties of local tensors. We develop a framework to describe on-site symmetries, as occurring in systems exhibiting symmetry-protected topological (SPT) quantum order, in terms of virtual symmetries of the local tensors expressed as a set of matrix product operators (MPOs) labeled by distinct group elements. These MPOs describe the possibly anomalous symmetry of the edge theory, whose local degrees of freedom are concretely identified in a PEPS. A classification of SPT phases is obtained by studying the obstructions to continuously deforming one set of MPOs into another, recovering the results derived for fixed-point models [X. Chen et al., Phys. Rev. B 87, 155114 (2013)]. Our formalism accommodates perturbations away from fixed point models, opening the possibility of studying phase transitions between different SPT phases. We also demonstrate that applying the recently developed quantum state gauging procedure to a SPT PEPS yields a PEPS with topological order determined by the initial symmetry MPOs. The MPO framework thus unifies the different approaches to classifying SPT phases, via fixed-points models, boundary anomalies, or gauging the symmetry, into the single problem of classifying inequivalent sets of matrix product operator symmetries that are defined purely in terms of a PEPS.

preprint2016arXiv

Matrix product state renormalization

The truncation or compression of the spectrum of Schmidt values is inherent to the matrix product state (MPS) approximation of one-dimensional quantum ground states. We provide a renormalization group picture by interpreting this compression as an application of Wilson's numerical renormalization group along the imaginary time direction appearing in the path integral representation of the state. The location of the physical index is considered as an impurity in the transfer matrix and static MPS correlation functions are reinterpreted as dynamical impurity correlations. Coarse-graining the transfer matrix is performed using a hybrid variational ansatz based on matrix product operators, combining ideas of MPS and the multi-scale entanglement renormalization ansatz. Through numerical comparison with conventional MPS algorithms, we explicitly verify the impurity interpretation of MPS compression, as put forward by [V. Zauner et al., New J. Phys. 17, 053002 (2015)] for the transverse-field Ising model. Additionally, we motivate the conceptual usefulness of endowing MPS with an internal layered structure by studying restricted variational subspaces to describe elementary excitations on top of the ground state, which serves to elucidate a transparent renormalization group structure ingrained in MPS descriptions of ground states.

preprint2016arXiv

Quasi-particle interactions in frustrated Heisenberg chains

Interactions between elementary excitations in quasi-one dimensional antiferromagnets are of experimental relevance and their quantitative theoretical treatment has been a theoretical challenge for many years. Using matrix product states, one can explicitly determine the wavefunctions of the one- and two-particle excitations, and, consequently, the contributions to dynamical correlations. We apply this framework to the (non integrable) frustrated dimerized spin-1/2 chain, a model for generic spin-Peierls systems, where low-energy quasi-particle excitations are bound states of topological solitons. The spin structure factor involving two quasi-particle scattering states is obtained in the thermodynamic limit with full momentum and frequency resolution. This allows very subtle features in the two-particle spectral function to be revealed which, we argue, could be seen, e.g., in inelastic neutron scattering of spin-Peierls compounds under a change of the external pressure.

preprint2015arXiv

Classical Non-equilibrium Phases from a Quantum Perspective: Quasi-adiabatic evolution and Local Ansatz for Solving The Equilibrium State of Local Stochastic Dynamics

Starting from a simple mapping of a generator of local stochastic dynamics to a quantum Hamiltonian, we derive a condition, which allows us to use the quasi-adiabatic evolution and so relate gapped quantum phases with non-equilibrium's. This leads us to a study of invertible matrix product operators. Finally, we present an ansatz for constructing local stochastic dynamics for which the Perron-Frobenius vector has a Matrix Product Representation. Additionally, we get for free that the dynamics satisfy a generalized form of detailed balance.

preprint2015arXiv

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

A natural way to generalise 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.

preprint2015arXiv

Efficient DMFT impurity solver using real-time dynamics with Matrix Product States

We propose to calculate spectral functions of quantum impurity models using the Time Evolving Block Decimation (TEBD) for Matrix Product States. The resolution of the spectral function is improved by a so-called linear prediction approach. We apply the method as an impurity solver within the Dynamical Mean Field Theory (DMFT) for the single- and two-band Hubbard model on the Bethe lattice. For the single-band model we observe sharp features at the inner edges of the Hubbard bands. A finite size scaling shows that they remain present in the thermodynamic limit. We analyze the real time-dependence of the double occupation after adding a single electron and observe oscillations at the same energy as the sharp feature in the Hubbard band, indicating a long-lived coherent superposition of states that correspond to the Kondo peak and the side peaks. For a two-band Hubbard model we observe an even richer structure in the Hubbard bands, which cannot be related to a multiplet structure of the impurity, in addition to sharp excitations at the band edges of a type similar to the single-band case.

preprint2015arXiv

Excitations and the tangent space of projected entangled-pair states

We develop tangent space methods for projected entangled-pair states (PEPS) that provide direct access to the low-energy sector of strongly-correlated two-dimensional quantum systems. More specifically, we construct a variational ansatz for elementary excitations on top of PEPS ground states that allows for computing gaps, dispersion relations, and spectral weights directly in the thermodynamic limit. Solving the corresponding variational problem requires the evaluation of momentum transformed two-point and three-point correlation functions on a PEPS background, which we can compute efficiently by using a contraction scheme. As an application we study the spectral properties of the magnons of the Affleck-Kennedy-Lieb-Tasaki model on the square lattice and the anyonic excitations in a perturbed version of Kitaev's toric code.

preprint2015arXiv

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.

preprint2015arXiv

Revisiting Many-body Localization with Random Networks of Tensors

We argue first that translational invariant Matrix Product can be interpreted as a stationary sea of particles. Next, rather than starting from some local Hamiltonian with random potentials, we consider fluctuations of the local tensors of a continuous one-parameter family of Matrix Product States. This leads to a mapping from 1+1+0 to 0+1+1, being time-space-$λ$. We argue that both pictures are equivalent. Finally, the holography principle and an operational argumentation is used to map the problem onto 1+0+1. Localization in 1-dimension, can be understood from a simple study spectral and mixing properties of finite dimensional quantum channels.

preprint2015arXiv

Scattering particles in quantum spin chains

A variational approach for constructing an effective particle description of the low-energy physics of one-dimensional quantum spin chains is presented. Based on the matrix product state formalism, we compute the one- and two-particle excitations as eigenstates of the full microscopic Hamiltonian. We interpret the excitations as particles on a strongly-correlated background with non-trivial dispersion relations, spectral weights and two-particle S matrices. Based on this information, we show how to describe a finite density of excitations as an interacting gas of bosons, using its approximate integrability at low densities. We apply our framework to the Heisenberg antiferromagnetic ladder: we compute the elementary excitation spectrum and the magnon-magnon S matrix, study the formation of bound states and determine both static and dynamic properties of the magnetized ladder.

preprint2015arXiv

Shadows of Anyons

The eigenvalue structure of the quantum transfer matrix is known to encode essential information about the elementary excitations. Here we study transfer matrices of quantum states in a topological phase using the tensor network formalism. We demonstrate that topological quantum order requires a particular type of `symmetry breaking' for the fixed point subspace of the transfer matrix, and relate physical anyon excitations to domain wall excitations at the level of the transfer matrix. A topological phase transition to a trivial phase triggers a change in the fixed point subspace to either a larger or smaller symmetry and we explain how this relates to a condensation or confinement of the corresponding anyon sectors. The tensor network formalism enables us to determine the structure of the topological sectors in two-dimensional gapped phases very efficiently, therefore opening novel avenues for studying fundamental open questions related to anyon condensation.

preprint2015arXiv

Tensor networks for gauge field theories

Over the last decade tensor network states (TNS) have emerged as a powerful tool for the study of quantum many body systems. The matrix product states (MPS) are one particular class of TNS and are used for the simulation of (1+1)-dimensional systems. In this proceeding we use MPS to determine the elementary excitations of the Schwinger model in the presence of an electric background field. We obtain an estimate for the value of the background field where the one-particle excitation with the largest energy becomes unstable and decays into two other elementary particles with smaller energy.

preprint2015arXiv

The entanglement of distillation for gauge theories

We study the entanglement structure of lattice gauge theories from the local operational point of view, and, similar to Soni and Trivedi (arXiv:1510.07455), we show that the usual entanglement entropy for a spatial bipartition can be written as the sum of an undistillable gauge part and of another part corresponding to the LOCC distillable entanglement, which is obtained by depolarizing the local superselection sectors. We demonstrate that the distillable entanglement is zero for pure abelian gauge theories in the weak coupling limit, while it is in general nonzero for the nonabelian case. We also consider gauge theories with matter, and show in a perturbative approach how area laws -- including a topological correction -- emerge for the distillable entanglement.

preprint2015arXiv

Truncating an exact Matrix Product State for the XY model: transfer matrix and its renormalisation

We discuss how to analytically obtain an -- essentially infinite -- Matrix Product State (MPS) representation of the ground state of the XY model. On the one hand this allows to illustrate how the Ornstein-Zernike form of the correlation function emerges in the exact case using standard MPS language. On the other hand we study the consequences of truncating the bond dimension of the exact MPS, which is also part of many tensor network algorithms, and analyze how the truncated MPS transfer matrix is representing the dominant part of the exact quantum transfer matrix. In the gapped phase we observe that the correlation length obtained from a truncated MPS approaches the exact value following a power law in effective bond dimension. In the gapless phase we find a good match between a state obtained numerically from standard MPS techniques with finite bond dimension, and a state obtained by effective finite imaginary time evolution in our framework. This provides a direct hint for a geometric interpretation of Finite Entanglement Scaling at the critical point in this case. Finally, by analyzing the spectra of transfer matrices, we support the interpretation put forward by [V. Zauner at. al., New J. Phys. 17, 053002 (2015)] that the MPS transfer matrix emerges from the quantum transfer matrix though the application of Wilson's Numerical Renormalisation Group along the imaginary-time direction.

preprint2015arXiv

Unifying time evolution and optimization with matrix product states

We show that the time-dependent variational principle provides a unifying framework for time-evolution methods and optimisation methods in the context of matrix product states. In particular, we introduce a new integration scheme for studying time-evolution, which can cope with arbitrary Hamiltonians, including those with long-range interactions. Rather than a Suzuki-Trotter splitting of the Hamiltonian, which is the idea behind the adaptive time-dependent density matrix renormalization group method or time-evolving block decimation, our method is based on splitting the projector onto the matrix product state tangent space as it appears in the Dirac-Frenkel time-dependent variational principle. We discuss how the resulting algorithm resembles the density matrix renormalization group (DMRG) algorithm for finding ground states so closely that it can be implemented by changing just a few lines of code and it inherits the same stability and efficiency. In particular, our method is compatible with any Hamiltonian for which DMRG can be implemented efficiently and DMRG is obtained as a special case of imaginary time evolution with infinite time step.

preprint2014arXiv

Approximating Gibbs states of local Hamiltonians efficiently with PEPS

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 ($β^{-1}$), and system size ($N$). First, we introduce a compression method in which the bond dimension scales as $D=e^{O(\log^2(N/ε))}$ if $β<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/ε\right)^{O(β)}$. 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.

preprint2014arXiv

Characterizing Topological Order with Matrix Product Operators

One of the most striking features of quantum phases that exhibit topological order is the presence of long range entanglement that cannot be detected by any local order parameter. The formalism of projected entangled-pair states is a natural framework for the parameterization of the corresponding ground state wavefunctions, in which the full wavefunction is encoded in terms of local tensors. Topological order is reflected in the symmetries of these tensors, and we give a characterization of those symmetries in terms of matrix product operators acting on the virtual level. This leads to a set of algebraic rules characterizing states with topological quantum order. The corresponding matrix product operators fully encode all topological features of the theory, and provide a systematic way of constructing topological states. We generalize the conditions of $\mathsf{G}$ and twisted injectivity to the matrix product operator case, and provide a complete picture of the ground state manifold on the torus. As an example, we show how all string-net models of Levin and Wen fit within this formalism, and in doing so provide a particularly intuitive interpretation of the pentagon equation for F-symbols as the pulling of certain matrix product operators through the string-net tensor network. Our approach paves the way to finding novel topological phases beyond string-nets, and elucidates the description of topological phases in terms of entanglement Hamiltonians and edge theories.

preprint2014arXiv

Chebyshev expansion for Impurity Models using Matrix Product States

We improve a recently developed expansion technique for calculating real frequency spectral functions of any one-dimensional model with short-range interactions, by postprocessing computed Chebyshev moments with linear prediction. This can be achieved at virtually no cost and, in sharp contrast to existing methods based on the dampening of the moments, improves the spectral resolution rather than lowering it. We validate the method for the exactly solvable resonating level model and the single impurity Anderson model. It is capable of resolving sharp Kondo resonances, as well as peaks within the Hubbard bands when employed as an impurity solver for dynamical mean-field theory (DMFT). Our method works at zero temperature and allows for arbitrary discretization of the bath spectrum. It achieves similar precision as the dynamical density matrix renormalization group (DDMRG), at lower cost. We also propose an alternative expansion, of 1-exp(-tau H) instead of the usual H, which opens the possibility of using established methods for the time evolution of matrix product states to calculate spectral functions directly.

preprint2014arXiv

Conformal Data from Finite Entanglement Scaling

In this paper we apply the formalism of translation invariant (continuous) matrix product states in the thermodynamic limit to $(1+1)$ dimensional critical models. Finite bond dimension bounds the entanglement entropy and introduces an effective finite correlation length, so that the state is perturbed away from criticality. The assumption that the scaling hypothesis holds for this kind of perturbation is known in the literature as finite entanglement scaling. We provide further evidence for the validity of finite entanglement scaling and based on this formulate a scaling algorithm to estimate the central charge and critical exponents of the conformally invariant field theories describing the critical models under investigation. The algorithm is applied to three exemplary models; the cMPS version to the non-relativistic Lieb-Liniger model and the relativistic massless boson, and MPS version to the one-dimensional quantum Ising model at the critical point. Another new aspect to our approach is that we directly use the (c)MPS induced correlation length rather than the bond dimension as scaling parameter. This choice is motivated by several theoretical arguments as well as by the remarkable accuracy of our results.

preprint2014arXiv

Entanglement Rates and the Stability of the Area Law for the Entanglement Entropy

We prove a conjecture by Bravyi on an upper bound on entanglement rates of local Hamiltonians. We then use this bound to prove the stability of the area law for the entanglement entropy of quantum spin systems under adiabatic and quasi-adiabatic evolutions.

preprint2014arXiv

Linear-optical generation of eigenstates of the two-site XY model

Much of the anticipation accompanying the development of a quantum computer relates to its application to simulating dynamics of another quantum system of interest. Here we study the building blocks for simulating quantum spin systems with linear optics. We experimentally generate the eigenstates of the XY Hamiltonian under an external magnetic field. The implemented quantum circuit consists of two CNOT gates, which are realized experimentally by harnessing entanglement from a photon source and by applying a CPhase gate. We tune the ratio of coupling constants and magnetic field by changing local parameters. This implementation of the XY model using linear quantum optics might open the door to the future studies of quenching dynamics using linear optics.

preprint2014arXiv

Matrix product states for gauge field theories

The matrix product state formalism is used to simulate Hamiltonian lattice gauge theories. To this end, we define matrix product state manifolds which are manifestly gauge invariant. As an application, we study 1+1 dimensional one flavour quantum electrodynamics, also known as the massive Schwinger model, and are able to determine very accurately the ground state properties and elementary one-particle excitations in the continuum limit. In particular, a novel particle excitation in the form of a heavy vector boson is uncovered, compatible with the strong coupling expansion in the continuum. We also study non-equilibrium dynamics by simulating the real-time evolution of the system induced by a quench in the form of a uniform background electric field.

preprint2014arXiv

Matrix product states for Hamiltonian lattice gauge theories

Over the last decade tensor network states (TNS) have emerged as a powerful tool for the study of quantum many body systems. The matrix product states (MPS) are one particular case of TNS and are used for the simulation of 1+1 dimensional systems. In [1] we considered the MPS formalism for the simulation of the Hamiltonian lattice gauge formulation of 1+1 dimensional one flavor quantum electrodynamics, also known as the massive Schwinger model. We deduced the ground state and lowest lying excitations. Furthermore, we performed a full quantum real-time simulation for a quench with a uniform background electric field. In this proceeding we continue our work on the Schwinger model. We demonstrate the advantage of working with gauge invariant MPS by comparing with MPS simulations on the full Hilbert space, that includes numerous non-physical gauge variant states. Furthermore, we compute the chiral condensate and recover the predicted UV-divergent behavior.

preprint2014arXiv

S matrix from matrix product states

We use the matrix product state formalism to construct stationary scattering states of elementary excitations in generic one-dimensional quantum lattice systems. Our method is applied to the spin-1 Heisenberg antiferromagnet, for which we calculate the full magnon-magnon S matrix for arbitrary momenta and spin, the two-particle contribution to the spectral function and the magnetization curve. As our method provides an accurate microscopic representation of the interaction between elementary excitations, we envisage the description of low-energy dynamics of one-dimensional spin chains in terms of these particlelike excitations.

preprint2014arXiv

Tensor product methods and entanglement optimization for ab initio quantum chemistry

The treatment of high-dimensional problems such as the Schrödinger equation can be approached by concepts of tensor product approximation. We present general techniques that can be used for the treatment of high-dimensional optimization tasks and time-dependent equations, and connect them to concepts already used in many-body quantum physics. Based on achievements from the past decade, entanglement-based methods, -- developed from different perspectives for different purposes in distinct communities already matured to provide a variety of tools -- can be combined to attack highly challenging problems in quantum chemistry. The aim of the present paper is to give a pedagogical introduction to the theoretical background of this novel field and demonstrate the underlying benefits through numerical applications on a text book example. Among the various optimization tasks we will discuss only those which are connected to a controlled manipulation of the entanglement which is in fact the key ingredient of the methods considered in the paper. The selected topics will be covered according to a series of lectures given on the topic "New wavefunction methods and entanglement optimizations in quantum chemistry" at the Workshop on Theoretical Chemistry, 18 - 21 February 2014, Mariapfarr, Austria.

preprint2014arXiv

The nested Algebraic Bethe Ansatz for the supersymmetric t-J and Tensor Networks

We consider a model of strongly correlated electrons in 1D called the t-J model, which was solved by graded algebraic Bethe ansatz. We use it to design graded tensor networks which can be contracted approximately to obtain a Matrix Product State. As a proof of principle, we calculate observables of ground states and excited states of finite lattices up to $18$ lattice sites.

preprint2013arXiv

An Information-Theoretic Proof of the Constructive Commutative Quantum Lovász Local Lemma

The Quantum Lovász Local Lemma (QLLL) [AKS12] establishes non-constructively that any quantum system constrained by a local Hamiltonian has a zero-energy ground state, if the local Hamiltonian terms overlap only in a certain restricted way. In this paper, we present an efficient quantum algorithm to prepare this ground state for the special case of commuting projector terms. The related classical problem has been open for more than 34 years. Our algorithm follows the breakthrough ideas of Moser's [Moser09] classical algorithm and lifts his information theoretic argument to the quantum setting. A similar result has been independently published by Arad and Sattath [AS13] recently.

preprint2013arXiv

Constructing gapless spin liquid state for the spin-1/2 J1-J2 Heisenberg model on a square lattice

We construct a class of projected entangled pair states (PEPS) which is exactly the resonating valence bond (RVB) wavefunctions endowed with both short range and long range valence bonds. With an energetically preferred RVB pattern, the wavefunction is simplified to live in a one parameter variational space. We tune this variational parameter to minimize the energy for the frustrated spin 1/2 J1-J2 antiferromagnetic Heisenberg model on the square lattice. Taking a cylindrical geometry, we are able to construct four topological sectors with even or odd number of fluxes penetrating the cylinder and even or odd number of spinons on the boundary. The energy splitting in different topological sectors is exponentially small with the cylinder perimeter. We find a power law decay of the dimer correlation function on a torus, and a lnL correction to the entanglement entropy, indicating a gapless spin liquid phase at the optimum parameter.

preprint2013arXiv

Elementary excitations in gapped quantum spin systems

For quantum lattice systems with local interactions, the Lieb-Robinson bound acts as an alternative for the strict causality of relativistic systems and allows to prove many interesting results, in particular when the energy spectrum exhibits an energy gap. In this Letter, we show that for translation invariant systems, simultaneous eigenstates of energy and momentum with an eigenvalue that is separated from the rest of the spectrum in that momentum sector, can be arbitrarily well approximated by building a momentum superposition of a local operator acting on the ground state. The error decreases in the size of the support of the local operator, with a rate that is set by the gap below and above the targeted eigenvalue. We show this explicitly for the AKLT model and discuss generalizations and applications of our result.

preprint2013arXiv

Entanglement rates and area laws

We prove an upper bound on the maximal rate at which a Hamiltonian interaction can generate entanglement in a bipartite system. The scaling of this bound as a function of the subsystem dimension on which the Hamiltonian acts nontrivially is optimal and is exponentially improved over previously known bounds. As an application, we show that a gapped quantum many-body spin system on an arbitrary lattice satisfies an area law for the entanglement entropy if and only if any other state with which it is adiabatically connected (i.e. any state in the same phase) also satisfies an area law.

preprint2013arXiv

Particles, holes and solitons: a matrix product state approach

We introduce a variational method for calculating dispersion relations of translation invariant (1+1)-dimensional quantum field theories. The method is based on continuous matrix product states and can be implemented efficiently. We study the critical Lieb-Liniger model as a benchmark and excelent agreement with the exact solution is found. Additionally, we observe solitonic signatures of Lieb's Type II excitation. In addition, a non-integrable model is introduced where a U(1)-symmetry breaking term is added to the Lieb-Liniger Hamiltonian. For this model we find evidence of a non-trivial bound-state excitation in the dispersion relation.

preprint2013arXiv

Post-Matrix Product State Methods: To tangent space and beyond

We develop in full detail the formalism of tangent states to the manifold of matrix product states, and show how they naturally appear in studying time-evolution, excitations and spectral functions. We focus on the case of systems with translation invariance in the thermodynamic limit, where momentum is a well defined quantum number. We present some new illustrative results and discuss analogous constructions for other variational classes. We also discuss generalizations and extensions beyond the tangent space, and give a general outlook towards post matrix product methods.

preprint2013arXiv

Tree tensor networks and entanglement spectra

A tree tensor network variational method is proposed to simulate quantum many-body systems with global symmetries where the optimization is reduced to individual charge configurations. A computational scheme is presented, how to extract the entanglement spectra in a bipartite splitting of a loopless tensor network across multiple links of the network, by constructing a matrix product operator for the reduced density operator and simulating its eigenstates. The entanglement spectra of 2 x L, 3 x L and 4 x L with either open or periodic boundary conditions on the rungs are studied using the presented methods, where it is found that the entanglement spectrum depends not only on the subsystem but also on the boundaries between the subsystems.

preprint2013arXiv

Variational Matrix Product Ansatz for Nonuniform Dynamics in the Thermodynamic Limit

We describe how to implement the time-dependent variational principle for matrix product states in the thermodynamic limit for nonuniform lattice systems. This is achieved by confining the nonuniformity to a (dynamically growable) finite region with fixed boundary conditions. The suppression of unphysical quasiparticle reflections from the boundary of the nonuniform region is also discussed. Using this algorithm we study the dynamics of localized excitations in infinite systems, which we illustrate in the case of the spin-1 anti-ferromagnetic Heisenberg model and the $ϕ^4$ model.

preprint2012arXiv

Calculus of continuous matrix product states

We discuss various properties of the variational class of continuous matrix product states, a class of ansatz states for one-dimensional quantum fields that was recently introduced as the direct continuum limit of the highly successful class of matrix product states. We discuss both attributes of the physical states, e.g. by showing in detail how to compute expectation values, as well as properties intrinsic to the representation itself, such as the gauge freedom. We consider general translation non-invariant systems made of several particle species and derive certain regularity properties that need to be satisfied by the variational parameters. We also devote a section to the translation invariant setting in the thermodynamic limit and show how continuous matrix product states possess an intrinsic ultraviolet cutoff. Finally, we introduce a new set of states which are tangent to the original set of continuous matrix product states. For the case of matrix product states, this construction has recently proven relevant in the development of new algorithms for studying time evolution and elementary excitations of quantum spin chains. We thus lay the foundation for similar developments for one-dimensional quantum fields.

preprint2012arXiv

Geometry of Matrix Product States: metric, parallel transport and curvature

We study the geometric properties of the manifold of states described as (uniform) matrix product states. Due to the parameter redundancy in the matrix product state representation, matrix product states have the mathematical structure of a (principal) fiber bundle. The total space or bundle space corresponds to the parameter space, i.e. the space of tensors associated to every physical site. The base manifold is embedded in Hilbert space and can be given the structure of a Kähler manifold by inducing the Hilbert space metric. Our main interest is in the states living in the tangent space to the base manifold, which have recently been shown to be interesting in relation to time dependence and elementary excitations. By lifting these tangent vectors to the (tangent space) of the bundle space using a well-chosen prescription (a principal bundle connection), we can define and efficiently compute an inverse metric, and introduce differential geometric concepts such as parallel transport (related to the Levi-Civita connection) and the Riemann curvature tensor.

preprint2012arXiv

Numerical Contraction of the Tensor Network generated by the Algebraic Bethe Ansatz

The algebraic Bethe Ansatz is a prosperous and well-established method for solving one-dimensional quantum models exactly. The solution of the complex eigenvalue problem is thereby reduced to the solution of a set of algebraic equations. Whereas the spectrum is usually obtained directly, the eigenstates are available only in terms of complex mathematical expressions. This makes it very hard in general to extract properties from the states, like, for example, correlation functions. In our work, we apply the tools of Tensor Network States to describe the eigenstates approximately as Matrix Product States. From the Matrix Product State expression, we then obtain observables like correlation functions directly.

preprint2012arXiv

Preparing projected entangled pair states on a quantum computer

We present a quantum algorithm to prepare injective PEPS on a quantum computer, a class of open tensor networks representing quantum states. The run-time of our algorithm scales polynomially with the inverse of the minimum condition number of the PEPS projectors and, essentially, with the inverse of the spectral gap of the PEPS' parent Hamiltonian.

preprint2012arXiv

Preparing topological PEPS on a quantum computer

Simulating of exotic phases of matter that are not amenable to classical techniques is one of the most important potential applications of quantum information processing. We present an efficient algorithm for preparing a large class of topological quantum states -- the G-injective Projected Entangled Pair States (PEPS) -- on a quantum computer. Important examples include the resonant valence bond (RVB) states, conjectured to be topological spin liquids. The runtime of the algorithm scales polynomially with the condition number of the PEPS projectors, and inverse-polynomially in the spectral gap of the PEPS parent Hamiltonian.

preprint2012arXiv

Quantum state discrimination bounds for finite sample size

In the problem of quantum state discrimination, one has to determine by measurements the state of a quantum system, based on the a priori side information that the true state is one of two given and completely known states, rho or sigma. In general, it is not possible to decide the identity of the true state with certainty, and the optimal measurement strategy depends on whether the two possible errors (mistaking rho for sigma, or the other way around) are treated as of equal importance or not. Results on the quantum Chernoff and Hoeffding bounds and the quantum Stein's lemma show that, if several copies of the system are available then the optimal error probabilities decay exponentially in the number of copies, and the decay rate is given by a certain statistical distance between rho and sigma (the Chernoff distance, the Hoeffding distances, and the relative entropy, respectively). While these results provide a complete solution to the asymptotic problem, they are not completely satisfying from a practical point of view. Indeed, in realistic scenarios one has access only to finitely many copies of a system, and therefore it is desirable to have bounds on the error probabilities for finite sample size. In this paper we provide finite-size bounds on the so-called Stein errors, the Chernoff errors, the Hoeffding errors and the mixed error probabilities related to the Chernoff and the Hoeffding errors.

preprint2012arXiv

The Algebraic Bethe Ansatz and Tensor Networks

We describe the Algebraic Bethe Ansatz for the spin-1/2 XXX and XXZ Heisenberg chains with open and periodic boundary conditions in terms of tensor networks. These Bethe eigenstates have the structure of Matrix Product States with a conserved number of down-spins. The tensor network formulation suggestes possible extensions of the Algebraic Bethe Ansatz to two dimensions.

preprint2012arXiv

The continuum limit of a tensor network: a path integral representation

We argue that the natural way to generalise a tensor network variational class to a continuous quantum system is to use the Feynman path integral to implement a continuous tensor contraction. This approach is illustrated for the case of a recently introduced class of quantum field states known as continuous matrix-product states (cMPS). As an example of the utility of the path-integral representation we argue that the state of a dynamically evolving quantum field admits a natural representation as a cMPS. An argument that all states in Fock space admit a cMPS representation when the number of variational parameters tends to infinity is also provided.

preprint2012arXiv

Variational matrix product ansatz for dispersion relations

A variational ansatz for momentum eigenstates of translation invariant quantum spin chains is formulated. The matrix product state ansatz works directly in the thermodynamic limit and allows for an efficient implementation (cubic scaling in the bond dimension) of the variational principle. Unlike previous approaches, the ansatz includes topologically non-trivial states (kinks, domain walls) for systems with symmetry breaking. The method is benchmarked using the spin-1/2 XXZ antiferromagnet and the spin-1 Heisenberg antiferromagnet and we obtain surprisingly accurate results.

preprint2012arXiv

Variational Numerical Renormalization Group: Bridging the gap between NRG and Density Matrix Renormalization Group

The numerical renormalization group (NRG) is rephrased as a variational method with the cost function given by the sum of all the energies of the effective low-energy Hamiltonian. This allows to systematically improve the spectrum obtained by NRG through sweeping. The ensuing algorithm has a lot of similarities to the density matrix renormalization group (DMRG) when targeting many states, and this synergy of NRG and DMRG combines the best of both worlds and extends their applicability. We illustrate this approach with simulations of a quantum spin chain and a single impurity Anderson model (SIAM) where the accuracy of the effective eigenstates is greatly enhanced as compared to the NRG, especially in the transition to the continuum limit.

preprint2011arXiv

A matrix product state based algorithm for determining dispersion relations of quantum spin chains with periodic boundary conditions

We study a matrix product state (MPS) algorithm to approximate excited states of translationally invariant quantum spin systems with periodic boundary conditions. By means of a momentum eigenstate ansatz generalizing the one of Östlund and Rommer [1], we separate the Hilbert space of the system into subspaces with different momentum. This gives rise to a direct sum of effective Hamiltonians, each one corresponding to a different momentum, and we determine their spectrum by solving a generalized eigenvalue equation. Surprisingly, many branches of the dispersion relation are approximated to a very good precision. We benchmark the accuracy of the algorithm by comparison with the exact solutions of the quantum Ising and the antiferromagnetic Heisenberg spin-1/2 model.

preprint2011arXiv

Applying the variational principle to (1+1)-dimensional quantum field theories

We extend the recently introduced continuous matrix product state (cMPS) variational class to the setting of (1+1)-dimensional relativistic quantum field theories. This allows one to overcome the difficulties highlighted by Feynman concerning the variational procedure applied to relativistic theories, and provides a new way to regularize quantum field theories. A fermionic version of the continuous matrix product state is introduced which is manifestly free from fermion doubling and sign problems. We illustrate the power of the formalism with the simulation of free massive Dirac fermions, the Gross-Neveu model, and the Casimir effect. We find that cMPS can capture chiral symmetry breaking with absolute scaling of the chiral parameter, and that boundary effects can be accommodated with modest computational effort.

preprint2011arXiv

Cluster update for tensor network states

We propose a novel recursive way of updating the tensors in projected entangled pair states by evolving the tensor in imaginary time evolution on clusters of different sizes. This generalizes the so- called simple update method of Jiang et al. [Phys. Rev. Lett. 101, 090603 (2008)] and the updating schemes in the single layer picture of Pižorn et al. [Phys. Rev. A 83, 052321 (2011)]. A finite-size scaling of the observables as a function of the cluster size provides a remarkable improvement in the accuracy as compared to the simple update scheme. We benchmark our results on the hand of the spin 1/2 staggered dimerized antiferromagnetic model on the square lattice, and accurate results for the magnetization and the critical exponents are determined.

preprint2011arXiv

Entanglement renormalization for quantum fields

It is shown how to construct renormalization group flows of quantum field theories in real space, as opposed to the usual Wilsonian approach in momentum space. This is achieved by generalizing the multiscale entanglement renormalization ansatz to continuum theories. The variational class of wavefunctions arising from this RG flow are translation invariant and exhibit an entropy-area law. We illustrate the construction for a free non-relativistic boson model, and argue that the full power of the construction should emerge in the case of interacting theories.

preprint2011arXiv

Entanglement spectrum and boundary theories with projected entangled-pair states

In many physical scenarios, close relations between the bulk properties of quantum systems and theories associated to their boundaries have been observed. In this work, we provide an exact duality mapping between the bulk of a quantum spin system and its boundary using Projected Entangled Pair States (PEPS). This duality associates to every region a Hamiltonian on its boundary, in such a way that the entanglement spectrum of the bulk corresponds to the excitation spectrum of the boundary Hamiltonian. We study various specific models, like a deformed AKLT [1], an Ising-type [2], and Kitaev's toric code [3], both in finite ladders and infinite square lattices. In the latter case, some of those models display quantum phase transitions. We find that a gapped bulk phase with local order corresponds to a boundary Hamiltonian with local interactions, whereas critical behavior in the bulk is reflected on a diverging interaction length of the boundary Hamiltonian. Furthermore, topologically ordered states yield non-local Hamiltonians. As our duality also associates a boundary operator to any operator in the bulk, it in fact provides a full holographic framework for the study of quantum many-body systems via their boundary.

preprint2011arXiv

Monte Carlo simulation with Tensor Network States

It is demonstrated that Monte Carlo sampling can be used to efficiently extract the expectation value of projected entangled pair states with large virtual bond dimension. We use the simple update rule introduced by Xiang et al. to obtain the tensors describing the ground state wavefunction of the antiferromagnetic Heisenberg model and evaluate the finite size energy and staggered magnetization for square lattices with periodic boundary conditions of sizes up to L=16 and virtual bond dimensions up to D=16. The finite size magnetization errors are 0.003(2) and 0.013(2) at D=16 for a system of size L=8,16 respectively. Finite D extrapolation provides exact finite size magnetization for L=8, and reduces the magnetization error to 0.005(3) for L=16, significantly improving the previous state of the art results.

preprint2011arXiv

Quantum simulation of time-dependent Hamiltonians and the convenient illusion of Hilbert space

We consider the manifold of all quantum many-body states that can be generated by arbitrary time-dependent local Hamiltonians in a time that scales polynomially in the system size, and show that it occupies an exponentially small volume in Hilbert space. This implies that the overwhelming majority of states in Hilbert space are not physical as they can only be produced after an exponentially long time. We establish this fact by making use of a time-dependent generalization of the Suzuki-Trotter expansion, followed by a counting argument. This also demonstrates that a computational model based on arbitrarily rapidly changing Hamiltonians is no more powerful than the standard quantum circuit model.

preprint2011arXiv

Time evolution of projected entangled pair states in the single layer picture

We propose an efficient algorithm for simulating quantum many-body systems in two spatial dimensions using projected entangled pair states. This is done by approximating the environment, arising in the context of updating tensors in the process of time evolution, using a single-layered tensor network structure. This significantly reduces the computational costs and allows simulations in a larger submanifold of the Hilbert space as bounded by the bond dimension of the tensor network. We present numerical evidence for stability of the method on an antiferromagnetic isotropic Heisenberg model where good agreement is found with the available accurate results.

preprint2011arXiv

Time-dependent variational principle for quantum lattices

We develop a new algorithm based on the time-dependent variational principle applied to matrix product states to efficiently simulate the real- and imaginary time dynamics for infinite one-dimensional quantum lattice systems. This procedure: (1) is argued to be optimal; (2) does not rely on the Trotter decomposition and thus has no Trotter error; (3) explicitly preserves all symmetries and conservation laws; and (4) has low computational complexity. The algorithm is illustrated using both imaginary time and real-time examples.

preprint2010arXiv

Complete-Graph Tensor Network States: A New Fermionic Wave Function Ansatz for Molecules

We present a new class of tensor network states that are specifically designed to capture the electron correlation of a molecule of arbitrary structure. In this ansatz, the electronic wave function is represented by a Complete-Graph Tensor Network (CGTN) ansatz which implements an efficient reduction of the number of variational parameters by breaking down the complexity of the high-dimensional coefficient tensor of a full-configuration-interaction (FCI) wave function. We demonstrate that CGTN states approximate ground states of molecules accurately by comparison of the CGTN and FCI expansion coefficients. The CGTN parametrization is not biased towards any reference configuration in contrast to many standard quantum chemical methods. This feature allows one to obtain accurate relative energies between CGTN states which is central to molecular physics and chemistry. We discuss the implications for quantum chemistry and focus on the spin-state problem. Our CGTN approach is applied to the energy splitting of states of different spin for methylene and the strongly correlated ozone molecule at a transition state structure. The parameters of the tensor network ansatz are variationally optimized by means of a parallel-tempering Monte Carlo algorithm.

preprint2010arXiv

Computational Complexity of interacting electrons and fundamental limitations of Density Functional Theory

One of the central problems in quantum mechanics is to determine the ground state properties of a system of electrons interacting via the Coulomb potential. Since its introduction by Hohenberg, Kohn, and Sham, Density Functional Theory (DFT) has become the most widely used and successful method for simulating systems of interacting electrons, making their original work one of the most cited in physics. In this letter, we show that the field of computational complexity imposes fundamental limitations on DFT, as an efficient description of the associated universal functional would allow to solve any problem in the class QMA (the quantum version of NP) and thus particularly any problem in NP in polynomial time. This follows from the fact that finding the ground state energy of the Hubbard model in an external magnetic field is a hard problem even for a quantum computer, while given the universal functional it can be computed efficiently using DFT. This provides a clear illustration how the field of quantum computing is useful even if quantum computers would never be built.

preprint2010arXiv

Fermionic Implementation of Projected Entangled Pair States Algorithm

We present and implement an efficient variational method to simulate two-dimensional finite size fermionic quantum systems by fermionic projected entangled pair states. The approach differs from the original one due to the fact that there is no need for an extra string-bond for contracting the tensor network. The method is tested on a bi-linear fermionic model on a square lattice for sizes up to ten by ten where good relative accuracy is achieved. Qualitatively good results are also obtained for an interacting fermionic system.

preprint2010arXiv

Fermionic Projected Entangled Pair States

We introduce a family of states, the fPEPS, which describes fermionic systems on lattices in arbitrary spatial dimensions. It constitutes the natural extension of another family of states, the PEPS, which efficiently approximate ground and thermal states of spin systems with short-range interactions. We give an explicit mapping between those families, which allows us to extend previous simulation methods to fermionic systems. We also show that fPEPS naturally arise as exact ground states of certain fermionic Hamiltonians. We give an example of such a Hamiltonian, exhibiting criticality while obeying an area law.

preprint2010arXiv

Grassmann tensor network states and its renormalization for strongly correlated fermionic and bosonic states

The projective construction (the slave-particle approach) has played an very important role in understanding strongly correlated systems, such as the emergence of fermions, anyons, and gauge theory in quantum spin liquids and quantum Hall states. Recently, fermionic Projected Entangled Pair States (fPEPS) have been introduced to effciently represent many-body fermionic states. In this paper, we show that the strongly correlated bosonic/fermionic states obtained both from the projective construction and the fPEPS approach can be represented systematically as Grassmann tensor product states. This construction can also be applied to all other tensor network states approaches. The Grassmann tensor product states allow us to encode many-body bosonic/fermionic states effciently with a polynomial number of parameters. We also generalize the tensor-entanglement renormalization group (TERG) method for complex tensor networks to Grassmann tensor networks. This allows us to approximate the norm and average local operators of Grassmann tensor product states in polynomial time, and hence leads to a variational approach for describing strongly correlated bosonic/fermionic systems in higher dimensions.

preprint2010arXiv

Holographic quantum states

We show how continuous matrix product states of quantum field theories can be described in terms of the dissipative non-equilibrium dynamics of a lower-dimensional auxiliary boundary field theory. We demonstrate that the spatial correlation functions of the bulk field can be brought into one-to-one correspondence with the temporal statistics of the quantum jumps of the boundary field. This equivalence: (1) illustrates an intimate connection between the theory of continuous quantum measurement and quantum field theory; (2) gives an explicit construction of the boundary field theory allowing the extension of real-space renormalization group methods to arbitrary dimensional quantum field theories without the introduction of a lattice parameter; and (3) yields a novel interpretation of recent cavity QED experiments in terms of quantum field theory, and hence paves the way toward observing genuine quantum phase transitions in such zero-dimensional driven quantum systems.

preprint2010arXiv

Mutual information in classical spin models

The total many-body correlations present in finite temperature classical spin systems are studied using the concept of mutual information. As opposed to zero-temperature quantum phase transitions, the total correlations are not maximal at the phase transition, but reach a maximum in the high temperature paramagnetic phase. The Shannon and Renyi mutual information in both Ising and Potts models in 2 dimensions are calculated numerically by combining matrix product states algorithms and Monte Carlo sampling techniques.

preprint2010arXiv

Simulating Strongly Correlated Quantum Systems with Tree Tensor Networks

We present a tree-tensor-network-based method to study strongly correlated systems with nonlocal interactions in higher dimensions. Although the momentum-space and quantum-chemistry versions of the density matrix renormalization group (DMRG) method have long been applied to such systems, the spatial topology of DMRG-based methods allows efficient optimizations to be carried out with respect to one spatial dimension only. Extending the matrix-product-state picture, we formulate a more general approach by allowing the local sites to be coupled to more than two neighboring auxiliary subspaces. Following Shi. et. al. [Phys. Rev. A, 74, 022320 (2006)], we treat a tree-like network ansatz with arbitrary coordination number z, where the z=2 case corresponds to the one-dimensional scheme. For this ansatz, the long-range correlation deviates from the mean-field value polynomially with distance, in contrast to the matrix-product ansatz, which deviates exponentially. The computational cost of the tree-tensor-network method is significantly smaller than that of previous DMRG-based attempts, which renormalize several blocks into a single block. In addition, we investigate the effect of unitary transformations on the local basis states and present a method for optimizing such transformations. For the 1-d interacting spinless fermion model, the optimized transformation interpolates smoothly between real space and momentum space. Calculations carried out on small quantum chemical systems support our approach.

preprint2010arXiv

Simulation of interacting fermions with entanglement renormalization

We propose and test an algorithm to simulate a lattice system of interacting fermions in two spatial dimensions. The approach is an extension of the entanglement renormalization technique [Phys. Rev. Lett. 99, 220405 (2007)] and the related multi-scale entanglement renormalization ansatz. Benchmark calculations for free and interacting fermions on a lattice of $6\times 6$ sites and periodic boundary conditions confirm the validity of this proposal.

preprint2010arXiv

Stochastic exclusion processes versus coherent transport

Stochastic exclusion processes play an integral role in the physics of non-equilibrium statistical mechanics. These models are Markovian processes, described by a classical master equation. In this paper a quantum mechanical version of a stochastic hopping process in one dimension is formulated in terms of a quantum master equation. This allows the investigation of coherent and stochastic evolution in the same formal framework. The focus lies on the non-equilibrium steady state. Two stochastic model systems are considered, the totally asymmetric exclusion process and the fully symmetric exclusion process. The steady state transport properties of these models is compared to the case with additional coherent evolution, generated by the $XX$-Hamiltonian.

preprint2010arXiv

Stochastic Matrix Product States

The concept of stochastic matrix product states is introduced and a natural form for the states is derived. This allows to define the analogue of Schmidt coefficients for steady states of non-equilibrium stochastic processes. We discuss a new measure for correlations which is analogous to the entanglement entropy, the entropy cost $S_C$, and show that this measure quantifies the bond dimension needed to represent a steady state as a matrix product state. We illustrate these concepts on the hand of the asymmetric exclusion process.

preprint2007arXiv

The computational complexity of PEPS

We determine the computational power of preparing Projected Entangled Pair States (PEPS), as well as the complexity of classically simulating them, and generally the complexity of contracting tensor networks. While creating PEPS allows to solve PP problems, the latter two tasks are both proven to be #P-complete. We further show how PEPS can be used to approximate ground states of gapped Hamiltonians, and that creating them is easier than creating arbitrary PEPS. The main tool for our proofs is a duality between PEPS and postselection which allows to use existing results from quantum compexity.

Frank Verstraete

What is connected

Connect this record

See the researcher in context

Building this map preview

88 published item(s)

Contrasting pseudo-criticality in the classical two-dimensional Heisenberg and $\mathrm{RP}^2$ models: zero-temperature phase transition versus finite-temperature crossover

Direct sampling of projected entangled-pair states

Efficient MPS methods for extracting spectral information on rings and cylinders

Generative modeling with projected entangled-pair states

Matrix quantum groups as matrix product operator representations of Lie groups

Spinon confinement and deconfinement in spin-1 chains

The Haldane gap in the SU(3) [3 0 0] Heisenberg chain

Variational methods for contracting projected entangled-pair states

Matrix product operator symmetries and intertwiners in string-nets with domain walls

Tangent-space methods for truncating uniform MPS

Variational optimization of continuous matrix product states

Real-time scattering of interacting quasiparticles in quantum spin chains

Restricted Boltzmann Machines for Quantum States with Nonabelian or Anyonic Symmetries

Solving frustrated Ising models using tensor networks

Uncertainty and Trade-offs in Quantum Multiparameter Estimation

Approaching the Kosterlitz-Thouless transition for the classical XY model with tensor networks

Galois conjugated tensor fusion categories and non-unitary CFT

Symmetric cluster expansions with tensor networks

A tensor product state approach to spin-1/2 square $J_1$-$J_2$ antiferromagnetic Heisenberg model: evidence for deconfined quantum criticality

Boundary-field-driven control of discontinuous phase transitions on hyperbolic lattices

Confinement and string breaking for QED$_2$ in the Hamiltonian picture

Gradient methods for variational optimization of projected entangled-pair states

Hamiltonian simulation of the Schwinger model at finite temperature

Matrix product operators for symmetry-protected topological phases: Gauging and edge theories

Matrix product state renormalization

Quasi-particle interactions in frustrated Heisenberg chains

Classical Non-equilibrium Phases from a Quantum Perspective: Quasi-adiabatic evolution and Local Ansatz for Solving The Equilibrium State of Local Stochastic Dynamics

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

Efficient DMFT impurity solver using real-time dynamics with Matrix Product States

Excitations and the tangent space of projected entangled-pair states

Gauging quantum states: from global to local symmetries in many-body systems

Revisiting Many-body Localization with Random Networks of Tensors

Scattering particles in quantum spin chains

Shadows of Anyons

Tensor networks for gauge field theories

The entanglement of distillation for gauge theories

Truncating an exact Matrix Product State for the XY model: transfer matrix and its renormalisation

Unifying time evolution and optimization with matrix product states

Approximating Gibbs states of local Hamiltonians efficiently with PEPS

Characterizing Topological Order with Matrix Product Operators

Chebyshev expansion for Impurity Models using Matrix Product States

Conformal Data from Finite Entanglement Scaling

Entanglement Rates and the Stability of the Area Law for the Entanglement Entropy

Linear-optical generation of eigenstates of the two-site XY model

Matrix product states for gauge field theories

Matrix product states for Hamiltonian lattice gauge theories

S matrix from matrix product states

Tensor product methods and entanglement optimization for ab initio quantum chemistry

The nested Algebraic Bethe Ansatz for the supersymmetric t-J and Tensor Networks

An Information-Theoretic Proof of the Constructive Commutative Quantum Lovász Local Lemma

Constructing gapless spin liquid state for the spin-1/2 J1-J2 Heisenberg model on a square lattice

Elementary excitations in gapped quantum spin systems

Entanglement rates and area laws

Particles, holes and solitons: a matrix product state approach

Post-Matrix Product State Methods: To tangent space and beyond

Tree tensor networks and entanglement spectra

Variational Matrix Product Ansatz for Nonuniform Dynamics in the Thermodynamic Limit

Calculus of continuous matrix product states

Geometry of Matrix Product States: metric, parallel transport and curvature

Numerical Contraction of the Tensor Network generated by the Algebraic Bethe Ansatz

Preparing projected entangled pair states on a quantum computer

Preparing topological PEPS on a quantum computer

Quantum state discrimination bounds for finite sample size

The Algebraic Bethe Ansatz and Tensor Networks

The continuum limit of a tensor network: a path integral representation

Variational matrix product ansatz for dispersion relations

Variational Numerical Renormalization Group: Bridging the gap between NRG and Density Matrix Renormalization Group

A matrix product state based algorithm for determining dispersion relations of quantum spin chains with periodic boundary conditions

Applying the variational principle to (1+1)-dimensional quantum field theories

Cluster update for tensor network states

Entanglement renormalization for quantum fields

Entanglement spectrum and boundary theories with projected entangled-pair states

Monte Carlo simulation with Tensor Network States

Quantum simulation of time-dependent Hamiltonians and the convenient illusion of Hilbert space