Celestial Matrix Model.

We construct a Hermitian random matrix model that provides a stable nonperturbative completion of Cangemi-Jackiw (CJ) gravity, a two-dimensional theory of flat spacetimes. The matrix model reproduces, to all orders in the topological expansion, the Euclidean partition function of CJ gravity with an arbitrary number of boundaries. The nonperturbative completion enables the exact computation of observables in flat space quantum gravity which we use to explicitly characterize the Bondi Hamiltonian spectrum. We discuss the implications of our results for the flat space S-matrix and black holes.