⚛️ high-energy theory

On the stabilizer complexity of Hawking radiation

This paper investigates the stabilizer complexity of Hawking radiation by calculating Wigner negativity in various black hole evaporation models, demonstrating that it remains low before the Page transition but grows exponentially afterward, and proposing a geometric formula linking this complexity to the presence of a "python's lunch" in the entanglement wedge.

Original authors: Ritam Basu, Onkar Parrikar, Suprakash Paul, Harshit Rajgadia

Published 2026-08-04

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Ritam Basu, Onkar Parrikar, Suprakash Paul, Harshit Rajgadia

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe as a giant, cosmic computer. For decades, physicists have been trying to figure out how this computer processes the most mysterious data of all: black holes. When a black hole "evaporates" (shrinks away by emitting radiation), it seems to delete information, which breaks the fundamental rules of quantum mechanics. This is the famous "black hole information paradox." Recently, scientists discovered that the information isn't actually lost; it's hidden in a very tricky way. The radiation coming out of the black hole is entangled with the inside of the black hole, but only after a specific moment in time called the "Page time." Before this moment, the radiation looks simple and random. After this moment, the radiation holds the secret code to the black hole's interior. But here's the catch: just because the information is there doesn't mean we can easily read it. It might be hidden behind a wall of computational difficulty, so complex that even a supercomputer couldn't crack the code in the lifetime of the universe.

This is where the concept of "complexity" comes in. In the world of quantum computing, there's a special class of problems that are easy for classical computers to solve, and others that are incredibly hard. The "hard" ones require something called "magic" (a technical term for non-classical resources) to solve. If a quantum state has a lot of "magic," it's like a locked safe that requires a very specific, difficult key to open. The paper you are about to read dives into exactly how "magical" (or complex) the radiation from an evaporating black hole is, using a tool called "Wigner negativity" to measure it.

The Paper: Measuring the "Magic" of Black Hole Smoke

In this study, the authors Ritam Basu, Onkar Parrikar, Suprakash Paul, and Harshit Rajgadia from the Tata Institute of Fundamental Research in India set out to answer a burning question: How hard is it to simulate the radiation coming from an evaporating black hole on a regular computer?

To do this, they use a concept from quantum computing called stabilizer complexity. Think of a quantum computer as a chef trying to cook a meal. Some ingredients are "stabilizers"—they are like basic salt and water. You can mix them together in any way, and a regular computer can easily predict the result. But to make a truly complex, delicious quantum dish, you need "magic" ingredients. These are special resources that make the recipe impossible for a standard computer to follow. The more "magic" you need, the harder the recipe is to simulate.

The authors use a mathematical tool called Wigner negativity to measure this "magic." Imagine the Wigner function as a map of a quantum state. In a simple, easy-to-simulate state, this map looks like a nice, positive hill. But in a complex, "magical" state, the map dips below zero, creating "negative valleys." The deeper and more numerous these negative valleys are, the more "magic" the state has, and the harder it is for a classical computer to simulate.

The Journey Through Time: Before and After the Page Transition

The paper explores two main scenarios to see how this "negativity" changes as a black hole evaporates.

1. The PSSY Model: A Static Snapshot
First, the authors look at a simplified model (the PSSY model) where a black hole is already in equilibrium with a bath of radiation. They use the "gravitational path integral"—a fancy way of summing up all possible shapes of spacetime—to calculate the negativity.

Before the Page Transition: When the black hole is young and hasn't evaporated much, the radiation is simple. The Wigner negativity is small (order 1). This means the radiation is like a basic soup; a classical computer can simulate it easily.
After the Page Transition: Once the black hole passes the "Page time" (the halfway point of its life), things get wild. The Wigner negativity explodes, becoming exponentially large. The formula they find is roughly:
$N \sim \sqrt{\frac{2}{\pi}} \exp\left[\frac{1}{2}(S_{\text{max}} - S_2)\right]$
Here, $S_{\text{max}}$ is the maximum possible entropy, and $S_2$ is a measure of how much information is already in the radiation. The result suggests that after the Page time, the radiation is so "magical" that simulating it would require a computer with resources that grow exponentially. It's like trying to simulate a universe inside a universe.

2. The Dynamical Model: A Real-Time Movie
Next, they move to a more realistic, dynamic model where the black hole and the radiation are interacting in real-time, like two dancers coupling up.

The "Negativity Shock": As soon as the black hole starts interacting with the radiation, there is a sharp spike in negativity. It's a sudden "shock" of complexity caused by the coupling.
Settling Down: As time goes on, the system settles. Surprisingly, the negativity eventually drops down and matches the exact same universal formula found in the static PSSY model. This confirms that the exponential complexity is a robust feature of the late-time radiation, regardless of how the system started.

The Geometric Connection: The Python's Lunch

Finally, the authors connect this quantum complexity to the shape of spacetime itself. In the holographic view of the universe, the interior of a black hole is encoded on the surface of the radiation. Sometimes, the geometry of this encoding looks like a "Python's Lunch"—a region where the space bulges out, creating a bottleneck.

The paper proposes a geometric formula for the Wigner negativity:
$N \sim \sqrt{\frac{2}{\pi}} \exp\left[\frac{1}{2} \frac{A(\gamma_{\text{out}}) - A(\gamma_{\text{min}})}{4G_N}\right]$
In this formula, $A(\gamma_{\text{out}})$ is the area of the outermost surface, and $A(\gamma_{\text{min}})$ is the area of the smallest surface (the "bottleneck"). The difference between these areas, divided by Newton's constant $G_N$ , determines the complexity.

What this means: If there is a "Python's Lunch" (a big gap between the outer and inner surfaces), the stabilizer complexity is exponentially large. This suggests that the reason we can't easily manipulate the black hole interior from the outside isn't just because the information is hidden, but because the "key" to unlock it is buried behind a computational wall that is exponentially high.

The Verdict

The authors don't claim to have solved the black hole information paradox (that was done by others using similar tools). Instead, they provide a new way to measure the difficulty of the problem. They show that:

Before the Page time, the radiation is simple and easy to simulate.
After the Page time, the radiation becomes exponentially complex, making it practically impossible for classical computers to simulate.
This complexity is directly linked to the geometry of spacetime, specifically the presence of "Python's Lunch" regions.

In short, the universe seems to have a built-in safety mechanism. While the information inside a black hole is never truly lost (it's just encoded in the radiation), the "magic" required to decode it grows so fast that it protects the semi-classical nature of spacetime from being easily manipulated by observers at infinity. The black hole isn't just a cosmic vacuum cleaner; it's a cosmic firewall made of pure computational complexity.

Technical Summary: On the Stabilizer Complexity of Hawking Radiation

Problem Statement
The black hole information paradox posits a conflict between the unitary evolution of quantum mechanics and the thermal nature of Hawking radiation, which suggests a loss of information. Recent resolutions involving replica wormholes and the emergence of the Page curve indicate that the black hole interior is encoded in the early radiation. However, this encoding implies a form of non-locality where an observer at infinity could, in principle, manipulate degrees of freedom deep inside the black hole. The Harlow-Hayden conjecture suggests that while such operations are information-theoretically possible, they are computationally intractable for observers with bounded resources. This paper investigates the "stabilizer complexity" of Hawking radiation—a measure of the difficulty of classically simulating the quantum state—to quantify this computational protection. Specifically, the authors utilize Wigner negativity as a magic monotone to measure the non-stabilizer content of the radiation state.

Methodology
The authors employ the framework of stabilizer quantum computation, where circuits composed of Clifford unitaries and stabilizer states can be efficiently simulated classically (Gottesman-Knill theorem). States requiring "magic" (non-stabilizer resources) exhibit Wigner negativity, which serves as a proxy for classical simulation complexity. The study proceeds through three distinct approaches:

Gravitational Path Integral in the PSSY Model: The authors analyze the Penington-Shenker-Stanford-Yang (PSSY) model of black hole evaporation in Jackiw-Teitelboim (JT) gravity coupled to an external bath. They compute the Wigner negativity of the reduced density matrix of the radiation ( $\rho_R$ ) directly using the gravitational path integral. This involves calculating ensemble averages over Hamiltonians and utilizing replica wormholes (specifically $n=2$ saddles) to handle the absolute value in the negativity definition via an integral representation and analytic continuation.
Dynamical Model with Random Matrices: A time-dependent model is constructed where a black hole system $B$ is coupled to a radiation bath $R$ via a joint Hamiltonian $H_{BR} = O_B \otimes O_R$ . The operator $O_R$ is modeled as a random matrix drawn from a unitarily invariant ensemble. The authors compute the ensemble-averaged Wigner negativity of the bath system over time using Haar integration techniques and diagrammatic expansions.
Geometric Conjecture for Holographic States: Drawing intuition from fixed-area states and random tensor networks, the authors propose a geometric formula for Wigner negativity in general holographic states satisfying the Ryu-Takayanagi formula.

Key Results

PSSY Model (Static/Equilibrium):
- Before the Page Transition: When the radiation system dimension $D$ is small compared to the black hole entropy ( $D \ll e^{S_0}$ ), the Wigner negativity is $O(1)$ . This indicates low stabilizer complexity, implying the radiation state is classically simulable.
- After the Page Transition: When $D \gg e^{S_0}$ , the Wigner negativity becomes exponentially large. The authors derive the asymptotic behavior:
  $N \sim \sqrt{\frac{2}{\pi}} \exp\left[\frac{1}{2}(S_{\text{max}} - S_2)\right]$
  where $S_{\text{max}} = \log D$ is the coarse-grained entropy and $S_2$ is the second Rényi entropy of the radiation.
- Universality: The derived formula is basis-independent and purely information-theoretic. The gravitational path integral naturally averages over the computational basis, yielding a result that depends only on entropic quantities.
Dynamical Model:
- Early Times: Upon coupling the black hole to the bath, the Wigner negativity exhibits a sharp spike (a "negativity shock") due to the initial interaction.
- Late Times: As the system decoheres and settles, the negativity decays and converges to a time-independent value. Crucially, this late-time value precisely matches the universal formula derived from the PSSY gravitational path integral, confirming the consistency between the dynamical random matrix approach and the semi-classical gravity calculation.
Geometric Formula and Python's Lunch:
- The authors propose a geometric formula for Wigner negativity in holographic states:
  $N \sim \sqrt{\frac{2}{\pi}} \exp\left[\frac{1}{2} \frac{A(\gamma_{\text{out}}) - A(\gamma_{\text{min}})}{4G_N}\right]$
  where $\gamma_{\text{out}}$ is the outermost extremal surface and $\gamma_{\text{min}}$ is the minimal extremal surface.
- When the gap between these areas vanishes, negativity is $O(1)$ . When a large gap exists (characteristic of a "python's lunch" region in the entanglement wedge), the negativity is exponentially large.

Significance and Claims
The paper claims to provide a concrete, computable measure of the complexity of Hawking radiation that aligns with the Harlow-Hayden conjecture. By demonstrating that Wigner negativity grows exponentially after the Page time, the authors argue that the state of the Hawking radiation becomes exponentially hard to simulate classically. This computational hardness serves as a protective mechanism for the semi-classical notion of spacetime, preventing observers at infinity from easily performing the unitary operations required to decode the black hole interior.

The work establishes a direct link between gravitational path integral calculations (specifically replica wormholes) and quantum information theoretic measures of complexity. The authors emphasize that their results are "universal" and "basis-independent," suggesting that the high complexity of the radiation is an intrinsic feature of the quantum state rather than an artifact of a specific basis choice. Furthermore, the proposed geometric formula suggests that the presence of a "python's lunch" in the entanglement wedge corresponds to an exponentially large stabilizer complexity, offering a new perspective on the computational cost of bulk reconstruction in AdS/CFT.

The Paper: Measuring the "Magic" of Black Hole Smoke

The Journey Through Time: Before and After the Page Transition

The Geometric Connection: The Python's Lunch

The Verdict

More like this