On the stabilizer complexity of Hawking radiation

We study the complexity of Hawking radiation for an evaporating black hole from the perspective of the stabilizer theory of quantum computation. Specifically, we calculate Wigner negativity - a magic monotone which can be interpreted as a measure of the stabilizer complexity, or equivalently, the complexity of classical simulation - in various toy models for evaporating black holes. We first calculate the Wigner negativity of Hawking radiation in the PSSY model directly using the gravitational path integral, and show that the negativity is O(1) before the Page transition, but becomes exponentially large past the Page transition. We also derive a universal, information theoretic formula for the negativity which interpolates between the two extremes. We then study the Wigner negativity of radiation in a dynamical model of black hole evaporation. In this case, the negativity shows a sharp spike at early times resulting from the coupling between the black hole and the radiation system, but at late times when the system settles down, we find that the negativity satisfies the same universal formula as in the PSSY model. Finally, we also propose a geometric formula for Wigner negativity in general holographic states using intuition from fixed area states and random tensor networks, and argue that a python's lunch in the entanglement wedge implies a stabilizer complexity which is exponentially large in frac{1}{8G_N} times the difference between the areas corresponding to the outermost and minimal extremal surfaces.