Efficient numerical integration of neutrino oscillations in matter

TL;DR

This paper introduces a Magnus expansion-based solver that significantly accelerates the numerical integration of three-neutrino oscillation equations in matter, facilitating faster data analysis in neutrino physics.

Contribution

The paper presents a specialized, efficient solver for neutrino oscillation equations that outperforms general integrators in speed by up to two orders of magnitude.

Findings

Speed improvement up to 100 times over general integrators

Detailed analysis of numerical procedure and computational costs

Potential to enhance large-scale neutrino data analysis

Abstract

A special purpose solver, based on the Magnus expansion, well suited for the integration of the linear three neutrino oscillations equations in matter is proposed. The computations are speeded up to two orders of magnitude with respect to a general numerical integrator, a fact that could smooth the way for massive numerical integration concomitant with experimental data analyses. Detailed illustrations about numerical procedure and computer time costs are provided.