We describe a scenario where the smallness of neutrino masses is related to a global symmetry that is only violated by quantum gravitational effects. The coupling of neutrinos to gauge singlet right-handed fermions is attributed to symmetry preserving gravitational operators suppressed by the Planck mass, in this framework. The proposed scenario leads to axion particles that decay into neutrinos, which could be probed through cosmological measurements and may help explain the Hubble parameter tension. Depending on the details of the implementation, the scenario could provide axion dark matter candidates.

Producing an acceptable thermal relic abundance of dark matter with masses $\gg 10^2$ TeV is a challenge. We propose a novel mechanism where GeV-scale states establish a tiny thermal relic abundance for dark matter, which is later promoted to ultra massive status by a very light scalar. We refer to this dark matter as a THermal Ultra Massive Particle (THUMP). Direct detection of THUMPs can be naturally expected due to large scattering cross sections mediated by low mass states that couple THUMPs to the Standard Model. Our model generically leads to signals for the associated GeV-scale states at accelerator experiments.

The initial data from the Event Horizon Telescope (EHT) on M87$^*$, the supermassive black hole at the center of the M87 galaxy, provide direct observational information on its mass, spin, and accretion disk properties. A combination of the EHT data and other constraints provide evidence that M87$^*$ has a mass $\sim 6.5 \times 10^9\,M_\odot$ and dimensionless spin parameter $|a^*|\gtrsim 0.5$. These determinations disfavor ultra light bosons of mass $\mu_b\sim 10^{-21}$ eV, within the range considered for fuzzy dark matter, invoked to explain dark matter distribution on $\sim$ kpc scales. Future observations of M87$^*$ could be expected to strengthen our conclusions.

We propose that $O(10 M_{\rm Sun})$ black holes observed by LIGO/Virgo originate from a first order QCD phase transition at a temperature $T_* \lesssim 100$~MeV. This is realized by keeping the quark masses small compared to confinement scale down to $T\sim T_*$, making QCD transition first order. We implement this scenario using a light scalar that could potentially be a good dark matter candidate.

We consider a long range scalar force that mainly couples to dark matter and unstable Standard Model states, like the muon, with tiny strength. Probing this type of force would present a challenge to observations. We point out that the dependence of the induced background scalar field on dark matter number density can cause the mass of the unstable particles to have spatial and temporal variations. These variations, in turn, leave an imprint on the value of the fine structure constant $\alpha$, through threshold corrections, that could be detected in astronomical and cosmological measurements. Our mechanism can accommodate the mild preference of the Planck data for such a deviation, $(\alpha_{_{\rm CMB}}-\alpha_{\rm present})/\alpha_{\rm present} = (-3.6\pm 3.7)\times 10^{-3}$. In this case, the requisite parameters typically imply that violations of Equivalence Principle may be within reach of future experiments.

We propose that neutrino masses can be zero $in ~vacuo$ and may be generated by the local distribution of dark matter through a feeble long range scalar force. We discuss potential phenomenological constraints and implications of this framework. Our model typically implies that the relic neutrino background left over from the Big Bang is mostly absent in our Galactic neighborhood. Hence, a positive detection signal from future proposed experiments, such as PTOLEMY, could in principle falsify our scenario.

We consider a feeble repulsive interaction between ordinary matter and dark matter, with a range similar to or larger than the size of the Earth. Dark matter can thus be repelled from the Earth, leading to null results in direct detection experiments, regardless of the strength of the short-distance interactions of dark matter with atoms. Generically, such a repulsive force would not allow trapping of dark matter inside astronomical bodies. In this scenario, accelerator-based experiments may furnish the only robust signals of asymmetric dark matter models, which typically lack indirect signals from self-annihilation. Some of the variants of our hypothesis are also briefly discussed.

A very light boson of mass $\mathcal{O}(10^{-22})$ eV may potentially be a viable dark matter (DM) candidate which can avoid phenomenological problems associated with cold DM. Such "fuzzy DM (FDM)" may naturally be an axion with a decay constant $f_a \sim 10^{16} \div 10^{18}$ GeV, and a mass $m_a \sim \mu^2/f_a$ with $\mu\sim 10^2$ eV. Here we propose a concrete model where $\mu$ arises as a dynamical scale from infrared confining dynamics, analogous to QCD. Our model is an alternative to the usual approach of generating $\mu$ through string theoretic instanton effects. We outline the features of this scenario that result from various cosmological constraints. We find that those constraints are suggestive of a period of mild of inflation, perhaps from a strong first order phase transition, that reheats the Standard Model (SM) sector only. A typical prediction of our scenario, broadly speaking, is a larger effective number of neutrinos compared to the SM value $N_{\text{eff}} \approx 3$, as inferred from precision measurements of the cosmic microwave background. Some of the new degrees of freedom may be identified as "sterile neutrinos," which may be required to explain certain neutrino oscillation anomalies. Hence, aspects of our scenario could be testable in terrestrial experiments, which is a novelty of our FDM model.

We entertain the possibility that primordialblack holes of mass $\sim (10^{26}$--$10^{29})$~g, with Schwarzschild radii of $\mathcal{O}{\text{cm}}$, constitute $\sim 10\%$ or more of cosmic dark matter, as allowed by various constraints. These black holes would typically originate from cosmological eras corresponding to temperatures $\mathcal{O}{10-100}$~GeV, and may be associated with first order phase transitions in the visible or hidden sectors. In case these small primordial black holes get captured in orbits around neutron stars or astrophysical black holes in our galactic neighborhood, gravitational waves from the resulting "David and Goliath (D\&G)" binaries could be detectable at Advanced LIGO or Advanced Virgo for hours or more, possibly over distances of $\mathcal{O}{10}$~Mpc encompassing the Local Supercluster of galaxies. The proposed Einstein Telescope would further expand the reach for these signals. A positive signal could be further corroborated by the discovery of new particles in the $\mathcal{O}{10-100}$~GeV mass range, and potentially also the detection of long wavelength gravitational waves originating from the first order phase transition era.

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.

We describe a general scenario, dubbed "Inflatable Dark Matter", in which the density of dark matter particles can be reduced through a short period of late-time inflation in the early universe. The overproduction of dark matter that is predicted within many otherwise well-motivated models of new physics can be elegantly remedied within this context, without the need to tune underlying parameters or to appeal to anthropic considerations. Thermal relics that would otherwise be disfavored can easily be accommodated within this class of scenarios, including dark matter candidates that are very heavy or very light. Furthermore, the non-thermal abundance of GUT or Planck scale axions can be brought to acceptable levels, without invoking anthropic tuning of initial conditions. A period of late-time inflation could have occurred over a wide range of scales from ~ MeV to the weak scale or above, and could have been triggered by physics within a hidden sector, with small but not necessarily negligible couplings to the Standard Model.

Dark sectors, consisting of new, light, weakly-coupled particles that do not interact with the known strong, weak, or electromagnetic forces, are a particularly compelling possibility for new physics. Nature may contain numerous dark sectors, each with their own beautiful structure, distinct particles, and forces. This review summarizes the physics motivation for dark sectors and the exciting opportunities for experimental exploration. It is the summary of the Intensity Frontier subgroup "New, Light, Weakly-coupled Particles" of the Community Summer Study 2013 (Snowmass). We discuss axions, which solve the strong CP problem and are an excellent dark matter candidate, and their generalization to axion-like particles. We also review dark photons and other dark-sector particles, including sub-GeV dark matter, which are theoretically natural, provide for dark matter candidates or new dark matter interactions, and could resolve outstanding puzzles in particle and astro-particle physics. In many cases, the exploration of dark sectors can proceed with existing facilities and comparatively modest experiments. A rich, diverse, and low-cost experimental program has been identified that has the potential for one or more game-changing discoveries. These physics opportunities should be vigorously pursued in the US and elsewhere.

We examine the possibility that dark matter may be the manifestation of dark forces of a hidden sector, i.e. "Dark Force = Dark Matter." As an illustrative and minimal example we consider the hidden SU(2)_h x U(1)_h gauge group. The hidden dynamics is indirectly coupled to the Standard Model (SM) through kinetic mixing of U(1)_h with the U(1)_Y of hypercharge. We assume a hidden symmetry breaking pattern analogous to that of the SM electroweak symmetry, augmented with an extra scalar that allows both the "hidden Z boson" Z_h and the "hidden photon" \gamma_h to be massive. The "hidden W" bosons W_h are dark matter in this scenario. This setup can readily accommodate a potential direct detection signal for dark matter at ~10 GeV from CDMSII-Si data. For some choices of parameters, the model can lead to signals both in "dark matter beam" experiments, from Z_h\to W_h W_h, as well as in experiments that look for visible signals of dark photons, mediated by \gamma_h. Other possible phenomenological consequences are also briefly discussed.

The "gravitational baryogenesis" scenario is extended to generate both baryon and dark matter asymmetries, in the matter dominated era corresponding to post-inflationary reheating. A minimal extension requires a singlet fermion X for dark matter and a singlet scalar S. With two or more hidden sector fermions, the scenario can lead to nucleon decay into dark matter with a lifetime of order 10^{34-36} yr, which is relevant for current or future experiments. The correct multi-component relic density can be obtained if dark matter fermions couple to a sub-GeV vector boson that weakly interacts with the Standard Model through mixing. The typical inflationary scale in the scenario is of order 10^{16} GeV which suggests that tensor mode perturbations could potentially be within observational reach.

We investigate new and unusual signals that arise in theories where dark matter is asymmetric and carries a net antibaryon number, as may occur when the dark matter abundance is linked to the baryon abundance. Antibaryonic dark matter can cause {\it induced nucleon decay} by annihilating visible baryons through inelastic scattering. These processes lead to an effective nucleon lifetime of 10^{29}-10^{32} years in terrestrial nucleon decay experiments, if baryon number transfer between visible and dark sectors arises through new physics at the weak scale. The possibility of induced nucleon decay motivates a novel approach for direct detection of cosmic dark matter in nucleon decay experiments. Monojet searches (and related signatures) at hadron colliders also provide a complementary probe of weak-scale dark-matter--induced baryon number violation. Finally, we discuss the effects of baryon-destroying dark matter on stellar systems and show that it can be consistent with existing observations.

We present a novel mechanism for generating both the baryon and dark matter densities of the Universe. A new Dirac fermion X carrying a conserved baryon number charge couples to the Standard Model quarks as well as a GeV-scale hidden sector. CP-violating decays of X, produced non-thermally in low-temperature reheating, sequester antibaryon number in the hidden sector, thereby leaving a baryon excess in the visible sector. The antibaryonic hidden states are stable dark matter. A spectacular signature of this mechanism is the baryon-destroying inelastic scattering of dark matter that can annihilate baryons at appreciable rates relevant for nucleon decay searches.

We examine under what circumstances the INTEGRAL/SPI 511 keV signal can originate from decays of MeV-scale composite states produced by: (A) thermonuclear (type Ia) or (B) core collapse supernovae (SNe). The requisite dynamical properties that would account for the observed data are quite distinct, for cases (A) and (B). We determine these requirements in simple hidden valley models, where the escape fraction problem is naturally addressed, due to the long lifetime of the new composite states. A novel feature of scenario (A) is that the dynamics of type Ia SNe, standard candles for cosmological measurements, might be affected by our mechanism. In case (A), the mass of the state mediating between the hidden sector and the SM $\ee$ could be a few hundred GeV and within the reach of a 500 GeV $\ee$ linear collider. We also note that kinetic mixing of the photon with a light vector state may provide an interesting alternate mediation mechanism in this case. Scenarios based on case (B) are challenged by the need for a mechanism to transport some of the produced positrons toward the Galactic bulge, due to the inferred distribution of core collapse sources. The mass of the mediator in case (B) is typically hundreds of TeV, leading to long-lived particles that could, under certain circumstances, include a viable dark matter candidate. The appearance of long-lived particles in typical models leads to cosmological constraints and we address how a consistent cosmic history may be achieved.

Two general problems arise when interpreting the recent cosmic ray data as signals of Dark Matter (DM) annihilation: (i) the required cross section is too large by O(100), and (ii) the annihilation products seem to be mostly leptonic. We propose to address these two problems by assuming that the couplings of DM to leptons grow with time. This can be achieved by a dynamic localization of DM in extra dimensions. A possible outcome of this proposal is a time (red-shift) dependent annihilation signal, in terms of strength and dominant final states.

We present a detailed computation of the expected rate for Geomagnetic Conversion of Solar Axions to X-rays (GECOSAX) along the orbit of an x-ray satellite. We use realistic satellite orbits and propagation in time. A realistic model for the Earth's magnetic field, which properly accounts for its spatial non-uniformity, is used. We also account for the effect of the Earth's atmosphere on the propagation of x-rays in our calculation of axion-photon conversion probability. To estimate possible sensitivities to the axion-photon coupling g_{a\gamma}, we use an actual measurement of the expected backgrounds by the SUZAKU satellite. Assuming a detector area of 10^3 cm^2 and about 10^6 s of data, we show that a 2 \sigma limit of g_{a\gamma} < (4.7-6.6) times 10^{-11} GeV^{-1} from GECOSAX is achievable, for axion masses m_a<10^{-4} eV. This significantly exceeds current laboratory sensitivities to g_{a\gamma}.

We estimate the production rate of axion-type particles in the core of the Earth, at a temperature T~5000K. We constrain thermal geo-axion emission by demanding a core-cooling rate less than 100K/Gyr, as suggested by geophysics. This yields a "quasi-vacuum" (unaffected by extreme stellar conditions) bound on the axion-electron fine structure constant \alpha_a^{QV} < 10^{-18}, stronger than the existing accelerator (vacuum) bound by 4 orders of magnitude. We consider the prospects for measuring the geo-axion flux through conversion into photons in a geoscope; such measurements can further constrain \alpha_a^{QV}.

It has recently been suggested that a scale invariant "unparticle" sector with a non-trivial infrared fixed point may couple to the Standard Model (SM) via higher dimensional operators. The weakness of such interactions hides the the unparticle phenomena at low energies. We demonstrate how cosmology and astrophysics can place significant bounds on the strength of unparticle-SM interactions. We also discuss the possibility of a having a non-negligible unparticle relic density today.

Theoretical axion models state that axions are very weakly interacting particles. In order to experimentally detect them, the use of colorful and inspired techniques becomes mandatory. There is a wide variety of experimental approaches that were developed during the last 30 years, most of them make use of the Primakoff effect, by which axions convert into photons in the presence of an electromagnetic field. We review the experimental techniques used to search for axions and will give an outlook on experiments planned for the near future.

We consider cosmological consequences of the spontaneous breaking of a global symmetry that is anomalous under technicolor interactions, leading to the emergence of a light axion-like particle. Avoiding overclosure of the universe by such axions yields the upper bound $f_a \lsim 10^{10}$ GeV on the symmetry breaking scale, corresponding to keV-scale axions. However, diffuse x-ray background data typically require larger values of $f_a$. The overclosure and x-ray bounds can be reconciled if the axion initial amplitude of oscillations $A_i \sim f_a/10$. In this case, a viable axionic dark matter candidate with a mass in the 50-100 eV range emerges. The detection of this type of dark matter may pose a challenge.

We show that solar axion conversion to photons in the Earth's magnetosphere can produce an x-ray flux, with average energy \sim 4 keV, which is measurable on the dark side of the Earth. The smallness of the Earth's magnetic field is compensated by a large magnetized volume. For axion masses < 10^{-4} eV, a low-Earth-orbit x-ray detector with an effective area of 10^4 cm^2, pointed at the solar core, can probe the photon-axion coupling down to 10^{-11} GeV^{-1}, in one year. Thus, the sensitivity of this new approach will be an order of magnitude beyond current laboratory limits.

We study type II supernova signatures of neutrino mass generation via symmetry breaking at a scale in the range from keV to MeV. The scalar responsible for symmetry breaking can be thermalized in the supernova core and restore the symmetry. The neutrinos from scalar decays have about half the average energy of thermal neutrinos and are Bose-Einstein distributed. We find that, even without a detailed knowledge of the supernova parameters, a discovery is well within reach at Super-Kamiokande.

We construct the New Minimal Standard Model that incorporates the new discoveries of physics beyond the Minimal Standard Model (MSM): Dark Energy, non-baryonic Dark Matter, neutrino masses, as well as baryon asymmetry and cosmic inflation, adopting the principle of minimal particle content and the most general renormalizable Lagrangian. We base the model purely on empirical facts rather than aesthetics. We need only six new degrees of freedom beyond the MSM. It is free from excessive flavor-changing effects, CP violation, too-rapid proton decay, problems with electroweak precision data, and unwanted cosmological relics. Any model of physics beyond the MSM should be measured against the phenomenological success of this model.

We show that a gravitational interaction between the derivative of the Ricci scalar curvature and the baryon-number current dynamically breaks CPT in an expanding universe and, combined with baryon-number-violating interactions, can drive the universe towards an equilibrium baryon asymmetry that is observationally acceptable.

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.

The phenomenological possibilities of the Randall-Sundrum non-compact extra dimension scenario with the AdS horizon increased to approximately a millimeter length, corresponding to an effective brane tension of TeV^4, are investigated. The corrections to the Newtonian potential are found to be the only observationally accessible probe of this scenario, as previously suggested in the literature. In particular, the presence of the continuum of KK modes does not lead to any observable collider signatures. The extent to which experimental tests of Newtonian gravity can distinguish this scenario from the scenario of Arkani-Hamed, Dimopoulos, and Dvali with one and two millimeter size extra dimensions is explicitly demonstrated.

The flux of cosmic rays beyond the GZK cutoff ($\sim 10^{20}$ eV) may be explained through their production by ultra high energy cosmic neutrinos, annihilating on the relic neutrino background, in the vicinity of our galaxy. This process is mediated through the production of a $Z$ boson at resonance, and is generally known as the $Z$-Burst mechanism. We show that a similar mechanism can also contribute to the super-GZK spectrum at even higher, ultra-GZK energies, where the particles produced at resonance are the Kaluza-Klein gravitons of weak scale mass and coupling from the Randall-Sundrum (RS) hierarchy model of localized gravity model. We call this mechanism Gravi-Burst. We discuss the parameter space of relevance to Gravi-Bursts, and comment on the possibility of its contribution to the present and future super-GZK cosmic ray data and place bounds on the RS model parameters. Under certain assumptions about the energy spectrum of the primary neutrinos we find that cosmic ray data could be potentially as powerful as the LHC in probing the RS model.

We describe a new mechanism for the generation of the baryon asymmetry of the universe during a first order electroweak phase transition. The mechanism requires the existence of two (or more) baryon number carrying scalar fields with masses and CP violating mixing which vary with the Higgs field expectation value. This mechanism can be implemented using squarks in supersymmetric theories or using leptoquarks. Our central observation is that reflection of these scalars from a bubble wall can yield a significant net baryon number flux into the symmetric phase, balanced by a flux of opposite sign into the broken phase. For generic parameter choices, scalars with incident energies in a specific, but not narrow, range yield order one reflection asymmetries (between the probability of reflection of the scalars and of their antiparticles). The interesting energies are those for which there are two propagating scalars in the symmetric phase but only one in the broken phase. Electroweak sphaleron processes drive the baryon number in the symmetric phase toward zero, but do not act in the broken phase. Our estimate of the resulting baryon asymmetry is consistent with cosmological observations for a range of mass parameters and CP violating phases in a supersymmetric implementation, as long as the bubble walls are not too fast and not too thick.