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:2112.14199 · Physical Review Research

Quantum Rényi entropy by optimal thermodynamic integration paths

Miha Srdinšek, Michele Casula, Rodolphe Vuilleumier

Published Jul 5, 2022• 7 citations
arXivHTMLPDF

Abstract

Despite being a well-established operational approach to quantify entanglement, Rényi entropy calculations have been plagued by their computational complexity.We introduce here a theoretical framework based on an optimal thermodynamic integration scheme, where the Rényi entropy can be efficiently evaluated using regularizing paths.This approach avoids slowly convergent fluctuating contributions and leads to low-variance estimates.In this way, large system sizes and high levels of entanglement in model or first-principles Hamiltonians are within our reach.We demonstrate this approach in the one-dimensional quantum Ising model and perform an evaluation of entanglement entropy in the formic acid dimer, by discovering that its two shared protons are entangled even above room temperature.

DOI
10.1103/physrevresearch.4.l032002

Related papers

Numerical evaluation of Rényi entropy with path integrals

Miha Srdinšek

Quantum phase transitions and quantum entanglement can be studied by quantum Rényi entropy, a generalization of Von Neumann entropy. However, Rényi entropy calculations have been plagued by their computational complexity. In this thesis we present a new numerical technique based on thermodynamic integration via regularized integration paths. The optimal path is constructed by minimizing the variance of the Rényi entropy estimator leading to an efficient evaluation of its quantity for complex quantum systems of distinguishable particles with unprecedented entanglement and size. We implement the method in path integral Monte Carlo (PIMC) and path integral molecular dynamics (PIMD) schemes. We demonstrate in the case of the one-dimensional Ising model in transverse magnetic field, that the method outperforms some of the existing state-of-the-art numerical methods. We then study the complex of formic acid dimer and discover that its two shared protons are entangled even above room temperature. By the means of hybrid PIMC-PIMD scheme we then fully resolve the phase diagram of the discrete ϕ4 model, and show that the peak of the Rényi entropy at finite temperature exactly coincides with the quantum critical point (QCP) and the freezing of the field in the imaginary time, interpreted as the local moments formation. We then use the Rényi entropy to determine the phase diagram of high-pressure water ice phases VIII, VII and X. The quantum Rényi entropy of protons sharply peaks at the pressure of 43 GPa, which is a signature of a QCP. The peak precisely coincides with the freezing of protons in the imaginary time, where phase X is lost, and the formation of a macroscopic polarization at lower temperature, characteristic of phase VIII. In this way we are able to scan the full phase diagram of high pressure water ice by looking at the changes in Rényi entropy behavior as a function of pressure and temperature.

•Sep 26, 2023•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

The Rényi entanglement entropy of a general quantum dimer model at the RK point: a highly efficient algorithm

Jiquan Pei, Steve Han, Haijun Liao +1 more

A highly efficient and simple-to-implement Monte Carlo algorithm is proposed for the evaluation of the Rényi entanglement entropy (REE) of the quantum dimer model (QDM) at the Rokhsar-Kivelson (RK) point. It makes possible the evaluation of REE at the RK point to the thermodynamic limit for a general QDM. We apply the algorithm to a QDM defined on the triangular and the square lattice in two dimensions and the simple and the face centered cubic (fcc) lattice in three dimensions. We find the REE on all these lattices follows perfect linear scaling in the thermodynamic limit, apart from an even-odd oscillation in the case of the square lattice. We also evaluate the topological entanglement entropy (TEE) with both a subtraction and an extrapolation procedure. We find the QDMs on both the triangular and the fcc lattice exhibit robust Z2 topological order. The expected TEE of ln2 is clearly demonstrated in both cases. Our large scale simulation also proves the recently proposed extrapolation procedure in cylindrical geometry to be a highly reliable way to extract the TEE of a topologically ordered system.

arXiv:1211.1504•Journal of Physics Condensed Matter•Dec 12, 2013•8 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

Rényi entropy of a quantum anharmonic chain at nonzero temperature

Miha Srdinšek, Michele Casula, Rodolphe Vuilleumier

The interplay of quantum and classical fluctuations in the vicinity of a quantum critical point (QCP) gives rise to various regimes or phases with distinct quantum character. In this work, we show that the R\'enyi entropy is a precious tool to characterize the phase diagram of critical systems not only around the QCP but also away from it, thanks to its capability to detect the emergence of local moments at finite temperature. For an efficient evaluation of the R\'enyi entropy, we introduce an algorithm based on a path-integral Langevin dynamics combined with a previously proposed thermodynamic integration method built on regularized paths. We apply this framework to study the critical behavior of a linear chain of anharmonic oscillators, a particular realization of the ${\ensuremath{\phi}}^{4}$ model. We fully resolved its phase diagram, as a function of both temperature and interaction strength. At finite temperature, we find a sequence of three regimes---para, disordered, and antiferro---met as the interaction is increased. The R\'enyi entropy divergence coincides with the crossover between the para and disordered regimes, which shows no temperature dependence. The occurrence of the antiferro regime, on the other hand, is temperature dependent. The two crossover lines merge in proximity of the QCP, at zero temperature, where the R\'enyi entropy is sharply peaked. Via its subsystem-size scaling, we confirm that the transition belongs to the two-dimensional Ising universality class. This phenomenology is expected to happen in all ${\ensuremath{\phi}}^{4}$-like systems, as well as in the elusive water-ice transition across phases VII, VIII, and X.

arXiv:2303.04768•Physical review. B./Physical review. B•Dec 7, 2023•4 citations
Despite being a well-established operational approach to quantify entanglement, Rényi entropy calculations have been plagued by their computational complexity.We introduce here a theoretical framework based on an optimal thermodynamic integration scheme, where the Rényi entropy can be efficiently evaluated using regularizing paths.This approach avoids slowly convergent fluctuating contributions and leads to low-variance estimates.In this way, large system sizes and high levels of entanglement in model or first-principles Hamiltonians are within our reach.We demonstrate this approach in the one-dimensional quantum Ising model and perform an evaluation of entanglement entropy in the formic acid dimer, by discovering that its two shared protons are entangled even above room temperature.
Quantum phase transitions and quantum entanglement can be studied by quantum Rényi entropy, a generalization of Von Neumann entropy. However, Rényi entropy calculations have been plagued by their computational complexity. In this thesis we present a new numerical technique based on thermodynamic integration via regularized integration paths. The optimal path is constructed by minimizing the variance of the Rényi entropy estimator leading to an efficient evaluation of its quantity for complex quantum systems of distinguishable particles with unprecedented entanglement and size. We implement the method in path integral Monte Carlo (PIMC) and path integral molecular dynamics (PIMD) schemes. We demonstrate in the case of the one-dimensional Ising model in transverse magnetic field, that the method outperforms some of the existing state-of-the-art numerical methods. We then study the complex of formic acid dimer and discover that its two shared protons are entangled even above room temperature. By the means of hybrid PIMC-PIMD scheme we then fully resolve the phase diagram of the discrete ϕ4 model, and show that the peak of the Rényi entropy at finite temperature exactly coincides with the quantum critical point (QCP) and the freezing of the field in the imaginary time, interpreted as the local moments formation. We then use the Rényi entropy to determine the phase diagram of high-pressure water ice phases VIII, VII and X. The quantum Rényi entropy of protons sharply peaks at the pressure of 43 GPa, which is a signature of a QCP. The peak precisely coincides with the freezing of protons in the imaginary time, where phase X is lost, and the formation of a macroscopic polarization at lower temperature, characteristic of phase VIII. In this way we are able to scan the full phase diagram of high pressure water ice by looking at the changes in Rényi entropy behavior as a function of pressure and temperature.
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.
A highly efficient and simple-to-implement Monte Carlo algorithm is proposed for the evaluation of the Rényi entanglement entropy (REE) of the quantum dimer model (QDM) at the Rokhsar-Kivelson (RK) point. It makes possible the evaluation of REE at the RK point to the thermodynamic limit for a general QDM. We apply the algorithm to a QDM defined on the triangular and the square lattice in two dimensions and the simple and the face centered cubic (fcc) lattice in three dimensions. We find the REE on all these lattices follows perfect linear scaling in the thermodynamic limit, apart from an even-odd oscillation in the case of the square lattice. We also evaluate the topological entanglement entropy (TEE) with both a subtraction and an extrapolation procedure. We find the QDMs on both the triangular and the fcc lattice exhibit robust Z2 topological order. The expected TEE of ln2 is clearly demonstrated in both cases. Our large scale simulation also proves the recently proposed extrapolation procedure in cylindrical geometry to be a highly reliable way to extract the TEE of a topologically ordered system.
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.
The interplay of quantum and classical fluctuations in the vicinity of a quantum critical point (QCP) gives rise to various regimes or phases with distinct quantum character. In this work, we show that the R\'enyi entropy is a precious tool to characterize the phase diagram of critical systems not only around the QCP but also away from it, thanks to its capability to detect the emergence of local moments at finite temperature. For an efficient evaluation of the R\'enyi entropy, we introduce an algorithm based on a path-integral Langevin dynamics combined with a previously proposed thermodynamic integration method built on regularized paths. We apply this framework to study the critical behavior of a linear chain of anharmonic oscillators, a particular realization of the ${\ensuremath{\phi}}^{4}$ model. We fully resolved its phase diagram, as a function of both temperature and interaction strength. At finite temperature, we find a sequence of three regimes---para, disordered, and antiferro---met as the interaction is increased. The R\'enyi entropy divergence coincides with the crossover between the para and disordered regimes, which shows no temperature dependence. The occurrence of the antiferro regime, on the other hand, is temperature dependent. The two crossover lines merge in proximity of the QCP, at zero temperature, where the R\'enyi entropy is sharply peaked. Via its subsystem-size scaling, we confirm that the transition belongs to the two-dimensional Ising universality class. This phenomenology is expected to happen in all ${\ensuremath{\phi}}^{4}$-like systems, as well as in the elusive water-ice transition across phases VII, VIII, and X.