Gravitational Wave Astronomy Research Articles

Asymmetric matter from a dark first-order phase transition

We introduce a model for matter genesis in which both the baryonic and dark matter asymmetries originate from a first-order phase transition in a dark sector with an $SU(3)\ifmmode\times\else\texttimes\fi{}SU(2)\ifmmode\times\else\texttimes\fi{}U(1)$ gauge group and minimal matter content. In the simplest scenario, we predict that dark matter is a dark antineutron with mass of either ${m}_{\overline{n}}=1.36\text{ }\text{ }\mathrm{GeV}$ or ${m}_{\overline{n}}=1.63\text{ }\text{ }\mathrm{GeV}$. Alternatively, dark matter may be comprised of equal numbers of dark antiprotons and pions. In either scenario, this model is highly discoverable through both dark matter direct detection and dark photon search experiments. The strong dark matter self-interactions may ameliorate small-scale structure problems, while the strongly first-order phase transition may be confirmed at future gravitational wave observatories.

Quiescent and Active Galactic Nuclei as Factories of Merging Compact Objects in the Era of Gravitational Wave Astronomy

Galactic nuclei harbouring a central supermassive black hole (SMBH), possibly surrounded by a dense nuclear cluster (NC), represent extreme environments that house a complex interplay of many physical processes that uniquely affect stellar formation, evolution, and dynamics. The discovery of gravitational waves (GWs) emitted by merging black holes (BHs) and neutron stars (NSs), funnelled a huge amount of work focused on understanding how compact object binaries (COBs) can pair up and merge together. Here, we review from a theoretical standpoint how different mechanisms concur with the formation, evolution, and merger of COBs around quiescent SMBHs and active galactic nuclei (AGNs), summarising the main predictions for current and future (GW) detections and outlining the possible features that can clearly mark a galactic nuclei origin.

Pseudoplane-wave gravitational calibrator for gravitational wave observatories

The precisions of existing gravitational calibrators for gravitational wave observatories are limited by their dependence on the relative position between the calibrators and the observatory's test masses. Here we present a novel geometry consisting of four quadrupole rotors placed at the vertices of a rectangle centered on the test mass. The phases and rotation directions are selected to produce a pseudo plane-wave sinusoidal gravitational acceleration with amplitude of ~ 100 fm/s^2. We show that this acceleration only has minimal dependence on the test mass position relative to the rotor array and can yield 0.15% acceleration amplitude uncertainty while tolerating a 1-cm test mass position uncertainty. The acceleration can be directed precisely along the optical axis of the interferometer arm and applies no torque on the test mass. In addition, the small size of the rotors has significant engineering and safety benefits.

Adaptive charge control for the space inertial sensor

Inertial sensors are key components of gravitational wave observations and Earth geodesy missions. An inertial sensor includes an isolated free-floating test mass (TM) surrounded by capacitive electrodes and a housing frame (EH) to perform the relative-position measurement and control the TM in six degrees of freedom. Owing to galactic cosmic rays and solar energetic particles, many additional accelerations are introduced through the Coulomb interaction between charged TMs and their surrounding conducting surfaces. Thus, the TM charge control is critical in space-based missions. A contact-free and ultraviolet light-based charge management system (CMS) was developed to reduce charge-induced noises acting on the TMs and minimize force disturbances that can perturb measurements or interrupt science tasks. However, the operating environment for space charge control is full of uncertainties and disturbances. Physical parameters in the discharging process are rarely measured and will vary owing to changes in solar activity, temperature, and so on. The unpredictability and variability of these parameters affects the CMS performance in long-term space missions and must be evaluated or eliminated. This paper presents a simplified physical model for the discharge process based on electron exchange between the TM and the opposing EH. Subsequently, a model reference adaptive control (MRAC) is proposed for the CMS with parametric uncertainties to maintain the TM charge below a certain level and improve its robustness. The simulation results show that the MRAC can automatically adjust control parameters to eliminate the effect of the variability of the aforementioned physical parameters, and the control precision can reach 0.1 mV under uncertainties, which is superior to that of a classic proportional–integral–derivative controller. This study demonstrated the effects of adaptive charge control and its potential for actual applications.

Effect of a high-precision semianalytical mass function on the merger rate of primordial black holes in dark matter halos

In this work, we study the effect of a high-precision semianalytical mass function on the merger rate of primordial black holes (PBHs) in dark matter halos. For this purpose, we first explain a theoretical framework for dark matter halo models and introduce relevant quantities such as the halo density profile, the concentration parameter, and a high-precision semianalytical function, namely Del Popolo (DP) mass function. In the following, we calculate the merger rate of PBHs in the framework of ellipsoidal-collapse dark matter halo models while considering the DP mass function, and compare it with our previous study for the Sheth-Tormen (ST) mass function. The results show that by taking the mass of PBHs as ${M}_{\mathrm{PBH}}=30{M}_{\ensuremath{\bigodot}}$, the DP mass function can potentially amplify the merger rate of PBHs. Moreover, we calculate the merger rate of PBHs for the DP mass function as a function of their mass and fraction and compare it with the black hole mergers recorded by the Laser Interferometer Gravitational-Wave Observatory (LIGO)-Virgo detectors during the latest observing run. Our findings show that the merger rate of PBHs will fall within the LIGO-Virgo band if ${f}_{\mathrm{PBH}}\ensuremath{\gtrsim}\mathcal{O}({10}^{\ensuremath{-}1})$. This implies that the DP mass function can be used to strengthen constraints on the fraction of PBHs.

The effects of Galactic model uncertainties on LISA observations of double neutron stars

ABSTRACT Observations of binaries containing pairs of neutron stars using the upcoming space-based gravitational wave observatory, LISA, have the potential to improve our understanding of neutron star physics and binary evolution. In this work, we assess the effect of changing the model of the Milky Way’s kinematics and star formation history on predictions of the population of double neutron stars that will be detected and resolved by LISA. We conclude that the spatial distribution of these binaries is insensitive to the choice of Galactic models, compared to the stochastic variation induced by the small sample size. In particular, the time-consuming computation of the binaries’ Galactic orbits is not necessary. The distributions of eccentricity and gravitational-wave frequency are, however, affected by the choice of star formation history. Binaries with eccentricities e > 0.1, which can be measured by LISA observations, are mostly younger than $100\, {\rm Myr}$. We caution that comparisons between different predictions for LISA observations need to use consistent star formation histories, and that the Galactic star formation history should be taken into account in the analysis of the observations themselves. The lack of strong dependence on Galactic models means that LISA detection of double neutron star binaries may provide a relatively clean probe of massive binary star evolution.

Accumulating Errors in Tests of General Relativity with Gravitational Waves: Overlapping Signals and Inaccurate Waveforms

Observations of gravitational waves (GWs) from compact binary coalescences provide powerful tests of general relativity (GR), but systematic errors in data analysis could lead to incorrect scientific conclusions. This issue is especially serious in the third-generation GW detectors in which the signal-to-noise ratio (S/N) is high and the number of detections is large. In this work, we investigate the impacts of overlapping signals and inaccurate waveform models on tests of GR. We simulate mock catalogs for Einstein Telescope and Cosmic Explorer and perform parametric tests of GR using waveform models with different levels of inaccuracy. We find that the systematic error in non-GR parameter estimates could accumulate toward a false deviation from GR when combining results from multiple events, although a Bayesian model selection analysis may not favor a deviation. Waveform inaccuracies contribute most to the systematic errors, but multiple overlapping signals could magnify the effects of systematics owing to the incorrect removal of signals. We also point out that testing GR using selected “golden binaries” with high S/N is even more vulnerable to false deviations from GR. The problem of error accumulation is universal; we emphasize that it must be addressed to fully exploit the data from third-generation GW detectors and that further investigations, particularly in waveform accuracy, will be essential.

A Systematic View of Ten New Black Hole Spins

The launch of NuSTAR and the increasing number of binary black hole (BBH) mergers detected through gravitational wave observations have exponentially advanced our understanding of BHs. Despite the simplicity owed to being fully described by their mass and angular momentum, BHs have remained mysterious laboratories that probe the most extreme environments in the universe. While significant progress has been made in the recent decade, the distribution of spin in BHs has not yet been understood. In this work, we provide a systematic analysis of all known BHs in X-ray binary systems (XBs) that have previously been observed by NuSTAR, but have not yet had a spin measurement made using the “relativistic reflection” method obtained from those data. By looking at all the available archival NuSTAR data of these sources, we measure 10 new BH spins: IGR J17454-2919 − GRS 1758-258 − MAXI J1727-203 − MAXI J0637-430 − a = 0.97 ± 0.02; Swift J1753.5-0127 − V4641 Sgr − 4U 1543-47 − 4U 1957+11 − H 1743-322 − and MAXI J1820+070 − (all uncertainties are at the 1σ confidence level). We discuss the implications of our measurements on the entire distribution of stellar-mass BH spins in XBs, and we compare them with the spin distribution in BBHs, finding that the two distributions are clearly in disagreement. Additionally, we discuss the implications of this work on our understanding of how the “relativistic reflection” spin measurement technique works, and discuss possible sources of systematic uncertainty that can bias our measurements.

The Missing Link between Black Holes in High-mass X-Ray Binaries and Gravitational-wave Sources: Observational Selection Effects

There are few observed high-mass X-ray binaries (HMXBs) that harbor massive black holes (BHs), and none are likely to result in a binary black hole (BBH) that merges within a Hubble time; however, we know that massive merging BBHs exist from gravitational-wave (GW) observations. We investigate the role that X-ray and GW observational selection effects play in determining the properties of their respective detected binary populations. We find that, as a result of selection effects, detectable HMXBs and detectable BBHs form at different redshifts and metallicities, with detectable HMXBs forming at much lower redshifts and higher metallicities than detectable BBHs. We also find disparities in the mass distributions of these populations, with detectable merging BBH progenitors pulling to higher component masses relative to the full detectable HMXB population. Fewer than 3% of detectable HMXBs host BHs >35M ⊙ in our simulated populations. Furthermore, we find the probability that a detectable HMXB will merge as a BBH system within a Hubble time is ≃0.6%. Thus, it is unsurprising that no currently observed HMXB is predicted to form a merging BBH with high probability.

Gravitational Wave Source Populations: Disentangling an AGN Component

The astrophysical origin of over 90 compact binary mergers discovered by the LIGO and Virgo gravitational wave observatories is an open question. While the unusual mass and spin of some of the discovered objects constrain progenitor scenarios, the observed mergers are consistent with multiple interpretations. A promising approach to solve this question is to consider the observed distributions of binary properties and compare them to expectations from different origin scenarios. Here we describe a new hierarchical population analysis framework to assess the relative contribution of different formation channels simultaneously. For this study we considered binary formation in active galactic nucleus (AGN) disks along with phenomenological models, but the same framework can be extended to other models. We find that high-mass and high-mass-ratio binaries appear more likely to have an AGN origin compared to having the same origin as lower-mass events. Future observations of high-mass black hole mergers could further disentangle the AGN component from other channels.

Testing gravitational wave propagation with multiband detections

Effective field theories (EFT) of dark energy (DE) — built to parameterise the properties of DE in an agnostic manner — are severely constrained by measurements of the propagation speed of gravitational waves (GW). However, GW frequencies probed by ground-based interferometers lie around the typical strong coupling scale of the EFT, and it is likely that the effective description breaks down before even reaching that scale. We discuss how this leaves the possibility that an appropriate ultraviolet completion of DE scenarios, valid at scales beyond an EFT description, can avoid present constraints on the GW speed. Instead, additional constraints in the lower frequency LISA band would be harder to escape, since the energies involved are orders of magnitude lower. By implementing a method based on GW multiband detections, we show indeed that a single joint observation of a GW150914-like event by LISA and a terrestrial interferometer would allow one to constrain the speed of light and gravitons to match to within 10-15. Multiband GW observations can therefore firmly constrain scenarios based on the EFT of DE, in a robust and unambiguous way.

Cover Your Basis: Comprehensive Data-driven Characterization of the Binary Black Hole Population

We introduce the first complete nonparametric model for the astrophysical distribution of the binary black hole (BBH) population. Constructed from basis splines, we use these models to conduct the most comprehensive data-driven investigation of the BBH population to date, simultaneously fitting nonparametric models for the BBH mass ratio, spin magnitude and misalignment, and redshift distributions. With GWTC-3, we report the same features previously recovered with similarly flexible models of the mass distribution, most notably the peaks in merger rates at primary masses of ∼10M ⊙ and ∼35M ⊙. Our model reports a suppressed merger rate at low primary masses and a mass-ratio distribution consistent with a power law. We infer a distribution for primary spin misalignments that peaks away from alignment, supporting conclusions of recent work. We find broad agreement with the previous inferences of the spin magnitude distribution: the majority of BBH spins are small (a < 0.5), the distribution peaks at a ∼ 0.2, and there is mild support for a nonspinning subpopulation, which may be resolved with larger catalogs. With a modulated power law describing the BBH merger rate’s evolution in redshift, we see hints of the rate evolution either flattening or decreasing at z ∼ 0.2–0.5, but the full distribution remains entirely consistent with a monotonically increasing power law. We conclude with a discussion of the astrophysical context of our new findings and how nonparametric methods in gravitational-wave population inference are uniquely poised to complement to the parametric approach as we enter the data-rich era of gravitational-wave astronomy.

Rapid hierarchical inference of neutron star equation of state from multiple gravitational wave observations of binary neutron star coalescences

Bayesian hierarchical inference of phenomenological parametrized neutron star equations of state (EoS) from multiple gravitational wave observations of binary neutron star mergers is of fundamental importance in improving our understanding of neutron star structure, the general properties of matter at supra nuclear densities, and the strong nuclear force. However, such an analysis is computationally costly, as it is unable to reuse single-event EoS agnostic parameter estimation runs that are carried out regardless for generating gravitational wave transient catalogs. With the number of events expected to be observable during the fourth observing run (O4) of LIGO/Virgo/KAGRA, this problem can only be expected to worsen. We develop a novel and robust algorithm for rapid and computationally cheap hierarchical inference of parametrized EoS from gravitational wave data which reuses single-event EoS agnostic parameter estimation samples to significantly reduce computational cost. We efficiently include a priori knowledge of neutron star physics as Bayesian priors on the EoS parameters. The high speed and low computational cost of our method allow for efficient recomputation of EoS inference every time a new binary neutron star event is discovered or whenever new observations and theoretical discoveries change the prior on EoS parameters. We test our method on both real and simulated gravitational wave data to demonstrate its accuracy. We show that our computationally cheap method produces EoS constraints that are completely consistent with existing analysis for real data, the chosen fiducial EoS for simulated data. Armed with our fast analysis scheme, we also study the variability of EoS constraints with binary neutron star properties for sets of simulated events drawn in different signal-to-noise ratio and mass ranges.

Bayesian modeling of scattered light in the LIGO interferometers

Excess noise from scattered light poses a persistent challenge in the analysis of data from gravitational wave detectors such as Laser Interferometer Gravitational-wave Observatory. We integrate a physically motivated model for the behavior of these “glitches” into a standard Bayesian analysis pipeline used in gravitational wave science. This allows for the inference of the free parameters in this model, and subtraction of these models to produce glitch-free versions of the data. We show that this inference is an effective discriminator of the presence of the features of these glitches, even when those features may not be discernible in standard visualizations of the data.

Measuring the primordial curvature perturbations from the scalar induced gravitational waves

The scalar induced gravitational waves are produced from primordial curvature perturbations in the second order of perturbations. We constrain the fractional energy density of scalar induced gravitational waves from gravitational waves observations. If there is no detection of the scalar induced gravitational waves, the fractional energy density of scalar induced gravitational waves is constrained by some upper limits. Depends on these upper limits, we can obtain the constraints on the power spectrum of the primordial curvature perturbations. For a power-law scalar power spectrum, the constraints on the power spectrum are affected by adding the upper limit of scalar induced gravitational waves from Square Kilometer Array (SKA). In the standard model, the mean values of the scalar amplitude and the spectral index shift to lower values when SKA is added to the combination of Cosmic Microwave Background (CMB) and Baryon Acoustic Oscillation (BAO) datasets, namely $\ln(10^{10}A_s)=3.038\pm0.013$ and $n_s=0.9589^{+0.0021}_{-0.0011}$ at $68\%$ confidence level. We also consider the effects of the existing ground-based gravitational-wave detectors, the existing Pulsar Timing Arrays (PTAs) and Five-hundred-meter Aperture Spherical radio Telescope (FAST), while the constraints from CMB+BAO datasets are totally within their upper limits of scalar induced gravitational waves. Furthermore, we characterize the scalar fluctuation spectrum in terms of the spectral index $n_s$ and its first two derivatives. We calculate corresponding power spectrum of scalar induced gravitational waves theoretically and give the constraints on the running of the spectral index and the running of the running of the spectral index.

Relaxing CMB bounds on primordial black holes: The role of ionization fronts

The accretion of matter onto primordial black holes (PBHs) during the dark ages and the subsequent energy injection in the medium should have left imprints on the cosmic microwave background (CMB) anisotropies. Recent works have claimed stringent CMB limits on the PBH abundance, hardly compatible with a PBH interpretation of the gravitational-wave observations of binary BH mergers. By using a more realistic accretion model based on hydrodynamical simulations and conservative assumptions for the emission efficiency, we show that CMB limits on the PBH abundance are up to two orders of magnitude less stringent than previously estimated between $10$ and $10^4$ M$_\odot$. This reopens the possibility that PBHs might explain at the same time (at least a fraction of) the dark matter, some of the LIGO-Virgo-KAGRA binary BH mergers and the existence of super-massive BHs. More generally, we emphasize that PBH accretion can be a rather complex physical process with velocity dependences that are hard to assess, which introduces large uncertainties in accretion-based limits on the PBH abundance.

Gravitational wave constraints on spatial covariant gravities

The direct discovery of gravitational waves (GWs) from the coalescence of compact binary components by the LIGO/Virgo/KAGRA Collaboration provides an unprecedented opportunity for exploring the underlying theory of gravity that drives the coalescence process in the strong and highly dynamical field regime of gravity. In this paper, we consider the observational effects of spatial covariant gravities on the propagation of GWs in the cosmological background and obtain the observational constraints on coupling coefficients in the action of spatial covariant gravities from GW observations. We first decompose the GWs into the left- and right-hand circular polarization modes and derive the effects of the spatial covariant gravities on the propagation equation of GWs. We find that these effects can be divided into three classes: 1) frequency-independent effects on GW speed and friction, 2) parity-violating amplitude and velocity birefringences, and 3) a Lorentz-violating damping rate and dispersion of GWs. With these effects, we calculate the corresponding modified waveform of GWs generated by the coalescence of compact binaries. By comparing these new effects with the publicly available posterior samples or results from various tests of gravities with LIGO/Virgo/KAGRA data in the literature, we derive the observational constraints on coupling coefficients of the spatial covariant gravities. These results represent the most comprehensive constraints on the spatial covariant gravities in the literature.

Data quality up to the third observing run of advanced LIGO: Gravity Spy glitch classifications

Understanding the noise in gravitational-wave detectors is central to detecting and interpreting gravitational-wave signals. Glitches are transient, non-Gaussian noise features that can have a range of environmental and instrumental origins. The Gravity Spy project uses a machine-learning algorithm to classify glitches based upon their time–frequency morphology. The resulting set of classified glitches can be used as input to detector-characterisation investigations of how to mitigate glitches, or data-analysis studies of how to ameliorate the impact of glitches. Here we present the results of the Gravity Spy analysis of data up to the end of the third observing run of advanced laser interferometric gravitational-wave observatory (LIGO). We classify 233981 glitches from LIGO Hanford and 379805 glitches from LIGO Livingston into morphological classes. We find that the distribution of glitches differs between the two LIGO sites. This highlights the potential need for studies of data quality to be individually tailored to each gravitational-wave observatory.

Primordial black holes and gravitational waves from nonminimally coupled supergravity inflation

We study the formation of primordial black holes and the generation of gravitational waves in a class of cosmological models that are direct supersymmetric analogs of the observationally favored nonminimally coupled Higgs inflation model. It is known that this type of model naturally includes multiple scalar fields which may be regarded as the inflaton. For the sake of simplicity we focus on the case where the inflaton field space is two dimensional. We analyze the multifield dynamics and find the region of parameters that gives copious production of primordial black holes that may comprise a significant part of the present dark matter abundance. We also compute the spectrum of the gravitational waves and discuss their detectability by means of future ground-based and space-borne gravitational wave observatories.

(An)isotropy measurement with gravitational wave observations

We constrain the distribution of merging compact binaries across the celestial sphere using the GWTC-3 catalog from the LIGO-Virgo-KAGRA Collaborations' (LVK) third observing run. With 63 confident detections from O3, we constrain the relative variability (standard deviation) of the rate density across the sky to be $\ensuremath{\lesssim}16%$ at 90% confidence assuming the logarithm of the rate density is described by a Gaussian random field with correlation length $\ensuremath{\ge}10\ifmmode^\circ\else\textdegree\fi{}$. This tightens to $\ensuremath{\lesssim}3.5%$ when the correlation length is $\ensuremath{\ge}20\ifmmode^\circ\else\textdegree\fi{}$. While the new data provides the tightest constraints on anisotropies available to date, we do not find overwhelming evidence in favor of isotropy, either. A simple counting experiment favors an isotropic distribution by a factor of ${\mathcal{B}}_{\mathrm{ani}}^{\mathrm{iso}}=3.7$, which is nonetheless an improvement of more than a factor of two compared to analogous analyses based on only the LVK's first and second observing runs.