A two-mode model for black hole evaporation and information flow

We develop and analyze a two-oscillator model for black hole evaporation in which an effective geometric degree of freedom and a representative Hawking radiation mode are described by coupled harmonic oscillators with opposite signs in their free Hamiltonians. The normal-mode structure is obtained analytically and the corresponding modal amplitudes determine the pattern of energy exchange between the two sectors. To bridge the discrete and semiclassical pictures, we introduce smooth envelope functions that provide a continuous effective description along the geometric variable. Numerical simulations in a truncated Fock space show that the two oscillators exchange quanta in an approximately out-of-phase manner, consistent with an effective conservation of $\langle n_x\rangle - \langle n_y\rangle$. The reduced entropy $S_x(t)$ exhibits periodic growth, indicating entanglement generation. These results demonstrate that even a minimal two-mode framework can capture key qualitative features of energy transfer and information flow during evaporation.

💡 Research Summary

The paper presents a minimalist yet insightful quantum‑mechanical model of black‑hole evaporation using two coupled harmonic oscillators. One oscillator (denoted x) represents an effective geometric degree of freedom of the black hole (e.g., its mass or radius), while the second oscillator (denoted y) stands for a representative Hawking‑radiation mode. The crucial ingredient is that the free Hamiltonians of the two oscillators carry opposite signs, so that energy lost by the x‑sector appears as energy gained by the y‑sector. The total Hamiltonian reads

 H = ½(pₓ² + ωₓ²x²) − ½(pᵧ² + ωᵧ²y²) + g x y,

where g is a linear coupling constant. This sign structure mimics the classical picture of a black hole radiating away its mass.

The authors first solve the classical equations of motion, casting them into a symplectic matrix form. By seeking solutions of the type ξ(t)∝e^{iΩt} they obtain the secular equation

 (Ω² − ωₓ²)(Ω² − ωᵧ²) + g² = 0,

which yields two normal‑mode frequencies Ω₁ and Ω₂. The discriminant Δ = (ωₓ² − ωᵧ²)² − 4g² determines the nature of the dynamics: Δ ≥ 0 gives real Ω’s and stable oscillatory exchange; Δ < 0 leads to one complex pair, signalling exponential growth/decay—interpreted as a “run‑away” regime where the geometric sector collapses rapidly. This analytic classification provides a clear criterion for stability in terms of the coupling strength and the bare frequencies.

Quantization proceeds by introducing ladder operators aₓ, aᵧ. The difference operator D = aₓ†aₓ − aᵧ†aᵧ is not conserved by the interaction, but for weak coupling it remains approximately constant, reflecting an effective conservation of ⟨nₓ⟩ − ⟨nᵧ⟩. The total excitation number Q = aₓ†aₓ + aᵧ†aᵧ is not conserved because of the opposite signs in the free Hamiltonians.

The initial quantum state is chosen as a product of a coherent state |α⟩ for the black‑hole oscillator and the vacuum |0⟩ for the radiation oscillator. The authors truncate the infinite‑dimensional Fock space to a finite cutoff N_cut, diagonalize the Hamiltonian exactly, and evolve the state via |Ψ(t)⟩ = e^{−iHt}|Ψ(0)⟩. Observable quantities—occupation numbers ⟨n_i⟩, quadrature variances ⟨x²⟩, ⟨p²⟩, and the reduced von‑Neumann entropy Sₓ(t)=−Tr