This Letter reports new results on muon neutrino disappearance from NOvA, using a 14 kton detector equivalent exposure of 6.05×10^{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 θ_{23} mixing (θ_{23}=π/4). Assuming the normal mass hierarchy, we find Δm_{32}^{2}=(2.67±0.11)×10^{-3}  eV^{2} and sin^{2}θ_{23} at the two statistically degenerate values 0.404_{-0.022}^{+0.030} and 0.624_{-0.030}^{+0.022}, both at the 68% confidence level. Our data disfavor the maximal mixing scenario with 2.6σ significance.

Using two years of data from the NOvA Near Detector at Fermilab, we report a seasonal variation of cosmic ray induced multiple-muon (Nμ≥2) event rates which has an opposite phase to the seasonal variation in the atmospheric temperature. The strength of the seasonal multiple-muon variation is shown to increase as a function of the muon multiplicity. However, no significant dependence of the strength of the seasonal variation of the multiple-muon variation is seen as a function of the muon zenith angle, or the spatial or angular separation between the correlated muons.

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×1020 protons-on-target in the NuMI beam at Fermilab, we observe 95 neutral-current candidates at the Far Detector compared with 83.5±9.7(stat)±9.4(syst) events predicted assuming mixing only occurs between active neutrino species. No evidence for νμ→νs transitions is found. Interpreting these results within a 3+1 model, we place constraints on the mixing angles θ24<20.8° and θ34<31.2° at the 90% C.L. for 0.05 eV2≤Δm412≤0.5 eV2, the range of mass splittings that produce no significant oscillations over the Near Detector baseline.

The NOvA experiment has seen a 4.4σ signal of ν[over ¯]_{e} appearance in a 2 GeV ν[over ¯]_{μ} beam at a distance of 810 km. Using 12.33×10^{20} protons on target delivered to the Fermilab NuMI neutrino beamline, the experiment recorded 27 ν[over ¯]_{μ}→ν[over ¯]_{e} candidates with a background of 10.3 and 102 ν[over ¯]_{μ}→ν[over ¯]_{μ} candidates. This new antineutrino data are combined with neutrino data to measure the parameters |Δm_{32}^{2}|=2.48_{-0.06}^{+0.11}×10^{-3}  eV^{2}/c^{4} and sin^{2}θ_{23} in the ranges from (0.53-0.60) and (0.45-0.48) in the normal neutrino mass hierarchy. The data exclude most values near δ_{CP}=π/2 for the inverted mass hierarchy by more than 3σ and favor the normal neutrino mass hierarchy by 1.9σ and θ_{23} values in the upper octant by 1.6σ.

The two-detector design of the NOvA neutrino oscillation experiment, in which two functionally identical detectors are exposed to an intense neutrino beam, aids in canceling leading order effects of cross-section uncertainties. However, limited knowledge of neutrino interaction cross sections still gives rise to some of the largest systematic uncertainties in current oscillation measurements. We show contemporary models of neutrino interactions to be discrepant with data from NOvA, consistent with discrepancies seen in other experiments. Adjustments to neutrino interaction models in GENIE are presented, creating an effective model that improves agreement with our data. We also describe systematic uncertainties on these models, including uncertainties on multi-nucleon interactions from a newly developed procedure using NOvA near detector data.

Results are reported from an improved measurement of ν_{μ}→ν_{e} transitions by the NOvA experiment. Using an exposure equivalent to 6.05×10^{20} protons on target, 33 ν_{e} candidates are observed with a background of 8.2±0.8 (syst.). Combined with the latest NOvA ν_{μ} disappearance data and external constraints from reactor experiments on sin^{2}2θ_{13}, the hypothesis of inverted mass hierarchy with θ_{23} in the lower octant is disfavored at greater than 93% C.L. for all values of δ_{CP}.

The cross section of neutrino-induced neutral-current coherent π0 production on a carbon-dominated target is measured in the NOvA near detector. This measurement uses a narrow-band neutrino beam with an average neutrino energy of 2.7 GeV, which is of interest to ongoing and future long-baseline neutrino oscillation experiments. The measured, flux-averaged cross section is σ=13.8±0.9(stat)±2.3(syst)×10-40 cm2/nucleus, consistent with model prediction. This result is the most precise measurement of neutral-current coherent π0 production in the few-GeV neutrino energy region.

Using the NOvA neutrino detectors, a broad search has been performed for any signal coincident with 28 gravitational wave events detected by the LIGO/Virgo Collaboration between September 2015 and July 2019. For all of these events, NOvA is sensitive to possible arrival of neutrinos and cosmic rays of GeV and higher energies. For five (seven) events in the NOvA Far (Near) Detector, timely public alerts from the LIGO/Virgo Collaboration allowed recording of MeV-scale events. No signal candidates were found.

Searches for electron antineutrino, muon neutrino, and muon antineutrino disappearance driven by sterile neutrino mixing have been carried out by the Daya Bay and MINOS+ collaborations. This Letter presents the combined results of these searches, along with exclusion results from the Bugey-3 reactor experiment, framed in a minimally extended four-neutrino scenario. Significantly improved constraints on the θ_{μe} mixing angle are derived that constitute the most constraining limits to date over five orders of magnitude in the mass-squared splitting Δm_{41}^{2}, excluding the 90% C.L. sterile-neutrino parameter space allowed by the LSND and MiniBooNE observations at 90% CL_{s} for Δm_{41}^{2}<13  eV^{2}. Furthermore, the LSND and MiniBooNE 99% C.L. allowed regions are excluded at 99% CL_{s} for Δm_{41}^{2}<1.6 eV^{2}.

This paper reports the first measurement using the NOvA detectors of νμ disappearance in a νμ beam. The analysis uses a 14 kton-equivalent exposure of 2.74×1020 protons-on-target from the Fermilab NuMI beam. Assuming the normal neutrino mass hierarchy, we measure Δm322=(2.52-0.18+0.20)×10-3 eV2 and sin2θ23 in the range 0.38-0.65, both at the 68% confidence level, with two statistically degenerate best-fit points at sin2θ23=0.43 and 0.60. Results for the inverted mass hierarchy are also presented.