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

Entanglement spectra of interacting fermions in quantum Monte Carlo simulations

Fakher F. Assaad, Thomas C. Lang, Francesco Parisen Toldin

Published Mar 24, 2014• 80 citations
arXivHTMLPDF

Abstract

The authors present a detailed method for calculating higher order Renyi entanglement entropies and the entanglement spectrum using auxiliary quantum Monte Carlo. These findings allow the controlled and uncomplicated access to the entanglement properties of correlated fermion systems.
DOI
10.1103/physrevb.89.125121

Related papers

Entanglement studies of interacting fermionic models

Francesco Parisen Toldin, Fakher F. Assaad

Recent advances in the field of strongly correlated electron systems allow to access the entanglement properties of interacting fermionic models, by means of Monte Carlo simulations. We briefly review the techniques used in this context to determine the entanglement entropies and correlations of the entanglement Hamiltonian. We further apply these methods to compute the spin two-point function of entanglement Hamiltonian for a stripe embedded into a correlated topological insulator. Further we discuss a recent method that allows an unbiased, numerically exact, direct determination of the entanglement Hamiltonian by means of auxiliary field quantum Monte Carlo simulations.

arXiv:1810.06595•Journal of Physics Conference Series•Feb 1, 2019•6 citations

Entanglement of Interacting Fermions in Quantum Monte Carlo Calculations

Tarun Grover

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.

arXiv:1307.1486•Physical Review Letters•Sep 24, 2013•119 citations

Stable quantum Monte Carlo simulations for entanglement spectra of interacting fermions

Fakher F. Assaad

We show that the two recently proposed methods to compute Renyi entanglement entropies in the realm of determinant quantum Monte Carlo methods for fermions are in principle equivalent, but differ in sampling strategies. The analogy allows us to formulate a numerically stable calculation of the entanglement spectrum at strong coupling. We demonstrate the approach by studying static and dynamical properties of the entanglement Hamiltonian across the interaction driven quantum phase transition between a topological insulator and quantum antiferromagnet in the Kane-Mele Hubbard model. The formulation is not limited to fermion systems and can readily be adapted to world-line-based simulations of bosonic systems.

arXiv:1501.01418•Physical Review B•Mar 31, 2015•51 citations

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

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

Entanglement Hamiltonian of Interacting Fermionic Models

Francesco Parisen Toldin, Fakher F. Assaad

Recent numerical advances in the field of strongly correlated electron systems allow the calculation of the entanglement spectrum and entropies for interacting fermionic systems. An explicit determination of the entanglement (modular) Hamiltonian has proven to be a considerably more difficult problem, and only a few results are available. We introduce a technique to directly determine the entanglement Hamiltonian of interacting fermionic models by means of auxiliary field quantum Monte Carlo simulations. We implement our method on the one-dimensional Hubbard chain partitioned into two segments and on the Hubbard two-leg ladder partitioned into two chains. In both cases, we study the evolution of the entanglement Hamiltonian as a function of the physical temperature.

arXiv:1804.03163•Physical Review Letters•Nov 14, 2018•68 citations
Recent advances in the field of strongly correlated electron systems allow to access the entanglement properties of interacting fermionic models, by means of Monte Carlo simulations. We briefly review the techniques used in this context to determine the entanglement entropies and correlations of the entanglement Hamiltonian. We further apply these methods to compute the spin two-point function of entanglement Hamiltonian for a stripe embedded into a correlated topological insulator. Further we discuss a recent method that allows an unbiased, numerically exact, direct determination of the entanglement Hamiltonian by means of auxiliary field quantum Monte Carlo simulations.
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.
We show that the two recently proposed methods to compute Renyi entanglement entropies in the realm of determinant quantum Monte Carlo methods for fermions are in principle equivalent, but differ in sampling strategies. The analogy allows us to formulate a numerically stable calculation of the entanglement spectrum at strong coupling. We demonstrate the approach by studying static and dynamical properties of the entanglement Hamiltonian across the interaction driven quantum phase transition between a topological insulator and quantum antiferromagnet in the Kane-Mele Hubbard model. The formulation is not limited to fermion systems and can readily be adapted to world-line-based simulations of bosonic systems.
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.
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.
Recent numerical advances in the field of strongly correlated electron systems allow the calculation of the entanglement spectrum and entropies for interacting fermionic systems. An explicit determination of the entanglement (modular) Hamiltonian has proven to be a considerably more difficult problem, and only a few results are available. We introduce a technique to directly determine the entanglement Hamiltonian of interacting fermionic models by means of auxiliary field quantum Monte Carlo simulations. We implement our method on the one-dimensional Hubbard chain partitioned into two segments and on the Hubbard two-leg ladder partitioned into two chains. In both cases, we study the evolution of the entanglement Hamiltonian as a function of the physical temperature.