⚛️ high-energy experiments

The Scalar MSW Effect: Compact analytical formulas for neutrino oscillations with large matter effects, including $\mu -\mu$ and $\tau -\tau$ dominant matter potentials

This paper extends the Jacobi diagonalization method to derive compact, accurate analytical formulas for three-flavor neutrino oscillations in matter, revealing that scalar non-standard interactions can induce an energy-independent resonant "scalar MSW" effect and providing effective parameter mappings for both standard and exotic matter potentials.

Original authors: Sandhya Choubey, Andreas Lund

Published 2026-10-02

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Sandhya Choubey, Andreas Lund

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Neutrinos are the most abundant massive particles in the universe, yet they remain among the most elusive. These ghostly particles zip through everything, from the core of the sun to the human body, rarely interacting with the matter they pass through. For decades, physicists have known that neutrinos change their identity, or "flavor," as they travel. A neutrino born as an electron type might arrive at a detector as a muon or tau type. This transformation, known as oscillation, depends on the energy of the particle and the distance it travels. However, when neutrinos move through dense matter, such as the interior of the Earth or the sun, their journey becomes more complex. The electrons in that matter interact with the neutrinos, effectively altering their path and the speed at which they switch flavors. This phenomenon, known as the Mikheyev-Smirnov-Wolfenstein effect, is a well-established part of modern physics, but it only accounts for the standard interactions we currently understand.

Scientists are now hunting for signs of new physics beyond the standard model, and one promising avenue involves looking for subtle deviations in how neutrinos behave in matter. If there are new, undiscovered forces or particles, they might create additional "potentials" that nudge the neutrinos in ways standard physics cannot explain. These new influences could appear in different parts of the mathematical description of the neutrino's journey, potentially affecting muons or taus just as strongly as they affect electrons. The challenge has been that the mathematical tools used to predict these behaviors were built for the standard case, where only electrons matter. When researchers tried to apply these old tools to new scenarios involving muons or taus, the equations became impossibly messy, or simply failed to work. Without clean, simple formulas to describe these new possibilities, it is difficult to design experiments that can spot them or to interpret the data when they do appear.

In a recent study, researchers at the KTH Royal Institute of Technology in Stockholm have developed a new way to solve these complex equations. They revisited a mathematical technique called Jacobi diagonalization, which is essentially a method for untangling a complicated set of interconnected variables until they stand alone and can be understood individually. The key insight of this new work was realizing that the standard way of writing down the neutrino mixing matrix—the mathematical object that describes how the flavors blend—was not the only valid choice. The standard method works perfectly when the new physics affects only the electron, but it breaks down when the new physics targets the muon or the tau. The authors realized that by rearranging the order of the mathematical rotations used to describe the neutrinos, they could align the math with the specific type of matter effect they were studying.

By choosing a custom arrangement for the muon and tau cases, the team was able to derive simple, compact formulas that describe how neutrinos oscillate even when these exotic matter effects are present. These formulas are powerful because they remain accurate even when the new effects are very strong, a regime where older approximation methods usually fail. The researchers applied this new method to a specific scenario involving scalar non-standard interactions, a theoretical model where a new type of field interacts with neutrinos. They found that these interactions could create a resonant enhancement, a sudden boost in the probability of flavor change, similar to the famous resonance seen in standard matter effects. However, they discovered a crucial difference: while the standard resonance depends heavily on the energy of the neutrino, this new "scalar" resonance can occur at a specific strength of interaction regardless of the neutrino's energy. This means the effect would look different in experiments depending on the density of the material the neutrinos are passing through, offering a unique signature to distinguish it from known physics.

The team did not just produce formulas; they rigorously tested how well these formulas matched reality. They compared their analytical results against precise numerical simulations for several major upcoming experiments, including JUNO in China and Hyper-Kamiokande in Japan. For the case where the new interaction affects electrons, their formulas were remarkably accurate, matching the simulations with less than one percent error across a wide range of energies. This high level of precision is vital for experiments like JUNO, which aim to measure neutrino properties with extreme care. The situation was more challenging for the muon and tau cases. Because the standard mathematical method for handling the Earth's ordinary matter effects does not fit neatly with their new approach for these specific particles, the initial formulas showed larger errors. To fix this, the researchers applied a layer of correction, treating the standard matter effects as a small addition to their new solution. This adjustment significantly improved the accuracy, bringing the error down to acceptable levels for many experiments, though it remained less precise for the longest-distance experiments where the standard matter effects are very strong.

The work provides a clear, analytical roadmap for understanding how neutrinos might behave if new forces exist. It demonstrates that by carefully choosing the mathematical language used to describe the problem, scientists can unlock simple solutions to problems that previously seemed intractable. The study confirms that if these scalar interactions exist, they would create a distinct, energy-independent resonance that could be spotted by comparing neutrino behavior in different environments. While the formulas for muon and tau interactions require careful handling of standard matter effects to remain precise, the overall framework offers a robust tool for the next generation of neutrino experiments. By providing these clear, testable predictions, the research helps physicists know exactly what to look for as they search for the hidden layers of the universe that govern these ghostly particles.

Technical Summary: The Scalar MSW Effect: Compact analytical formulas for neutrino oscillations with large matter effects

Problem Statement
Deriving simple and accurate analytical formulas for three-flavor neutrino oscillation probabilities in the presence of matter is a critical challenge in neutrino phenomenology. While exact solutions for the three-generation Hamiltonian with matter effects are analytically intractable, approximate methods are required to interpret data from current and upcoming experiments. Standard perturbation theory often struggles with high matter densities. Conversely, the Jacobi diagonalization method has proven effective for standard matter effects (where the potential is dominated by the electron-electron, $e$ - $e$ , element) but has historically been limited by its reliance on a specific matrix structure. It fails to elegantly handle matter potentials where the dominant contribution arises from other diagonal elements (e.g., $\mu$ - $\mu$ or $\tau$ - $\tau$ ), a scenario relevant for certain Beyond the Standard Model (BSM) physics, specifically Scalar Non-Standard Interactions (SNSI).

Methodology
The authors extend the Jacobi diagonalization method to handle matter potential matrices with a dominant contribution in any single diagonal element ( $e$ - $e$ , $\mu$ - $\mu$ , or $\tau$ - $\tau$ ). The core innovation lies in the strategic parameterization of the neutrino mixing matrix, $U_{\text{PMNS}}$ .

Adaptive Parameterization: The authors argue that the standard Particle Data Group (PDG) parameterization ( $U_e = R_{23}U_\delta R_{13}U_\delta^\dagger R_{12}$ ) is optimal only when the matter potential is symmetric under $R_{23}$ rotations (the $e$ - $e$ case). For potentials dominated by the $\mu$ - $\mu$ element, they utilize a parameterization ( $U_\mu$ ) where $R_{13}$ is the leftmost rotation, aligning with the symmetry of the potential. For $\tau$ - $\tau$ dominance, they employ a parameterization ( $U_\tau$ ) where a modified $R_{12}$ (denoted $\tilde{R}_{12}$ ) is the leftmost factor.
Jacobi Diagonalization: In each case, the Hamiltonian is rotated into an auxiliary basis where the matter potential's symmetry is preserved. The method then applies a sequence of Jacobi rotations to diagonalize the Hamiltonian, systematically eliminating off-diagonal elements.
Approximations: The method relies on neglecting specific off-diagonal elements proportional to $\Delta m^2_{21}$ (the solar mass splitting), which are significantly smaller than terms proportional to $\Delta m^2_{31}$ . This approximation holds to high accuracy ( $O(\epsilon^2)$ ) for relevant neutrino energies.
Perturbative Corrections: For cases where the standard charged-current matter potential ( $A_{CC}$ ) cannot be neglected alongside the new physics potential (specifically for $\delta m_{\mu\mu}$ and $\delta m_{\tau\tau}$ ), the authors apply first-order perturbative corrections to the derived effective parameters.

Key Contributions

Generalization of Jacobi Diagonalization: The paper successfully generalizes the Jacobi method to handle matter potentials of the form $A \propto \text{diag}(0, 1, 0)$ and $A \propto \text{diag}(0, 0, 1)$ , which were previously inaccessible to this specific analytical technique.
Derivation of Effective Parameter Mappings: The authors derive compact, analytical formulas for effective mixing angles ( $\theta^{\text{eff}}_{ij}$ ) and mass-squared differences ( $\Delta m^2_{ij, \text{eff}}$ ) for three distinct SNSI scenarios: $\delta m_{ee}$ , $\delta m_{\mu\mu}$ , and $\delta m_{\tau\tau}$ .
Identification of the Scalar MSW (SMSW) Effect: The derived formulas reveal that diagonal SNSI couplings can induce a resonant enhancement of mixing angles analogous to the standard Mikheyev-Smirnov-Wolfenstein (MSW) effect. Crucially, unlike the standard MSW effect which is energy-dependent, the SMSW resonance condition for $\delta m_{\alpha\alpha}$ can be energy-independent (in the limit of negligible $A_{CC}$ ), offering a distinct phenomenological signature.
Resonance Conditions: Explicit resonance conditions are derived for all three diagonal SNSI cases, showing how they depend on the intrinsic neutrino masses and mixing parameters. The authors note that these resonances can occur for both neutrinos and antineutrinos depending on the sign of the coupling and the specific mixing angle sector.

Results and Accuracy
The authors validate their analytical formulas by comparing them against numerically exact diagonalization results for several experiments: JUNO, Hyper-Kamiokande, ESSnuSB, and DUNE.

$\delta m_{ee}$ Case: The analytical formulas are remarkably accurate (relative error $< 1\%$ ) across all considered experiments and energy ranges. This is because the standard $A_{CC}$ and the SNSI term $\delta m_{ee}$ both enter the $e$ - $e$ element, allowing them to be treated simultaneously within the Jacobi framework.
$\delta m_{\mu\mu}$ and $\delta m_{\tau\tau}$ Cases: The accuracy of the zeroth-order formulas is significantly lower for these cases because the derivation assumes $A_{CC} = 0$ . In long-baseline experiments where $A_{CC}$ is significant, the error can be substantial (up to 100% in DUNE for certain energies).
Improvement via Perturbation: The authors demonstrate that applying first-order perturbative corrections to the effective parameters (or combining their effective parameters with perturbative probability formulas from literature) significantly improves accuracy. For the JUNO experiment, this reduces the error to the sub-1% level. For Hyper-Kamiokande and ESSnuSB, the error is reduced to $\sim 1\%$ . However, for DUNE at high energies ( $E > 5$ GeV), the perturbative approach breaks down as $A_{CC}$ becomes too large relative to the mass splittings.

Significance and Claims
The paper claims that by carefully selecting the parameterization of the mixing matrix based on the symmetry of the matter potential, the Jacobi diagonalization method can be elegantly extended to a broader class of BSM scenarios. The primary significance lies in providing compact, analytical expressions for neutrino oscillations in the presence of large matter effects that are valid for the full three-generation framework.

The authors emphasize that their results for $\delta m_{\mu\mu}$ and $\delta m_{\tau\tau}$ are the first analytical expressions derived using the Jacobi method for these specific parameters. They highlight the discovery of the "Scalar MSW" effect, which offers a unique, potentially energy-independent resonant signature that distinguishes scalar non-standard interactions from standard matter effects. The work provides a necessary toolkit for interpreting data from future high-precision experiments, particularly for distinguishing between standard oscillations and new physics scenarios involving scalar mediators.

More like this