Total: 1
Computing neutrino oscillation probabilities is a solved problem; computing their derivatives is not. We present MANGO (MANGO: A Neutrino Gradient Oscillator), a composable oscillation engine in which every computed quantity is differentiable with respect to all inputs, including propagation geometry, detector depth, and individual Earth-shell densities. None of these quantities appear in traditional analytic probability formulas. The engine supports vacuum, constant-density, layered-PREM, arbitrary-profile, and adiabatic solar propagation, alongside front-ends for non-standard interactions, 3+N sterile states, decoherence, and non-unitary mixing. Because reverse-mode cost depends on output rather than input dimension, evaluating sensitivities incurs a constant 2.5-3x forward-pass overhead. As a result, calculating sensitivities for all 369 density, electron-fraction, and shell-radius parameters of a layered Earth costs no more than the 6 standard oscillation parameters. By maintaining exact sensitivity signals, MANGO allows gradients to flow continuously past the probability stage and through detector response, event weighting, binning, and likelihoods. We demonstrate this on a stylized Earth-tomography analysis. In a single pass, MANGO computes the marginalized uncertainty on a six-zone radial density model and, by differentiating through the inverse Fisher matrix, evaluates its sensitivity with respect to detector angular resolution. This provides an experimental-design metric unreachable using traditional analytic probability formulas. Three-flavor and layered-Earth probabilities match external benchmarks (OscProb, NuFast-Earth) to within 10^-9 to 10^-5, while all forward models, BSM limits, and differentiation paths are verified against exact analytic solutions and finite differences.