Probing Quadratically Coupled Ultralight Dark Matter with Pulsar Timing Arrays

This paper characterizes the coherent and stochastic signals produced by ultralight dark matter with quadratic couplings in pulsar timing arrays, finding that while current PTA observations can compete with other probes for coherent signals, they remain less sensitive than equivalence principle tests for stochastic signals.

The Gist

The Big Picture: Listening to the Universe's "Hum"

Imagine the universe is filled with an invisible, ultra-light fog called Ultralight Dark Matter (ULDM). We can't see it, but the authors of this paper suggest it might be "wiggling" or oscillating like a giant, cosmic drumhead.

Usually, scientists think this dark matter interacts with normal matter (like us, the Sun, or stars) only through gravity, like a gentle, invisible hand pushing everything. But this paper asks: What if this dark matter also has a "chemical" connection to normal matter? Specifically, what if it interacts in a way that depends on the square of its strength (a "quadratic" interaction)?

To find out, the authors used Pulsar Timing Arrays (PTAs). Think of a PTA as a galaxy-sized drum set. Pulsars are dead stars that spin incredibly fast and send out radio beams like lighthouses. They are so regular that they act as the universe's most precise clocks. By listening to the "ticks" of these cosmic clocks, scientists can detect if something is messing with the time or the rhythm.

The Two Types of "Wiggles"

The paper explains that if this dark matter interacts quadratically, it creates two very different kinds of signals in the pulsar data:

1. The Coherent Signal (The "Fast Beat"):

   • Analogy: Imagine a single, pure musical note played by a violin. It's steady, rhythmic, and happens at a specific, fast speed.
   • What it is: This is a fast, regular oscillation. The dark matter field wiggles back and forth, causing fundamental constants (like the mass of particles or the strength of forces) to oscillate rapidly.
   • The Result: This creates a predictable "beat" in the pulsar timing data. The authors found that current pulsar data is already very good at hearing this beat, sometimes even better than other Earth-based experiments like atomic clocks.

2. The Stochastic Signal (The "Static Noise"):

   • Analogy: Imagine standing in a crowded room where everyone is whispering randomly. You don't hear a single note; you hear a chaotic, low-frequency "hiss" or static.
   • What it is: This is a slow, messy fluctuation caused by the dark matter waves interfering with each other. It's not a steady beat; it's a random, low-frequency rumble.
   • The Result: The authors found that pulsar arrays are currently not very good at hearing this "static" compared to other methods (like testing the Equivalence Principle). The signal gets lost in the noise.

The "Matter Effect": The Invisible Shield

One of the most important discoveries in the paper is the "Matter Effect."

• The Analogy: Imagine you are trying to hear a whisper from a person standing outside a thick, soundproof concrete wall. If the wall is thin, you hear the whisper clearly. But if the wall is thick and dense, the sound waves get absorbed or blocked before they reach you.
• The Reality: When this dark matter passes near very dense objects like the Sun, the Earth, or a Pulsar, the "quadratic" interaction changes the behavior of the dark matter field. It's as if the dense object creates a "shield" or a "screen" that distorts or suppresses the dark matter signal.
• The Consequence:
  • For the fast beat (coherent signal), the Earth's shield doesn't block it much, so we can still hear it.
  • For the slow rumble (stochastic signal), the shield is very effective. It dampens the signal so much that our current pulsar data can't reliably detect it. The authors had to draw "transparent" lines on their graphs to show where their measurements are no longer trustworthy because of this shielding.

The "Clock" vs. The "Spin"

The paper breaks down exactly how the dark matter messes with the pulsars:

1. The Clock Signal: The dark matter changes the "ticking speed" of the atomic clocks on Earth used to measure the pulsars. If the dark matter makes the clock tick faster, the pulsar seems to spin slower.
2. The Spin Signal: The dark matter changes the actual mass and size of the pulsar itself, causing it to physically speed up or slow down (like a figure skater pulling in their arms).
3. The Doppler Signal: The dark matter pushes the Earth and the Sun slightly, changing their speed relative to the pulsar.

The authors found that for the fast signals, the "Clock" effect is the loudest. For the slow signals, the "Doppler" (pushing) effect is the loudest, but it gets blocked by the Matter Effect.

The "Light QCD Axion" Case Study

The paper also applies this logic to a specific type of dark matter candidate called the QCD Axion (a hypothetical particle proposed to solve a problem in particle physics).

• The Finding: If axions exist and interact this way, they would create the same two signals (fast beat and slow noise).
• The Limit: The authors mapped out where we can rule out these axions. They found that for the "slow noise" part, the "Matter Effect" (the shielding by Earth) is so strong that our current telescopes can't see it yet. However, for the "fast beat," pulsar arrays are competitive with other experiments and can already tell us that these axions don't exist in certain mass ranges.

Summary of Conclusions

• We are listening: Pulsar Timing Arrays are powerful tools for hunting this specific type of dark matter.
• We hear the beat, but not the noise: We are very good at detecting the fast, rhythmic signals (coherent), but the slow, random signals (stochastic) are currently too weak or too blocked by the "Matter Effect" to be detected by pulsars alone.
• The Shield is real: Dense objects like the Sun and Earth act as filters that can hide the dark matter signal, a factor that must be accounted for to avoid false conclusions.
• Competition: For the fast signals, pulsar data is now a top-tier detective, competing with the best atomic clocks and gravity tests on Earth.

In short, the paper teaches us how to listen to the universe's hidden "wiggles," warns us that dense planets can act like noise-canceling headphones, and tells us exactly which frequencies we are currently good at hearing.

Technical Summary

Technical Summary: Probing Quadratically Coupled Ultralight Dark Matter with Pulsar Timing Arrays

Problem Statement
Ultralight dark matter (ULDM), with masses $m_\phi \ll 1$ eV, behaves as a classical oscillating field. While gravitational interactions of ULDM with Pulsar Timing Arrays (PTAs) have been extensively studied, non-gravitational interactions remain less explored. Specifically, this work addresses the phenomenology of ULDM that couples quadratically to Standard Model (SM) particles. Unlike linear couplings, which induce oscillations in fundamental constants at frequency $\omega \approx m_\phi$, quadratic couplings (e.g., $\phi^2$) generate signals at both coherent frequencies ($\omega = 2m_\phi$) and stochastic low-frequency modes ($\omega \lesssim m_\phi \sigma^2$). Furthermore, the propagation of such fields through finite-density astrophysical objects (Sun, Earth, pulsars) introduces "matter effects" that can significantly distort the field profile and suppress observable signals. The paper aims to characterize these signals, assess the sensitivity of current and future PTAs, and determine the viability of these probes against other experimental constraints.

Methodology
The authors develop a comprehensive framework to model PTA timing residuals induced by quadratically coupled scalar fields.

1. Signal Characterization:

   • Lagrangian: The interaction is defined by a scalar field $\phi$ coupled to SM operators (gluons, photons, fermion masses) via a dimensionless quadratic operator $\phi^2/2M_{pl}^2$. This induces fractional variations in fundamental constants ($\Lambda_{QCD}, \alpha, m_\psi$).
   • Signal Types: Three distinct PTA signals are derived:
     • Doppler Signal: Arises from the force on Solar System objects and pulsars due to spatial gradients of the field ($\nabla \phi^2$).
     • Clock Signal: Arises from the modulation of the Terrestrial Time (TT) frequency standard (atomic clocks) due to variations in fundamental constants.
     • Pulsar Spin Signal: Arises from variations in the pulsar's moment of inertia due to changes in constituent particle masses.
   • Statistical Properties: The power spectra for these signals are decomposed into fast modes (coherent oscillations at $\omega = 2m_\phi$) and slow modes (stochastic fluctuations at $\omega \lesssim m_\phi \sigma^2$). The authors demonstrate that while the underlying field is Gaussian, the quadratic signal is non-Gaussian; however, for the stochastic regime relevant to PTAs, the deviation from Gaussianity (kurtosis) is negligible, justifying the use of standard Gaussian likelihoods.

2. Matter Effects:

   • The authors solve the Klein-Gordon equation with a spatially dependent mass term induced by the ambient matter density. This leads to a screening effect where the field amplitude inside massive objects is suppressed.
   • Analytical form factors ($A_{dop}, A_{clk}, A_{psr}$) are derived to quantify the suppression of the Doppler, Clock, and Pulsar Spin signals as a function of the coupling strength and object density.

3. Analysis Framework:

   • Bayesian Inference: The analysis utilizes the PTArcade package within the enterprise framework.
   • Data Sets: Constraints are derived using the NANOGrav 12.5-year data set and projections for a mock 30-year data set (simulating 166 pulsars).
   • Likelihood: The timing residuals are modeled as a combination of white noise, intrinsic red noise, a gravitational wave background (GWB), and the specific ULDM signal models (coherent sinusoidal or stochastic power-law).

Key Contributions

• Unified Signal Model: The paper provides the first complete characterization of both coherent and stochastic PTA signals arising from arbitrary quadratic couplings to the SM, extending previous work that focused primarily on conformal couplings or linear interactions.
• Matter Effect Formalism: It rigorously quantifies how matter effects screen quadratic interactions. The authors derive specific critical coupling thresholds above which PTA signals are suppressed, rendering constraints physically meaningless in those regimes.
• Gaussianity Justification: The work explicitly validates the assumption of Gaussian statistics for stochastic ULDM signals in the PTA frequency band, despite the quadratic nature of the interaction.
• Axion-Specific Constraints: The framework is applied to the Light QCD Axion, deriving constraints in the mass-decay constant plane that account for the specific couplings generated by strong-sector dynamics.

Results

• Coherent Signals ($10^{-24} \text{ eV} \lesssim m_\phi \lesssim 10^{-22} \text{ eV}$):
  • PTA sensitivities for coherent signals (dominated by the Clock signal) compete with and, in certain mass ranges, exceed constraints from equivalence principle tests (MICROSCOPE) and atomic clocks.
  • For the Light QCD Axion, PTA limits are competitive with terrestrial bounds in the low-mass regime.
• Stochastic Signals ($10^{-18} \text{ eV} \lesssim m_\phi \lesssim 10^{-16} \text{ eV}$):
  • The stochastic signal is dominated by the Doppler effect.
  • Crucial Finding: Due to matter effects, the Doppler signal is strongly screened across the entire stochastic mass range for physically relevant couplings. Consequently, current and projected PTA data do not provide physically meaningful constraints on stochastic quadratic interactions in this regime, as the signal is suppressed before it can be detected.
  • In contrast, equivalence principle tests (MICROSCOPE) remain sensitive in this regime, particularly when attractive potentials induce phase transitions inside the Earth, leading to static field profiles and enhanced forces.
• Dilaton-like Scalars: Constraints on general dilaton couplings ($d_g, d_\gamma, d_m$) are presented. The semi-transparent regions in the exclusion plots indicate where matter effects suppress the signal, highlighting the importance of accounting for screening.

Significance and Claims
The paper claims that while PTAs are powerful probes for coherent quadratic interactions, their utility for stochastic quadratic signals is severely limited by matter effects. The authors emphasize that ignoring these screening effects leads to overestimated sensitivity.

• For Coherent Signals: PTAs offer a competitive, independent probe of quadratic couplings, particularly for the Light QCD Axion, complementing terrestrial experiments.
• For Stochastic Signals: The paper concludes that PTA sensitivities underperform compared to equivalence principle constraints for existing and upcoming data sets. The matter effect effectively "screens" the stochastic signal, preventing PTAs from placing meaningful limits on the stochastic abundance of quadratically coupled ULDM in the $10^{-18} - 10^{-16}$ eV range.
• Methodological Impact: The work establishes a necessary correction (form factors) for any future analysis of quadratic couplings in dense environments, ensuring that exclusion limits are not reported for parameter spaces where the signal is physically suppressed.

In summary, the paper refines the search strategy for quadratically coupled ULDM, demonstrating that while PTAs are valuable for coherent searches, their potential for stochastic searches is currently subdominant to other experimental probes due to the unavoidable influence of matter effects.