6533b82ffe1ef96bd12951d3

RESEARCH PRODUCT

Efficient numerical integration of neutrino oscillations in matter

Juan Carlos D'olivoJ A OteoFernando Casas

subject

Physics010308 nuclear & particles physicsComputationNumerical analysisFOS: Physical sciencesNumerical Analysis (math.NA)65L05 65L20Computational Physics (physics.comp-ph)Solver01 natural sciencesNumerical integrationHigh Energy Physics - PhenomenologyHigh Energy Physics - Phenomenology (hep-ph)Classical mechanicsIntegratorMagnus expansion0103 physical sciencesFOS: MathematicsApplied mathematicsMathematics - Numerical Analysis010306 general physicsNeutrino oscillationPhysics - Computational PhysicsNumerical stability

description

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.

10.1103/physrevd.94.113008http://arxiv.org/abs/1611.06814