Regularizing dual-frame generalized harmonic gauge at null infinity

The dual-frame formalism leads to an approach to extend numerical relativity simulations in generalized harmonic gauge (GHG) all the way to null infinity. A major setback is that without care, even simple choices of initial data give rise to logarithmically divergent terms that would result in irregular variables and equations on the compactified domain, which would in turn prevent accurate numerical approximation. It has been shown, however, that a suitable choice of gauge and constraint addition can be used to prevent their appearance. Presently we give a first order symmetric hyperbolic reduction of general relativity in GHG on compactified hyperboloidal slices that exploits this knowledge and eradicates these log-terms at leading orders. Because of their effect on the asymptotic solution space, specific formally singular terms are systematically chosen to remain. Such formally singular terms have been successfully treated numerically in toy models and result in a formulation with the desirable property that unphysical radiation content near infinity is suppressed.