Abstract We prove that the Hilbert–Schmidt norm of k -particle reduced density matrices of N -body fermionic states is bounded from above by $$C_kN^{k/2}$$ C k N k / 2 - matching the scaling behaviour of Slater determinant states. This generalises a recent result of Christiansen [3] on 2-particle reduced density matrices to higher order density matrices. Moreover, our estimate directly yields a lower bound on the von Neumann entropy and the 2-Rényi entropy of reduced density matrices, thereby providing further insight into conjectures of Carlen–Lieb–Reuvers [2, 8].
Abstract We prove that the Hilbert–Schmidt norm of k -particle reduced density matrices of N -body fermionic states is bounded from above by CkNk/2 C k N k / 2 - matching the scaling behaviour of Slater determinant states. This generalises a recent result of Christiansen [3] on 2-particle reduced density matrices to higher order density matrices. Moreover, our estimate directly yields a lower bound on the von Neumann entropy and the 2-Rényi entropy of reduced density matrices, thereby providing further insight into conjectures of Carlen–Lieb–Reuvers [2, 8].
Abstract We prove that the 2-body operator γ2Ψ γ 2 Ψ of a fermionic N-particle state Ψ Ψ obeys ∥γ2Ψ∥HS≤5N ‖ γ 2 Ψ ‖ HS ≤ 5 N , which complements the bound of Yang (Rev Mod Phys 34:694, 1962) that ∥γ2Ψ∥op≤N ‖ γ 2 Ψ ‖ op ≤ N . This estimate furthermore resolves a conjecture of Carlen–Lieb–Reuvers (Commun Math Phys 344:655–671, 2016) concerning the entropy of the normalized 2-body operator. We also prove that the Hilbert–Schmidt norm of the truncated 2-body operator γ2Ψ,T γ 2 Ψ , T obeys the inequality ∥γ2Ψ,T∥HS≤5Ntr(γ1Ψ(1−γ1Ψ)) ‖ γ 2 Ψ , T ‖ HS ≤ 5 N tr ( γ 1 Ψ ( 1 - γ 1 Ψ ) ) .
We discuss a general bipartitelike representation and Schmidt decomposition of an arbitrary pure state of N indistinguishable fermions, based on states of M<N and $(N\ensuremath{-}M)$ fermions. It is directly connected with the reduced M- and $(N\ensuremath{-}M)$-body density matrices (DMs), which have the same spectrum in such states. The concept of M-body entanglement emerges naturally in this scenario, generalizing that of one-body entanglement. Rigorous majorization relations satisfied by the normalized M-body DM are then derived, which imply that the associated entropy will not increase, on average, under a class of operations which have these DMs as postmeasurement states. Moreover, such entropy is an upper bound to the bipartite entanglement entropy generated by a class of operations which map the original state to a bipartite state of M and $N\ensuremath{-}M$ effectively distinguishable fermions. Analytic evaluation of the spectrum of M-body DMs in some strongly correlated fermionic states is also provided.
Fermionic Hamiltonians play a critical role in quantum chemistry, one of the most promising use cases for near-term quantum computers. However, since encoding nonlocal fermionic statistics using conventional qubits results in significant computational overhead, fermionic quantum hardware, such as fermion atom arrays, were proposed as a more efficient platform. In this context, we here study the many-body entanglement structure of fermionic N-particle states by concentrating on M-body reduced density matrices (DMs) across various bipartitions in Fock space. The von Neumann entropy of the reduced DM is a basis independent entanglement measure which generalizes the traditional quantum chemistry concept of the one-particle DM entanglement, which characterizes how a single fermion is entangled with the rest. We carefully examine the upper bounds on the M-body entanglement, which are analogous to the volume law of conventional entanglement measures. To this end, we establish a connection between M-body reduced DM and the mathematical structure of hypergraphs. Specifically, we show that a special class of hypergraphs, known as t-designs, corresponds to maximally entangled fermionic states. We also explore fermionic many-body entanglement in random states. We semianalytically demonstrate that the distribution of reduced DMs associated with random fermionic states corresponds to the trace-fixed Wishart-Laguerre random matrix ensemble. In the limit of large single-particle dimension D and a non-zero filling fraction, random states asymptotically become absolutely maximally entangled. We also analyze many-body entanglement in interacting fermionic systems by using the M-body DM. We demonstrate that if a particle number conserving fermionic Hamiltonian contains only up to M-body interaction terms, then its N-particle ground state cannot be maximally M-body entangled. As a key step in the proof, we show that the energy expectation value of a maximally M-body mixed state is equal to the spectral mean of the Hamiltonian on the corresponding N-particle subspace. We further demonstrate that the many-body entanglement structure of a ground state can place quantitative constraint on the interaction strength of its parent Hamiltonian. We illustrate the theorem and its implications in Hubbard and extended SYK models. Going beyond ground states, we analyze entanglement generation under unitary dynamics from Slater-determinant initial states in these models. We determine early-time growth and estimate entanglement saturation times. Finally, we derive explicit symmetry-refined saturation upper bounds for M-body entanglement in the presence of an Abelian symmetry.
Abstract We consider entanglement of first-quantized identical particles by adopting an algebraic approach. In particular, we investigate fermions whose wave functions are given by the Slater determinants, as for singlet sectors of one-matrix models. We show that the upper bounds of the general Rényi entropies are N log 2 for N particles or an N × N matrix. We compute the target space entanglement entropy and the mutual information in a free one-matrix model. We confirm the area law: the single-interval entropy for the ground state scales as 31 1 3 log N in the large N model. We obtain an analytical O(N0) O N 0 expression of the mutual information for two intervals in the large N expansion.
The following thesis focuses on the scaling of entanglement entropy in lower dimensions and is divided into three main parts. Chapter 2 studies the thermal reduced density matrices in fermion and spin systems on ladders. Chapter 3 studies the many-body localization phase transition in a Rokhsar-Kivelson type wave function. Chapter 4 studies the subleading correction term of entanglement entropy in 2+1 dimensional scale invariant systems. Chapter 5 studies the bulk-boundary correspondence in 3 + 1 dimensional topological phases and its entanglement entropy. In chapter 2, we investigate the reduced density matrices for a model of free fermions on a two-leg ladder (gapped by the inter-chain tunneling operator) and in 1/2 spin systems on a ladder with a gapped ground state using exact solutions and several controlled approximations. We calculate the reduced density matrix and the entanglement entropy for a leg of the ladder (i.e. cut made between the chains). In the fermionic system we find the exact form of the reduced density matrix for one of the chains and determine the entanglement spectrum explicitly. Here we find that in the weak tunneling limit of the ladder the entanglement entropy of one chain of the gapped ladder has a simple and universal form dictated by conformal invariance. In the case of the spin system, we consider the strong coupling limit by using perturbation theory and get the reduced density matrix by the Schmidt decomposition. The entanglement entropies of a general gapped system of two coupled conformal field theories (in 1+1 dimensions) is discussed using the replica trick and scaling arguments. We show that 1) for a system with a bulk gap, the reduced density matrix has the form of a thermal density matrix, 2) the long-wavelength modes of one subsystem (a chain) of a gapped coupled system are always thermal, 3) the von Neumann entropy equals to the thermodynamic entropy of one chain, and 4) the bulk gap plays the role of effective temperature. In chapter 3, we construct a family of many-body wave functions to study the many-body localization phase transition. The wave functions have a Rokhsar-Kivelson form, in which the weight for the configurations are chosen from the Gibbs weights of a classical spin glass model, known as the Random Energy Model, multiplied by a random sign structure to represent a highly excited state. These wave functions show a phase transition into an MBL phase. In addition, we see three regimes of entanglement scaling with subsystem size: scaling with entanglement corresponding to an infinite temperature thermal phase, constant scaling, and a sub-extensive scaling between these limits. Near the phase transition point, the fluctuations of the Renyi entropies are non-Gaussian. We find that Renyi entropies with different Renyi index transition into the MBL phase at different points and have different scaling behavior, suggesting a multifractal behavior. In chapter 4, we study the universal scaling behavior of the entanglement entropy of critical theories in 2+1 dimensions. We specially consider two fermionic scale-invariant models, free massless Dirac fermions and a model of fermions with quadratic band touching, and numerically study the two-cylinder entanglement entropy of the models on the torus. We find that in both cases the entanglement entropy satisfies the area law and has the subleading term which is a scaling function of the aspect ratios of the cylindrical regions. We test the scaling of entanglement in both the free fermion models using three possible scaling functions for the subleading term derived from a) the quasi-one-dimensional conformal field theory, b) the bosonic quantum Lifshitz model, and c) the holographic AdS/CFT correspondence. For the later case we construct an analytic scaling function using holography, appropriate for critical theories with a gravitational dual description. We find that the subleading term in the fermionic models is well described, for a range of aspect ratios, by the scaling form derived from the quantum Lifshitz model as well as that derived using the AdS/CFT correspondence (in this case only for the Dirac model). For the case where the fermionic models are placed on a square torus we find the fit to the different scaling forms is in agreement to surprisingly high precision. In chapter 5, we discuss (2+1)-dimensional gapless surface theories of bulk (3+1)-dimensional topological phases, such as the BF theory at level K, and its generalization. In particular, we put these theories on a flat (2+1) dimensional torus T3 parameterized by its modular parameters, and compute the partition functions obeying various twisted boundary conditions. We show the partition functions are transformed into each other under SL(3, Z) modular transformations, and furthermore establish the bulk-boundary correspondence in (3+1) dimensions by matching the modular S and T matrices computed from the boundary field theories with those computed in the bulk. We propose the three-loop braiding statistics can be studied by constructing the modular S and T matrices from an appropriate boundary field theory. We also study the EE for 3 + 1- dimensional topological phase with or without three-loop braiding phase.