Hierarchical analysis of the binary black hole (BBH) detections by the Advanced LIGO and Virgo detectors has offered an increasingly clear picture of their mass, spin, and redshift distributions. Fully understanding the formation and evolution of BBH
mergers will require not just the characterization of these marginal distributions, though, but the discovery of any correlations that exist between the properties of BBHs. Here, we hierarchically analyze the ensemble of BBHs discovered by the LIGO and Virgo with a model that allows for intrinsic correlations between their mass ratios $q$ and effective inspiral spins $chi_mathrm{eff}$. At $98.7%$ credibility, we find that the mean of the $chi_mathrm{eff}$ distribution varies as a function of $q$, such that more unequal-mass BBHs exhibit systematically larger $chi_mathrm{eff}$. We find Bayesian odds ratio of $10.5$ in favor of a model that allows for such a correlation over one that does not. Finally, we use simulated signals to verify that our results are robust against degeneracies in the measurements of $q$ and $chi_mathrm{eff}$ for individual events. While many proposed astrophysical formation channels predict some degree correlation between spins and mass ratio, these predicted correlations typically act in an opposite sense to the trend we observationally identify in the data.
Intermediate mass ratio inspiral (IMRI) binaries -- containing stellar-mass black holes coalescing into intermediate-mass black holes ($M>100M_{odot}$) -- are a highly anticipated source of gravitational waves (GWs) for Advanced LIGO/Virgo. Their det
ection and source characterization would provide a unique probe of strong-field gravity and stellar evolution. Due to the asymmetric component masses and the large primary, these systems generically excite subdominant modes while reducing the importance of the dominant quadrupole mode. Including higher order harmonics can also result in a $10%-25%$ increase in signal-to-noise ratio for IMRIs, which may help to detect these systems. We show that by including subdominant GW modes into the analysis we can achieve a precise characterization of IMRI source properties. For example, we find that the source properties for IMRIs can be measured to within $2%-15%$ accuracy at a fiducial signal-to-noise ratio of 25 if subdominant modes are included. When subdominant modes are neglected, the accuracy degrades to $9%-44%$ and significant biases are seen in chirp mass, mass ratio, primary spin and luminosity distances. We further demonstrate that including subdominant modes in the waveform model can enable an informative measurement of both individual spin components and improve the source localization by a factor of $sim$10. We discuss some important astrophysical implications of high-precision source characterization enabled by subdominant modes such as constraining the mass gap and probing formation channels.
Gravitational waves in general relativity contain two polarization degrees of freedom, commonly labeled plus and cross. Besides those two tensor modes, generic theories of gravity predict up to four additional polarization modes: two scalar and two v
ector. Detection of nontensorial modes in gravitational wave data would constitute a clean signature of physics beyond general relativity. Previous measurements have pointed to the unambiguous presence of tensor modes in gravitational waves, but the presence of additional generic nontensorial modes has not been directly tested. We propose a model-independent analysis capable of detecting and characterizing mixed tensor and nontensor components in transient gravitational wave signals, including those from compact binary coalescences. This infrastructure can constrain the presence of scalar or vector polarization modes on top of the tensor modes predicted by general relativity. Our analysis is morphology-independent (as it does not rely on a waveform templates), phase-coherent, and agnostic about the source sky location. We apply our analysis to data from GW190521 and simulated data and demonstrate that it is capable of placing upper limits on the strength of nontensorial modes when none are present, or characterizing their morphology in the case of a positive detection. Tests of the polarization content of a transient gravitational wave signal hinge on an extended detector network, wherein each detector observes a different linear combination of polarization modes. We therefore anticipate that our analysis will yield precise polarization constraints in the coming years, as the current ground-based detectors LIGO Hanford, LIGO Livingston, and Virgo are joined by KAGRA and LIGO India.
The recent observation of GW190412, the first high-mass ratio binary black-hole (BBH) merger, by the LIGO-Virgo Collaboration (LVC) provides a unique opportunity to probe the impact of subdominant harmonics and precession effects encoded in a gravita
tional wave signal. We present refined estimates of source parameters for GW190412 using texttt{NRSur7dq4}, a recently developed numerical relativity waveform surrogate model that includes all $ell leq 4$ spin-weighted spherical harmonic modes as well as the full physical effects of precession. We compare our results with two different variants of phenomenological precessing BBH waveform models, texttt{IMRPhenomPv3HM} and texttt{IMRPhenomXPHM}, as well as to the LVC results. Our results are broadly in agreement with texttt{IMRPhenomXPHM} results and the reported LVC analysis compiled with the texttt{SEOBNRv4PHM} waveform model, but in tension with texttt{IMRPhenomPv3HM}. Using the texttt{NRSur7dq4} model, we provide a tighter constraint on the mass-ratio ($0.26^{+0.08}_{-0.06}$) as compared to the LVC estimate of $0.28^{+0.13}_{-0.07}$ (both reported as median values withs 90% credible intervals). We also constrain the binary to be more face-on, and find a broader posterior for the spin precession parameter. We further find that even though $ell=4$ harmonic modes have negligible signal-to-noise ratio, omission of these modes will influence the estimated posterior distribution of several source parameters including chirp mass, effective inspiral spin, luminosity distance, and inclination. We also find that commonly used model approximations, such as neglecting the asymmetric modes (which are generically excited during precession), have negligible impact on parameter recovery for moderate SNR-events similar to GW190412.
The identification of the electromagnetic counterpart candidate ZTF19abanrhr to the binary black hole merger GW190521 opens the possibility to infer cosmological parameters from this standard siren with a uniquely identified host galaxy. The distant
merger allows for cosmological inference beyond the Hubble constant. Here we show that the three-dimensional spatial location of ZTF19abanrhr calculated from the electromagnetic data remains consistent with the updated sky localization of GW190521 provided by the LIGO-Virgo Collaboration. If ZTF19abanrhr is associated with the GW190521 merger and assuming a flat wCDM model we find that $H_0 =48^{+24}_{-10}$ km/s/Mpc, $Omega_m =0.39^{+0.38}_{-0.29}$, and $w_0 = -1.29^{+0.63}_{-0.50}$ (median and 68% credible interval). If we use the Hubble constant value inferred from another gravitational-wave event, GW170817, as a prior for our analysis, together with assumption of a flat ${Lambda}$CDM and the model-independent constraint on the physical matter density ${omega}_m$ from Planck, we find $H_0 = 69.1^{8.7}_{-6.0}$ km/s/Mpc.
The data from ground based gravitational-wave detectors such as Advanced LIGO and Virgo must be calibrated to convert the digital output of photodetectors into a relative displacement of the test masses in the detectors, producing the quantity of int
erest for inference of astrophysical gravitational wave sources. Both statistical uncertainties and systematic errors are associated with the calibration process, which would in turn affect the analysis of detected sources, if not accounted for. Currently, source characterization algorithms either entirely neglect the possibility of calibration uncertainties or account for them in a way that does not use knowledge of the calibration process itself. We present physiCal, a new approach to account for calibration errors during the source characterization step, which directly uses all the information available about the instrument calibration process. Rather than modeling the overall detectors response function, we consider the individual components that contribute to the response. We implement this method and apply it to the compact binaries detected by LIGO and Virgo during the second observation run, as well as to simulated binary neutron stars for which the sky position and distance are known exactly. We find that the physiCal model performs as well as the method currently used within the LIGO-Virgo collaboration, but additionally it enables improving the measurement of specific components of the instrument control through astrophysical calibration.
Gravitational waves may be one of the few direct observables produced by ultralight bosons, conjectured dark matter candidates that could be the key to several problems in particle theory, high-energy physics and cosmology. These axionlike particles
could spontaneously form clouds around astrophysical black holes, leading to potent emission of continuous gravitational waves that could be detected by instruments on the ground and in space. Although this scenario has been thoroughly studied, it has not been yet appreciated that both types of detector may be used in tandem (a practice known as multibanding). In this paper, we show that future gravitational-wave detectors on the ground and in space will be able to work together to detect ultralight bosons with masses $25 lesssim mu/left(10^{-15}, mathrm{eV}right)lesssim 500$. In detecting binary-black-hole inspirals, the LISA space mission will provide crucial information enabling future ground-based detectors, like Cosmic Explorer or Einstein Telescope, to search for signals from boson clouds around the individual black holes in the observed binaries. We lay out the detection strategy and, focusing on scalar bosons, chart the suitable parameter space. We study the impact of ignorance about the systems history, including cloud age and black hole spin. We also consider the tidal resonances that may destroy the boson cloud before its gravitational signal becomes detectable by a ground-based follow-up. Finally, we show how to take all of these factors into account, together with uncertainties in the LISA measurement, to obtain boson mass constraints from the ground-based observation facilitated by LISA.
Gravitational waves emitted by neutron star black hole mergers encode key properties of neutron stars - such as their size, maximum mass and spins - and black holes. However, the presence of matter and the high mass ratio makes generating long and ac
curate waveforms from these systems hard to do with numerical relativity, and not much is known about systematic uncertainties due to waveform modeling. We simulate gravitational waves from neutron star black hole mergers by hybridizing numerical relativity waveforms produced with the SpEC code with a recent numerical relativity surrogate NRHybSur3dq8Tidal. These signals are analyzed using a range of available waveform families, and statistical and systematic errors are reported. We find that at a network signal-to-noise ratio (SNR) of 30, statistical uncertainties are usually larger than systematic offsets, while at an SNR of 70 the two become comparable. The individual black hole and neutron star masses, as well as the mass ratios, are typically measured very precisely, though not always accurately at high SNR. At a SNR of 30 the neutron star tidal deformability can only be bound from above, while for louder sources it can be measured and constrained away from zero. All neutron stars in our simulations are non-spinning, but in no case we can constrain the neutron star spin to be smaller than $sim0.4$ (90% credible interval). Waveform families whose late inspiral has been tuned specifically for neutron star black hole signals typically yield the most accurate characterization of the source parameters. Their measurements are in tension with those obtained using waveform families tuned against binary neutron stars, even for mass ratios that could be relevant for both binary neutron stars and neutron star black holes mergers.
Our understanding of observed Gravitational Waves (GWs) comes from matching data to known signal models describing General Relativity (GR). These models, expressed in the post-Newtonian formalism, contain the mathematical constant $pi$. Allowing $pi$
to vary thus enables a strong, universal and generalisable null test of GR. From a population of 22 GW observations, we make an astrophysical measurement of $pi=3.115^{+0.048}_{-0.088}$, and prefer GR as the correct theory of gravity with a Bayes factor of 321. We find the variable $pi$ test robust against simulated beyond-GR effects.
Finite-size effects on the gravitational wave signal from a neutron star merger typically manifest at high frequencies where detector sensitivity decreases. Proposed sensitivity improvements can give us access both to stronger signals and to a myriad
of weak signals from cosmological distances. The latter will outnumber the former and the relevant part of signal will be redshifted towards the detector most sensitive band. We study the redshift dependence of information about neutron star matter and find that single-scale properties, such as the star radius or the post-merger frequency, are better measured from the distant weak sources from $zsim 1$.
