BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jutho Haegeman

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

ResearcherUnclaimed source record
quant-phcond-mat.str-elcond-mat.stat-mechhep-lathep-thmath-phmath.MPcond-mat.quant-gascond-mat.dis-nnmath.QA

Catalog footprint

What is connected

46works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

46 published item(s)

preprint2025arXiv

Numerical study of boson mixtures with multi-component continuous matrix product states

The continuous matrix product state (cMPS) ansatz is a promising numerical tool for studying quantum many-body systems in continuous space. Although it provides a clean framework that allows one to directly simulate continuous systems, the optimization of cMPS is known to be a very challenging task, especially in the case of multi-component systems. In this work, we have developed an improved optimization scheme for multi-component cMPS that enables simulations of bosonic quantum mixtures with substantially larger bond dimensions than previous works. We benchmark our method on the two-component Lieb-Liniger model, obtaining numerical results that agree well with analytical predictions. Our work paves the way for further numerical studies of quantum mixture systems using the cMPS ansatz.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Riemannian optimization of isometric tensor networks

Several tensor networks are built of isometric tensors, i.e. tensors satisfying $W^\dagger W = \mathrm{I}$. Prominent examples include matrix product states (MPS) in canonical form, the multiscale entanglement renormalization ansatz (MERA), and quantum circuits in general, such as those needed in state preparation and quantum variational eigensolvers. We show how gradient-based optimization methods on Riemannian manifolds can be used to optimize tensor networks of isometries to represent e.g. ground states of 1D quantum Hamiltonians. We discuss the geometry of Grassmann and Stiefel manifolds, the Riemannian manifolds of isometric tensors, and review how state-of-the-art optimization methods like nonlinear conjugate gradient and quasi-Newton algorithms can be implemented in this context. We apply these methods in the context of infinite MPS and MERA, and show benchmark results in which they outperform the best previously-known optimization methods, which are tailor-made for those specific variational classes. We also provide open-source implementations of our algorithms.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The Chiral Gross-Neveu model on the lattice via a Landau-forbidden phase transition

We study the phase diagram of the $(1+1)$-dimensional Gross-Neveu model with both $g_x^2(\barψψ)^2$ and $g_y^2(\barψiγ_5ψ)^2$ interaction terms on a spatial lattice. The continuous chiral symmetry, which is present in the continuum model when $g_x^2=g_y^2$, has a mixed 't~Hooft anomaly with the charge conservation symmetry, which guarantees the existence of a massless mode. However, the same 't~Hooft anomaly also implies that the continuous chiral symmetry is broken explicitly in our lattice model. Nevertheless, from numerical matrix product state simulations we find that, for certain parameters of the lattice model, the continuous chiral symmetry reemerges in the infrared fixed point theory, even at strong coupling. We argue that, to understand this phenomenon, it is crucial to go beyond mean-field theory (or, equivalently, beyond the leading order term in a $1/N$ expansion). Interestingly, on the lattice, the chiral Gross-Neveu model appears at a Landau-forbidden second order phase transition separating two distinct and unrelated symmetry-breaking orders. We point out the crucial role of two different 't Hooft anomalies or Lieb-Schultz-Mattis obstructions for this Landau-forbidden phase transition to occur.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Matrix product states and variational methods applied to critical quantum field theory

We study the second-order quantum phase-transition of massive real scalar field theory with a quartic interaction ($ϕ^4$ theory) in (1+1) dimensions on an infinite spatial lattice using matrix product states (MPS). We introduce and apply a naive variational conjugate gradient method, based on the time-dependent variational principle (TDVP) for imaginary time, to obtain approximate ground states, using a related ansatz for excitations to calculate the particle and soliton masses and to obtain the spectral density. We also estimate the central charge using finite-entanglement scaling. Our value for the critical parameter agrees well with recent Monte Carlo results, improving on an earlier study which used the related DMRG method, verifying that these techniques are well-suited to studying critical field systems. We also obtain critical exponents that agree, as expected, with those of the transverse Ising model. Additionally, we treat the special case of uniform product states (mean field theory) separately, showing that they may be used to investigate non-critical quantum field theories under certain conditions.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An order parameter for symmetry-protected phases in one dimension

We introduce an order parameter for symmetry-protected phases in one dimension which allows to directly identify those phases. The order parameter consists of string-like operators and swaps, but differs from conventional string order operators in that it only depends on the symmetry but not on the state. We verify our framework through numerical simulations for the SO(3) invariant spin-1 bilinear-biquadratic model which exhibits a dimerized and a Haldane phase, and find that the order parameter not only works very well for the dimerized and the Haldane phase, but it also returns a distinct signature for gapless phases. Finally, we discuss possible ways to measure the order parameter in experiments with cold atoms.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Minimally doubled chiral fermions with C, P and T symmetry on the staggered lattice

Recently, the interest in local lattice actions for chiral fermions has revived, with the proposition of new local actions in which only the minimal number of doublers appear. The trigger role of graphene having a minimally doubled, chirally invariant, Dirac-like excitation spectrum can not be neglected. The challenge is to construct an action which preserves enough symmetries to be useful in lattice gauge calculations. We present a new approach to obtain local lattice actions for fermions using a reinterpretation of the staggered lattice approach of Kogut and Susskind. This interpretation is based on the similarity with the staggered lattice approach in FDTD simulations of acoustics and electromagnetism. It allows us to construct a local action for chiral fermions which has all discrete symmetries and the minimal number of fermion flavors, but which is non-Hermitian in real space. However, we argue that this will not pose a threat to the usability of the theory.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access