We present a comprehensive review of keV-scale sterile neutrino Dark Matter, collecting views and insights from all disciplines involved - cosmology, astrophysics, nuclear, and particle physics - in each case viewed from both theoretical and experimental/observational perspectives. After reviewing the role of active neutrinos in particle physics, astrophysics, and cosmology, we focus on sterile neutrinos in the context of the Dark Matter puzzle. Here, we first review the physics motivation for sterile neutrino Dark Matter, based on challenges and tensions in purely cold Dark Matter scenarios. We then round out the discussion by critically summarizing all known constraints on sterile neutrino Dark Matter arising from astrophysical observations, laboratory experiments, and theoretical considerations. In this context, we provide a balanced discourse on the possibly positive signal from X-ray observations. Another focus of the paper concerns the construction of particle physics models, aiming to explain how sterile neutrinos of keV-scale masses could arise in concrete settings beyond the Standard Model of elementary particle physics. The paper ends with an extensive review of current and future astrophysical and laboratory searches, highlighting new ideas and their experimental challenges, as well as future perspectives for the discovery of sterile neutrinos.

This white paper addresses the hypothesis of light sterile neutrinos based on recent anomalies observed in neutrino experiments and the latest astrophysical data.

Fractionally charged massive particles (FCHAMPs) appear in extensions of the standard model, especially those with superstring constructions. The lightest FCHAMP would be absolutely stable and any produced during the early evolution of the Universe would be present today. The production, annihilation, and survival of L, a lepton with electroweak but no strong interactions, of mass m_L and charge Q_L (in units of the positron charge) are explored. Since massive charged particles behave like baryons, primordial nucleosynthesis and the cosmic background radiation temperature anisotropies limit the FCHAMP relic density leading to constraints on the Q_L - m_L relation. Further constraints are provided by the invisible width of the Z and by accelerator searches for massive, charged particles. We exploit the fact that in the early Universe the negatively charged L will combine with alpha particles and protons forming tightly bound, positively charged states. The Coulomb barriers between these positively charged bound states and the free L+ suppress late time annihilation in the Galaxy and on Earth, limiting any late-time reduction of relic FCHAMP pairs. The surviving FCHAMP abundance on Earth is orders of magnitude higher than the limits from terrestrial searches for fractionally charged particles, appearing to close the window on FCHAMPs. However, as Q_L approaches an integer these searches become increasingly insensitive, leaving some "islands" in the Q_L - m_L plane which may be explored by searching for FCHAMPs in the cosmic rays.

We investigate the correlated predictions of singlet extended MSSM models for direct detection of the lightest neutralino with its cosmological relic density. To illustrate the general effects of the singlet, we take heavy sleptons and squarks. We apply LEP, $(g-2)_\mu$ and perturbativity constraints. We find that the WMAP upper bound on the cold dark matter density limits much of the parameter space to regions where the lightest neutralino can be discovered in recoil experiments. The results for the NMSSM and UMSSM are typically similar to the MSSM since their light neutralinos have similar compositions and masses. In the n/sMSSM the neutralino is often very light and its recoil detection is within the reach of the CDMS II experiment. In general, most points in the parameter spaces of the singlet models we consider are accessible to the WARP experiment.

We construct an anomaly free supersymmetric U(1)' model with a secluded U(1)'-breaking sector. We study the one-loop effective potential at finite temperature, and show that there exists a strong enough first order electroweak phase transition for electroweak baryogenesis (EWBG) because of the large trilinear term $A_h h S H_d H_u$ in the tree-level Higgs potential. Unlike in the MSSM, the lightest stop can be very heavy. We consider the non-local EWBG mechanism in the thin wall regime, and find that within uncertainties the observed baryon number can be generated from the $\tau$ lepton contribution, with the secluded sector playing an essential role. The chargino and neutralino contributions and the implications for the Z' mass and electric dipole moments are briefly discussed.

We quantify the extent to which extra relativistic energy density can be concealed by a neutrino asymmetry without conflicting with the baryon asymmetry measured by the Wilkinson Microwave Anisotropy Probe (WMAP). In the presence of a large electron neutrino asymmetry, slightly more than seven effective neutrinos are allowed by Big Bang Nucleosynthesis (BBN) and WMAP at 2\sigma. The same electron neutrino degeneracy that reconciles the BBN prediction for the primordial helium abundance with the observationally inferred value also reconciles the LSND neutrino with BBN by suppressing its thermalization prior to BBN.

Time variation of fundamental constants would not be surprising in the framework of theories involving extra dimensions. The variation of any one constant is likely to be correlated with variations of others in a pattern that is diagnostic of the underlying physics.

Alternatives to the traditional grand unified theory seesaw for neutrino masses are briefly described. These include the possibility of large extra dimensions and various possibilities for models involving an extra U(1)' gauge symmetry. The difficulty of observing Majorana phases in neutrinoless double beta decay is also briefly commented on.

In models involving new TeV-scale Z' gauge bosons, the new U(1)' symmetry often prevents the generation of Majorana masses needed for a conventional neutrino seesaw, leading to three superweakly interacting ``right-handed'' neutrinos nu_R, the Dirac partners of the ordinary neutrinos. These can be produced prior to big bang nucleosynthesis by the Z' interactions, leading to a faster expansion rate and too much ^4He. We quantify the constraints on the Z' properties from nucleosynthesis for Z' couplings motivated by a class of E_6 models parametrized by an angle theta_E6. The rate for the annihilation of three approximately massless right-handed neutrinos into other particle pairs through the Z' channel is calculated. The decoupling temperature, which is higher than that of ordinary left-handed neutrinos due to the large Z' mass, is evaluated, and the equivalent number of new doublet neutrinos Delta N_nu is obtained numerically as a function of the Z' mass and couplings for a variety of assumptions concerning the Z-Z' mixing angle and the quark-hadron transition temperature T_c. Except near the values of theta_E6 for which the Z' decouples from the right-handed neutrinos, the Z' mass and mixing constraints from nucleosynthesis are much more stringent than the existing laboratory limits from searches for direct production or from precision electroweak data, and are comparable to the ranges that may ultimately be probed at proposed colliders. For the case T_c = 150 MeV with the theoretically favored range of Z-Z' mixings, Delta N_nu < 0.3 for M_Z' > 4.3 TeV for any value of theta_E6. Larger mixing or larger T_c often lead to unacceptably large Delta N_nu except near the nu_R decoupling limit.

The existence of bulk sterile neutrinos in theories with large extra dimensions can naturally explain small 4-dimensional Dirac masses for the active neutrinos. We study a model with 3 bulk neutrinos and derive various constraints on the size of the extra dimensions from neutrino oscillation experiments. Our analysis includes recent solar and atmospheric data from SNO and Super-Kamiokande, respectively, as well as various reactor and accelerator experiments. In addition, we comment on possible extensions of the model that could accommodate the LSND results, using natural parameters only. Our most conservative bound, obtained from the atmospheric data in the hierarchical mass scheme, constrains the largest extra dimension to have a radius R < 0.82 microns. Thus, in the context of the model studied here, future gravitational experiments are unlikely to observe the effects of extra dimensions.

Unification of the gauge couplings would imply that time variations of the fine structure constant are accompanied by calculable and very significant time variations in the QCD scale parameter $\Lambda_{QCD}$. Since $\Lambda_{QCD}$ is the dominant factor in setting the hadron masses, estimates made by simple variations of the fine structure constant may not provide meaningful limits. There may also be related variations in Yukawa couplings and the electroweak scale. Implications for the 21 cm hyperfine transition, big bang nucleosynthesis, and the triple alpha process are discussed. We find that the first of these already provides strong constraints on the underlying theory. It is emphasized more generally that time (and space) variations of fundamental couplings and their correlations may be a significant probe of ultra-high-energy physics.

The NOW2000 Workshop summarized the present status and future possibilities for all aspects of neutrino physics and astrophysics. Neutrino oscillation physics has truly come of age in the last few years. It is now data driven (analogous to cosmology and much of high energy physics in general). Experimental techniques and the theoretical interpretation have developed dramatically. An example of the latter includes the more realistic analysis of neutrino oscillations in the framework of 3 or more mass eigenstates. The phenomenological emphasis has shifted from degenerate to hierarchical neutrino spectra.

Efficient oscillations between ordinary (active) and sterile neutrinos can occur only if Dirac and Majorana mass terms exist which are both small and comparable. It is shown that this can occur naturally in a class of string models, in which higher-dimensional operators in the superpotential lead to an intermediate scale expectation value for a scalar field and to suppressed Dirac and Majorana fermion masses.

We present an updated analysis of astrophysical solutions, two-flavor MSW solutions, and vacuum oscillation solutions to the solar neutrino anomaly. The recent results of each of the five solar neutrino experiments are incorporated, including both the zenith angle (day-night) and spectral information from the Kamiokande experiment, and the preliminary Super-Kamiokande results. New theoretical developments include the use of the most recent Bahcall-Pinsonneault flux predictions (and uncertainties) and density and production profiles, the radiative corrections to the neutrino-electron scattering cross section, and new constraints on the Ga absorption cross section inferred from the gallium source experiments. From a model independent analysis, arbitrary astrophysical solutions are excluded at more that 98% C.L. even if one ignores any one of the three classes of experiment, relaxes the luminosity constraint, or allows more suppression of the 7Be than 8B flux. The data is well described by large and small mixing angle two-flavor MSW conversions, MSW conversions into a sterile neutrino with small mixing, or vacuum oscillations. We also present MSW fits for nonstandard solar models parameterized by an arbitrary solar core temperature or arbitrary 8B flux.

Constraints on big bang nucleosynthesis (BBN) and on cosmological parameters from conflicting deuterium observations in different high red-shift QSO systems are discussed. The high deuterium observations by Carswell et al., Songaila et al., and Rugers & Hogan is consistent with 4He and 7Li observations and Standard BBN (N_nu =3) and allows N_nu \leq 3.6 at 95% C.L., but is inconsistent with local observations of D and 3He in the context of conventional theories of stellar and Galactic evolution. In contrast, the low deuterium observations by Tytler, Fan & Burles and Burles & Tytler are consistent with the constraints from local Galactic observations, but require N_nu = 1.9 +/- 0.3 at 68% C.L., excluding Standard BBN at 99.9% C.L., unless the systematic uncertainties in the 4He observations have been underestimated by a large amount. The high and low primordial deuterium abundances imply, respectively, Omega_B h^2 = 0.005 - 0.01 and Omega_B h^2 = 0.02 - 0.03 at 95% C.L. When combined with the high baryon fraction inferred from x-ray observations of rich clusters, the corresponding total mass densities (for 50 \le H_0 \le 90) are Omega_M = 0.05 - 0.20 and Omega_M = 0.2 - 0.7, respectively (95\% C.L.) The range of Omega_M corresponding to high D is in conflict with dynamical constraints (Omega_M > 0.2 - 0.3) and with the shape parameter constraint (Gamma = Omega_M h = 0.25 +/- 0.05) from large scale structure formation in CDM and Lambda-CDM models.

A new evaluation of the constraint on the number of light neutrino species (N_nu) from big bang nucleosynthesis suggests a discrepancy between the predicted light element abundances and those inferred from observations, unless the inferred primordial 4He abundance has been underestimated by 0.014 +/- 0.004 (1 sigma) or less than 10% (95%C.L.) of 3He survives stellar processing. With the quoted systematic errors in the observed abundances and a conservative chemical evolution parameterization, the best fit to the combined data is N_nu = 2.1 +/- 0.3 (1 sigma) and the upper limit is N_nu < 2.6 (95% C.L.). The data are inconsistent with the Standard Model (N_nu = 3) at the 98.6% C.L.

The status of the solar neutrino problem is reviewed. Attempts to explain the observed deficit and spectral distortion, both by astrophysical and particle physics methods, are described. It is argued that the comparison of all experiments strongly prefers the particle physics solutions.

Besides the opportunity for discovering new neutrino physics, solar neutrino measurements provide a sensitive probe of the solar interior, and thus a rigorous test of solar model predictions. We present model independent determinations of the neutrino spectrum by using relevant flux components as free parameters subject only to the luminosity constraint. (1) Without the Mikheyev-Smirnov-Wolfenstein (MSW) effect, the best fit for the combined data is poor. Furthermore, the data indicate a severe suppression of the $^7$Be flux relative to the $^8$B, contradicting both standard and nonstandard solar models in general; the $pp$ flux takes its maximum value allowed by the luminosity constraint. This pathology consistently appears even if we ignore any one of the three data. (2) In the presence of the two-flavor MSW effect, the current constraint on the initial $^8$B flux is weak, but consistent with the SSM and sufficient to exclude nonstandard models with small $^8$B fluxes. No meaningful constraint is obtained for the other fluxes. In the future, even allowing MSW, the $^8$B and $^7$Be fluxes can be determined at the $\pm$(15 -- 20)\% level, making competing solar models distinguishable. We emphasize that the neutral current sensitivity for $^7$Be neutrinos in BOREXINO, HELLAZ, and HERON is essential for determining the initial fluxes. The constraints on the MSW parameters in the model independent analysis are also discussed.

Various theoretical uncertainties in the standard solar model and in the Mikheyev-Smirnov-Wolfenstein (MSW) analysis are discussed. It is shown that two methods of estimating the solar neutrino flux uncertainties are equivalent: (a) a simple parametrization of the uncertainties using the core temperature and the nuclear production cross sections; (b) the Monte Carlo method of Bahcall and Ulrich. In the MSW analysis, we emphasize proper treatments of correlation of theoretical uncertainties between flux components and between different detectors, the Earth effect, and multiple solutions in a combined $\chi^2$ procedure. The MSW solutions for various standard and nonstandard solar models are also shown. The MSW predictions of the global solutions for the future solar neutrino experiments are given, emphasizing the measurement of the energy spectrum and the day-night effect in Sudbury Neutrino Observatory and Super-Kamiokande to distinguish the two solutions.

We update the analysis of the MSW and general astrophysical solutions to the combined solar neutrino observations by including the GALLEX II result. We also show that our parametrized flux uncertainties are equivalent to the Monte-Carlo results of Bahcall and Ulrich.

We consider the most general solar model, using the neutrino fluxes as free parameters constrained only by the solar luminosity, and show that the combined solar neutrino data exclude any astrophysical solution at 98\% C.L.\ Our best fit to the $^7$Be and $^8$B fluxes is respectively $<$7\% and 37$\pm$4\% of the standard solar model prediction, but only with a large $\chi^2$ (5.6 for 1 d.f.). This best fit to the fluxes contradicts explicit nonstandard solar models, which generally reduce the $^8$B flux more than the $^7$Be. Those models are well parameterized by a single parameter, the central temperature.

We consider the Earth effect in the MSW analysis of the Homestake, Kamiokande, GALLEX, and SAGE solar neutrino experiments. Using the time-averaged data and assuming two-flavor oscillations, the large-angle region of the combined fit extends to much smaller angles (to $\sin^22\theta \simeq 0.1$) than when the Earth effect is ignored. However, the additional constraint from the Kamiokande II day-night data excludes most of the parameter space sensitive to the Earth effect independent of astrophysical uncertainties, and leaves only a small large-angle region close to maximal mixing at 90\% C.L. The nonadiabatic solution remains unaffected by the Earth effect and is still preferred. Both theoretical and experimental uncertainties are included in the analysis.

Constraints on the core temperature (T_c) of the Sun and on neutrino- oscillation parameters are obtained from the existing solar neutrino data, including the recent GALLEX and Kamiokande III results. (1) A purely astrophysical solution to the solar neutrino problem is strongly disfavored by the data: the best fit in a cooler Sun model requires an 8% reduction in T_c, but the chi-sqaured test rejects this hypothesis at 99.99% C.L., suggesting new neutrino physics. (2) Assuming the Standard Solar Model (SSM) and MSW oscillations, the MSW parameters are constrained to two small regions: one in the non-adiabatic region and the other in the large-mixing region. The non-adiabatic solution gives a considerably better fit. For nu_e oscillations into sterile neutrinos, the allowed region (90%) is constrained to non- adiabatic oscillations. As long as the SSM is assumed, the neutrino mixing angles are at least four times larger, or considerably smaller, than the corresponding quark mixing angles. (3) Allowing both MSW oscillations and a non-standard core temperature, a) the experiments determine the core temperature at the 5% level, yielding a value consistent with the SSM prediction. b) When T_c is used as a free parameter, the allowed MSW region is broadened: a cooler Sun (T_c=0.95) allows mass and mixing implied by the SUSY SO(10) GUT, while a warmer Sun (T_c=1.05) allows parameter space suggested by intermediate-scale SO(10) GUTs. Superstring-inspired models are consistent with all solutions. (4) From the narrowed parameter space, we predict the neutrino spectral shape which should be observed in SNO. Throughout the calculation we use the latest Bahcall-Pinsonneault SSM, and include nuclear and astrophysical uncertainties in a simplified, but physically transparent way.