$\texttt{SNuDD}$: Solar Neutrinos for Direct Detection

Technical Summary: SNuDD – Solar Neutrinos for Direct Detection

Problem Statement
Direct detection experiments have reached a sensitivity regime where solar neutrinos, previously considered negligible backgrounds, now constitute a significant "neutrino fog" that mimics dark matter signals. However, this background also presents an opportunity: these experiments are evolving into sensitive neutrino observatories capable of studying Standard Model (SM) properties and searching for Beyond Standard Model (BSM) physics. A critical gap exists in the theoretical toolkit required to analyze this data. While existing codes (e.g., NuSQuIDS, NuFAST, PEANUTS) can compute solar neutrino propagation, no public package currently exists that can simultaneously compute the Coherent Elastic Neutrino-Nucleus Scattering (CE$\nu$NS) and Elastic Neutrino-Electron Scattering (E$\nu$ES) rate spectra in the presence of Non-Standard Interactions (NSI). Crucially, such a tool must account for the full phase-correlation of different solar neutrino flavor states and the interplay between propagation effects (in the Sun and Earth) and scattering cross-section modifications.

Methodology: The SNuDD Framework
The authors introduce SNuDD (Solar Neutrinos for Direct Detection), an open-source Python package designed to compute solar neutrino recoil spectra at direct detection experiments. The methodology is built upon a comprehensive formalism for neutrino NSI:

1. NSI Formalism: The code utilizes an effective Lagrangian approach where NSI are parameterized by coupling coefficients $\varepsilon_{\alpha\beta}^{fP}$. The framework focuses on vector interactions with stable matter building blocks (electrons, protons, neutrons), defined by a magnitude $\varepsilon_{\eta,\phi}^{\alpha\beta}$ and relative strength angles $\eta$ (proton-to-neutron) and $\phi$ (proton-to-electron).
2. Propagation Modeling:
   • Solar Propagation: The code solves the time-dependent Schrödinger equation for neutrinos traversing the Sun. It employs an adiabatic approximation (validated for the relevant parameter space, including non-zero $\delta_{CP}$) to evolve the neutrino density matrix from the production point to the solar surface. It accounts for spatial averaging over the solar production volume and decoherence during the vacuum propagation to Earth.
   • Earth Propagation: To model the day-night effect and matter regeneration, SNuDD implements a "slab method," approximating the neutrino path through the Earth as a series of constant-density slabs (typically $\sim$50) based on the PREM Earth model. It numerically diagonalizes the Hamiltonian in each slab to compute the evolution S-matrix. A Magnus expansion is also provided for fast analytical approximations at lower energies.
3. Scattering Cross Sections:
   • CE$\nu$NS: The package computes the generalized cross-section by summing amplitudes over final neutrino flavors, incorporating the Helm form factor and NSI-induced modifications to the coherence factor.
   • E$\nu$ES: The code calculates the neutrino-electron scattering cross-section, including SM $Z$- and $W$-exchange and NSI contact terms. It models atomic binding effects in xenon detectors using a stepping approximation and Relativistic Random Phase Approximation (RRPA) scaling to correct for the reduction in event rates due to bound electrons.
4. Detector Integration: SNuDD includes a Convolver class to apply detector-specific effects, such as energy-dependent efficiency, resolution functions (modeled as Gaussian), and quenching factors (Lindhard) for nuclear recoils. It supports configurations for LZ, XENONnT, and PandaX-4T, as well as projections for future observatories.

Key Contributions

• Software Release: The primary contribution is the release of SNuDD, the first public tool capable of computing both CE$\nu$NS and E$\nu$ES spectra with full flavor coherence and NSI effects for direct detection experiments.
• Unified Analysis: It provides a consistent framework to treat neutrino propagation (Sun + Earth) and interaction physics simultaneously, ensuring that flavor correlations are correctly handled in the trace formalism for the scattering rate.
• Validation: The code's Earth propagation module was cross-validated against NuFAST, showing good agreement. Detector modeling was validated by reproducing expected event counts for LZ, XENONnT, and PandaX-4T using SM predictions.

Results: Constraints on NSI
The authors applied SNuDD to derive new constraints on NSI parameters using data from LZ, XENONnT, and PandaX-4T, and projected sensitivities for future 200 tonne-year xenon detectors (XLZD/PandaX-xT).

• Nucleon NSI Plane ($\phi=0$): The analysis reveals that current direct detection limits (particularly from LZ and XENONnT) are rapidly approaching the sensitivity of global fits combining dedicated neutrino experiments (e.g., COHERENT). The study identifies "blind directions" in the parameter space where NSI effects cancel due to the specific proton-to-neutron ratio of xenon isotopes ($\eta \approx -\pi/5$).
• Charged NSI Plane ($\eta=0$): By combining Nuclear Recoil (NR) and Electron Recoil (ER) channels, the authors demonstrate that direct detection experiments can probe NSI with electrons and protons simultaneously. While NR searches are most sensitive to proton couplings ($\phi=0$) and lose sensitivity for pure electron couplings ($\phi \to \pm \pi/2$), ER searches exhibit the opposite behavior. The combination of both channels effectively covers the blind spots of single-channel analyses.
• Future Sensitivity: Projections indicate that future xenon detectors could improve constraints on electron-NSI parameters ($\varepsilon_{ee}, \varepsilon_{e\mu}, \varepsilon_{e\tau}$) by approximately an order of magnitude, potentially surpassing current global fit intervals.

Significance and Claims
The paper claims that SNuDD enables the direct detection community to move beyond treating solar neutrinos merely as a background. Instead, these experiments can now function as leading observatories for neutrino physics. The authors emphasize that:

1. Direct detection experiments are reaching sensitivities comparable to dedicated neutrino experiments for specific NSI parameters.
2. The combination of NR and ER channels in xenon detectors provides a powerful, complementary probe that mitigates the blind directions inherent in single-channel analyses.
3. SNuDD is essential for integrating direct detection data into future global fits of neutrino physics, thereby placing these experiments within the broader landscape of neutrino phenomenology.

The work concludes that with the availability of SNuDD, the inclusion of direct detection data in global NSI analyses is not only feasible but necessary to fully exploit the physics potential of current and next-generation dark matter detectors.