Compact Binary Systems Research Articles

Extreme mass ratio inspirals triggered by massive black hole binaries: from relativistic dynamics to cosmological rates

ABSTRACT Extreme mass ratio inspirals (EMRIs) are compact binary systems characterized by a mass ratio q = m/M in the range 10−9–10−4 and represent primary gravitational wave (GW) sources for the forthcoming Laser Interferometer Space Antenna (LISA). While their standard formation channel involves relaxation processes deflecting compact objects on very low angular momentum orbits around the central massive black hole, a number of alternative formation channels has been proposed, including binary tidal break-up, migration in accretion discs and secular and chaotic dynamics around a massive black hole binary (MBHB). In this work, we take an extensive closer look at this latter scenario, investigating how EMRIs can be triggered by MBHBs, formed in the aftermath of galaxy mergers. By employing a suite of relativistic three-body simulations, we evaluate the efficiency of EMRI formation for different parameters of the MBHB, assessing the importance of both secular and chaotic dynamics. By modelling the distribution of compact objects in galaxy nuclei, we estimate the resulting EMRI formation rate, finding that EMRI are produced in a sharp burst, with peak rates that are 10–100 times higher than the standard two-body relaxation channel, lasting for 106–108 yr. By coupling our results with an estimate of the cosmic MBHB merger rate, we finally forecast that LISA could observe ${\cal O}(10)$ EMRIs per year formed by this channel.

Parameterised population models of transient non-Gaussian noise in the LIGO gravitational-wave detectors

The two interferometric LIGO gravitational-wave observatories provide the most sensitive data to date to study the gravitational-wave universe. As part of a global network, they have completed their third observing run in which they observed many tens of signals from merging compact binary systems. It has long been known that a limiting factor in identifying transient gravitational-wave signals is the presence of transient non-Gaussian noise, which reduce the ability of astrophysical searches to detect signals confidently. Significant efforts are taken to identify and mitigate this noise at the source, but its presence persists, leading to the need for software solutions. Taking a set of transient noise artefacts categorised by the GravitySpy software during the O3a observing era, we produce parameterised population models of the noise projected into the space of astrophysical model parameters of merging binary systems. We compare the inferred population properties of transient noise artefacts with observed astrophysical systems from the GWTC2.1 catalogue. We find that while the population of astrophysical systems tend to have near equal masses and moderate spins, transient noise artefacts are typically characterised by extreme mass ratios and large spins. This work provides a new method to calculate the consistency of an observed candidate with a given class of noise artefacts. This approach could be used in assessing the consistency of candidates found by astrophysical searches (i.e. determining if they are consistent with a known glitch class). Furthermore, the approach could be incorporated into astrophysical searches directly, potentially improving the reach of the detectors, though only a detailed study would verify this.

Gravitational waves in scalar-tensor theory to one-and-a-half post-Newtonian order

We compute the gravitational waves generated by compact binary systems in a class of massless scalar-tensor (ST) theories to the 1.5 post-Newtonian (1.5PN) order beyond the standard quadrupole radiation in general relativity (GR). Using and adapting to ST theories the multipolar-post-Minkowskian and post-Newtonian formalisms originally defined in GR, we obtain the tail and non-linear memory terms associated with the dipole radiation in ST theory. The multipole moments and GW flux of compact binaries are derived for general orbits including the new 1.5PN contribution, and comparison is made with previous results in the literature. In the case of quasi-circular orbits, we present ready-to-use templates for the data analysis of detectors, and for the first time the scalar GW modes for comparisons with numerical relativity results.

Cosmology application of the astrophysics originated gravitational wave

We review the current status of the astrophysics-originated gravitational wave astronomy and several cosmology-related science topics of the future spaceborne/ground-based gravitational wave observatory. Specifically, we highlight the gravitational lensing and gravitational wave sirens-related science cases. We first introduce several astrophysical gravitational wave sources, including the compact binary system (stellar black hole binary, and neutron star binary), continuous-wave sources (fast-spinning neutron star, and extreme mass ratio binary system), and the galactic white dwarf binary. Second, we review gravitational wave time-delay strong lensing, lensing magnification, wave-optics and the diffraction effect. Third, we mention the theoretical development of the luminosity distance space distortion effect. Finally, we discuss the electromagnetic counterpart, multi-messenger astronomy, and its cosmological application.

Global Kinetic Modeling of the Intrabinary Shock in Spider Pulsars

Spider pulsars are compact binary systems composed of a millisecond pulsar and a low-mass companion. The relativistic magnetically dominated pulsar wind impacts onto the companion, ablating it and slowly consuming its atmosphere. The interaction forms an intrabinary shock, a proposed site of particle acceleration. We perform global fully kinetic particle-in-cell simulations of the intrabinary shock, assuming that the pulsar wind consists of plane-parallel stripes of alternating polarity and that the shock wraps around the companion. We find that particles are efficiently accelerated via shock-driven reconnection. We extract first-principles synchrotron spectra and light curves, which are in good agreement with X-ray observations: (1) the synchrotron spectrum is nearly flat, F ν ∝ const; (2) when the pulsar spin axis is nearly aligned with the orbital angular momentum, the light curve displays two peaks, just before and after the pulsar eclipse (pulsar superior conjunction), separated in phase by ∼0.8 rad; (3) the peak flux exceeds the one at the inferior conjunction by a factor of 10.

The quadrupole moment of compact binaries to the fourth post-Newtonian order: I. Non-locality in time and infra-red divergencies

With the aim of providing high accuracy post-Newtonian (PN) templates for the analysis of gravitational waves generated by compact binary systems, we complete the analytical derivation of the source type mass quadrupole moment of compact binaries (without spins) at the fourth PN order of general relativity. Similarly to the case of the conservative 4PN equations of motion, we show that the quadrupole moment at that order contains a non-local (in time) contribution, arising from the tail-transported interaction entering the conservative part of the dynamics. Furthermore, we investigate the infra-red (IR) divergences of the quadrupole moment. In a previous work, this moment has been computed using a Hadamard partie finie procedure for the IR divergences, but the knowledge of the conservative equations of motion indicates that those divergences have to be dealt with by means of dimensional regularization. This work thus derives the difference between the two regularization schemes, which has to be added on top of the previous result. We show that unphysical IR poles start to appear at the 3PN order, and we determine all of these up to the 4PN order. In particular, the non-local tail term comes in along with a specific pole at the 4PN order. It will be proven in a companion paper that the poles in the source-type quadrupole are cancelled in the physical radiative type quadrupole moment measured at future null infinity.

Detection of gravitational wave mixed polarization with single space-based detectors

General Relativity predicts only two tensor polarization modes for gravitational waves while at most six possible polarization modes are allowed in the general metric theory of gravity. The number of polarization modes is determined by the specific modified theory of gravity. Therefore, the determination of polarization modes can be used to test gravitational theory. We introduce a concrete data analysis pipeline for a space-based detector such as LISA to detect the polarization modes of gravitational waves. This method can be used for monochromatic gravitational waves emitted from any compact binary system with a known sky position and frequency to detect mixtures of tensor and extra polarization modes. We use the source $\mathrm{J}0806.3+1527$ with one year of simulation data as an example to show that this approach is capable of probing pure and mixed polarizations without knowing the exact polarization modes. We also find that the ability of detection of extra polarization depends on the gravitational-wave source location and the amplitude of nontensorial components.

Detection of strong scattering close to the eclipse region of PSR B1957+20

ABSTRACT We present the first measurement of pulse scattering close to the eclipse region of PSR B1957+20, which is in a compact binary system with a low-mass star. We measured pulse scattering time-scales up to 0.2 ms close to the eclipse and showed that it scales with the dispersion measure (DM) excess roughly as τ ∝ ΔDM2. Our observations provide the first evidence of strong scattering due to multipath propagation effects in the eclipsing material. We show that Kolmogorov turbulence in the eclipsing material with an inner scale of ∼100 m and an outer scale of the size of the eclipse region can naturally explain the observation. Our results show that the eclipsing material in such systems can be highly turbulent and suggest that scattering is one of the main eclipsing mechanisms at around 1.4 GHz.

High-accuracy simulations of highly spinning binary neutron star systems

With an increasing number of expected gravitational-wave detections of binary neutron star mergers, it is essential that gravitational-wave models employed for the analysis of observational data are able to describe generic compact binary systems. This includes systems in which the individual neutron stars are millisecond pulsars for which spin effects become essential. In this work, we perform numerical-relativity simulations of binary neutron stars with aligned and anti-aligned spins within a range of dimensionless spins of $\chi \sim [-0.28,0.58]$. The simulations are performed with multiple resolutions, show a clear convergence order and, consequently, can be used to test existing waveform approximants. We find that for very high spins gravitational-wave models that have been employed for the interpretation of GW170817 and GW190425 are not capable of describing our numerical-relativity dataset. We verify through a full parameter estimation study in which clear biases in the estimate of the tidal deformability and effective spin are present. We hope that in preparation of the next gravitational-wave observing run of the Advanced LIGO and Advanced Virgo detectors our new set of numerical-relativity data can be used to support future developments of new gravitational-wave models.

Linear momentum flux from inspiralling compact binaries in quasielliptical orbits at 2.5 post-Newtonian order

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.

Dark photon bursts from compact binary systems and constraints

In this work, we consider the burst signal of the dark photon, the hypothetical vector boson of the $U(1)_B$ or $U(1)_{B-L}$ gauge group, generated by a compact binary star system. The absence of the signal in the laser interferometer puts bounds on the coupling constant $\epsilon$ to the ordinary matter. It turns out that if the dark photon is massless, $\epsilon^2$ is on the order of $10^{-37}-10^{-33}$ at most; in the massive case, the upper bound of $\epsilon^2$ is about $10^{-38}-10^{-31}$ in the mass range from $10^{-19}$ eV to $10^{-11}$ eV. These are the first bounds derived from the interferometer observations independent of the assumption of dark photons being dark matter.

Training strategies for deep learning gravitational-wave searches

Compact binary systems emit gravitational radiation which is potentially detectable by current Earth bound detectors. Extracting these signals from the instruments' background noise is a complex problem and the computational cost of most current searches depends on the complexity of the source model. Deep learning may be capable of finding signals where current algorithms hit computational limits. Here we restrict our analysis to signals from non-spinning binary black holes and systematically test different strategies by which training data is presented to the networks. To assess the impact of the training strategies, we re-analyze the first published networks and directly compare them to an equivalent matched-filter search. We find that the deep learning algorithms can generalize low signal-to-noise ratio (SNR) signals to high SNR ones but not vice versa. As such, it is not beneficial to provide high SNR signals during training, and fastest convergence is achieved when low SNR samples are provided early on. During testing we found that the networks are sometimes unable to recover any signals when a false alarm probability $<10^{-3}$ is required. We resolve this restriction by applying a modification we call unbounded Softmax replacement (USR) after training. With this alteration we find that the machine learning search retains $\geq 91.5\%$ of the sensitivity of the matched-filter search down to a false-alarm rate of 1 per month.

Constraining the orbital eccentricity of inspiralling compact binary systems with Advanced LIGO

The detection of ~50 coalescing compact binaries with the Advanced LIGO and Virgo detectors has allowed us to test general relativity, constrain merger rates, and look for evidence of tidal effects, compact object spins, higher waveform modes, and black hole ringdowns. An effect that has not yet been confidently detected is binary eccentricity, which might be present in a small fraction of binaries formed dynamically. Here we discuss general limits on eccentricity that can, in-principle, be placed on all types of compact object binaries by a detector operating at the design sensitivity of Advanced LIGO. Using a post-Newtonian model for gravitational-wave phasing valid in the small eccentricity regime, we assess the relative measurement error for eccentricity for a variety of spinning and non-spinning binaries. Errors and correlations involving the mass and spin parameters are also investigated. We find that decreasing the low frequency limit of a detector's observational frequency band is one of the key design factors for increasing the odds of measuring binary eccentricity. We also introduce and analytically explore the eccentric chirp mass parameter, which replaces the chirp mass as the key measurable parameter combination in eccentric gravitational waveform models. The eccentric chirp mass parameter explains a degeneracy between the chirp mass and the eccentricity. This degeneracy leads to a bias in the standard chirp mass parameter. We also investigate the systematic parameter bias that arises when eccentric systems are recovered using circular waveform templates. We use both Fisher matrix and Bayesian-inference-based Markov Chain Monte Carlo (MCMC) methods to investigate these parameter estimation issues, and we find good agreement between the two approaches (for both statistical and systematic errors) in the appropriate signal-to-noise ratio regime. (abridged)

Leveraging gravitational-wave memory to distinguish neutron star-black hole binaries from black hole binaries

In the observation of gravitational waves (GWs) from a compact binary coalescence system where the mass of one of the companions is $\leq 5~M_{\odot}$ the nature of the object is ambiguous until the measurements of tidal effects give evidence for the presence of a neutron star (NS) or a low mass black hole (BH). The relevance of tidal effects in a neutron star--black hole (NSBH) binary system depends crucially upon the mass and the spin of the companion BH. These effects become important predominantly when the binary system is of comparable mass and/or has large aligned spins. Depending upon the masses and spins the NS can even get tidally disrupted before the innermost stable circular orbit (ISCO) is reached. The gravitational-wave signatures of various tidal effects are encoded in the phasing of the signal and in the case of tidal disruption an abrupt cutoff of the signal amplitude occurs. In this work we show that tidal effects can also be captured by the nonlinear memory of the GW signal. Although small in amplitude, nonlinear memory is present at low frequency in contrast to the oscillatory GW signal. We introduce nonlinear memory in the NSBH and binary black hole (BBH) waveform models and show how the addition of memory aids in distinguishing NSBH systems from BBH systems for a larger part of the parameter space. We discuss the recently detected events of interest by LIGO-Virgo and provide the future prospects for the third generation detectors where nonlinear memory can play a crucial role in inferring the nature of the coalescence as BBH or NSBH from its GW signal alone.

Non-principal axis rotation in binary asteroid systems and how it weakens the BYORP effect

Using viscoelastic mass/spring model simulations, we explore tidal evolution and migration of compact binary asteroid systems. We find that after the secondary is captured into a spin-synchronous state, non-principal axis rotation in the secondary can be long-lived. The secondary’s long axis can remain approximately aligned along the vector connecting secondary to primary while the secondary rocks back and forth about its long axis. Inward orbital semi-major axis migration can also resonantly excite non-principal axis rotation. By estimating solar radiation forces on triangular surface meshes, we show that the magnitude of the BYORP effect induced torque is sensitive to the secondary’s spin state. Non-principal axis rotation within the 1:1 spin–orbit resonance can reduce the BYORP torque or cause frequent reversals in its direction.

Effective chirp mass in the inspiral frequency evolution of the nonspinning eccentric compact binary

Compact binary systems with black holes are the primary sources of interferometric advanced gravitational wave detectors. Astrophysical models suggest finite possibility of binary systems with appreciable non-zero eccentricity in the dense stellar environment like globular clusters and galactic nuclei. While most of the compact binary measurements have not shown appreciable eccentricity, constrains on the eccentricity have been placed on few detected events. With increasing sensitivity of the gravitational wave detectors, the eccentric binaries are plausible sources in the gravitational wave window. Nevertheless, the challenges in the waveform modelling for high values of eccentricity constrain both the search methodologies of eccentric system as well as parameter estimation of the eccentricity. Waveform independent approaches are being investigated. In this work, we demonstrate that a new effective chirp mass $\mathcal{M}_e$ parameter can be defined which governs the dynamics of the non-spinning eccentric compact binary system. We obtain the phenomenological model of the same for low to moderate eccentricity values. One direct implication is its application in the unmodelled searches and as an another implication, we demonstrate that this model can be used to constrain the eccentricity of the non-spinning eccentric binaries in the advanced detector era in the waveform model-independent way.

YY Hya and its interstellar environment

Context. During a search for previously unknown Galactic emission nebulae, we discovered a faint 36′ diameter Hα emission nebula centered around the periodic variable YY Hya. Although this star has been classified as an RR-Lyr variable, such a classification is inconsistent with a Gaia distance of ≃450 pc. The GALEX image data also show YY Hya as having a strong UV excess, suggesting the existence of a hot and compact binary companion. Aims. We aim to clarify the nature of YY Hya and its nebula. Methods. In addition to our discovery image data, we obtained a 2.°5 × 2.°5 image mosaic of the whole region with CHILESCOPE facilities and time-series spectroscopy at MDM observatory. Also, we used data from various space missions to derive an exact orbital period and a spectral energy distribution. Results. We find that YY Hya is a compact binary system containing a K dwarf star that is strongly irradiated by a hot white dwarf companion. The spectral characteristics of the emission lines that are visible only during the maximum light of the perfectly sinusoidal optical light curve show signatures resembling those of members of the BE UMa variable family. These are post-common-envelope pre-cataclysmic variables. However, the companion star here is more massive than that found in other group members and, thus, the progenitor of the white dwarf must have been a star between 3 and 4 M⊙. The nebula appears to be an ejected common-envelope shell with a mass on the order of one M⊙ and an age of 500 000 yr. This makes it the biggest such shell known thus far. The alignment of neighboring nebulosities some 45′ to the northeast and southwest of YY Hya suggests that the system has had strong bipolar outflows. We also briefly speculate that it might be related to the 1065 BP “guest-star” reported in ancient Chinese records.

Following up the afterglow: strategy for X-ray observation triggered by gravitational wave events

The multi-messenger observation of coalescing compact binary systems promises great scientific treasure. However, synthesising observations from both gravitational wave and electromagnetic channels remains challenging. In the context of the day-to-week long emission from a macronova, the binary neutron star merger GW170817 remains the only event with successful electromagnetic followup. In this manuscript, we explore the possibility of using the early stage X-ray afterglow to search for the electromagnetic counterpart of a gravitational wave event. Two algorithms, the simple and straightforward sequential observation (SO) and the step-wise optimizing local optimization are considered and applied to some simulated events. We consider the WXT from the proposed Einstein Probe as a candidate X-ray telescope, which has a very wide field of view of 3600 deg2. Benefiting from the large field of view and high sensitivity, we find that the SO algorithm not only is easy to implement, but also promises a good chance of actual detection.

Ultralight dark matter: constraints from gravitational waves and other astrophysical observations

The orbital period loss of the compact binary systems is the first indirect evidence of gravitational waves which agrees well with Einstein’s general theory of relativity to a very good accuracy. However, there is less than one percent uncertainty in the measurement of orbital period loss from the general reltivistic prediction. Perihelion precession of planets, Gravitational light bending and Shapiro delay are three other successful tests of general relativity theory. Though there are uncertainties in the measurements of those observations from the general reltivistic predictions as well. To resolve these uncertainties, we assume radiation of ultralight axions and light gauge boson particles of L i – Lj type from those systems which can be a possible candidate of fuzzy dark matter. In this article, we obtain bounds on new physics parameters from those astrophysical observations.

Estimating the Chirp Mass of Eccentric Inspiraling Binary Systems from Time-Frequency Representations of their Gravitational Radiation

Using the simplest yet meaningful Peters-Mathews model describing the orbital damping of a compact binary system under the emission of gravitatonal radiation, we show that the chirp-mass of an eccentric inspiraling binary, and its (Keplerian) orbital eccentricity at some reference time, can be estimated from the time-frequency skeleton of its gravitational wave signal. The estimation algorithm is nicely simple, and is robust against the non-ideal (non Gaussian, non stationary) features of detector noise.