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:1203.5752 · Physical Review B

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

Stephan Humeniuk, Tommaso Roscilde

Published Dec 12, 2012• 125 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevb.86.235116

Related papers

Wang-Landau method for calculating Rényi entropies in finite-temperature quantum Monte Carlo simulations

Stephen Inglis, Roger G. Melko

We implement a Wang-Landau sampling technique in quantum Monte Carlo (QMC) simulations for the purpose of calculating the Rényi entanglement entropies and associated mutual information. The algorithm converges an estimate for an analog to the density of states for stochastic series expansion QMC, allowing a direct calculation of Rényi entropies without explicit thermodynamic integration. We benchmark results for the mutual information on two-dimensional (2D) isotropic and anisotropic Heisenberg models, a 2D transverse field Ising model, and a three-dimensional Heisenberg model, confirming a critical scaling of the mutual information in cases with a finite-temperature transition. We discuss the benefits and limitations of broad sampling techniques compared to standard importance sampling methods.

arXiv:1207.5052•Physical Review E•Jan 22, 2013•49 citations

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

Nonequilibrium quantum Monte Carlo algorithm for stabilizer Rényi entropy in spin systems

Zejun Liu, Bryan K. Clark

Quantum magic, or nonstabilizerness, provides a crucial characterization of quantum systems, regarding the classical simulability with stabilizer states. In this work, we propose an alternative and efficient algorithm for computing stabilizer R\'enyi entropy, one of the measures for quantum magic, in spin systems with sign-problem free Hamiltonians. This algorithm is based on the quantum Monte Carlo simulation of the path integral of the work between two partition function ensembles and it applies to all spatial dimensions and temperatures. We demonstrate this algorithm on the one- and two-dimensional transverse field Ising model at both finite and zero temperatures and show the quantitative agreements with tensor-network based algorithms. We analyze the computational cost and provide analytical and numerical evidences for it to be polynomial in system size. This work also suggests a unifying framework for calculating various types of entropy quantities including entanglement R\'enyi entropy and entanglement R\'enyi negativity.

arXiv:2405.19577•Physical review. B./Physical review. B•Feb 20, 2025•17 citations

Implementations and applications of Renyi entanglement in Monte Carlo simulations of spin models

Stephen Inglis

Although entanglement is a well studied property in the context of quantum systems, the ability to measure it in Monte Carlo methods is relatively new. \nThrough measures of the Renyi entanglement entropy and mutual information one is able to examine and characterize criticality, pinpoint phase transitions, and probe universality. \nWe describe the most basic algorithms for calculating these quantities in straightforward Monte Carlo methods and state of the art techniques used in high performance computing. \nThis description emphasizes the core principal of these measurements and allows one to both build an intuition for these quantities and how they are useful in numerical studies. \n \nUsing the Renyi entanglement entropy we demonstrate the ability to detect thermal phase transitions in the Ising model and XY model without use of an order parameter. \nThe scaling near the critical point also shows signatures identifying the universality class of the model. \nImproved methods are explored using extended ensemble techniques that can increase calculation efficiency, and show good agreement with the standard approach. \nWe explore the "ratio trick" at finite temperature and use it to explore the quantum critical fan of the one dimensional transverse field Ising model, showing agreement with finite temperature and finite size scaling from field theory. \nThis same technique is used at zero temperature to explore the geometric dependence of the entanglement entropy and examine the universal scaling functions in the two dimensional transverse field Ising model. \nAll of this shows the multitude of ways in which the study of the Renyi entanglement entropy can be efficiently and practically used in conventional and exotic condensed matter systems, and should serve as a reference for those wishing to use it as a tool.

•Sep 27, 2013•8 citations

Measuring Rényi entanglement entropy with high efficiency and precision in quantum Monte Carlo simulations

Jiarui Zhao, Bin-Bin Chen, Yan-Cheng Wang +3 more

Abstract We develop a nonequilibrium increment method in quantum Monte Carlo simulations to obtain the Rényi entanglement entropy of various quantum many-body systems with high efficiency and precision. To demonstrate its power, we show the results on a few important yet difficult (2 + 1)dquantum lattice models, ranging from the Heisenberg quantum antiferromagnet with spontaneous symmetry breaking, the quantum critical point with O(3) conformal field theory (CFT) to the toric code $${{\mathbb{Z}}}_{2}$$ Z2 topological ordered state and the Kagome $${{\mathbb{Z}}}_{2}$$ Z2 quantum spin liquid model with frustration and multi-spin interactions. In all these cases, our method either reveals the precise CFT data from the logarithmic correction or extracts the quantum dimension in topological order, from the dominant area law in finite-size scaling, with very large system sizes, controlled errorbars, and minimal computational costs. Our method, therefore, establishes a controlled and practical computation paradigm to obtain the difficult yet important universal properties in highly entangled quantum matter.

arXiv:2112.15178•npj Quantum Materials•Jun 29, 2022•61 citations

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

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

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.

arXiv:2207.11301•SciPost Physics•Mar 29, 2023•6 citations
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.
We implement a Wang-Landau sampling technique in quantum Monte Carlo (QMC) simulations for the purpose of calculating the Rényi entanglement entropies and associated mutual information. The algorithm converges an estimate for an analog to the density of states for stochastic series expansion QMC, allowing a direct calculation of Rényi entropies without explicit thermodynamic integration. We benchmark results for the mutual information on two-dimensional (2D) isotropic and anisotropic Heisenberg models, a 2D transverse field Ising model, and a three-dimensional Heisenberg model, confirming a critical scaling of the mutual information in cases with a finite-temperature transition. We discuss the benefits and limitations of broad sampling techniques compared to standard importance sampling methods.
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.
Quantum magic, or nonstabilizerness, provides a crucial characterization of quantum systems, regarding the classical simulability with stabilizer states. In this work, we propose an alternative and efficient algorithm for computing stabilizer R\'enyi entropy, one of the measures for quantum magic, in spin systems with sign-problem free Hamiltonians. This algorithm is based on the quantum Monte Carlo simulation of the path integral of the work between two partition function ensembles and it applies to all spatial dimensions and temperatures. We demonstrate this algorithm on the one- and two-dimensional transverse field Ising model at both finite and zero temperatures and show the quantitative agreements with tensor-network based algorithms. We analyze the computational cost and provide analytical and numerical evidences for it to be polynomial in system size. This work also suggests a unifying framework for calculating various types of entropy quantities including entanglement R\'enyi entropy and entanglement R\'enyi negativity.
Although entanglement is a well studied property in the context of quantum systems, the ability to measure it in Monte Carlo methods is relatively new. \nThrough measures of the Renyi entanglement entropy and mutual information one is able to examine and characterize criticality, pinpoint phase transitions, and probe universality. \nWe describe the most basic algorithms for calculating these quantities in straightforward Monte Carlo methods and state of the art techniques used in high performance computing. \nThis description emphasizes the core principal of these measurements and allows one to both build an intuition for these quantities and how they are useful in numerical studies. \n \nUsing the Renyi entanglement entropy we demonstrate the ability to detect thermal phase transitions in the Ising model and XY model without use of an order parameter. \nThe scaling near the critical point also shows signatures identifying the universality class of the model. \nImproved methods are explored using extended ensemble techniques that can increase calculation efficiency, and show good agreement with the standard approach. \nWe explore the "ratio trick" at finite temperature and use it to explore the quantum critical fan of the one dimensional transverse field Ising model, showing agreement with finite temperature and finite size scaling from field theory. \nThis same technique is used at zero temperature to explore the geometric dependence of the entanglement entropy and examine the universal scaling functions in the two dimensional transverse field Ising model. \nAll of this shows the multitude of ways in which the study of the Renyi entanglement entropy can be efficiently and practically used in conventional and exotic condensed matter systems, and should serve as a reference for those wishing to use it as a tool.
Abstract We develop a nonequilibrium increment method in quantum Monte Carlo simulations to obtain the Rényi entanglement entropy of various quantum many-body systems with high efficiency and precision. To demonstrate its power, we show the results on a few important yet difficult (2 + 1)dquantum lattice models, ranging from the Heisenberg quantum antiferromagnet with spontaneous symmetry breaking, the quantum critical point with O(3) conformal field theory (CFT) to the toric code Z2{{\mathbb{Z}}}_{2}Z2​ Z2 topological ordered state and the Kagome Z2{{\mathbb{Z}}}_{2}Z2​ Z2 quantum spin liquid model with frustration and multi-spin interactions. In all these cases, our method either reveals the precise CFT data from the logarithmic correction or extracts the quantum dimension in topological order, from the dominant area law in finite-size scaling, with very large system sizes, controlled errorbars, and minimal computational costs. Our method, therefore, establishes a controlled and practical computation paradigm to obtain the difficult yet important universal properties in highly entangled quantum matter.
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.