Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userHigh-accuracy comparison of numerical relativity simulations with post-Newtonian expansions
150
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Numerical simulations of 15 orbits of an equal-mass binary black hole system are presented. Gravitational waveforms from these simulations, covering more than 30 cycles and ending about 1.5 cycles before merger, are compared with those from quasi-circular zero-spin post-Newtonian (PN) formulae. The cumulative phase uncertainty of these comparisons is about 0.05 radians, dominated by effects arising from the small residual spins of the black holes and the small residual orbital eccentricity in the simulations. Matching numerical results to PN waveforms early in the run yields excellent agreement (within 0.05 radians) over the first $sim 15$ cycles, thus validating the numerical simulation and establishing a regime where PN theory is accurate. In the last 15 cycles to merger, however, {em generic} time-domain Taylor approximants build up phase differences of several radians. But, apparently by coincidence, one specific post-Newtonian approximant, TaylorT4 at 3.5PN order, agrees much better with the numerical simulations, with accumulated phase differences of less than 0.05 radians over the 30-cycle waveform. Gravitational-wave amplitude comparisons are also done between numerical simulations and post-Newtonian, and the agreement depends on the post-Newtonian order of the amplitude expansion: the amplitude difference is about 6--7% for zeroth order and becomes smaller for increasing order. A newly derived 3.0PN amplitude correction improves agreement significantly ($<1%$ amplitude difference throughout most of the run, increasing to 4% near merger) over the previously known 2.5PN amplitude terms.
rate research
Read More
We study the effectiveness of stationary-phase approximated post-Newtonian waveforms currently used by ground-based gravitational-wave detectors to search for the coalescence of binary black holes by comparing them to an accurate waveform obtained from numerical simulation of an equal-mass non-spinning binary black hole inspiral, merger and ringdown. We perform this study for the Initial- and Advanced-LIGO detectors. We find that overlaps between the templates and signal can be improved by integrating the match filter to higher frequencies than used currently. We propose simple analytic frequency cutoffs for both Initial and Advanced LIGO, which achieve nearly optimal matches, and can easily be extended to unequal-mass, spinning systems. We also find that templates that include terms in the phase evolution up to 3.5 pN order are nearly always better, and rarely significantly worse, than 2.0 pN templates currently in use. For Initial LIGO we recommend a strategy using templates that include a recently introduced pseudo-4.0 pN term in the low-mass ($M leq 35 MSun$) region, and 3.5 pN templates allowing unphysical values of the symmetric reduced mass $eta$ above this. This strategy always achieves overlaps within 0.3% of the optimum, for the data used here. For Advanced LIGO we recommend a strategy using 3.5 pN templates up to $M=12 MSun$, 2.0 pN templates up to $M=21 MSun$, pseudo-4.0 pN templates up to $65 MSun$, and 3.5 pN templates with unphysical $eta$ for higher masses. This strategy always achieves overlaps within 0.7% of the optimum for Advanced LIGO.
Gravitational waveforms which describe the inspiral, merger and ringdown of coalescing binaries are usually constructed by synthesising information from perturbative descriptions, in particular post-Newtonian theory and black-hole perturbation theory, with numerical solutions of the full Einstein equations. In this paper we discuss the glueing of numerical and post-Newtonian waveforms to produce hybrid waveforms which include subdominant spherical harmonics (higher order modes), and focus in particular on the process of consistently aligning the waveforms, which requires a comparison of both descriptions and a discussion of their imprecisions. We restrict to the non-precessing case, and illustrate the process using numerical waveforms of up to mass ratio $q=18$ produced with the BAM code, and publicly available waveforms from the SXS catalogue. The results also suggest new ways of analysing finite radius errors in numerical simulations.
We present analytic computations of gauge invariant quantities for a point mass in a circular orbit around a Schwarzschild black hole, giving results up to 15.5 post-Newtonian order in this paper and up to 21.5 post-Newtonian order in an online repository. Our calculation is based on the functional series method of Mano, Suzuki and Takasugi (MST) and a recent series of results by Bini and Damour. We develop an optimised method for generating post-Newtonian expansions of the MST series, enabling significantly faster computations. We also clarify the structure of the expansions for large values of $ell$, and in doing so develop an efficient new method for generating the MST renormalised angular momentum, $ u$.
We present an analytic computation of Detweilers redshift invariant for a point mass in a circular orbit around a Kerr black hole, giving results up to 8.5 post-Newtonian order while making no assumptions on the magnitude of the spin of the black hole. Our calculation is based on the functional series method of Mano, Suzuki and Takasugi, and employs a rigorous mode-sum regularization prescription based on the Detweiler-Whiting singular-regular decomposition. The approximations used in our approach are minimal; we use the standard self-force expansion to linear order in the mass ratio, and the standard post-Newtonian expansion in the separation of the binary. A key advantage of this approach is that it produces expressions that include contributions at all orders in the spin of the Kerr black hole. While this work applies the method to the specific case of Detweilers redshift invariant, it can be readily extended to other gauge invariant quantities and to higher post-Newtonian orders.
The recent detection of gravitational waves and electromagnetic counterparts emitted during and after the collision of two neutron stars marks a breakthrough in the field of multi-messenger astronomy. Numerical relativity simulations are the only tool to describe the binarys merger dynamics in the regime when speeds are largest and gravity is strongest. In this work we report state-of-the-art binary neutron star simulations for irrotational (non-spinning) and spinning configurations. The main use of these simulations is to model the gravitational-wave signal. Key numerical requirements are the understanding of the convergence properties of the numerical data and a detailed error budget. The simulations have been performed on different HPC clusters, they use multiple grid resolutions, and are based on eccentricity reduced quasi-circular initial data. We obtain convergent waveforms with phase errors of 0.5-1.5 rad accumulated over approximately 12 orbits to merger. The waveforms have been used for the construction of a phenomenological waveform model which has been applied for the analysis of the recent binary neutron star detection. Additionally, we show that the data can also be used to test other state-of-the-art semi-analytical waveform models.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
