The existence of dark matter particles that carry phenomenologically relevant self-interaction cross sections mediated by light dark sector states is considered to be severely constrained through a combination of experimental and observational data. The conclusion is based on the assumption of specific dark matter production mechanisms such as thermal freeze-out together with an extrapolation of a standard cosmological history beyond the epoch of primordial nucleosynthesis. In this work, we drop these assumptions and examine the scenario from the perspective of the current firm knowledge we have: results from direct and indirect dark matter searches and cosmological and astrophysical observations, without additional assumptions on dark matter genesis or the thermal state of the very early universe. We show that even in the minimal set-up, where dark matter particles self-interact via a kinetically mixed vector mediator, a significant amount of parameter space remains allowed. Interestingly, however, these parameter regions imply a meta-stable, light mediator, which in turn calls for modified search strategies.

Photons traveling cosmological distances through the inhomogeneous Universe experience a great variation in their in-medium induced effective mass. Using the EAGLE suite of hydrodynamical simulations, we infer the free electron distribution and thereby the effective photon mass after reionization. We use this data to study the inter-conversion of kinetically mixed photons and dark photons, which may occur at a great number of resonance redshifts, and obtain the "optical depth" against conversion along random lines-of-sight. Using COBE/FIRAS, Planck, and SPT measurements, we constrain the dark photon parameter space from the depletion of CMB photons into dark photons that causes both spectral distortions and additional anisotropies in the CMB. Finally, we also consider the conversion of sub-eV dark radiation into ordinary photons. We make the line-of-sight distributions of both, free electrons and dark matter, publicly available.

In the inhomogeneous Universe, the cosmological conversion of dark photons into ordinary photons (and vice versa) may happen at a great number of resonance redshifts. This alters the CMB observed energy spectrum and degree of small-scale anisotropies. We utilize results from the EAGLE simulation to obtain the conversion probability along random line-of-sights to quantify these effects. We then apply our results to the case where dark photons are sourced by dark matter decay and their high-redshift conversion into ordinary photons modify the global 21 cm signal expected from the cosmic dawn era. Concretely, we show that a significant portion of the parameter space for which a converted population of photons in the Rayleigh-Jeans tail of the CMB explains the absorption strength observed by EDGES, is ruled out from the brightness temperature measurements of COBE/FIRAS and the CMB anisotropy measurements of Planck and SPT.

Electromagnetically neutral dark sector particles may directly couple to the photon through higher dimensional effective operators. Considering electric and magnetic dipole moment, anapole moment, and charge radius interactions, we derive constraints from stellar energy loss in the Sun, horizontal branch and red giant stars, as well as from cooling of the proto-neutron star of SN1987A. We provide the exact formula for in-medium photon-mediated pair production to leading order in the dark coupling, and compute the energy loss rates explicitly for the most important processes, including a careful discussion on resonances and potential double counting between the processes. Stringent limits for dark states with masses below $3\,$keV ($40\,$MeV) arise from red giant stars (SN1987A), implying an effective lower mass-scale of approximately $10^9\,$GeV ($10^7\,$GeV) for mass-dimension five, and $100\,$GeV ($2.5\,$TeV) for mass-dimension six operators as long as dark states stream freely; for the proto-neutron star, the trapping of dark states is also evaluated. Together with direct limits previously derived by us in Chu et al. (2018), this provides the first comprehensive overview of the viability of effective electromagnetic dark-state interactions below the GeV mass-scale.

New particles $\chi$ that are electrically neutral but couple to the electromagnetic current via higher-dimensional operators and that are sufficiently light, at or below the GeV-mass scale, can be produced in pairs in a number of dedicated high-intensity experiments. In this work we consider the production of $\chi$ through magnetic- and electric-dipole moments as well as through anapole moment and charge radius interactions in electron beams. We derive new constraints from BaBar, NA64 and mQ and forecast the future sensitivity on the existence of such states, from Belle-II, LDMX and BDX. We present for the first time a detailed treatment of the off-shell production of photons in electron beams with subsequent decay into a $\chi\bar\chi$ pair in a 2-to-4 process. These direct limits are then compared to the effects on SM precision observables, as well as to bounds from flavor physics and high energy colliders. Finally, we consider the possibility that $\chi$ is dark matter and study ensuing astrophysical and cosmological constraints. We find that a combination of all considered probes rule out $\chi$ particles with mass-dimension five and six photon interactions as dark matter when assuming a standard freeze-out abundance.

Precision cosmology provides a sensitive probe of extremely weakly coupled states due to thermal freeze-in production, with subsequent decays impacting physics during well-tested cosmological epochs. We explore the cosmological implications of the freeze-in production of a new scalar $S$ via the super-renormalizable Higgs portal. If the mass of $S$ is at or below the electroweak scale, peak freeze-in production occurs during the electroweak epoch. We improve the calculation of the freeze-in abundance by including all relevant QCD and electroweak production channels. The resulting abundance and subsequent decay of $S$ is constrained by a combination of X-ray data, cosmic microwave background anisotropies and spectral distortions, $N_{\rm eff}$, and the consistency of BBN with observations. These probes constrain technically natural couplings for such scalars from $m_S \sim$ keV all the way to $m_S \sim 100$ GeV. The ensuing constraints are similar in spirit to typical beam bump limits, but extend to much smaller couplings, down to mixing angles as small as $\theta_{Sh} \sim 10^{-16}$, and to masses all the way to the electroweak scale.

We show that despite stringent constraints on the shape of the main part of the CMB spectrum, there is considerable room for its modification within its Rayleigh-Jeans (RJ) end, $\omega \ll T_{\rm CMB}$. We construct explicit New Physics models that give an order one (or larger) increase of photon count in the RJ tail, which can be tested by existing and upcoming experiments aiming to detect the cosmological 21 cm emission/absorption signal. This class of models stipulates the decay of unstable particles to dark photons, $A'$, that have a small mass, $m_{A'} \sim 10^{-14} - 10^{-9}$ eV, non-vanishing mixing angle $\epsilon$ with electromagnetism, and energies much smaller than $T_{\rm CMB}$. The non-thermal number density of dark photons can be many orders of magnitude above the number density of CMB photons, and even a small probability of $A'\to A$ oscillations, for values of as small as $\epsilon \sim 10^{-9}$, can significantly increase the number of RJ photons. In particular, we show that resonant oscillations of dark photons into regular photons in the interval of redshifts $20 < z < 1700$ can be invoked as an explanation of the recent tentative observation of a stronger-than-expected absorption signal of 21 cm photons. We present a model that realizes this possibility, where milli-eV mass dark matter decays to dark photons, with a lifetime longer than the age of the Universe.

We consider the generic possibility that the Universe's energy budget includes some form of relativistic or semi-relativistic dark radiation (DR) with non-gravitational interactions with Standard Model (SM) particles. Such dark radiation may consist of SM singlets or a non-thermal, energetic component of neutrinos. If such DR is created at a relatively recent epoch, it can carry sufficient energy to leave a detectable imprint in experiments designed to search for very weakly interacting particles: dark matter and underground neutrino experiments. We analyze this possibility in some generality, assuming that the interactive dark radiation is sourced by late decays of an unstable particle, potentially a component of dark matter, and considering a variety of possible interactions between the dark radiation and SM particles. Concentrating on the sub-GeV energy region, we derive constraints on different forms of DR using the results of the most sensitive neutrino and dark matter direct detection experiments. In particular, for interacting dark radiation carrying a typical momentum of $\sim30$~MeV$/c$, both types of experiments provide competitive constraints. This study also demonstrates that non-standard sources of neutrino emission (e.g. via dark matter decay) are capable of creating a "neutrino floor" for dark matter direct detection that is closer to current bounds than is expected from standard neutrino sources.

This white paper summarizes the workshop "U.S. Cosmic Visions: New Ideas in Dark Matter" held at University of Maryland on March 23-25, 2017.

Inspired by the string axiverse idea, it has been suggested that the recent transition from decelerated to accelerated cosmic expansion is driven by an axion-like quintessence field with a sub-Planckian decay constant. The scenario requires that the axion field be rather near the maximum of its potential, but is less finely tuned than other explanations of cosmic acceleration. The model is parametrized by an axion decay constant $f$, the axion mass $m$, and an initial misalignment angle $|\theta_i|$ which is close to $\pi$. In order to determine the $m$ and $\theta_{i}$ values consistent with observations, these parameters are mapped onto observables: the Hubble parameter $H(z)$ at and angular diameter distance $d_{A}(z)$ to redshift $z= 0.57$, as well as the angular sound horizon of the cosmic microwave background (CMB). Measurements of the baryon acoustic oscillation (BAO) scale at $z\simeq 0.57$ by the BOSS survey and Planck measurements of CMB temperature anisotropies are then used to probe the $\left\{m,f,\theta_i\right\}$ parameter space. With current data, CMB constraints are the most powerful, allowing a fraction of only $\sim 0.2$ of the parameter-space volume. Measurements of the BAO scale made using the SPHEREx or SKA experiments could go further, observationally distinguishing all but $\sim 10^{-2}$ or $\sim 10^{-5}$ of the parameter-space volume (allowed by simple priors) from the $\Lambda$CDM model.

We point out that the cosmological abundance of ${}^7$Li can be reduced down to observed values if during its formation Big Bang Nucleosynthesis is modified by the presence of light electrically neutral particles $X$ that have substantial interactions with nucleons. We find that the lithium problem can be solved without affecting the precisely measured abundances of deuterium and helium if the following conditions are satisfied: the mass and lifetimes of such particles are bounded by $ 1.6~{\rm MeV}\leq m_X \leq 20~{\rm MeV}$ and $ {\rm few}~100~{\rm s} \lesssim \tau_X \lesssim 10^4~{\rm s}$, and the abundance times the absorption cross section by either deuterium or ${}^7$Be are comparable to the Hubble rate, $n_X \sigma_{\rm abs} v \sim H$, at the time of ${}^7$Be formation. We include $X$-initiated reactions into the primordial nucleosynthesis framework, observe that it leads to a substantial reduction of the freeze-out abundances of ${}^7$Li+${}^7$Be, and find specific model realizations of this scenario. Concentrating on the axion-like-particle case, $X=a$, we show that all these conditions can be satisifed if the coupling to $d$-quarks is in the range of $f_d^{-1} \sim {\rm TeV}^{-1}$, which can be probed at intensity frontier experiments.

New Abelian U(1)' gauge bosons $V_{\mu}$ can couple to the Standard Model through mixing of the associated field strength tensor $V_{\mu\nu}$ with the one from hypercharge, $F_{\mu\nu}^Y$. Here we consider early Universe sensitivity to this vector portal and show that the effective mixing parameter with the photon, $\kappa$, is being probed for vector masses in the GeV ballpark down to values $10^{-10} \lesssim \kappa \lesssim 10^{-14}$ where no terrestrial probes exist. The ensuing constraints are based on a detailed calculation of the vector relic abundance and an in-depth analysis of relevant nucleosynthesis processes.

String theories suggest the existence of a plethora of axion-like fields with masses spread over a huge number of decades. Here we show that these ideas lend themselves to a model of quintessence with no super-Planckian field excursions and in which all dimensionless numbers are order unity. The scenario addresses the "why now" problem, i.e., why has accelerated expansion begun only recently, by suggesting that the onset of dark-energy domination occurs randomly with a slowly decreasing probability per unit logarithmic interval in cosmic time. The standard axion potential requires us to postulate a rapid decay of most of the axion fields that do no become dark energy. The need for these decays is averted, though, with the introduction of a slightly modified axion potential. In either case, a Universe like ours arises in roughly 1 in 100 universes. The scenario may have a host of observable consequences.

We explore the cosmological consequences of kinetically mixed dark photons with a mass between 1 MeV and 10 GeV, and an effective electromagnetic fine structure constant as small as $10^{-38}$. We calculate the freeze-in abundance of these dark photons in the early Universe and explore the impact of late decays on BBN and the CMB. This leads to new constraints on the parameter space of mass $m_V$ vs kinetic mixing parameter $\kappa$.

We study the dissipation of small-scale adiabatic perturbations at early times when the Universe is hotter than T ~ 0.5 keV. When the wavelength falls below the damping scale 1/kD, the acoustic modes diffuse and thermalize, causing entropy production. Before neutrino decoupling, kD is primarily set by the neutrino shear viscosity, and we study the effect of acoustic damping on the relic neutrino number, primordial nucleosynthesis, dark-matter freeze-out, and baryogenesis. This sets a new limit on the amplitude of primordial fluctuations of DeltaR^2 < 0.007 at 10^4/Mpc< k < 10^5/Mpc and a model dependent limit of DeltaR^2 < 0.3 at k < 10^{20-25}/Mpc.

Light new particles with masses below 10 keV, often considered as a plausible extension of the Standard Model, will be emitted from the solar interior, and can be detected on the Earth with a variety of experimental tools. Here we analyze the new "dark" vector state V, a massive vector boson mixed with the photon via an angle kappa, that in the limit of the small mass m_V has its emission spectrum strongly peaked at low energies. Thus, we utilize the constraints on the atomic ionization rate imposed by the results of the XENON10 experiment to set the limit on the parameters of this model: kappa times m_V< 3 times10^{-12} eV. This makes low-threshold Dark Matter experiments the most sensitive dark vector helioscopes, as our result not only improves current experimental bounds from other searches by several orders of magnitude, but also surpasses even the most stringent astrophysical and cosmological limits in a seven-decade-wide interval of m_V. We generalize this approach to other light exotic particles, and set the most stringent direct constraints on "mini-charged" particles.

We consider the stellar production of vector states V within the minimal model of "dark photons". We show that when the Stuckelberg mass of the dark vector becomes smaller than plasma frequency, the emission rate is dominated by the production of the longitudinal modes of V, and scales as \kappa^2 m_V^2, where \kappa and m_V are the mixing angle with the photon and the mass of the dark state. This is in contrast with earlier erroneous claims in the literature that the emission rate decouples as the forth power of the mass. We derive ensuing constraints on the (\kappa, m_V) parameter space by calculating the cooling rates for the Sun and red giant stars. We find that stellar bounds for m_V < 10 eV are significantly strengthened, to the extent that all current "light-shining-through-wall" experiments find themselves within deeply excluded regions.

The rate of the direct decay of 40K to the ground state of 40Ar through electron capture has not been experimentally reported. Aside from its inherent importance for the theory of electron capture as the only such decay known of its type (unique third-forbidden), this decay presents an irreducible background in the DAMA experiment. We find that the presence of this background, as well as others, poses a challenge to any interpretation of the DAMA results in terms of a Dark Matter model with a small modulation fraction. A 10ppb contamination of natural potassium requires a 20% modulation fraction or more. A 20ppb contamination, which is reported as an upper limit by DAMA, disfavors any Dark Matter origin of the signal. This conclusion is based on the efficiency of detecting 40K decays as inferred from simulation. We propose measures to help clarify the situation.

We reply to the critiques of our paper arXiv:1210.5501 by the DAMA collaboration which appeared in arXiv:1210.6199 and arXiv:1211.6346. Our original claim that the observed background levels are likely to require a large modulation fraction of any putative signal holds. In fact, in light of DAMA's recent comment our claim is further corroborated. We identify the source of the discrepancy between our own analysis and DAMA's claimed levels of unmodulated background. Our analysis indicates that the background in the signal region as reported by DAMA is indeed likely underestimated.

If the neutral component of weak-scale dark matter is accompanied by a charged excitation separated by a mass gap of less than ~20 MeV, WIMPs can form stable bound states with nuclei. We show that the recent progress in experiments searching for neutrinoless double-beta decay sets the first direct constraint on the exoergic reaction of WIMP-nucleus bound state formation. We calculate the rate for such process in representative models and show that the double-beta decay experiments provide unique sensitivity to a large fraction of parameter space of the WIMP doublet model, complementary to constraints imposed by cosmology and direct collider searches.

The DAMA dark matter search experiment observes a statistically significant percent-level variation of its low-energy count rate with a period of one year. In this note we recall some of the arguments which challenge the hypothesis that the cosmic ray induced underground muon flux can be the origin of the modulation. In addition, we provide new comments on recent works on this subject.

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

The aim of the current work is a detailed time-series analysis of the data from Dark Matter direct detection experiments as well as related datasets. We examine recent claims that the cosmic ray muon flux can be responsible for generating the modulation signals seen by DAMA and, more recently, by the CoGeNT collaboration. We find no evidence for such a strong correlation and show that the two phenomena differ in their power spectrum, phase, and possibly in amplitude. In addition, we investigate in more detail, the time dependence of Dark Matter scattering. Since the signal is periodic with period of a year (due to the Earth's motion around the Sun), the presence of higher harmonics can be expected. We show that the higher harmonics generically have similar phase to the annual modulation and the biannual mode in particular could provide another handle in searching for Dark Matter in the laboratory.

Many models of new physics including variants of supersymmetry predict metastable long-lived particles that can decay during or after primordial nucleosynthesis, releasing significant amounts of non-thermal energy. The hadronic energy injection in these decays leads to the formation of ^9Be via the chain of non-equilibrium transformations: Energy_h -> T, ^3He -> ^6He, ^6Li -> ^9Be. We calculate the efficiency of this transformation and show that if the injection happens at cosmic times of a few hours, the release of 10 MeV per baryon can be sufficient for obtaining a sizable ^9Be abundance. The absence of a plateau-structure in the ^9Be/H abundance down to a 10^{-14} level allows one to use beryllium as a robust constraint on new physics models with decaying or annihilating particles.

We consider the minimal model in which dark matter is stabilized by a non-Abelian discrete symmetry. The symmetry group is taken to be D_3, which is the smallest non-Abelian finite group. The minimal model contains (nontrivial) singlet and doublet scalar representations of D_3 which couple to the Standard Model fields via the Higgs portal. This construction predicts two species of dark matter over much of the parameter space. Nontrivial interactions under D_3 lead to a novel thermal history of dark matter, while the multi-component nature of dark matter can be tested by future direct detection experiments.

Big bang nucleosynthesis (BBN), an epoch of primordial nuclear transformations in the expanding Universe, has left an observable imprint in the abundances of light elements. Precision observations of such abundances, combined with high-accuracy predictions, provide a nontrivial test of the hot big bang and probe non-standard cosmological and particle physics scenarios. We give an overview of BBN sensitivity to different classes of new physics: new particle or field degrees of freedom, time-varying couplings, decaying or annihilating massive particles leading to non-thermal processes, and catalysis of BBN by charged relics.

The persistent discrepancy between observations of 7Li with putative primordial origin and its abundance prediction in Big Bang Nucleosynthesis (BBN) has become a challenge for the standard cosmological and astrophysical picture. We point out that the decay of GeV-scale metastable particles X may significantly reduce the BBN value down to a level at which it is reconciled with observations. The most efficient reduction occurs when the decay happens to charged pions and kaons, followed by their charge exchange reactions with protons. Similarly, if X decays to muons, secondary electron antineutrinos produce a similar effect. We consider the viability of these mechanisms in different classes of new GeV-scale sectors, and find that several minimal extensions of the Standard Model with metastable vector and/or scalar particles are capable of solving the cosmological lithium problem. Such light states can be a key to the explanation of recent cosmic ray anomalies and can be searched for in a variety of high-intensity medium-energy experiments.

The catalysis of nuclear reactions by negatively charged relics leads to increased outputs of primordial ^6Li and ^9Be. In combination with observational constraints on the primordial fractions of ^6Li and ^9Be, this imposes strong restrictions on the primordial abundance and the lifetime of charged relics. We analyze the constraints from the catalysis of ^9Be on supersymmetric models in which the gravitino is the lightest supersymmetric particle and a charged slepton--such as the lighter stau--the next-to-lightest supersymmetric particle (NLSP). Barring the special cases in which the primordial fraction of the slepton NLSP is significantly depleted, we find that the ^9Be data require a slepton NLSP lifetime of less than 6x10^3 seconds. We also address the issue of the catalytic destruction of ^6Li and ^9Be by late forming bound states of protons with negatively charged relics finding that it does not lead to any significant modification of the limit on the slepton lifetime.

In supersymmetric models with a long-lived stau being the lightest Standard Model superpartner, the stau abundance during primordial nucleosynthesis is tightly constrained. Considering the complete set of stau annihilation channels in the minimal supersymmetric Standard Model (MSSM) with real parameters for scenarios in which sparticle coannihilations are negligible, we calculate the decoupling of the lighter stau from the primordial plasma and identify processes which are capable to deplete the resulting stau abundance significantly. We find particularly efficient stau annihilation at the resonance of the heavy CP-even Higgs boson and for a lighter stau with a sizeable left-right mixing due to enhanced stau-Higgs couplings. Even within the constrained MSSM, we encounter both effects leading to exceptionally small values of the resulting stau abundance. Prospects for collider phenomenology are discussed and possible implications of our findings are addressed with emphasis on gravitino dark matter scenarios.

We investigate gravitino dark matter scenarios in which the primordial ^6Li production is catalyzed by bound-state formation of long-lived negatively charged particles X^- with ^4He. In the constrained minimal supersymmetric Standard Model (CMSSM) with the stau^- as the X^-, the observationally inferred bound on the primordial ^6Li abundance allows us to derive a rigid lower limit on the gaugino mass parameter for a standard cosmological history. This limit can have severe implications for supersymmetry searches at the Large Hadron Collider and for the reheating temperature after inflation.

Considering gravitino dark matter scenarios, we obtain the full gauge-invariant result for the relic density of thermally produced gravitinos to leading order in the Standard Model gauge couplings. For the temperatures required by thermal leptogenesis, we find gaugino mass bounds which will be probed at future colliders. We show that a conceivable determination of the gravitino mass will allow for a unique test of the viability of thermal leptogenesis in the laboratory.

Considering gravitino dark matter scenarios, we study constraints on the reheating temperature of inflation. We present the gauge-invariant result for the thermally produced gravitino yield to leading order in the Standard Model gauge couplings. Within the framework of the constrained minimal supersymmetric Standard Model (CMSSM), we find a maximum reheating temperature of about 10^7 GeV taking into account bound-state effects on the primordial $^6$Li abundance. We show that late-time entropy production can relax this constraint significantly. Only with a substantial entropy release after the decoupling of the lightest Standard Model superpartner, thermal leptogenesis remains a viable explanation of the cosmic baryon asymmetry within the CMSSM.