Providing a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using data collected from December 2011 to August 2017, a search has been performed for electron-antineutrino signals coinciding with detected GW events, including GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817. We used three time windows of $\mathrm{\pm 10~s}$, $\mathrm{\pm 500~s}$, and $\mathrm{\pm 1000~s}$ relative to the occurrence of the GW events, and a neutrino energy range of 1.8 to 100 MeV to search for correlated neutrino candidates. The detected electron-antineutrino candidates are consistent with the expected background rates for all the three time windows. Assuming monochromatic spectra, we found upper limits (90% confidence level) on electron-antineutrino fluence of $(1.13~-~2.44) \times 10^{11}~\rm{cm^{-2}}$ at 5 MeV to $8.0 \times 10^{7}~\rm{cm^{-2}}$ at 100 MeV for the three time windows. Under the assumption of a Fermi-Dirac spectrum, the upper limits were found to be $(5.4~-~7.0)\times 10^{9}~\rm{cm^{-2}}$ for the three time windows.

The 4.8$\sigma$ anomaly in MiniBooNE data cannot be reconciled with MINOS+ and IceCube data within the vanilla framework of neutrino oscillations involving an eV-mass sterile neutrino. We show that an apparently consistent picture can be drawn if charged-current and neutral-current nonstandard neutrino interactions are at work in the 3+1 neutrino scheme. It appears that either the neutrino sector is more elaborate than usually envisioned, or one or more datasets needs revision.

Couplings between light scalar dark matter (DM) and neutrinos induce a perturbation to the neutrino mass matrix. If the DM oscillation period is smaller than ten minutes (or equivalently, the DM particle is heavier than $0.69\times10^{-17}$ eV), the fast-averaging over an oscillation cycle leads to a modification of the measured oscillation parameters. We present a specific $\mu-\tau$ symmetric model in which the measured value of $\theta_{13}$ is entirely generated by the DM interaction, and which reproduces the other measured oscillation parameters. For a scalar DM particle lighter than $10^{-15}$ eV, adiabatic solar neutrino propagation is maintained. A suppression of the sensitivity to CP violation at long-baseline neutrino experiments is predicted in this model. We find that DUNE cannot exclude the DM scenario at more than $3\sigma$ C.L. for bimaximal, tribimaximal and hexagonal mixing, while JUNO can rule it out at more than $6\sigma$ C.L. by precisely measuring both $\theta_{12}$ and $\theta_{13}$.

Measurements of the solar neutrino mass-squared difference from KamLAND and solar neutrino data are somewhat discrepant, perhaps due to nonstandard neutrino interactions in matter. We show that the zenith angle distribution of solar neutrinos at Hyper-Kamiokande and the energy spectrum of reactor antineutrinos at JUNO can conclusively confirm the discrepancy and detect new neutrino interactions.

The NSF has chosen the site for the Deep Underground Science and Engineering Laboratory (DUSEL) to be in Lead, South Dakota. In fact, the state of South Dakota has already stepped up to the plate and contributed its own funding for the proposed lab, see http://www.sanfordlaboratoryathomestake.org/index.html. The final decision by NSF for funding the Initial Suite of Experiments for DUSEL will be made early in 2009. At that time the NSF Science Board must make a decision. Of order 200 experimentalists have already expressed an interest in performing experiments at DUSEL. In order to assess the interest of the theoretical community, the Center for Cosmology and Astro-Particle Physics (CCAPP) at The Ohio State University (OSU) organized a 3-day DUSEL Theory Workshop in Columbus, Ohio from April 4 - 6, 2008. The workshop focused on the scientific case for six proposed experiments for DUSEL: long baseline neutrino oscillations, proton decay, dark matter, astrophysical neutrinos, neutrinoless double beta decay and N-Nbar oscillations. The outcome of this workshop is the DUSEL Theory White paper addressing the scientific case at a level which may be useful in the decision making process for policy makers at the NSF and in the U.S. Congress. In order to assess the physics interest in the DUSEL project we have posted the DUSEL Theory White paper on the following CCAPP link http://ccapp.osu.edu/whitepaper.html . Please read the white paper and, if you are interested, use the link to show your support by co-signing the white paper.

We study the proposal that mass-varying neutrinos could provide an explanation for the LSND signal for \bar\nu_mu to \bar\nu_e oscillations. We first point out that all positive oscillation signals occur in matter and that three active mass-varying neutrinos are insufficient to describe all existing neutrino data including LSND. We then examine the possibility that a model with four mass-varying neutrinos (three active and one sterile) can explain the LSND effect and remain consistent with all other neutrino data. We find that such models with a 3+1 mass structure in the neutrino sector may explain the LSND data and a null MiniBooNE result for 0.10 < \sin^2 2\theta_x < 0.30. Predictions of the model include a null result at Double-CHOOZ, but positive signals for underground reactor experiments and for \nu_\mu to \nu_e oscillations in long-baseline experiments.

We perform a comparative study of two methods of determining the survival probabilities of low, intermediate, and high energy solar neutrinos that emphasizes the general agreement between the Large Mixing Angle (LMA) solution and extant solar neutrino data. The first analysis is oscillation parameter-independent and the second analysis involves an approximate calculation of the survival probabilities in the three energy ranges that depends only on oscillation parameters. We show that future experiments like BOREXino, CLEAN, Heron, LENS and MOON, that measure $pp$ and $^7$Be neutrinos, will facilitate a stringent test of the LMA solution independently of the Standard Solar Model (SSM), without recourse to earth-matter effects. Throughout, we describe the role of SSM assumptions on our results. If the LMA solution passes the test without needing to be modified, it may be possible to establish that $\theta_x$ is nonzero at more than $2\sigma$ assuming the SSM prediction for the $pp$ flux is correct.

The current status of the physics of massive neutrinos is reviewed with a forward-looking emphasis. The article begins with the general phenomenology of neutrino oscillations in vacuum and matter and documents the experimental evidence for oscillations of solar, reactor, atmospheric and accelerator neutrinos. Both active and sterile oscillation possibilities are considered. The impact of cosmology (BBN, CMB, leptogenesis) and astrophysics (supernovae, highest energy cosmic rays) on neutrino observables and vice versa, is evaluated. The predictions of grand unified, radiative and other models of neutrino mass are discussed. Ways of determining the unknown parameters of three-neutrino oscillations are assessed, taking into account eight-fold degeneracies in parameters that yield the same oscillation probabilities, as well as ways to determine the absolute neutrino mass scale (from beta-decay, neutrinoless double-beta decay, large scale structure and Z-bursts). Critical unknowns at present are the amplitude of \nu_\mu to \nu_e oscillations and the hierarchy of the neutrino mass spectrum; the detection of CP violation in the neutrino sector depends on these and on an unknown phase. The estimated neutrino parameter sensitivities at future facilities (reactors, superbeams, neutrino factories) are given. The overall agenda of a future neutrino physics program to construct a bottom-up understanding of the lepton sector is presented.

We perform a global analysis in the framework of two active neutrino oscillations of all solar neutrino data, including the recent SNO day and night spectra (comprised of the charged current (CC), elastic scattering (ES) and neutral current (NC) events), the Super-Kamiokande (SK) day and night spectra (from 1496 days) and the updated SAGE results. We find that the Large Mixing Angle (LMA) solution is selected at the 99% C.L.; the best-fit parameters are \Delta m^2=5.6 \times 10^{-5} eV^2 and \theta=32^{\circ}. No solutions with \theta\geq \pi/4 are allowed at the 5\sigma C.L. Oscillations to a pure sterile state are excluded at 5.3\sigma, but a sizeable sterile neutrino component could still be present in the solar flux.

We examine how constraints can be placed on the neutrino component of dark matter by an accurate measurement of neutrinoless double beta ($0\nu\beta\beta$) decay and the solar oscillation amplitude. We comment on the alleged evidence for $0\nu\beta\beta$ decay.

We analyze the 1258-day Super-Kamiokande day and night solar neutrino energy spectra with various $\chi^2$ definitions. The best-fit lies in the LMA region at $(\Delta m^2, tan^2 \theta)=(5.01\times 10^{-5} eV^2, 0.60)$, independently of whether systematic errors are included in the $\chi^2$-definition. We compare the exclusion and allowed regions from the different definitions and choose the most suitable definition to predict the regions from SNO at the end of three years of data accumulation. We first work under the assumption that Super-Kamiokande sees a flux-suppressed flat energy spectrum. Then, we consider the possibility of each one of the three MSW regions being the solution to the solar neutrino problem. We find that the exclusion and allowed regions for the flat spectrum hypothesis and the LMA and LOW solutions are alike. In three years, we expect SNO to find very similar regions to that obtained by Super-Kamiokande. We evaluate whether the zenith angle distribution at SNO with optimum binning will add anything to the analysis of the day and night spectra; for comparison, we show the results of our analysis of the 1258-day zenith angle distribution from Super-Kamiokande, for which the best-fit parameters are $(\Delta m^2, tan^2 \theta)=(5.01\times 10^{-5} eV^2, 0.56)$.

A revolution in our understanding of the neutrino sector is underway, driven by observations that are interpreted in terms of changes in neutrino flavors as they propagate. Since neutrino oscillations occur only if neutrinos are massive, these phenomena indicate physics beyond the Standard Model. With the present evidence for oscillations from atmospheric, solar and accelerator data, we are already able to begin to make strong inferences about the mass spectrum and mixings of neutrinos. Theoretical efforts to achieve a synthesis have produced a variety of models with differing testable consequences. A combination of particle physics, nuclear physics and astrophysics is needed for a full determination of the fundamental properties of neutrinos. This article reviews what has been achieved thus far and the future prospects for understanding the nature of neutrino masses and mixing.

When three Majorana neutrinos describe the solar and atmospheric neutrino data via oscillations, a nonzero measurement of neutrinoless double beta ($0\nu\beta\beta$) decay can determine the sum of neutrino masses $\sum m_\nu$ if the solar solution has small-angle mixing, and place a lower bound on $\sum m_\nu$ for large-angle solar mixing. If in addition a nonzero $\sum m_\nu$ is deduced from cosmology, the neutrino mass spectrum may be uniquely specified for some ranges of neutrino parameters. For $\sum m_\nu > 0.75$ eV, the small-angle solar solution is excluded by the current upper limit on neutrinoless double beta decay. In models with maximal solar mixing the $CP$ phases of the neutrinos may be strongly constrained by stringent upper bounds on $0\nu\beta\beta$ decay.

We determine the three-neutrino mixing and mass parameters that are allowed by the solar and atmospheric neutrino data when vacuum oscillations are responsible for both phenomena. The global fit does not appreciably change the allowed regions for the parameters obtained from effective two-neutrino fits. We discuss how measurements of the solar electron energy spectrum below 6.5 GeV in Super-Kamiokande and seasonal variations in the Super-Kamiokande, $^{71}$Ga, and BOREXINO experiments can distinguish the different solar vacuum solutions.

We examine the prospects for determining the neutrino mixing matrix and for observing CP and T violation in neutrino oscillations in four-neutrino models. We focus on a general class of four-neutrino models with two pairs of nearly degenerate mass eigenstates separated by approximately 1 eV, which can describe the solar, atmospheric and LSND neutrino data. We present a general parametrization of these models and discuss in detail the determination of the mixing parameters and the mass matrix texture from current and future neutrino data in the case where $\nu_e$ and $\nu_\mu$ each mix primarily with one other neutrino. We find that measurable CP/T violation in long-baseline experiments, with amplitude at the level of the LSND signal, is possible given current experimental constraints. Also, additional oscillation effects in short- and long-baseline experiments may be measurable in many cases. We point out that, given separate scales for the mass-squared differences of the solar and atmospheric oscillations, observable CP/T violation effects in neutrino oscillations signals the existence of a sterile neutrino. We examine several textures of the neutrino mass matrix and determine which textures can have measurable CP/T violation in neutrino oscillations in long-baseline experiments. We also briefly discuss some possible origins of the neutrino mass terms in straightforward extensions of the Standard Model.

By combining experimental constraints from atmospheric and solar neutrino oscillations and the tritium beta decay endpoint, we infer upper limits of 4.4 eV on the $\nu_\mu$ and $\nu_\tau$ masses, if the universe consists of three neutrinos. For hierarchical mass spectra $m_3 \gg m_1, m_2$ or $m_3 \simeq m_2 \gg m_1$ we infer that $m_{\nu_\alpha} \lsim 0.08$ eV for $\alpha = e, \mu, \tau$. In a four neutrino universe, and assuming neutrino oscillations also account for the LSND experimental results, $m_{\nu_\alpha} \lsim 5.4$ eV. We also obtain lower limits on the masses.

We show that if the solar and atmospheric data are both described by maximal vacuum oscillations at the relevant mass scales then there exists a unique mixing matrix for three neutrino flavors. The solution necessarily conserves CP and automatically implies that there is no disappearance of atmospheric $\nu_e$, consistent with indications from the Super-Kamiokande experiment. We also investigate the consequences for three-neutrino mixing if the solar and atmospheric oscillations exhibit mixing that is large but not maximal. For non-maximal mixing $\nu_e\leftrightarrow\nu_\tau$ and $\nu_e\leftrightarrow\nu_\mu$ oscillations are predicted that may be observable in future long-baseline experiments.

We make a model-independent analysis of all available data that indicate neutrino oscillations. Using probability diagrams, we confirm that a mass spectrum with two nearly degenerate pairs of neutrinos separated by a mass gap of $\simeq1$ eV is preferred over a spectrum with one mass eigenstate separated from the others. We derive some new relations among the four-neutrino mixing matrix elements. We design four-neutrino mass matrices with three active neutrinos and one sterile neutrino that naturally incorporate maximal oscillations of atmospheric $\nu_\mu$ and explain the solar neutrino and LSND results. The models allow either a large or small angle MSW or vacuum oscillation description of the solar neutrino deficit. The models predict (i) oscillations of either $\nu_e \to \nu_\tau$ or $\nu_e \to\nu_s$ in long-baseline experiments at $L/E \gg 1$ km/GeV, with amplitude determined by the LSND oscillation amplitude and argument given by the atmospheric $\delta m^2$, and (ii) the equality of the $\nu_e$ disappearance probability, the $\nu_\mu$ disappearance probability, and the LSND $\nu_\mu\to\nu_e$ appearance probability in short-baseline experiments.

We present a four-neutrino model with three active neutrinos and one sterile neutrino which naturally has maximal $\nu_\mu\to\nu_\tau$ oscillations of atmospheric neutrinos and can explain the solar neutrino and LSND results. The model predicts $\nu_e \to \nu_\tau$ and $\nu_e \to \nu_\mu$ oscillations in long-baseline experiments with $L/E \gg 1$ km/GeV with amplitudes that are determined by the LSND oscillation amplitude and argument controlled by the atmospheric $\delta m^2$.