Nigel T. Bishop

Characteristic extraction in numerical relativity: binary black hole merger waveforms at null infinity

The accurate modeling of gravitational radiation is a key issue for gravitational wave astronomy. As simulation codes reach higher accuracy, systematic errors inherent in current numerical relativity wave extraction methods become evident, and may lead to a wrong astrophysical interpretation of the data. In this paper, we give a detailed description of the Cauchy-characteristic extraction technique applied to binary black hole inspiral and merger evolutions to obtain gravitational waveforms that are defined unambiguously, that is, at future null infinity. By this method, we remove finite-radius approximations and the need to extrapolate data from the near zone. Further, we demonstrate that the method is free of gauge effects and thus is affected only by numerical error. Various consistency checks reveal that energy and angular momentum are conserved to high precision and agree very well with extrapolated data. In addition, we revisit the computation of the gravitational recoil and find that finite-radius extrapolation very well approximates the result at J^+. However, the (non-convergent) systematic differences in the extrapolated data are of the same order of magnitude as the (convergent) discretization error of the Cauchy evolution, thus highlighting the need for correct wave extraction.

Strategies for the characteristic extraction of gravitational waveforms

We develop, test, and compare new numerical and geometrical methods for improving the accuracy of extracting waveforms using characteristic evolution. The new numerical method involves use of circular boundaries to the stereographic grid patches which cover the spherical cross sections of the outgoing null cones. We show how an angular version of numerical dissipation can be introduced into the characteristic code to damp the high frequency error arising form the irregular way the circular patch boundary cuts through the grid. The new geometric method involves use of the Weyl tensor component Psi4 to extract the waveform as opposed to the original approach via the Bondi news function. We develop the necessary analytic and computational formula to compute the O(1/r) radiative part of Psi4 in terms of a conformally compactified treatment of null infinity. These methods are compared and calibrated in test problems based upon linearized waves.

Extraction of gravitational waves in numerical relativity

A numerical-relativity calculation yields in general a solution of the Einstein equations including also a radiative part, which is in practice computed in a region of finite extent. Since gravitational radiation is properly defined only at null infinity and in an appropriate coordinate system, the accurate estimation of the emitted gravitational waves represents an old and non-trivial problem in numerical relativity. A number of methods have been developed over the years to “extract” the radiative part of the solution from a numerical simulation and these include: quadrupole formulas, gauge-invariant metric perturbations, Weyl scalars, and characteristic extraction. We review and discuss each method, in terms of both its theoretical background as well as its implementation. Finally, we provide a brief comparison of the various methods in terms of their inherent advantages and disadvantages.

General relativistic null-cone evolutions with a high-order scheme

We present a high-order scheme for solving the full non-linear Einstein equations on characteristic null hypersurfaces using the framework established by Bondi and Sachs. This formalism allows asymptotically flat spaces to be represented on a finite, compactified grid, and is thus ideal for far-field studies of gravitational radiation. We have designed an algorithm based on 4th-order radial integration and finite differencing, and a spectral representation of angular components. Consequently the scheme offers more accuracy at a given computational cost compared to previous methods which are second-order accurate. Based on a newly implemented code, we show that the new numerical scheme remains stable and is convergent at the expected order of accuracy.

Towards the geometry of the Universe from data

We present a new algorithm that can reconstruct the full distributions of metric components within the class of spherically symmetric dust universes that may include a cosmological constant. The algorithm is capable of confronting this class of solutions with arbitrary data and opens a new observational window to determine the value of the cosmological constant. In this work we use luminosity and age data to constrain the geometry of the universe up to a redshift of

Improving fold resistance prediction of HIV-1 against protease and reverse transcriptase inhibitors using artificial neural networks

z = 1.75

The gravitational wave strain in the characteristic formalism of numerical relativity

. We show that, although current data are perfectly compatible with homogeneous models of the universe, simple radially inhomogeneous void models that are sometimes used as alternative explanations for the apparent acceleration of the late time universe cannot yet be ruled out. In doing so we reconstruct the density of cold dark matter out to

Characterizing drug resistance using geometric ensembles from HIV protease dynamics

z = 1.75

Numerically reconstructing the geometry of the Universe from data

and derive constraints on the metric components when the universe was 10.5 Gyr old within a comoving volume of approximately 1 Gpc

Extraction of GWs from a Numerical Simulation

^{3}