Scalar fields around black hole binaries in LIGO-Virgo-KAGRA

This paper presents a validated semi-analytic waveform model for black hole binaries in scalar environments, which is applied to LIGO-Virgo-KAGRA data to set upper limits on scalar densities and identify tentative evidence for a light scalar field around the GW190728 event.

The Gist

Imagine two black holes dancing around each other, spiraling closer and closer until they crash together. This cosmic waltz creates ripples in space-time called gravitational waves, which detectors like LIGO, Virgo, and KAGRA can "hear."

This paper asks a simple but profound question: What if the black holes aren't dancing in empty space, but are actually wading through a thick, invisible fog?

The Invisible Fog

The authors are looking for a specific type of "fog" made of light scalar particles. Think of these particles as the "ghosts" of the universe. They are a leading candidate for Dark Matter, the mysterious stuff that holds galaxies together but never touches us.

Usually, we think of dark matter as a thin, diffuse gas spread out across the galaxy. But near a spinning black hole, gravity can act like a vacuum cleaner, pulling these particles in and piling them up into a dense, swirling cloud. The paper suggests that if a black hole binary (two black holes orbiting each other) is surrounded by this cloud, the dance changes.

The Dance Floor Analogy

Imagine two ice skaters spinning on a perfectly smooth rink (this is a black hole binary in a vacuum). They spin faster and faster until they collide.

Now, imagine that same rink is covered in a layer of thick, sticky syrup (the scalar field cloud).

• The Drag: As the skaters spin, they have to push through the syrup. This creates friction.
• The Effect: The syrup steals energy from their spin. They lose speed and spiral inward faster than they would on the empty ice.
• The Sound: If you recorded their spin, the "chirp" (the rising pitch of the sound) would change. It would sound slightly different because the syrup is altering their rhythm.

The authors built a mathematical model (a "soundboard") to predict exactly how this syrup changes the gravitational wave signal. They then tested this model against supercomputer simulations (which act like a "wind tunnel" for black holes) to make sure their math was right.

The Detective Work

With their new "soundboard" ready, the team went to the police station (the LIGO-Virgo-KAGRA data catalog) to look at 28 recent black hole collisions. They asked: "Does any of this data sound like it happened in syrup?"

For most of the events, the answer was no. The data looked like the skaters were on clean ice. The team set strict upper limits, saying, "If there was syrup, it couldn't have been thicker than X amount."

However, two cases stood out: GW190728 and GW190814.

• For these two events, the "clean ice" explanation didn't fit the data perfectly.
• The data hinted that the skaters might have been wading through a bit of syrup.
• Specifically, for GW190728, the evidence was "tentative" but intriguing. The statistical tools suggested there was a roughly 3.5 times higher chance that the event happened in a scalar field environment than in a vacuum.

The "Goldilocks" Particle

If this "syrup" is real, what is it made of? The paper suggests it could be a new type of particle with a very specific weight: about $10^{-12}$ electron-volts.

• To put that in perspective, this is incredibly light—billions of times lighter than an electron.
• The authors call this a "light scalar." If it exists, it solves a puzzle in physics and explains where some of the universe's missing dark matter might be hiding.

The Caveats

The authors are careful not to shout "Eureka!" just yet.

• It's just a hint: The evidence is "tentative," meaning it's a strong whisper, not a shout.
• Other possibilities: The "syrup" could be something else entirely, like gas or a different astrophysical effect, though the authors checked and didn't find strong evidence for those.
• The "Syrup" might be thin: The cloud might be very sparse, or it might have been partially eaten away by the black holes before they merged.

The Bottom Line

This paper is a new way of listening to the universe. Instead of just looking at the black holes themselves, the authors are listening to the "air" around them. They found a few events where the air might be a little thicker than we thought, possibly hinting at the existence of a new, ultra-light particle that makes up the dark matter of our cosmos. If confirmed, it would be a massive discovery, proving that black holes can grow "hairs" of dark matter that we can detect from Earth.

Technical Summary

Title: Scalar fields around black hole binaries in LIGO-Virgo-KAGRA

Problem Statement
Light scalar particles are compelling dark matter (DM) candidates arising in extensions of the Standard Model. Gravitational interactions near black holes (BHs) can trigger the growth of dense scalar configurations (e.g., via superradiance), which may persist during the inspiral of compact binaries. If present, these environments alter binary dynamics and imprint signatures on gravitational-wave (GW) signals. While previous studies suggested that close-to-equal-mass binaries (the primary targets of LIGO-Virgo-KAGRA, or LVK) would disrupt such DM structures, recent numerical relativity (NR) simulations and analytic approaches indicate that significant fractions of these clouds can survive, and additional mass may be captured. However, a dedicated search for these effects in LVK GW events using a validated waveform model has been lacking.

Methodology
The authors developed a fast, semi-analytic waveform model for black hole binaries embedded in scalar environments and applied it to the GWTC-3 catalog.

1. Theoretical Framework: The model assumes a real scalar field $\phi$ minimally coupled to gravity. In the non-relativistic limit, the Einstein-Klein-Gordon system reduces to a Poisson-Schrödinger system. The authors employ a perturbative scheme tailored for near-equal-mass binaries ($q \sim 1$), where the scalar field is expressed as a sum of global "gravitational molecule" states of the binary monopole, perturbed by higher multipoles. This approach captures the angular momentum exchange between the binary and the scalar field, leading to a dephasing of the GW signal relative to the vacuum chirp.
2. Waveform Construction: The model incorporates environmental corrections into the co-precessing frame of the IMRPhenomXPHM phenomenological waveform model. It accounts for higher-order modes to improve constraints on asymmetric systems. The phase correction is computed using the stationary phase approximation (SPA) based on the conservation of angular momentum, including the torque exerted by the scalar field.
3. Validation: The waveform model was validated against numerical relativity (NR) simulations using the GRChombo code. These simulations evolved equal-mass BH binaries in scalar field environments. The authors performed Bayesian parameter estimation on injected NR waveforms, comparing the recovery of binary parameters (chirp mass, mass ratio, spins) and the average scalar density ($\bar{\rho}_\phi$) using both the standard vacuum model (IMRPhenomXP) and the new scalar-environment model (IMRPhenomXP_Scalar).
4. Data Analysis: A Bayesian analysis was performed on 28 selected LVK events from the GWTC-3 catalog (satisfying false-alarm rate and inspiral SNR criteria). The analysis used the Bilby package with the Dynesty nested sampling algorithm. Priors were set for the average scalar density and the coupling parameter $\alpha$.
5. Superradiance Interpretation: For events showing potential evidence, the authors re-analyzed the data under physically motivated priors derived from superradiance theory. This imposed constraints on the particle mass $m_\phi$ and the delay time $\tau_d$ between BH formation and merger, ensuring the scalar cloud could survive and reach the late inspiral.

Key Results

• Model Validation: The semi-analytic model successfully reproduced the frequency chirp observed in NR simulations. When applied to injected NR waveforms, the vacuum model exhibited significant biases (e.g., overestimating chirp mass to compensate for accelerated inspiral), whereas the scalar-environment model correctly recovered the binary parameters and inferred the order of magnitude of the scalar density. The Bayes factor favored the scalar model ($\ln B_{vac}^{env} \approx 3.8$) for the injected signal.
• Catalog Analysis: Applying the model to GWTC-3 events yielded upper limits on scalar-field densities for most events, with many bounds lying orders of magnitude below typical superradiant cloud densities. For the majority of events, the posterior distributions were compatible with a vacuum environment.
• GW190728 and GW190814: These two events were exceptions where the vacuum hypothesis lay outside the 90% credible region.
  • GW190728: When analyzed with superradiance-motivated priors, this event showed tentative evidence for a scalar environment with a Bayes factor of $\ln B_{vac}^{env} \approx 3.5$ (for delay times $\tau_d \sim 10^6$ years). This evidence is consistent with a light scalar of mass $m_\phi \sim 10^{-12}$ eV. The environmental model also inferred a lower chirp mass and an effective spin consistent with zero, contrasting with the vacuum model's bias toward higher mass and non-zero spin.
  • GW190814: This event showed only weak evidence for a scalar environment ($\ln B_{vac}^{env} \lesssim 1$) under the same priors.
• Parameter Biases: The analysis demonstrated that neglecting environmental effects can lead to systematic biases in inferred binary parameters, specifically an overestimation of the chirp mass and effective spin.

Significance and Claims
The paper claims to present the first semi-analytic waveform model specifically designed for equal-mass black hole binaries in scalar environments, validated against full numerical relativity. By applying this model to the LVK catalog, the authors provide the first upper bounds on scalar-field environments around compact binaries.

The primary significance lies in the tentative detection of a scalar environment in GW190728. If confirmed, this would point to the existence of a new light scalar particle with a mass of $\sim 10^{-12}$ eV. The authors emphasize that this result is "tentative" and that the vacuum hypothesis cannot be ruled out, particularly given the model's extrapolation into the relativistic regime and potential degeneracies with other astrophysical effects (though eccentricity and tidal deformability were checked and found consistent with vacuum). The work offers a novel, complementary approach to constraining ultralight bosons, distinct from previous methods relying on black hole spin measurements or searches for monochromatic GWs. The authors note that future observing runs and next-generation detectors (Einstein Telescope, Cosmic Explorer) will provide louder signals to enable more stringent tests.