Abstract. We present a spectral algorithm for solving the full nonlinear vacuum Einstein field equations in the Bondi framework. Developed within the Spectral Einstein Code (SpEC), we demonstrate spectral characteristic evolution as a technical precursor to Cauchy Characteristic Extraction (CCE), a rigorous method for obtaining gauge-invariant gravitational waveforms from existing and future astrophysical simulations. We demonstrate the new algorithm's stability, convergence, and agreement with existing evolution methods. We explain how an innovative spectral approach enables a two orders of magnitude improvement in computational efficiency.
What is Characteristic Evolution, and why?As an international network of gravitational wave observatories come online, the race to the first direct detection of gravitational waves is expected to herald the beginning of gravitational wave astronomy. Detectors such as Advanced LIGO, VIRGO, GEO, and KAGRA aim for strain sensitivities approaching 10 −24 [1,2,3,4]. Nevertheless, signal candidates from compact binaries or supernovae will be on the cusp of detectability, with very poor signal to noise ratios requiring matched filtering for detection [5]. Matched filtering requires a comprehensive template bank, the generation of which has been a primary goal of the field of numerical relativity. These templates cover a range of expected astronomical phenomena, and are generated by a variety of numerical codes, including the Spectral Einstein Code (SpEC) [6]. Filling out the template bank requires a balance of numerical relativity and analytic waveforms and effective-one-body[8]), with waveform models often calibrated using numerical results [9,10,11].One technical challenge facing the construction of a large template bank is extraction of gauge invariant waveforms from simulations. In general, large computationally intensive simulations are needed to describe the physics of events such as supernovae and compact object binary inspirals. While a waveform-like signal can readily be extracted from arXiv:1406.7029v2 [gr-qc] 24 Sep 2014 anywhere in the computational domain, waveforms are only rigorously defined at future null infinity. Finite radii waveform approximations are universally contaminated by coordinate system dynamics, or gauge effects, which are poorly understood and nearly impossible to remove or even quantify. Comparison with Cauchy Characteristic Extraction, or CCE, an alternative which applies Characteristic Evolution to enable waveform computation at future null infinity, suggests that extrapolation gauge errors could dominate the global error [12]. Characteristic Evolution has been previously implemented at up to 4 th order radial accuracy [13], while complete extraction has been achieved with finite difference/volume methods up to 2 nd order [14,15]. Here we implement inner boundary extraction and evolution. This must be combined with an appropriate outer boundary algorithm to generate the gauge-invariant news at future null infinity. In a hypothetical complete ...