GitHub
LinkedIn
ORCIDFlow-Based Surrogates for High-Dimensional Likelihoods in Experimental Neutrino Physics
Precision long-baseline neutrino experiments use near-detector data to constrain systematic uncertainties on the unoscillated neutrino flux, a prerequisite for accurate oscillation parameter measurements at the far detector. When the constrained likelihood is high-dimensional and non-Gaussian, this procedure demands advanced statistical treatment. Here we show that normalizing flows provide faithful and portable likelihood models for this problem. Leveraging an initial Gaussian approximation of the likelihood, we train a hybrid architecture combining coupling transformations and autoregressive spline flows. We demonstrate the method on a representative near-detector likelihood replica with 110 systematic uncertainty parameters, 10 of which explicitly introduce non-Gaussianities in the posterior. The trained model achieves a relative effective sample size of 98%, compared with about 5% for the Gaussian approximation, and reproduces a Markov chain Monte Carlo reference while remaining closed-form, samplable, and pointwise evaluable, making it suited to downstream uncertainty propagation and future near-detector to far-detector fits.
