Source author record

Laurens Vanderstraeten

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

22works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

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

22 published item(s)

preprint2026arXiv

Interface roughening in the 3-D Ising model with tensor networks

Interfaces in three-dimensional many-body systems can exhibit rich phenomena beyond the corresponding bulk properties. In particular, they can fluctuate and give rise to massless low energy degrees of freedom even in the presence of a gapped bulk. In this work, we present the first tensor-network study of the paradigmatic interface roughening transition of the 3-D Ising model using highly asymmetric lattices that are infinite in the $(xy)$ direction and finite in $z$. By reducing the problem to an effective 2-D tensor network, we study how truncating the $z$ direction reshapes the physics of the interface. For a truncation based on open boundary conditions, we demonstrate that varying the interface width gives rise to either a $\mathbb{Z}_2$ symmetry breaking transition (for odd $L_z$) or a smooth crossover(for even $L_z$). For antiperiodic boundary conditions, we obtain an effective $\mathbb{Z}_q$ clock model description with $q=2L_z$ that exhibits an intermediate Luttinger liquid phase with an emergent $\U(1)$ symmetry.

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

Momentum-resolved time evolution with matrix product states

We introduce a method based on matrix product states (MPS) for computing spectral functions of (quasi) one-dimensional spin chains, working directly in momentum space in the thermodynamic limit. We simulate the time evolution after applying a momentum operator to an MPS ground state by working with the momentum superposition of a window MPS. We show explicitly for the spin-1 Heisenberg chain that the growth of entanglement is smaller in momentum space, even inside a two-particle continuum, such that we can attain very accurate spectral functions with relatively small bond dimension. We apply our method to compute spectral lineshapes of the gapless XXZ chain and the square-lattice J1-J2 Heisenberg model on a six-leg cylinder.

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

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.

preprint2020arXiv

Quasiparticle origin of dynamical quantum phase transitions

Considering nonintegrable quantum Ising chains with exponentially decaying interactions, we present matrix product state results that establish a connection between low-energy quasiparticle excitations and the kind of nonanalyticities in the Loschmidt return rate. When domain walls in the spectrum of the quench Hamiltonian are energetically favored to be bound rather than freely propagating, anomalous cusps appear in the return rate regardless of the initial state. In the nearest-neighbor limit, domain walls are always freely propagating, and anomalous cusps never appear. As a consequence, our work illustrates that models in the same equilibrium universality class can still exhibit fundamentally distinct out-of-equilibrium criticality. Our results are accessible to current ultracold-atom and ion-trap experiments.

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

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.

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

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

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

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

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

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.

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.

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.