Population-coherence routes to purity in Page-type models of black-hole evaporation

1. Problem Statement

The black-hole information paradox centers on the apparent conflict between the unitary evolution of quantum mechanics (which preserves information and purity) and Hawking's semiclassical prediction that black holes emit thermal radiation, leading to a mixed state and information loss. While recent developments (e.g., the "island formula" and Page curves) suggest that the entanglement entropy of radiation follows a unitary trajectory—rising to a maximum at the Page time and then decreasing back to zero—these analyses primarily rely on entropic measures (von Neumann entropy).

The author argues that entropy alone is insufficient to describe how information is encoded in the radiation. A density matrix $\rho$ contains more structural information than a single entropy number. Specifically, the paper asks: When the radiation recovers purity at late times, is this purity stored in the level populations (diagonal elements) or in the coherences (off-diagonal elements)?

The central hypothesis is that if the radiation remains locally thermal (i.e., its energy populations remain close to a thermal distribution), then the recovery of global purity cannot be driven by population anisotropy. Instead, it must be driven by the build-up of quantum coherences.

2. Methodology

The paper employs a population–coherence decomposition of density-matrix purity, adapted from polarization optics and general density-matrix analysis.

A. Theoretical Framework

1. Intrinsic Basis Decomposition:
   For an $n$-dimensional density matrix $\rho$, the author defines an "intrinsic basis" via a real orthogonal transformation $Q$ that diagonalizes the real part of the density matrix, $\Re(\rho)$.
   $$\rho_O = Q \rho Q^\dagger = A + iN$$
   Where:

   • $A = \Re(\rho_O)$ is a real diagonal matrix containing intrinsic populations $\{a_k\}$.
   • $N = \Im(\rho_O)$ is a real antisymmetric matrix containing intrinsic coherences $\{n_{ij}\}$.

2. Purity Indices:
   The standard global purity $P(n)$ is normalized to range from 0 (maximally mixed) to 1 (pure). The author decomposes the squared global purity into two orthogonal components:
   $$P^2(n) = P_p^2 + \alpha_n P_c^2$$
   Where:

   • $P_p$ (Population Index): Measures the anisotropy of the intrinsic populations relative to a uniform distribution.
   • $P_c$ (Coherence Index): Measures the strength of the intrinsic off-diagonal correlations.
   • $\alpha_n = \frac{n}{2(n-1)}$ is a dimension-dependent scaling factor.

3. Geometric Representation:
   By defining a rescaled coherence index $\tilde{P}_c = \sqrt{\alpha_n} P_c$, the decomposition becomes a Euclidean relation:
   $$P^2(n) = P_p^2 + \tilde{P}_c^2$$
   This allows states to be plotted in a Population–Coherence Plane $(P_p, \tilde{P}_c)$.

   • Population-dominated route: Trajectory moves toward $(1, 0)$. Purity comes from non-thermal diagonal populations.
   • Coherence-dominated route: Trajectory moves toward $(0, 1)$. Purity comes from off-diagonal correlations while populations remain uniform/thermal.

B. Application to Black Hole Evaporation

The framework is applied to a Page-type model of black-hole evaporation:

• System: A bipartite system consisting of a Black Hole (BH) and Radiation (Rad) with Hilbert space dimensions $d_{BH}(k)$ and $d_{rad}(k)$ evolving as $k$ qubits are emitted.
• Assumption: The global state is a typical pure state drawn from the Haar measure (representing chaotic scrambling).
• Constraint: The analysis assumes the radiation populations in the physical energy basis remain locally thermal (nearly uniform in the simplified model, or close to a Gibbs distribution in realistic extensions). This implies $P_p$ must remain small throughout the process.

3. Key Contributions and Results

A. Two-Level Illustration (The Prototype)

The author constructs two explicit families of 2-level states with identical eigenvalues (and thus identical global purity $P$) but opposite internal structures:

1. Route A (Population-Dominated): Purely diagonal states where $P_c = 0$ and $P_p = P$. Information is encoded in the deviation of emission probabilities from a thermal spectrum.
2. Route B (Coherence-Dominated): States where intrinsic populations are exactly equiprobable ($P_p = 0$) and $P_c = P$. Information is encoded entirely in off-diagonal coherences.

• Result: This proves that global purity and spectrum do not uniquely determine the encoding of information.

B. Large-Dimension Page Model Analysis

Applying the decomposition to the Page model with $N$ qubits:

1. Typical Behavior: For typical reduced states in large dimensions, the diagonal elements (populations) of the radiation density matrix fluctuate very little around the uniform distribution ($1/n$).
2. Constraint on $P_p$: Under the assumption of locally thermal populations, the population index $P_p(k)$ remains small ($P_p \approx 0$) throughout the evaporation process.
3. Necessity of Coherence: Since the global purity $P(k)$ must recover to 1 at late times (unitary evaporation), and $P_p$ is constrained to be small, the decomposition $P^2 = P_p^2 + \tilde{P}_c^2$ forces $\tilde{P}_c$ to approach 1.
4. Geometric Trajectory: The trajectory of the radiation state in the $(P_p, \tilde{P}_c)$ plane starts at the origin (maximally mixed), stays near the origin around the Page time, and then moves vertically to the point $(0, 1)$.
   • Conclusion: The late-time recovery of purity in unitary evaporation models with locally thermal populations is necessarily coherence-dominated.

C. Robustness and Generalization

• Dynamical Models: The author argues this result is robust beyond the static Haar-random ensemble. In chaotic dynamical models (random unitary circuits, fast scramblers), the reduced state is statistically indistinguishable from a Haar-random state, preserving the concentration of measure that keeps $P_p$ small.
• Island Formula Connection: The paper suggests that "islands" and replica wormholes in gravitational path integrals generate the necessary nonlocal correlations. These manifest as intrinsic coherences ($P_c$) in the reduced radiation density matrix, purifying the state without altering the coarse-grained thermal spectrum.

4. Significance and Implications

1. Refining the Information Paradox: The paper shifts the focus from whether information is recovered (entropy) to how it is recovered (structural decomposition). It provides a precise mathematical language to distinguish between "thermal deviations" and "quantum correlations."
2. Diagnostic Tool: The $(P_p, \tilde{P}_c)$ plane serves as a phase diagram for information recovery routes. It allows researchers to classify evaporation models:
   • Models ending near $(1,0)$ imply information escapes via spectral deviations (non-thermal emission).
   • Models ending near $(0,1)$ imply information escapes via subtle correlations (coherence) while the spectrum remains thermal.
3. Validation of Semiclassical Intuition: The result supports the intuition that Hawking radiation can appear thermal locally (small $P_p$) while the global state remains pure, provided the "missing" information is stored in high-order correlations (large $P_c$).
4. Future Directions: The framework suggests extensions to "thermal reference states" (deviations from Gibbs states rather than uniform distributions) and higher-order statistical descriptors (beyond second-order purity) to further characterize the fine structure of quantum gravity states.

Summary Conclusion

José J. Gil demonstrates that in standard Page-type models of black-hole evaporation, if the radiation maintains locally thermal populations (a physically motivated constraint), the restoration of global purity at late times must be driven by the growth of intrinsic coherences. The information is not recovered by the radiation becoming spectrally non-thermal; rather, it is recovered through the build-up of complex, non-local quantum correlations that are invisible to diagonal measurements but essential for unitarity. This "coherence-dominated route" is a robust kinematic feature of unitary evaporation in chaotic systems.