This article reviews the concept of Lorentz invariant relative velocity that is often misunderstood or unknown in high energy physics literature. The properties of the relative velocity allow to formulate the invariant flux and cross section without recurring to non--physical velocities or any assumption about the reference frame. Applications such as the luminosity of a collider, the use as kinematic variable, and the statistical theory of collisions in a relativistic classical gas are reviewed. It is emphasized how the hyperbolic properties of the velocity space explain the peculiarities of relativistic scattering.

We find an exact formula for the thermally averaged cross section times the relative velocity $\langle \sigma v_{\text{rel}} \rangle$ with relativistic Maxwell-Boltzmann statistics. The formula is valid in the effective field theory approach when the masses of the annihilation products can be neglected compared with the dark matter mass and cut-off scale. The expansion at $x=m/T\gg 1$ directly gives the nonrelativistic limit of $\langle \sigma v_{\text{rel}}\rangle$ which is usually used to compute the relic abundance for heavy particles that decouple when they are nonrelativistic. We compare this expansion with the one obtained by expanding the total cross section $\sigma(s)$ in powers of the nonrelativistic relative velocity $v_r$. We show the correct invariant procedure that gives the nonrelativistic average $\langle \sigma_{nr} v_r \rangle_{nr}$ coinciding with the large $x$ expansion of $\langle \sigma v_{\text{rel}}\rangle$ in the comoving frame. We explicitly formulate flux, cross section, thermal average, collision integral of the Boltzmann equation in an invariant way using the true relativistic relative $v_\text{rel}$, showing the uselessness of the M\o{}ller velocity and further elucidating the conceptual and numerical inconsistencies related with its use.

We show that the standard theory of thermal production and chemical decoupling of WIMPs is incomplete. The hypothesis that WIMPs are produced and decouple from a thermal bath implies that the rate equation the bath particles interacting with the WIMPs is an algebraic equation that constraints the actual WIMPs abundance to have a precise analytical form down to the temperature $x_\ast=m_\chi /T_\ast$. The point $x_\ast$, which coincides with the stationary point of the equation for the quantity $\Delta= Y-Y_0$, is where the maximum departure of the WIMPs abundance $Y$ from the thermal value $Y_0$ is reached. For each mass $m_\chi$ and total annihilation cross section $\langle \sigma_\text{ann}v_\text{r}\rangle$, the temperature $x_\ast$ and the actual WIMPs abundance $Y(x_\ast)$ are exactly known. This value provides the true initial condition for the usual differential equation that have to be integrated in the interval $x\ge x_\ast$. The matching of the two abundances at $x_\ast$ is continuous and differentiable. The dependence of the present relic abundance on the abundance at an intermediate temperature is an exact result. The exact theory suggests a new analytical approximation that furnishes the relic abundance accurate at the level of $1\%-2\%$ in the case of $S$-wave and $P$-wave scattering cross sections. We conclude the paper studying the evolution of the WIMPs chemical potential and the entropy production using methods of non equilibrium thermodynamics.

We find the probability density function $\mathcal{P}(V_{\texttt{r}})$ of the relativistic relative velocity for two colliding particles in a non-degenerate relativistic gas. The distribution reduces to Maxwell distribution for the relative velocity in the non-relativistic limit. We find an exact formula for the mean value $\langle V_{\texttt{r}}\rangle$. The mean velocity tends to the Maxwell's value in the non-relativistic limit and to the velocity of light in the ultra-relativistic limit. At a given temperature $T$, when at least for one of the two particles the ratio of the rest energy over the thermal energy $m c^2/k_B T$ is smaller than 40 the Maxwell distribution is inadequate.

In this paper we clarify the relation between the invariant relativistic relative velocity $V_\text{r}$, the M\o{}ller velocity $\bar{v}$, and the non-relativistic relative velocity $v_r$. Adopting $V_{\text{r}}$ as the true physical relative velocity for pair-collisions in a non-degenerate relativistic gas, we show that in the frame co-moving with the gas (i) the thermally averaged cross section times relative velocity $<\sigma v_\text{rel}>$ that appears in the density evolution equation for thermal relics is reformulated only in terms of $V_\text{r}$ and $\mathcal{P}(V_{\text{r}})$ in a manifestly Lorentz invariant form; (ii) the frame-dependent issues of the standard formulation in terms of the M\o{}ller velocity, as well as "superluminal" relative velocities, are not present in this formulation. Furthermore, considering the annihilation of dark matter into a particle-antiparticle pair $f\bar{f}$, in the cases $m_f=0$, $m_f=m$ and $m_f \gg m$, we find that the coefficients of the low velocity expansion of $<\sigma V_{\text{r}}>$ admit an exact analytical representation in terms of the Meijer $G$ functions that can be reduced to combinations of modified Bessel functions of the second kind.

We show how to include in the existing calculations for nuclei other than $^{129}$Xe, $^{131}$Xe, the corrections to the isovector coupling arising in chiral effective field theory recently found in Ref. \cite{Menendez1}. The dominant, momentum independent, 2-body currents effect can be taken into account by formally redefining the static spin matrix elements $< \mathbf{S}_{p,n} >$. By further using the normalized form factor at $q\neq 0$ built with the 1-body level structure functions, we show that the WIMP-nucleus cross section and the upper limits on the WIMP-nucleon cross sections coincide with the ones derived using the exact functions at the 2-body level. We explicitly show it in the case of XENON100 limits on the WIMP-neutron cross section and we recalculate the limits on the WIMP-proton spin dependent cross section set by COUPP. We also give practical formulas to obtain $< \mathbf{S}_{p,n} >$ given the structure functions in the various formalisms/notations existing in literature. We argue that the standard treatment of the spin-dependent cross section in terms of three independent isospin functions, $S_{00}(q)$, $S_{11}(q)$, $S_{01}(q)$, is redundant in the sense that the interference function $S_{01}(q)$ is the double product $|S_{01}(q)|=2\sqrt{S_{00}(q)}\sqrt{S_{11}(q)}$ even when including the new effective field theory corrections.

We discuss on the possibility that colliding dark matter particles in the form of neutralinos may be gravitationally boosted near the super-massive black hole at the galactic center so that they can have enough collision energy to annihilate into a stau pair. Since in some phenomenologically favored supersymmetric models the mass splitting between the neutralino and the lightest stau, one of the two scalar superpartners of the tau lepton, is a few GeVs, this channel may be allowed. In addition, staus can only decay into a tau lepton and another neutralino. We calculate the gamma-ray spectrum and flux generated by the tau pair discussing the observability of the obtained features.

In this paper, we compare dwarf galaxies and galaxy clusters in order to elucidate which object class is the best target for gamma-ray DM searches with imaging atmospheric Cherenkov telescopes (IACTs). We have built a mixed dwarfs+clusters sample containing some of the most promising nearby dwarf galaxies (Draco, Ursa Minor, Wilman 1 and Segue 1) and local galaxy clusters (Perseus, Coma, Ophiuchus, Virgo, Fornax, NGC5813 and NGC5846), and then compute their DM annihilation flux profiles by making use of the latest modeling of their DM density profiles. We also include in our calculations the effect of DM substructure. Willman 1 appears as the best candidate in the sample. However, its mass modeling is still rather uncertain, so probably other candidates with less uncertainties and quite similar fluxes, namely Ursa Minor and Segue 1, might be better options. As for galaxy clusters, Virgo represents the one with the highest flux. However, its large spatial extension can be a serious handicap for IACT observations and posterior data analysis. Yet, other local galaxy cluster candidates with more moderate emission regions, such as Perseus, may represent good alternatives. After comparing dwarfs and clusters, we found that the former exhibit annihilation flux profiles that, at the center, are roughly one order of magnitude higher than those of clusters, although galaxy clusters can yield similar, or even higher, integrated fluxes for the whole object once substructure is taken into account. Even when any of these objects are strictly point-like according to the properties of their annihilation signals, we conclude that dwarf galaxies are best suited for observational strategies based on the search of point-like sources, while galaxy clusters represent best targets for analyses that can deal with rather extended emissions. Finally, we study the detection prospects for IACTs [ABRIDGED]

We revise the spin-dependent neutralino-nucleus elastic scattering comparing the formalisms and approximations found in literature for the momentum transfer dependent structure functions. We argue that one of the normalized structure functions of Divari, Kosmas, Vergados and Skouras is all that one needs to correctly take into account the detailed nuclear physics information provided by shell-model calculations. The factorization of the particle physics degrees of freedom from the nuclear physics momentum dependent structure functions implied by this formalism allows for a better understanding of the so-called model independent method for setting upper limits. We further discuss the possibility of experiments with spin-dependent sensitivity like COUPP to test or set limits on the proton spin-dependent cross section in the framework of the stau co-annihilation region of the constrained minimal supersymmetric standard model. For this model with $A_0 =0$, we provide a fitting formula by which it is possible to convert an upper limit on the spin-independent cross section as a function of the neutralino mass directly into an exclusion plot in the ($m_{1/2}$, $\tan\beta$) plane.

We propose a scheme that allows to analytically determine the three elementary cross sections and connect the solutions to the relative sign between the proton and the neutron spin scattering amplitudes once the measurements of total event rate from three appropriate targets become available. In this way it is thus possible to extract the maximum information on the supersymmetric parameter space obtainable from direct detection experiments, in the case that the dark matter particle is the lightest neutralino. Our scheme is based on suitably normalized form of the isospin momentum dependent structure functions entering in the spin-dependent elastic neutralino-nucleus cross section. We compare these functions with the commonly used ones and discuss their advantages: in particular, these allow in the spin-dependent cross section to factorize the particle physics degrees of freedom from the momentum transfer dependent nuclear structure functions as it happens in the spin-independent cross section with the nuclear form factor.

We present a detailed study of the effect of internal bremsstrahlung photons in the context of the minimal supersymmetric standard models and their impact on gamma-ray dark matter annihilation searches. We find that although this effect has to be included for the correct evaluation of fluxes of high energy photons from neutralino annihilation, its contribution is relevant only in models and at energies where the lines contribution is dominant over the secondary photons. Therefore, we find that the most optimistic supersymmetric scenarios for dark matter detection do not change significantly when including the internal bremsstrahlung. As an example, we review the gamma-ray dark matter detection prospects of the Draco dwarf spheroidal galaxy for the MAGIC stereoscopic system and the CTA project. Though the flux of high energy photons is enhanced by an order of magnitude in some regions of the parameter space, the expected fluxes are still much below the sensitivity of the instruments.