We study spherically symmetric gravitational collapse in cubic Horndeski theories of gravity. By varying the coupling constants and the initial amplitude of the scalar field, we determine the region in the space of couplings and amplitudes for which it is possible to construct global solutions to the Horndeski theories. Furthermore, we identify the regime of validity of effective field theory as the sub-region for which a certain weak field condition remains small at all times. We evolve the initial data using the CCZ4 formulation of the Einstein equations and horizon penetrating coordinates without assuming spherical symmetry.

We demonstrate the flexibility and utility of the Berger-Rigoutsos Adaptive Mesh Refinement (AMR) algorithm used in the open-source numerical relativity code GRChombo for generating gravitational waveforms from binary black-hole inspirals, and for studying other problems involving non-trivial matter configurations. We show that GRChombo can produce high quality binary black-hole waveforms through a code comparison with the established numerical relativity code \lean. We also discuss some of the technical challenges involved in making use of full AMR (as opposed to, e.g. moving box mesh refinement), including the numerical effects caused by using various refinement criteria when regridding. We suggest several ``rules of thumb'' for when to use different tagging criteria for simulating a variety of physical phenomena. We demonstrate the use of these different criteria through example evolutions of a scalar field theory. Finally, we also review the current status and general capabilities of GRChombo.

No abstract

No abstract