Quasielastic electromagnetic scattering cross sections and world data comparisons in the GENIE Monte Carlo event generator

The usage of Monte Carlo neutrino event generators (MC nu EGs) is a norm within the high-energy nu scattering community. The relevance of quasielastic (QE) energy regimes to nu oscillation experiments implies that accurate calculations of nu A cross sections in this regime will be a key contributor to reducing the systematic uncertainties affecting the extraction of oscillation parameters. In spite of this, many MC nu EGs utilize highly phenomenological, parametrized models of QE scattering cross sections. Moreover, a culture of validation of MC nu EGs against prolific electron (e) scattering data has been historically lacking. In this work, we implement new eA cross sections obtained from nuclear ab initio approaches in GENIE, the primary MC nu EG utilized by the FNAL community. In particular, we utilize results from quantum MC methods which solve the many-body nuclear problem in the short-time approximation (STA), allowing consistent retention of two-nucleon dynamics which are crucial to explain available nuclear electromagnetic (electroweak) data over a wide range of energy and momentum transfers. This new implementation in GENIE is fully tested against the world QE electromagnetic data, finding agreement with available data below similar to 2 GeV of beam energy with the aid of a scaling function formalism. The STA is currently limited to study A <= 12 nuclei, however, its semi-inclusive multibody identity components are exportable to other many-body computational techniques such as auxiliary field diffusion MC which can reach A <= 40 systems while continuing to realize the factorization contained within the STA's multinucleon dynamics. Together, these developments promise to make future experiments such as DUNE more accurate in their assessment of MC nu EG systematics, nu properties, and potentially empower the discovery of physics beyond the Standard Model.