Abstract

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.Electronic supplementary materialThe online version of this article (doi:10.1007/s41114-016-0001-9) contains supplementary material, which is available to authorized users.

Highlights

  • With the commissioning of the second generation of laser interferometric gravitationalwave detectors, and the recent detection of gravitational waves (Abbott 2016), there is considerable interest in gravitational-wave astronomy

  • The key difficulty is that gravitational waves are unambiguously defined only at future null infinity (J +), whereas in practice the domain of numerical simulations is a region of finite extent using a “3+1” foliation of the spacetime

  • The relative difference in amplitude was of 50 % at most, which is probably acceptable given that these formulas are usually employed in complex astrophysical calculations in which the systematic errors coming from the microphysical modelling are often much larger

Read more

Summary

Introduction

With the commissioning of the second generation of laser interferometric gravitationalwave detectors, and the recent detection of gravitational waves (Abbott 2016), there is considerable interest in gravitational-wave astronomy This is a huge field, covering the diverse topics of: detector hardware construction and design; data analysis; astrophysical source modeling; approximate methods for gravitational-wave calculation; and, when the weak field approach is not valid, numerical relativity. The key difficulty is that gravitational waves are unambiguously defined only at future null infinity (J +), whereas in practice the domain of numerical simulations is a region of finite extent using a “3+1” foliation of the spacetime. This is true for most of the numerical codes, but there are notable exceptions. These attempts include the hyperboloidal method (Frauendiener 2004), Cauchy characteristic matching (Winicour 2005), and

Page 4 of 117
Linearized Einstein equations
Page 6 of 117
Page 8 of 117
Making sense of the TT gauge
Page 10 of 117
The quadrupole formula
Page 14 of 117
Extensions of the quadrupole formula
Page 16 of 117
Basic numerical approaches
Page 18 of 117
Page 20 of 117
The ADM formalism
Page 22 of 117
Page 24 of 117
Page 26 of 117
Extracting gravitational waves using ψ4 in practice: the extrapolation method
Page 28 of 117
Page 30 of 117
Gravitational waves in the Cauchy-perturbative approach
Gauge-invariant metric perturbations
Page 34 of 117
Multipolar expansion of metric perturbations
Page 36 of 117
Gauge-invariant odd-parity perturbations
Λ dΩ sin θ i ta2
Gauge-invariant even-parity perturbations
Xsin m θ
Numerical implementations of the Cauchy-perturbative approach
Asymptotic expressions from odd-parity perturbations
Page 44 of 117
Asymptotic expressions from even-parity perturbations
Page 46 of 117
Asymptotic general expressions
Energy and angular momentum losses
Page 48 of 117
A commonly used convention
Page 50 of 117
Implementation summary
Page 54 of 117
Perturbative expansion
Page 56 of 117
Extraction
Perturbative evolution
Page 58 of 117
Reconstruction
Gravitational waves in the characteristic approach
Page 60 of 117
The Einstein equations in Bondi–Sachs coordinates
Page 62 of 117
The Bondi gauge
Page 64 of 117
General gauge
Page 66 of 117
Page 68 of 117
The gravitational-wave strain
Conformal compactification
Page 70 of 117
The news N
Page 72 of 117
“Appendix
The Newman–Penrose quantity ψ40
Linearized case
Page 74 of 117
Numerical implementations of the characteristic approach
Page 78 of 117
Worldtube boundary data
Reconstruction from spectral modes
Transformation to null affine coordinates
Page 82 of 117
Null affine metric
Metric in Bondi–Sachs coordinates
Starting up the null code at the worldtube
Taylor series method
Page 86 of 117
Initial data
Setting worldtube boundary data for the characteristic code
Estimation of gravitational waves
A comparison among different methods
Page 90 of 117
Comparisons of the accuracy of extraction methods
Page 94 of 117
Page 96 of 117
Page 98 of 117
Page 100 of 117
Page 102 of 117
Page 104 of 117
Page 106 of 117
Page 108 of 117
Page 110 of 117
Page 112 of 117
Page 114 of 117
Page 116 of 117
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call