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:1307.1486 · Physical Review Letters

Entanglement of Interacting Fermions in Quantum Monte Carlo Calculations

Tarun Grover

Published Sep 24, 2013• 119 citations
arXivHTMLPDF

Abstract

Given a specific interacting quantum Hamiltonian in a general spatial dimension, can one access its entanglement properties, such as the entanglement entropy corresponding to the ground state wave function? Even though progress has been made in addressing this question for interacting bosons and quantum spins, as yet there exist no corresponding methods for interacting fermions. Here we show that the entanglement structure of interacting fermionic Hamiltonians has a particularly simple form-the interacting reduced density matrix can be written as a sum of operators that describe free fermions. This decomposition allows one to calculate the Renyi entropies for Hamiltonians which can be simulated via determinantal quantum Monte Carlo calculations, while employing the efficient techniques hitherto available only for free fermions. The method presented works for the ground state, as well as for the thermally averaged reduced density matrix.

DOI
10.1103/physrevlett.111.130402

Related papers

Rényi entropies of interacting fermions from determinantal quantum Monte Carlo simulations

Peter Broecker, Simon Trebst

Entanglement measures such as the entanglement entropy have become an indispensable tool to identify the fundamental character of ground states of interacting quantum many-body systems. For systems of interacting spin or bosonic degrees of freedom much recent progress has been made not only in the analytical description of their respective entanglement entropies but also in their numerical classification. Systems of interacting fermionic degrees of freedom, however, have proved to be more difficult to control, in particular with regard to the numerical understanding of their entanglement properties. Here we report a generalization of the replica technique for the calculation of Rényi entropies to the framework of determinantal Quantum Monte Carlo simulations—the numerical method of choice for unbiased, large-scale simulations of interacting fermionic systems. We demonstrate the strength of this approach over a recent alternative proposal based on a decomposition in free fermion Green's functions by studying the entanglement entropy of one-dimensional Hubbard systems both at zero and finite temperatures.

arXiv:1404.3027•Journal of Statistical Mechanics Theory and Experiment•Aug 22, 2014•56 citations

Entanglement Rényi negativity of interacting fermions from quantum Monte Carlo simulations

Fo-Hong Wang, Xiao Yan Xu

Many-body entanglement unveils additional aspects of quantum matter and offers insights into strongly correlated physics. While ground-state entanglement has received much attention in the past decade, the study of mixed-state quantum entanglement using negativity in interacting fermionic systems remains largely unexplored. We demonstrate that the partially transposed density matrix of interacting fermions, similar to their reduced density matrix, can be expressed as a weighted sum of Gaussian states describing free fermions, enabling the calculation of rank-n Rényi negativity within the determinant quantum Monte Carlo framework. We calculate the rank-two Rényi negativity for the half-filled Hubbard model and the spinless t-V model. Our calculation reveals that the area law coefficient of the Rényi negativity for the spinless t-V model has a logarithmic finite-size scaling at the finite-temperature transition point. Our work contributes to the calculation of entanglement and sets the stage for future studies on quantum entanglement in various fermionic many-body mixed states.

arXiv:2312.14155•Nature Communications•Mar 17, 2025•14 citations

Renyi Entanglement Entropy of Interacting Fermions Calculated Using the Continuous-Time Quantum Monte Carlo Method

Lei Wang, Matthias Troyer

We present a new algorithm for calculating the Renyi entanglement entropy of interacting fermions using the continuous-time quantum Monte Carlo method. The algorithm only samples the interaction correction of the entanglement entropy, which by design ensures the efficient calculation of weakly interacting systems. Combined with Monte Carlo reweighting, the algorithm also performs well for systems with strong interactions. We demonstrate the potential of this method by studying the quantum entanglement signatures of the charge-density-wave transition of interacting fermions on a square lattice.

arXiv:1407.0707•Physical Review Letters•Sep 10, 2014•51 citations

A scaling function for the particle entanglement entropy of fermions

Harini Radhakrishnan, Matthias Thamm, Hatem Barghathi +2 more

Abstract Entanglement entropy under a particle bipartition generates complementary information to mode entanglement because it is sensitive to interactions and particle statistics at leading order and does not depend on any externally imposed length scale. In this paper, we investigate the particle entanglement entropy in a system of N interacting spinless lattice fermions in one spatial dimension by combining bosonization techniques with exact and approximate numerical methods. We introduce a general scaling form for the fermionic particle entanglement entropy captured by a shape function that enters as an extensive interaction-induced correction to a known free fermion result. A general asymptotic expansion in the total number of particles demonstrates that its form is robust for different values of the Rényi index and highlights how quantum correlations are encoded in the n-particle density matrix of a pure many-body quantum state.

arXiv:2302.09093•Journal of Statistical Mechanics Theory and Experiment•Aug 1, 2023•1 citation

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

Spatial Entanglement of Fermions in One-Dimensional Quantum Dots

Ivan P. Christov

The time-dependent quantum Monte Carlo method for fermions is introduced and applied in the calculation of the entanglement of electrons in one-dimensional quantum dots with several spin-polarized and spin-compensated electron configurations. The rich statistics of wave functions provided by this method allow one to build reduced density matrices for each electron, and to quantify the spatial entanglement using measures such as quantum entropy by treating the electrons as identical or distinguishable particles. Our results indicate that the spatial entanglement in parallel-spin configurations is rather small, and is determined mostly by the spatial quantum nonlocality introduced by the ground state. By contrast, in the spin-compensated case, the outermost opposite-spin electrons interact like bosons, which prevails their entanglement, while the inner-shell electrons remain largely at their Hartree-Fock geometry. Our findings are in close correspondence with the numerically exact results, wherever such comparison is possible.

arXiv:2107.03834•Entropy•Jul 7, 2021•6 citations
Given a specific interacting quantum Hamiltonian in a general spatial dimension, can one access its entanglement properties, such as the entanglement entropy corresponding to the ground state wave function? Even though progress has been made in addressing this question for interacting bosons and quantum spins, as yet there exist no corresponding methods for interacting fermions. Here we show that the entanglement structure of interacting fermionic Hamiltonians has a particularly simple form-the interacting reduced density matrix can be written as a sum of operators that describe free fermions. This decomposition allows one to calculate the Renyi entropies for Hamiltonians which can be simulated via determinantal quantum Monte Carlo calculations, while employing the efficient techniques hitherto available only for free fermions. The method presented works for the ground state, as well as for the thermally averaged reduced density matrix.
Entanglement measures such as the entanglement entropy have become an indispensable tool to identify the fundamental character of ground states of interacting quantum many-body systems. For systems of interacting spin or bosonic degrees of freedom much recent progress has been made not only in the analytical description of their respective entanglement entropies but also in their numerical classification. Systems of interacting fermionic degrees of freedom, however, have proved to be more difficult to control, in particular with regard to the numerical understanding of their entanglement properties. Here we report a generalization of the replica technique for the calculation of Rényi entropies to the framework of determinantal Quantum Monte Carlo simulations—the numerical method of choice for unbiased, large-scale simulations of interacting fermionic systems. We demonstrate the strength of this approach over a recent alternative proposal based on a decomposition in free fermion Green's functions by studying the entanglement entropy of one-dimensional Hubbard systems both at zero and finite temperatures.
Many-body entanglement unveils additional aspects of quantum matter and offers insights into strongly correlated physics. While ground-state entanglement has received much attention in the past decade, the study of mixed-state quantum entanglement using negativity in interacting fermionic systems remains largely unexplored. We demonstrate that the partially transposed density matrix of interacting fermions, similar to their reduced density matrix, can be expressed as a weighted sum of Gaussian states describing free fermions, enabling the calculation of rank-n Rényi negativity within the determinant quantum Monte Carlo framework. We calculate the rank-two Rényi negativity for the half-filled Hubbard model and the spinless t-V model. Our calculation reveals that the area law coefficient of the Rényi negativity for the spinless t-V model has a logarithmic finite-size scaling at the finite-temperature transition point. Our work contributes to the calculation of entanglement and sets the stage for future studies on quantum entanglement in various fermionic many-body mixed states.
We present a new algorithm for calculating the Renyi entanglement entropy of interacting fermions using the continuous-time quantum Monte Carlo method. The algorithm only samples the interaction correction of the entanglement entropy, which by design ensures the efficient calculation of weakly interacting systems. Combined with Monte Carlo reweighting, the algorithm also performs well for systems with strong interactions. We demonstrate the potential of this method by studying the quantum entanglement signatures of the charge-density-wave transition of interacting fermions on a square lattice.
Abstract Entanglement entropy under a particle bipartition generates complementary information to mode entanglement because it is sensitive to interactions and particle statistics at leading order and does not depend on any externally imposed length scale. In this paper, we investigate the particle entanglement entropy in a system of N interacting spinless lattice fermions in one spatial dimension by combining bosonization techniques with exact and approximate numerical methods. We introduce a general scaling form for the fermionic particle entanglement entropy captured by a shape function that enters as an extensive interaction-induced correction to a known free fermion result. A general asymptotic expansion in the total number of particles demonstrates that its form is robust for different values of the Rényi index and highlights how quantum correlations are encoded in the n-particle density matrix of a pure many-body quantum state.
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.
The time-dependent quantum Monte Carlo method for fermions is introduced and applied in the calculation of the entanglement of electrons in one-dimensional quantum dots with several spin-polarized and spin-compensated electron configurations. The rich statistics of wave functions provided by this method allow one to build reduced density matrices for each electron, and to quantify the spatial entanglement using measures such as quantum entropy by treating the electrons as identical or distinguishable particles. Our results indicate that the spatial entanglement in parallel-spin configurations is rather small, and is determined mostly by the spatial quantum nonlocality introduced by the ground state. By contrast, in the spin-compensated case, the outermost opposite-spin electrons interact like bosons, which prevails their entanglement, while the inner-shell electrons remain largely at their Hartree-Fock geometry. Our findings are in close correspondence with the numerically exact results, wherever such comparison is possible.