Ultralight dark matter in long-baseline accelerator neutrino oscillations

This paper systematically analyzes ultralight dark matter effects on neutrino oscillations using T2K and NO$\nu$A data, revealing that stochastic fluctuations in the low-mass regime significantly relax coupling constraints while failing to resolve the existing tension in the CP-violating phase $\delta_{CP}$ between the two experiments.

The Gist

Imagine the universe is filled with a mysterious, invisible ocean called Dark Matter. For decades, scientists thought this ocean was made of heavy, slow-moving "fish" (particles like WIMPs). But recently, a new theory suggests the ocean might actually be made of ultralight, ghostly waves (Ultralight Dark Matter, or ULDM) that ripple through everything, including the Earth.

This paper is like a team of detectives (physicists) using two massive underwater microphones (the T2K and NOνA experiments in Japan and the USA) to listen for these ghostly waves. They are specifically listening to neutrinos—tiny, nearly massless particles that zip through the Earth at near light speed.

Here is the breakdown of their investigation, explained simply:

1. The Setup: Neutrinos as Messengers

Neutrinos are like invisible messengers that change their "costume" (flavor) as they travel. Sometimes they are "electron neutrinos," sometimes "muon neutrinos." This changing is called oscillation.

Scientists have been measuring these changes very precisely. However, the two main experiments (T2K and NOνA) are currently disagreeing on a specific detail: the "CP-violating phase" (let's call it the Neutrino Twist). It's like two people looking at the same clock but disagreeing on the exact time. This disagreement is a "tension" in physics that needs solving.

2. The Suspect: The Ghostly Wave (ULDM)

The researchers asked: Could these ghostly dark matter waves be messing with the neutrinos' costumes?

If ULDM exists, it creates a background field that oscillates. Imagine the neutrinos are swimmers in a pool. Usually, the water is calm. But if ULDM exists, the water is gently sloshing back and forth. This sloshing could nudge the swimmers, changing how they move and potentially explaining why the two experiments disagree.

3. The Twist: The "Flickering" Lightbulb

Here is where the paper gets really clever. Previous studies treated this dark matter wave like a steady, constant hum. But the authors realized that because the wave is so light and the universe is so big, the wave isn't steady—it flickers.

• The Analogy: Imagine you are trying to measure the brightness of a lightbulb.
  • High Mass (Heavy Wave): The lightbulb flickers super fast (thousands of times a second). If you take a photo with a slow camera, the flickering blurs out, and you just see a steady average brightness.
  • Low Mass (Light Wave): The lightbulb flickers very slowly. If your camera takes a picture for 10 years, you might catch the bulb in a "dim" phase or a "bright" phase. You can't just average it out; the randomness matters.

The authors realized that for very light dark matter, this random flickering (stochastic fluctuations) is huge. They built a new statistical model to account for this "flickering" rather than assuming a steady hum.

4. The Investigation Results

The team ran the numbers using the latest data from T2K and NOνA. Here is what they found:

• The "Flicker" Matters: In the low-mass range (where the waves flicker slowly), the rules change. Because of the randomness, the scientists had to loosen their "safety net." They found that the limits on how strong the dark matter interaction can be are about 10 times weaker (more relaxed) than previously thought. It's like realizing you can't catch a thief as easily when they are wearing a mask that changes color randomly.
• No Smoking Gun: Despite looking for this ghostly wave, they did not find it. The data still fits the "standard model" (no dark matter interference) just as well as it fits the new theory.
• The Tension Remains: The big hope was that this dark matter wave would fix the disagreement between T2K and NOνA regarding the "Neutrino Twist." It didn't. The two experiments still disagree, and the dark matter wave didn't act as the magic glue to fix it.

5. The Conclusion: Keep Looking

The paper concludes that while this specific type of ultralight dark matter is an interesting idea, the current data from T2K and NOνA doesn't prove it exists.

• The Good News: The scientists have set very strict new rules (constraints) on how this dark matter could behave, especially accounting for that "flickering" effect.
• The Future: To really solve the mystery of the "Neutrino Twist" or catch the ghostly wave, we need even better, more precise microphones. The authors are pointing the finger at future experiments like DUNE and JUNO, which will have the sensitivity to finally see if these ultralight waves are real.

In a nutshell: The scientists checked if invisible, ghostly dark matter waves were messing with neutrino experiments. They realized the waves might be "flickering" randomly, which changes how we look for them. They didn't find the waves, and the waves didn't fix the disagreement between the experiments. But they've drawn a much better map for where to look next.

Technical Summary

1. Problem Statement

The paper addresses the potential impact of Ultralight Dark Matter (ULDM) on neutrino oscillation experiments, specifically focusing on the tension between the T2K and NO$\nu$A experiments regarding the measurement of the CP-violating phase ($\delta_{CP}$) and the atmospheric mixing angle ($\theta_{23}$).

While previous studies have explored ULDM effects on neutrinos, this work identifies two critical gaps in existing analytical frameworks:

1. Coherence and Stochasticity: Most studies model the ULDM field as a single-mode background with static amplitude and phase, averaging out time-dependent effects. However, due to the macroscopic de Broglie wavelength of ULDM, the field is a superposition of many near-degenerate modes. When the experimental exposure time ($T_{exp}$) is shorter than the DM coherence time ($\tau_c$), the field exhibits stochastic fluctuations in amplitude and phase that cannot be averaged out.
2. Neutrino Mass Dependence: In scalar interactions, the induced potential is coupled to the intrinsic neutrino mass matrix. Previous phenomenological studies often neglected the systematic dependence of exclusion limits on the absolute neutrino mass scale ($m_\nu$).

The authors aim to perform a rigorous exclusion analysis of ULDM parameter space using the latest joint datasets from T2K and NO$\nu$A, explicitly incorporating stochastic fluctuations and examining the dependence on the absolute neutrino mass.

2. Methodology

Theoretical Framework

The authors model the ULDM field $\phi(t)$ as a classical oscillating background:
$$\phi(t) = \phi_0 \cos(m_\phi t + \alpha)$$
where $\phi_0$ follows a Rayleigh distribution due to the random walk of many field modes. The analysis covers:

• Scalar Interactions: Flavor-universal ($y_{ee}=y_{\mu\mu}=y_{\tau\tau}$) and flavor-general couplings. The effective potential includes terms proportional to the neutrino mass matrix ($m_\nu$) and the coupling ($y$).
• Vector Interactions: Couplings to $L_e - L_\mu$ and $L_\mu - L_\tau$ gauge symmetries. The potential depends on the charge matrix $Q$ and coupling $g$.

Statistical Treatment of Stochasticity

A key methodological innovation is the Bayesian treatment of stochastic fluctuations.

• The number of independent samples is defined as $N_{DM} = \lfloor T_{exp}/\tau_c \rfloor + 1$.
• Instead of a deterministic $\chi^2$, the authors calculate a marginalized likelihood by integrating over the probability distribution of the field amplitude (Gamma distribution).
• Regimes:
  • Low Mass ($m_\phi \lesssim 10^{-17}$ eV): $N_{DM} \approx 1$. Stochastic fluctuations are maximal; constraints are relaxed.
  • High Mass ($m_\phi \gtrsim 10^{-15}$ eV): $N_{DM} \gg 1$. Fluctuations average out, recovering the deterministic limit.

Data Analysis

• Datasets: Combined latest data from T2K (1.97 $\times$ 10$^{21}$ POT neutrino mode, 1.63 $\times$ 10$^{21}$ POT antineutrino mode) and NO$\nu$A (26.6 $\times$ 10$^{20}$ POT neutrino mode, 12.5 $\times$ 10$^{20}$ POT antineutrino mode).
• Tools: Modified GLoBES (General Long-Baseline Experiment Simulator) to calculate oscillation probabilities including scalar and vector DM potentials.
• Parameters: Solar parameters ($\theta_{12}, \Delta m^2_{21}$) are fixed. The analysis varies $\theta_{23}, \Delta m^2_{31}, \delta_{CP}$, and the DM parameters ($m_\phi$, coupling strength).
• Mass Ordering: Assumes Normal Ordering (NO) with a benchmark lightest neutrino mass $m_1 = 0$ (and tests $m_1 = 10^{-2}$ eV).

3. Key Contributions

1. Explicit Stochastic Modeling: The paper rigorously incorporates the stochastic nature of the ULDM field into the statistical analysis, demonstrating that ignoring these fluctuations leads to overly stringent (and potentially incorrect) constraints in the low-mass regime.
2. Systematic Mass Dependence: The authors quantify how the absolute neutrino mass scale affects scalar ULDM constraints, finding that while the effect exists, it is marginal compared to the stochastic effects.
3. Comprehensive Parameter Scan: The study covers a broad mass range ($10^{-23}$ to $10^{-12}$ eV) for both scalar and vector (specifically $L_e-L_\mu$ and $L_\mu-L_\tau$) interactions, utilizing the most recent joint T2K/NO$\nu$A datasets.
4. Resolution of Tension: The analysis specifically tests whether ULDM interactions can resolve the existing tension between T2K and NO$\nu$A regarding $\delta_{CP}$ and $\theta_{23}$.

4. Key Results

Constraints on Couplings

The inclusion of stochastic fluctuations significantly relaxes constraints in the low-mass regime ($m_\phi \lesssim 10^{-17}$ eV) compared to the deterministic high-mass regime.

• Relaxation Factor: Constraints are relaxed by approximately one order of magnitude (factors of $\sim 8.0 - 8.6$) in the stochastic-dominated regime.
• 3$\sigma$ Exclusion Limits (Low Mass Regime):
  • Scalar (Flavor-Universal): $y_0/m_\phi \lesssim 19.4 \text{ eV}^{-1}$.
  • Vector ($L_e - L_\mu$): $g/m_\phi \lesssim 7.41 \times 10^{-10} \text{ eV}^{-1}$.
  • Vector ($L_\mu - L_\tau$): $g'/m_\phi \lesssim 4.08 \times 10^{-10} \text{ eV}^{-1}$.

Best-Fit Points and Model Preference

Global scans identified best-fit points for all scenarios, but the statistical improvement over the Standard Model (SM) is negligible:

• Scalar: $m_\phi \simeq 6.31 \times 10^{-13}$ eV, $y_0 \simeq 7.93 \times 10^{-13}$. ($\Delta\chi^2 \sim 0.4$)
• Vector ($L_e - L_\mu$): $m_\phi \simeq 1.58 \times 10^{-13}$ eV, $g \simeq 7.93 \times 10^{-24}$. ($\Delta\chi^2 \sim 0.5$)
• Vector ($L_\mu - L_\tau$): $m_\phi \simeq 3.98 \times 10^{-13}$ eV, $g' \simeq 7.93 \times 10^{-24}$. ($\Delta\chi^2 \sim 0.2$)

Conclusion on Tension: The inclusion of ULDM leads to only marginal shifts in the preferred values of $\sin^2\theta_{23}$ and $\delta_{CP}$. It does not provide a statistically significant solution to the tension between T2K and NO$\nu$A. The parameter spaces remain largely consistent with the SM predictions.

Impact of Neutrino Mass

Varying the lightest neutrino mass from $m_1=0$ to $m_1=10^{-2}$ eV resulted in only a marginal relaxation of the excluded regions, confirming that the absolute mass scale is not the dominant factor in the current constraints.

5. Significance

• Methodological Advancement: This work establishes a new standard for analyzing ULDM in neutrino experiments by moving beyond static, deterministic approximations to a stochastic framework that accurately reflects the physical nature of the DM field in the low-mass regime.
• Experimental Guidance: The results indicate that current long-baseline experiments (T2K, NO$\nu$A) are not yet sensitive enough to detect ULDM or resolve the $\delta_{CP}$ tension via these mechanisms.
• Future Prospects: The authors emphasize that future high-precision facilities (such as DUNE and JUNO) are essential. These experiments will have the statistical power to either definitively probe the relaxed ULDM parameter space or achieve a precise measurement of $\delta_{CP}$ independent of BSM physics, thereby clarifying the nature of the current experimental tension.
• Re-evaluation of Limits: The finding that stochastic effects relax constraints by an order of magnitude suggests that previous exclusion limits in the low-mass regime may need to be re-evaluated if they did not account for these fluctuations.