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:2209.07103 · SciPost Physics

Many versus one: The disorder operator and entanglement entropy in fermionic quantum matter

Weilun Jiang, Bin-Bin Chen, Zi Hong Liu, Junchen Rong, Fakher F. Assaad, Meng Cheng, Kai Sun, Zi Yang Meng

Published Sep 7, 2023• 26 citations
arXivHTMLPDF

Abstract

Motivated by recent development of the concept of the disorder operator and its relation with entanglement entropy in bosonic systems, here we show the disorder operator successfully probes many aspects of quantum entanglement in fermionic many-body systems. From both analytical and numerical computations in free and interacting fermion systems in 1D and 2D, we find the disorder operator and the entanglement entropy exhibit similar universal scaling behavior, as a function of the boundary length of the subsystem, but with subtle yet important differences. In 1D they both follow the \log{L} log L scaling behavior with the coefficient determined by the Luttinger parameter for disorder operator, and the conformal central charge for entanglement entropy. In 2D they both show the universal L\log L L log L scaling behavior in free and interacting Fermi liquid states, with the coefficients depending on the geometry of the Fermi surfaces. However at a 2D quantum critical point with non-Fermi-liquid state, extra symmetry information is needed in the design of the disorder operator, so as to reveal the critical fluctuations as does the entanglement entropy. Our results demonstrate the fermion disorder operator can be used to probe quantum many-body entanglement related to global symmetry, and provide new tools to explore the still largely unknown territory of highly entangled fermion quantum matter in 2 or higher dimensions.

DOI
10.21468/scipostphys.15.3.082

Related papers

Disorder operators in two-dimensional Fermi and non-Fermi liquids through multidimensional bosonization

Kang-Le Cai, Meng Cheng

Disorder operators are a type of non-local observables for quantum many-body systems, measuring the fluctuations of symmetry charges inside a region. It has been shown that disorder operators can reveal global aspects of many-body states that are otherwise difficult to access through local measurements. We study the disorder operator for U(1) (charge or spin) symmetry in two-dimensional Fermi and non-Fermi liquid states, using the multidimensional bosonization formalism. For a region $A$, the logarithm of the charge disorder parameter in a Fermi liquid with isotropic interactions scales asymptotically as $l_A\ln l_A$, with $l_A$ being the linear size of the region $A$. We calculate the proportionality coefficient in terms of Landau parameters of the Fermi liquid theory. We then study models of the Fermi surface coupled to gapless bosonic fields realizing non-Fermi liquid states. In a simple spinless model, where the fermion density is coupled to a critical scalar, we find that at the quantum critical point the scaling behavior of the charge disorder operators is drastically modified to $l_A \ln^2 l_A$. We also consider the composite Fermi liquid state and argue that the charge disorder operator scales as $l_A$.

arXiv:2404.04334•Physical review. B./Physical review. B•Sep 22, 2025•3 citations

Entanglement entropy of fermions in a strange metal

Santanu Singh, Surajit Bera, Chenyuan Li +2 more

The subsystem-size dependence of ground-state entanglement entropy and its crossover to thermal entropy as a function of temperature are well understood for one-dimensional (1D) gapless systems described by conformal field theory (CFT), and for free fermions with a Fermi surface in any dimension. However, little is known about the entanglement entropy for gapless fermionic systems without quasi-particles, such as a strange metal. Here we study the entanglement entropy of fermions in a solvable large-$N$ 1D lattice model akin to the Yukawa-Sachdev-Ye-Kitaev (Yukawa SYK) model. In this model, two Fermi points are coupled to scalar bosons via spatially random Yukawa interactions, providing a solvable model of a strange metal when the bosons become critical at a quantum critical point. We exactly compute the second Rényi entropy of fermions in a spatial subregion in this model. Our results unravel crucial role of intra-subregion entanglement between fermionic and bosonic degrees of freedom along with the inter-subregion entanglement in understanding the ground states of such strongly coupled fermion-boson systems. We show that the crossover from thermal entropy to entanglement entropy, is captured by a single scaling ansatz, that collapses the second Rényi entropy of fermions for different subregion sizes and temperatures into a single universal curve. We further show that the universal scaling curve for the critical strange metal is well described by the standard CFT formula, albeit with an effective central charge substantially larger than the non-interacting value.

cond-mat.str-elCross-listedarXiv:2608.04098•Aug 6, 2026•0 citations

Fermion Disorder Operator at Gross-Neveu and Deconfined Quantum Criticalities

Zi Hong Liu, Weilun Jiang, Bin-Bin Chen +5 more

The fermion disorder operator has been shown to reveal the entanglement information in 1D Luttinger liquids and 2D free and interacting Fermi and non-Fermi liquids emerging at quantum critical points (QCPs) [W. Jiang et al., arXiv:2209.07103]. Here we study, by means of large-scale quantum Monte Carlo simulation, the scaling behavior of the disorder operator in correlated Dirac systems. We first demonstrate the logarithmic scaling behavior of the disorder operator at the Gross-Neveu (GN) chiral Ising and Heisenberg QCPs, where consistent conformal field theory (CFT) content of the GN-QCP in its coefficient is found. Then we study a 2D monopole-free deconfined quantum critical point (DQCP) realized between a quantum-spin Hall insulator and a superconductor. Our data point to negative values of the logarithmic coefficients such that the DQCP does not correspond to a unitary CFT. Density matrix renormalization group calculations of the disorder operator on a 1D DQCP model also detect emergent continuous symmetries.

arXiv:2212.11821•Physical Review Letters•Jun 28, 2023•37 citations

Universal entanglement signatures of quantum liquids as a guide to fermionic criticality

Siddhartha Patra, Anirban Mukherjee, Siddhartha Lal

Abstract An outstanding challenge involves understanding the many-particle entanglement of liquid states of quantum matter that arise in systems of interacting electrons. The Fermi liquid (FL) shows a violation of the area-law in real-space entanglement entropy of a subsystem, believed to be a signature of the ground state of a gapless quantum critical system of interacting fermions. Here, we apply a T = 0 renormalization group approach to the FL, unveiling the long-wavelength quantum fluctuations from which long-range entanglement arises. A similar analysis of non-Fermi liquids such as the 2D marginal Fermi liquid (MFL) and the 1D Tomonaga–Luttinger liquid reveals a universal logarithmic violation of the area-law in gapless electronic liquids, with a proportionality constant that depends on the nature of the underlying Fermi surface. We extend this analysis to classify the gapped quantum liquids emergent from the destabilisation of the Fermi surface by renormalisation group relevant quantum fluctuations arising from backscattering processes.

arXiv:2205.11123•New Journal of Physics•May 25, 2023•2 citations

Entanglement entropies and fermion signs of critical metals

Nikolaos K. Kaplis, Frank Krüger, Jan Zaanen

The fermion sign problem is often viewed as a sheer inconvenience that plagues numerical studies of strongly interacting electron systems. Only recently has it been suggested that fermion signs are fundamental for the universal behavior of critical metallic systems and crucially enhance their degree of quantum entanglement. In this work we explore potential connections between emergent scale invariance of fermion sign structures and scaling properties of bipartite entanglement entropies. Our analysis is based on a wave-function Ansatz that incorporates collective, long-range backflow correlations into fermionic Slater determinants. Such wave functions mimic the collapse of a Fermi liquid at a quantum critical point. Their nodal surfaces, a representation of the fermion sign structure in many-particle configurations space, show fractal behavior up to a length scale $\ensuremath{\xi}$ that diverges at a critical backflow strength. We show that the Hausdorff dimension of the fractal nodal surface depends on $\ensuremath{\xi}$, the number of fermions and the exponent of the backflow. For the same wave functions we numerically calculate the second R\'enyi entanglement entropy ${S}_{2}$. Our results show a crossover from volume scaling, ${S}_{2}\ensuremath{\sim}{\ensuremath{\ell}}^{\ensuremath{\theta}}$ $(\ensuremath{\theta}=2$ in $d=2$ dimensions), to the characteristic Fermi-liquid behavior ${S}_{2}\ensuremath{\sim}\ensuremath{\ell}ln\ensuremath{\ell}$ on scales larger than $\ensuremath{\xi}$. We find that volume scaling of the entanglement entropy is a robust feature of critical backflow fermions, independent of the backflow exponent and hence the fractal dimension of the scale invariant sign structure.

arXiv:1605.02477•Physical review. B./Physical review. B•Apr 3, 2017•6 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
Motivated by recent development of the concept of the disorder operator and its relation with entanglement entropy in bosonic systems, here we show the disorder operator successfully probes many aspects of quantum entanglement in fermionic many-body systems. From both analytical and numerical computations in free and interacting fermion systems in 1D and 2D, we find the disorder operator and the entanglement entropy exhibit similar universal scaling behavior, as a function of the boundary length of the subsystem, but with subtle yet important differences. In 1D they both follow the \log{L} log L scaling behavior with the coefficient determined by the Luttinger parameter for disorder operator, and the conformal central charge for entanglement entropy. In 2D they both show the universal L\log L L log L scaling behavior in free and interacting Fermi liquid states, with the coefficients depending on the geometry of the Fermi surfaces. However at a 2D quantum critical point with non-Fermi-liquid state, extra symmetry information is needed in the design of the disorder operator, so as to reveal the critical fluctuations as does the entanglement entropy. Our results demonstrate the fermion disorder operator can be used to probe quantum many-body entanglement related to global symmetry, and provide new tools to explore the still largely unknown territory of highly entangled fermion quantum matter in 2 or higher dimensions.
Disorder operators are a type of non-local observables for quantum many-body systems, measuring the fluctuations of symmetry charges inside a region. It has been shown that disorder operators can reveal global aspects of many-body states that are otherwise difficult to access through local measurements. We study the disorder operator for U(1) (charge or spin) symmetry in two-dimensional Fermi and non-Fermi liquid states, using the multidimensional bosonization formalism. For a region AAA, the logarithm of the charge disorder parameter in a Fermi liquid with isotropic interactions scales asymptotically as lAln⁡lAl_A\ln l_AlA​lnlA​, with lAl_AlA​ being the linear size of the region AAA. We calculate the proportionality coefficient in terms of Landau parameters of the Fermi liquid theory. We then study models of the Fermi surface coupled to gapless bosonic fields realizing non-Fermi liquid states. In a simple spinless model, where the fermion density is coupled to a critical scalar, we find that at the quantum critical point the scaling behavior of the charge disorder operators is drastically modified to lAln⁡2lAl_A \ln^2 l_AlA​ln2lA​. We also consider the composite Fermi liquid state and argue that the charge disorder operator scales as lAl_AlA​.
The subsystem-size dependence of ground-state entanglement entropy and its crossover to thermal entropy as a function of temperature are well understood for one-dimensional (1D) gapless systems described by conformal field theory (CFT), and for free fermions with a Fermi surface in any dimension. However, little is known about the entanglement entropy for gapless fermionic systems without quasi-particles, such as a strange metal. Here we study the entanglement entropy of fermions in a solvable large-NNN 1D lattice model akin to the Yukawa-Sachdev-Ye-Kitaev (Yukawa SYK) model. In this model, two Fermi points are coupled to scalar bosons via spatially random Yukawa interactions, providing a solvable model of a strange metal when the bosons become critical at a quantum critical point. We exactly compute the second Rényi entropy of fermions in a spatial subregion in this model. Our results unravel crucial role of intra-subregion entanglement between fermionic and bosonic degrees of freedom along with the inter-subregion entanglement in understanding the ground states of such strongly coupled fermion-boson systems. We show that the crossover from thermal entropy to entanglement entropy, is captured by a single scaling ansatz, that collapses the second Rényi entropy of fermions for different subregion sizes and temperatures into a single universal curve. We further show that the universal scaling curve for the critical strange metal is well described by the standard CFT formula, albeit with an effective central charge substantially larger than the non-interacting value.
The fermion disorder operator has been shown to reveal the entanglement information in 1D Luttinger liquids and 2D free and interacting Fermi and non-Fermi liquids emerging at quantum critical points (QCPs) [W. Jiang et al., arXiv:2209.07103]. Here we study, by means of large-scale quantum Monte Carlo simulation, the scaling behavior of the disorder operator in correlated Dirac systems. We first demonstrate the logarithmic scaling behavior of the disorder operator at the Gross-Neveu (GN) chiral Ising and Heisenberg QCPs, where consistent conformal field theory (CFT) content of the GN-QCP in its coefficient is found. Then we study a 2D monopole-free deconfined quantum critical point (DQCP) realized between a quantum-spin Hall insulator and a superconductor. Our data point to negative values of the logarithmic coefficients such that the DQCP does not correspond to a unitary CFT. Density matrix renormalization group calculations of the disorder operator on a 1D DQCP model also detect emergent continuous symmetries.
Abstract An outstanding challenge involves understanding the many-particle entanglement of liquid states of quantum matter that arise in systems of interacting electrons. The Fermi liquid (FL) shows a violation of the area-law in real-space entanglement entropy of a subsystem, believed to be a signature of the ground state of a gapless quantum critical system of interacting fermions. Here, we apply a T = 0 renormalization group approach to the FL, unveiling the long-wavelength quantum fluctuations from which long-range entanglement arises. A similar analysis of non-Fermi liquids such as the 2D marginal Fermi liquid (MFL) and the 1D Tomonaga–Luttinger liquid reveals a universal logarithmic violation of the area-law in gapless electronic liquids, with a proportionality constant that depends on the nature of the underlying Fermi surface. We extend this analysis to classify the gapped quantum liquids emergent from the destabilisation of the Fermi surface by renormalisation group relevant quantum fluctuations arising from backscattering processes.
The fermion sign problem is often viewed as a sheer inconvenience that plagues numerical studies of strongly interacting electron systems. Only recently has it been suggested that fermion signs are fundamental for the universal behavior of critical metallic systems and crucially enhance their degree of quantum entanglement. In this work we explore potential connections between emergent scale invariance of fermion sign structures and scaling properties of bipartite entanglement entropies. Our analysis is based on a wave-function Ansatz that incorporates collective, long-range backflow correlations into fermionic Slater determinants. Such wave functions mimic the collapse of a Fermi liquid at a quantum critical point. Their nodal surfaces, a representation of the fermion sign structure in many-particle configurations space, show fractal behavior up to a length scale $\ensuremath{\xi}$ that diverges at a critical backflow strength. We show that the Hausdorff dimension of the fractal nodal surface depends on $\ensuremath{\xi}$, the number of fermions and the exponent of the backflow. For the same wave functions we numerically calculate the second R\'enyi entanglement entropy S2{S}_{2}S2​. Our results show a crossover from volume scaling, ${S}_{2}\ensuremath{\sim}{\ensuremath{\ell}}^{\ensuremath{\theta}}$ $(\ensuremath{\theta}=2$ in d=2d=2d=2 dimensions), to the characteristic Fermi-liquid behavior ${S}_{2}\ensuremath{\sim}\ensuremath{\ell}ln\ensuremath{\ell}$ on scales larger than $\ensuremath{\xi}$. We find that volume scaling of the entanglement entropy is a robust feature of critical backflow fermions, independent of the backflow exponent and hence the fractal dimension of the scale invariant sign structure.
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.