A Neutrino Factory where neutrinos of all species are produced in equal quantities by muon decay is described as a facility at the intensity frontier for exquisite precision providing ideal conditions for ultimate neutrino studies and the ideal complement to Long Baseline Facilities like LBNF at Fermilab. It is foreseen to be built in stages with progressively increasing complexity and performance, taking advantage of existing or proposed facilities at an existing laboratory like Fermilab. A tentative layout based on a recirculating linac providing opportunities for considerable saving is discussed as well as its possible evolution toward a muon collider if and when requested by Physics. Tentative parameters of the various stages are presented as well as the necessary R&D to address the technological issues and demonstrate their feasibility.

Neutrons produced by cosmic ray muons are an important background for underground experiments studying neutrino oscillations, neutrinoless double beta decay, dark matter, and other rare-event signals. A measurement of the neutron yield in the three different experimental halls of the Daya Bay Reactor Neutrino Experiment at varying depth is reported. The neutron yield in Daya Bay's liquid scintillator is measured to be $Y_n=(10.26\pm 0.86)\times 10^{-5}$, $(10.22\pm 0.87)\times 10^{-5}$, and $(17.03\pm 1.22)\times 10^{-5}~\mu^{-1}~$g$^{-1}~$cm$^2$ at depths of 250, 265, and 860 meters-water-equivalent. These results are compared to other measurements and the simulated neutron yield in Fluka and Geant4. A global fit including the Daya Bay measurements yields a power law coefficient of $0.77 \pm 0.03$ for the dependence of the neutron yield on muon energy.

The influence of a supersonic projectile on a three-dimensional complex plasma is studied. Micron sized particles in a low-temperature plasma formed a large undisturbed system in the new 'Zyflex' chamber during microgravity conditions. A supersonic probe particle excited a Mach cone with Mach number M $\approx$ 1.5 - 2 and double Mach cone structure in the large weakly damped particle cloud. The speed of sound is measured with different methods and particle charge estimations are compared to calculations from standard theories. The high image resolution enables the study of Mach cones in microgravity on the single particle level of a three-dimensional complex plasma and gives insight to the dynamics. A heating of the microparticles is discovered behind the supersonic projectile but not in the flanks of the Mach cone.

The Daya Bay Experiment consists of eight identically designed detectors located in three underground experimental halls named as EH1, EH2, EH3, with 250, 265 and 860 meters of water equivalent vertical overburden, respectively. Cosmic muon events have been recorded over a two-year period. The underground muon rate is observed to be positively correlated with the effective atmospheric temperature and to follow a seasonal modulation pattern. The correlation coefficient $\alpha$, describing how a variation in the muon rate relates to a variation in the effective atmospheric temperature, is found to be $\alpha_{\text{EH1}} = 0.362\pm0.031$, $\alpha_{\text{EH2}} = 0.433\pm0.038$ and $\alpha_{\text{EH3}} = 0.641\pm0.057$ for each experimental hall.

The Daya Bay experiment has observed correlations between reactor core fuel evolution and changes in the reactor antineutrino flux and energy spectrum. Four antineutrino detectors in two experimental halls were used to identify 2.2 million inverse beta decays (IBDs) over 1230 days spanning multiple fuel cycles for each of six 2.9 GW$_{\textrm{th}}$ reactor cores at the Daya Bay and Ling Ao nuclear power plants. Using detector data spanning effective $^{239}$Pu fission fractions, $F_{239}$, from 0.25 to 0.35, Daya Bay measures an average IBD yield, $\bar{\sigma}_f$, of $(5.90 \pm 0.13) \times 10^{-43}$ cm$^2$/fission and a fuel-dependent variation in the IBD yield, $d\sigma_f/dF_{239}$, of $(-1.86 \pm 0.18) \times 10^{-43}$ cm$^2$/fission. This observation rejects the hypothesis of a constant antineutrino flux as a function of the $^{239}$Pu fission fraction at 10 standard deviations. The variation in IBD yield was found to be energy-dependent, rejecting the hypothesis of a constant antineutrino energy spectrum at 5.1 standard deviations. While measurements of the evolution in the IBD spectrum show general agreement with predictions from recent reactor models, the measured evolution in total IBD yield disagrees with recent predictions at 3.1$\sigma$. This discrepancy indicates that an overall deficit in measured flux with respect to predictions does not result from equal fractional deficits from the primary fission isotopes $^{235}$U, $^{239}$Pu, $^{238}$U, and $^{241}$Pu. Based on measured IBD yield variations, yields of $(6.17 \pm 0.17)$ and $(4.27 \pm 0.26) \times 10^{-43}$ cm$^2$/fission have been determined for the two dominant fission parent isotopes $^{235}$U and $^{239}$Pu. A 7.8% discrepancy between the observed and predicted $^{235}$U yield suggests that this isotope may be the primary contributor to the reactor antineutrino anomaly.

The precise measurement of neutrino properties is among the highest priorities in fundamental particle physics, involving many experiments worldwide. Since the experiments rely on the interactions of neutrinos with bound nucleons inside atomic nuclei, the planned advances in the scope and precision of these experiments requires a commensurate effort in the understanding and modeling of the hadronic and nuclear physics of these interactions, which is incorporated as a nuclear model in neutrino event generators. This model is essential to every phase of experimental analyses and its theoretical uncertainties play an important role in interpreting every result. In this White Paper we discuss in detail the impact of neutrino-nucleus interactions, especially the nuclear effects, on the measurement of neutrino properties using the determination of oscillation parameters as a central example. After an Executive Summary and a concise Overview of the issues, we explain how the neutrino event generators work, what can be learned from electron-nucleus interactions and how each underlying physics process - from quasi-elastic to deep inelastic scattering - is understood today. We then emphasize how our understanding must improve to meet the demands of future experiments. With every topic we find that the challenges can be met only with the active support and collaboration among specialists in strong interactions and electroweak physics that include theorists and experimentalists from both the nuclear and high energy physics communities.

A new measurement of the reactor antineutrino flux and energy spectrum by the Daya Bay reactor neutrino experiment is reported. The antineutrinos were generated by six 2.9~GW$_{\mathrm{th}}$ nuclear reactors and detected by eight antineutrino detectors deployed in two near (560~m and 600~m flux-weighted baselines) and one far (1640~m flux-weighted baseline) underground experimental halls. With 621 days of data, more than 1.2 million inverse beta decay (IBD) candidates were detected. The IBD yield in the eight detectors was measured, and the ratio of measured to predicted flux was found to be $0.946\pm0.020$ ($0.992\pm0.021$) for the Huber+Mueller (ILL+Vogel) model. A 2.9~$\sigma$ deviation was found in the measured IBD positron energy spectrum compared to the predictions. In particular, an excess of events in the region of 4-6~MeV was found in the measured spectrum, with a local significance of 4.4~$\sigma$. A reactor antineutrino spectrum weighted by the IBD cross section is extracted for model-independent predictions.

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

A measurement of electron antineutrino oscillation by the Daya Bay Reactor Neutrino Experiment is described in detail. Six 2.9-GW$_{\rm th}$ nuclear power reactors of the Daya Bay and Ling Ao nuclear power facilities served as intense sources of $\overline{\nu}_{e}$'s. Comparison of the $\overline{\nu}_{e}$ rate and energy spectrum measured by antineutrino detectors far from the nuclear reactors ($\sim$1500-1950 m) relative to detectors near the reactors ($\sim$350-600 m) allowed a precise measurement of $\overline{\nu}_{e}$ disappearance. More than 2.5 million $\overline{\nu}_{e}$ inverse beta decay interactions were observed, based on the combination of 217 days of operation of six antineutrino detectors (Dec. 2011--Jul. 2012) with a subsequent 1013 days using the complete configuration of eight detectors (Oct. 2012--Jul. 2015). The $\overline{\nu}_{e}$ rate observed at the far detectors relative to the near detectors showed a significant deficit, $R=0.949 \pm 0.002(\mathrm{stat.}) \pm 0.002(\mathrm{syst.})$. The energy dependence of $\overline{\nu}_{e}$ disappearance showed the distinct variation predicted by neutrino oscillation. Analysis using an approximation for the three-flavor oscillation probability yielded the flavor-mixing angle $\sin^22\theta_{13}=0.0841 \pm 0.0027(\mathrm{stat.}) \pm 0.0019(\mathrm{syst.})$ and the effective neutrino mass-squared difference of $\left|{\Delta}m^2_{\mathrm{ee}}\right|=(2.50 \pm 0.06(\mathrm{stat.}) \pm 0.06(\mathrm{syst.})) \times 10^{-3}\ {\rm eV}^2$. Analysis using the exact three-flavor probability found ${\Delta}m^2_{32}=(2.45 \pm 0.06(\mathrm{stat.}) \pm 0.06(\mathrm{syst.})) \times 10^{-3}\ {\rm eV}^2$ assuming the normal neutrino mass hierarchy and ${\Delta}m^2_{32}=(-2.56 \pm 0.06(\mathrm{stat.}) \pm 0.06(\mathrm{syst.})) \times 10^{-3}\ {\rm eV}^2$ for the inverted hierarchy.

This Letter reports an improved search for light sterile neutrino mixing in the electron antineutrino disappearance channel with the full configuration of the Daya Bay Reactor Neutrino Experiment. With an additional 404 days of data collected in eight antineutrino detectors, this search benefits from 3.6 times the statistics available to the previous publication, as well as from improvements in energy calibration and background reduction. A relative comparison of the rate and energy spectrum of reactor antineutrinos in the three experimental halls yields no evidence of sterile neutrino mixing in the $2\times10^{-4} \lesssim |\Delta m^{2}_{41}| \lesssim 0.3$ eV$^{2}$ mass range. The resulting limits on $\sin^{2}2\theta_{14}$ are improved by approximately a factor of 2 over previous results and constitute the most stringent constraints to date in the $|\Delta m^{2}_{41}| \lesssim 0.2$ eV$^{2}$ region.

The disappearance of reactor $\bar{\nu}_e$ observed by the Daya Bay experiment is examined in the framework of a model in which the neutrino is described by a wave packet with a relative intrinsic momentum dispersion $\sigma_\text{rel}$. Three pairs of nuclear reactors and eight antineutrino detectors, each with good energy resolution, distributed among three experimental halls, supply a high-statistics sample of $\bar{\nu}_e$ acquired at nine different baselines. This provides a unique platform to test the effects which arise from the wave packet treatment of neutrino oscillation. The modified survival probability formula was used to fit Daya Bay data, providing the first experimental limits: $2.38 \cdot 10^{-17} < \sigma_{\rm rel} < 0.23$. Treating the dimensions of the reactor cores and detectors as constraints, the limits are improved: $10^{-14} \lesssim \sigma_{\rm rel} < 0.23$, and an upper limit of $\sigma_{\rm rel} <0.20$ is obtained. All limits correspond to a 95\% C.L. Furthermore, the effect due to the wave packet nature of neutrino oscillation is found to be insignificant for reactor antineutrinos detected by the Daya Bay experiment thus ensuring an unbiased measurement of the oscillation parameters $\sin^22\theta_{13}$ and $\Delta m^2_{32}$ within the plane wave model.

This article reports an improved independent measurement of neutrino mixing angle $\theta_{13}$ at the Daya Bay Reactor Neutrino Experiment. Electron antineutrinos were identified by inverse $\beta$-decays with the emitted neutron captured by hydrogen, yielding a data-set with principally distinct uncertainties from that with neutrons captured by gadolinium. With the final two of eight antineutrino detectors installed, this study used 621 days of data including the previously reported 217-day data set with six detectors. The dominant statistical uncertainty was reduced by 49%. Intensive studies of the cosmogenic muon-induced $^9$Li and fast neutron backgrounds and the neutron-capture energy selection efficiency, resulted in a reduction of the systematic uncertainty by 26%. The deficit in the detected number of antineutrinos at the far detectors relative to the expected number based on the near detectors yielded $\sin^22\theta_{13} = 0.071 \pm 0.011$ in the three-neutrino-oscillation framework. The combination of this result with the gadolinium-capture result is also reported.

This is a brief summary of the 11th Applied Antineutrino Physics 2015 workshop held at the Virginia Tech Arlington Research Facility from December 7-8, 2015.

The Daya Bay experiment was the first to report simultaneous measurements of reactor antineutrinos at multiple baselines leading to the discovery of $\bar{\nu}_e$ oscillations over km-baselines. Subsequent data has provided the world's most precise measurement of $\rm{sin}^22\theta_{13}$ and the effective mass splitting $\Delta m_{ee}^2$. The experiment is located in Daya Bay, China where the cluster of six nuclear reactors is among the world's most prolific sources of electron antineutrinos. Multiple antineutrino detectors are deployed in three underground water pools at different distances from the reactor cores to search for deviations in the antineutrino rate and energy spectrum due to neutrino mixing. Instrumented with photomultiplier tubes (PMTs), the water pools serve as shielding against natural radioactivity from the surrounding rock and provide efficient muon tagging. Arrays of resistive plate chambers over the top of each pool provide additional muon detection. The antineutrino detectors were specifically designed for measurements of the antineutrino flux with minimal systematic uncertainty. Relative detector efficiencies between the near and far detectors are known to better than 0.2%. With the unblinding of the final two detectors' baselines and target masses, a complete description and comparison of the eight antineutrino detectors can now be presented. This paper describes the Daya Bay detector systems, consisting of eight antineutrino detectors in three instrumented water pools in three underground halls, and their operation through the first year of eight detector data-taking.

We report a new measurement of electron antineutrino disappearance using the fully-constructed Daya Bay Reactor Neutrino Experiment. The final two of eight antineutrino detectors were installed in the summer of 2012. Including the 404 days of data collected from October 2012 to November 2013 resulted in a total exposure of 6.9$\times$10$^5$ GW$_{\rm th}$-ton-days, a 3.6 times increase over our previous results. Improvements in energy calibration limited variations between detectors to 0.2%. Removal of six $^{241}$Am-$^{13}$C radioactive calibration sources reduced the background by a factor of two for the detectors in the experimental hall furthest from the reactors. Direct prediction of the antineutrino signal in the far detectors based on the measurements in the near detectors explicitly minimized the dependence of the measurement on models of reactor antineutrino emission. The uncertainties in our estimates of $\sin^{2}2\theta_{13}$ and $|\Delta m^2_{ee}|$ were halved as a result of these improvements. Analysis of the relative antineutrino rates and energy spectra between detectors gave $\sin^{2}2\theta_{13} = 0.084\pm0.005$ and $|\Delta m^{2}_{ee}|= (2.42\pm0.11) \times 10^{-3}$ eV$^2$ in the three-neutrino framework.

This Letter reports a measurement of the flux and energy spectrum of electron antineutrinos from six 2.9~GW$_{th}$ nuclear reactors with six detectors deployed in two near (effective baselines 512~m and 561~m) and one far (1,579~m) underground experimental halls in the Daya Bay experiment. Using 217 days of data, 296,721 and 41,589 inverse beta decay (IBD) candidates were detected in the near and far halls, respectively. The measured IBD yield is (1.55 $\pm$ 0.04) $\times$ 10$^{-18}$~cm$^2$/GW/day or (5.92 $\pm$ 0.14) $\times$ 10$^{-43}$~cm$^2$/fission. This flux measurement is consistent with previous short-baseline reactor antineutrino experiments and is $0.946\pm0.022$ ($0.991\pm0.023$) relative to the flux predicted with the Huber+Mueller (ILL+Vogel) fissile antineutrino model. The measured IBD positron energy spectrum deviates from both spectral predictions by more than 2$\sigma$ over the full energy range with a local significance of up to $\sim$4$\sigma$ between 4-6 MeV. A reactor antineutrino spectrum of IBD reactions is extracted from the measured positron energy spectrum for model-independent predictions.

The US neutrino community gathered at the Workshop on the Intermediate Neutrino Program (WINP) at Brookhaven National Laboratory February 4-6, 2015 to explore opportunities in neutrino physics over the next five to ten years. Scientists from particle, astroparticle and nuclear physics participated in the workshop. The workshop examined promising opportunities for neutrino physics in the intermediate term, including possible new small to mid-scale experiments, US contributions to large experiments, upgrades to existing experiments, R&D plans and theory. The workshop was organized into two sets of parallel working group sessions, divided by physics topics and technology. Physics working groups covered topics on Sterile Neutrinos, Neutrino Mixing, Neutrino Interactions, Neutrino Properties and Astrophysical Neutrinos. Technology sessions were organized into Theory, Short-Baseline Accelerator Neutrinos, Reactor Neutrinos, Detector R&D and Source, Cyclotron and Meson Decay at Rest sessions.This report summarizes discussion and conclusions from the workshop.

We describe a new detector, called NuLat, to study electron anti-neutrinos a few meters from a nuclear reactor, and search for anomalous neutrino oscillations. Such oscillations could be caused by sterile neutrinos, and might explain the "Reactor Antineutrino Anomaly". NuLat, is made possible by a natural synergy between the miniTimeCube and mini-LENS programs described in this paper. It features a "Raghavan Optical Lattice" (ROL) consisting of 3375 boron or $^6$Li loaded plastic scintillator cubical cells 6.3\,cm (2.500") on a side. Cell boundaries have a 0.127\,mm (0.005") air gap, resulting in total internal reflection guiding most of the light down the 3 cardinal directions. The ROL detector technology for NuLat gives excellent spatial and energy resolution and allows for in-depth event topology studies. These features allow us to discern inverse beta decay (IBD) signals and the putative oscillation pattern, even in the presence of other backgrounds. We discuss here test venues, efficiency, sensitivity and project status.

Current models of antineutrino production in nuclear reactors predict detection rates and spectra at odds with the existing body of direct reactor antineutrino measurements. High-resolution antineutrino detectors operated close to compact research reactor cores can produce new precision measurements useful in testing explanations for these observed discrepancies involving underlying nuclear or new physics. Absolute measurement of the 235U-produced antineutrino spectrum can provide additional constraints for evaluating the accuracy of current and future reactor models, while relative measurements of spectral distortion between differing baselines can be used to search for oscillations arising from the existence of eV-scale sterile neutrinos. Such a measurement can be performed in the United States at several highly-enriched uranium fueled research reactors using near-surface segmented liquid scintillator detectors. We describe here the conceptual design and physics potential of the PROSPECT experiment, a U.S.-based, multi-phase experiment with reactor-detector baselines of 7-20 meters capable of addressing these and other physics and detector development goals. Current R&D status and future plans for PROSPECT detector deployment and data-taking at the High Flux Isotope Reactor at Oak Ridge National Laboratory will be discussed.

A search for light sterile neutrino mixing was performed with the first 217 days of data from the Daya Bay Reactor Antineutrino Experiment. The experiment's unique configuration of multiple baselines from six 2.9~GW$_{\rm th}$ nuclear reactors to six antineutrino detectors deployed in two near (effective baselines 512~m and 561~m) and one far (1579~m) underground experimental halls makes it possible to test for oscillations to a fourth (sterile) neutrino in the $10^{\rm -3}~{\rm eV}^{2} < |\Delta m_{41}^{2}| < 0.3~{\rm eV}^{2}$ range. The relative spectral distortion due to electron antineutrino disappearance was found to be consistent with that of the three-flavor oscillation model. The derived limits on $\sin^22\theta_{14}$ cover the $10^{-3}~{\rm eV}^{2} \lesssim |\Delta m^{2}_{41}| \lesssim 0.1~{\rm eV}^{2}$ region, which was largely unexplored.

A new measurement of the $\theta_{13}$ mixing angle has been obtained at the Daya Bay Reactor Neutrino Experiment via the detection of inverse beta decays tagged by neutron capture on hydrogen. The antineutrino events for hydrogen capture are distinct from those for gadolinium capture with largely different systematic uncertainties, allowing a determination independent of the gadolinium-capture result and an improvement on the precision of $\theta_{13}$ measurement. With a 217-day antineutrino data set obtained with six antineutrino detectors and from six 2.9 GW$_{th}$ reactors, the rate deficit observed at the far hall is interpreted as $\sin^22\theta_{13}=0.083\pm0.018$ in the three-flavor oscillation model. When combined with the gadolinium-capture result from Daya Bay, we obtain $\sin^22\theta_{13}=0.089\pm0.008$ as the final result for the six-antineutrino-detector configuration of the Daya Bay experiment.

A facility that can deliver beams of electron and muon neutrinos from the decay of a stored muon beam has the potential to unambiguously resolve the issue of the evidence for light sterile neutrinos that arises in short-baseline neutrino oscillation experiments and from estimates of the effective number of neutrino flavors from fits to cosmological data. In this paper, we show that the nuSTORM facility, with stored muons of 3.8 GeV/c $\pm$ 10%, will be able to carry out a conclusive muon neutrino appearance search for sterile neutrinos and test the LSND and MiniBooNE experimental signals with 10$\sigma$ sensitivity, even assuming conservative estimates for the systematic uncertainties. This experiment would add greatly to our knowledge of the contribution of light sterile neutrinos to the number of effective neutrino flavors from the abundance of primordial helium production and from constraints on neutrino energy density from the cosmic microwave background. The appearance search is complemented by a simultaneous muon neutrino disappearance analysis that will facilitate tests of various sterile neutrino models.

A measurement of the energy dependence of antineutrino disappearance at the Daya Bay Reactor Neutrino Experiment is reported. Electron antineutrinos ($\overline{\nu}_{e}$) from six $2.9$ GW$_{\rm th}$ reactors were detected with six detectors deployed in two near (effective baselines 512 m and 561 m) and one far (1579 m) underground experimental halls. Using 217 days of data, 41589 (203809 and 92912) antineutrino candidates were detected in the far hall (near halls). An improved measurement of the oscillation amplitude $\sin^{2}2\theta_{13} = 0.090^{+0.008}_{-0.009} $ and the first direct measurement of the $\overline{\nu}_{e}$ mass-squared difference $|\Delta m^{2}_{ee}|= (2.59_{-0.20}^{+0.19}) \times 10^{-3}\ {\rm eV}^2 $ is obtained using the observed $\overline{\nu}_{e}$ rates and energy spectra in a three-neutrino framework. This value of $|\Delta m^{2}_{ee}|$ is consistent with $|\Delta m^{2}_{\mu\mu}|$ measured by muon neutrino disappearance, supporting the three-flavor oscillation model.

A staged approach towards muon based facilities for Intensity and Energy Frontier science, building upon existing and proposed facilities at Fermilab, is presented. At each stage, a facility exploring new physics also provides an R&D platform to validate the technology needed for subsequent stages. The envisioned program begins with nuSTORM, a sensitive sterile neutrino search which also provides precision neutrino cross-section measurements while developing the technology of using and cooling muons. A staged Neutrino Factory based upon Project X, sending beams towards the Sanford Underground Research Facility (SURF), which will house the LBNE detector, could follow for detailed exploration of neutrino properties at the Intensity Frontier, while also establishing the technology of using intense bunched muon beams. The complex could then evolve towards Muon Colliders, starting at 126 GeV with measurements of the Higgs resonance to sub-MeV precision, and continuing to multi-TeV colliders for the exploration of physics beyond the Standard Model at the Energy Frontier. An Appendix addresses specific questions raised by the Lepton Colliders subgroup of the CSS2013 Frontier Capabilities Study Group.

We study the possible impact of nuclear effects and final state interactions on the determination of the oscillation parameters due to mis-reconstruction of non-quasi-elastic events as quasi-elastic events at low energies. We analyze a $\nu_\mu$ disappearance experiment using a water \v{C}erenkov detector. We find that, if completely ignored in the fit, nuclear effects can induce a significant bias in the determination of atmospheric oscillation parameters, particularly for the atmospheric mixing angle. Even after inclusion of a near detector a bias in the determination of the atmospheric mixing angle comparable to the statistical error remains.