Average entanglement entropy of a small subsystem in a constrained pure Gaussian state ensemble

This paper demonstrates that for constrained ensembles of pure Gaussian states modeling thermalization and Hawking radiation, the average entanglement entropy of a small subsystem equals the von Neumann entropy of a corresponding uncorrelated mixed Gaussian state, thereby providing a framework to analyze the Page curve under unitary evolution.

The Gist

The Big Picture: The Black Hole Mystery

Imagine a black hole as a giant, cosmic vacuum cleaner. For decades, physicists have been puzzled by a paradox:

1. The Thermal Clue: When a black hole "eats" matter and eventually evaporates, it spits out radiation (Hawking radiation) that looks exactly like heat from a toaster. It's random, messy, and "thermal."
2. The Quantum Clue: In the quantum world, information can never be truly destroyed. If the black hole started as a neat, organized object (a "pure state"), the radiation it spits out must also be a neat, organized object, even if it looks messy on the surface.

The Problem: How can something look like random noise (thermal) to a small observer, yet be a perfectly organized, information-rich whole (pure) to the universe?

The Authors' Solution: The "Gaussian Ensemble"

The authors (Erik Aurell, Lucas Hackl, and Mario Kieburg) decided to stop guessing and start building a model. They asked: "What if we create a giant library of possible quantum states that all look thermal on the surface, but are actually pure inside?"

They built a mathematical "ensemble" (a collection) of these states. Think of it like a massive bag of marbles.

• The Constraint: Every marble in the bag must have a specific color pattern on its surface (the "marginals"). This ensures that if you look at just one marble, it looks exactly like the thermal radiation Hawking predicted.
• The Freedom: Inside the bag, the marbles can be arranged in any way, as long as they fit the surface pattern.

The Discovery: The "Small Subsystem" Surprise

The researchers wanted to know: If I pick a tiny handful of these marbles (a small subsystem) and look at them, what do I see?

The Analogy: The Orchestra and the Soloist
Imagine a massive orchestra (the whole universe/black hole) playing a symphony.

• The Global View: The whole orchestra is playing a complex, perfectly coordinated piece of music (a "pure state"). Every instrument is listening to every other instrument.
• The Local View: You put on noise-canceling headphones and only listen to the violin section (a small subsystem).

The Result:
The authors proved that if the orchestra is huge and you only listen to a tiny section, the violin section sounds exactly like a group of musicians playing random, uncoordinated jazz.

• No Correlations: The violins aren't "talking" to each other in a complex way. They aren't entangled with their immediate neighbors.
• Maximal Entanglement: However, the violins are maximally entangled with the rest of the orchestra (the drums, the brass, the choir). The "randomness" you hear is actually because the violins are perfectly synchronized with the entire rest of the universe, which you can't hear.

In simple terms: If you look at a small piece of the black hole's radiation, it looks like a mixed-up, thermal mess with no internal connections. But that's only because it's secretly holding hands with everything else in the universe.

Why Does This Matter? (The Page Curve)

This connects to the famous Page Curve, which is a graph showing how much information is "lost" or "entangled" as a black hole evaporates over time.

• The Early Days: When the black hole is young, the radiation it emits is small compared to the black hole itself. The authors' math shows that in this phase, the radiation modes are essentially independent of each other. They are just thermal noise.
• The Implication: This supports the idea that the "Information Paradox" might not be a glitch in physics, but a consequence of geometry. The information isn't gone; it's just hidden in the massive, complex web of connections between the tiny radiation particles and the rest of the black hole.

The "Gaussian" Assumption

The paper uses a specific type of math called "Gaussian states."

• The Metaphor: Imagine a bell curve (the classic "hump" shape in statistics). Most things in nature follow this shape. The authors assume the black hole's radiation follows this simple, smooth bell curve.
• The Caveat: They admit that near the very end of a black hole's life (when it's tiny and gravity is crazy), the math might get weird and non-Gaussian. But for the early and middle stages of a black hole's life, their "bell curve" model works perfectly.

The Takeaway

This paper is a mathematical proof that you don't need complex, exotic physics to explain why black hole radiation looks thermal.

If you have a giant, pure quantum system (the black hole) and you only look at a tiny, random slice of it, that slice will automatically look like a hot, messy, thermal soup. The "messiness" is an illusion caused by your limited view. The information is safe; it's just spread out across the entire universe, connected in a way that makes the small pieces look like they have no secrets of their own.

In a nutshell: The universe is a giant, perfectly organized puzzle. If you only look at a few pieces, they look like random junk. But if you step back, you realize they are part of a masterpiece.

Technical Summary

1. Problem Statement and Motivation

The paper addresses two fundamental problems in theoretical physics:

1. The Black Hole Information Paradox: Specifically, the nature of entanglement in Hawking radiation. While Hawking radiation appears thermal (mixed) locally, unitary evolution requires the global state of an evaporating black hole to be pure. The paradox lies in understanding how a globally pure state can exhibit local thermal properties and what the entanglement structure between radiation modes looks like.
2. Eigenstate Thermalization Hypothesis (ETH): Understanding how isolated quantum systems (pure states) can thermalize locally, appearing as mixed thermal states to local observers.

The authors aim to construct a rigorous statistical ensemble of pure Gaussian states that are constrained to have specific single-mode marginals (thermal states) but are otherwise random. The goal is to calculate the average entanglement entropy of a small subsystem within this ensemble to determine if typical correlations exist between modes or if the system behaves as if it were a product of independent thermal states.

2. Methodology

The authors employ a combination of Random Matrix Theory (RMT), Symplectic Geometry, and the Replica Trick.

• Ensemble Construction:

  • They consider a system of $N$ bosonic modes with a global pure Gaussian state $|\psi\rangle$.
  • The state is characterized by a covariance matrix $C = SS^\top$, where $S$ is a symplectic transformation ($S\Omega S^\top = \Omega$).
  • Constraints: The ensemble is constrained such that the single-mode marginals (reduced density matrices for individual modes) are fixed to specific thermal states. This fixes the diagonal blocks of the covariance matrix.
  • Scenarios:
    • Scenario I: Arbitrary correlations between modes are allowed, provided marginals are fixed.
    • Scenario II: Correlations between modes within specific subsets (e.g., time windows in Hawking radiation) are constrained to vanish, mimicking the "Wald" condition where modes in the same epoch are uncorrelated.
  • The ensemble measure is the natural invariant volume form of the symplectic group restricted to the sub-manifold defined by these marginal constraints.

• Analytical Tools:

  • Coherent State Representation: The reduced density matrix $\rho_A$ of a subsystem $A$ is expressed using the characteristic function and displacement operators.
  • Replica Trick: To compute the average entanglement entropy $S_A = -\text{Tr}(\rho_A \ln \rho_A)$, the authors calculate the average moments $\langle \text{Tr}(\rho_A^x) \rangle$ for integer $x$, and then analytically continue to $x \to 1$.
  • Saddle Point Approximation: The integrals over the symplectic group are evaluated in the limit where the total number of modes $N$ is much larger than the subsystem size $N_A$ ($N_A/N \to 0$). This allows the use of a mean-field approximation where fluctuations around the saddle point are negligible.

3. Key Contributions and Results

A. Equivalence of Average Purity and Mixed State Purity

The central technical result is that for a small subsystem $A$ ($N_A \ll N$), the average purity of the reduced state in the constrained pure Gaussian ensemble is identical to the purity of a mixed Gaussian state constructed solely from the constrained marginals, assuming no correlations between the modes in $A$.

Mathematically, if $\hat{C}_A$ is the block-diagonal constraint matrix (containing only the fixed single-mode covariances) and $C_A$ is the actual reduced covariance matrix of the random pure state:
$$\langle \text{Tr}(\rho_A^2) \rangle = \text{Tr}(\hat{\rho}_A^2)$$
where $\hat{\rho}_A = \bigotimes_{i \in A} \rho_i$ is the product state of the individual thermal modes.

B. Generalization to Higher Moments and Entropy

The authors extend this result to higher moments $\langle \text{Tr}(\rho_A^x) \rangle$ for $x > 2$. They demonstrate that:
$$\langle \text{Tr}(\rho_A^x) \rangle = \text{Tr}(\hat{\rho}_A^x)$$
By analytically continuing $x \to 1$, they derive the average von Neumann entropy:
$$\langle S_A \rangle = \sum_{i \in A} S(\rho_i)$$
Physical Interpretation: The average entanglement entropy of a small subsystem in this ensemble is maximal. It is exactly equal to the sum of the entropies of the individual modes. This implies that, on average, the small subsystem is maximally entangled with the rest of the system (the complement), and there are vanishingly small correlations (entanglement) between the modes within the small subsystem itself.

C. Application to Hawking Radiation

Applying this to the Hawking radiation model:

• Local Thermalism: The model reproduces the thermal nature of individual Hawking modes.
• Absence of Mode-Mode Correlations: The result implies that for a typical state in this ensemble, two distinct Hawking modes (even from the same time window) are almost surely not entangled.
• Global Entanglement: While local modes are uncorrelated, the small set of modes is maximally entangled with the collective of all other modes. This supports the idea that information is encoded in the global correlations between the early radiation and the late radiation (or the black hole interior), rather than in pairwise correlations between early modes.

D. The Page Curve

The authors discuss the implications for the Page curve (entanglement entropy vs. time):

• Their results strictly apply to the initial slope of the Page curve (early times, $N_A \ll N$).
• In this regime, the entropy grows linearly with the number of modes, consistent with the Page curve's initial rise.
• They acknowledge that their mean-field approximation breaks down when $N_A$ becomes comparable to $N$ (the "Page time"), where fluctuations become significant. However, they argue that the initial behavior is robust and likely independent of the specific Gaussian vs. non-Gaussian nature of the state.

4. Significance and Discussion

• Resolution of the "Paradox" in Simple Models: The paper suggests that the apparent paradox of a pure state looking thermal locally is a generic feature of high-dimensional random states with fixed marginals. It does not require exotic physics; simple random matrix theory on the symplectic group suffices to reproduce the thermal properties and the lack of local mode-mode correlations.
• Information Paradox Implications: The findings support the view that the black hole information paradox is primarily a geometric problem related to high-dimensional entanglement structures. If the total state is a random pure Gaussian state consistent with Hawking's marginals, then:
  1. Local observers see thermal radiation (consistent with Hawking's calculation).
  2. No information is lost because the global state is pure.
  3. Information retrieval requires measuring correlations between a small subsystem and the entire rest of the radiation, not just pairwise correlations between modes.
• Limitations: The authors explicitly note that their results are asymptotic ($N \to \infty, N_A/N \to 0$). They do not fully resolve the behavior at the Page time (where $N_A \approx N/2$) or the final evaporation stage, where non-Gaussian effects (due to strong gravity) might become relevant. However, they argue that the initial slope of the Page curve is likely universal for both Gaussian and non-Gaussian random states.

Conclusion

The paper provides a rigorous analytical proof that in an ensemble of pure Gaussian states constrained to have thermal single-mode marginals, the average entanglement entropy of a small subsystem is maximal and equal to the sum of the entropies of its constituent modes. This implies that typical correlations between different modes are negligible, while the subsystem is maximally entangled with the environment. This offers a compelling statistical mechanical explanation for the thermal appearance of Hawking radiation and the structure of the Page curve in the early stages of black hole evaporation.