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}$.

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.

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.