No Arabic abstract
We initiate the development of a horizon-based initial (or rather final) value formalism to describe the geometry and physics of the near-horizon spacetime: data specified on the horizon and a future ingoing null boundary determine the near-horizon geometry. In this initial paper we restrict our attention to spherically symmetric spacetimes made dynamic by matter fields. We illustrate the formalism by considering a black hole interacting with a) inward-falling, null matter (with no outward flux) and b) a massless scalar field. The inward-falling case can be exactly solved from horizon data. For the more involved case of the scalar field we analytically investigate the near slowly evolving horizon regime and propose a numerical integration for the general case.
We present a set of inner boundary conditions for the numerical construction of dynamical black hole space-times, when employing a 3+1 constrained evolution scheme and an excision technique. These inner boundary conditions are heuristically motivated by the dynamical trapping horizon framework and are enforced in an elliptic subsystem of the full Einstein equation. In the stationary limit they reduce to existing isolated horizon boundary conditions. A characteristic analysis completes the discussion of inner boundary conditions for the radiative modes.
A maximally rotating Kerr black hole is said to be extremal. In this paper we introduce the corresponding restrictions for isolated and dynamical horizons. These reduce to the standard notions for Kerr but in general do not require the horizon to be either stationary or rotationally symmetric. We consider physical implications and applications of these results. In particular we introduce a parameter e which characterizes how close a horizon is to extremality and should be calculable in numerical simulations.
We investigate the initial-boundary value problem for linearized gravitational theory in harmonic coordinates. Rigorous techniques for hyperbolic systems are applied to establish well-posedness for various reductions of the system into a set of six wave equations. The results are used to formulate computational algorithms for Cauchy evolution in a 3-dimensional bounded domain. Numerical codes based upon these algorithms are shown to satisfy tests of robust stability for random constraint violating initial data and random boundary data; and shown to give excellent performance for the evolution of typical physical data. The results are obtained for plane boundaries as well as piecewise cubic spherical boundaries cut out of a Cartesian grid.
We present the first numerically stable nonlinear evolution for the leading-order gravitational effective field theory (Quadratic Gravity) in the spherically-symmetric sector. The formulation relies on (i) harmonic gauge to cast the evolution system into quasi-linear form (ii) the Cartoon method to reduce to spherical symmetry in keeping with harmonic gauge, and (iii) order-reduction to 1st-order (in time) by means of introducing auxiliary variables. Well-posedness of the respective initial-value problem is numerically confirmed by evolving randomly perturbed flat-space and black-hole initial data. Our study serves as a proof-of-principle for the possibility of stable numerical evolution in the presence of higher derivatives.
We report on a numerical investigation of black hole evolution in an Einstein dilaton Gauss-Bonnet (EdGB) gravity theory where the Gauss-Bonnet coupling and scalar (dilaton) field potential are symmetric under a global change in sign of the scalar field (a Z2 symmetry). We find that for sufficiently small Gauss-Bonnet couplings Schwarzschild black holes are stable to radial scalar field perturbations, and are unstable to such perturbations for sufficiently large couplings. For the latter case, we provide numerical evidence that there is a band of coupling parameters and black hole masses where the end states are stable scalarized black hole solutions, in general agreement with the results of Macedo et al (2019 Phys. Rev. D 99 104041). For Gauss-Bonnet couplings larger than those in the stable band, we find that an elliptic region forms outside of the black hole horizon, indicating the theory does not possess a well-posed initial value formulation in that regime.