A search for mixing between active neutrinos and light sterile neutrinos has been performed by looking for muon neutrino disappearance in two detectors at baselines of 1.04 km and 735 km, using a combined MINOS and MINOS+ exposure of $16.36\times10^{20}$ protons-on-target. A simultaneous fit to the charged-current muon neutrino and neutral-current neutrino energy spectra in the two detectors yields no evidence for sterile neutrino mixing using a 3+1 model. The most stringent limit to date is set on the mixing parameter $\sin^2\theta_{24}$ for most values of the sterile neutrino mass-splitting $\Delta m^2_{41} > 10^{-4}$ eV$^2$.

We report results from the first search for sterile neutrinos mixing with active neutrinos through a reduction in the rate of neutral-current interactions over a baseline of 810\,km between the NOvA detectors. Analyzing a 14-kton detector equivalent exposure of 6.05$\times$10$^{20}$ protons-on-target in the NuMI beam at Fermilab, we observe 95 neutral-current candidates at the Far Detector compared with $83.5 \pm 9.7 \mbox{(stat.)} \pm 9.4 \mbox{(syst.)}$ events predicted assuming mixing only occurs between active neutrino species. No evidence for $\nu_{\mu} \rightarrow \nu_{s}$ transitions is found. Interpreting these results within a 3+1 model, we place constraints on the mixing angles $\theta_{24}<20.8^{\circ}$ and $\theta_{34}<31.2^{\circ}$ at the 90% C.L. for $0.05~eV^2\leq \Delta m^2_{41}\leq 0.5~eV^2$, the range of mass splittings that produce no significant oscillations over the Near Detector baseline.

Results are reported from an improved measurement of $\nu_\mu \rightarrow \nu_e$ transitions by the NOvA experiment. Using an exposure equivalent to $6.05\times10^{20}$ protons-on-target 33 $\nu_e$ candidates were observed with a background of $8.2\pm0.8$ (syst.). Combined with the latest NOvA $\nu_\mu$ disappearance data and external constraints from reactor experiments on $\sin^22\theta_{13}$, the hypothesis of inverted mass hierarchy with $\theta_{23}$ in the lower octant is disfavored at greater than $93\%$ C.L. for all values of $\delta_{CP}$.

This Letter reports new results on muon neutrino disappearance from NOvA, using a 14 kton detector equivalent exposure of $6.05\times10^{20}$ protons-on-target from the NuMI beam at the Fermi National Accelerator Laboratory. The measurement probes the muon-tau symmetry hypothesis that requires maximal mixing ($\theta_{23} = \pi/4$). Assuming the normal mass hierarchy, we find $\Delta m^2 = (2.67 \pm 0.11)\times 10^{-3}$ eV$^2$ and $\sin^2 \theta_{23}$ at the two statistically degenerate values $0.404^{+0.030}_{-0.022}$ and $0.624^{+0.022}_{-0.030}$, both at the 68% confidence level. Our data disfavor the maximal mixing scenario with 2.6 $\sigma$ significance.

We report new constraints on the size of large extra dimensions from data collected by the MINOS experiment between 2005 and 2012. Our analysis employs a model in which sterile neutrinos arise as Kaluza-Klein states in large extra dimensions and thus modify the neutrino oscillation probabilities due to mixing between active and sterile neutrino states. Using Fermilab's NuMI beam exposure of $10.56 \times 10^{20}$ protons-on-target, we combine muon neutrino charged current and neutral current data sets from the Near and Far Detectors and observe no evidence for deviations from standard three-flavor neutrino oscillations. The ratios of reconstructed energy spectra in the two detectors constrain the size of large extra dimensions to be smaller than $0.45\,\mu\text{m}$ at 90% C.L. in the limit of a vanishing lightest active neutrino mass. Stronger limits are obtained for non-vanishing masses.

We report new constraints on flavor-changing non-standard neutrino interactions from the MINOS long-baseline experiment using $\nu_{e}$ and $\bar{\nu}_{e}$ appearance candidate events from predominantly $\nu_{\mu}$ and $\bar{\nu}_{\mu}$ beams. We used a statistical selection algorithm to separate $\nu_{e}$ candidates from background events, enabling an analysis of the combined MINOS neutrino and antineutrino data. We observe no deviations from standard neutrino mixing, and thus place constraints on the non-standard interaction matter effect, $|\varepsilon_{e\tau}|$, and phase, $(\delta_{CP}+\delta_{e\tau})$, using a thirty-bin likelihood fit.

Data from the MINOS experiment has been used to search for mixing between muon neutrinos and muon antineutrinos using a time-independent Lorentz-violating formalism derived from the Standard-Model Extension (SME). MINOS is uniquely capable of searching for muon neutrino-antineutrino mixing given its long baseline and ability to distinguish between neutrinos and antineutrinos on an event-by-event basis. Neutrino and antineutrino interactions were observed in the MINOS Near and Far Detectors from an exposure of 10.56$\times10^{20}$ protons-on-target from the NuMI neutrino-optimized beam. No evidence was found for such transitions and new, highly stringent limits were placed on the SME coefficients governing them. We place the first limits on the SME parameters $(c_{L})^{TT}_{\mu\mu} $ and $(c_{L})^{TT}_{\tau\tau}$ at $-8.4\times10^{-23} < (c_{L})^{TT}_{\mu\mu} < 8.0\times10^{-23}$ and $-8.0\times10^{-23} < (c_{L})^{TT}_{\tau\tau} < 8.4\times10^{-23}$, and the world's best limits on the $\tilde{g}^{ZT}_{\mu\overline{\mu}}$ and $\tilde{g}^{ZT}_{\tau\overline{\tau}}$ parameters at $|\tilde{g}^{ZT}_{\mu\overline{\mu}}| < 3.3\times 10^{-23}$ and $|\tilde{g}^{ZT}_{\tau\overline{\tau}}| < 3.3\times 10^{-23}$, all limits quoted at $3\sigma$.

Searches for a light sterile neutrino have been performed independently by the MINOS and the Daya Bay experiments using the muon (anti)neutrino and electron antineutrino disappearance channels, respectively. In this Letter, results from both experiments are combined with those from the Bugey-3 reactor neutrino experiment to constrain oscillations into light sterile neutrinos. The three experiments are sensitive to complementary regions of parameter space, enabling the combined analysis to probe regions allowed by the LSND and MiniBooNE experiments in a minimally extended four-neutrino flavor framework. Stringent limits on $\sin^2 2\theta_{\mu e}$ are set over 6 orders of magnitude in the sterile mass-squared splitting $\Delta m^2_{41}$. The sterile-neutrino mixing phase space allowed by the LSND and MiniBooNE experiments is excluded for $\Delta m^2_{41} < 0.8$ eV$^2$ at 95% CL$_s$.

Forward single $\pi^0$ production by coherent neutral-current interactions, $\nu \mathcal{A} \to \nu \mathcal{A} \pi^0$, is investigated using a 2.8$\times 10^{20}$ protons-on-target exposure of the MINOS Near Detector. For single-shower topologies, the event distribution in production angle exhibits a clear excess above the estimated background at very forward angles for visible energy in the range~1-8 GeV. Cross sections are obtained for the detector medium comprised of 80% iron and 20% carbon nuclei with $\langle \mathcal{A} \rangle = 48$, the highest-$\langle \mathcal{A} \rangle$ target used to date in the study of this coherent reaction. The total cross section for coherent neutral-current single-$\pi^0$ production initiated by the $\nu_\mu$ flux of the NuMI low-energy beam with mean (mode) $E_{\nu}$ of 4.9 GeV (3.0 GeV), is $77.6\pm5.0\,(\text{stat}) ^{+15.0}_{-16.8}\,(\text{syst})\times10^{-40}\,\text{cm}^2~\text{per nucleus}$. The results are in good agreement with predictions of the Berger-Sehgal model.

We report results of a search for oscillations involving a light sterile neutrino over distances of 1.04 and $735\,\mathrm{km}$ in a $\nu_{\mu}$-dominated beam with a peak energy of $3\,\mathrm{GeV}$. The data, from an exposure of $10.56\times 10^{20}\,\textrm{protons on target}$, are analyzed using a phenomenological model with one sterile neutrino. We constrain the mixing parameters $\theta_{24}$ and $\Delta m^{2}_{41}$ and set limits on parameters of the four-dimensional Pontecorvo-Maki-Nakagawa-Sakata matrix, $|U_{\mu 4}|^{2}$ and $|U_{\tau 4}|^{2}$, under the assumption that mixing between $\nu_{e}$ and $\nu_{s}$ is negligible ($|U_{e4}|^{2}=0$). No evidence for $\nu_{\mu} \to \nu_{s}$ transitions is found and we set a world-leading limit on $\theta_{24}$ for values of $\Delta m^{2}_{41} \lesssim 1\,\mathrm{eV}^{2}$.

Convolutional neural networks (CNNs) have been widely applied in the computer vision community to solve complex problems in image recognition and analysis. We describe an application of the CNN technology to the problem of identifying particle interactions in sampling calorimeters used commonly in high energy physics and high energy neutrino physics in particular. Following a discussion of the core concepts of CNNs and recent innovations in CNN architectures related to the field of deep learning, we outline a specific application to the NOvA neutrino detector. This algorithm, CVN (Convolutional Visual Network) identifies neutrino interactions based on their topology without the need for detailed reconstruction and outperforms algorithms currently in use by the NOvA collaboration.

We report results from the first search for $\nu_\mu\to\nu_e$ transitions by the NOvA experiment. In an exposure equivalent to $2.74\times10^{20}$ protons-on-target in the upgraded NuMI beam at Fermilab, we observe 6 events in the Far Detector, compared to a background expectation of $0.99\pm0.11$ (syst.) events based on the Near Detector measurement. A secondary analysis observes 11 events with a background of $1.07\pm0.14$ (syst.). The $3.3\sigma$ excess of events observed in the primary analysis disfavors $0.1\pi < \delta_{CP} < 0.5\pi$ in the inverted mass hierarchy at the 90% C.L.

The charge ratio, $R_\mu = N_{\mu^+}/N_{\mu^-}$, for cosmogenic multiple-muon events observed at an under- ground depth of 2070 mwe has been measured using the magnetized MINOS Far Detector. The multiple-muon events, recorded nearly continuously from August 2003 until April 2012, comprise two independent data sets imaged with opposite magnetic field polarities, the comparison of which allows the systematic uncertainties of the measurement to be minimized. The multiple-muon charge ratio is determined to be $R_\mu = 1.104 \pm 0.006 {\rm \,(stat.)} ^{+0.009}_{-0.010} {\rm \,(syst.)} $. This measurement complements previous determinations of single-muon and multiple-muon charge ratios at underground sites and serves to constrain models of cosmic ray interactions at TeV energies.

This paper reports the first measurement using the NOvA detectors of $\nu_\mu$ disappearance in a $\nu_\mu$ beam. The analysis uses a 14 kton-equivalent exposure of $2.74 \times 10^{20}$ protons-on-target from the Fermilab NuMI beam. Assuming the normal neutrino mass hierarchy, we measure $\Delta m^{2}_{32}=(2.52^{+0.20}_{-0.18})\times 10^{-3}$ eV$^{2}$ and $\sin^2\theta_{23}$ in the range 0.38-0.65, both at the 68% confidence level, with two statistically-degenerate best fit points at $\sin^2\theta_{23} = $ 0.43 and 0.60. Results for the inverted mass hierarchy are also presented.

This paper describes the hardware and operations of the Neutrinos at the Main Injector (NuMI) beam at Fermilab. It elaborates on the design considerations for the beam as a whole and for individual elements. The most important design details of individual components are described. Beam monitoring systems and procedures, including the tuning and alignment of the beam and NuMI long-term performance, are also discussed.

We report the first observation of seasonal modulations in the rates of cosmic ray multiple-muon events at two underground sites, the MINOS Near Detector with an overburden of 225 mwe, and the MINOS Far Detector site at 2100 mwe. At the deeper site, multiple-muon events with muons separated by more than 8 m exhibit a seasonal rate that peaks during the summer, similar to that of single-muon events. In contrast and unexpectedly, the rate of multiple-muon events with muons separated by less than 5-8 m, and the rate of multiple-muon events in the smaller, shallower Near Detector, exhibit a seasonal rate modulation that peaks in the winter.

Kinematic distributions from an inclusive sample of 1.41 x 10^6 charged-current nu_mu interactions on iron, obtained using the MINOS Near Detector exposed to a wide-band beam with peak flux at 3 GeV, are compared to a conventional treatment of neutrino scattering within a Fermi gas nucleus. Results are used to guide the selection of a subsample enriched in quasielastic nu_mu Fe interactions, containing an estimated 123,000 quasielastic events of incident energies 1 < E_nu < 8 GeV, with <E_nu> = 2.79 GeV. Four additional subsamples representing topological and kinematic sideband regions to quasielastic scattering are also selected for the purpose of evaluating backgrounds. Comparisons using subsample distributions in four-momentum transfer Q^2 show the Monte Carlo model to be inadequate at low Q^2. Its shortcomings are remedied via inclusion of a Q^2-dependent suppression function for baryon resonance production, developed from the data. A chi-square fit of the resulting Monte Carlo simulation to the shape of the Q^2 distribution for the quasielastic-enriched sample is carried out with the axial-vector mass M_A of the dipole axial-vector form factor of the neutron as a free parameter. The effective M_A which best describes the data is 1.23 +0.13/-0.09 (fit) +0.12/-0.15 (syst.) GeV.

A sample of 1.53$\times$10$^{9}$ cosmic-ray-induced single muon events has been recorded at 225 meters-water-equivalent using the MINOS Near Detector. The underground muon rate is observed to be highly correlated with the effective atmospheric temperature. The coefficient $\alpha_{T}$, relating the change in the muon rate to the change in the vertical effective temperature, is determined to be 0.428$\pm$0.003(stat.)$\pm$0.059(syst.). An alternative description is provided by the weighted effective temperature, introduced to account for the differences in the temperature profile and muon flux as a function of zenith angle. Using the latter estimation of temperature, the coefficient is determined to be 0.352$\pm$0.003(stat.)$\pm$0.046(syst.).

We report on a new analysis of neutrino oscillations in MINOS using the complete set of accelerator and atmospheric data. The analysis combines the $\nu_{\mu}$ disappearance and $\nu_{e}$ appearance data using the three-flavor formalism. We measure $|\Delta m^{2}_{32}|=[2.28-2.46]\times10^{-3}\mbox{\,eV}^{2}$ (68% C.L.) and $\sin^{2}\theta_{23}=0.35-0.65$ (90% C.L.) in the normal hierarchy, and $|\Delta m^{2}_{32}|=[2.32-2.53]\times10^{-3}\mbox{\,eV}^{2}$ (68% C.L.) and $\sin^{2}\theta_{23}=0.34-0.67$ (90% C.L.) in the inverted hierarchy. The data also constrain $\delta_{CP}$, the $\theta_{23}$ octant degeneracy and the mass hierarchy; we disfavor 36% (11%) of this three-parameter space at 68% (90%) C.L.

The fixed-target MIPP experiment, Fermilab E907, was designed to measure the production of hadrons from the collisions of hadrons of momenta ranging from 5 to 120 GeV/c on a variety of nuclei. These data will generally improve the simulation of particle detectors and predictions of particle beam fluxes at accelerators. The spectrometer momentum resolution is between 3 and 4%, and particle identification is performed for particles ranging between 0.3 and 80 GeV/c using $dE/dx$, time-of-flight and Cherenkov radiation measurements. MIPP collected $1.42 \times10^6$ events of 120 GeV Main Injector protons striking a target used in the NuMI facility at Fermilab. The data have been analyzed and we present here charged pion yields per proton-on-target determined in bins of longitudinal and transverse momentum between 0.5 and 80 GeV/c, with combined statistical and systematic relative uncertainties between 5 and 10%.

We report measurements of oscillation parameters from $\nu_{mu}$ and $\bar{\nu}_{\mu}$ disappearance using beam and atmospheric data from MINOS. The data comprise exposures of \unit[$10.71 \times 10^{20}$]{protons on target (POT)} in the $\nu_{\mu}$-dominated beam, $\unit[3.36\times10^{20}]{POT}}$ in the $\bar{\nu}_{\mu}$-enhanced beam, and 37.88 kton-years of atmospheric neutrinos. Assuming identical $\nu$ and $\bar{\nu}$ oscillation parameters, we measure \mbox{$|\Delta m^2}| = \unit[2.41^{+0.09}_{-0.10}) \times 10^{-3}]{eV^{2}}$} and $\sin^{2}/!/left(2\theta \right) = 0.950^{+0.035}_{-0.036}$. Allowing independent $\nu$ and $\bar{\nu}$ oscillations, we measure antineutrino parameters of $|\Delta \bar{m}^2| = \unit[(2.50 ^{+0.23}_{-0.25}) \times 10^{-3}]{eV^{2}}$ and $\sin^{2}/!/left(2\bar{\theta} \right) = 0.97^{+0.03}_{-0.08}$, with minimal change to the neutrino parameters.

We report on $\nu_e$ and $\bar{\nu}_e$ appearance in $\nu_\mu$ and $\bar{\nu}_\mu$ beams using the full MINOS data sample. The comparison of these $\nu_e$ and $\bar{\nu}_e$ appearance data at a 735 km baseline with $\theta_{13}$ measurements by reactor experiments probes $\delta$, the $\theta_{23}$ octant degeneracy, and the mass hierarchy. This analysis is the first use of this technique and includes the first accelerator long-baseline search for $\bar{\nu}_\mu\rightarrow\bar{\nu}_e$. Our data disfavor 31% (5%) of the three-parameter space defined by $\delta$, the octant of the $\theta_{23}$, and the mass hierarchy at the 68% (90%) C.L. We measure a value of 2sin$^2(2\theta_{13})$sin$^2(\theta_{23})$ that is consistent with reactor experiments.

We report new constraints on flavor-changing non-standard neutrino interactions (NSI) using data from the MINOS experiment. We analyzed a combined set of beam neutrino and antineutrino data, and found no evidence for deviations from standard neutrino mixing. The observed energy spectra constrain the NSI parameter to the range $-0.20 < \varepsilon_{\mu\tau} < 0.07\;\text{(90% C.L.)}$

This paper reports measurements of atmospheric neutrino and antineutrino interactions in the MINOS Far Detector, based on 2553 live-days (37.9 kton-years) of data. A total of 2072 candidate events are observed. These are separated into 905 contained-vertex muons and 466 neutrino-induced rock-muons, both produced by charged-current $\nu_{\mu}$ and $\bar{\nu}_{\mu}$ interactions, and 701 contained-vertex showers, composed mainly of charged-current $\nu_{e}$ and $\bar{\nu}_{e}$ interactions and neutral-current interactions. The curvature of muon tracks in the magnetic field of the MINOS Far Detector is used to select separate samples of $\nu_{\mu}$ and $\bar{\nu}_{\mu}$ events. The observed ratio of $\bar{\nu}_{\mu}$ to $\nu_{\mu}$ events is compared with the Monte Carlo simulation, giving a double ratio of $R^{data}_{\bar{\nu}/\nu}/R^{MC}_{\bar{\nu}/\nu} = 1.03 \pm 0.08 (stat.) \pm 0.08 (syst.)$. The $\nu_{\mu}$ and $\bar{\nu}_{\mu}$ data are separated into bins of $L/E$ resolution, based on the reconstructed energy and direction of each event, and a maximum likelihood fit to the observed $L/E$ distributions is used to determine the atmospheric neutrino oscillation parameters. This fit returns 90% confidence limits of $|\Delta m^{2}| = (1.9 \pm 0.4) \times 10^{-3} eV^{2}$ and $sin^{2} 2\theta > 0.86$. The fit is extended to incorporate separate $\nu_{\mu}$ and $\bar{\nu}_{\mu}$ oscillation parameters, returning 90% confidence limits of $|\Delta m^{2}|-|\Delta \bar{m}^{2}| = 0.6^{+2.4}_{-0.8} \times 10^{-3} eV^{2}$ on the difference between the squared-mass splittings for neutrinos and antineutrinos.

We report an improved measurement of muon anti-neutrino disappearance over a distance of 735km using the MINOS detectors and the Fermilab Main Injector neutrino beam in a muon anti-neutrino enhanced configuration. From a total exposure of 2.95e20 protons on target, of which 42% have not been previously analyzed, we make the most precise measurement of the anti-neutrino "atmospheric" delta-m squared = 2.62 +0.31/-0.28 (stat.) +/- 0.09 (syst.) and constrain the anti-neutrino atmospheric mixing angle >0.75 (90%CL). These values are in agreement with those measured for muon neutrinos, removing the tension reported previously.