Paper detail

Free evolution of the hyperboloidal initial value problem in spherical symmetry

The hyperboloidal initial value problem is addressed in the context of Numerical Relativity, motivated by its use of hyperboloidal slices - smooth spacelike slices that reach future null infinity, the "place" in spacetime where radiation is to be extracted. This is beneficial for studying the global properties of isolated systems and unambiguously extracting their gravitational radiation. The present approach implements the Einstein equations as a free evolution, using the BSSN and Z4 formulations (standard in current codes) expressed in terms of a conformally rescaled metric as suggested by Penrose, and with a time-independent conformal factor. The main difficulty is that the resulting system of PDEs includes formally divergent terms at null infinity that require a special treatment. The numerical simulations in this thesis are restricted to spherical symmetry, although the regularization in the radial direction is expected to also apply to the full 3D case up to some extent. A critical ingredient are the gauge conditions, which rely on well-chosen source functions and damping terms and control the treatment of future null infinity by means of the scri-fixing condition. Once the numerical implementation was stabilized, stable numerical evolutions of a massless scalar field coupled to the Einstein equations could be performed with regular and black hole trumpet initial data on a hyperboloidal slice. The signal of the scalar field has been successfully extracted at future null infinity. Small perturbations of regular initial data give stationary data that are stable forever, while larger scalar field perturbations result in the formation of a black hole. Schwarzschild trumpet initial data have been found to slowly drift away from the expected stationary values, but for small perturbations the effect is slow enough to allow the observation of the power-law decay tails of the scalar field.