Density of States at Fixed Entanglement

• Density of states at fixed entanglement is the measure of quantum states whose reduced density matrices maintain a prescribed entanglement, analogous to fixed-energy ensembles in statistical mechanics.
• Methodologies including random matrix theory and Coulomb gas models are employed to derive equilibrium distributions and identify phase transitions in the entanglement spectrum.
• Applications span quantum information, many-body physics, and holography, providing practical insights into entanglement fluctuations and the microstate counting of complex systems.

The density of states at fixed entanglement quantifies the number or distribution of quantum states whose reduced density matrices possess prescribed entanglement properties, often measured via the entanglement entropy or purity. This object serves as the entanglement-analog of the microcanonical density of states at fixed energy in statistical mechanics. It has significant roles in quantum information theory, random matrix theory, many-body physics, and holography, providing insight into both the typical and atypical structure of entanglement, phase transitions in entanglement spectrum statistics, and the microcanonical ensembles associated to fixed entanglement entropy.

1. Foundations: Entanglement Spectrum and Entropy

For a pure quantum state $|\Psi\rangle$ defined on a bipartite system $\mathcal{H}_A \otimes \mathcal{H}_{A^\mathrm{c}}$, the reduced density matrix on subsystem $A$ is

$$\rho_A = \mathrm{Tr}_{A^\mathrm{c}} |\Psi\rangle\langle\Psi|\,,$$

which can be diagonalized as

$$\rho_A = \sum_{i}\lambda_i |i\rangle_A\langle i|_A\,,$$

where the collection $\{\lambda_i\}$ is called the entanglement spectrum. The von Neumann entanglement entropy is

$S(\rho_A) = -\sum_i \lambda_i \ln \lambda_i\,,$

and alternative measures include purity $P = \sum_i \lambda_i^2$ or Rényi entropies. The distribution and typicality of the entanglement spectrum for reduced subregions—especially at fixed values of $S(\rho_A)$ or $P$—forms the subject of the density of states at fixed entanglement (Facchi et al., 2013, Cunden et al., 2013, Guo, 2020, Shekhar et al., 2 Feb 2024, Dong et al., 2018).

2. Microcanonical Ensembles and Density of States at Fixed Entanglement

The density of states at fixed entanglement refers to the number (or measure) of pure states whose reduced density matrices possess a prescribed value (or narrow range) of an entanglement functional, e.g.,

$$g(S_A)\,dS_A = \frac{\textrm{number of pure states with } S_A \in [S_A, S_A + dS_A]}{\textrm{total number of states}}\,,$$

with analogous definitions for purity, Rényi entropy, or other entanglement measures (Shekhar et al., 2 Feb 2024, Cunden et al., 2013).

Formally, the microcanonical ensemble at fixed entanglement entropy $S_A$ is the uniform ensemble over all states whose reduced density matrices have $S(\rho_A) = S_A$, and the corresponding density of states is the measure of this hypersurface.

Example: Purity Microcanonical Ensemble

For fixed purity $P$, the microcanonical density of states $\Omega(P)$ is given (in the large-$N$ limit for a subsystem of dimension $N$) by a large deviation form

$$\Omega(P) \simeq \exp\bigl[N^2\,S(P)\bigr],\qquad S(P) = \frac{1}{N P - 1} + \frac{1}{2}\ln\frac{1}{N P - 1} - 1\,.$$

This expression describes the exponential rarity (or typicality) of states at atypical (resp. typical) purity (Cunden et al., 2013).

3. Random Matrix and Coulomb Gas Approaches

Random matrix theory provides a powerful statistical framework to characterize the ensemble of reduced density matrices, particularly for Haar-random pure states. The eigenvalues $\{\lambda_i\}$ are distributed according to a joint PDF reflecting the induced Haar measure and constraints such as fixed trace, and possibly fixed entanglement or purity. This situation maps directly to a Coulomb gas model with logarithmic repulsion and external fields: $$p_N(\vec{\lambda}) \propto \prod_{i<j} (\lambda_i - \lambda_j)^2\,,$$ subject to $\sum_i \lambda_i = 1$ and $\lambda_i \geq 0$ (Facchi et al., 2013).

Imposing a fixed von Neumann entropy $S_\mathrm{vN} = \ln N - u$ introduces a further constraint, analogous to a microcanonical ensemble at that entropy. The most probable entanglement spectrum is then obtained by a variational principle, leading (in the large-$N$ limit) to an equilibrium density $\sigma(x)$ supported on $[a,b]$, given by saddle-point and Tricomi techniques (Facchi et al., 2013).

Phases of the Entanglement Spectrum

Analysis yields distinct entanglement-spectrum regimes:

• Gapped (high-entropy) phase: $\beta > 3/2$. The spectrum avoids zero, and $\sigma(x)$ is a deformed semicircle.
• Touching (intermediate-entropy) phase: $0 \leq \beta < 3/2$. The spectrum touches zero, generating a square-root singularity.
• Separable-tail phase: $\beta < 0$ ($u > 1/2$). The largest eigenvalue detaches, signaling a nearly separable state (Facchi et al., 2013).

4. Holographic States and Flat Entanglement Spectra at Fixed Area

In holographic theories, states prepared by Euclidean path integrals can be projected onto subspaces of fixed area $A$ of the HRT (Ryu-Takayanagi/Hubeny-Rangamani-Takayanagi) surface. The resulting reduced density matrix $\rho_{R,A}$, to leading order in $1/G$ (semiclassical gravity), becomes

$$\rho_{R,A} \propto \frac{1}{d}\,I\,,\qquad d = \exp\left(\frac{A}{4G}\right),$$

implying a flat entanglement spectrum and all Rényi entropies $S_n = A/(4G)$ independent of $n$ (Dong et al., 2018).

The corresponding density of states at fixed entropy $S_A = A/(4G)$ is

$g(S_A) = \exp(S_A)\,,$

matching the Bekenstein-Hawking microcanonical counting of geometrical edge modes. The flat spectrum arises from maximal mixing of center/edge degrees of freedom on the entanglement wedge boundary (Dong et al., 2018).

Restoring area fluctuations (i.e., not fixing $A$) and performing a saddle-point integral recovers the standard $n$-dependence of Rényi entropies, linking the $n$-dependence directly to the HRT-area fluctuations around the saddle point (Dong et al., 2018).

5. Exact Distributions and Fluctuations of Entanglement at Fixed Value

Shekhar & Shukla derived the full probability density $f_v(R_1;Y)$ for the von Neumann entropy $R_1 = -\sum_{n=1}^{N_A}\lambda_n \ln \lambda_n$ in non-ergodic (or ergodic) states, modeled by reduced density matrices from the generalized Wishart ensemble with unit trace. The time-like parameter $Y$ interpolates between ensembles with typical separable and highly entangled states (Shekhar et al., 2 Feb 2024).

As $Y \to \infty$, $f_v(R_1;Y)$ approaches a delta function at $R_1 \approx \ln N_A$, corresponding to the microcanonical density of states at maximum entanglement. For finite $Y$, the distribution broadens, and $f_v(R_1;Y)\,dR_1$ is the normalized number of states with entanglement entropy in $[R_1, R_1 + dR_1]$ (Shekhar et al., 2 Feb 2024).

6. Applications: Quantum Information, Many-Body Systems, and Holography

Density of states at fixed entanglement underpins classification of random pure states, large deviation studies of entanglement, and phase transitions in the entanglement spectrum, as well as microcanonical quantum information protocols. In holographic models, it formalizes the microstate counting implicit in the AdS/CFT correspondence and supports the structure of quantum error-correcting tensor network models, providing a continuum realization of maximally entangled EPR links (Dong et al., 2018).

In conformal field theory, the density of eigenstates at fixed entanglement is extracted by Laplace-inverting the replica-trick formula for $\operatorname{Tr} \rho_A^n$, leading to explicit expressions for the density (e.g., in terms of Bessel functions for single-interval states) (Guo, 2020). Microcanonical ensemble states constructed in this way have expectation values for non-extensive local probes that match those of the original reduced density matrix for small regions (Guo, 2020).

7. Open Problems and Physical Implications

Several unresolved questions remain. The identification of geometric duals for measures such as the Holevo information in holographic contexts remains open. The precise "critical region" demarcating distinguishable and indistinguishable microcanonical eigenstates near subsystem boundaries is yet to be sharply formulated (Guo, 2020). Connections between microcanonical states at fixed entanglement and operator-algebra quantum error correction structures continue to be elucidated (Dong et al., 2018). Exploring large deviation principles and full entanglement spectrum statistics in more general non-ergodic or out-of-equilibrium ensembles presents further fertile ground (Shekhar et al., 2 Feb 2024).

---

References:

• (Facchi et al., 2013) Facchi, G., Florio, G., Marzolino, U., Parisi, G., & Pascazio, S., "Entropy-Driven Phase Transitions of Entanglement."
• (Cunden et al., 2013) Facchi, G., Florio, G., Marzolino, U., Parisi, G., & Pascazio, S., "Typical Entanglement."
• (Dong et al., 2018) Dong, X., Harlow, D., & Marolf, D., "Flat entanglement spectra in fixed-area states of quantum gravity."
• (Guo, 2020) Guo, Q., "Entanglement spectrum of geometric states."
• (Shekhar et al., 2 Feb 2024) Shekhar, Y. & Shukla, P., "Distribution of the entanglement entropy of a non-ergodic quantum state."