Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2207.11301 · SciPost Physics

A path integral ground state Monte Carlo algorithm for entanglement of lattice bosons

Emanuel Casiano-Diaz, C. M. Herdman, Adrian Del Maestro

Published Mar 29, 2023• 6 citations
arXivHTMLPDF

Abstract

A ground state path integral quantum Monte Carlo algorithm is introduced that allows for the study of entanglement in lattice bosons at zero temperature. The Rényi entanglement entropy between spatial subregions is explored across the phase diagram of the one dimensional Bose-Hubbard model for systems consisting of up to L=256 L = 256 sites at unit-filling without any restrictions on site occupancy, far beyond the reach of exact diagonalization. The favorable scaling of the algorithm is demonstrated through a further measurement of the Rényi entanglement entropy at the two dimensional superfluid-insulator critical point for large system sizes, confirming the existence of the expected entanglement boundary law in the ground state. The Rényi estimator is extended to measure the symmetry resolved entanglement that is operationally accessible as a resource for experimentally relevant lattice gases with fixed total particle number.

DOI
10.21468/scipostphys.14.3.054

Related papers

Path-integral Monte Carlo method for Rényi entanglement entropies

C. M. Herdman, Stephen Inglis, Pierre–Nicholas Roy +2 more

We introduce a quantum Monte Carlo algorithm to measure the Rényi entanglement entropies in systems of interacting bosons in the continuum. This approach is based on a path-integral ground state method that can be applied to interacting itinerant bosons in any spatial dimension with direct relevance to experimental systems of quantum fluids. We demonstrate how it may be used to compute spatial mode entanglement, particle partitioned entanglement, and the entanglement of particles, providing insights into quantum correlations generated by fluctuations, indistinguishability, and interactions. We present proof-of-principle calculations and benchmark against an exactly soluble model of interacting bosons in one spatial dimension. As this algorithm retains the fundamental polynomial scaling of quantum Monte Carlo when applied to sign-problem-free models, future applications should allow for the study of entanglement entropy in large-scale many-body systems of interacting bosons.

arXiv:1404.7104•Physical Review E•Jul 25, 2014•44 citations

Measuring Renyi Entanglement Entropy in Quantum Monte Carlo Simulations

Matthew B. Hastings, Iván González, Ann B. Kallin +1 more

We develop a quantum Monte Carlo procedure, in the valence bond basis, to measure the Renyi entanglement entropy of a many-body ground state as the expectation value of a unitary Swap operator acting on two copies of the system. An improved estimator involving the ratio of Swap operators for different subregions enables convergence of the entropy in a simulation time polynomial in the system size. We demonstrate convergence of the Renyi entropy to exact results for a Heisenberg chain. Finally, we calculate the scaling of the Renyi entropy in the two-dimensional Heisenberg model and confirm that the Néel ground state obeys the expected area law for systems up to linear size L=32.

arXiv:1001.2335•Physical Review Letters•Apr 14, 2010•427 citations

Probing universality with entanglement entropy via quantum Monte Carlo

Bohdan Kulchytskyy

Our understanding of physical phenomena hinges on finding universal core mechanisms that unite them. The concept of universality is deeply ingrained in the study of quantum many-body systems. At zero temperature, microscopically different systems with long-range order collapse into their universal state described by a handful of universal parameters. Establishing those parameters implies identification of the universal theory that is effectively describing the system and ultimately providing the desired understanding. We take up this task with the help of Renyi entanglement entropy. Defined with respect to a system bipartition, this measure quantifies information shared between the subsystems. Understanding what insights into universality are encoded within this information-theoretic quantity as well as developing numerical tools to efficiently estimate the Renyi entanglement entropy are the subjects of this thesis. \n \nOn the computational end of this far-reaching goal, we develop a novel theoretical framework for constructing improved Renyi entanglement entropy estimators in the context of d+1 quantum Monte Carlo methods. The discovery of a connection of this methodology to the well-established Kandel-Domany formalism provides a clear path towards generalization. Additionally, we embrace a data-driven approach towards learning the ground state wavefunction. We demonstrate how a restricted Boltzmann machine can be used to reconstruct the Renyi entanglement entropy from projective measurements of a quantum ground state. Furthermore, we extend this classical generative architecture to a quantum analogue that we call the quantum Boltzmann machine. \n \nOn the theoretical side of this endeavour, we study the Renyi entanglement entropy scaling terms for two quantum lattice models embedded in two dimensional space via extensive quantum Monte Carlo simulations. For the ground state of the XY model, we provide conclusive numerical evidence for a logarithmic contribution that uniquely characterizes the continuous symmetry of the emerging order parameter. Moreover, we confirm the form of the subleading universal geometric contribution arising due to the bosonic nature of low-energy degrees of freedom in this model. For the critical ground state of the transverse field Ising model, we develop a novel scaling procedure to extract a universal number κ revealed via a cylindrical entangling bipartition in the thin-slice limit. The combined product of our work sheds new light on the entanglement-based classification of universality and brings a suite of new powerful numerical tools to continue illuminating this theoretical program in the future.

•Aug 30, 2019•9 citations

Spatial entanglement entropy in the ground state of the Lieb-Liniger model

C. M. Herdman, Pierre–Nicholas Roy, Roger G. Melko +1 more

We consider the entanglement between two spatial subregions in the Lieb-Liniger model of bosons in one spatial dimension interacting via a contact interaction. Using ground-state path integral quantum Monte Carlo we numerically compute the R\'enyi entropy of the reduced density matrix of the subsystem as a measure of entanglement. Our numerical algorithm is based on a replica method previously introduced by the authors, which we extend to efficiently study the entanglement of spatial subsystems of itinerant bosons. We confirm a logarithmic scaling of the R\'enyi entropy with subsystem size that is expected from conformal field theory, and compute the nonuniversal subleading constant for interaction strengths ranging over two orders of magnitude. In the strongly interacting limit, we find agreement with the known free fermion result.

arXiv:1606.03396•Physical review. B./Physical review. B•Aug 31, 2016•13 citations

Detecting symmetry-resolved entanglement: A quantum Monte Carlo approach

Kuangjie Chen, Weizhen Jia, Xiaopeng Li +2 more

Symmetry and entanglement are two fundamental concepts in quantum many-body physics. Their interplay is captured by symmetry-resolved entanglement, which decomposes the total entanglement into contributions from different symmetry sectors. Computing symmetry-resolved entanglement in strongly interacting higher-dimensional quantum systems remains challenging. Here, we formulate and implement an estimator-based quantum Monte Carlo (QMC) framework for computing symmetry-resolved Rényi entropies (SRRE) in sign-problem-free interacting lattice systems by measuring disorder (symmetry-twisted) operators in ordinary and replica ensembles and reconstructing SRRE from the corresponding charged moments. We validate the framework in two controlled one-dimensional settings: the transverse-field Ising model (TFIM), for which exact conformal-field-theory predictions are available, and the interacting Heisenberg chain, which tests the U(1) U ( 1 ) symmetry-sector reconstruction and its finite-size behavior. We then apply the method to the two-dimensional TFIM. Within the accessible system sizes and a phenomenological finite-size extrapolation, our data provide numerical evidence consistent with entanglement equipartition at the (2+1) ( 2 + 1 ) D Ising critical point. Our work establishes a practical numerical route to symmetry-resolved entanglement in interacting lattice models and provides a framework for future studies beyond one dimension.

arXiv:2604.02307•SciPost Physics•Oct 21, 2026•0 citations

Quantum Monte Carlo calculation of entanglement Rényi entropies for generic quantum systems

Stephan Humeniuk, Tommaso Roscilde

We present a general scheme for the calculation of the R\'enyi entropy of a subsystem in quantum many-body models that can be efficiently simulated via quantum Monte Carlo. When the simulation is performed at very low temperature, the above approach delivers the entanglement R\'enyi entropy of the subsystem, and it allows us to explore the crossover to the thermal R\'enyi entropy as the temperature is increased. We implement this scheme explicitly within the stochastic series expansion as well as within path-integral Monte Carlo, and apply it to quantum spin and quantum rotor models. In the case of quantum spins, we show that relevant models in two dimensions with reduced symmetry ($XX$ model or hard-core bosons, transverse-field Ising model at the quantum critical point) exhibit an area law for the scaling of the entanglement entropy.

arXiv:1203.5752•Physical Review B•Dec 12, 2012•125 citations
A ground state path integral quantum Monte Carlo algorithm is introduced that allows for the study of entanglement in lattice bosons at zero temperature. The Rényi entanglement entropy between spatial subregions is explored across the phase diagram of the one dimensional Bose-Hubbard model for systems consisting of up to L=256 L = 256 sites at unit-filling without any restrictions on site occupancy, far beyond the reach of exact diagonalization. The favorable scaling of the algorithm is demonstrated through a further measurement of the Rényi entanglement entropy at the two dimensional superfluid-insulator critical point for large system sizes, confirming the existence of the expected entanglement boundary law in the ground state. The Rényi estimator is extended to measure the symmetry resolved entanglement that is operationally accessible as a resource for experimentally relevant lattice gases with fixed total particle number.
We introduce a quantum Monte Carlo algorithm to measure the Rényi entanglement entropies in systems of interacting bosons in the continuum. This approach is based on a path-integral ground state method that can be applied to interacting itinerant bosons in any spatial dimension with direct relevance to experimental systems of quantum fluids. We demonstrate how it may be used to compute spatial mode entanglement, particle partitioned entanglement, and the entanglement of particles, providing insights into quantum correlations generated by fluctuations, indistinguishability, and interactions. We present proof-of-principle calculations and benchmark against an exactly soluble model of interacting bosons in one spatial dimension. As this algorithm retains the fundamental polynomial scaling of quantum Monte Carlo when applied to sign-problem-free models, future applications should allow for the study of entanglement entropy in large-scale many-body systems of interacting bosons.
We develop a quantum Monte Carlo procedure, in the valence bond basis, to measure the Renyi entanglement entropy of a many-body ground state as the expectation value of a unitary Swap operator acting on two copies of the system. An improved estimator involving the ratio of Swap operators for different subregions enables convergence of the entropy in a simulation time polynomial in the system size. We demonstrate convergence of the Renyi entropy to exact results for a Heisenberg chain. Finally, we calculate the scaling of the Renyi entropy in the two-dimensional Heisenberg model and confirm that the Néel ground state obeys the expected area law for systems up to linear size L=32.
Our understanding of physical phenomena hinges on finding universal core mechanisms that unite them. The concept of universality is deeply ingrained in the study of quantum many-body systems. At zero temperature, microscopically different systems with long-range order collapse into their universal state described by a handful of universal parameters. Establishing those parameters implies identification of the universal theory that is effectively describing the system and ultimately providing the desired understanding. We take up this task with the help of Renyi entanglement entropy. Defined with respect to a system bipartition, this measure quantifies information shared between the subsystems. Understanding what insights into universality are encoded within this information-theoretic quantity as well as developing numerical tools to efficiently estimate the Renyi entanglement entropy are the subjects of this thesis. \n \nOn the computational end of this far-reaching goal, we develop a novel theoretical framework for constructing improved Renyi entanglement entropy estimators in the context of d+1 quantum Monte Carlo methods. The discovery of a connection of this methodology to the well-established Kandel-Domany formalism provides a clear path towards generalization. Additionally, we embrace a data-driven approach towards learning the ground state wavefunction. We demonstrate how a restricted Boltzmann machine can be used to reconstruct the Renyi entanglement entropy from projective measurements of a quantum ground state. Furthermore, we extend this classical generative architecture to a quantum analogue that we call the quantum Boltzmann machine. \n \nOn the theoretical side of this endeavour, we study the Renyi entanglement entropy scaling terms for two quantum lattice models embedded in two dimensional space via extensive quantum Monte Carlo simulations. For the ground state of the XY model, we provide conclusive numerical evidence for a logarithmic contribution that uniquely characterizes the continuous symmetry of the emerging order parameter. Moreover, we confirm the form of the subleading universal geometric contribution arising due to the bosonic nature of low-energy degrees of freedom in this model. For the critical ground state of the transverse field Ising model, we develop a novel scaling procedure to extract a universal number κ revealed via a cylindrical entangling bipartition in the thin-slice limit. The combined product of our work sheds new light on the entanglement-based classification of universality and brings a suite of new powerful numerical tools to continue illuminating this theoretical program in the future.
We consider the entanglement between two spatial subregions in the Lieb-Liniger model of bosons in one spatial dimension interacting via a contact interaction. Using ground-state path integral quantum Monte Carlo we numerically compute the R\'enyi entropy of the reduced density matrix of the subsystem as a measure of entanglement. Our numerical algorithm is based on a replica method previously introduced by the authors, which we extend to efficiently study the entanglement of spatial subsystems of itinerant bosons. We confirm a logarithmic scaling of the R\'enyi entropy with subsystem size that is expected from conformal field theory, and compute the nonuniversal subleading constant for interaction strengths ranging over two orders of magnitude. In the strongly interacting limit, we find agreement with the known free fermion result.
Symmetry and entanglement are two fundamental concepts in quantum many-body physics. Their interplay is captured by symmetry-resolved entanglement, which decomposes the total entanglement into contributions from different symmetry sectors. Computing symmetry-resolved entanglement in strongly interacting higher-dimensional quantum systems remains challenging. Here, we formulate and implement an estimator-based quantum Monte Carlo (QMC) framework for computing symmetry-resolved Rényi entropies (SRRE) in sign-problem-free interacting lattice systems by measuring disorder (symmetry-twisted) operators in ordinary and replica ensembles and reconstructing SRRE from the corresponding charged moments. We validate the framework in two controlled one-dimensional settings: the transverse-field Ising model (TFIM), for which exact conformal-field-theory predictions are available, and the interacting Heisenberg chain, which tests the U(1) U ( 1 ) symmetry-sector reconstruction and its finite-size behavior. We then apply the method to the two-dimensional TFIM. Within the accessible system sizes and a phenomenological finite-size extrapolation, our data provide numerical evidence consistent with entanglement equipartition at the (2+1) ( 2 + 1 ) D Ising critical point. Our work establishes a practical numerical route to symmetry-resolved entanglement in interacting lattice models and provides a framework for future studies beyond one dimension.
We present a general scheme for the calculation of the R\'enyi entropy of a subsystem in quantum many-body models that can be efficiently simulated via quantum Monte Carlo. When the simulation is performed at very low temperature, the above approach delivers the entanglement R\'enyi entropy of the subsystem, and it allows us to explore the crossover to the thermal R\'enyi entropy as the temperature is increased. We implement this scheme explicitly within the stochastic series expansion as well as within path-integral Monte Carlo, and apply it to quantum spin and quantum rotor models. In the case of quantum spins, we show that relevant models in two dimensions with reduced symmetry (XXXXXX model or hard-core bosons, transverse-field Ising model at the quantum critical point) exhibit an area law for the scaling of the entanglement entropy.