Search

We introduce a Bayesian protocol based on artificial neural networks that is suitable for modeling inclusive electron-nucleus scattering on a variety of nuclear targets with quantified uncertainties. Unlike previous applications in the field, which directly parameterize the cross sections, our approach employs artificial neural networks to represent the longitudinal and transverse response functions. In contrast to cross sections, which depend on the incoming energy, scattering angle, and energy transfer, the response functions are determined solely by the energy and momentum transfer to the system, allowing the angular component to be treated analytically. We assess the accuracy and predictive power of our framework against the extensive data in the quasielastic inclusive electron-scattering database. Additionally, we present novel extractions of the longitudinal and transverse response functions and compare them with previous experimental analysis and nuclear ab-initio calculations., Comment: 9 pages, 6 figures

We compute ground-state and dynamical properties of $^4$He and $^{16}$O nuclei using as input high-resolution, phenomenological nucleon-nucleon and three-nucleon forces that are local in coordinate space. The nuclear Schr\"odinger equation for both nuclei is accurately solved employing the auxiliary-field diffusion Monte Carlo approach. For the $^4$He nucleus, detailed benchmarks are carried out with the hyperspherical harmonics method. In addition to presenting results for the binding energies and radii, we also analyze the momentum distributions of these nuclei and their Euclidean response function corresponding to the isoscalar density transition. The latter quantity is particularly relevant for lepton-nucleus scattering experiments, as it paves the way to quantum Monte Carlo calculations of electroweak response functions of $^{16}$O., Comment: 14 pages, 9 figures

Neutrino physics is advancing into a precision era with the construction of new experiments, particularly in the few GeV energy range. Within this energy range, neutrinos exhibit diverse interactions with nucleons and nuclei. This study delves in particular into neutrino--nucleus quasi-elastic cross sections, taking into account both standard and, for the first time, non-standard interactions, all within the framework of effective field theory (EFT). The main uncertainties in these cross sections stem from uncertainties in the nucleon-level form factors, and from the approximations necessary to solve the nuclear many-body problem. We explore how these uncertainties influence the potential of neutrino experiments to probe new physics introduced by left-handed, right-handed, scalar, pseudoscalar, and tensor interactions. For some of these interactions the cross section is enhanced, making long-baseline experiments an excellent place to search for them. Our results, including tabulated cross sections for all interaction types and all neutrino flavors, can serve as the foundation for such searches., Comment: 27 pages, 3 figures; v2: Corrected typos, as appears in JHEP

Modeling lepton-nucleus scattering with the accuracy required to extract neutrino-oscillation parameters from long- and short-baseline experiments necessitates retaining most quantum-mechanical effects. One such effect is the interference between one- and two-body current operators in the transition currents, which has been known to enhance the cross-sections, especially in transverse kinematics. In this work, we incorporate such interference in the spectral function formalism, which combines relativistic currents and kinematics with an accurate description of the initial target state. Our analysis of lepton-scattering off $^{12}$C demonstrates that interference effects appreciably enhance the transverse electromagnetic response functions and the flux-folded neutrino-nucleus cross section, in both cases, improving the agreement with experimental data. We discuss the impact on the neutrino-oscillation program and the determination of nucleon axial form factors., Comment: 11 pages, 5 figures

Progress in neutrino-nucleus cross section models is being driven by the need for highly accurate predictions for the neutrino oscillation community. These sophisticated models are being developed within a microscopic description of the nucleus with the goal of encompassing all reaction modes relevant for the accelerator neutrino program. The disconnect between these microscopic models and the event generators that will be used in the next generation of experiments represents a critical obstacle that must be overcome in order to precisely measure the neutrino oscillation parameters. To this end we have developed a Fortran wrapper for lepton-nucleus quasielastic (QE) scattering within the GENIE event generator as a proof of principle, with the broader goal of creating an efficient pipeline for incorporating advanced theoretical models in event generators. As a demonstration of this interface, we have implemented the Spectral Function model into GENIE, offering a more complete description of the nuclear ground state, as well as the ability to provide quantifiable theoretical uncertainties. We validate this implementation and compare its predictions against data and against QE models already available in GENIE., Comment: 10 pages, 9 figures

This review paper emphasizes the significance of microscopic calculations with quantified theoretical error estimates in studying lepton-nucleus interactions and their implications for electron-scattering and accelerator neutrino-oscillation measurements. We investigate two approaches: Green's Function Monte Carlo and the extended factorization scheme, utilizing realistic nuclear target spectral functions. In our study, we include relativistic effects in Green's Function Monte Carlo and validate the inclusive electron-scattering cross section on carbon using available data. We compare the flux folded cross sections for neutrino-Carbon scattering with T2K and MINER$\nu$A experiments, noting the substantial impact of relativistic effects in reducing the theoretical curve strength when compared to MINER$\nu$A data. Additionally, we demonstrate that quantum Monte Carlo-based spectral functions accurately reproduce the quasi-elastic region in electron-scattering data and T2K flux folded cross sections. By comparing results from Green's Function Monte Carlo and the spectral function approach, which share a similar initial target state description, we quantify errors associated with approximations in the factorization scheme and the relativistic treatment of kinematics in Green's Function Monte Carlo., Comment: 30 pages, 9 figures

Microscopic calculations of neutrino-nucleus scattering cross sections are critical for the success of the neutrino-oscillation program. In addition to retaining nuclear correlations in the initial and final state of the reaction, they are based on consistent nuclear interactions and transition current operators, thereby enabling robust uncertainty quantification. In this work, we address a significant limitation of these microscopic methods, which arises from their nonrelativistic nature. By performing the calculations in a reference frame that minimizes nucleon momenta and utilizing the so-called ``two-fragment'' model, we extend the applicability of Green's function Monte Carlo calculations of neutrino-nucleus scattering to higher momenta than currently possible. To validate this approach, we compare our theoretical predictions against inclusive data measured by the MiniBooNE, T2K, and MINER$\nu$A experiments., Comment: 12 pages, 7 figures

The Green's function Monte Carlo (GFMC) method provides accurate solutions to the nuclear many-body problem and predicts properties of light nuclei starting from realistic two- and three-body interactions. Controlling the GFMC fermion-sign problem is crucial, as the signal-to-noise ratio decreases exponentially with Euclidean time, requiring significant computing resources. Inspired by similar scenarios in lattice quantum field theory and spin systems, in this work, we employ integration contour deformations to improve the GFMC signal-to-noise ratio. Machine learning techniques are used to select optimal contours with minimal variance from parameterized families of deformations. As a proof of principle, we consider the deuteron binding energies and Euclidean density response functions. We only observe mild signal-to-noise improvement for the binding energy case. On the other hand, we achieve an order of magnitude reduction of the variance for Euclidean density response functions, paving the way for computing electron- and neutrino-nucleus cross-sections of larger nuclei., Comment: 21 pages, 10 figures

To achieve its design goals, the next generation of neutrino-oscillation accelerator experiments requires percent-level predictions of neutrino-nucleus cross sections supplemented by robust estimates of the theoretical uncertainties involved. The latter arise from both approximations in solving the nuclear many-body problem and in the determination of the single- and few-nucleon quantities taken as input by many-body methods. To quantify both types of uncertainty, we compute flux-averaged double-differential cross sections using the Green's function Monte Carlo and spectral function methods as well as different parameterizations of the nucleon axial form factors based on either deuterium bubble-chamber data or lattice quantum chromodynamics calculations. The cross-section results are compared with available experimental data from the MiniBooNE and T2K collaborations. We also discuss the uncertainties associated with $N\rightarrow \Delta$ transition form factors that enter the two-body current operator. We quantify the relations between neutrino-nucleus cross section and nucleon form factor uncertainties. These relations enable us to determine the form factor precision targets required to achieve a given cross-section precision., Comment: Minor changes to text and figure labels

We present a novel lepton-nucleus event generator: ACHILLES, A CHIcagoLand Lepton Event Simulator. The generator factorizes the primary interaction from the propagation of hadrons in the nucleus, which allows for a great deal of modularity, facilitating further improvements and interfaces with existing codes. We validate our generator against high quality electron-carbon scattering data in the quasielastic regime, including the recent CLAS/e4v reanalysis of existing data. We find good agreement in both inclusive and exclusive distributions. By varying the assumptions on the propagation of knocked out nucleons throughout the nucleus, we estimate a component of theoretical uncertainties. We also propose novel observables that will allow for further testing of lepton-nucleus scattering models. ACHILLES is readily extendable to generate neutrino-nucleus scattering events., Comment: 19 pages, 13 figures

The nuclei $^{151}$Sm and $^{171}$Tm have been identified as attractive candidates for the detection of the cosmic neutrino background. Both isotopes undergo first-forbidden non-unique beta decays which inhibits a prediction of their spectral shape using symmetries alone and this has, so far, obstructed a definitive prediction of their neutrino capture cross sections. In this work we point out that for both elements the so-called "$\xi$-approximation" is applicable and this effectively reduces the spectral shape to deviate by at most $1\%$ from the one that would arise if beta decays were of the allowed type. Using measured half-lives we extract the relevant nuclear matrix element and predict the neutrino capture cross sections for both isotopes at $1\%$ level, accounting for a number of relevant effects including radiative corrections and the finite size of the nuclei. We obtained $(1.12\pm 0.01)\times 10^{-46}{\rm cm}^2$ for $^{171}$Tm and $(4.77\pm 0.01)\times 10^{-48}{\rm cm}^2$ for $^{151}$Sm. This method is robust as it does not rely on the data points near the end-point of the beta spectrum which may be contaminated by atomic physics effects, namely shake-up and shake-off. Finally, we calculate the target mass which is necessary for cosmic neutrino discovery and discuss several bottlenecks and respective solutions associated to the experimental program. We conclude that the detection of cosmic neutrino background by neutrino capture on $^{151}$Sm and $^{171}$Tm is achievable and free from theoretical limitations but still subject to technical issues that should be further investigated by the experimentalists in the context of the proposed PTOLEMY project., Comment: published version; 10 pages, 2 figures

An event generation framework is presented that enables the automatic simulation of events for next-generation neutrino experiments in the Standard Model or extensions thereof. The new generator combines the calculation of the leptonic current based on an automated matrix element generator, and the computation of the hadronic current based on a state-of-the-art nuclear physics model. The approach is validated in Standard-Model simulations for electron scattering and neutrino scattering. Furthermore, the first fully-differential neutrino trident production results are shown in the quasielastic region., Comment: 13 pages, 8 figures

We perform first-principle calculations of electron-nucleus scattering on $^3$He and $^3$H using the Green's function Monte Carlo method and two approaches based on the factorization of the final hadronic state: the spectral-function formalism and the short-time approximation. These three methods are benchmarked among each other and compared to the experimental data for the longitudinal and transverse electromagnetic response functions of $^3$He, and the inclusive cross sections of both $^3$He and $^3$H. Since these three approaches are based on the same description of nuclear dynamics of the initial target state, comparing their results enables a precise quantification of the uncertainties inherent to factorization schemes. At sufficiently large values of the momentum transfer, we find an excellent agreement of the Green's function Monte Carlo calculation with experimental data and with both the spectral-function formalism and the short-time approximation. We also analyze the relevance of relativistic effects, whose inclusion becomes crucial to explain data at high momentum and energy transfer., Comment: 16 pages, 8 figures

The ground-breaking works of Weinberg have opened the way to calculations of atomic nuclei that are based on systematically improvable Hamiltonians. Solving the associated many-body Schr\"odinger equation involves non-trivial difficulties, due to the non-perturbative nature and strong spin-isospin dependence of nuclear interactions. Artificial neural networks have proven to be able to compactly represent the wave functions of nuclei with up to $A=4$ nucleons. In this work, we extend this approach to $^6$Li and $^6$He nuclei, using as input a leading-order pionless effective field theory Hamiltonian. We successfully benchmark their binding energies, point-nucleon densities, and radii with the highly accurate hyperspherical harmonics method., Comment: 15 pages, 3 figures, invited contribution to the special issue in Few-Body Systems "Celebrating 30 years of Steven Weinberg's papers on Nuclear Forces from Chiral Lagrangians"

A variety of new physics scenarios allow 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., Comment: 17 pages, 6 figures. v2: minor changes, references added and typos fixed. Matches version accepted for publication in JHEP

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., Comment: 12 pages, 10 figures

A neutrino community workshop was held at Fermilab in Jan 2020, with the aim of developing an implementation plan for a set of common interfaces to Neutrino Event Generators. This white paper summarizes discussions at the workshop and the resulting plan.

We propose a novel approach to intranuclear cascades which takes as input quantum MonteCarlo nuclear configurations and uses a semi-classical, impact-parameter based algorithm to modelthe propagation of protons and neutrons in the nuclear medium. We successfully compare oursimulations to available proton-carbon scattering data and nuclear-transparency measurements. Byanalyzing the dependence of the simulated observables upon the ingredients entering our intranuclearcascade algorithm, we provide a quantitative understanding of their impact. Particular emphasisis devoted to the role played by nuclear correlations, the Pauli exclusion principle, and interactionprobability distributions., Comment: 14 pages, 10 figures

The complexity of many-body quantum wave functions is a central aspect of several fields of physics and chemistry where non-perturbative interactions are prominent. Artificial neural networks (ANNs) have proven to be a flexible tool to approximate quantum many-body states in condensed matter and chemistry problems. In this work we introduce a neural-network quantum state ansatz to model the ground-state wave function of light nuclei, and approximately solve the nuclear many-body Schr\"odinger equation. Using efficient stochastic sampling and optimization schemes, our approach extends pioneering applications of ANNs in the field, which present exponentially-scaling algorithmic complexity. We compute the binding energies and point-nucleon densities of $A\leq 4$ nuclei as emerging from a leading-order pionless effective field theory Hamiltonian. We successfully benchmark the ANN wave function against more conventional parametrizations based on two- and three-body Jastrow functions, and virtually-exact Green's function Monte Carlo results., Comment: 4 pages, 2 figures

We have applied the extended factorization scheme to investigate the electroweak pion production on nuclei. The ANL-Osaka model, which was obtained by analyzing the data of $\pi N$, $\gamma N$, $N(e,e'\pi)$ and $N(\nu,\mu\,\pi) $ reactions up to invariant mass $W=$ 2 GeV, is used to generate the matrix elements of current operators relevant to pion-production off the nucleon. Medium effects on the $\Delta$ (1232) component of meson-exchange current are included by using a $\Delta$-nucleus potential determined from the previous $\Delta$-hole model studies of pion-nucleus reactions. Nuclear correlations in the initial target state and in the spectator system(s) are modeled using realistic hole spectral functions. As a first step, we show that the data of $^{12}$C$(e,e')$ up to the $\Delta$ (1232) region can be described reasonably well. The interplay between the pion production and two-body meson-exchange mechanisms is shown to be essential in improving the agreement with the data in the ``dip'' region, between the quasielastic and the $\Delta$ (1232) peaks. Predictions for $^{12}$C$(\nu,\mu\,\pi)$ have also been made. They can be used to estimate pion-emission rates in neutrino-nucleus cross section, which constitutes an important systematic uncertainty to the reconstructed neutrino energy. With further improvements of the Metropolis Monte-Carlo techniques to account for final states comprised of more than two particles, our approach can be employed up to $W=$ 2 GeV, where two-pion production and higher mass nucleon resonances must be included for analyzing the data from accelerator-based neutrino-oscillation experiments., Comment: 16 pages, 11 figures

The factorization scheme, based on the impulse approximation and the spectral function formalism, has been recently generalized to allow the description of electromagnetic nuclear interactions driven by two-nucleon currents. We have extended this framework to the case of weak charged and neutral currents and carried out calculations of the double-differential neutrino-carbon and neutrino-oxygen cross sections using two different models of the target spectral functions. The results, showing a moderate dependence on the input spectral function, confirm that our approach provides a consistent treatment of all reaction mechanisms contributing to the signals detected by accelerator-based neutrino experiments., Comment: 18 pages, 7 figures

The electromagnetic responses obtained from Green's function Monte Carlo (GFMC) calculations are based on realistic treatments of nuclear interactions and currents. The main limitations of this method comes from its nonrelativistic nature and its computational cost, the latter hampering the direct evaluation of the inclusive cross sections as measured by experiments. We extend the applicability of GFMC in the quasielastic region to intermediate momentum transfers by performing the calculations in a reference frame that minimizes nucleon momenta. Additional relativistic effects in the kinematics are accounted for employing the two-fragment model. In addition, we developed a novel algorithm, based on the concept of first-kind scaling, to compute the inclusive electromagnetic cross section of $^4$He through an accurate and reliable interpolation of the response functions. A very good agreement is obtained between theoretical and experimental cross sections for a variety of kinematical setups. This offers a promising prospect for the data analysis of neutrino-oscillation experiments that requires an accurate description of nuclear dynamics in which relativistic effects are fully accounted for., Comment: 10 pages, 7 figures

The electromagnetic responses of carbon obtained from the Green's function Monte Carlo and spectral function approaches using the same dynamical input are compared in the kinematical region corresponding to momentum transfer in the range 300-570 MeV. The results of our analysis, aimed at pinning down the limits of applicability of the approximations involved in the two schemes, indicate that the factorization ansatz underlying the spectral function formalism provides remarkably accurate results down to momentum transfer as low as 300 MeV. On the other hand, it appears that at 570 MeV relativistic corrections to the electromagnetic current not included in the Monte Carlo calculations may play a significant role in the transverse channel., Comment: 8 pages, 6 figures

We developed a novel approach based on a generalization of factorization and nuclear spectral functions, allowing for a consistent treatment of the amplitudes involving one- and 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 quasi elastic and $\Delta$-production peaks. The implications for the analysis of neutrino-nucleus cross sections are discussed.

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 quasi elastic sector. We analyse 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 non relativistic approximation breaks down, is outlined., Comment: Physical Review C, in press. arXiv admin note: text overlap with arXiv:1312.1210

The excitation of two particle-two hole final states in neutrino-nucleus scattering, not taken into account in Monte Carlo simulations for data analysis, has been advocated by many authors as the source of the excess cross section observed by the MiniBooNE collaboration in the quasi elastic channel. We analyse the mechanisms leading to the appearance of two particle-two hole states, and show that interference between the amplitudes involving one- and two-nucleon currents, not consistently included in existing calculations, plays a major role. A novel approach allowing to treat one- and two-nucleon current contributions on the same footing is outlined.

The quantitative description of the effects of nuclear dynamics on the measured neutrino-nucleus cross sections -- needed to reduce the systematic uncertainty of long baseline neutrino oscillation experiments -- involves severe difficulties. Owing to the uncertainty on the incoming neutrino energy, different reaction mechanisms contribute to the cross section measured at fixed energy and scattering angle of the outgoing lepton, and must therefore be consistently taken into account within a unified model. We review the theoretical approach based on the impulse approximation and the use of realistic nucleon spectral functions, allowing one to describe a variety of reaction mechanisms active in the broad kinematical range covered by neutrino experiments. The extension of this scheme to include more complex mechanisms involving the two-nucleon current, which are believed to be important, is also outlined. The impact of nuclear effects on the determination of neutrino oscillation parameters is illustrated by analyzing the problem of neutrino energy reconstruction., Comment: To appear in the special issue of Advances in High Energy Physics on "Neutrino Masses and Oscillations"