The thermal decoupling description of dark matter (DM) and co-annihilating partners is reconsidered. If DM is realized at around the TeV-mass region or above, even the heaviest electroweak force carriers could act as long-range forces, leading to the existence of meta-stable DM bound states. The formation and subsequent decay of the latter further deplete the relic density during the freeze-out process on top of the Sommerfeld enhancement, allowing for larger DM masses. While so far the bound-state formation was described via the emission of an on-shell mediator ($W^{\pm}$, $Z$, $H$, $g$, photon or exotic), we point out that this particular process does not have to be the dominant scattering-bound state conversion channel in general. If the mediator is coupled in a direct way to any relativistic species present in the Early Universe, the bound-state formation can efficiently occur through particle scattering, where a mediator is exchanged virtually. To demonstrate that such a virtually stimulated conversion process can dominate the on-shell emission even for all temperatures, we analyze a simplified model where DM is coupled to only one relativistic species in the primordial plasma through an electroweak-scale mediator. We find that the bound-state formation cross section via particle scattering can exceed the on-shell emission by up to several orders of magnitude.

We point out that current and planned telescopes have the potential of probing annihilating Dark Matter (DM) with a mass of O(100) TeV and beyond. As a target for such searches, we propose models where DM annihilates into lighter mediators, themselves decaying into Standard Model (SM) particles. These models allow to reliably compute the energy spectra of the SM final states, and to naturally evade the unitarity bound on the DM mass. Indeed, long-lived mediators may cause an early matter-dominated phase in the evolution of the Universe and, upon decaying, dilute the density of preexisting relics thus allowing for very large DM masses. We compute this dilution in detail and provide results in a ready-to-use form. Considering for concreteness a model of dark U(1) DM, we then study both dilution and the signals at various high energy telescopes observing gamma rays, neutrinos and charged cosmic rays. This study enriches the physics case of these experiments, and opens a new observational window on heavy new physics sectors.

The unexpectedly strong 21cm absorption signal detected by the EDGES experiment suggests that the baryonic gas was colder at redshift $z\sim 17$ than predicted in the standard scenario. We discuss a mechanism to lower the baryon temperature after recombination. We introduce a stable, negatively-charged particle with a non-negligible cosmological abundance, such that the universe remains charge-neutral but the electron and proton numbers are no longer equal. The deficit of electrons during recombination results in an earlier decoupling of the baryon gas temperature from that of the CMB. This implies a smaller ratio of the gas and CMB temperature at $z\sim 17$. The parameter space of the mechanism where the 21 cm absorption signal is significantly enhanced is probed by the CMB spectrum, cooling of stars and supernovae, and colliders. Nevertheless, we find viable regions corresponding to sub-eV or MeV-scale milli-charged particles, or to TeV-scale multi-charged particles.

Dark matter (DM) charged under a dark U(1) force appears in many extensions of the Standard Model, and has been invoked to explain anomalies in cosmic-ray data, as well as a self-interacting DM candidate. In this paper, we perform a comprehensive phenomenological analysis of such a model, assuming that the DM abundance arises from the thermal freeze-out of the dark interactions. We include, for the first time, bound-state effects both in the DM production and in the indirect detection signals, and quantify their importance for Fermi, AMS, and CMB experiments. We find that DM in the mass range 1 GeV to 100 TeV, annihilating into dark photons of MeV to GeV mass, is in conflict with observations. Instead, DM annihilation into heavier dark photons is viable. We point out that the late decays of multi-GeV dark photons can produce significant entropy and thus dilute the DM density. This can lower considerably the dark coupling needed to obtain the DM abundance, and in turn relax the existing constraints.

We consider indirect detection signals of atomic dark matter, with a massive dark photon which mixes kinetically with hypercharge. In significant regions of parameter space, dark matter remains at least partially ionized today, and dark atom formation can occur efficiently in dense regions, such as the centers of galactic halos. The formation of dark atoms is accompanied by emission of a dark photon, which can subsequently decay into Standard Model particles. We discuss the expected signal strength and compare it to that of annihilating dark matter. As a case study, we explore the possibility that dark atom formation can account for the observed 511 keV line and outline the relevant parameter space.

We show that the relic abundance of thermal dark matter annihilating via a long-range interaction, is significantly affected by the formation and decay of dark matter bound states in the early universe, if the dark matter mass is above a few TeV. We determine the coupling required to obtain the observed dark matter density, taking into account both the direct 2-to-2 annihilations and the formation of bound states, and provide an analytical fit. We argue that the unitarity limit on the inelastic cross-section is realized only if dark matter annihilates via a long-range interaction, and we determine the upper bound on the mass of thermal-relic dark matter to be about 197 (139) TeV for (non)-self-conjugate dark matter.

Dark matter (DM) with sizeable self-interactions mediated by a light species offers a compelling explanation of the observed galactic substructure; furthermore, the direct coupling between DM and a light particle contributes to the DM annihilation in the early universe. If the DM abundance is due to a dark particle-antiparticle asymmetry, the DM annihilation cross-section can be arbitrarily large, and the coupling of DM to the light species can be significant. We consider the case of asymmetric DM interacting via a light (but not necessarily massless) Abelian gauge vector boson, a dark photon. In the massless dark photon limit, gauge invariance mandates that DM be multicomponent, consisting of positive and negative dark ions of different species which partially bind in neutral dark atoms. We argue that a similar conclusion holds for light dark photons; in particular, we establish that the multi-component and atomic character of DM persists in much of the parameter space where the dark photon is sufficiently light to mediate sizeable DM self-interactions. We discuss the cosmological sequence of events in this scenario, including the dark asymmetry generation, the freeze-out of annihilations, the dark recombination and the phase transition which gives mass to the dark photon. We estimate the effect of self-interactions in DM haloes, taking into account this cosmological history. We place constraints based on the observed ellipticity of large haloes, and identify the regimes where DM self-scattering can affect the dynamics of smaller haloes, bringing theory in better agreement with observations. Moreover, we estimate the cosmological abundance of dark photons in various regimes, and derive pertinent bounds.

Asymmetric dark matter models are based on the hypothesis that the present-day abundance of dark matter has the same origin as the abundance of ordinary or visible matter: an asymmetry in the number densities of particles and antiparticles. They are largely motivated by the observed similarity in the mass densities of dark and visible matter, with the former observed to be about five times the latter. This review discusses the construction of asymmetric dark matter models, summarizes cosmological and astrophysical implications and bounds, and touches on direct detection prospects and collider signatures.

It has been argued that the existence of old neutron stars excludes the possibility of non-annihilating light bosonic dark matter, such as that arising in asymmetric dark matter scenarios. If non-annihilating dark matter is captured by neutron stars, the density will eventually become sufficient for black hole formation. However, the dynamics of collapse is highly sensitive to dark-matter self-interactions. Repulsive self-interactions, even if extremely weak, can prevent black hole formation. We argue that self-interactions will necessarily be present, and estimate their strength in representative models. We also consider co-annihilation of dark matter with nucleons, which arises naturally in many asymmetric dark matter models, and which again acts to prevent black hole formation. We demonstrate how the excluded region of the dark-matter parameter space shrinks as the strength of such interactions is increased, and conclude that neutron star observations do not exclude most realistic bosonic asymmetric dark matter models.

Pangenesis is the mechanism for jointly producing the visible and dark matter asymmetries via Affleck-Dine dynamics in a baryon-symmetric universe. The baryon-symmetric feature means that the dark asymmetry cancels the visible baryon asymmetry and thus enforces a tight relationship between the visible and dark matter number densities. The purpose of this paper is to analyse the general dynamics of this scenario in more detail and to construct specific models. After reviewing the simple symmetry structure that underpins all baryon-symmetric models, we turn to a detailed analysis of the required Affleck-Dine dynamics. Both gravity-mediated and gauge-mediated supersymmetry breaking are considered, with the messenger scale left arbitrary in the latter, and the viable regions of parameter space are determined. In the gauge-mediated case where gravitinos are light and stable, the regime where they constitute a small fraction of the dark matter density is identified. We discuss the formation of Q-balls, and delineate various regimes in the parameter space of the Affleck-Dine potential with respect to their stability or lifetime and their decay modes. We outline the regions in which Q-ball formation and decay is consistent with successful pangenesis. Examples of viable dark sectors are presented, and constraints are derived from big bang nucleosynthesis, large scale structure formation and the Bullet cluster. Collider signatures and implications for direct dark matter detection experiments are briefly discussed. The following would constitute evidence for pangenesis: supersymmetry, GeV-scale dark matter mass(es) and a Z' boson with a significant invisible width into the dark sector.

The similar cosmological abundances observed for visible and dark matter suggest a common origin for both. By viewing the dark matter density as a dark-sector asymmetry, mirroring the situation in the visible sector, we show that the visible and dark matter asymmetries may have arisen simultaneously through a first-order phase transition in the early universe. The dark asymmetry can then be equal and opposite to the usual visible matter asymmetry, leading to a universe that is symmetric with respect to a generalised baryon number. We present both a general structure, and a precisely defined example of a viable model of this type. In that example, the dark matter is atomic as well as asymmetric, and various cosmological and astrophysical constraints are derived. Testable consequences for colliders include a Z' boson that couples through the B-L charge to the visible sector, but also decays invisibly to dark sector particles. The additional scalar particles in the theory can mix with the standard Higgs boson and provide other striking signatures.

The similarity of the visible and dark matter abundances indicates that they may originate via the same mechanism. If both the dark and the visible matter are charged under a generalized baryon number which remains always conserved, then the asymmetry of the visible sector may be compensated by an asymmetry in the dark sector. We show how the separation of baryonic and antibaryonic charge can originate in the vacuum, via the Affleck-Dine mechanism, due to the breaking of a symmetry orthogonal to the baryon number. Symmetry restoration in the current epoch guarantees the individual stability of the two sectors.

We calculate the neutrino signal resulting from annihilation of secluded dark matter in the Sun. In this class of models, dark matter annihilates first into metastable mediators, which subsequently decay into Standard Model particles. If the mediators are long lived, they will propagate out from the dense solar core before decaying. High energy neutrinos undergo absorption in the Sun. In the standard scenario in which neutrinos are produced directly in the centre of the Sun, absorption is relevant for E > 100 GeV, resulting in a significant suppression of the neutrino spectrum beyond E ~ 1 TeV. In the secluded dark matter scenario, the neutrino signal is greatly enhanced because neutrinos are injected away from the core, at lower density. Since the solar density falls exponentially with radius, metastable mediators have a significant effect on the neutrino flux, even for decay lengths which are small compared to the solar radius. Moreover, since neutrino detection cross sections grow with energy, this enhancement of the high energy region of the neutrino spectrum would have a large effect on overall event rates.

We consider a class of late-decaying dark-matter models, in which a dark matter particle decays to a heavy stable daughter of approximately the same mass, together with one or more relativistic particles which carry away only a small fraction of the parent rest mass. Such decays can affect galactic halo structure and evolution, and have been invoked as a remedy to some of the small scale structure-formation problems of cold dark matter. There are existing stringent limits on the dark matter lifetime if the decays produce photons. By considering examples in which the relativistic decay products instead consist of neutrinos or electron-position pairs, we derive stringent limits on these scenarios for a wide range of dark matter masses. We thus eliminate a sizable portion of the parameter space for these late decay models if the dominant decay channel involves Standard Model final states.

We consider sterile neutrinos with rest masses ~0.2 GeV. Such sterile neutrinos could augment core collapse supernova shock energies by enhancing energy transport from the core to the vicinity of the shock front. The decay of these neutrinos could produce a flux of very energetic active neutrinos, detectable by future neutrino observations from a galactic supernova. The relevant range of sterile neutrino masses and mixing angles can be probed in future laboratory experiments.

Sterile neutrino with mass of several keV can be the cosmological dark matter, can explain the observed velocities of pulsars, and can play an important role in the formation of the first stars. We describe the production of sterile neutrinos in a model with an extended Higgs sector, in which the Majorana mass term is generated by the vacuum expectation value of a gauge-singlet Higgs boson. In this model the relic abundance of sterile neutrinos does not necessarily depend on their mixing angles, the free-streaming length can be much smaller than in the case of warm dark matter produced by neutrino oscillations, and, therefore, some of the previously quoted bounds do not apply. The presence of the gauge singlet in the Higgs sector has important implications for the electroweak phase transition, baryogenesis, and the upcoming experiments at the Large Hadron Collider and a Linear Collider.

We propose a model for the acceleration of charged particles in interplanetary space that appear during quiet time periods, that is, not associated with solar activity events like intense flares or coronal mass ejections. The interaction of charged particles with modeled turbulent electromagnetic fields, which mimic the fields observed in the interplanetary medium, is studied. The turbulence is modeled by means of a dynamical system,the Gledzer-Ohkitani-Yamada (GOY) shell model, which describes the gross features of the Navier-Stokes equations. The GOY model is used to build a 3-D velocity field, which in turn is used to numerically solve the ideal magneto-hydrodynamic (MHD) induction equation, while the electric field is calculated from the ideal Ohm's law. Particle acceleration in such an environment is investigated by test particle simulations, and the resulting energy distributions are discussed and compared to observations of suprathermal electrons and ions during quiet periods in interplanetary space.