Coexistence Test of Primordial Black Holes and Particle Dark Matter from Diffractive Lensing

If dark matter (DM) consists of primordial black holes (PBHs) and particles simultaneously, PBHs are generically embedded within particle DM halos. Such "dressed PBHs" (dPBHs) are subject to modified constraints compared to PBHs and can contribute to significant DM abundance in the mass range 10^{-1}-10^{2}M_{⊙}. We show that diffractive lensing of chirping gravitational waves from binary mergers can not only discover, but can also identify dPBH lenses and discriminate them from bare PBHs on the event-by-event basis, with potential to definitively establish the coexistence of subdominant PBHs and particle DM.

Resolving the Stellar-Collapse and Hierarchical-Merger Origins of the Coalescing Black Holes

Spin and mass properties provide essential clues in distinguishing the origins of coalescing black holes (BHs). With a dedicated semiparametric population model for the coalescing binary black holes (BBHs), we identify two distinct categories of BHs among the GWTC-3 events, which is favored over the one population scenario by a logarithmic Bayes factor (lnB) of 7.5. One category, with a mass ranging from ∼25M_{⊙} to ∼80M_{⊙}, is distinguished by the high spin magnitudes (∼0.75) and consistent with the hierarchical merger origin. The other category, characterized by low spins, has a sharp mass cutoff at ∼40M_{⊙}, which is natural for the stellar-collapse origin and in particular the pair-instability explosion of massive stars. We infer the local hierarchical merger rate density as 0.46_{-0.24}^{+0.61}  Gpc^{-3} yr^{-1}. Additionally, we find that a fraction of the BBHs has a cosine-spin-tilt-angle distribution concentrated preferentially around 1, and the fully isotropic distribution for spin orientation is disfavored by a lnB of -6.3, suggesting that the isolated field evolution channels are contributing to the total population.

First Constraints on Compact Binary Environments from LIGO-Virgo Data

The LIGO-Virgo analyses of signals from compact binary mergers observed so far have assumed isolated binary systems in a vacuum, neglecting the potential presence of astrophysical environments. We present here the first investigation of environmental effects on each of the events of GWTC-1 and two low-mass events from GWTC-2. We find no evidence for the presence of environmental effects. Most of the events decisively exclude the scenario of dynamical fragmentation of massive stars as their formation channel. GW170817 results in the most stringent upper bound on the medium density (≲21  g/cm^{3}). We find that environmental effects can substantially bias the recovered parameters in the vacuum model, even when these effects are not detectable. We forecast that the Einstein Telescope and B-DECIGO will be able to probe the environmental effects of accretion disks and superradiant boson clouds on compact binaries.

Detection of early-universe gravitational-wave signatures and fundamental physics

Detection of a gravitational-wave signal of non-astrophysical origin would be a landmark discovery, potentially providing a significant clue to some of our most basic, big-picture scientific questions about the Universe. In this white paper, we survey the leading early-Universe mechanisms that may produce a detectable signal-including inflation, phase transitions, topological defects, as well as primordial black holes-and highlight the connections to fundamental physics. We review the complementarity with collider searches for new physics, and multimessenger probes of the large-scale structure of the Universe.

Spectral Sirens: Cosmology from the Full Mass Distribution of Compact Binaries

We explore the use of the mass spectrum of neutron stars and black holes in gravitational-wave compact binary sources as a cosmological probe. These standard siren sources provide direct measurements of luminosity distance. In addition, features in the mass distribution, such as mass gaps or peaks, will redshift and thus provide independent constraints on their redshift distribution. We argue that the entire mass spectrum should be utilized to provide cosmological constraints. For example, we find that the mass spectrum of LIGO-Virgo-KAGRA events introduces at least five independent mass "features": the upper and lower edges of the pair instability supernova (PISN) gap, the upper and lower edges of the neutron star-black hole gap, and the minimum neutron star mass. We find that although the PISN gap dominates the cosmological inference with current detectors (second generation, 2G), as shown in previous work, it is the lower mass gap that will provide the most powerful constraints in the era of Cosmic Explorer and Einstein Telescope (third generation, 3G). By using the full mass distribution, we demonstrate that degeneracies between mass evolution and cosmological evolution can be broken, unless an astrophysical conspiracy shifts all features of the full mass distribution simultaneously following the (nontrivial) Hubble diagram evolution. We find that this self-calibrating "spectral siren" method has the potential to provide precision constraints of both cosmology and the evolution of the mass distribution, with 2G achieving better than 10% precision on H(z) at z≲1 within a year and 3G reaching ≲1% at z≳2 within one month.

Hints of Spin-Orbit Resonances in the Binary Black Hole Population

Binary black hole spin measurements from gravitational wave observations can reveal the binary's evolutionary history. In particular, the spin orientations of the component black holes within the orbital plane, ϕ_{1} and ϕ_{2}, can be used to identify binaries caught in the so-called spin-orbit resonances. In a companion paper, we demonstrate that ϕ_{1} and ϕ_{2} are best measured near the merger of the two black holes. In this work, we use these spin measurements to provide the first constraints on the full six-dimensional spin distribution of merging binary black holes. In particular, we find that there is a preference for Δϕ=ϕ_{1}-ϕ_{2}∼±π in the population, which can be a signature of spin-orbit resonances. We also find a preference for ϕ_{1}∼-π/4 with respect to the line of separation near merger, which has not been predicted for any astrophysical formation channel. However, the strength of these preferences depends on our prior choices, and we are unable to constrain the widths of the ϕ_{1} and Δϕ distributions. Therefore, more observations are necessary to confirm the features we find. Finally, we derive constraints on the distribution of recoil kicks in the population and use this to estimate the fraction of merger remnants retained by globular and nuclear star clusters. We make our spin and kick population constraints publicly available.

Gravitational Bremsstrahlung and Hidden Supersymmetry of Spinning Bodies

The recently established formalism of a worldline quantum field theory, which describes the classical scattering of massive bodies (black holes, neutron stars, or stars) in Einstein gravity, is generalized up to quadratic order in spin, revealing an alternative N=2 supersymmetric description of the symmetries inherent in spinning bodies. The far-field time domain waveform of the gravitational waves produced in such a spinning encounter is computed at leading order in the post-Minkowskian (weak field, but generic velocity) expansion, and exhibits this supersymmetry. From the waveform we extract the leading-order total radiated angular momentum in a generic reference frame, and the total radiated energy in the center-of-mass frame to leading order in a low-velocity approximation.

Lensing Magnification Seen by Gravitational Wave Detectors

In this paper, we studied the gravitational lensing of gravitational wave events. The probability that an observed gravitational wave source has been (de-)amplified by a given amount is a detector-dependent quantity which depends on different ingredients: the lens distribution, the underlying distribution of sources and the detector sensitivity. The main objective of the present work was to introduce a semi-analytic approach to study the distribution of the magnification of a given source population observed with a given detector. The advantage of this approach is that each ingredient can be individually varied and tested. We computed the expected magnification as both a function of redshift and of the observedsource luminosity distance, which is the only quantity one can access via observation in the absence of an electromagnetic counterpart. As a case study, we then focus on the LIGO/Virgo network and on strong lensing (μ>1).

Toward optomechanical parametric instability prediction in ground-based gravitational wave detectors

Increasing laser power is essential to improve the sensitivity of interferometric gravitational wave detectors. However, optomechanical parametric instabilities can set a limit to that power. It is of major importance to understand and characterize the many parameters and effects that influence these instabilities. Here, we model with a high degree of precision the optical and mechanical modes involved in these parametric instabilities, such that our model can become predictive. As an example, we perform simulations for the Advanced Virgo interferometer (O3 configuration). In particular, we compute mechanical mode losses by combining both on-site measurements and finite element analysis with unprecedented levels of detail and accuracy. We also study the influence on optical modes and parametric gains of mirror finite size effects, and mirror deformations due to thermal absorption. We show that these effects play an important role if transverse optical modes of orders higher than four are involved in the instability process.

Mass and Rate of Hierarchical Black Hole Mergers in Young, Globular and Nuclear Star Clusters

Hierarchical mergers are one of the distinctive signatures of binary black hole (BBH) formation through dynamical evolution. Here, we present a fast semi-analytic approach to simulate hierarchical mergers in nuclear star clusters (NSCs), globular clusters (GCs) and young star clusters (YSCs). Hierarchical mergers are more common in NSCs than they are in both GCs and YSCs because of the different escape velocity. The mass distribution of hierarchical BBHs strongly depends on the properties of first-generation BBHs, such as their progenitor’s metallicity. In our fiducial model, we form black holes (BHs) with masses up to ∼103 M⊙ in NSCs and up to ∼102 M⊙ in both GCs and YSCs. When escape velocities in excess of 100 km s−1 are considered, BHs with mass >103 M⊙ are allowed to form in NSCs. Hierarchical mergers lead to the formation of BHs in the pair instability mass gap and intermediate-mass BHs, but only in metal-poor environments. The local BBH merger rate in our models ranges from ∼10 to ∼60 Gpc−3 yr−1; hierarchical BBHs in NSCs account for ∼10−2–0.2 Gpc−3 yr−1, with a strong upper limit of ∼10 Gpc−3 yr−1. When comparing our models with the second gravitational-wave transient catalog, we find that multiple formation channels are favored to reproduce the observed BBH population.

The first 5 years of gravitational-wave astrophysics

Gravitational waves are ripples in spacetime generated by the acceleration of astrophysical objects; a direct consequence of general relativity, they were first directly observed in 2015. Here, I review the first 5 years of gravitational-wave detections. More than 50 gravitational-wave events have been found, emitted by pairs of merging compact objects such as neutron stars and black holes. These signals yield insights into the formation of compact objects and their progenitor stars, enable stringent tests of general relativity, and constrain the behavior of matter at densities higher than that of an atomic nucleus. Mergers that emit both gravitational and electromagnetic waves probe the formation of short gamma-ray bursts and the nucleosynthesis of heavy elements, and they measure the local expansion rate of the Universe.

GW190521 from the Merger of Ultradwarf Galaxies

We present an alternative formation scenario for the gravitational wave event GW190521 that can be explained as the merger of central black holes (BHs) from two ultradwarf galaxies of stellar mass ∼10^{5}-10^{6}  M_{⊙}, which had themselves previously undergone a merger. The GW190521 components' masses of 85_{-14}^{+21}  M_{⊙} and 66_{-18}^{+17}  M_{⊙} challenge standard stellar evolution models, as they fall in the so-called mass gap. We demonstrate that the merger history of ultradwarf galaxies at high redshifts (1≲z≲2) matches well the LIGO-Virgo inferred merger rate for BHs within the mass range of the GW190521 components, resulting in a likely time delay of ≲4  Gyr considering the redshift of this event. We further demonstrate that the predicted timescales are consistent with expectations for central BH mergers, although with large uncertainties due to the lack of high-resolution simulations in low-mass dwarf galaxies. Our findings show that this BH production and merging channel is viable and extremely interesting as a new way to explore galaxies' BH seeds and galaxy formation. We recommend this scenario be investigated in detail with simulations and observations.

Deep Learning for Gravitational-Wave Data Analysis: A Resampling White-Box Approach

In this work, we apply Convolutional Neural Networks (CNNs) to detect gravitational wave (GW) signals of compact binary coalescences, using single-interferometer data from real LIGO detectors. Here, we adopted a resampling white-box approach to advance towards a statistical understanding of uncertainties intrinsic to CNNs in GW data analysis. We used Morlet wavelets to convert strain time series to time-frequency images. Moreover, we only worked with data of non-Gaussian noise and hardware injections, removing freedom to set signal-to-noise ratio (SNR) values in GW templates by hand, in order to reproduce more realistic experimental conditions. After hyperparameter adjustments, we found that resampling through repeated k-fold cross-validation smooths the stochasticity of mini-batch stochastic gradient descent present in accuracy perturbations by a factor of 3.6. CNNs are quite precise to detect noise, 0.952 for H1 data and 0.932 for L1 data; but, not sensitive enough to recall GW signals, 0.858 for H1 data and 0.768 for L1 data-although recall values are dependent on expected SNR. Our predictions are transparently understood by exploring tthe distribution of probabilistic scores outputted by the softmax layer, and they are strengthened by a receiving operating characteristic analysis and a paired-sample t-test to compare with a random classifier.

Prospects for Measuring the Hubble Constant with Neutron-Star-Black-Hole Mergers

Gravitational wave (GW) and electromagnetic (EM) observations of neutron-star-black-hole (NSBH) mergers can provide precise local measurements of the Hubble constant (H_{0}), ideal for resolving the current H_{0} tension. We perform end-to-end analyses of realistic populations of simulated NSBHs, incorporating both GW and EM selection for the first time. We show that NSBHs could achieve unbiased 1.5%-2.4% precision H_{0} estimates by 2030. The achievable precision is strongly affected by the details of spin precession and tidal disruption, highlighting the need for improved modeling of NSBH mergers.

New Spin on LIGO-Virgo Binary Black Holes

Gravitational waves from binary black holes have the potential to yield information on both of the intrinsic parameters that characterize the compact objects: their masses and spins. While the component masses are usually resolvable, the component spins have proven difficult to measure. This limitation stems in great part from our choice to inquire about the spins of the most and least massive objects in each binary, a question that becomes ill defined when the masses are equal. In this Letter, we show that one can ask a different question of the data: what are the spins of the objects with the highest and lowest dimensionless spins in the binary? We show that this can significantly improve estimates of the individual spins, especially for binary systems with comparable masses. When applying this parametrization to the first 13 gravitational-wave events detected by the LIGO-Virgo Collaboration (LVC), we find that the highest-spinning object is constrained to have nonzero spin for most sources and to have significant support at the Kerr limit for GW151226 and GW170729. A joint analysis of all the confident binary black hole detections by the LVC finds that, unlike with the traditional parametrization, the distribution of spin magnitude for the highest-spinning object has negligible support at zero spin. Regardless of the parametrization used, the configuration where all of the spins in the population are aligned with the orbital angular momentum is excluded from the 90% credible interval for the first ten events and from the 99% credible interval for all current confident detections.

Cygnus X-1 contains a 21-solar mass black hole-Implications for massive star winds

The evolution of massive stars is influenced by the mass lost to stellar winds over their lifetimes. These winds limit the masses of the stellar remnants (such as black holes) that the stars ultimately produce. We used radio astrometry to refine the distance to the black hole x-ray binary Cygnus X-1, which we found to be [Formula: see text] kiloparsecs. When combined with archival optical data, this implies a black hole mass of 21.2 ± 2.2 solar masses, which is higher than previous measurements. The formation of such a high-mass black hole in a high-metallicity system (within the Milky Way) constrains wind mass loss from massive stars.

Optical loss study of the cryogenic molecular layer using a folded cavity for future gravitational-wave detectors

In order to increase the number of detectable gravitational-wave sources, future gravitational-wave detectors will operate with cryogenically cooled mirrors. However, recent studies showed that cryogenic mirrors can suffer from the molecular layer formation, which introduces an additional optical loss, and the detector's performance degrades. In order to evaluate the impact of the molecular layer on future cryogenic gravitational-wave detectors, we built a cryogenic folded-cavity setup and developed an ellipsometric measurement method. The optical loss induced by the cryogenic molecular layer shows a large value even at a few nanometer thickness and can deteriorate the performance of the future cryogenic gravitational-wave detectors.

Hierarchical black hole mergers in multiple systems: constrain the formation of GW190412-, GW190814-, and GW190521-like events

The merging black hole (BH) binaries GW190412, GW190814, and GW190521 from the third LIGO/VIRGO observing run exhibit some extraordinary properties, including highly asymmetric masses, significant spin, and component mass in the ‘mass gap’. These features can be explained if one or both components of the binary are the remnants of previous mergers. In this paper, we explore hierarchical mergers in multiple stellar systems, taking into account the natal kick and mass-loss due to the supernova explosion (SN) on each component, as well as the merger kick received by the merger remnant. The binaries that have survived the SNe and kicks generally have too wide orbital separations to merge by themselves, but can merge with the aid of an external companion that gives rise to Lidov–Kozai oscillations. The BH binaries that consist of second-generation BHs can also be assembled in dense star clusters through binary interactions. We characterize the parameter space of these BH binaries by merger fractions in an analytical approach. Combining the distributions of the survived binaries, we further constrain the parameters of the external companion, using the analytically formulated tertiary perturbation strength. We find that to produce the three LIGO/VIRGO O3 events, the external companions must be at least a few hundreds M⊙, and fall in the intermediate-mass BH and supermassive BH range. We suggest that GW190412, GW190814, and GW190521 could all be produced via hierarchical mergers in multiples, likely in a nuclear star cluster, with the final merger induced by a massive BH.

Consistency of Primordial Black Hole Dark Matter with LIGO/Virgo Merger Rates

The LIGO/Virgo Collaboration has by now observed or constrained the gravitational merger rates of different classes of compact objects. We consider the possibility that the bulk of these mergers are primordial black hole (PBH) mergers of PBHs formed during the QCD epoch making up the entirety of the dark matter. Having shown in a companion paper that mergers due to the initial binary population formed in the early Universe are likely negligible, we compute current merger rates in PBH clusters in which the typical PBH resides. We consider two scenarios: (i) the PBH mass function dictated by the QCD equation of state and (ii) the PBH mass function dictated by the existence of a peak in the inflationary perturbation spectrum. In the first scenario, which is essentially parameter-free, we reproduce very well the merger rates for heavy BHs, the merger rate of mass-asymmetric BHs such as GW190814, a recently discovered merger of a 23  M_{⊙} black hole with a 2.6  M_{⊙} object, and can naturally explain why LIGO/Virgo has not yet observed mergers of two light PBHs from the dominant ∼1  M_{⊙} PBH population. In the second scenario, which has some parameter freedom, we match well the observed rate of heavy PBHs but can currently not explain the rate for mass-asymmetric events. In either case, it is tantalizing that in both scenarios PBH merger rates made with a minimum of assumptions match most LIGO/Virgo observed rates very well.

Spin-Induced Scalarized Black Holes

It was recently shown that a scalar field suitably coupled to the Gauss-Bonnet invariant G can undergo a spin-induced linear tachyonic instability near a Kerr black hole. This instability appears only once the dimensionless spin j is sufficiently large, that is, j≳0.5. A tachyonic instability is the hallmark of spontaneous scalarization. Focusing, for illustrative purposes, on a class of theories that do exhibit this instability, we show that stationary, rotating black hole solutions do indeed have scalar hair once the spin-induced instability threshold is exceeded, while black holes that lie below the threshold are described by the Kerr solution. Our results provide strong support for spin-induced black hole scalarization.

Multiparameter Tests of General Relativity Using Multiband Gravitational-Wave Observations

In this Letter, we show that multiband observations of stellar-mass binary black holes by the next generation of ground-based observatories (3G) and the space-based Laser Interferometer Space Antenna (LISA) would facilitate a comprehensive test of general relativity by simultaneously measuring all the post-Newtonian coefficients. Multiband observations would measure most of the known post-Newtonian phasing coefficients to an accuracy below a few percent-2 orders-of-magnitude better than the best bounds achievable from even "golden" binaries in the 3G or LISA bands. Such multiparameter bounds would play a pivotal role in constraining the parameter space of modified theories of gravity beyond general relativity.

Intermediate Mass-Ratio Black Hole Binaries: Applicability of Small Mass-Ratio Perturbation Theory

The inspiral phasing of binary black holes at intermediate mass ratios (m_{2}/m_{1}∼10^{-3}) is important for gravitational wave observations, but not accessible to standard modeling techniques: The accuracy of the small mass-ratio (SMR) expansion is unknown at intermediate mass ratios, whereas numerical relativity simulations cannot reach this regime. This article assesses the accuracy of the SMR expansion by extracting the first three terms of the SMR expansion from numerical relativity data for nonspinning, quasicircular binaries. We recover the leading term predicted by SMR theory and obtain a robust prediction of the next-to-leading term. The influence of higher-order terms is bounded to be small, indicating that the SMR series truncated at next-to-leading order is quite accurate at intermediate mass ratios and even at nearly comparable mass binaries. We estimate the range of applicability for SMR and post-Newtonian series for nonspinning, quasicircular inspirals.

Constraining the Lensing of Binary Black Holes from Their Stochastic Background

Gravitational waves (GWs) are subject to gravitational lensing in the same way as electromagnetic radiation. However, to date, no unequivocal observation of a lensed GW transient has been reported. Independently, GW observatories continue to search for the stochastic GW signal that is produced by many transient events at high redshift. We exploit a surprising connection between the lensing of individual transients and limits to the background radiation produced by the unresolved population of binary back hole mergers: we show that it constrains the fraction of individually resolvable lensed binary black holes to less than ∼4×10^{-5} at present sensitivity. We clarify the interpretation of existing, low redshift GW observations (obtained assuming no lensing) in terms of their apparent lensed redshifts and masses and explore constraints from GW observatories at future sensitivity. Based on our results, recent claims of observations of lensed events are statistically disfavored.

Neutrino and Positron Constraints on Spinning Primordial Black Hole Dark Matter

Primordial black holes can have substantial spin-a fundamental property that has a strong effect on its evaporation rate. We conduct a comprehensive study of the detectability of primordial black holes with non-negligible spin, via the searches for the neutrinos and positrons in the MeV energy range. Diffuse supernova neutrino background searches and observation of the 511 keV gamma-ray line from positrons in the Galactic center set competitive constraints. Spinning primordial black holes are probed up to a slightly higher mass range compared to nonspinning ones. Our constraint using neutrinos is slightly weaker than that due to the diffuse gamma-ray background, but complementary and robust. Our positron constraints are typically weaker in the lower mass range and stronger in the higher mass range for the spinning primordial black holes compared to the nonspinning ones. They are generally stronger than those derived from the diffuse gamma-ray measurements for primordial black holes having masses greater than a few ×10^{16}  g.

Astrophysical Implications of GW190412 as a Remnant of a Previous Black-Hole Merger

Two of the dominant channels to produce merging stellar-mass black-hole binaries are believed to be the isolated evolution of binary stars in the field and dynamical formation in star clusters. The first reported black-hole binary event from the third LIGO/Virgo observing run (GW190412) is unusual in that it has unequal masses, nonzero effective spin, and nonzero primary spin at 90% confidence interval. We show that this event should be exceedingly rare in the context of both the field and cluster formation scenarios. Interpreting GW190412 as a remnant of a previous black-hole merger provides a promising route to explain its features. If GW190412 indeed formed hierarchically, we show that the region of the parameter space that is best motivated from an astrophysical standpoint (low natal spins and light clusters) cannot accommodate the observation. We analyze public GW190412 LIGO/Virgo data with a Bayesian prior where the more massive black hole resulted from a previous merger and find that this interpretation is equally supported by the data. If the heavier component of GW190412 is indeed a merger remnant, then its spin magnitude is χ_{1}=0.56_{-0.21}^{+0.19}, which is higher than the value previously reported by the LIGO/Virgo collaboration.

GW190521: A Binary Black Hole Merger with a Total Mass of 150 M_{⊙}

On May 21, 2019 at 03:02:29 UTC Advanced LIGO and Advanced Virgo observed a short duration gravitational-wave signal, GW190521, with a three-detector network signal-to-noise ratio of 14.7, and an estimated false-alarm rate of 1 in 4900 yr using a search sensitive to generic transients. If GW190521 is from a quasicircular binary inspiral, then the detected signal is consistent with the merger of two black holes with masses of 85_{-14}^{+21}  M_{⊙} and 66_{-18}^{+17}  M_{⊙} (90% credible intervals). We infer that the primary black hole mass lies within the gap produced by (pulsational) pair-instability supernova processes, with only a 0.32% probability of being below 65  M_{⊙}. We calculate the mass of the remnant to be 142_{-16}^{+28}  M_{⊙}, which can be considered an intermediate mass black hole (IMBH). The luminosity distance of the source is 5.3_{-2.6}^{+2.4}  Gpc, corresponding to a redshift of 0.82_{-0.34}^{+0.28}. The inferred rate of mergers similar to GW190521 is 0.13_{-0.11}^{+0.30}  Gpc^{-3} yr^{-1}.