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

Entanglement Hamiltonian of Interacting Fermionic Models

Francesco Parisen Toldin, Fakher F. Assaad

Published Nov 14, 2018• 68 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevlett.121.200602

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

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

Hybrid Monte Carlo approach to the entanglement entropy of interacting fermions

Joaquín E. Drut, William J. Porter

The Monte Carlo calculation of R\'enyi entanglement entropies ${S}_{n}^{}$ of interacting fermions suffers from a well-known signal-to-noise problem, even for a large number of situations in which the infamous sign problem is absent. A few methods have been proposed to overcome this issue, such as ensemble switching and the use of auxiliary partition-function ratios. Here, we present an approach that builds on the recently proposed free-fermion decomposition method; it incorporates entanglement in the probability measure in a natural way; it takes advantage of the hybrid Monte Carlo algorithm (an essential tool in lattice quantum chromodynamics and other gauge theories with dynamical fermions); and it does not suffer from noise problems. This method displays no sign problem for the same cases as other approaches and is therefore useful for a wide variety of systems. As a proof of principle, we calculate ${S}_{2}^{}$ for the one-dimensional, half-filled Hubbard model and compare with results from exact diagonalization and the free-fermion decomposition method.

arXiv:1506.06654•Physical Review B•Sep 14, 2015•31 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

High-efficiency quantum Monte Carlo algorithm for extracting entanglement entropy in interacting fermion systems

Weilun Jiang, Gaopei Pan, Wang, Zhe +3 more

The entanglement entropy probing novel phases and phase transitions numerically via quantum Monte Carlo has made great achievements in large-scale interacting spin/boson systems. In contrast, the numerical exploration in interacting fermion systems is rare, even though fermion systems attract more attentions in condensed matter. The fundamental restrictions is that the computational cost of fermion quantum Monte Carlo ($\sim βN^3$) is much higher than that of spin/boson ($\sim βN$). Here, $N$ is the total number of sites and $β$ is the inverse temperature or projection length. To tackle this problem, we propose a fermionic quantum Monte Carlo algorithm based on the incremental technique along physical parameters, which greatly improves the efficiency of extracting entanglement entropy. We benchmark the developed algorithm by calculating the scaling behavior of the entanglement entropy in a two-dimensional square lattice Hubbard model. The obtained phase diagram including Fermi surface and Goldstone modes validates the correctness of the algorithm. Remarkably, our method shows the high-efficiency with respect to the existing algorithms, while keeping the high computation precision. We proceed to apply this algorithm to explore the scaling behavior of the entanglement entropy and particularly its derivative at Gross-Neveu criticality. Our results elucidate that such critical behavior can be quantified by the correlation length exponent.

arXiv:2409.20009•Sep 30, 2024•0 citations

Entanglement spectra of interacting fermions in quantum Monte Carlo simulations

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

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.

arXiv:1311.5851•Physical Review B•Mar 24, 2014•80 citations
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.
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.
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.
The Monte Carlo calculation of R\'enyi entanglement entropies Sn{S}_{n}^{}Sn​ of interacting fermions suffers from a well-known signal-to-noise problem, even for a large number of situations in which the infamous sign problem is absent. A few methods have been proposed to overcome this issue, such as ensemble switching and the use of auxiliary partition-function ratios. Here, we present an approach that builds on the recently proposed free-fermion decomposition method; it incorporates entanglement in the probability measure in a natural way; it takes advantage of the hybrid Monte Carlo algorithm (an essential tool in lattice quantum chromodynamics and other gauge theories with dynamical fermions); and it does not suffer from noise problems. This method displays no sign problem for the same cases as other approaches and is therefore useful for a wide variety of systems. As a proof of principle, we calculate S2{S}_{2}^{}S2​ for the one-dimensional, half-filled Hubbard model and compare with results from exact diagonalization and the free-fermion decomposition method.
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.
The entanglement entropy probing novel phases and phase transitions numerically via quantum Monte Carlo has made great achievements in large-scale interacting spin/boson systems. In contrast, the numerical exploration in interacting fermion systems is rare, even though fermion systems attract more attentions in condensed matter. The fundamental restrictions is that the computational cost of fermion quantum Monte Carlo (∼βN3\sim βN^3∼βN3) is much higher than that of spin/boson (∼βN\sim βN∼βN). Here, NNN is the total number of sites and βββ is the inverse temperature or projection length. To tackle this problem, we propose a fermionic quantum Monte Carlo algorithm based on the incremental technique along physical parameters, which greatly improves the efficiency of extracting entanglement entropy. We benchmark the developed algorithm by calculating the scaling behavior of the entanglement entropy in a two-dimensional square lattice Hubbard model. The obtained phase diagram including Fermi surface and Goldstone modes validates the correctness of the algorithm. Remarkably, our method shows the high-efficiency with respect to the existing algorithms, while keeping the high computation precision. We proceed to apply this algorithm to explore the scaling behavior of the entanglement entropy and particularly its derivative at Gross-Neveu criticality. Our results elucidate that such critical behavior can be quantified by the correlation length exponent.