Researcher profile

Jutho Haegeman

Jutho Haegeman contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

13 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.

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

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

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.

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

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.

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.

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.

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.