The formalism based on factorization and nuclear spectral functions has been generalized to treat transition matrix elements involving two-nucleon currents, whose contribution to the nuclear electromagnetic response in the transverse channel is known to be significant. We report the results of calculations of the inclusive electron-carbon cross section, showing that the inclusion of processes involving two-nucleon currents appreciably improves the agreement between theory and data in the dip region, between the quasielastic and $\mathrm{\Delta }$-production peaks. The relation to approaches based on the independent particle of the nucleus and the implications for the analysis of flux-integrated neutrino-nucleus cross sections are discussed.

We carry out an ab initio calculation of the neutrino flux-folded inclusive cross sections measured on ${}^{12}\mathrm{C}$ by the MiniBooNE and T2K Collaborations in the charged-current quasielastic regime. The calculation is based on realistic two- and three-nucleon interactions, and on a realistic nuclear electroweak current with one- and two-nucleon terms that are constructed consistently with these interactions and reproduce low-energy electroweak transitions. Numerically exact quantum Monte Carlo methods are utilized to compute the nuclear weak response functions, by fully retaining many-body correlations in the initial and final states and interference effects between one- and two-body current contributions. We employ a nucleon axial form factor of the dipole form with ${\mathrm{\Lambda }}_A=1.0$ or 1.15 GeV, the latter more in line with a very recent lattice QCD determination. The calculated cross sections are found to be in good agreement with the neutrino data of MiniBooNE and T2K, and antineutrino MiniBooNE data, yielding a consistent picture of nuclei and their electroweak properties across a wide regime of energy and momenta.

The excitation of two-particle–two-hole final states in neutrino-nucleus scattering has been advocated by many authors as the source of the excess cross section observed by the MiniBooNE Collaboration in the quasielastic sector. We analyze the mechanisms leading to the appearance of these final states and illustrate their significance through the results of accurate calculations of the nuclear electromagnetic response in the transverse channel. A novel approach, allowing for a consistent treatment of the amplitudes involving one- and two-nucleon currents in the kinematical region in which the nonrelativistic approximation breaks down, is outlined, and its preliminary results are reported.

A variety of new physics scenarios allows for neutrinos to up-scatter into a heavy neutral lepton state. For a range of couplings and neutrino energies, the heavy neutrino may travel some distance before decaying to visible final states. When both the up-scattering and decay occur within the detector volume, these “double bang” events produce distinctive phenomenology with very low background. In this work, we first consider the current sensitivity at Super-Kamiokande via the atmospheric neutrino flux, and find current data may already provide new constraints. We then examine projected future sensitivity at DUNE and Hyper-Kamiokande, including both atmospheric and beam flux contributions to double-bang signals.

A microscopic description of the interaction of atomic nuclei with external electroweak probes is required for elucidating aspects of short-range nuclear dynamics and for the correct interpretation of neutrino oscillation experiments. Nuclear quantum Monte Carlo methods infer the nuclear electroweak response functions from their Laplace transforms. Inverting the Laplace transform is a notoriously ill-posed problem; and Bayesian techniques, such as maximum entropy, are typically used to reconstruct the original response functions in the quasielastic region. In this work, we present a physics-informed artificial neural network architecture suitable for approximating the inverse of the Laplace transform. Utilizing simulated, albeit realistic, electromagnetic response functions, we show that this physics-informed artificial neural network outperforms maximum entropy in both the low-energy transfer and the quasielastic regions, thereby allowing for robust calculations of electron scattering and neutrino scattering on nuclei and inclusive muon capture rates.