How Bright in Gravitational Waves are Millisecond Pulsars for the Galactic Center GeV Gamma-Ray Excess? A Systematic Study and Implications for Dark Matter

This study systematically investigates the gravitational wave emission from a population of millisecond pulsars proposed to explain the Galactic Center GeV gamma-ray excess, concluding that while current detectors cannot observe them, next-generation instruments like the Einstein Telescope and Cosmic Explorer may detect these signals to definitively test the pulsar versus dark matter interpretations of the excess.

The Gist

The Big Mystery: What is the "Galactic Glitch"?

Imagine looking at the center of our galaxy, the Milky Way, through a gamma-ray telescope. You see a massive, glowing "glitch" or excess of energy that doesn't quite fit the background noise. Scientists call this the Galactic Center GeV Excess (GCE).

For years, there have been two main suspects for what's causing this glow:

1. The Dark Matter Culprit: Invisible particles (Dark Matter) crashing into each other and exploding into gamma rays.
2. The Millisecond Pulsar Suspects: A hidden crowd of dead stars called Millisecond Pulsars (MSPs). These are neutron stars that spin hundreds of times per second, like cosmic lighthouses beaming radiation.

The Problem: Trying to spot these pulsars with light (gamma rays, radio waves, etc.) is like trying to count individual fireflies in a dense, foggy forest at night. The light from the stars gets blurred, blocked by dust, or confused with other bright sources. We can't see them clearly enough to say, "Yes, it's definitely a crowd of pulsars."

The New Detective Tool: Listening to Gravity

This paper proposes a new way to solve the mystery: Gravitational Waves (GWs).

Think of gravitational waves as ripples in the fabric of space-time, like the ripples you make when you skip a stone in a pond. Usually, we detect these ripples from massive collisions (like black holes smashing together). But this paper asks: Can we hear the ripples from a single, spinning neutron star?

If a neutron star isn't perfectly round (if it has a tiny "mountain" or bump on it), its rapid spin creates a steady, rhythmic hum in gravitational waves. It's like a spinning top that wobbles slightly, creating a constant vibration.

The Investigation: How the Authors Solved It

The authors, Ming-Yu Lei, Bei Zhou, and Xiaoyuan Huang, decided to run a simulation to see if these "wobbly" pulsars would be loud enough for our detectors to hear.

1. Building the Crowd (The Population Models)
Since we can't see the pulsars, they had to guess how many there are and what they are like. They built two different "crowd" models:

• Model A (The Copycat): They assumed the pulsars in the galactic center are just like the ones we already know in the rest of the galaxy.
• Model B (The Evolutionary): They assumed these are ancient, older pulsars that have lived in the center for billions of years, changing their properties over time.

2. The "Wobble" Factor (Ellipticity)
For a star to make gravitational waves, it needs a "wobble" (called ellipticity). The authors tested three reasons why a star might wobble:

• The Magnetic Stretch: Imagine the star's internal magnetic field is so strong it pulls the star into an egg shape.
• The Crustal Mountain: Imagine the star's solid crust has a tiny, frozen mountain on it (maybe formed when the star was "reborn" by eating a companion star).
• The Energy Leak: A theoretical scenario where the star loses energy only by making gravitational waves (the "maximum possible" wobble).

3. The Search Strategy (Coherent vs. Incoherent)
They looked at how we would listen for these signals:

• Coherent Search (The Perfect Ear): This assumes we know exactly where the star is and can correct for the Earth's movement. It's like listening for a specific singer in a choir with perfect noise-canceling headphones.
• Incoherent Search (The Group Ear): This assumes we don't know exactly where the stars are or if they are in binary systems (dancing with a partner). The signal gets smeared out. It's like trying to hear a group of people shouting in a crowded room; you can't pick out one voice, but you can hear the total volume of the crowd getting louder.

The Results: Will We Hear Them?

The team ran the numbers and found some exciting, but cautious, results:

• Current Detectors (LIGO/Virgo): With our current technology, the "hum" from these pulsars is too quiet to hear. It's like trying to hear a whisper in a hurricane. Even in the most optimistic "maximum wobble" scenario, we might only catch a few of the loudest ones.
• Future Detectors (Einstein Telescope & Cosmic Explorer): This is where it gets interesting. The next generation of gravitational wave detectors (planned for the 2030s) will be like super-sensitive microphones.
  • The Verdict: These future machines might actually be able to hear a fraction of these pulsars! If they do, it proves the "Pulsar Suspect" is guilty, and the "Dark Matter Culprit" might be innocent.
  • The Twist: If the future detectors don't hear anything, it puts a very strict limit on how "wobbly" these stars can be. If they aren't wobbly enough to make waves, then maybe they aren't the source of the gamma-ray glow after all, which would make the Dark Matter theory more likely.

The Takeaway

This paper is a roadmap for the future. It tells us that while we can't see the hidden crowd of pulsars with our eyes (telescopes), we might be able to hear them with our ears (gravitational wave detectors).

• If we hear them: We solve the mystery of the Galactic Center glow and confirm it's just a crowd of spinning dead stars.
• If we don't: We learn that these stars are perfectly smooth, which forces us to look harder at Dark Matter as the cause.

Either way, listening to the "gravity song" of the Milky Way's center is a brand new way to solve one of astronomy's oldest riddles.

Technical Summary

1. Problem Statement

The Galactic Center GeV Gamma-Ray Excess (GCE) is an unexplained surplus of gamma rays centered on the Milky Way's core. Two primary competing hypotheses exist:

1. Dark Matter (DM): The signal arises from the annihilation of Weakly Interacting Massive Particles (WIMPs).
2. Astrophysical (Millisecond Pulsars - MSPs): The signal arises from an unresolved population of recycled neutron stars (MSPs) in the Galactic bulge.

The Challenge: Distinguishing between these two interpretations via electromagnetic (EM) observations is extremely difficult due to source confusion, pulse broadening from interstellar scattering, and high extinction toward the Galactic Center (GC). Only one candidate MSP has been tentatively identified.
The Proposed Solution: Gravitational Waves (GWs) offer a complementary probe. A steadily rotating, non-axisymmetric MSP emits a continuous, nearly monochromatic GW signal at $f_{GW} \approx 2f_{rot}$. Unlike EM signals, GWs are unaffected by dust or scattering. This paper investigates whether the hypothetical MSP population responsible for the GCE produces detectable GW signals with current and future detectors.

2. Methodology

The authors conducted a systematic study involving three main components: modeling the ellipticity of MSPs, modeling the population properties, and calculating detection strategies.

A. Origins of Ellipticity ($\epsilon$)

The GW strain amplitude depends on the neutron star's ellipticity ($\epsilon$). The authors modeled three distinct physical scenarios:

1. Magnetic-Field-Induced Deformation: Ellipticity caused by strong internal magnetic fields ($B_{int}$) misaligned with the rotation axis. They assumed a conservative ratio $B_{int} = 150 B_{surf}$ and derived $\epsilon$ based on observed spin parameters ($P, \dot{P}$).
2. Crustal Mountains: Ellipticity arising from "mountains" formed by accretion-driven temperature gradients, limited by the crust's breaking strain. They adopted a constant long-term residual ellipticity of $\epsilon \sim 10^{-9}$.
3. Spin-Down Energy Fraction: A model-independent parametrization where a fraction $\eta$ of the spin-down power is converted into GWs. They used a benchmark of $\eta = 1\%$ and a theoretical maximum of $\eta = 100\%$ (the spin-down limit).

B. Population Models

Two frameworks were used to model the MSP population in the Galactic bulge:

1. Empirical Disk-Proxy Model (Pop. I): Assumes bulge MSPs share properties with the well-characterized Galactic disk population. Derived from the ATNF pulsar catalog (filtered for $f_{rot} > 100$ Hz and excluding spin-ups). Estimated total population $N_{MSP} \approx 2.9 \times 10^5$.
2. Evolutionary Model (Pop. II): A theoretical model where bulge MSPs evolve from universal birth parameters over the age of the bulge (following Ploeg et al.). This yields an older, dimmer population with $N_{MSP} \approx 1.6 \times 10^4$.

Both models were normalized to reproduce the observed GCE luminosity and spatially distributed using a "boxy bulge" morphology (consistent with stellar mass distributions).

C. Detection Strategies

The authors evaluated two search strategies to account for Doppler shifts caused by Earth's motion and potential binary orbits:

1. Coherent Search: Assumes Doppler shifts can be fully corrected (requires precise source location and high computational power). Sensitivity scales as $T_{obs}^{-1/2}$.
2. Incoherent Search: Assumes Doppler shifts cannot be fully corrected (due to unknown binary orbits or computational limits). Signals are broadened in frequency. Sensitivity is degraded by a factor of $(T_{obs}/T_{stack})^{1/4}$. The effective strain is calculated as the quadrature sum of all MSPs falling within a Doppler-broadened frequency bin.

3. Key Contributions

• Systematic Comparison: Unlike previous studies that often focused on single scenarios or specific pulsars, this work systematically compares multiple ellipticity origins and two distinct population models.
• Realistic Population Constraints: The study rigorously normalizes the MSP population size to the GCE flux, moving beyond arbitrary assumptions about the number of sources.
• Doppler Treatment: The paper explicitly addresses the impact of binary orbital motion and Earth's revolution on GW detectability, distinguishing between coherent and incoherent search limits.
• Correction of Previous Literature: The authors clarify that previous optimistic estimates (e.g., Miller & Zhao) were inflated by including PSR J0537-6910, a young pulsar with properties inconsistent with the MSP population, and that other studies (Bartel & Profumo) overestimated collective strain by using linear sums instead of quadrature sums.

4. Results

The results are presented in terms of GW strain amplitude ($h_0$) and effective strain ($h_{eff}$) compared against detector sensitivity curves (LVK O4, Einstein Telescope, Cosmic Explorer).

• Current Detectors (LVK O4):

  • Under conservative assumptions (Magnetic deformation, Crustal mountains, $\eta=1\%$), the predicted signals are below the sensitivity of current LIGO-Virgo-KAGRA detectors.
  • Theoretical Maximum: If $\eta = 100\%$ (all spin-down power goes to GWs), a fraction of the high-ellipticity MSPs in the Pop. I model could be detected by current detectors. A null detection in O4 would place stringent limits on this optimistic scenario.

• Next-Generation Detectors (Einstein Telescope & Cosmic Explorer):

  • Coherent Strategy: The predicted strain amplitudes for the high-ellipticity tail of the population (especially in the Magnetic deformation and $\eta=1\%$ scenarios) fall within the reach of ET and CE.
  • Incoherent Strategy: Only the theoretical maximum ($\eta=100\%$) from the larger Pop. I model marginally exceeds the sensitivity of ET/CE.

• Morphology Independence: The results are robust against changes in the assumed bulge morphology (boxy vs. spherical), as the distance distribution differences are negligible compared to the broad ellipticity distribution.

5. Significance and Implications

• Testing the GCE Origin: This work establishes GW astronomy as a viable, independent method to test the MSP hypothesis for the GCE.
  • Positive Detection: A detection of continuous GWs from the Galactic Center would strongly support the MSP interpretation, implying the GCE is astrophysical, not Dark Matter.
  • Null Detection: If next-generation detectors (ET/CE) fail to detect these signals, it would place severe constraints on the ellipticity of MSPs and their abundance. This would challenge the MSP explanation for the GCE, indirectly strengthening the case for Dark Matter annihilation.
• Astrophysical Constraints: Even without a detection, the study provides upper limits on the ellipticity of recycled neutron stars, informing our understanding of neutron star interiors, magnetic field structures, and crustal physics.
• Future Pathways: The paper outlines a clear roadmap for directed continuous GW searches toward the Galactic Center, optimized for the specific frequency bands and Doppler characteristics of the bulge MSP population.

In summary, the paper concludes that while current detectors cannot see this signal, next-generation GW observatories will either uncover the long-sought population of bulge MSPs or rule out their ability to explain the GCE via GW emission, providing a definitive test for one of the most significant anomalies in high-energy astrophysics.