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

Entanglement spectroscopy using quantum Monte Carlo

Chia-Min Chung, Lars Bonnes, Pochung Chen, Andreas Martin Läuchli

Published May 30, 2014• 36 citations
arXivHTMLPDF

Abstract

We present a numerical scheme to reconstruct a subset of the entanglement spectrum of quantum many body systems using quantum Monte Carlo. The approach builds on the replica trick to evaluate particle number resolved traces of the first $n$ of powers of a reduced density matrix. From this information we reconstruct $n$ entanglement spectrum levels using a polynomial root solver. We illustrate the power and limitations of the method by an application to the extended Bose-Hubbard model in one dimension where we are able to resolve the quasidegeneracy of the entanglement spectrum in the Haldane-insulator phase. In general, the method is able to reconstruct the largest few eigenvalues in each symmetry sector and typically performs better when the eigenvalues are not too different.

DOI
10.1103/physrevb.89.195147

Related papers

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

Sampling reduced density matrix to extract fine levels of entanglement spectrum and restore entanglement Hamiltonian

Bin-Bin Mao, Yi-Ming Ding, Zhe Wang +2 more

The reduced density matrix (RDM) plays a key role in quantum entanglement and measurement, as it allows the extraction of almost all physical quantities related to the reduced degrees of freedom. However, restricted by the degrees of freedom in the environment, the total system size is often limited, let alone the subsystem. To address this challenge, we propose a quantum Monte Carlo scheme with a low technical barrier, enabling precise extraction of the RDM. To demonstrate the power of the method, we present the fine levels of the entanglement spectrum (ES), which is the logarithmic eigenvalues of the RDM. We clearly show the ES for a 1D ladder with a long entangled boundary, and that for the 2D Heisenberg model with a tower of states. Furthermore, we put forward an efficient way to restore the entanglement Hamiltonian in operator-form from the sampled RDM data. Our simulation results, utilizing unprecedentedly large system sizes, establish a practical computational framework for determining entanglement quantities based on the RDM, such as the ES, particularly in scenarios where the environment has a huge number of degrees of freedom. The reduced density matrix enables extraction of entanglement-related quantities but computing it for large quantum systems is challenging. Mao et al. combine quantum Monte Carlo with exact diagonalization to efficiently extract the reduced density matrix and entanglement spectrum in large 1D and 2D systems.

arXiv:2310.16709•Nature Communications•Mar 24, 2025•15 citations

Exploring the limit of the Lattice-Bisognano-Wichmann form describing the Entanglement Hamiltonian: A quantum Monte Carlo study

Siyi Yang, Yi-Ming Ding, Zheng Yan

As a powerful theoretical construct, the entanglement Hamiltonian (EH) encapsulates the essential entanglement properties of a quantum many-body system. From the EH, one can extract a variety of entanglement quantities, such as entanglement entropies, negativity, and the entanglement spectrum. However, its general analytical form remains largely unknown. While the Bisognano-Wichmann theorem gives an exact EH form for Lorentz-invariant field theories, its validity on lattice systems is limited, especially when Lorentz invariance is absent. In this work, we propose a general scheme based on the lattice-Bisognano-Wichmann (LBW) ansatz and multi-replica-trick quantum Monte Carlo methods to numerically reconstruct the entanglement Hamiltonian in two-dimensional systems and systematically explore its applicability to systems without translational invariance, going beyond the original scope of the primordial Bisognano-Wichmann theorem. Various quantum phases--including gapped and gapless phases, critical points, and phases with either discrete or continuous symmetry breaking--are investigated, demonstrating the versatility of our method in reconstructing entanglement Hamiltonians. Furthermore, we find that when the entanglement boundary of a system is ordinary (i.e., free from surface anomalies), the LBW ansatz provides an accurate approximation well beyond Lorentz-invariant cases. Our work thus establishes a general framework for investigating the analytical structure of entanglement in the complex quantum many-body systems.

cond-mat.str-elCross-listedarXiv:2511.00950•Physical review. B./Physical review. B•Nov 4, 2025•0 citations

Calculating the entanglement spectrum in quantum Monte Carlo with application toab initioHamiltonians

Norm M. Tubman, David ChangMo Yang

Several algorithms have been proposed to calculate the spatial entanglement spectrum from high order R\'enyi entropies and maximum entropy techniques. We present an alternative approach for computing the entanglement spectrum with quantum Monte Carlo for both continuum and lattice Hamiltonians. This method provides direct access to the matrix elements of the spatially reduced density matrix and we determine an estimator that can be used in variational Monte Carlo as well as other Monte Carlo methods. The algorithm is based on using a generalization of the swap operator, which can be extended to calculate a general class of density matrices that can include combinations of spin, space, particle, and even momentum degrees of freedom. We demonstrate the method by applying it to the ${\mathrm{H}}_{2}$ and ${\mathrm{N}}_{2}$ molecules and describe how the spatial entanglement spectrum encodes a covalent bond that includes all the many body correlations.

arXiv:1402.0503•Physical Review B•Aug 28, 2014•19 citations

Quantum Monte Carlo in the Age of Many-Body Quantum Information

Yi-Ming Ding, Bin-Bin Mao, Zheng Yan

Quantum Monte Carlo (QMC) methods are among the central numerical tools for studying strongly correlated quantum many-body systems, particularly in higher dimensions. As quantum information has introduced new information-theoretic perspectives and diagnostics into many-body physics, QMC methods have accordingly been extended beyond the measurement of conventional linear observables. This review summarizes recent progress in adapting QMC to many-body quantum-information, focusing on qubit or spin-$1/2$ systems as a concrete setting while keeping the discussion broadly applicable to qudit and bosonic systems. We present a unified perspective on the extraction of nonlinear diagnostics, including entanglement entropies and entanglement spectra, Rényi negativities for mixed-state entanglement, stabilizer entropies for quantum magic, and decoherence-driven phenomena such as the interplay between imaginary-time evolution and decoherence and strong-to-weak spontaneous symmetry breaking.

quant-pharXiv:2608.23231•Aug 25, 2026•0 citations

Entanglement negativity via the replica trick: A quantum Monte Carlo approach

Chia-Min Chung, Vincenzo Alba, Lars Bonnes +2 more

Motivated by recent developments in conformal field theory (CFT), we devise a quantum Monte Carlo (QMC) method to calculate the moments of the partially transposed reduced density matrix at finite temperature. These are used to construct scale invariant combinations that are related to the negativity, a true measure of entanglement for two intervals embedded in a chain. These quantities can serve as witnesses of criticality. In particular, we study several scale invariant combinations of the moments for the one-dimensional (1D) hard-core boson model. For two adjacent intervals unusual finite-size corrections are present, showing parity effects that oscillate with a filling dependent period. These are more pronounced in the presence of boundaries. For large chains we find perfect agreement with CFT. Oppositely, for disjoint intervals corrections are more severe and CFT is recovered only asymptotically. Furthermore, we provide evidence that their exponent is the same as that governing the corrections of the mutual information. Additionally we study the 1D Bose-Hubbard model in the superfluid phase. Remarkably, the finite-size effects are smaller and QMC data are already in impressive agreement with CFT at moderately large sizes.

arXiv:1312.1168•Physical Review B•Aug 4, 2014•76 citations
We present a numerical scheme to reconstruct a subset of the entanglement spectrum of quantum many body systems using quantum Monte Carlo. The approach builds on the replica trick to evaluate particle number resolved traces of the first nnn of powers of a reduced density matrix. From this information we reconstruct nnn entanglement spectrum levels using a polynomial root solver. We illustrate the power and limitations of the method by an application to the extended Bose-Hubbard model in one dimension where we are able to resolve the quasidegeneracy of the entanglement spectrum in the Haldane-insulator phase. In general, the method is able to reconstruct the largest few eigenvalues in each symmetry sector and typically performs better when the eigenvalues are not too different.
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.
The reduced density matrix (RDM) plays a key role in quantum entanglement and measurement, as it allows the extraction of almost all physical quantities related to the reduced degrees of freedom. However, restricted by the degrees of freedom in the environment, the total system size is often limited, let alone the subsystem. To address this challenge, we propose a quantum Monte Carlo scheme with a low technical barrier, enabling precise extraction of the RDM. To demonstrate the power of the method, we present the fine levels of the entanglement spectrum (ES), which is the logarithmic eigenvalues of the RDM. We clearly show the ES for a 1D ladder with a long entangled boundary, and that for the 2D Heisenberg model with a tower of states. Furthermore, we put forward an efficient way to restore the entanglement Hamiltonian in operator-form from the sampled RDM data. Our simulation results, utilizing unprecedentedly large system sizes, establish a practical computational framework for determining entanglement quantities based on the RDM, such as the ES, particularly in scenarios where the environment has a huge number of degrees of freedom. The reduced density matrix enables extraction of entanglement-related quantities but computing it for large quantum systems is challenging. Mao et al. combine quantum Monte Carlo with exact diagonalization to efficiently extract the reduced density matrix and entanglement spectrum in large 1D and 2D systems.
As a powerful theoretical construct, the entanglement Hamiltonian (EH) encapsulates the essential entanglement properties of a quantum many-body system. From the EH, one can extract a variety of entanglement quantities, such as entanglement entropies, negativity, and the entanglement spectrum. However, its general analytical form remains largely unknown. While the Bisognano-Wichmann theorem gives an exact EH form for Lorentz-invariant field theories, its validity on lattice systems is limited, especially when Lorentz invariance is absent. In this work, we propose a general scheme based on the lattice-Bisognano-Wichmann (LBW) ansatz and multi-replica-trick quantum Monte Carlo methods to numerically reconstruct the entanglement Hamiltonian in two-dimensional systems and systematically explore its applicability to systems without translational invariance, going beyond the original scope of the primordial Bisognano-Wichmann theorem. Various quantum phases--including gapped and gapless phases, critical points, and phases with either discrete or continuous symmetry breaking--are investigated, demonstrating the versatility of our method in reconstructing entanglement Hamiltonians. Furthermore, we find that when the entanglement boundary of a system is ordinary (i.e., free from surface anomalies), the LBW ansatz provides an accurate approximation well beyond Lorentz-invariant cases. Our work thus establishes a general framework for investigating the analytical structure of entanglement in the complex quantum many-body systems.
Several algorithms have been proposed to calculate the spatial entanglement spectrum from high order R\'enyi entropies and maximum entropy techniques. We present an alternative approach for computing the entanglement spectrum with quantum Monte Carlo for both continuum and lattice Hamiltonians. This method provides direct access to the matrix elements of the spatially reduced density matrix and we determine an estimator that can be used in variational Monte Carlo as well as other Monte Carlo methods. The algorithm is based on using a generalization of the swap operator, which can be extended to calculate a general class of density matrices that can include combinations of spin, space, particle, and even momentum degrees of freedom. We demonstrate the method by applying it to the H2{\mathrm{H}}_{2}H2​ and N2{\mathrm{N}}_{2}N2​ molecules and describe how the spatial entanglement spectrum encodes a covalent bond that includes all the many body correlations.
Quantum Monte Carlo (QMC) methods are among the central numerical tools for studying strongly correlated quantum many-body systems, particularly in higher dimensions. As quantum information has introduced new information-theoretic perspectives and diagnostics into many-body physics, QMC methods have accordingly been extended beyond the measurement of conventional linear observables. This review summarizes recent progress in adapting QMC to many-body quantum-information, focusing on qubit or spin-1/21/21/2 systems as a concrete setting while keeping the discussion broadly applicable to qudit and bosonic systems. We present a unified perspective on the extraction of nonlinear diagnostics, including entanglement entropies and entanglement spectra, Rényi negativities for mixed-state entanglement, stabilizer entropies for quantum magic, and decoherence-driven phenomena such as the interplay between imaginary-time evolution and decoherence and strong-to-weak spontaneous symmetry breaking.
Motivated by recent developments in conformal field theory (CFT), we devise a quantum Monte Carlo (QMC) method to calculate the moments of the partially transposed reduced density matrix at finite temperature. These are used to construct scale invariant combinations that are related to the negativity, a true measure of entanglement for two intervals embedded in a chain. These quantities can serve as witnesses of criticality. In particular, we study several scale invariant combinations of the moments for the one-dimensional (1D) hard-core boson model. For two adjacent intervals unusual finite-size corrections are present, showing parity effects that oscillate with a filling dependent period. These are more pronounced in the presence of boundaries. For large chains we find perfect agreement with CFT. Oppositely, for disjoint intervals corrections are more severe and CFT is recovered only asymptotically. Furthermore, we provide evidence that their exponent is the same as that governing the corrections of the mutual information. Additionally we study the 1D Bose-Hubbard model in the superfluid phase. Remarkably, the finite-size effects are smaller and QMC data are already in impressive agreement with CFT at moderately large sizes.