This report, based on the Dark Sectors workshop at SLAC in April 2016, summarizes the scientific importance of searches for dark sector dark matter and forces at masses beneath the weak-scale, the status of this broad international field, the important milestones motivating future exploration, and promising experimental opportunities to reach these milestones over the next 5-10 years.

The DUNE/LBNF program aims to address key questions in neutrino physics and astroparticle physics. Realizing DUNE's potential to reconstruct low-energy particles in the 10-100 MeV energy range will bring significant benefits for all DUNE's science goals. In neutrino physics, low-energy sensitivity will improve neutrino energy reconstruction in the GeV range relevant for the kinematics of DUNE's long-baseline oscillation program. In astroparticle physics, low-energy capabilities will make DUNE's far detectors the world's best apparatus for studying the electron-neutrino flux from a supernova. This will open a new window to unrivaled studies of the dynamics and neutronization of a star's central core in real time, the potential discovery of the neutrino mass hierarchy, provide new sensitivity to physics beyond the Standard Model, and evidence of neutrino quantum-coherence effects. The same capabilities will also provide new sensitivity to `boosted dark matter' models that are not observable in traditional direct dark matter detectors.

We examine a framework with light new physics, which couples to the Standard Model only via neutrino mixing. Taking the hints from the short-baseline anomalies seriously and combining them with modern cosmological data and recent IceCube measurements, we obtain surprisingly effective constraints on the hidden force: keV $\lesssim M \lesssim0.3$ GeV for the mediator mass and $g_{h}>10^{-6}-10^{-3}$ for the coupling constant. Flavor equilibration between the hidden and active neutrinos can be delayed until temperatures of $\sim 1$ MeV, but not below $\sim 100$ keV. This scenario can be tested with next-generation Cosmic Microwave Background, IceCube, and oscillation experiments.

It has been suggested that the baseline scenario of collisionless cold dark matter over-predicts the numbers of satellite galaxies, as well as the dark matter (DM) densities in galactic centers. This apparent lack of structure at small scales can be accounted for if one postulates neutrino-DM and DM-DM interactions mediated by light O(MeV) force carriers. In this letter, we consider a simple, consistent model of neutrinophilic DM with these features where DM and a "secluded" SM-singlet neutrino species are charged under a new $U(1)$ gauge symmetry. An important ingredient of this model is that the secluded sector couples to the Standard Model fields only through neutrino mixing. We observe that the secluded and active neutrinos recouple, leading to a large relic secluded neutrino population. This relic population can prevent small-scale halos from collapsing, while at same time significantly modifying the optical depth of ultra-high-energy neutrinos recently observed at Icecube. We find that the bulk of the parameter space accommodating an (a)symmetric thermal relic has potentially observable consequences for the IceCube high energy signal, with some of the parameter ranges already ruled out by the existing data. Future data may confirm this mechanism if either spectral absorption features or correlations with nearby sources are observed.

Accretion disks arising from neutron star- neutron star mergers or black hole- neutron star mergers produce large numbers of neutrinos and antineutrinos. In contrast to other astrophysical scenarios, like supernovae, in mergers the antineutrinos outnumber the neutrinos. This antineutrino dominance gives neutrinos from merger disks the opportunity to exhibit new oscillation physics, specifically a matter-neutrino resonance. We explore this resonance, finding that consequences can be a large transition of $\nu_e$ to other flavors, while the $\bar{\nu}_e$s return to their initial state. We present numerical calculations of neutrinos from merger disks and compare with a single energy model. We explain both the basic features and the conditions for a transition.

We present a brief overview of the ongoing searches for the axion particle via its coupling to photons. Both the classical QCD axions and more recently proposed Axion-Like-Particles are considered. Astrophysical bounds on the axion-photon coupling come from considerations of stellar energy loss during Helium burning, in both low- and high-mass stars. Helioscopes look for back-conversion of solar axions into x-ray photons in strong laboratory magnetic fields. Finally, haloscopes aim to detect dark matter axions in our galactic halo. Both types of searches are expecting significant advances in the future, which will enable them to probe large, well-motivated parts of the parameter space below the stellar cooling bounds.

We study the capabilities of the MAJORANA DEMONSTRATOR, a neutrinoless double-beta decay experiment currently under construction at the Sanford Underground Laboratory, as a light WIMP detector. For a cross section near the current experimental bound, the MAJORANA DEMONSTRATOR should collect hundreds or even thousands of recoil events. This opens up the possibility of simultaneously determining the physical properties of the dark matter and its local velocity distribution, directly from the data. We analyze this possibility and find that allowing the dark matter velocity distribution to float considerably worsens the WIMP mass determination. This result is traced to a previously unexplored degeneracy between the WIMP mass and the velocity dispersion. We simulate spectra using both isothermal and Via Lactea II velocity distributions and comment on the possible impact of streams. We conclude that knowledge of the dark matter velocity distribution will greatly facilitate the mass and cross section determination for a light WIMP.

We give the first self-consistent calculation of the effect of the scattered neutrino halo on flavor evolution in supernovae. Our example case is an O-Ne-Mg core collapse supernova neutronization neutrino burst. We find that the addition of the halo neutrinos produces qualitative and quantitative changes in the final flavor states of neutrinos. We also find that the halo neutrinos produce a novel distortion of the neutrino flavor swap. Our results provide strong motivation for tackling the full multidimensional and composition-dependent aspects of this problem in the future.

Isotropy is a key assumption in many models of cosmic-ray electrons and positrons. We find that simulation results imply a critical energy of ~10-1000 GeV above which electrons and positrons can spend their entire lives in streams threading magnetic fields, due to energy losses. This would restrict the number of electron/positron sources contributing at Earth, likely leading to smooth electron and positron spectra, as is observed. For positrons, this could be as few as one, with an enhanced flux that would ease energetics concerns of a pulsar origin of the positron excess, or even zero, bringing dark matter into play. We conclude that ideas about electron/positron propagation based on either isotropic diffusion or turbulent fields must be changed.

We point out that stars in the mass window ~ 8-12 Msun can serve as sensitive probes of the axion-photon interaction, g_{A\gamma\gamma}. Specifically, for these stars axion energy losses from the helium-burning core would shorten and eventually eliminate the blue loop phase of the evolution. This would contradict observational data, since the blue loops are required, e.g., to account for the existence of Cepheid stars. Using the MESA stellar evolution code, modified to include the extra cooling, we conservatively find g_{A\gamma\gamma} <~ 0.8 * 10^{-10} GeV^{-1}, which compares favorably with the existing bounds.

We argue that the small fraction of neutrinos that undergo direction-changing scattering outside of the neutrinosphere could have significant influence on neutrino flavor transformation in core-collapse supernova environments. We show that the standard treatment for collective neutrino flavor transformation is adequate at late times, but could be inadequate in the crucial shock revival/explosion epoch of core-collapse supernovae, where the potentials that govern neutrino flavor evolution are affected by the scattered neutrinos. Taking account of this effect, and the way it couples to entropy and composition, will require a new paradigm in supernova modeling.

We investigate collective flavor oscillations of supernova neutrinos at late stages of the explosion. We first show that the frequently used single-angle (averaged coupling) approximation predicts oscillations close to, or perhaps even inside, the neutrinosphere, potentially invalidating the basic neutrino transport paradigm. Fortunately, we also find that the single-angle approximation breaks down in this regime; in the full multiangle calculation, the oscillations start safely outside the transport region. The new suppression effect is traced to the interplay between the dispersion in the neutrino-neutrino interactions and the vacuum oscillation term.

Cosmological and astrophysical measurements provide powerful constraints on neutrino masses complementary to those from accelerators and reactors. Here we provide a guide to these different probes, for each explaining its physical basis, underlying assumptions, current and future reach.

Recently, it has been demonstrated that neutrinos in a supernova oscillate collectively. This process occurs much deeper than the conventional matter-induced MSW effect and hence may have an impact on nucleosynthesis. In this paper we explore the effects of collective neutrino oscillations on the r-process, using representative late-time neutrino spectra and outflow models. We find that accurate modeling of the collective oscillations is essential for this analysis. As an illustration, the often-used "single-angle" approximation makes grossly inaccurate predictions for the yields in our setup. With the proper multiangle treatment, the effect of the oscillations is found to be less dramatic, but still significant. Since the oscillation patterns are sensitive to the details of the emitted fluxes and the sign of the neutrino mass hierarchy, so are the r-process yields. The magnitude of the effect also depends sensitively on the astrophysical conditions - in particular on the interplay between the time when nuclei begin to exist in significant numbers and the time when the collective oscillation begins. A more definitive understanding of the astrophysical conditions, and accurate modeling of the collective oscillations for those conditions, is necessary.

Neutrinos in a core-collapse supernova undergo coherent flavor transformations in their own background. We explore this phenomenon during the cooling stage of the explosion. Our three-flavor calculations reveal qualitatively new effects compared to a two-flavor analysis. These effects are especially clearly seen for the inverted mass hierarchy: we find a different pattern of spectral "swaps" in the neutrino spectrum and a novel "mixed" spectrum for the antineutrinos. A brief discussion of the relevant physics is presented, including the instability of the two-flavor evolution trajectory, the 3-flavor pattern of spectral "swaps," and partial nonadiabaticity of the evolution.

We explore the sensitivity of massive stars to neutrino magnetic moments. We find that the additional cooling due to the neutrino magnetic moments bring about qualitative changes to the structure and evolution of stars in the mass window 7 Msun < M < 18 Msun, rather than simply changing the time scales for the burning. We describe some of the consequences of this modified evolution: the shifts in the threshold masses for creating core-collapse supernovae and oxygen-neon-magnesium white dwarfs and the appearance of a new type of supernova in which a partial carbon-oxygen core explodes within a massive star. The resulting sensitivity to the magnetic moment is at the level of (2-4) * 10^{-11} \mu_B.

In a class of extra dimensional models with a warped metric and a single brane the photon can be localized on the brane by gravity only. An intriguing feature of these models is the possibility of the photon escaping into the extra dimensions. The search for this effect has motivated the present round of precision orthopositronium decay experiments. We point out that in this framework a photon in plasma should be metastable. We consider the astrophysical consequences of this observation, in particular, what it implies for the plasmon decay rate in globular cluster stars and for the core-collapse supernova cooling rate. The resulting bounds on the model parameter exceed the possible reach of orthopositronium experiments by many orders of magnitude.

In several classes of particle physics models -- ranging from the classical Majoron models, to the more recent scenarios of late neutrino masses or Mass-Varying Neutrinos -- one or more of the neutrinos are postulated to couple to a new light scalar field. As a result of this coupling, neutrinos in the early universe instead of streaming freely could form a self-coupled fluid, with potentially observable signatures in the Cosmic Microwave Background and the large scale structure of the universe. We re-examine the constraints on this scenario from the presently available cosmological data and investigate the sensitivity expected from the Planck satellite. In the first case, we find that the sensitivity strongly depends on which piece of data is used. The SDSS Main sample data, combined with WMAP and other data, disfavors the scenario of three coupled neutrinos at about the 3.5$\sigma$ confidence level, but also favors a high number of freely streaming neutrinos, with the best fit at 5.2. If the matter power spectrum is instead taken from the SDSS Large Red Galaxy sample, best fit point has 2.5 freely streaming neutrinos, but the scenario with three coupled neutrinos becomes allowed at $2\sigma$. In contrast, Planck alone will exclude even a single self-coupled neutrino at the $4.2\sigma$ confidence level, and will determine the total radiation at CMB epoch to $\Delta N_\nu^{eff} = ^{+0.5}_{-0.3}$ ($1\sigma$ errors). We investigate the robustness of this result with respect to the details of Planck's detector. This sensitivity to neutrino free-streaming implies that Planck will be capable of probing a large region of the Mass-Varying Neutrino parameter space. Planck may also be sensitive to a scale of neutrino mass generation as high as 1 TeV.

Convection that develops behind the shock front during the first second of a core-collapse supernova explosion is believed to play a crucial role in the explosion mechanism. We demonstrate that the resulting turbulent density fluctuations may be directly observable in the neutrino signal starting at t>~ 3-4 s after the onset of the explosion. The effect comes from the modulation of the MSW flavor transformations by the turbulent density fluctuations. We derive a simple and general criterion for neutrino flavor depolarization in a Kolmogorov-type turbulence and apply it to the turbulence seen in modern numerical simulations. The turbulence casts a "shadow", by making other features, such as the shock front, unobservable in the density range covered by the turbulence.

It is of great interest whether the recent KamLAND bound on the flux of electron antineutrinos from the Sun constrains the electromagnetic properties of the neutrino. We examine the efficiency of the electron antineutrino production in the solar magnetic fields, assuming the neutrinos are Majorana particles with a relatively large transition moment. We consider fields both in the radiative and convective zones of the Sun, with physically plausible strengths, and take into account the recently established values of the oscillation parameters. Our analysis shows that the production rate in question is presently unobservable. In the radiative zone, it is suppressed by the large measured value of the flavor mixing angle which eliminates the resonant level crossing. A corresponding general resonance condition, valid for large as well as small values of the mixing angle, is derived. Likewise, in the convective zone, the strength of the small-scale magnetic field is likely insufficient. Thus, no useful bound on the neutrino transition moment can be derived from the published KamLAND bound. KamLAND may be, however, on the edge of probing an ``optimistic'' scenario, making further improvements of its sensitivity desirable.

The possibility that the energy density of the Universe is dominated by a network of low-tension domain walls provides an alternative to the commonly discussed cosmological constant and scalar-field quintessence models of dark energy. We quantify the lower bound on the number density of the domain walls that follows from the observed near-isotropy of the cosmic microwave background radiation. This bound can be satisfied by a strongly frustrated domain wall network. No fine-tuning of the parameters of the underlying field theory model is required. We briefly outline the observational consequences of this model.

We reexamine the transition magnetic moment solution to the solar neutrino problem. We argue that the absence of large time variations in the Super-Kamiokande rate provides strong evidence against spin-flavor flip in the solar convective zone. Spin-flavor flip could, however, occur in the primordial magnetic field in the radiative zone. We compute the longest lived toroidal mode for this field and show that spin-flavor flip in the radiative zone can account for all available solar data.

We discuss bounds on the strength of the magnetic fields that could be buried in the radiative zone of the Sun. The field profiles and decay times are computed for all axisymmetric toroidal Ohmic decay eigenmodes with lifetimes exceeding the age of the Sun. The measurements of the solar oblateness yield a bound <~ 7 MG on the strength of the field. A comparable bound is expected to come from the analysis of the splitting of the solar oscillation frequencies. The theoretical analysis of the double diffusive instability also yields a similar bound. The oblateness measurements at their present level of sensitivity are therefore not expected to measure a toroidal field contribution.

The traditional physical description of neutrino flavor conversion in the Sun focuses on the notion of resonance. However, the resonance picture is valid only in the limit of small mixing angles theta. For large values of theta, the resonance picture leads to seemingly paradoxical results. This observation is important for understanding the physics of neutrino flavor conversion in the Sun, since the latest solar neutrino data seems to prefer large mixing angles. Here we review the basic arguments and in particular show that the resonance does not in general coincide with either the point of maximal violation of adiabaticity in the nonadiabatic case or the point of maximal flavor conversion in the adiabatic case. We also discuss a modified adiabaticity criterion.

Recently, several interesting proposals were made modifying the law of gravity on large scales, within a sensible relativistic formulation. This allows a precise formulation of the idea that such a modification might account for galaxy rotation curves, instead of the usual interpretation of these curves as evidence for dark matter. We here summarize several observational constraints which any such modification must satisfy, and which we believe make more challenging any interpretation of galaxy rotation curves in terms of new gravitational physics.

We point out that for solar neutrino oscillations with the mass-squared difference of Delta m^2 ~ 10^-10 - 10^-9 eV^2, traditionally known as ``vacuum oscillation'' range, the solar matter effects are non-negligible, particularly for the low energy pp neutrinos. One consequence of this is that the values of the mixing angle theta and pi/2-theta are not equivalent, leading to the need to consider the entire physical range of the mixing angle 0<=theta<=pi/2 when determining the allowed values of the neutrino oscillation parameters.

Results of neutrino oscillation experiments have always been presented on the (sin^2 2theta, Delta m^2) parameter space for the case of two-flavor oscillations. We point out, however, that this parameterization misses the half of the parameter space pi/4 < theta <= pi/2 (``the dark side''), which is physically inequivalent to the region 0 <= theta <= pi/4 (``the light side'') in the presence of matter effects. The MSW solutions to the solar neutrino problem can extend to the dark side, especially if we take the conservative attitude to allow higher confidence levels, ignore some of the experimental results in the fits, or relax theoretical predictions. Furthermore even the so-called ``vacuum oscillation'' solution distinguishes the dark and the light sides. We urge experimental collaborations to present their results on the entire parameter space.

Measuring the 7Be solar neutrino flux is crucial towards solving the solar neutrino puzzle. The Borexino experiment, and possibly the KamLAND experiment, will be capable of studying the 7Be neutrinos in the near future. We discuss (1) how the seasonal variation of the Borexino and KamLAND data can be used to measure the 7Be solar neutrino flux in a background independent way and (2) how anomalous seasonal variations might be used to discover vacuum neutrino oscillations, independent of the solar model and the measurement of the background. In particular, we find that, after three years of Borexino or KamLAND running, vacuum neutrino oscillations can be either established or excluded for almost all values of (sin^2 2 theta, Delta m^2) preferred by the Homestake, GALLEX, SAGE, and Super-Kamiokande data. We also discuss how well seasonal variations of the data can be used to measure (sin^2 2 theta, Delta m^2) in the case of vacuum oscillations.