Spherical evolution of the generalized harmonic gauge formulation of general relativity on compactified hyperboloidal slices

We report on the successful numerical evolution of the compactified hyperboloidal initial value problem in general relativity using generalized harmonic gauge. We work in spherical symmetry, using a massless scalar field to drive dynamics. Our treatment is based on the dual-foliation approach, proceeding either by using a height function or by solving the eikonal equation to map between frames. Both are tested here with a naive implementation and with hyperboloidal layers. We present a broad suite of numerical evolutions, including pure gauge perturbations, constraint violating, and satisfying data with and without scalar field matter. We present calculations of spacetimes with a regular center. For black hole spacetimes, we use excision to remove part of the black hole interior. We demonstrate both pointwise and norm convergence at the expected rate of our discretization. We present evolutions in which the scalar field collapses to form a black hole. Evolving nonlinear scalar field perturbations of the Schwarzschild spacetime, we recover the expected quasinormal frequencies and tail decay rates from linear theory.