Crowdsourcing Gravitational Waves from Superradiant Axions

This paper proposes a crowdsourcing approach to detect ultralight axions by analyzing the collective gravitational wave signals from axion clouds around the entire population of black holes in the Milky Way and the universe, demonstrating that current and future detectors like LIGO, Einstein Telescope, and Cosmic Explorer can robustly probe axion masses in the range of $10^{-13}$ to $10^{-10}$ eV.

The Gist

Imagine the universe is filled with invisible, ghostly particles called axions. These are hypothetical particles that scientists have been hunting for decades because they might explain why the universe looks the way it does and what "dark matter" actually is.

This paper is a new, massive "crowdsourcing" project to find these ghosts. Instead of looking at one or two specific black holes, the authors are asking: "What if we listen to the entire galaxy (and the whole universe) at once?"

Here is the story of their discovery, explained with simple analogies.

1. The Setup: The Black Hole as a Spin-Doctor

Imagine a spinning black hole as a giant, cosmic ice skater spinning on a frozen lake.

• The Axion Cloud: If axions exist, they don't just float by; they get caught in the black hole's "spin." They form a swirling cloud around the black hole, like a gravitational atom.
• The Steal: As the black hole spins, it accidentally "steals" energy from the axion cloud, making the cloud grow bigger and bigger. It's like the skater spinning faster and faster, pulling more and more air into their vortex.
• The Scream: Eventually, this cloud gets so crowded with axions that they start bumping into each other and annihilating. When they do, they burst into Gravitational Waves (ripples in space-time). It's like the cloud finally screaming out a single, pure musical note that we can hear with our detectors.

2. The Old Way vs. The New Way

The Old Way (Looking at Individuals):
Previously, scientists looked at a few famous black holes (like Cygnus X-1) and said, "Hey, you're spinning too fast! If axions existed, you would have slowed down by now." This is like checking the speedometer of one specific car to see if the engine is broken. It works, but it only checks a tiny sample.

The New Way (The Crowdsource):
This paper says, "Let's stop looking at just a few cars. Let's listen to the traffic jam of 100 million black holes in our galaxy, plus billions more in the rest of the universe."

• The Forest of Signals: Instead of one loud scream, we might hear a "forest" of millions of faint, distinct notes (one from each black hole).
• The White Noise: For black holes far away in the rest of the universe, we can't hear them individually. Instead, their combined whispers create a constant, low-level "hiss" or static (a stochastic background) that fills the whole universe.

3. The Detective Work: Listening with LIGO

The authors used the LIGO detector (and future super-detectors) as their "ear." They ran a massive simulation:

• They created a fake universe with 100 million black holes, giving them random masses, spins, ages, and locations.
• They asked: "If axions exist with a specific weight (mass), how many of these black holes would be screaming right now?"
• They calculated exactly what LIGO would hear.

The Results:

• The Sweet Spot: They found that LIGO is perfectly tuned to hear axions with a mass between roughly $10^{-13}$ and $4 \times 10^{-12}$ electron-volts. It's like finding the exact radio frequency where the signal is strongest.
• The "What If" Scenarios: They also looked at "optimistic" scenarios. What if there are lighter black holes than we thought? What if we listen to higher-pitched notes? In these best-case scenarios, they could potentially hear axions up to $10^{-10}$ eV. This is a huge deal because it pushes into a range that other experiments (like microwave ovens designed to catch axions) can't reach yet.

4. The Challenges (The "Systematics")

The paper is very honest about the uncertainties. It's like trying to hear a whisper in a storm.

• The Storm: We don't know the exact age, spin, or location of every black hole. If we guess wrong about how many young, fast-spinning black holes exist, our prediction changes.
• The Analogy: Imagine trying to predict how many birds are singing in a forest. If you guess there are 100 birds but there are actually 1,000, your prediction of the noise level will be way off. The authors tested dozens of different "guesses" (models) to see how much their results would change. They found that while the exact numbers wiggle, the general ability to find axions is very robust.

5. The Future: High-Frequency Detectors

The paper also looks at a new type of detector called a Magnetic Weber Bar.

• The Analogy: LIGO is like a giant drum that hears low, deep booms. The Weber Bar is like a tiny, high-pitched bell.
• Why it matters: Heavier axions would sing at higher pitches. If we can build these high-frequency "bells," we might be able to hear axions that are too heavy for LIGO to detect. The authors show that this could be the key to finding axions that are the "Goldilocks" size for dark matter theories.

The Bottom Line

This paper is a roadmap for the next decade of axion hunting.
It tells us: "Don't just look at the famous black holes. Turn up the volume on the whole universe. If axions exist, the universe is likely screaming with their signal, and our current and future detectors are finally sensitive enough to hear it."

It turns the search for dark matter from a game of "hide and seek" with a few players into a massive "musical chairs" game where the whole universe is the music.

Technical Summary

1. Problem Statement

The existence of ultralight axions (bosons with mass $\mu \sim 10^{-13} - 10^{-10}$ eV) is a leading candidate for solving the Strong CP problem and potentially constituting dark matter. A primary mechanism for detecting such particles is black hole (BH) superradiance. If an axion exists, it can extract rotational energy and angular momentum from spinning black holes, forming a "gravitational atom" (an axion cloud). This cloud eventually decays via axion-axion annihilation ($aa \to g$), emitting nearly monochromatic gravitational waves (GWs).

The Gap: Previous constraints on axion masses have relied primarily on individual black hole observations (e.g., spin measurements of Cygnus X-1 via X-rays or spins inferred from GW mergers like GW231123). These methods are limited by the small number of known BHs and significant systematic uncertainties associated with modeling individual objects. The paper addresses the need to quantify the population-level signal arising from the entire ensemble of black holes in the Milky Way and the universe, which could provide a much more robust and sensitive probe.

2. Methodology

The authors perform a comprehensive population-level analysis of the GW signal generated by axion superradiance, dividing the search into two distinct channels: Galactic (resolved) and Extragalactic (stochastic).

A. Theoretical Framework

• Superradiance Condition: They utilize the condition $\omega_R < m\chi / 2M\bar{r}_+$ (where $\chi$ is spin, $M$ is mass, and $m$ is the azimuthal quantum number) to determine which BHs can support an axion cloud.
• Cloud Evolution: They model the exponential growth of the cloud (timescale $T_c$) and its subsequent decay via GW emission (timescale $T_h$). The analysis focuses primarily on the fastest-growing mode ($n=2, \ell=1, m=1$), while systematically exploring the impact of higher modes ($m > 1$).
• Self-Interactions: The study generally neglects axion self-interactions (dependent on the decay constant $f_a$) for the fiducial results, assuming they are negligible for QCD axions in the lowest modes, though this is noted as a key theoretical uncertainty for higher masses.

B. Galactic Analysis (Resolved Signals)

• Population Modeling: The authors simulate an ensemble of $10^8$ BHs in the Milky Way. They model distributions for:
  • Mass: Fiducial Salpeter-like power law ($5-20 M_\odot$), with variations including Gaussian, exponential, and LIGO-inferred distributions (extending down to $<5 M_\odot$).
  • Spin: Uniform distribution $\chi \in [0, 1]$, with conservative ($[0, 0.5]$) and optimistic ($[0.5, 1]$) variants.
  • Age: A linear uniform distribution based on the Milky Way's star formation history (SFH), contrasting with previous log-flat assumptions.
  • Spatial Distribution: Thin disk, thick disk, and bulge profiles.
• Detection Strategy: They calculate the Signal-to-Noise Ratio (SNR) for each BH. They employ a semi-coherent search strategy (dividing a 1-year observation into 4-hour coherent segments) to account for Earth's rotation and frequency drift.
• Sensitivity Metric: They define sensitivity based on the expected number of signals with SNR $\ge 5$. A mass is considered detectable if the expected count $\lambda \ge 3$ (95% confidence).

C. Extragalactic Analysis (Stochastic Background)

• Population Modeling: They model the cosmic BH population using the Star Formation Rate (SFR), Initial Mass Function (IMF), stellar lifetimes, and metallicity evolution to derive the BH formation rate $R_{BH}(z, M)$.
• Signal Calculation: Instead of resolving individual sources, they compute the stochastic GW energy density spectrum $\Omega_h(f)$ by integrating the GW emission from all BHs throughout cosmic history, accounting for redshift and the non-instantaneous nature of cloud decay.
• Sensitivity Metric: They calculate the SNR for the stochastic background using a 1-year integration time, setting the detection threshold at SNR $\ge 8$.

D. Systematic Uncertainties

The paper rigorously tests the robustness of results against variations in:

• BH mass and spin distributions (including LIGO-inferred distributions).
• BH age and spatial distributions.
• Inclusion of higher superradiant modes.
• Alternative SFR models (e.g., Gamma-Ray Burst inferred).
• The existence of sub-solar mass BHs (via dark matter capture scenarios).

3. Key Contributions

1. Population-Level Sensitivity: Demonstrates that population-based searches are complementary to individual BH studies and can probe a wider range of axion masses with different systematic uncertainties.
2. Systematic Quantification: Provides the first detailed mapping of how astrophysical uncertainties (mass, spin, age, spatial distribution) impact axion sensitivity limits.
3. High-Frequency Probes: Explores the potential of high-frequency detectors (specifically the Magnetic Weber Bar, MWB) to probe axion masses up to $10^{-10}$ eV, a regime where interferometers like LIGO lose sensitivity due to noise floors.
4. Sub-Solar Mass BHs: Investigates the "optimistic" scenario where sub-solar mass BHs exist (formed via dark matter capture), showing this could extend sensitivity to axion masses $\sim 10^{-10}$ eV.

4. Key Results

A. Galactic Sensitivity (LIGO O5)

• Robust Range: LIGO can robustly probe axion masses from $10^{-13}$ eV to $4 \times 10^{-12}$ eV.
• Optimistic Scenarios:
  • If the BH population extends to masses slightly below $5 M_\odot$ (as hinted by recent LIGO inspiral observations) and higher modes ($m>1$) are included, sensitivity could reach $\sim 10^{-11}$ eV.
  • If sub-solar mass BHs ($M \sim 0.075 M_\odot$) exist, sensitivity could push toward $10^{-10}$ eV.
• Systematics: The lower mass bound ($\sim 10^{-13}$ eV) is robust against systematic variations. The upper bound is highly sensitive to the BH mass distribution (lighter BHs allow higher mass axions) and the inclusion of higher modes.

B. Extragalactic Sensitivity

• Stochastic Background: LIGO O5 is projected to be competitive with current constraints in just 1 year of observation for the range $1 \times 10^{-13}$ eV $\lesssim \mu \lesssim 4 \times 10^{-12}$ eV.
• Future Detectors: The Einstein Telescope (ET) and Cosmic Explorer (CE) will significantly extend this reach, with CE potentially probing down to $7 \times 10^{-14}$ eV.
• Robustness: Extragalactic projections are generally robust across systematic variations, particularly for pushing to lower masses.

C. High-Frequency Detectors

• The Magnetic Weber Bar (MWB) offers a unique path to probe axion masses $> 10^{-10}$ eV.
• In optimistic scenarios (sub-solar BHs + higher modes), MWB could reach masses approaching $10^{-10}$ eV, overlapping with the projected sensitivity of axion dark matter searches (haloscopes) and string theory motivated parameter spaces.

5. Significance

• Complementarity: This work establishes that GW population searches offer a distinct and powerful avenue for axion detection, complementary to X-ray spin-down limits and direct dark matter searches. It accesses a different set of systematics, making it a crucial cross-check.
• Discovery Potential: The results suggest that even current and near-future GW detectors (LIGO O5) could discover axions in the $10^{-13} - 10^{-12}$ eV range.
• Theoretical Implications: The potential to probe masses up to $10^{-10}$ eV is significant because this range corresponds to axion decay constants ($f_a \sim 10^{15}-10^{17}$ GeV) predicted by string theory and Grand Unified Theories (GUTs).
• Future Directions: The paper highlights the need for better astrophysical constraints on the BH mass and spin distribution (especially at the low-mass end) and the development of high-frequency GW detectors to fully exploit the superradiance window.

In summary, the paper argues that "crowdsourcing" the gravitational wave signal from the entire population of black holes transforms superradiance from a niche theoretical possibility into a robust, high-priority experimental program for discovering ultralight axions.