We compare numerical evolutions performed with the BSSNOK formulation and a conformal decomposition of a Z4-like formulation of General Relativity. The important difference between the two formulations is that the Z4 formulation has a propagating Hamiltonian constraint, whereas BSSNOK has a zero-speed characteristic variable in the constraint subsystem. In spherical symmetry we evolve both puncture and neutron star initial data. We demonstrate that the propagating nature of the Z4 constraints leads to results that compare favorably with BSSNOK evolutions, especially when matter is present in the spacetime. From the point of view of implementation the new system is a simple modification of BSSNOK.

Numerical relativity simulations of compact binaries with the Z4c and BSSNOK formulations are compared. The Z4c formulation is advantageous in every case considered. In simulations of non-vacuum spacetimes the constraint violations due to truncation errors are between one and three orders of magnitude lower in the Z4c evolutions. Improvements are also found in the accuracy of the computed gravitational radiation. For equal-mass irrotational binary neutron star evolutions we find that the absolute errors in phase and amplitude of the waveforms can be up to a factor of four smaller. The quality of the Z4c numerical data is also demonstrated by a remarkably accurate computation of the ADM mass from surface integrals. For equal-mass non-spinning binary puncture black hole evolutions we find that the absolute errors in phase and amplitude of the waveforms can be up to a factor of two smaller. In the same evolutions we find that away from the punctures the Hamiltonian constraint violation is reduced by between one and two orders of magnitude. Furthermore, the utility of gravitational radiation controlling, constraint preserving boundary conditions for the Z4c formulation is demonstrated. The evolution of spacetimes containing a single compact object confirm earlier results in spherical symmetry. The boundary conditions avoid spurious and non-convergent effects present in high resolution runs with either formulation with a more naive boundary treatment. We conclude that Z4c is preferable to BSSNOK for the numerical solution of the 3+1 Einstein equations with the puncture gauge

We derive an initial value formulation for dynamical Chem-Simons gravity, a modification of general relativity involving parity-violating higher derivative terms. We investigate the structure of the resulting system of partial differential equations thinking about linearization around arbitrary backgrounds. This type of consideration is necessary if we are to establish well-posedness of the Cauchy problem. Treating the field equations as an effective field theory we find that weak necessary conditions for hyperbolicity are satisfied. For the full field equations we find that there are states from which subsequent evolution is not determined. Generically the evolution system closes, but is not hyperbolic in any sense that requires a first order pseudodifferential reduction. In a cursory mode analysis we find that the equations of motion contain terms that may cause ill-posedness of the initial value problem.

We present a new pseudospectral code, bamps, for numerical relativity written with the evolution of collapsing gravitational waves in mind. We employ the first order generalized harmonic gauge formulation. The relevant theory is reviewed and the numerical method is critically examined and specialized for the task at hand. In particular we investigate formulation parameters, gauge and constraint preserving boundary conditions well-suited to non-vanishing gauge source functions. Different types of axisymmetric twist-free moment of time symmetry gravitational wave initial data are discussed. A treatment of the axisymmetric apparent horizon condition is presented with careful attention to regularity on axis. Our apparent horizon finder is then evaluated in a number of test cases. Moving on to evolutions, we investigate modifications to the generalized harmonic gauge constraint damping scheme to improve conservation in the strong field regime. We demonstrate strong-scaling of our pseudospectral penalty code. We employ the Cartoon method to efficiently evolve axisymmetric data in our 3+1 dimensional code. We perform test evolutions of Schwarzschild perturbed by gravitational waves and by gauge pulses, both to demonstrate the use of our blackhole excision scheme and for comparison with earlier results. Finally numerical evolutions of supercritical Brill waves are presented to demonstrate durability of the excision scheme for the dynamical formation of a blackhole.

These lecture notes accompany two classes given at the NRHEP2 school. In the first lecture I introduce the basic concepts used for analyzing well-posedness, that is the existence of a unique solution depending continuously on given data, of evolution partial differential equations. I show how strong hyperbolicity guarantees well-posedness of the initial value problem. Symmetric hyperbolic systems are shown to render the initial boundary value problem well-posed with maximally dissipative boundary conditions. I discuss the Laplace-Fourier method for analyzing the initial boundary value problem. Finally I state how these notions extend to systems that are first order in time and second order in space. In the second lecture I discuss the effect that the gauge freedom of electromagnetism has on the PDE status of the initial value problem. I focus on gauge choices, strong-hyperbolicity and the construction of constraint preserving boundary conditions. I show that strongly hyperbolic pure gauges can be used to build strongly hyperbolic formulations. I examine which of these formulations is additionally symmetric hyperbolic and finally demonstrate that the system can be made boundary stable.