Two-Oscillator Black Hole Evaporation

• The paper introduces a quantum two-oscillator model that encapsulates energy exchange and entanglement dynamics, reproducing the qualitative features of the Page curve.
• It employs Hamiltonian constructions from optical formalisms and matrix mechanics to represent black hole geometry and Hawking radiation with clear energy conservation and unitary evolution.
• The framework offers analytically tractable insights into reversible information flow while addressing key limitations such as neglecting an infinite radiation bath and fully dynamic backreaction.

A two-oscillator model for black hole evaporation is a quantum mechanical framework in which the complex process of evaporation and information transfer is distilled into a system of two coupled harmonic oscillators that encode, respectively, certain geometric properties of the black hole and a representative degree of freedom of Hawking radiation. This class of toy models aims to capture the core phenomenology of quantum black hole evaporation, such as energy exchange, entanglement growth, the qualitative structure of the Page curve, and operational aspects of information flow, within a tractable and analytically transparent dynamical system. Recent developments extend these models to include optical formalism, Hamiltonian constructions with atomic-like microstate degeneracy, and matrix mechanics reductions; these approaches are motivated by both operational simplicity and the need for models that rigorously encode unitarity, entanglement, and effective macroscopic behavior (Alsing, 14 Jan 2026, Bayenat et al., 27 Jan 2026, Zeng, 2021, Berenstein et al., 2021).

1. Model Hamiltonians and Mode Assignments

Canonical two-oscillator models introduce two quantum harmonic oscillators:

• One oscillator represents a coarse-grained geometric or pump variable associated with the black hole (often a degree of freedom analogous to internal energy, mass, collective geometry, or matrix eigenvalue sector).
• The second oscillator models a single representative mode of Hawking radiation, or, in generalizations, a selected sector of the exterior radiation field.

In the model developed in "A two-mode model for black hole evaporation and information flow" (Bayenat et al., 27 Jan 2026), the Hamiltonian is

$$H = \hbar\omega_x\, a^\dagger_x a_x - \hbar\omega_y\, a^\dagger_y a_y + \frac{g}{2}(a_x + a_x^\dagger)(a_y + a_y^\dagger),$$

where $a_x, a_x^\dagger$ and $a_y, a_y^\dagger$ are annihilation/creation operators for the black hole ("geometry") and radiation modes, respectively. The negative sign in the radiation sector ensures that quanta lost by $x$ are gained by $y$, reflecting effective energy conservation between subsystems.

In matrix-inspired models (Berenstein et al., 2021), two degrees of freedom emerge from Yang-Mills reductions: a collective coordinate $R(t)$ (e.g., brane separation) and a fast off-diagonal oscillator $x(t)$ encoding light string excitations, with an effective Hamiltonian

$$H_{\mathrm{eff}}(R, x) = \frac{1}{2}p_R^2 + \frac{p_\phi^2}{2R^2} - aR + \frac{p_x^2}{2m_x} + \frac{1}{2}m_x\omega^2(R)x^2$$

where the oscillator frequency depends on $R$.

Gaussian optics-inspired models (Alsing, 14 Jan 2026) represent the black hole as a single-mode squeezed state and the radiation as a vacuum mode, interacting via beam splitters and squeezers, entirely within the symplectic (covariance-matrix) formalism.

2. Normal-Mode Structure and Energy Transfer

The quadratic form of the two-oscillator Hamiltonians allows for exact diagonalization via a Bogoliubov or symplectic transformation. With

$$X_\pm = \cos\theta\,x \pm \sin\theta\,y,\qquad P_\pm = \cos\theta\,p_x \pm \sin\theta\,p_y,$$

the mixing angle $\theta$ satisfies $\tan2\theta = 2g/(\omega_x^2 - \omega_y^2)$, and the decoupled Hamiltonian is

$$H = \frac{1}{2}\left(P_+^2 + \Omega_+^2 X_+^2\right) - \frac{1}{2}(P_-^2 + \Omega_-^2 X_-^2),$$

with squared normal-mode frequencies specified by

$$(\Omega^2 - \omega_x^2)(\Omega^2 - \omega_y^2) + g^2 = 0.$$

Time evolution yields out-of-phase energy exchange:

$$\langle n_x(t)\rangle \approx |\alpha|^2 \cos^2\left(\frac{1}{2}\Delta\Omega\, t\right), \quad \langle n_y(t)\rangle \approx |\alpha|^2 \sin^2\left(\frac{1}{2}\Delta\Omega\, t\right),$$

where $\Delta\Omega = |\Omega_+ - \Omega_-|$, indicating that occupation transfer between black hole and radiation oscillators is coherent and oscillatory in this minimal setup. The difference $$D(t) = \langle n_x(t)\rangle - \langle n_y(t)\rangle$$ remains approximately conserved over a beat period, a consequence of the quadratic and unitary evolution (Bayenat et al., 27 Jan 2026).

3. Entanglement Generation and the Page Curve

The full quantum state of the two-oscillator system remains pure, but subsystems (such as the black hole oscillator or the radiation oscillator) become entangled due to the interaction. The reduced density matrix for, e.g., the $x$-oscillator subsystem is constructed by tracing over the radiation sector:

$$\rho_x(t) = \mathrm{Tr}_y\left[|\Psi(t)\rangle\langle\Psi(t)|\right].$$

For Gaussian initial states and quadratic evolutions, $\rho_x(t)$ remains Gaussian, and the symplectic eigenvalue $\nu(t)$ of its covariance matrix determines the subsystem entropy by

$$S_x(t) = \left(\nu + \frac{1}{2}\right) \ln\left(\nu + \frac{1}{2}\right) - \left(\nu - \frac{1}{2}\right) \ln\left(\nu - \frac{1}{2}\right).$$

Numerical and analytic results show that $S_x(t)$ (and similarly $S_y(t)$) exhibits periodic growth, peaking when oscillator occupations equalize, directly analogous to the rising limb of the Page curve. Subsequent decline in $S_x(t)$ as energy recoheres into a single oscillator mimics the late-time descent of the Page curve, providing dynamical evidence for unitary information recovery in the toy evaporation process (Bayenat et al., 27 Jan 2026, Alsing, 14 Jan 2026).

In the Gaussian optics model, this Page-curve-like behavior is achieved by dynamically tuning the beam-splitter angle such that the reflectivity $R(t) = \sin^2\theta(t)$ evolves from no leakage, to half-occupancy (maximal entanglement), then back to total transfer (final purity), with the entropy $S_A(t)$ displaying the characteristic rise and fall (Alsing, 14 Jan 2026).

4. Physical Interpretations and Limitations

The two-oscillator paradigm captures central qualitative phenomena:

• Energy Flow: Oscillatory and beat-frequency transfer of occupation between black hole and radiation modes; reflective of effective energy conservation and quantum backreaction.
• Entanglement: Generation and periodic reduction of entanglement entropy, reproducing the Page curve profile, though only strictly in the minimal two-mode system.
• Information Flow: Reversible unitary dynamics ensure that no information is irretrievably lost—a central requirement in candidate resolutions of the black hole information paradox.

Limitations of the two-mode models include:

• Real Hawking evaporation involves an infinite bath of radiation modes, while most two-oscillator models capture only a single mode, so do not exhibit genuine irreversibility or continuous thermal emission.
• The quadratic time-independent coupling cannot fully describe long-time evaporation, late-time decay to zero mass, or realistic decoherence.
• Backreaction and geometric response are only approximately encoded or modeled via envelope functions or semiclassical estimation, not from first principles (Bayenat et al., 27 Jan 2026).

As such, these models serve as analytically tractable laboratories rather than full descriptions.

5. Extensions: Multi-Level Oscillator and Atomic Analogies

More sophisticated approaches generalize these models to multi-level, degenerate systems, inspired by spontaneous emission in quantum optics. In (Zeng, 2021), the black hole is modeled as a multi-level "oscillator" (with degeneracy given by the Bekenstein-Hawking entropy), coupled by monopole transitions to the radiation field:

$H = H_{BH} + H_{rad} + H_{int},$

with $H_{BH}$ diagonal in the black hole microbasis and $H_{int}$ encoding transitions mediated by quantum emission.

Using the Wigner–Weisskopf approximation, the first emission reproduces the exact Hawking spectrum. For multiple emissions, numerical integration of the full Schrödinger equation yields a Page curve for the reduced entropy that rises and then decreases, matching the expected behavior from unitarity.

This formalism also supports an "atomic-like" inner structure: Families of solutions to Einstein's equations represent quantized spherically symmetric matter configurations, with the total number of microstates growing exponentially with the area, thereby recovering the Bekenstein-Hawking law within a quantized dust shell model.

Crucially, the framework demonstrates that information can escape via spontaneous monopole transitions (not pair creation at a sharp horizon), and that the unitary Page curve is a consequence of macroscopic quantum superpositions over distinct mass and inner-structure branches. This suggests a concrete mechanism to evade the small-correction no-go theorems and provides a non-replica-based resolution of the information paradox (Zeng, 2021).

6. Matrix-Mechanics and Adiabatic Invariants in Evaporation

Reduction of matrix quantum mechanics—specifically, 2×2 Yang-Mills models—yields an effective two-oscillator system composed of a "slow" radial variable $R(t)$ (the brane separation) and a "fast" off-diagonal string excitation $x(t)$ (Berenstein et al., 2021). At large $R$, the dynamics become adiabatic; the effective occupation number $N = H_{\mathrm{osc}} / \omega(R)$ is approximately conserved, and serves as an adiabatic invariant regulating evaporation.

Classically, the average return (lifetime before evaporation) diverges logarithmically due to rare, long excursions at threshold $N\to N_{\min}$. Quantum mechanically, discretization stabilizes the lifetime at $\tau \sim \log(1/\hbar)$. This model thus captures both the slow, quasi-irreversible leakage of black hole mass in semiclassical evaporation and the essential quantum regulation of evaporation timescales. A plausible implication is that adiabatic invariants play a critical role in determining long-time quantum decay properties of gravitational bound states (Berenstein et al., 2021).

7. Operational and Conceptual Significance

Two-oscillator models have established an impactful methodological niche:

• They connect quantum optics, information theory, and black hole thermodynamics within a non-perturbative, exactly solvable formalism (Alsing, 14 Jan 2026, Zeng, 2021).
• They provide explicit realizations of qualitative features—energy transfer, entanglement dynamics, and partial recoherence—that are otherwise obscured in infinite-mode, path-integral, or semiclassical analyses.
• By distilling the evaporation process to minimal degrees of freedom, these models clarify which features are generic (e.g., Page-curve entanglement) and which require greater microscopic detail (e.g., genuine thermalization and long-term decay).

Although limitations exist—primarily regarding the inclusion of many radiation modes and dynamical gravity—the two-oscillator framework continues to be an essential tool for testing and developing ideas concerning quantum gravity, black hole microstructure, and information recovery mechanisms.