Abstract
Scattering of high energy particles from nucleons probes their structure, as was done in the experiments that established the non-zero size of the proton using electron beams1. The use of charged leptons as scattering probes enables measuring the distribution of electric charges, which is encoded in the vector form factors of the nucleon2. Scattering weakly interacting neutrinos gives the opportunity to measure both vector and axial vector form factors of the nucleon, providing an additional, complementary probe of their structure. The nucleon transition axial form factor, FA, can be measured from neutrino scattering from free nucleons, νμn → μ−p and \({\bar{\nu }}_{\mu }p\to {\mu }^{+}n\), as a function of the negative four-momentum transfer squared (Q2). Up to now, FA(Q2) has been extracted from the bound nucleons in neutrino–deuterium scattering3,4,5,6,7,8,9, which requires uncertain nuclear corrections10. Here we report the first high-statistics measurement, to our knowledge, of the \({\bar{\nu }}_{\mu }\,p\to {\mu }^{+}n\) cross-section from the hydrogen atom, using the plastic scintillator target of the MINERvA11 experiment, extracting FA from free proton targets and measuring the nucleon axial charge radius, rA, to be 0.73 ± 0.17 fm. The antineutrino–hydrogen scattering presented here can access the axial form factor without the need for nuclear theory corrections, and enables direct comparisons with the increasingly precise lattice quantum chromodynamics computations12,13,14,15. Finally, the tools developed for this analysis and the result presented are substantial advancements in our capabilities to understand the nucleon structure in the weak sector, and also help the current and future neutrino oscillation experiments16,17,18,19,20 to better constrain neutrino interaction models.
Similar content being viewed by others
Main
Form factors measured in scattering processes describe the structure of composite objects. They have been thought to be the Fourier transform of charge distributions in the non-relativistic limit of low negative four-momentum transfer squared, Q2. Although this interpretation is not strictly true21, the slopes of the form factors at Q2 = 0 provide a measure of the mean-squared radius \(\left\langle {r}_{{\rm{A}}}^{2}\right\rangle \) for the particle in the charge species described. Nucleon electric (\({G}_{{\rm{E}}}^{{\rm{N}}}\)) and magnetic (\({G}_{{\rm{M}}}^{{\rm{N}}}\)) form factors are precisely measured in electron–nucleon elastic scattering experiments, enabling the radius of the nucleon22,23 to be inferred. Neutrino scattering measurements yield the analogous axial vector form factor, FA, which characterizes the weak charge distribution. FA is also a key input to neutrino oscillation experiments to precisely measure the neutrino oscillation parameters, including CP violation, and to establish mass hierarchy.
Previous measurements of FA in neutrino scattering were performed by measuring dσ/dQ2 in the reaction νμD → μ−pp in deuterium bubble chambers3,4,5,6,7,8,9. Even in deuterium nuclei, theoretical assumptions10 about the Fermi motion of the bound nucleons, the application of the Pauli exclusion principle in the proton-proton (pp) final state and the nuclear wave function are required to extract FA from these measurements24.
Previous extractions assumed FA to follow the dipole form factor, although more flexible models for the form factors are also available24,25,26. Although there are many efforts12,13,15 to calculate FA for Q2 > 0 from lattice quantum chromodynamics (QCD) with increasing precision, calculations in the Q2 region above 1 GeV2 remain imprecise. Pion electroproduction27 and muon capture measurements on hydrogen28 may also be interpreted as constraints on FA, but again with theoretical uncertainties.
In the charged-current elastic (CCE) reaction on the free nucleon \({\bar{\nu }}_{\mu }{\rm{H}}\to {\mu }^{+}n\), a muon antineutrino elastically scatters off the free proton from the hydrogen atom, turning the neutrino into the more massive positively charged muon μ+ and the proton into a neutron (see section 'Terminology of the signal process'). This reaction is free from the nuclear theory corrections described above in scattering from deuterium (D) and provides a direct measurement of FA. It is also a two-body reaction with a nucleon at rest; therefore, the neutrino direction and the final-state μ+ momentum fully specify the interacting system. The antineutrino energy (\({E}_{\bar{\nu }}\)) and Q2 are reconstructed under the CCE hypothesis (\({Q}_{{\rm{QE}}}^{2}\)) of equation (1) and use only the muon energy (Eμ), momentum (pμ) and angle (θμ) of the muon with the neutrino beam. The free proton in the initial state, and the neutron and muon in the final states have masses denoted as Mp, Mn and mμ, respectively.
The neutron momentum in components parallel and transverse to the neutrino direction are
The neutron in non-CCE background reactions will deviate from the predicted neutron directions due to nuclear effects, different initial state assumptions and final-state mass, making available physics-driven selections to reduce the backgrounds. No statistically significant measurement of the process has been performed so far; the only measurement previous to this work recorded 13 events in a hydrogen bubble chamber 29.
The free nucleon cross-section is given by
where M = (Mn + Mp)/2 is the average nucleon mass, GF is the Fermi coupling constant, θc is the Cabibbo angle, Eν is the neutrino energy, (s − u) = 4MEν − m2 − Q2 and m is the charged lepton mass30. The parameters A, B and C are functions of FA, the two vector form factors \({F}_{{\rm{V}}}^{1}\) and \(\xi {F}_{{\rm{V}}}^{2}\) derived from the proton and neutron electric and magnetic form factors, and a pseudoscalar form factor FP whose effect is suppressed by a factor of m2/M2 in A. The form factors are real assuming T invariance30 and charge symmetry. The vector form factors have been precisely parameterized from results of electron scattering experiments25,26,31,32. The pseudoscalar form factor is predicted from the axial form factor33. The cross-section reported in this paper is a convolution between the free nucleon cross-section and the wide-band neutrino flux34, with restrictions on the muon kinematics due to the detector geometry. FA can be derived from the restricted and flux-convolved cross-section given the other form factors.
The nucleon axial radius may be found from the slope of the small Q2 expansion of FA at Q2 = 0:
where FA(0) = −1.2723 ± 0.0023 is derived from neutron decay measurements35. We define \({r}_{{\rm{A}}}\equiv \sqrt{\left\langle {{r}_{{\rm{A}}}}^{2}\right\rangle }\).
Experiment description
The CCE process on hydrogen was measured with the MINERvA11 detector in a \({\bar{\nu }}_{\mu }\) beam produced at the NuMI neutrino beamline34 at Fermilab with an average energy of 5.4 GeV. Antineutrino interactions are selected by requiring a μ+ and detected neutron signatures in the MINERvA detector.
MINERvA is a segmented scintillator detector with hexagonal planes constructed from strips of triangular cross-section assembled into planes perpendicular to the \({\bar{\nu }}_{\mu }\) beam, and is described in more detail in the Methods. This analysis reconstructs neutrino interactions in the active tracker region of the detector, which consists of the scintillator. The scintillator strips point in either the vertical (X) or one of the ±60° (U,V) directions. This region is fully active, consisting of 128 tracker planes stacked in alternating patterns of XUXV. The alternating orientation enables extraction of a three-dimensional position from the strips when charged particles traverse two or more consecutive planes. Muons produced from charged-current neutrino interactions in the MINERvA detector may exit from the rear and enter the MINOS near detector (ND)36, which is located immediately downstream of the MINERvA detector. The MINOS ND is a fully magnetized scintillator and steel detector that determines the muon’s charge and momentum by measuring its curvature and range.
Only muons in the energy range 1.5 GeV < Eμ < 20 GeV with an opening angle θμ < 20° with respect to the neutrino direction are selected because they can be efficiently measured by the MINOS ND. The vertex is defined to be the beginning of the muon track. Energy deposits from other charged particles, such as protons and π±, can be reconstructed into tracks if they span at least four planes. Photon pairs from π0s can be reconstructed from their electromagnetic showers.
Although neutrons are not directly observable from ionization as charged particles, they produce secondary particles with observable energy deposits when they elastically, quasi-elastically or inelastically scatter in the detector. The dominant interactions produce low-energy protons, which can be observed37. Neutrons also scatter undetectably, for example by inelastically knocking out neutrons from carbon nuclei or elastically scattering from carbon, in ways that change the neutron direction and energy. Monte Carlo simulation studies of single-neutron transport in the MINERvA detector show that the angle between the reconstructed and true neutron directions follows the sum of two exponential distributions, with 68% of the candidates within 12°.
Cross-section extraction
The CCE cross-section is measured in bins of Q2. Control samples, events with neutrons pointing away from the predicted direction, provide data-driven constraints on the background models as a function of Q2. The reconstructed neutron directions from the signal events centre around the predicted direction from the μ+ reconstruction with deviations due to the angular resolution. Additional nuclear effects alter the initial neutron directions when the neutron is produced by the charged-current quasi-elastic (CCQE) neutrino interaction on a bound proton in the carbon nucleus. Fermi motion in the carbon nucleus imparts each bound nucleon with a random initial momentum resulting in the neutron direction further deviating from the two-body calculation. Although the CCQE cross-section is a function of FA as well, measurements of electromagnetic form factors even in nuclei as light as 4He have shown that nuclear effects obscure the relationship and make extractions of FA from measurements CCQE on 12C (refs. 38,39,40) susceptible to uncertain nuclear physics. Additionally, multinucleon knockout, such as two-particle-two-hole (2p2h) reactions, and secondary interactions of outgoing neutrons before exiting the nucleus are also present in carbon. The latter phenomena, collectively termed the final-state interactions (FSI), can change the direction and energy of the neutron, and can produce additional final-state particles, including pions produced through the excitation and subsequent decay of nucleons in the nucleus.
On the basis of the detected final-state particles, events are divided into those with only nucleons in the final state (QELike) and those with mesons present (non-QELike). Although the CCE signal is an exclusive subset of QELike, the carbon CCQE, 2p2h and resonant pion production events may experience FSI and land in either category. Both the signal and the background processes are simulated using a realistic Monte Carlo simulation of the detector based on the GEANT4 (ref. 41) simulation toolkit. The input model for neutrino-nucleus interactions is based on GENIE42 with theory and data-driven modifications43,44,45,46,47,48. Finally, the Fermi gas initial state model used by GENIE has been reweighted into a spectral function (SF)49,50 for a more realistic description of nucleus initial states51.
Samples with predominantly QELike(non-QELike) events are selected by requiring ≤1(>1) energetic shower(s) in the detector. Each sample is subdivided according to the neutron candidates’ opening angles to the predicted neutron direction on hydrogen to constrain the QELike and non-QELike backgrounds due to carbon. An alternative coordinate system can be defined as follows:
where the predicted neutron direction \(\hat{t}\), together with the neutrino direction \({\hat{p}}_{\nu }\) and measured muon direction \({\hat{p}}_{\mu }\), set up an alternative coordinate system \(({\hat{x}}^{{\prime} },{\hat{y}}^{{\prime} },{\hat{z}}^{{\prime} })\). The angular variables δθR and δθP describe the projected angles between the measured neutron candidate direction \(\hat{n}\) and \(\hat{t}\) in the neutrino-muon reaction plane and the orthogonal plane intersecting \({\hat{z}}^{{\prime} }\).
Figure 1 shows the relative contributions of the signal CCE and various backgrounds in each δθP–δθR bin in the QELike sample. Although most neutrons from the CCE interaction are reconstructed in the signal region, some undergo rescattering before they are detected, and are measured in non-signal regions. CCQE interaction is the dominant background in the signal region and the regions immediately surrounding it. Fermi momentum and FSI give CCQE events a broader distribution across δθR and δθP. Other interactions have a much broader distribution and are the main background at large δθP and δθR, but still have a small presence in the signal region.
To isolate CCE interactions, the background events are subtracted on the basis of a Monte Carlo simulation. The angular distributions are predicted by the Monte Carlo, but the data determine the Q2-dependent normalization. The normalization is determined by fitting the ‘QE’ region, which is dominated by CCQE, together with the ‘non-QE’ and the ‘non-QE and mesons’ regions, which have a large fraction of non-CCQE interactions. The normalization as a function of Q2 of each interaction channel is adjusted to provide the best fit for the three regions. The normalizations are validated by comparing to event rates in the ‘QE validation’ and ‘non-QE validation’ regions, which have different admixtures of each process. This comparison has an acceptable χ2 = 13.5 for 9 degrees of freedom. The normalizations are then extrapolated to the signal region, and the excess of events is attributed to CCE scattering. By using the data to normalize the background, any modifications in the CCQE cross-section due to changes in FA from the default model are taken into account. A cross-check of the fit strategy was performed on a proton sample in the beam of neutrinos rather than antineutrinos. Reactions of neutrinos are free of CCE events. The strategy yields good agreement at Q2 > 0.2 GeV/c2, the threshold above which CCQE proton tracks can be reconstructed by MINERvA. The antineutrino sample, on the other hand, reaches lower Q2 because the energy deposited by neutron secondary interactions need only span two planes. Figure 2 shows the antineutrino data with the signal and the predicted background levels after the fit, and the ratio to the post-fit model, as a function of Q2 computed under the hydrogen hypothesis (\({Q}_{{\rm{QE}}}^{2}\)) in the signal region. Figure 3 shows the same distributions in the fit regions. CCQE is the dominant background in the signal and QE fit regions.
The effect of detector resolution on the signal events is corrected using an iterative unfolding method52,53, and the detector efficiency under the muon phase space cut in each bin is assessed using the simulation. The cross-section is calculated using the fully integrated antineutrino flux54,55,56 and the total number of hydrogen atoms in the tracker. This measurement does not correct the efficiency loss as a result of the muon phase space cut, which manifests as a restriction on the range of antineutrino energies allowed into the signal sample at each Q2. Theoretical cross-section predictions need to account for the restricted energy range for comparison with the measurement.
Discussion
This is the first statistically significant measurement, to our knowledge, of the antineutrino CCE scattering on the free proton. We observe 5,580(180) signal events over the estimated background of 12,500 events. The measured cross-section is shown on the left in Fig. 4 in black data points, as a ratio to the cross-section prediction assuming a dipole FA described below. The analysis is dominated by statistical uncertainties at all Q2. Systematic uncertainties arise from the small remaining differences, due in part to the regularization, between the post-fit background prediction in each systematic variation of the input model. The dominant systematic uncertainties in this measurement are the neutron secondary interaction in the detector (4.8%), the normalization in the CCQE cross-section (4.5%), the muon energy scale (4.2% from MINOS and 3.1% from MINERvA), the flux (3.9%), neutron FSI (approximately 3%) and the 2p2h process (2.3%).
Theory prediction of the measured cross-section requires input from the electromagnetic vector form factors, the axial form factor, the muon momentum and angle restrictions described above, and convolution between the free nucleon cross-section with the antineutrino flux. The electromagnetic form factor used in this study assumes the BBBA2005 (ref. 31) parameterization. The axial form factor used by most neutrino experiments and generators36,42,49,57,58 assumes a dipole form, \({F}_{{\rm{A}}}({Q}^{2})={F}_{{\rm{A}}}(0){(1+{Q}^{2}/{M}_{{\rm{A}}}^{2})}^{-2}\), which is an approximation derived from the Fourier transform of an exponential charge distribution. In this ansatz, the shape of FA depends only on the axial mass term MA. A more general form, consistent with QCD, is the z expansion formalism59, which maps the one-dimensional variable t = −Q2 onto a unit circle bounded by \({t}_{{\rm{cut}}}=9{m}_{{\rm{\pi }}}^{2}\), the threshold of three-pion production allowed by the axial current24:
The hydrogen cross-section is fitted using FA from the z expansion with t0 = −0.75 (GeV/c)2, \({k}_{\max }=8\). t0 is chosen so that the Q2 bins with precise cross-section measurements are distributed symmetrically around z = 0. Small variations in t0 have no impact on the fit result. \({k}_{\max }\) was chosen to be as small as possible while still enabling the fit to describe the data, as tested by a χ2 statistic. The fit to data includes a bound on the higher order terms24, such that ∣ak/a0∣ ≲ 5 and, for k > 5, ∣ak/a0∣ ≲ 25/k. This bound is treated as a Gaussian regularization term during the χ2 minimization process with a strength parameter λ. The optimal λ of 0.13 was determined by an L-curve study comparing the minimum χ2 separated into the comparison to the data and the regularization. The behaviour of FA at low Q2 is constrained by FA(0) = −1.2723 ± 0.0023, the axial vector coupling as measured in beta decay. A more detailed discussion of the fitting method can be found in the Methods.
The resulting cross-section fit (in red) is shown on the left of Fig. 4 as ratio to a predicted dipole cross-section with MA = 1.014 GeV/c2, together with the predicted cross-section using FA from the Meyer24 fit (in yellow) on deuterium data and a fit derived jointly from deuterium and pion electroproduction data (BBBA2007, in dotted blue)25. The resulting form factor as a ratio to the dipole form factor is shown on the right. The cross-section ratio scales approximately linearly with FA ratios due to suppression of the A term in equations (4) and (5). The nucleon axial radius from the fit to this result is \({r}_{{\rm{A}}}\equiv \sqrt{\left\langle {{r}_{{\rm{A}}}}^{2}\right\rangle }=0.73(17)\ {\rm{fm}}\).
This result is the first statistically significant measurement, as far as we are aware, of the axial vector form factor on free protons without nuclear corrections or other theoretical assumptions. Theoretical uncertainties from the carbon background have been minimized by data-driven methods. By providing a precise and reliable prediction for the charged-current elastic scattering from nucleons, neutrino measurements on higher Z nuclei can benefit from better constrained nucleon effects to expose the nuclear effects. The method developed in this study will enable future experiments with hydrogen content in the target18,19 to make further measurements of the axial form factor. Future experiments with intrinsic three-dimensional capability would be able to observe the directions of low-energy neutron candidates, and improve the low Q2 measurement with more statistics.
Methods
Terminology of the signal process
This article refers to the \(\bar{\nu }+H\to {\mu }^{+}n\) scattering process as the CCE process to distinguish it from the CCQE process that experiments have measured in scattering from bound nucleons in nuclei. The signal process is an isolated, two-body reaction that conserves the invariant mass, and is ‘elastic’ in the sense that all energy from the neutrino and the target proton is transferred to the final-state muon and neutron.
Detector and neutrino beam
The MINERvA11 detector, shown in Extended Data Fig. 1, is located just upstream of the MINOS ND36 in the NuMI beamline34 at Fermilab. The main injector supplies a beam of 120 GeV protons striking a graphite target 1.04 km upstream of the detector. Secondary mesons produced at the target are charge selected by two magnetic horns to enhance the νμ production in the forward horn current mode and the \({\bar{\nu }}_{\mu }\) production in the reverse horn current (RHC) mode. The resulting neutrino or antineutrino beam travels at a downward angle of 0.059 radian and has a broad energy profile. In particular, the antineutrino beam in the RHC mode peaks at 5.5 GeV (ref. 60). Extended Data Fig. 8 (right) shows the flux distribution.
Extended Data Fig. 2 describes the detector components. The MINERvA detector is capable of measuring \(\bar{\nu }p\to {\mu }^{+}n\) events over a large range of Q2, up to 7 (GeV/c)2 by tracking the muon through the scintillator target and reconstructing the secondary interaction of the produced neutrons. Each plane in the active tracker region consists of scintillator strips with a triangular 33 mm (transverse) by 17 mm (along beam) cross-section interlocked together. Particles typically deposit energy in multiple adjacent strips, which gives position resolution smaller than the size of the strips. To be reconstructed for this analysis, muons produced in the scintillator target must exit from the back of the MINERvA detector and be measured by the MINOS ND36. The magnetized steel in MINOS provides an average 1.3 T field, and the interspersed scintillator planes measure the length and curvature of the muon tracks to determine the muon’s charge and momentum as it enters the MINOS ND. Muon tracks in MINERvA are reconstructed and fitted with a Kalman filter resulting in 3.1 mm tracking resolution. The total muon energy and angle resolutions are 6% and 0.06°, respectively61. Muons well tracked by the MINOS ND usually have accepted angles to the beam θμ < 20° and momenta in the range 1.5 GeV/c < pμ < 20 GeV/c, respectively.
Signal definition
The CCE, \(\bar{\nu }H\to {\mu }^{+}n\), signal is a subset of the quasi-elastic-like (QELike) sample in antineutrino scattering events. Precisely, ‘QELike’ events are those with a μ+ and only nucleons and remnant nuclei in the final state. In addition, events with energetic protons that produce reconstructable tracks, corresponding to a kinetic energy threshold of about 120 MeV, are excluded from the ‘QElike’ category in the antineutrino sample. The remaining charged-current antineutrino events, with other final-state particles or with sufficiently energetic protons, are labelled as ‘non-QELike’. The QELike category contains, apart from the signal events, CCQE interactions on protons bound in carbon nuclei, interactions in which mesons are produced but later absorbed in the nucleus, and multinucleon knockout events (2p2h), in which two target nucleons participate in the scattering process. The antineutrino QELike category is naturally neutron-rich because of the net negative charge transferred to the target. Neutron detection is inefficient and does not give a strong constraint on neutron energy; therefore, it is not possible to exclude classes of events with multiple neutrons by counting detectable neutrons or measuring their energies. Signal events are instead constrained by the consistency of the direction of the neutron with the kinematics of the charged-current elastic scattering signal process.
Reference interaction and detector model
A complete Monte Carlo simulation of the MINERvA detector, data acquisition and event reconstruction chain was developed to realistically predict the detector efficiency and resolution. The detector was simulated with the GEANT4 toolkit41. Signal and background events are generated by the GENIE42 v.2.12.6 neutrino event generator with data-driven modifications. GENIE describes the carbon nucleus with a relativistic Fermi gas (RFG) model that assumes the nucleons are free, with the Fermi momentum maximum of KF = 221 MeV. The single nucleon momentum is extended beyond the Fermi momentum cutoff to account for short-range correlations62. For the CCQE reactions, the RFG model is reweighted to the prediction of a spectral function model50,63 used in NuWro49, which incorporates both the shell structure of the nucleus and the nucleon correlation. Without the spectral function features, the initial state model is not sufficiently consistent with the carbon CCQE background to provide an accurate constraint. The reference CCE and CCQE cross-sections are modelled with a dipole axial form factor with MA = 0.99 GeV/c2 and the BBBA2005 (ref. 31) vector form factors. The Valencia 2p2h model44 used to describe the knockout of multiple nucleons is tuned by adding a large enhancement to match other MINERvA measurements43. The FSI is simulated by the GENIE INTRANUKE-hA42 model. The model emulates a full-scale intranucleus cascade64 in a way that is easily varied for systematic uncertainty studies. Some of the FSI reactions that change the final-state particle content in the detector are the inelastic FSI, pion production, pion absorption and pion and nucleon charge exchange. The elastic component of the FSI simulation in the CCQE channel is replaced with the no-FSI component through reweighting to circumvent an error in the INTRANUKE-hA kinematic implementation65. Non-CCQE interactions, such as resonant pion production, are based on the Rein–Seghal66 model, and the deep inelastic scattering is based on the Bodek–Yang67 model.
Neutron interaction and reconstruction
Neutrons produced in neutrino interactions can interact in the detector through elastic or inelastic interactions on hydrogen or carbon. Neutron capture rates are negligible at MeV energies68. Elastic scattering on hydrogen n1H → np or inelastic scattering on carbon, such as n12C → np11B, can produce protons with enough energy to be observed in the detector. Neutron interactions are modelled using GEANT4. However, GEANT4 predictions of some inelastic neutron-C cross-sections are not in agreement with measurements of neutron scattering on scintillators at energies similar to our signal reaction69; therefore, we reweight the probability of neutron scattering from the GEANT4 prediction to be consistent with the recent models70 that do agree with these data.
The MINERvA experiment has developed algorithms to reconstruct energy deposits from neutrons that produce ionizing particles by scattering off the material in the detector37. The daughter particles, primarily protons, are usually observed as isolated hits in the detector, far from the interaction vertex. If the particle is energetic enough to pass through two or more planes, a three-dimensional position of the interaction can be reconstructed and the direction from the vertex measured. The angle between the reconstructed and true neutron directions was assessed with the detector simulation. The distribution can be modelled as the sum of two exponential distributions \(N(\theta )={A}_{1}/{\Theta }_{1}\exp (-\theta /{\Theta }_{1})+{A}_{2}/{\Theta }_{2}\exp (-\theta /{\Theta }_{2})\), where the values of A1 and A2 are such that 40% of the neutrons have a narrow Θ1 = 6° exponential slope and the remainder have a Θ2 = 20° exponential slope. The heavy tail of this distribution probably results from the undetected scatters of neutrons in the detector before a detected scatter is observed. The intrinsic neutron angular resolution is convoluted with physics-driven neutron direction to enable constraint on the carbon CCQE events.
Extended Data Fig. 2 illustrates the schematics of the angular variables δθR and δθP, and shows raw data and fitted model distributions. The δθP distribution is symmetric about δθP = 0, largely due to Fermi momentum in the CCQE background, whereas for the CCE signal the smearing is dominated by secondary interaction effects. There are bulk shifts in the peak positions of δθR distribution from the non-CCQE components due to the neutron assumption required in the direction prediction. Whereas the transverse component of the predicted neutron momentum is fixed by the muon’s transverse momentum (equation (3)), the component parallel to the neutrino beam depends on the invariant hadron mass of the assumed final state (\({M}_{{\rm{n}}}^{2}\) in equation (2)), resulting in a systematic shift when the final-state hadronic system is not a single neutron. Extended Data Fig. 3 shows the Monte Carlo event rates for a few interaction types at different Q2.
Event selection
Events with a muon track starting in the scintillator tracker region and ending in the MINOS ND with a positive charge identification are selected for analysis. Events with additional tracks formed at the vertex are rejected because they indicate the presence of charged mesons or energetic protons. The recoil energy is defined as the total energy outside a 100 mm sphere from the vertex and not associated with the muon track. A Q2-dependent maximum recoil energy cut was applied to reduce the fraction of non-QELike background (Extended Data Table 1).
As described above, neutrons become detectable when they inelastically interact in the detector to create secondary particles. Neutron candidates in the detector are reconstructed by an algorithm that combines detector hits not attached to an interaction vertex and known particle tracks. Neutron candidates whose hits span two views or more produce three-dimensional position information, and only events with such candidates are selected for this analysis. Even at this early stage, the sample is already dominated by events with a final-state neutron from an antineutrino interaction.
In the CCE \(\bar{\nu }{\rm{H}}\) reaction, the energy transfer q0 from the neutrino to the initial state proton is proportional to the Q2:
Leading neutron candidates with detectable energy above this physics limit cannot originate from the signal reaction, so such events are rejected. The QELike selection also restricts the number of isolated energy deposits outside the vertex region. Events with a single deposit whose energy resembles proton energy deposits with at least 10 MeV per strip are defined to be ‘QELike’. Because of the inefficient neutron selection, events with multiple energy deposits are predominately events with particles other than neutrons, such as photons from π0 and mis-reconstructed π±. Events with multiple candidates are defined to be in the ‘non-QELike’ sideband sample.
Control regions are formulated on the basis of opening angular separations between the reconstructed neutron candidate and the expected neutron direction under the hydrogen hypothesis. The angular variables are defined according to equation (7), with regions formulated with varying fractions of QELike backgrounds, as shown in Fig. 1. The signal selection occupies the central −10° < δθP/δθR < 10° region, whereas the ‘non-QE and mesons’ region encompasses all angular space outside the named regions. A bin-by-bin background constraint is performed for both the QELike and non-QELike events in bins of \({Q}_{{\rm{QE}}}^{2}\), which is the Q2 computed under the hydrogen hypothesis. Equation (10) shows the χ2 function minimized in the fit:
where \({N}_{C,S,i}^{{\rm{mc}}}\) is the Monte Carlo event rate in the ith \({Q}_{{\rm{QE}}}^{2}\) bin in the interaction type category C, and in the angular control sample S. \({N}_{S,i}^{{\rm{data}}}\) is the data rate in the ith \({Q}_{{\rm{QE}}}^{2}\) bin in a control sample S. Each category used in the fit receives a weight wC,i; categories not fitted receive a constant weight of 1. A regularization term is added to the χ2 function to ensure the weights across the bins change smoothly. λS is the regularization strength and is a meta parameter tuned to ensure consistency between the fit and validation regions. To do so, the fit is performed on the central value model for a range of λS. An average χ2 for the validation regions, weighted by their event rates, is then calculated. The χ2 from the QE validation region dominates this calculation due to its high statistics. The value of λS is chosen through an L-curve study. A separate fit using these parameters is done for every source of systematic uncertainty to evaluate its effect on the fit.
The QELike categories constrained are CCQE on carbon, 2p2h and resonant events whose pions are absorbed through FSI. The CCE signal category is not tuned because its contribution is small in the control regions. QE carbon is the primary background type in the QELike sample. Events with mesons in the final states, such as single π0, π± and multiple pions, are primarily constrained with the non-QELike sample. The control regions used for fitting are the ‘QE’, ‘non-QE’ and ‘non-QE and mesons’, in the QELike sample, and all regions combined in the non-QELike sample. The QE and non-QE validation regions in the QELike sample are used to assess the goodness-of-fit. The results of the fit and the ratio to the fitted models in the validation regions are shown in Extended Data Fig. 4. The total event rate in the non-QE validation region is about one-fifth of the QE validation region.
Validation of CCQE model and background subtraction technique using ν μ events
The validity of the background constraint techniques with the above model was tested using a neutrino CCQE sample with no CCE scattering from the free hydrogen. Events with a negatively charged muon matched to MINOS and a reconstructed proton track are selected. As with the CCE reconstruction, we compute the theoretical nucleon direction, this time a final-state proton, on an event-by-event basis and compare it with the proton track direction to classify the event according to the same angular region scheme. The track reconstruction efficiency at MINERvA for protons scales with their momenta, and below 450 MeV/c (100 MeV in kinetic energy) the protons are no longer reconstructed. Therefore, QE protons with Q2 < 0.2 (GeV/c)2 are very rare, with events dominated by 2p2h background events whose reconstructed Q2 based on the QE hypothesis are small.
A fit was performed using the control samples. Extended Data Fig. 5 shows both the fitted neutrino event ratio to the simulation in the signal and QE regions analogous to Fig. 2 and the δθR in this sample. The fitted Monte Carlo adequately describes data when proton tracks from CCQE are reliably reconstructed. The small discrepancy at \({Q}_{{\rm{QE}}}^{2}\) between 0.2 (GeV/c)2 and 0.4 (GeV/c)2 is covered by a 100% uncertainty in the 2p2h component in that restricted region, as shown by the δθR plot at the bottom of the figure, and is the probable source of the discrepancy, given that the 2p2h component only begins to contribute significantly in this \({Q}_{{\rm{QE}}}^{2}\) range. Such a systematic deviation in the carbon model would be well within the systematic uncertainties applied to our CCE background analysis in the antineutrino beam.
Unfolding
The signal event rate obtained after background subtraction contains smearing effects from detector resolution. This analysis, similar to many recent results from the MINERvA Collaboration, undoes the smearing by means of an iterative unfolding algorithm developed by D’Agostini52. The inputs to the algorithm are the smearing matrix obtained through the Monte Carlo study of the true and reconstructed Q2 in the signal sample and the number of iteration steps, Ni, as a regularization parameter. When Ni is too small, the unfolded distribution will not have converged to the true kinematic distribution, but higher values of Ni increase statistical variances of the reconstructed distribution. To find the value of Ni that balances model bias and statistical uncertainties, we form statistical universes of a toy model obtained by reweighting the Monte Carlo signal sample. The toy model differs from the CCE model by about 10% at lower Q2, and 20% at higher Q2, to mimic the shape of the discrepancy between data and the model. An ensemble of pseudo-experiments is created by changing the value in each bin of the toy model to a random number generated by a Poisson distribution, with the mean set at the central value of the toy model in that bin. A total of 1,000 statistical universes have been made and Ni = 4 minimizes the median and mean of the χ2 of the result of the pseudo-experiments compared to the hypothesis of the toy model.
The unfolded data are corrected for the predicted efficiency obtained from the Monte Carlo study. The efficiency, shown in Extended Data Fig. 6, increases with higher Q2 until 0.5 (GeV/c)2 as the rate of inelastic interactions that produce protons over the reconstruction threshold increases. For Q2 greater than 0.5 (GeV/c)2, however, the signal efficiency drops due to two factors. First, the neutron angle with respect to the neutrino direction increases, and therefore the neutron passes through less detector material on average before exiting the detector. Second, the muon is produced at a wider angle with respect to the neutrino and therefore is less likely to be reconstructed in the downstream MINOS ND.
Uncertainties
Extended Data Fig. 7 shows the breakdown of the fractional systematic and statistical uncertainties in our measurement. We assess systematic uncertainties in six broad categories. Modelling uncertainties in GENIE are broken into the FSI and cross-section model uncertainties. The muon reconstruction uncertainties account for resolutions in the muon energy, direction and the reconstruction efficiency. The low recoil fit category is due to uncertainties in the tuned 2p2h model. The ‘other’ category includes target mass uncertainties and a set of particle response uncertainties. The largest particle response uncertainties are related to the neutron interactions in the detector. The ‘neutron reweight’ uncertainty is due to the neutron elastic and inelastic cross-sections assumed in the detector; it is assessed for each neutron and assigned a larger value for lower kinetic energy. The neutron interaction uncertainty accounts for the discrepancy in the energy deposited in the neutron candidates we observed in a previous study37. These are the largest uncertainties in the ‘other’ category, and they mainly affect the low Q2 cross-section. This analysis is dominated by the statistical uncertainty, in part because of the large background subtraction.
Cross-section calculation and the z expansion fit of F A(Q 2)
The cross-section measured in this analysis is the flux-integrated cross-section for the geometry of the MINERvA detector. The measured event and predicted background rates are shown in Supplementary Table 1. The efficiency corrected event rates are divided by the total number of hydrogen atoms and protons on target (Supplementary Table 2). The measured cross-section is shown in tabulated form in Supplementary Table 3, and the total covariance and statistical-only covariance matrices of the measurement are shown in Supplementary Tables 6 and 7, respectively.
Calculating the integrated cross-section in each bin needs to account for the muon phase space cuts described above. The effect of the phase space restriction is manifested as a restricted range of neutrino energy available in the selected sample at each Q2 point, ultimately reducing the differential cross-section at each Q2 bin because the acceptance-corrected event rate is divided by the fully integrated neutrino flux for this measurement. Extended Data Fig. 8 (right) illustrates the accepted neutrino energy at each Q2. Therefore, obtaining a fit of FA requires the convolution of equation (4) with the antineutrino flux from the RHC configuration of the NuMI34 neutrino beam with the Q2-dependent energy cutoff. The vector form factors used in this analysis are parameterized from electron scattering data31 used commonly by neutrino Monte Carlo generators42,49,58.
The axial vector form factor is also fitted using the z expansion59 formalism. We adopt the procedure of a fit to deuterium data by Meyer et al.24. The z expansion for FA is a polynomial of z with coefficients ak reproduced here for convenience.
Values of the coefficients ak are constrained by FA(0) = −1.2723, and a set of four coefficient sum rules (equation (13)) determined by the high Q2 behaviour required by QCD71,72. For a given value of \({k}_{\max }\), \({N}_{a}={k}_{\max }-4\) coefficients are free parameters: \(({a}_{1},...,{a}_{{k}_{\max }-4})\). The error function involved in the fit is a χ2 function with a regularization term of the form \(\lambda \left[{\sum }_{k=1}^{5}{\left(\frac{{a}_{k}}{5{a}_{0}}\right)}^{2}+{\sum }_{k=5}^{{k}_{\max }}{\left(\frac{k{a}_{k}}{25{a}_{0}}\right)}^{2}\right]\) based on the expected fall off of ak ∝ k4 at large k (ref. 24). We choose t0 = −0.75 (GeV/c)2 in our analysis to centre our data points around z = 0. Small variations in t0 have no effect on the final fit results. \({k}_{\max }\) and λ are parameters that affect the functional form of FA and the uncertainties of the form factor. Larger \({k}_{\max }\) increases allow for a more ‘curvy’ FA that might overfit and follow statistical variations of the data, whereas a non-zero λ suppresses this behaviour for larger \({k}_{\max }\). At \({k}_{\max }=6\), the value of λ has little effect on the low Q2 behaviour of the fit but prevents FA from adequately describing the data at larger Q2. The fit at \({k}_{\max }=8\) enables higher Q2 fit to follow the data points, but also exposes the low-Q2 region to a small dependence on λ, which manifests as a few per cent variation in the calculation of nucleon axial radius measured on a proton, sufficiently covered by the uncertainties from the fit. The behaviour of the regularization term and the data χ2 term were studied and the fit at λ = 0.13 was chosen according to the maximum curvature, or L-curve, criterion73. Extended Data Fig. 8 (left) shows the final fit to data. In general, varying the z expansion parameters, such as \({k}_{\max }\),t0 and λ, results in changes below 10% of the total uncertainty. Supplementary Tables 4 and 5 show the fit result and the correlation matrix using \({k}_{\max }=8,\lambda =0.13\) and another fit with \({k}_{\max }=6,\lambda =0\).
Data availability
The extracted cross-section data reported in this study are available in the extended figures and tables section, and will be made available for public access on https://minerva.fnal.gov/data-release-page/.
Code availability
The MINERvA Collaboration develops and maintains the code used for the simulation of the experimental apparatus and statistical analysis of the raw data used in this result. This code is shared among the collaboration, but not publicly distributed. Enquiries regarding the algorithms and methods used in this result may be directed to the corresponding author.
References
Chambers, E. E. & Hofstadter, R. Structure of the proton. Phys. Rev. 103, 1454–1463 (1956).
Ernst, F. J., Sachs, R. G. & Wali, K. C. Electromagnetic form factors of the nucleon. Phys. Rev. 119, 1105–1114 (1960).
Mann, W. A. et al. Study of the reaction ν + n → μ− + p. Phys. Rev. Lett. 31, 844–847 (1973).
Barish, S. J. et al. Study of neutrino interactions in hydrogen and deuterium. Description of the experiment and study of the reaction neutrino ν + d → μ− + p + ps. Phys. Rev. D 16, 3103 (1977).
Baker, N. J. et al. Quasielastic neutrino scattering: a measurement of the weak nucleon axial vector form-factor. Phys. Rev. D 23, 2499–2505 (1981).
Miller, K. L. et al. Study of the reaction νμd → μ−pps. Phys. Rev. D 26, 537–542 (1982).
Kitagaki, T. et al. High-energy quasielastic νμn → μ−p scattering in deuterium. Phys. Rev. D 28, 436–442 (1983).
Allasia, D. et al. Investigation of exclusive channels in neutrino/anti-neutrino deuteron charged current interactions. Nucl. Phys. B 343, 285–309 (1990).
Kitagaki, T. et al. Study of νd → μ−pps and νd → μ−Δ++(1232)ns using the BNL 7-foot deuterium filled bubble chamber. Phys. Rev. D 42, 1331–1338 (1990).
Singh, S. K. The effect of final state interactions and deuteron binding in νd → μ−pp. Nucl. Phys. B 36, 419–435 (1972).
Aliaga, L. et al. Design, calibration, and performance of the MINERνA detector. Nucl. Instrum. Meth. Phys. Res. A 743, 130–159 (2014).
Kronfeld, A. S. et al. Lattice QCD and neutrino-nucleus scattering. Eur. Phys. J. A 55, 196 (2019).
Alexandrou, C. et al. Nucleon axial and pseudoscalar form factors from lattice QCD at the physical point. Phys. Rev. D 103, 034509 (2021).
Lin, H.-W. Nucleon helicity generalized parton distribution at physical pion mass from lattice QCD. Phys. Lett. B 824, 136821 (2022).
Meyer, A. S., Walker-Loud, A. & Wilkinson, C. Status of lattice QCD determination of nucleon form factors and their relevance for the few-GeV neutrino program. Annu. Rev. Nucl. Part. Sci. 72, 205–232 (2022).
Ayres, D. S. et al. The NOvA Technical Design Report https://doi.org/10.2172/935497 (2007).
Abe, K. et al. The T2K experiment. Nucl. Instrum. Meth. Phys. Res. A 659, 106–135 (2011).
Abi, B. et al. Deep Underground Neutrino Experiment (DUNE) Far Detector Technical Design Report, Volume II: DUNE Physics. Preprint at https://arxiv.org/abs/2002.03005 (2020)
H.-K. Proto-Collaboration et al. Hyper-Kamiokande Design Report. Preprint at https://arxiv.org/abs/1805.04163 (2018).
An, F. et al. Neutrino physics with JUNO. J. Phys. G 43, 030401 (2016).
Miller, G. A. Defining the proton radius: a unified treatment. Phys. Rev. C 99, 035202 (2019).
Xiong, W. et al. A small proton charge radius from an electron–proton scattering experiment. Nature 575, 147–150 (2019).
Atac, H., Constantinou, M., Meziani, Z.-E., Paolone, M. & Sparveris, N. Measurement of the neutron charge radius and the role of its constituents. Nat. Commun. 12, 1759 (2021).
Meyer, A. S., Betancourt, M., Gran, R. & Hill, R. J. Deuterium target data for precision neutrino-nucleus cross sections. Phys. Rev. D 93, 113015 (2016).
Bodek, A., Avvakumov, S., Bradford, R. & Budd, H. S. Vector and axial nucleon form factors: a duality constrained parameterization. Eur. Phys. J. C 53, 349–354 (2008).
Megias, G. D., Bolognesi, S., Barbaro, M. B. & Tomasi-Gustafsson, E. New evaluation of the axial nucleon form factor from electron- and neutrino-scattering data and impact on neutrino-nucleus cross sections. Phys. Rev. C 101, 025501 (2020).
Bernard, V., Elouadrhiri, L. & Meißner, U.-G. Axial structure of the nucleon. J. Phys. G 28, R1 (2001).
Andreev, V. A. et al. Measurement of muon capture on the proton to 1% precision and determination of the pseudoscalar coupling gP. Phys. Rev. Lett. 110, 012504 (2013).
Fanourakis, G. et al. Study of low-energy anti-neutrino interactions on protons. Phys. Rev. D 21, 562 (1980).
Llewellyn Smith, C. H. Neutrino reactions at accelerator energies. Phys. Rep. 3, 261–379 (1972).
Bradford, R., Bodek, A., Budd, H. S. & Arrington, J. A new parameterization of the nucleon elastic form factors. Nucl. Phys. B Proc. Suppl. 159, 127–132 (2006).
Borah, K., Hill, R. J., Lee, G. & Tomalak, O. Parametrization and applications of the low-Q2 nucleon vector form factors. Phys. Rev. D 102, 074012 (2020).
Goldberger, M. L. & Treiman, S. B. Decay of the pi meson. Phys. Rev. 110, 1178–1184 (1958).
Adamson, P. et al. The NuMI neutrino beam. Nucl. Instrum. Meth. Phys. Res. A 806, 279–306 (2016).
Olive, K. A. Review of particle physics. Chin. Phys. C 40, 100001 (2016).
Michael, D. G. et al. The magnetized steel and scintillator calorimeters of the MINOS experiment. Nucl. Instrum. Meth. Phys. Res. A 596, 190–228 (2008).
Elkins, M. et al. Neutron measurements from antineutrino hydrocarbon reactions. Phys. Rev. D 100, 052002 (2019).
Schiavilla, R., Benhar, O., Kievsky, A., Marcucci, L. E. & Viviani, M. Polarization transfer in \({}^{4}{\rm{H}}{\rm{e}}{(\overrightarrow{e},{e}^{{\prime} }\overrightarrow{p})}^{3}{\rm{H}}\): is the ratio GEp/GMp modified in the nuclear medium? Phys. Rev. Lett. 94, 072303 (2005).
Strauch, S. et al. Polarization transfer in the \({}^{4}{\rm{H}}{\rm{e}}{(\overrightarrow{e},{e}^{{\prime} }\overrightarrow{p})}^{3}{\rm{H}}\) reaction up to Q2 = 2.6(GeV/c)2. Phys. Rev. Lett. 91, 052301 (2003).
Thomas, A. W. Role of Quarks in Nuclear Structure (Oxford Univ. Press, 2020); https://oxfordre.com/physics/view/10.1093/acrefore/9780190871994.001.0001/acrefore-9780190871994-e-1.
Agostinelli, S. et al. Geant4–a simulation toolkit. Nucl. Instrum. Meth. Phys. Res. A 506, 250–303 (2003).
Andreopoulos, C. The GENIE neutrino Monte Carlo generator. Acta Phys. Polon. B 40, 2461–2475 (2009).
Rodrigues, P. A., Demgen, J. & Miltenberger, E. et al. Identification of nuclear effects in neutrino-carbon interactions at low three-momentum transfer. Phys. Rev. Lett. 116, 071802 (2016).
Nieves, J., Ruiz Simo, I. & Vicente Vacas, M. J. Inclusive charged-current neutrino-nucleus reactions. Phys. Rev. C 83, 045501 (2011).
Wilkinson, C., Rodrigues, P., Cartwright, S., Thompson, L. & McFarland, K. Reanalysis of bubble chamber measurements of muon-neutrino induced single pion production. Phys. Rev. D 90, 112017 (2014).
Rodrigues, P., Wilkinson, C. & McFarland, K. Constraining the GENIE model of neutrino-induced single pion production using reanalyzed bubble chamber data. Eur. Phys. J. C 76, 474 (2016).
Nieves, J., Amaro, J. E. & Valverde, M. Inclusive quasi-elastic neutrino reactions. Phys. Rev. C 70, 055503 (2004).
Stowell, P. et al. Tuning the GENIE pion production model with MINERvA data. Phy. Rev. D 100, 072005 (2019).
Juszczak, C., Nowak, J. A. & Sobczyk, J. T. Simulations from a new neutrino event generator. Nucl. Phys. B 159, 211–216 (2006).
Ankowski, A. M. & Sobczyk, J. T. Construction of spectral functions for medium-mass nuclei. Phys. Rev. C 77, 044311 (2008).
Benhar, O., Farina, N., Nakamura, H., Sakuda, M. & Seki, R. Electron- and neutrino-nucleus scattering in the impulse approximation regime. Phys. Rev. D 72, 053005 (2005).
D’Agostini, G. A multidimensional unfolding method based on Bayes’ theorem. Nucl. Instrum. Meth. Phys. Res. A 362, 487–498 (1995).
D’Agostini, G. Improved iterative Bayesian unfolding. Preprint at https://doi.org/10.48550/arXiv.1010.0632 (2010).
Valencia, E. et al. Constraint of the MINERνA medium energy neutrino flux using neutrino-electron elastic scattering. Phys. Rev. D 100, 092001 (2019).
Ruterbories, D. et al. Constraining the NuMI neutrino flux using inverse muon decay reactions in MINERvA. Phys. Rev. D 104, 092010 (2021).
Zazueta, L. et al. Improved constraint on the MINERvA medium energy neutrino flux using \(\bar{\nu }{e}^{-}\to \bar{\nu }{e}^{-}\) data. Preprint at https://doi.org/10.48550/arXiv.2209.05540 (2022).
Abe, K. et al. Constraint on the matter–antimatter symmetry-violating phase in neutrino oscillations. Nature 580, 339–344 (2020).
Hayato, Y. A neutrino interaction simulation program library NEUT. Acta Phys. Polon. B 40, 2477–2489 (2009).
Hill, R.J. & Paz, G. Model-independent extraction of the proton charge radius from electron scattering. Phys. Rev. D https://doi.org/10.1103/physrevd.82.113005 (2010).
Aliaga, L. et al. Neutrino flux predictions for the NuMI beam. Phys. Rev. D 94, 092005 (2016).
Bashyal, A. et al. Use of neutrino scattering events with low hadronic recoil to inform neutrino flux and detector energy scale. J. Instrum. 16, P08068 (2021).
Bodek, A. & Ritchie, J. L. Fermi-motion effects in deep-inelastic lepton scattering from nuclear targets. Phys. Rev. D 23, 1070–1091 (1981).
Benhar, O., Day, D. & Sick, I. An archive for quasi-elastic electron-nucleus scattering data. Preprint at https://arxiv.org/abs/nucl-ex/0603032 (2006).
Dytman, S. Final state interaction models in neutrino nucleus cross sections. AIP Conf. Proc. 1189, 51–59 (2009).
Harewood, L. & Gran, R. Elastic hadron-nucleus scattering in neutrino-nucleus reactions and transverse kinematics measurements. Preprint at https://doi.org/10.48550/arXiv.1906.10576 (2019).
Rein, D. & Sehgal, L. M. Neutrino-excitation of baryon resonances and single pion production. Ann. Phys. 133, 79–153 (1981).
Bodek, A. & Yang, U. K. Modeling deep inelastic cross-sections in the few GeV region. Nucl. Phys. Proc. Suppl. 112, 70–76 (2002).
Brown, D. A. et al. ENDF/B-VIII.0: the 8th major release of the nuclear reaction data library with CIELO-project cross sections, new standards and thermal scattering data. Nucl. Data Sheets 148, 1–142 (2018) .
Del Guerra, A. A compilation of n-p and n-C cross sections and their use in a Monte Carlo program to calculate the neutron detection efficiency in plastic scintillator in the energy range 1–300 MeV. Nucl. Instrum. Meth. 135, 337–352 (1976).
Kohley, Z. et al. Modeling interactions of intermediate-energy neutrons in a plastic scintillator array with GEANT4. Nucl. Instrum. Meth. Phys. Res. A 682, 59–65 (2012).
Lepage, G. P. & Brodsky, S. J. Exclusive processes in perturbative quantum chromodynamics. Phys. Rev. D 22, 2157 (1980).
Chernyak, V. L. & Zhitnitsky, A. R. Asymptotic behavior of exclusive processes in QCD. Phys. Rep. 112, 173–318 (1984).
Hansen, P. C. & O’Leary, D. The use of the L-curve in the regularization of discrete ill-posed problems. SIAM J. Sci. Comput. 14, 1487–1503 (1993).
Acknowledgements
All authors are members of the MINERvA Collaboration. This document was prepared by members of the MINERvA Collaboration using the resources of the Fermi National Accelerator Laboratory (Fermilab), a US Department of Energy, Office of Science, HEP User Facility. Fermilab is managed by Fermi Research Alliance, LLC (FRA), acting under contract no. DE-AC02-07CH11359. These resources included support for the MINERvA construction project, and support for construction also was granted by the United States National Science Foundation under award no. PHY-0619727 and by the University of Rochester. Support for participating scientists was provided by NSF and DOE (USA); by NSERC (Canada); by CAPES and CNPq (Brazil); by CoNaCyT (Mexico); by Proyecto Basal FB 0821, CONICYT PIA ACT1413, and Fondecyt 3170845 and 11130133 (Chile); by CONCYTEC (Consejo Nacional de Ciencia, Tecnología e Innovación Tecnológica), DGI-PUCP (Dirección de Gestión de la Investigación-Pontificia Universidad Católica del Peru) and VRI-UNI (Vice-Rectorate for Research of National University of Engineering) (Peru); NCN Opus grant no. 2016/21/B/ST2/01092 (Poland); by Science and Technology Facilities Council (UK); by EU Horizon 2020 Marie Skłodowska-Curie Action. We thank the MINOS Collaboration for use of its near detector data. Finally, we thank the staff of Fermilab for support of the beamline, the detector and computing infrastructure.
Author information
Authors and Affiliations
Contributions
The operation, Monte Carlo simulation and data analysis of the MINERvA Experiment are carried out by the MINERvA Collaboration with contributions from all collaborators listed as authors on this manuscript. The scientific results presented here have been presented to and discussed by the full collaboration, and all authors have approved the final version of the manuscript. T.C. and A.B. first proposed this analysis to the authors. D.R. developed the analysis tools with input from T.C. and A.O. T.C. performed the analysis of data and simulation described in this work. T.C. and M.L.M. performed the fits to the data. T.C. and K.S.M. wrote the manuscript.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature thanks the anonymous reviewers for their contribution to the peer review of this work. Peer reviewer reports are available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data figures and tables
Extended Data Fig. 1 Illustration of MINERvA detector.
A 3D illustration of the MINERvA detector (left), a flat cross-sectional portrayal (reproduced from ref .11, right). The detector comprises 120 hexagonal modules, each consists of an Inner Detector (ID), side Electromagnetic Calorimeter (ECAL), and side Hadronic Calorimeter (HCAL) arranged outwards from the center. Horizontally along the beam direction in the ID region are four sub-detectors: nuclear target detector, active scintillator detector, the ECAL and the HCAL. Each scintillator plane in the detector is arranged in one of “X”, “U”, or “V” (bottom) orientations. A liquid helium tank is located upstream of the nuclear target detector. Downstream of the MINERvA detector is the magnetized MINOS36 near detector acting as MINERvA’s muon spectrometer.
Extended Data Fig. 2 Definition of δθP and δθR, and the event distributions.
(top) Angular variable definitions. The variables δθR and δθP are defined according to the rotated reference frame (\((\hat{x},{\hat{y}}^{{\prime} },{\hat{z}}^{{\prime} })\). Distributions of (bottom left) δθP and (bottom right) δθR in the QELike sample normalized to bin width, after sideband fits. The vertical error bars around the data points represent 1 standard deviation due to statistical uncertainty.
Extended Data Fig. 3 Simulated event rate in δθP − δθR plane in selected Q2 analysis bins.
Heat map showing Monte Carlo event rate for CCE, CCQE, QELike 2p2h, and QELike Resonant interaction models in the δθR-δθP plane at a few slices of Q2. CCE events are concentrated around origin, while the CCQE events have broader spread. Both QELike 2p2h and resonant events show diffused structure going out to larger δθR and δθP regions. The color scales in the heat maps are different because of different event rates for each subsample.
Extended Data Fig. 4 Post-fit event rates and ratios in the validation regions.
The “Non-QE validation” region is shown separated into two sub-regions. “Non-QE Val. 1” spans ∣δθP∣ < 20°, “Non-QE Val. 2” occupies 20° < ∣δθP∣ < 55°. The vertical error bars around the data points and the error band around the model prediction account for 1 standard deviation due to statistical uncertainty. The CCE signal and the regions used in the background fits are shown in Figs. 2 and 3, respectively.
Extended Data Fig. 5 Background constraint method tested on a neutrino sample.
(Top left) Signal and (top right) QE region distributions of the neutrino sample that looks for proton final states after fit. The vertical error bars around the data points and the error band around the model prediction account for 1 standard deviation due to statistical and systematic uncertainties. (Bottom) Application of ± 100% shift in 2p2h (gray band) on δθR for events in a − 10° < δθP < 10° and 0.2 (GeV/c)2 < Q2 < 0.4 (GeV/c)2 slice. The analog to the CCE selection selects between − 10° and 10° in δθR.
Extended Data Fig. 6 Signal event selection efficiency.
CCE signal efficiency as a function of Q2. The inner error bar on each data point accounts for the statistical effect of a Poisson standard deviation, while the full error bar account for all sources of systematic uncertainties.
Extended Data Fig. 7 Fractional uncertainties.
Fractional statistical and systematic uncertainties, as a function of Q2 uncertainties. Systematic uncertainties in the “other” category, including the neutron and proton interaction uncertainties, are shown on the right. The neutron interaction systematic accounts for the neutron secondary interaction uncertainties in detector. The leading interaction channels, such as (nC, Bnp),(nC, 3α), and \((nC,{n}^{{\prime} }C\gamma )\), are assigned 10% to 15% uncertainties below a kinetic energy of 100 MeV.
Extended Data Fig. 8 Cross section and the Q2-dependent flux cut.
(left) Measured cross-section and theory predictions, (right) Regions of neutrino energy and flux in signal selection at each Q2. The inner error bars on the data points account for 1 standard deviation due to statistical uncertainty only, and the full error bars include contribution from all sources of systematic uncertainties.
Supplementary information
Supplementary Information
Supplementary Tables 1–7.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Cai, T., Moore, M.L., Olivier, A. et al. Measurement of the axial vector form factor from antineutrino–proton scattering. Nature 614, 48–53 (2023). https://doi.org/10.1038/s41586-022-05478-3
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41586-022-05478-3
- Springer Nature Limited