Defining eccentricity for spin-precessing binaries

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.

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 report the results of a first study that uses numerical simulations to estimate the accuracy with which one can use gravitational wave observations of double neutron star inspiral to measure parameters of the neutron-star equation of state. The simulations use the evolution and initial-data codes of Shibata and Uryu to compute the last several orbits and the merger of neutron stars, with matter described by a parametrized equation of state. Previous work suggested the use of an effective cutoff frequency to place constraints on the equation of state. We find, however, that greater accuracy is obtained by measuring departures from the point-particle limit of the gravitational waveform produced during the late inspiral. As the stars approach their final plunge and merger, the gravitational wave phase accumulates more rapidly for smaller values of the neutron star compactness (the ratio of the mass of the neutron star to its radius). We estimate that realistic equations of state will lead to gravitational waveforms that are distinguishable from point particle inspirals at an effective distance (the distance to an optimally oriented and located system that would produce an equivalent waveform amplitude) of 100 Mpc or less. As Lattimer and Prakash observed, neutron-star radius is closely tied to the pressure at density not far above nuclear. Our results suggest that broadband gravitational wave observations at frequencies between 500 and 1000 Hz will constrain this pressure, and we estimate the accuracy with which it can be measured. Related first estimates of radius measurability show that the radius can be determined to an accuracy of ~1 km at 100 Mpc.

By reaching through shrouding blastwaves, efficiently discovering off-axis events and probing the central engine at work, gravitational wave (GW) observations will soon revolutionize the study of gamma-ray bursts. Already, analyses of GW data targeting gamma-ray bursts have helped constrain the central engines of selected events. Advanced GW detectors with significantly improved sensitivities are under construction. After outlining the GW emission mechanisms from gamma-ray burst progenitors (binary coalescences, stellar core collapses, magnetars and others) that may be detectable with advanced detectors, we review how GWs will improve our understanding of gamma-ray burst central engines, their astrophysical formation channels and the prospects and methods for different search strategies. We place special emphasis on multimessenger searches. To achieve the most scientific benefit, GW, electromagnetic and neutrino observations should be combined to provide greater discriminating power and science reach.

We investigate the impact of dark fluid accretion on gravitational waveforms emitted by a compact binary system consisting of a supermassive black hole and a stellar-mass black hole. Using a Lagrangian framework with 1 PN and 2.5 PN corrections, we analyze the effects of the spherically symmetric accretion of a fluid with steady-state flow, including those characterized by an equation of state parameter resembling dark energy, on the binary’s dynamics. We validate our approach by comparing it with previous studies in the common region of validity and extend the analysis to include both local effects, such as dynamical friction, and global gravitational interactions with the stellar-mass black hole, focusing on their dependence on the fluid’s properties. Our analysis reveals that these interactions induce de-phasing in gravitational waveforms, with the phase shift influenced by the fluid’s equation of state and energy density. We also extend the study to sudden cosmological singularities, finding that, although they can deform the binary’s orbit from initially circular to elliptical, their effect on de-phasing is negligible for cosmologically relevant energy densities. By incorporating both the local and global gravitational interactions of a fluid on a two-body system into the equations of motion, this preliminary study provides a framework for understanding the interplay between fluid dynamics and gravitational wave emissions in astrophysical systems. It further reinforces the potential for probing the properties of astrophysically relevant fluids through gravitational wave observations.

The detection of primordial black holes (PBHs) would mark a major breakthrough, with far-reaching implications for early universe cosmology, fundamental physics, and the nature of dark matter. Gravitational wave observations have recently emerged as a powerful tool to test the existence and properties of PBHs, as these objects leave distinctive imprints on the gravitational waveform. Notably, there are no known astrophysical processes that can form sub-solar mass black holes, making their discovery a compelling signal of new physics.In addition to PBHs, we consider other exotic compact object (ECO) candidates — such as strange quark stars and boson stars — which can produce similar gravitational signatures and potentially mimic PBHs.In this work, we employ the Fisher matrix formalism to explore a broad parameter space, including binary masses, spins, and different nuclear and quark matter equations of state. Our goal is to assess the ability of next-generation gravitational wave detectors — specifically Cosmic Explorer and the Einstein Telescope — to distinguish PBHs from ECOs, stellar BHs and neutron stars. We compute the maximum luminosity distances (LDs) at which confident (≥ 3σ) detections of sub-solar masses or tidal effects are possible, providing quantitative benchmarks for PBH identification or exclusion under various observational scenarios.Our results indicate that next-generation detectors will be capable of probing sub-solar mass PBHs out to cosmological distances of z ∼ 3. For heavier objects with masses up to M ≲ 2 M ⊙, we show that PBHs can be distinguished from neutron stars via their lack of tidal effects up to redshifts of z ∼ 0.2.

Gravitational wave (GW) observations provide insight into the gravity regime. These observations relies on the comparison of the data to the GWs model, which depends on the orbital evolution of the binary systems. In this paper, we study the orbital evolution of eccentric compact binaries within screened modified gravity, which is a kind of scalar-tensor theory with screening mechanisms. According to the Noether's theorem, from the theory's GW action, we compute the angular momentum flux both as a function of the fields in the theory and as a function of the multipole moments of a N-body system. We specialize to the orbital parameters of eccentric compact binaries to express its orbital energy and angular momentum decay, and then derive the decay rates of its orbital eccentricity and semimajor axis. We find that these quantities decay faster than in General Relativity due to the presence of dipole radiation.

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.

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.

To enable detection and maximise the physics output of gravitational wave observations from compact binary systems, it is crucial the availability of accurate waveform models. The present work aims at giving an overview for non-experts of the (inspiral) waveforms used in the gravitational wave data analysis for compact binary coalescence. We first provide the essential elements of gravitational radiation physics within a simple Newtonian orbital dynamics and the linearized gravity theory, describing the adiabatic approximation applied to binary systems: the key element to construct the theoretical gravitational waveforms in practice. We next lay out the gravitational waveforms in the post-Newtonian approximation to General Relativity, and highlight the basic input for the inspiral waveform of the slowly evolving, spinning, nonprecessing, quasicircular binary black holes, namely, post-Newtonian energy, fluxes and the (absorption-corrected) balance equation. The post-Newtonian inspiral templates are then presented both in the time and frequency domain. Finally, including the merger and subsequent ringdown phase, we briefly survey the two families of the full waveform models of compact binary mergers currently implemented in LSC Algorithm Library Simulation: the effective-one-body approach and the phenomenological frequency domain model.

Bayesian inference is the workhorse of gravitational-wave astronomy, for example, determining the mass and spins of merging black holes, revealing the neutron star equation of state, and unveiling the population properties of compact binaries. The science enabled by these inferences comes with a computational cost that can limit the questions we are able to answer. This cost is expected to grow. As detectors improve, the detection rate will go up, allowing less time to analyze each event. Improvement in low-frequency sensitivity will yield longer signals, increasing the number of computations per event. The growing number of entries in the transient catalog will drive up the cost of population studies. While Bayesian inference calculations are not entirely parallelizable, key components are embarrassingly parallel: calculating the gravitational waveform and evaluating the likelihood function. Graphical processor units (GPUs) are adept at such parallel calculations. We report on progress porting gravitational-wave inference calculations to GPUs. Using a single code - which takes advantage of GPU architecture if it is available - we compare computation times using modern GPUs (NVIDIA P100) and CPUs (Intel Gold 6140). We demonstrate speed-ups of $\sim 50 \times$ for compact binary coalescence gravitational waveform generation and likelihood evaluation and more than $100\times$ for population inference within the lifetime of current detectors. Further improvement is likely with continued development. Our python-based code is publicly available and can be used without familiarity with the parallel computing platform, CUDA.

The observation of gravitational waves (GWs) from merging black holes and neutron stars provides a unique opportunity to discern information about their astrophysical environment. Two signatures that are considered powerful tracers to distinguish between different binary formation channels are general relativistic spin-induced orbital precession and orbital eccentricity. Both effects leave characteristic imprints in the GW signal that can be extracted from observations. To date, neither precession nor eccentricity has been confidently discerned in merging neutron star–black hole binaries. Here we report the measurement of orbital eccentricity in a neutron star–black hole merger. Using, for the first time, a waveform model that incorporates precession and eccentricity, we perform Bayesian inference on the GW event GW200105 and infer a median orbital eccentricity of e 20 ∼ 0.145 at an orbital period of 0.1 s, ruling out eccentricities smaller than 0.028 with 99.5% confidence. We find inconclusive evidence for the presence of precession, consistent with previous, noneccentric results, but a more unequal mass ratio. Our result implies a fraction of these binaries will exhibit orbital eccentricity even at small separations, suggesting formation through mechanisms involving dynamical interactions beyond isolated binary evolution. Future observations will reveal the contribution of eccentric neutron star–black hole binaries to the total merger rate across cosmic time.

Determining the masses of neutron stars (NSs) accurately improves our understanding of the NS interior and complicated binary evolution. However, the masses of the systems are degenerate with the orbital inclination angle when using solely gravitational waves (GWs) or electromagnetic measurements, especially for face-on binaries. Taking advantage of both GWs and optical observations for LISA NS–white dwarf (WD) binaries, we propose a mass determination method utilising multimessenger observational information. By combining the binary mass function obtained from optical observations and a GW mass function, which we introduce, derived from GW observations, we demonstrate how we can set improved constraints on the NS mass and break the degeneracy in the mass and viewing inclination determination. We further comment on the universal relation of the error bar of the GW mass function versus the GW signal-to-noise ratio and propose a simple method for estimating the capability of using GW observations for mass determination with LISA. We show that for ultracompact NS–WD binaries within our Galaxy, the mass of the NS can be constrained to within an accuracy of ±0.2 M ⊙ with the proposed method.

Gravitational-wave (GW) observations from neutron stars (NSs) in a binary system provide an excellent scenario to constrain the nuclear parameters. The investigation of Pratten et al. has shown that the ignorance of f-mode dynamical tidal correction in the GW waveform model of the binary NS system can lead to substantial bias in the measurement of NS properties and NS equations of state. In this work, we investigate the bias in the nuclear parameters resulting from the ignorance of dynamical tidal correction. In addition, this work demonstrates the sensitivity of the nuclear parameters and the estimated constraints on nuclear parameters and NS properties from future GW observations. We infer the nuclear parameters from GW observations by describing the NS matter within the relativistic mean field model. For a population of GW events, we notice that the ignorance of dynamical tide predicts a lower median for nucleon effective mass (m*) by ∼6% compared to the scenario when dynamical tidal correction is considered. Whereas, at a 90% credible interval, m* gets constrained up to ∼5% and ∼3% in A+ (the LIGO-Virgo detectors with a sensitivity of the fifth observing run) and Cosmic Explorer, respectively. We also discuss the resulting constraints on all other nuclear parameters, including compressibility, symmetry energy, and slope of symmetry energy, considering an ensemble of GW events. We do not notice any significant impact in analyzing nuclear parameters other than m* due to the ignorance of f-mode dynamical tides.

It has long been argued that better timing precision allowed by satelites like Rossi X-ray Timing Experiments (RXTE) will allow us to measure the orbital eccentricity and the angle of periastron of some of the bright persistent high mass X-ray binaries (HMXBs) and hence a possible measurement of apsidal motion in these system. Measuring the rate of apsidal motion allows one to estimate the apsidal motion constant of the mass losing companion star and hence allows for the direct testing of the stellar structure models for these giant stars present in the HMXBs. In the present paper we use the archival RXTE data of two bright persistent sources, namely Cen X-3 and SMC X-1, to measure the very small orbital eccentricity and the angle of periastron. We find that the small variations in the pulse profiles of these sources rather than the intrinsic timing accuracy provided by RXTE, limit the accuracy with which we can measure arrival time of the pulses from these sources. This influences the accuracy with which one can measure the orbital parameters, especially the very small eccentricity and the angle of periastron in these sources. The observations of SMC X-1 in the year 2000 were taken during the high flux state of the source and we could determine the orbital eccentricity and $\omega$ using this data set.