PSSY Model in Quantum Gravity

• The PSSY Model is a theoretical framework in low-dimensional JT gravity utilizing EOW branes to analyze black hole evaporation and entanglement dynamics.
• It employs gravitational replica tricks and the island rule to compute von Neumann and Rényi entropies, revealing transitions between disconnected and connected saddle points.
• The model connects diagrammatic resummation, free probability, and quantum complexity—linking random matrix theory with stabilizer properties and corrections to holographic entropy.

The term "PSSY model" refers to the Penington–Shenker–Stanford–Yang framework, a class of analytically tractable toy models for quantum gravity in low dimensions, particularly designed to elucidate quantum aspects of black hole evaporation, entanglement entropy, and the Page curve. Initially formulated to address the microstate dynamics and entanglement structure in two-dimensional Jackiw–Teitelboim (JT) gravity with end-of-the-world (EOW) branes, PSSY constitutes the standard setup for replica wormhole analyses and fine-grained gravitational entropy calculations. It offers a complete semiclassical and quantum path-integral realization of the "island rule" and associated quantum extremal surface (QES) prescriptions, providing canonical cases for the study of dynamical entanglement transitions, quantum error correction in AdS/CFT, and corrections to leading-order holographic entropy formulas.

1. Mathematical Definition and Core Structure

The PSSY model comprises JT gravity—with action

$$I[g,\phi] = -S_0 \chi - \frac{1}{4\pi}\left[\int_M \sqrt{g}\,\phi (R+2) + \int_{\partial M} \sqrt{h}\,\phi K\right] + \mu \int_{\rm EOW} ds$$

where $S_0$ is extremal black hole entropy and $\mu$ is the brane tension—coupled to a non-dynamical quantum system ("bath" or "reference" $R$) via maximally entangled EOW brane degrees of freedom of rank $k$. The joint Hilbert space is

$$\mathcal{H}_{\rm full} = \mathcal{H}_{\mathrm{BH}} \otimes \mathcal{H}_R$$

with microstates $|\psi_i\rangle_{\mathrm{BH}} \otimes |i\rangle_R$, and the total pure state

$$|\Psi\rangle = \frac{1}{\sqrt{k}} \sum_{i=1}^k |\psi_i\rangle_{\mathrm{BH}} \otimes |i\rangle_R$$

prepared by a JT gravitational path integral with boundary EOW brane insertions and appropriate ensemble averages over interior microstates. This architecture situates the radiation sector as a subsystem entangled across a QES with the gravitationally-dressed black hole region.

2. Gravitational Replica Trick and Island Rule

The primary computational machinery in the PSSY model involves gravitational path integrals implementing the $n$-replica trick, allowing direct and exact calculation of von Neumann and Rényi entropies for various radiation subsystems. The partition function for $n$ boundaries is given, microcanonically, by

$Z_n(E) = e^{S_0} y(E)^n$

with $y(E)$ encapsulating brane factors and spectral density $\rho(E)$. Genus-zero path integrals sum over all possible genus-zero topologies (non-crossing brane diagrams), each associated with permutation group structure.

A key structural feature is the emergence of replica wormholes: topologically nontrivial solutions connecting multiple boundaries, yielding new saddle points not present in fixed spacetime background calculations. The island formula for fine-grained entropy reads

$$S(\rho_R) = \min_{\text{extremal}}\left\{ \frac{\phi(\partial I)}{4G_N} + S_{\mathrm{matter}}(R \cup I) \right\}$$

where $I$ is a candidate island region, $\partial I$ the QES. The formula is realized explicitly in the PSSY via a competition between disconnected ("Hawking") and connected (wormhole) path-integral geometries (Anderson et al., 2021). The dominance shifts at a critical $k \sim e^{S_0}$, producing a nonperturbative Hawking–Page-like entropy transition.

3. Statistical and Probabilistic Frameworks

The combinatorics of the PSSY replica path integral naturally dictates a correspondence with free probabilistic multiplicative convolution. For joint radiation-black hole spectral density, the moment expansion matches

$\nu_{\tilde{r}} = \nu_r \boxtimes \mu_{\nu_b}$

where $\nu_r$ is the reservoir (bath) spectrum and $\mu_{\nu_b}$ a compound Poisson law for the gravitational sector (Wang, 2022). The associated Cauchy transform, $S$-transform, and $R$-transform enable exact evaluation of the joint entropy and spectral density through harmonic analysis and random matrix theory, generalizing Page's formula and justifying solution uniqueness via the Hausdorff moment problem.

This interpretation extends to random tensor network representations and fixed-area state decompositions, situating the PSSY setup as a canonical model for understanding free random variables and convolution algebra in emergent von Neumann algebras induced by quantum gravity path-integral topology.

4. Diagrammatics and Matrix Model Correspondence

The diagrammatics of the PSSY model—planar and non-planar Feynman diagrams—encode the gravitational topologies contributing to entanglement entropy. The planar (disk) diagrams produce the leading-order Page curve, while non-planar crossing diagrams introduce $1/k^2$ corrections. Through an explicit mapping to the Iizuka–Okuda–Polchinski (IOP) matrix model, these non-planar crossing diagrams and their resummation can be matched to matrix model two-point functions, with the result

$$S_R = \ln k - \frac{1}{2} + \frac{1}{12k^2} + O(k^{-4})$$

when $k = e^{S_0}$, paralleling Page's original formula but missing additional $1/2k^2$ corrections contributed by bulk handles (Iizuka et al., 2024). This relationship situates gravitational topological transitions as diagrammatic phenomena in matrix integrals, clarifying the division between "crossings" and "handles" in contributing quantum gravity loop corrections.

5. Entropic Phase Structure and Corrections

Transition phenomena in the PSSY model are classified according to the competing path-integral saddles:

• Regime I (Low Entanglement, $k \ll e^{S_0}$): Disconnected diagrams dominate (empty islands), entropy determined by the boundary bath rank.
• Regime II (High Entanglement, $k \gg e^{S_0}$): Euclidean wormhole saddles dominate, corresponding to a two-sided black hole with an ER bridge—nonperturbative topology change induced by entanglement.
• Intermediate Corrections: Subleading genus-zero permutations ($\mathcal{O}(e^{-S_0})$), area fluctuations ($\mathcal{O}(G_N^{-1/2})$), incompressible bulk state corrections ($\mathcal{O}(G_N^{-1})$), and infinite violations due to rank-deficiency in bulk entanglement spectrum (Kudler-Flam et al., 2022).

Diagrammatic resummation and random matrix analogs precisely enumerate these corrections, with extra-handle diagrams yielding quantum gravity loop effects of $\mathcal{O}(1/k^2)$ asymptotically. Modified cosmic brane prescriptions further refine the treatment of Rényi entropies, particularly highlighting new diagonal “tie” saddles for $n < 1$ not visible in minimax approaches (Dong et al., 2023).

6. Quantum Complexity and Stabilizer Properties

Beyond entropic measures, the PSSY model enables quantitative evaluation of quantum complexity in terms of stabilizer theory, specifically Wigner negativity as a monotone for magic and simulation difficulty. Gravitational path-integral computation of negativity reveals $O(1)$ before the Page transition (maximally mixed phase) and exponentially large values, $\sim \exp[(S_{\rm max} - S_2)/4]$, after the transition (complex, non-stabilizer phase) (Basu et al., 21 Oct 2025). The universal formula,

$$\mathcal{N} = \text{erf}(\sqrt{r}) + \frac{e^{-r}}{\sqrt{\pi r}}, \quad r = \frac{e^{S_2}}{2e^{S_{\rm max}}}$$

interpolates between these regimes, with a geometric interpretation via the area difference of competing extremal surfaces ($\mathcal{N} \sim \exp[\Delta A/(8G_N)]$). This supports the "Python’s lunch" conjecture: entanglement wedge reconstruction and quantum gravity protection scale exponentially with QES area gaps.

Summary Table: Key Aspects of the PSSY Model

Aspect	PSSY Realization	Reference
Black hole evaporative setup	JT gravity + EOW brane	(Kudler-Flam et al., 2022)
Entropy calculation	Replica wormholes, island rule	(Anderson et al., 2021Wang, 2022)
Random matrix theory	Voiculescu free convolution	(Wang, 2022Iizuka et al., 2024)
Entropic corrections	Non-planar diagrams, bulk handles	(Iizuka et al., 2024)
Rényi entropy phases	Maximin cosmic brane saddles	(Dong et al., 2023)
Quantum complexity	Wigner negativity / stabilizer	(Basu et al., 21 Oct 2025)

The PSSY model stands as the paradigmatic analytic framework for resolving quantum gravitational puzzles in low-dimensional evaporating black holes, enabling precise entropic, diagrammatic, probabilistic, and complexity-theoretic analyses with implications extending to emergent spacetime, holographic reconstruction, and statistical properties of quantum gravity solutions.