Researcher profile

Laurens Vanderstraeten

Laurens Vanderstraeten contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
9topics
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

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