Post-Newtonian Templates for Gravitational Waves from Compact Binary Inspirals

The advent of gravitational wave (GW) astronomy has provided a new window through which to view and understand the Universe. To fully exploit the potential of GW astronomy, an understanding of all the potential electromagnetic counterparts to a gravitational wave detected source will help maximise the science returns. Here I present a study of the electromagnetic emission from relativistic jets that accompany the merger of binary neutron stars or black hole-neutron star systems. These counterparts provide a probe for the structure and dynamics of these relativistic outflows. Binary neutron star, or neutron star-black hole, mergers are thought to be the dominant progenitor of short gamma-ray bursts (GRBs). Here we investigate the possibility that there is a hidden population of low-Lorentz factor jets resulting in failed GRBs, on-axis orphan afterglows, and what kind of counterparts can be expected given a merger-jet population dominated by these failed-GRB jets. I find that for GW detected mergers, ∼ 80% of the population of on-axis events may result in a failed GRB afterglow. The afterglow of a failed GRB is characterised by the lack of any prompt emission; where the γ-rays are emitted within an optically thick region of the low-Lorentz factor (Γ) outflow and significant suppression via pair production and a high opacity results in the photons coupled to the pair plasma. This plasma will undergo adiabatic expansion, and the photons will decouple at the photospheric radius. The energy in the prompt photons, for a sufficiently low-Γ outflow, will have been significantly suppressed. GW detected mergers have a Malmquist bias towards on-axis events (i.e. the rotational axis of the system), where the peak of the probability distribution is an inclination ∼ 30⁰. If the jets from these mergers have an intrinsic structure out to wider angles, then the majority of mergers will be accompanied by electromagnetic counterparts from these various jet structures. By making some simple assumptions about the energetic structure of a jet outside of a bright core region, the various temporal features that result from a given jet structure can be predicted. Where the population of merger jets is dominated by a single structure model, I show the expected fractions of optical counterparts brighter than m_AB = 21. On 17 August 2017, the Light Interferometer Gravitational Wave Observatory (LIGO) in collaboration with Virgo detected the merger of a binary neutron star system. Various electromagnetic counterparts were detected: the GRB 170817A by Fermi/GBM and INTEGRAL; an optical, blue to red, macro/kilo-nova from ∼ 1/2 day post merger to ∼ 5 − 10 days; and a brightening radio, and X-ray counterpart from ∼ 10 days. Optical detection of this counterpart at a magnitude ∼ 26 was made at ∼ 100 days post-merger. Analysis of this counterpart is consistent with the afterglow of a Gaussian structured jet viewed at the system inclination, ∼ 18 ± 8⁰. If all short GRB jets have a similar jet structure, then the rates of orphan afterglows in deep drilling blind surveys e.g. the Large Synoptic Survey Telescope (LSST), will be higher than those expected from a homogeneous, or ‘top-hat’ jet, population. The rates for the various jet structures for orphan afterglows from mergers are discussed, showing that for a population of failed GRBs, or an intrinsic Gaussian structure, an excess in the orphan rate may be apparent. Understanding the dynamics and structure for the jets from black-hole systems born at the merger of a compact binary can help give clues as to the nature of jets from black holes on all scales. As an aside, I show empirically that regardless of black hole mass or system phenomenology, the relativistic jets from such systems share a universal scaling for the jet power and emitted γ-ray luminosity. This scaling could be due to the similar efficiencies of various processes, or alternatively, the scaling may be able to give insights into the emission and physical processes that are responsible for high-energy photons from these outflows. GW astronomy offers a probe of the most extreme relativistic outflows in the Universe, GRBs. The predicted electromagnetic counterparts from these outflows, in association with GW detections, provides a way to probe the Lorentz-factor distribution for merger-jets. Additionally, the phenomenological shape of the afterglows, at various inclinations, gives an indication of the intrinsic structure of these jets. An understanding of these dynamical and structural qualities can be used to constrain the parent population, merger rates, and binary evolution models for compact binary systems.

We provide ready-to-use time-domain gravitational waveforms for spinning compact binaries with precession effects through 1.5PN order in amplitude and compute their mode decomposition using spin-weighted -2 spherical harmonics. In the presence of precession, the gravitational-wave modes (l,m) contain harmonics originating from combinations of the orbital frequency and precession frequencies. We find that the gravitational radiation from binary systems with large mass asymmetry and large inclination angle can be distributed among several modes. For example, during the last stages of inspiral, for some maximally spinning configurations, the amplitude of the (2,0) and (2,1) modes can be comparable to the amplitude of the (2,2) mode. If the mass ratio is not too extreme, the l=3 and l=4 modes are generally one or two orders of magnitude smaller than the l = 2 modes. Restricting ourselves to spinning, non-precessing compact binaries, we apply the stationary-phase approximation and derive the frequency-domain gravitational waveforms including spin-orbit and spin(1)- spin(2) effects through 1.5PN and 2PN order respectively in amplitude, and 2.5PN order in phase. Since spin effects in the amplitude through 2PN order affect only the first and second harmonics of the orbital phase, they do not extend the mass reach of gravitational-wave detectors. However, they can interfere with other harmonics and lower or raise the signal-to-noise ratio depending on the spin orientation. These ready-to-use waveforms could be employed in the data-analysis of the spinning, inspiraling binaries as well as in comparison studies at the interface between analytical and numerical relativity.

The combined observation of gravitational and electromagnetic waves from the coalescence of two neutron stars marks the beginning of multi-messenger astronomy with gravitational waves (GWs). The development of accurate gravitational waveform models is a crucial prerequisite to extract information about the properties of the binary system that generated a detected GW signal. In binary neutron star systems (BNS), tidal effects also need to be incorporated in the modeling for an accurate waveform representation. Building on previous work [Phys.Rev.D96 121501], we explore the performance of inspiral-merger waveform models that are obtained by adding a numerical relativity (NR) based approximant for the tidal part of the phasing (NRTidal) to existing models for nonprecessing and precessing binary black hole systems (SEOBNRv4, PhenomD and PhenomPv2), as implemented in the LSC Algorithm Library Suite. The resulting BNS waveforms are compared and contrasted to target waveforms hybridizing NR waveforms, covering the last approx. 10 orbits up to merger and extending through the postmerger phase, with inspiral waveforms calculated from 30Hz obtained with TEOBResumS. The latter is a state-of-the-art effective-one-body waveform model that blends together tidal and spin effects. We probe that the combination of the PN-based self-spin terms and of the NRTidal description is necessary to obtain minimal mismatches (< 0.01) and phase differences (< 1 rad) with respect to the target waveforms. However, we also discuss possible improvements and drawbacks of the NRTidal approximant in its current form, since we find that it tends to overestimate the tidal interaction with respect to the TEOBResumS model during the inspiral.

Using our new numerical-relativity code SACRA, long-term simulations for inspiral and merger of black hole (BH)-neutron star (NS) binaries are performed, focusing particularly on gravitational waveforms. As the initial conditions, BH-NS binaries in a quasiequilibrium state are prepared in a modified version of the moving-puncture approach. The BH is modeled by a nonspinning moving puncture and, for the NS, a polytropic equation of state with $\ensuremath{\Gamma}=2$ and the irrotational velocity field are employed. The mass ratio of the BH to the NS, $Q={M}_{\mathrm{BH}}/{M}_{\mathrm{NS}}$, is chosen in the range between 1.5 and 5. The compactness of the NS, defined by $\mathcal{C}=G{M}_{\mathrm{NS}}/{c}^{2}{R}_{\mathrm{NS}}$, is chosen to be between 0.145 and 0.178. For a large value of $Q$ for which the NS is not tidally disrupted and is simply swallowed by the BH, gravitational waves are characterized by inspiral, merger and ringdown waveforms. In this case, the waveforms are qualitatively the same as that from BH-BH binaries. For a sufficiently small value of $Q\ensuremath{\lesssim}2$, the NS may be tidally disrupted before it is swallowed by the BH. In this case, the amplitude of the merger and ringdown waveforms is very low, and thus, gravitational waves are characterized by the inspiral waveform and subsequent quick damping. The difference in the merger and ringdown waveforms is clearly reflected in the spectrum shape and in the ``cutoff'' frequency above which the spectrum amplitude steeply decreases. When a NS is not tidally disrupted (e.g., for $Q=5$), kick velocity, induced by asymmetric gravitational-wave emission, agrees approximately with that derived for the merger of BH-BH binaries, whereas for the case when the tidal disruption occurs, the kick velocity is significantly suppressed.

Inspiraling compact binaries have been identified as one of the most promising sources of gravitational waves for interferometric detectors. Most of these binaries are expected to have circularized by the time their gravitational waves enter the instrument's frequency band. However, the possibility that some of the binaries might still possess a significant eccentricity is not excluded. We imagine a situation in which eccentric signals are received by the detector but not explicitly searched for in the data analysis, which uses exclusively circular waveforms as matched filters. We ascertain the likelihood that these filters, though not optimal, will nevertheless be successful at capturing the eccentric signals. We do this by computing the loss in signal-to-noise ratio incurred when searching for eccentric signals with those nonoptimal filters. We show that for a binary system of a given total mass, this loss increases with increasing eccentricity. We show also that for a given eccentricity, the loss decreases as the total mass is increased.

Gravitational waves (GWs) from binary neutron stars encode unique information about ultra-dense matter through characterisic signatures associated with a variety of phenomena including tidal effects during the inspiral. The main tidal signature depends predominantly on the equation of state (EoS)-related tidal deformability parameter Λ, but at late times is also characterised by the frequency of the star’s fundamental oscillation mode (f-mode). In General Relativity and for nuclear matter, Λ and the f-modes are related by universal relations which may not hold for alternative theories of gravity or exotic matter. Independently measuring Λ and the f-mode frequency enables tests of gravity and the nature of compact binaries. Here we present directly measured constraints on the f-mode frequencies of the companions of GW170817. We also show that future GW detector networks will measure f-mode frequencies to within tens of Hz, enabling precision GW asteroseismology with binary inspiral signals alone.

Matched-filter searches for gravitational waves from coalescing compact binaries by the LIGO Scientific Collaboration use the FINDCHIRP algorithm: an implementation of the optimal filter with innovations to account for unknown signal parameters and to improve performance on detector data that has nonstationary and non-Gaussian artifacts. We provide details on the FINDCHIRP algorithm as used in the search for subsolar mass binaries, binary neutron stars, neutron star--black hole binaries, and binary black holes.

Waveform models are important to gravitational wave data analysis. People recently pay much attention to the waveform model construction for eccentric binary black hole (BBH) coalescence. Several effective-one-body (EOB) Numerical-Relativity waveform models of eccentric BBH coalescence have been constructed. But none of them can treat orbit eccentricity and spin-precessing simultaneously. The current paper focuses on this problem. The authors previously have constructed waveform model for spin-aligned eccentric BBH coalescence SEOBNRE. Here we extend such waveform model to describe eccentric spin-precessing BBH coalescence. We calculate the 2PN orbital radiation-reaction forces and the instantaneous part of the decomposed waveform for a general spinning precessing BBH system in EOB coordinates. We implement these results based on our previous SEOBNRE waveform model. We have also compared our model waveforms to both SXS and RIT numerical relativity waveforms. We find good consistency between our model and numerical relativity. Based on our new waveform model, we analyze the impact of the non-perpendicular spin contributions on waveform accuracy. We find that the non-perpendicular spin contributions primarily affect the phase of the gravitational waveforms. For the current gravitational wave detectors, this contribution is not significant. The future detectors may be affected by such non-perpendicular spin contributions. More importantly our SEOBNRE waveform model, as the first theoretical waveform model to describe eccentric spin-precessing BBH coalescence, can help people to analyze orbit eccentricity and spin precession simultaneously for gravitational wave detection data.

Gravitational wave (GW) observations of coalescing compact binaries will be unique probes of strong-field, dynamical aspects of relativistic gravity. We present a short review of various schemes proposed in the literature to test general relativity (GR) and alternative theories of gravity using inspiral waveforms. Broadly these schemes may be classified into two types: model dependent and model independent. In the model dependent category, GW observations are compared against a specific waveform model representative of a particular theory or a class of theories such as scalar-tensor theories, dynamical Chern–Simons theory and massive graviton theories. Model independent tests are attempts to write down a parametrized gravitational waveform where the free parameters take different values for different theories and (at least some of) which can be constrained by GW observations. We revisit some of the proposed bounds in the case of downscaled LISA configuration (eLISA) and compare them with the original LISA configuration. We also compare the expected bounds on alternative theories of gravity from ground-based and space-based detectors and find that space-based GW detectors can test GR and other theories of gravity with unprecedented accuracies. We then focus on a recent proposal to use singular value decomposition of the Fisher information matrix to improve the accuracies with which post-Newtonian theory can be tested. We extend those results to the case of space-based detector eLISA and discuss its implications.

Emission of anisotropic gravitational radiation from compact binary system leads to a flux of linear momentum. This results in the recoil of the system. We investigate the rate of loss of Linear momentum flux in the far zone of the source using various mass type and current type multipole moments for inspiralling compact binary mergers in quasi-elliptical orbits at 2.5 Post Newtonian order. We compute the linear momentum flux accurate up to $\mathcal{O}(e_t)$ in harmonic coordinate. A 2.5 Post Newtonian Quasi-Keplarian representation of the parametric solution to the Post Newtonian equation of motion for the compact binary system has been adopted here. We also provide a closed-form expression for the accumulated linear momentum from the remote past through the binary evolution.

We review the formation and evolution of compact binary stars consisting of white dwarfs (WDs), neutron stars (NSs), and black holes (BHs). Binary NSs and BHs are thought to be the primary astrophysical sources of gravitational waves (GWs) within the frequency band of ground-based detectors, while compact binaries of WDs are important sources of GWs at lower frequencies to be covered by space interferometers (LISA). Major uncertainties in the current understanding of properties of NSs and BHs most relevant to the GW studies are discussed, including the treatment of the natal kicks which compact stellar remnants acquire during the core collapse of massive stars and the common envelope phase of binary evolution. We discuss the coalescence rates of binary NSs and BHs and prospects for their detections, the formation and evolution of binary WDs and their observational manifestations. Special attention is given to AM CVn-stars — compact binaries in which the Roche lobe is filled by another WD or a low-mass partially degenerate helium-star, as these stars are thought to be the best LISA verification binary GW sources.

Inspiraling compact-object binary systems are promising gravitational wave sources for ground and space-based detectors. The time-dependent signature of these sources is a well-characterized function of a relatively small number of parameters; thus, the favored analysis technique makes use of matched filtering and maximum likelihood methods. As the parameters that characterize the source model vary, so do the templates against which the detector data are compared in the matched filter. For small variations in the parameters, the filter responses are closely correlated. Current analysis methodology samples a bank of filters whose parameter values are chosen so that the correlation between successive samples from successive filters in the bank is 97%. Correspondingly, the additional information available with each successive template evaluation is, in a real sense, only 3% of that already provided by the nearby templates. The reason for such a dense coverage of parameter space is to minimize the chance that a real signal, near the detection threshold, will be missed by the parameter space sampling. Here we investigate the use of Chebyshev interpolation for reducing the number of templates that must be evaluated to obtain the same analysis sensitivity. Additionally, rather than focus on the ``loss'' of signal-to-noise associated with the finite number of filters in the template bank, we evaluate the receiver operating characteristic (ROC) as a measure of the effectiveness of an analysis technique. The ROC relates the false alarm probability to the false dismissal probability of an analysis, which are the quantities that bear most directly on the effectiveness of an analysis scheme. As a demonstration, we compare the present ``dense sampling'' analysis methodology with the ``interpolation'' methodology using Chebyshev polynomials, restricted to one dimension of the multidimensional analysis problem by plotting the ROC curves. We find that the interpolated search can be arranged to have the same false alarm and false dismissal probabilities as the dense sampling strategy using 25% fewer templates. Generalized to the two-dimensional space used in the computationally limited current analyses, this suggests a factor of 2 increase in computational efficiency; generalized to the full seven-dimensional parameter space that characterizes the signal associated with an eccentric binary system of spinning neutron stars or black holes, it suggests an order of magnitude increase in computational efficiency. Since the computational cost of the analysis is driven almost exclusively by the matched filter evaluations, a reduction in the number of template evaluations translates directly into an increase in computational efficiency; additionally, since the computational cost of the analysis is large, the increased efficiency translates also into an increase in the size of the parameter space that can be analyzed and, thus, the science that can be accomplished with the data.

Gravitational waves (GWs) from astrophysical sources can interact with background electromagnetic fields, giving rise to distinctive and potentially detectable electromagnetic signatures. In this paper, we study such interactions for far-field gravitational radiation using the 1 + 3 approach to relativity. Linearized equations for the electromagnetic field on perturbed Minkowski space are derived and solved analytically. The inverse Gertsenshteĭn conversion of GWs in a static electromagnetic field is rederived, and the resultant electromagnetic radiation is shown to be significant for highly magnetized pulsars in compact binary systems. We also obtain a variety of nonlinear interference effects for interacting gravitational and electromagnetic waves, although wave–wave resonances previously described in the literature are absent when the electric–magnetic self-interaction is taken into account. The fluctuation and amplification of electromagnetic energy flux as the GW strength increases towards the gravitational–electromagnetic frequency ratio is a possible signature of gravitational radiation from extended astrophysical sources.

Pinpointing the progenitors of long duration gamma-ray bursts (LGRBs) remains an extremely important question, although it is now clear that at least a fraction of LGRBs originate in the core collapse of massive stars in type Ic supernovae, the pathways to the production of these stars, and their initial masses, remain uncertain. Rotation is thought to be vital in the creation of LGRBs, and it is likely that black hole creation is also necessary. We suggest that these two constraints can be met if the GRB progenitors are very massive stars (>20 solar masses) and are formed in tight binary systems. Using simple models we compare the predictions of this scenario with observations and find that the location of GRBs on their host galaxies are suggestive of main-sequence masses in excess of 20 solar masses, while 50% of the known compact binary systems may have been sufficiently close to have had the necessary rotation rates for GRB creation. Thus, massive stars in compact binaries are a likely channel for at least some fraction of LGRBs.

Variation of mass supply rate from the companion can be smeared out by viscous processes inside an accretion disk. By the time the flow reaches the inner edge, the variation in X-rays needs not reflect the true variation of the rate at the outer edge. However, if the viscosity fluctuates around a mean value, one would expect the viscous time scale also to spread around a mean value. In HMXBs, the size of the viscous Keplerian disk is smaller & thus such a spread could be lower as compared to the LMXBs. If there is an increasing or decreasing trend in viscosity, the interval between enhanced emission would be modified systematically. In the absence of a full knowledge about the variation of mass supply rates at the outer edge, we study ideal circumstances where modulation must take place exactly in orbital time scales when there is an ellipticity in the orbit. We study a few compact binaries using long term RXTE/ASM(1.5-12 keV) & Swift/BAT(15-50keV) data to look for such effects & to infer what these results can tell us about the viscous processes inside the respective disks. We employ three different methods to seek imprints of periodicity on the X-ray variation & found that in all the cases, the location of the peak in the power density spectra is consistent with the orbital frequencies. Interestingly, in HMXBs the peaks are sharp with high rms values, consistent with a small Keplerian disk in a wind fed system. However, in LMXBs with larger Keplerian disk, the peaks are spread out with much lower rms values. X-ray reflections, or superhump phenomena which may also cause such X-ray modulations would not be affected by the size of the Keplerian disk. Our result confirms different sizes of Keplerian disks in 2 classes of binaries. If the orbital period of a binary system is not known, it may be obtained with reasonable accuracy for HMXBs & with lesser accuracy for LMXBs by our method.