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

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

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.

We report a two-detector measurement of the propagation speed of neutrinos over a baseline of 734 km. The measurement was made with the NuMI beam at Fermilab between the near and far MINOS detectors. The fractional difference between the neutrino speed and the speed of light is determined to be $(v/c-1) = (1.0 \pm 1.1) \times 10^{-6}$, consistent with relativistic neutrinos.

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.