We consider the impact of quantum diffusion on inflationary dynamics during an ultra-slow-roll phase, which can be of particular significance for the formation of primordial black holes. We show, by means of a fully analytical approach, that the power spectrum of comoving curvature perturbations computed in stochastic inflation matches precisely, at the linear level, the result obtained by solving the Mukhanov-Sasaki equation, even in the presence of an ultra-slow-roll phase. We confirm this result numerically in a model in which the inflaton has a polynomial potential and is coupled quadratically to the Ricci scalar. En route, we assess the role that quantum noise plays in the presence of an ultra-slow-roll phase, and clarify the issue of the quantum-to-classical transition in this scenario.

We consider the possibility that the majority of dark matter in our Universe consists of black holes of primordial origin. We determine the conditions under which such black holes may have originated from a single-field model of inflation characterized by a quartic polynomial potential. We also explore the effect of higher-dimensional operators. The large power spectrum of curvature perturbations that is needed for a large black hole abundance sources sizable second order tensor perturbations. The resulting stochastic background of primordial gravitational waves could be detected by the future space-based observatories LISA and DECIGO or --as long as we give up on the dark matter connection--by the ground-based Advanced LIGO-Virgo detector network.

Core-collapse supernovae presumably explode because trapped neutrinos push the material out of the stellar envelope. This process is directly controlled by the weak scale $v$: we argue that supernova explosions happen only if fundamental constants are tuned within a factor of few as $v \sim \Lambda_{\rm QCD}^{3/4} M_{\rm Pl}^{1/4}$, such that neutrinos are trapped in supernovae for a time comparable to the gravitational time-scale. We provide analytic arguments and simulations in spherical approximation, that need to be validated by more comprehensive simulations. The above result can be important for fundamental physics, because core-collapse supernova explosions seem anthropically needed, as they spread intermediate-mass nuclei presumably necessary for `life'. We also study stellar burning, finding that it does not provide anthropic boundaries on $v$.

At the dawn of a golden age for gravitational wave astronomy, we must leave no stone unturned in our quest for new phenomena beyond our current understanding of General Relativity (GR), particle physics and nuclear physics. In this paper we discuss gravitational echoes from ultracompact stars. We restrict our analysis to exact solutions of Einstein field equations in GR that are supported by physically motivated equations of state (EoS), and in particular we impose the constraint of causality. Our main conclusion is that ultracompact objects supported by physical EoS are not able to generate gravitational echoes like those that characterize the relaxation phase of a putative black hole mimicker. Nevertheless, we identify a class of physical exotic objects that are compact enough to accommodate the presence of an external unstable light ring, thus opening the possibility of trapping gravitational radiation and affecting the ringdown phase of a merger event. Most importantly, we show that once rotation is included these stars -- contrary to what usually expected for ultracompact objects -- are not plagued by any ergoregion instability. We extend our analysis for arbitrary values of angular velocity up to the Keplerian limit, and we comment about potential signals relevant for gravitational wave interferometers.

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 critically reexamine two possible Dark Matter candidate within the Standard Model. First, we consider the $uuddss$ exa-quark. Its QCD binding energy could be large enough to make it (quasi) stable. We show that the cosmological Dark Matter abundance is reproduced thermally if its mass is 1.2 GeV. However, we also find that such mass is excluded by the stability of Oxygen nuclei. Second, we consider the possibility that the instability in the Higgs potential leads to the formation of primordial black holes while avoiding vacuum decay during inflation. We show that the non-minimal Higgs coupling to gravity must be as small as allowed by quantum corrections, $|\xi_H| < 0.01$. Even so, one must assume that the Universe survived in $e^{120}$ independent regions to fluctuations that lead to vacuum decay with probability 1/2 each.

An inhomogeneous pseudo-scalar field configuration behaves like an optically active medium. Consequently, if a light ray passes through an axion cloud surrounding a Kerr black hole, it may experience a polarization-dependent bending. We explore the size and relevance of such effect considering both the QCD axion and a generic axion-like particle.

We present a novel interpretation of the gamma-ray diffuse emission measured by H.E.S.S. in the Galactic Center (GC) region and the Galactic ridge. Our starting base is an updated analysis of PASS8 Fermi-LAT data, which allows to extend down to few GeV the spectra measured by H.E.S.S. and to infer the primary CR radial distribution above 100 GeV. We compare those results with a CR transport model assuming a harder scaling of the diffusion coefficient with rigidity in the inner Galaxy. Such a behavior reproduces the radial dependence of the CR spectral index recently inferred from Fermi-LAT measurements in the inner GP. We find that, in this scenario, the bulk of the Galactic ridge emission can be naturally explained by the interaction of the diffuse, steady-state Galactic CR sea interacting with the gas present in the Central molecular zone. The evidence of a GC PeVatron is significantly weaker than that inferred adopting a conventional (softer) CR sea.

There exist non-trivial stationary points of the Euclidean action for an axion particle minimally coupled to Einstein gravity, dubbed wormholes. They explicitly break the continuos global shift symmetry of the axion in a non-perturbative way, and generate an effective potential that may compete with QCD depending on the value of the axion decay constant. In this paper, we explore both theoretical and phenomenological aspects of this issue. On the theory side, we address the problem of stability of the wormhole solutions, and we show that the spectrum of the quadratic action features only positive eigenvalues. On the phenomenological side, we discuss, beside the obvious application to the QCD axion, relevant consequences for models with ultralight dark matter, black hole superradiance, and the relaxation of the electroweak scale. We conclude discussing wormhole solutions for a generic coset and the potential they generate.

We present a novel interpretation of the $\gamma$-ray diffuse emission measured by Fermi-LAT and H.E.S.S. in the Galactic center (GC) region and the Galactic ridge (GR). In the first part we perform a data-driven analysis based on PASS8 Fermi-LAT data: we extend down to few GeV the spectra measured by H.E.S.S. and infer the primary cosmic-ray (CR) radial distribution between 0.1 and 3 TeV. In the second part we adopt a CR transport model based on a position-dependent diffusion coefficient. Such behavior reproduces the radial dependence of the CR spectral index recently inferred from the Fermi-LAT observations. We find that the bulk of the GR emission can be naturally explained by the interaction of the diffuse steady-state Galactic CR sea with the gas present in the Central Molecular Zone. Although our results leave room for a residual radial-dependent emission associated with a central source, the relevance of the large-scale background prevents from a solid evidence of a GC Pevatron.

The LIGO observation of gravitational waves from a binary black hole merger has begun a new era in fundamental physics. If new dark sector particles, be they bosons or fermions, can coalesce into exotic compact objects (ECOs) of astronomical size, then the first evidence for such objects, and their underlying microphysical description, may arise in gravitational wave observations. In this work we study how the macroscopic properties of ECOs are related to their microscopic properties, such as dark particle mass and couplings. We then demonstrate the smoking gun exotic signatures that would provide observational evidence for ECOs, and hence new particles, in terrestrial gravitational wave observatories. Finally, we discuss how gravitational waves can test a core concept in general relativity: Hawking's area theorem.

The last IceCube catalog of High Energy Starting Events (HESE) obtained with a livetime of 1347 days comprises 54 neutrino events equally-distributed between the three families with energies between 25 TeV and few PeVs. Considering the homogeneous flavors distribution (1:1:1) and the spectral features of these neutrinos the IceCube collaboration claims the astrophysical origin of these events with more than $5\sigma$. The spatial distribution of cited events does not show a clear correlation with known astrophysical accelerators leaving opened both the Galactic and the extra-Galactic origin interpretations. Here, we compute the neutrino diffuse emission of our Galaxy on the basis of a recently proposed phenomenological model characterized by radially-dependent cosmic-ray (CR) transport properties. We show that the astrophysical spectrum measured by IceCube experiment can be well explained adding to the diffuse Galactic neutrino flux (obtained with this new model) a extra-Galactic component derived from the astrophysical muonic neutrinos reconstructed in the Northern hemisphere. A good agreement between the expected astrophysical neutrino flux and the IceCube data is found for the full sky as well as for the Galactic plane region.

Recent years have seen increased theoretical and experimental effort towards the first-ever detection of cosmic-ray antideuterons, in particular as an indirect signature of dark matter annihilation or decay. In contrast to indirect dark matter searches using positrons, antiprotons, or gamma-rays, which suffer from relatively high and uncertain astrophysical backgrounds, searches with antideuterons benefit from very suppressed conventional backgrounds, offering a potential breakthrough in unexplored phase space for dark matter. This article is based on the first dedicated cosmic-ray antideuteron workshop, which was held at UCLA in June 2014. It reviews broad classes of dark matter candidates that result in detectable cosmic-ray antideuteron fluxes, as well as the status and prospects of current experimental searches. The coalescence model of antideuteron production and the influence of antideuteron measurements at particle colliders are discussed. This is followed by a review of the modeling of antideuteron propagation through the magnetic fields, plasma currents, and molecular material of our Galaxy, the solar system, the Earth's geomagnetic field, and the atmosphere. Finally, the three ongoing or planned experiments that are sensitive to cosmic-ray antideuterons, BESS, AMS-02, and GAPS, are detailed. As cosmic-ray antideuteron detection is a rare event search, multiple experiments with orthogonal techniques and backgrounds are essential. Many theoretical and experimental groups have contributed to these studies over the last decade, this review aims to provide the first coherent discussion of the relevant dark matter theories that antideuterons probe, the challenges to predictions and interpretations of antideuteron signals, and the experimental efforts toward cosmic antideuteron detection.

A spherical-symmetric gamma-ray emission from the central region of the Galaxy has been recently identified in Fermi-LAT data, and initially associated to dark matter particle annihilations. Guided by the evidence for a high gas density in the inner kpc of the Galaxy correlated with a very large Supernova rate, and hence with ongoing cosmic-ray acceleration, we investigate instead the possibility of addressing this excess in terms of ordinary cosmic-ray sources and standard steady-state diffusion. We alter the source term, and consistently the correlated gamma-ray emissions, in the context of a template-fitting analysis. We focus on a region of interest (ROI) defined as: $|l|<20^{\circ}; \hspace{5mm} 2^{\circ}<|b|<20^{\circ},$ with $l$ and $b$ the Galactic longitude and latitude coordinates. We analyze in detail the overall goodness of the fit of our framework, and perform a detailed direct comparison against data examining profiles in different directions. Remarkably, the test statistic of the fit related to our scenario turns out to be as good as the Dark Matter one in the ROI here considered.

We compute the gamma-ray and neutrino diffuse emission of the Galaxy on the basis of a recently proposed phenomenological model characterized by radially dependent cosmic-ray (CR) transport properties. We show how this model, designed to reproduce both Fermi-LAT gamma-ray data and local CR observables, naturally reproduces the anomalous TeV diffuse emission observed by Milagro in the inner Galactic plane. Above 100 TeV our picture predicts a neutrino flux that is about five (two) times larger than the neutrino flux computed with conventional models in the Galactic Center region (full-sky). Explaining in that way up to $\sim 25 \%$ of the flux measured by IceCube, we reproduce the full-sky IceCube spectrum adding an extra-Galactic component derived from the muonic neutrino flux in the northern hemisphere. We also present precise predictions for the Galactic plane region where the flux is dominated by the Galactic emission.

We show how the Higgs boson mass is protected from the potentially large corrections due to the introduction of minimal dark matter if the new physics sector is made supersymmetric. The fermionic dark matter candidate (a 5-plet of $SU(2)_L$) is accompanied by a scalar state. The weak gauge sector is made supersymmetric and the Higgs boson is embedded in a supersymmetric multiplet. The remaining standard model states are non-supersymmetric. Non vanishing corrections to the Higgs boson mass only appear at three-loop level and the model is natural for dark matter masses up to 15 TeV--a value larger than the one required by the cosmological relic density. The construction presented stands as an example of a general approach to naturalness that solves the little hierarchy problem which arises when new physics is added beyond the standard model at an energy scale around 10 TeV.

Several observations from Fermi-LAT, up to few hundred GeV, and from H.E.S.S., up to $\sim$ 10 TeV, reported an intense $\gamma$-ray emission from the inner part of the Galactic plane. After the subtraction of point-like contributions, the remaining $\gamma$-ray spectrum can provide important hints about the cosmic-ray (CR) population in that region. In particular, the diffuse spectrum measured by both Fermi-LAT and H.E.S.S. in the Galactic Ridge is significantly harder with respect to the rest of the Galaxy. These results were recently interpreted in terms of a comprehensive CR transport model which, adopting a spatial dependent diffusion coefficient and convective velocity, reproduces Fermi-LAT results on the whole sky as well as local CR spectra. We showed as that model predicts a significantly harder neutrino diffuse emission compared to conventional scenarios: The predicted signal is able to account for a significant fraction of the astrophysical flux measured by IceCube. In this contribution, we use the same setup to calculate the expected neutrino flux from several windows in the inner Galactic plane and compare the results with IceCube observations and the sensitivities of Mediterranean neutrino telescopes. In particular, for the ANTARES experiment, we compare the model expectations with the upper limits obtained from a recent unblinded data-analysis focused on the galactic ridge region. Moreover, we also show the expectations from the galactic ridge for the future KM3NeT observatory, whose position is optimal to observe this portion of the sky.

As recently shown, Fermi-LAT measurements of the diffuse gamma-ray emission from the Galaxy favor the presence of a smooth softening in the primary cosmic-ray spectrum with increasing Galactocentric distance. This result can be interpreted in terms of a spatial-dependent rigidity scaling of the diffusion coefficient. The DRAGON code was used to build a model based on such feature. That scenario correctly reproduces the latest Fermi-LAT results as well as local cosmic-ray measurements from PAMELA, AMS-02 and CREAM. Here we show that the model, if extrapolated at larger energies, grasps both the gamma-ray flux measured by MILAGRO at 15 TeV and the H.E.S.S. data from the Galactic ridge, assuming that the cosmic-ray spectral hardening found by those experiments at about 250 GeV/n is present in the whole inner Galactic plane region. Moreover, we show as that model also predicts a neutrino emission which may account for a significant fraction, as well as for the correct spectral shape, of the astrophysical flux measured by IceCube above 25 TeV.

The precise measurement of cosmic-ray antiparticles serves as important means for identifying the nature of dark matter. Recent years showed that identifying the nature of dark matter with cosmic-ray positrons and higher energy antiprotons is difficult, and has lead to a significantly increased interest in cosmic-ray antideuteron searches. Antideuterons may also be generated in dark matter annihilations or decays, offering a potential breakthrough in unexplored phase space for dark matter. Low-energy antideuterons are an important approach because the flux from dark matter interactions exceeds the background flux by more than two orders of magnitude in the low-energy range for a wide variety of models. This review is based on the "dbar14 - dedicated cosmic-ray antideuteron workshop", which brought together theorists and experimentalists in the field to discuss the current status, perspectives, and challenges for cosmic-ray antideuteron searches and discusses the motivation for antideuteron searches, the theoretical and experimental uncertainties of antideuteron production and propagation in our Galaxy, as well as give an experimental cosmic-ray antideuteron search status update. This report is a condensed summary of the article "Review of the theoretical and experimental status of dark matter identification with cosmic-ray antideuteron" (arXiv:1505.07785).

The standard approach to cosmic-ray (CR) propagation in the Galaxy is based on the assumption that local transport properties can be extrapolated to the whole CR confining volume. Such models tend to underestimate the gamma-ray flux above few GeV measured by the Fermi Large Area Telescope towards the inner Galactic plane. We consider here for the first time a phenomenological scenario allowing for both the rigidity scaling of the diffusion coefficient and convective effects to be position-dependent. We show that within this approach we can reproduce the observed gamma-ray spectra at both low and mid Galactic latitudes - including the Galactic center - without spoiling any local CR observable.

In this paper we analyze a simple scenario in which Dark Matter (DM) consists of free fermions with mass $m_f$. We assume that on galactic scales these fermions are capable of forming a degenerate Fermi gas, in which stability against gravitational collapse is ensured by the Pauli exclusion principle. The mass density of the resulting configuration is governed by a non-relativistic Lane-Emden equation, thus leading to a universal cored profile that depends only on one free parameter in addition to $m_f$. After reviewing the basic formalism, we test this scenario against experimental data describing the velocity dispersion of the eight classical dwarf spheroidal galaxies of the Milky Way. We find that, despite its extreme simplicity, the model exhibits a good fit to the data and realistic predictions for the size of DM halos providing that $m_f\simeq 200$ eV. Furthermore, we show that in this setup larger galaxies correspond to the non-degenerate limit of the gas. We propose a concrete realization of this model in which DM is produced non-thermally via inflaton decay. We show that imposing the correct relic abundance and the bound on the free-streaming length constrains the inflation model in terms of inflaton mass, its branching ratio into DM and the reheating temperature.

The scalar Higgs portal is a compelling model of dark matter (DM) in which a renormalizable coupling with the Higgs boson provides the connection between the visible world and the dark sector. In this paper we investigate the constraint placed on the parameter space of this model by the antiproton data. Due to the fact that the antiproton-to-proton ratio has relative less systematic uncertainties than the antiproton absolute flux, we propose and explore the possibility to combine all the available $\bar{p}/p$ data. Following this approach, we are able to obtain stronger limits if compared with the existing literature. In particular, we show that most of the parameter space close to the Higgs resonance is ruled out by our analysis. Furthermore, by studying the reach of the future AMS-02 antiproton and antideuteron data, we argue that a DM mass of $\mathcal{O}(150)$ GeV offers a promising discovery potential. The method of combining all the antiproton-to-proton ratio data proposed in this paper is quite general, and can be straightforwardly applied to other models.

A GeV gamma-ray excess has possibly been individuated in Fermi-LAT data from the Galactic Center and interpreted in terms of Dark Matter (DM) annihilations, either in hadronic (essentially $b\bar{b}$) or leptonic channels. In order to test this tantalizing interpretation, we address two issues: (i) we improve the computation of secondary emission from DM (Inverse Compton and Bremsstrahlung) with respect to previous works, confirming it to be very relevant for determining the DM spectrum in the leptonic channels, so that any conclusion on the DM nature of the signal critically depends on this contribution; (ii) we consider the constraints from antiprotons on the DM hadronic channel, finding that the uncertainties on the propagation model, and in particular on the halo height, play a major role. Moreover, we discuss the role of solar modulation, taking into account possible charge dependent effects whose importance is estimated exploiting detailed numerical tools. The limits that we obtain severely constrain the DM interpretation of the excess in the hadronic channel, for standard assumptions on the Galactic propagation parameters and solar modulation. However, they considerably relax if more conservative choices are adopted.

The actual realization of the electroweak symmetry breaking in the context of a natural extension of the Standard Model (SM) and the nature of Dark Matter (DM) are two of the most compelling questions in high-energy particle physics. Composite Higgs models may provide a unified picture in which both the Higgs boson and the DM particle arise as pseudo Nambu-Goldstone bosons of a spontaneously broken global symmetry at a scale $f\sim$ TeV. In this paper we analyze a general class of these models based on the coset SO(6)/SO(5). Assuming the existence of light and weakly coupled spin-1 and spin-1/2 resonances which mix linearly with the elementary SM particles, we are able to compute the effective potential of the theory by means of some generalized Weinberg sum rules. The properties of the Higgs boson, DM, top quark and the above resonances are thus calculable and tightly connected. We perform a wide phenomenological analysis, considering both collider physics at the LHC and astrophysical observables. We find that these models are tightly constrained by present experimental data, which are able to completely exclude the most natural setup with $f\simeq 800$ GeV. Upon increasing the value of $f$, an allowed region appears. In particular for $f\simeq 1.1$ TeV we find a concrete realization that predicts $m_{DM}\simeq 200$ GeV for the DM mass. This DM candidate lies close to the present sensitivity of direct detection experiments and will be ruled out - or discovered - in the near future.

The analysis of the gamma-ray photons collected by the Fermi Large Area Telescope reveals, after removal of astrophysical background, the existence of an excess towards the Galactic center. This excess peaks around few GeV, and its origin is compatible with the gamma-ray flux originating from Dark Matter annihilation. In this work we take a closer look on this interpretation; we investigate which kind of Dark Matter, and which type of interactions with the Standard Model fields are able to reproduce the observed signal. The structure of the paper is twofold. In the first part, we follow an effective field theory approach considering both fermionic and scalar Dark Matter. The computation of the relic density, the constraint imposed from the null result of direct searches, and the reliability of the effective field theory description allow us to single out only two viable dim-6 operators in the case of fermionic Dark Matter. In the second part, we analyze some concrete models. In particular, we find that the scalar Higgs portal can provide a simple, concrete and realistic scenario able to explain the GeV excess under scrutiny.

The quest for Dark Matter signals in the gamma-ray sky is one of the most intriguing and exciting challenges in astrophysics. In this paper we perform the analysis of the energy spectrum of the \textit{Fermi bubbles} at different latitudes, making use of the gamma-ray data collected by the Fermi Large Area Telescope. By exploring various setups for the full-sky analysis we achieve stable results in all the analyzed latitudes. At high latitude, $|b|=20^{\circ}-50^{\circ}$, the \textit{Fermi bubbles} energy spectrum can be reproduced by gamma-ray photons generated by inverse Compton scattering processes, assuming the existence of a population of high-energy electrons. At low latitude, $|b|=10^{\circ}-20^{\circ}$, the presence of a bump at $E_{\gamma}\sim 1-4$ GeV, reveals the existence of an extra component compatible with Dark Matter annihilation. Our best-fit candidate corresponds to annihilation into $b\bar{b}$ with mass $M_{\rm DM}= 61.8^{+6.9}_{-4.9}$ GeV and cross section $<\sigma v> = 3.30^{+0.69}_{-0.49}\times 10^{-26}$ cm$^{3}$s$^{-1}$. In addition, using the energy spectrum of the \textit{Fermi bubbles}, we derive new conservative but stringent upper limits on the Dark Matter annihilation cross section.

We analyze the validity of the theorems concerning the cancellation of the infrared and collinar divergences in the case of dark matter freeze-out in the early universe. In particular, we compute the electroweak logarithmic corrections of infrared origin to the annihilation cross section of a dark matter particle being the neutral component of a SU(2)_L multiplet. The inclusion of processes with final state W can modify significantly the cross sections computed with only virtual W exchange. Our results show that the inclusion of infrared logs is necessary for a precise computation of the dark matter relic abundance.

We show that the dark matter (DM) could be a light composite scalar $\eta$, emerging from a TeV-scale strongly-coupled sector as a pseudo Nambu-Goldstone boson (pNGB). Such state arises naturally in scenarios where the Higgs is also a composite pNGB, as in $O(6)/O(5)$ models, which are particularly predictive, since the low-energy interactions of $\eta$ are determined by symmetry considerations. We identify the region of parameters where $\eta$ has the required DM relic density, satisfying at the same time the constraints from Higgs searches at the LHC, as well as DM direct searches. Compositeness, in addition to justify the lightness of the scalars, can enhance the DM scattering rates and lead to an excellent discovery prospect for the near future. For a Higgs mass $m_h\simeq 125$ GeV and a pNGB characteristic scale $f \lesssim 1$ TeV, we find that the DM mass is either $m_\eta \simeq 50-70$ GeV, with DM annihilations driven by the Higgs resonance, or in the range 100-500 GeV, where the DM derivative interaction with the Higgs becomes dominant. In the former case the invisible Higgs decay to two DM particles could weaken the LHC Higgs signal.

Recent analyses have shown that the inclusion of electroweak corrections can alter significantly the energy spectra of Standard Model particles originated from dark matter annihilations. We investigate the important situation where the radiation of electroweak gauge bosons has a substantial influence: a Majorana dark matter particle annihilating into two light fermions. This process is in p-wave and hence suppressed by the small value of the relative velocity of the annihilating particles. The inclusion of electroweak radiation eludes this suppression and opens up a potentially sizeable s-wave contribution to the annihilation cross section. We study this effect in detail and explore its impact on the fluxes of stable particles resulting from the dark matter annihilations, which are relevant for dark matter indirect searches. We also discuss the effective field theory approach, pointing out that the opening of the s-wave is missed at the level of dimension-six operators and only encoded by higher orders.

The computation of the energy spectra of Standard Model particles originated from the annihilation/decay of dark matter particles is of primary importance in indirect searches of dark matter. We compute how the inclusion of electroweak corrections significantly alter such spectra when the mass M of dark matter particles is larger than the electroweak scale: soft electroweak gauge bosons are copiously radiated opening new channels in the final states which otherwise would be forbidden if such corrections are neglected. All stable particles are therefore present in the final spectrum, independently of the primary channel of dark matter annihilation/decay. Such corrections are model independent.