Semiclassical Dynamics of Hawking Radiation*

The problem of solving the equations for Hawking radiation, with full backreaction included, is cast as an initial value problem suitable for future numerical simulations. Quantum entanglement of the fields is captured using bilocal collective fields in a large N limit. Corrections to the Einstein equations contain up to four spacetime derivatives. These corrections are treated systematically using perturbative constraints.