Entanglement in a pure state of a many-body system can be characterized by the Rényi entropies S^{(\alpha)}=\ln\textrm{tr}(\rho^\alpha)/(1-\alpha) S ( α ) = ln tr ( ρ α ) / ( 1 − α ) of the reduced density matrix \rho ρ of a subsystem. These entropies are, however, difficult to access experimentally and can typically be determined for small systems only. Here we show that for free fermionic systems in a Gaussian state and with particle number conservation, S^{(2)} S ( 2 ) can be tightly bound—from above and below—by the much easier accessible Rényi number entropy S^{(2)}_N=-\ln \sum_n p^2(n) S N ( 2 ) = − ln ∑ n p 2 ( n ) which is a function of the probability distribution p(n) p ( n ) of the total particle number in the considered subsystem only. A dynamical growth in entanglement, in particular, is therefore always accompanied by a growth—albeit logarithmically slower—of the number entropy. We illustrate this relation by presenting numerical results for quenches in non-interacting one-dimensional lattice models including disorder-free, Anderson-localized, and critical systems with off-diagonal (bond) disorder.
Entanglement in a pure state of a many-body system can be characterized by the Rényi entropies S^{(\alpha)}=\ln\textrm{tr}(\rho^\alpha)/(1-\alpha) S ( α ) = ln tr ( ρ α ) / ( 1 − α ) of the reduced density matrix \rho ρ of a subsystem. These entropies are, however, difficult to access experimentally and can typically be determined for small systems only. Here we show that for free fermionic systems in a Gaussian state and with particle number conservation, S^{(2)} S ( 2 ) can be tightly bound—from above and below—by the much easier accessible Rényi number entropy S^{(2)}_N=-\ln \sum_n p^2(n) S N ( 2 ) = − ln ∑ n p 2 ( n ) which is a function of the probability distribution p(n) p ( n ) of the total particle number in the considered subsystem only. A dynamical growth in entanglement, in particular, is therefore always accompanied by a growth—albeit logarithmically slower—of the number entropy. We illustrate this relation by presenting numerical results for quenches in non-interacting one-dimensional lattice models including disorder-free, Anderson-localized, and critical systems with off-diagonal (bond) disorder.
Abstract The symmetry-resolved Rényi entanglement entropy is the Rényi entanglement entropy of each symmetry sector of a density matrix ρ ρ . This experimentally relevant quantity is known to have rich theoretical connections to conformal field theory (CFT). For a family of critical free-fermion chains, we present a rigorous lattice-based derivation of its scaling properties using the theory of Toeplitz determinants. We consider a class of critical quantum chains with a microscopic U(1) symmetry; each chain has a low energy description given byNmassless Dirac fermions. For the density matrix, ρA ρA , of subsystems ofLneighbouring sites we calculate the leading terms in the largeLasymptotic expansion of the symmetry-resolved Rényi entanglement entropies. This follows from a largeLexpansion of the charged moments of ρA ρA ; we derive tr(eiαQAρAn)=aeiα⟨QA⟩(σL)−x(1+O(L−μ)) tr(eiαQAρAn)=aeiα⟨QA⟩(σL)-x(1+O(L-μ)) , wherea, xand μ μ are universal and σ σ depends only on theNFermi momenta. We show that the exponentxcorresponds to the expectation from CFT analysis. The error term O(L−μ) O(L-μ) is consistent with but weaker than the field theory prediction O(L−2μ) O(L-2μ) . However, using further results and conjectures for the relevant Toeplitz determinant, we find excellent agreement with the expansion over CFT operators.
Fermionic Gaussian states are a fundamental tool in many-body physics, faithfully representing non-interacting quantum systems and allowing for efficient numerical simulations. Given a many-body wave function, it is therefore interesting to ask how much it differs from that of a Gaussian state, as quantified by the notion of non-Gaussianity. In this work, we relate measures of non-Gaussianity with the Shannon entropy of the particle-number distribution, coinciding with the particle-number asymmetry for pure states. We derive a lower bound on the relative entropy of non-Gaussianity in terms of the exponential of the Shannon entropy, and study numerically its tightness for large system sizes. Our bound is non-trivial for large values of the asymmetry and relies on the concentration of the particle-number distribution of (mixed) fermionic Gaussian states. Since the Shannon entropy of the particle-number distribution is often efficient to compute or experimentally measure, our results can be viewed as a practical way to lower bound non-Gaussianity, highlighting a non-trivial interplay with particle-number asymmetry.
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 prove rigorous bounds on the growth of $\ensuremath{\alpha}$-R\'enyi entropies ${S}_{\ensuremath{\alpha}}(t)$ (the von Neumann entropy being the special case $\ensuremath{\alpha}=1$) associated with any subsystem A of a general lattice quantum many-body system with finite on-site Hilbert space dimension. For completely nonlocal Hamiltonians, we show that the instantaneous growth rates $|{S}_{\ensuremath{\alpha}}^{\ensuremath{'}}(t)|$ (with $\ensuremath{\alpha}\ensuremath{\ne}1$) can be exponentially larger than $|{S}_{1}^{\ensuremath{'}}(t)|$ as a function of the subsystem size ∣A∣. For D-dimensional systems with geometric locality, we prove bounds on $|{S}_{\ensuremath{\alpha}}^{\ensuremath{'}}(t)|$ that depend on the decay rate of interactions with distance. When $\ensuremath{\alpha}=1$, the bound is ∣A∣ independent for all power-law decaying interactions $V(r)\ensuremath{\sim}{r}^{\ensuremath{-}w}$ with w>2D+1. However, for $\ensuremath{\alpha}>1$, the bound is ∣A∣ independent only when the interactions are finite-range or decay faster than $V(r)\ensuremath{\sim}{e}^{\ensuremath{-}c\phantom{\rule{0.16em}{0ex}}{r}^{D}}$ for some c depending on the local Hilbert space dimension. Using similar arguments, we also prove bounds on k-local systems with or without geometric locality. A central theme of this work is that the value of $\ensuremath{\alpha}$ strongly influences the interplay between locality and instantaneous entanglement growth. In other words, the von Neumann entropy and the $\ensuremath{\alpha}$-R\'enyi entropies cannot be regarded as proxies for each other in studies of entanglement dynamics. We compare these bounds with analytic and numerical results on Hamiltonians with varying degrees of locality and find concrete examples that almost saturate the bound for nonlocal dynamics.
We have recently shown that the logarithmic growth of the entanglement entropy following a quantum quench in a many-body localized phase is accompanied by a slow growth of the number entropy ${S}_{N}\ensuremath{\sim}ln\phantom{\rule{0.16em}{0ex}}ln\phantom{\rule{0.16em}{0ex}}t$. Here we provide an in-depth numerical study of SN(t) for the disordered Heisenberg chain and show that this behavior is not transient and persists even for very strong disorder. Calculating the truncated R\'enyi number entropy ${S}_{N}^{(\ensuremath{\alpha})}(t)={(1\ensuremath{-}\ensuremath{\alpha})}^{\ensuremath{-}1}ln{\ensuremath{\sum}}_{n}{p}^{\ensuremath{\alpha}}(n)$ for $\ensuremath{\alpha}\ensuremath{\ll}1$ and p(n)>pc---which is sensitive to large number fluctuations occurring with low probability---we demonstrate that the particle number distribution p(n) in one half of the system has a continuously growing tail. This indicates a slow but steady increase in the number of particles crossing between the partitions in the interacting case and is in sharp contrast to Anderson localization for which we show that ${S}_{N}^{(\ensuremath{\alpha}\ensuremath{\rightarrow}0)}(t)$ saturates for any cutoff pc>0. We show, furthermore, that the growth of SN is not the consequence of rare states or rare regions but rather represents typical behavior. These findings indicate that the interacting system is never fully localized even for very strong but finite disorder.
We calculate the typical bipartite entanglement entropy ${\ensuremath{\langle}{S}_{A}\ensuremath{\rangle}}_{N}$ in systems containing indistinguishable particles of any kind as a function of the total particle number N, the volume V, and the subsystem fraction f=VA/V, where VA is the volume of the subsystem. We expand our result as a power series ${\ensuremath{\langle}{S}_{A}\ensuremath{\rangle}}_{N}=afV+b\sqrt{V}+c+o(1)$, and find that c is universal (i.e., independent of the system type), while a and b can be obtained from a generating function characterizing the local Hilbert space dimension. We illustrate the generality of our findings by studying a wide range of different systems, e.g., bosons, fermions, spins, and mixtures thereof. We provide evidence that our analytical results describe the entanglement entropy of highly excited eigenstates of quantum-chaotic spin and boson systems, which is distinct from that of integrable counterparts.