Assuming that Dark Matter (DM) is made of fermions in the sub-GeV mass range with interactions dominated by electromagnetic moments of higher order, such as the electric and magnetic dipoles or the anapole moment, we show that direct detection experiments searching for atomic ionisation events in xenon targets can shed light on whether DM is a Dirac or Majorana particle. Specifically, we find that between about 45 (120) and 610 (1700) signal events are required to reject Majorana DM in favour of Dirac DM with a statistical significance corresponding to 3 (5) standard deviations. The exact number of DM signal events corresponding to a given significance depends on the relative size of the anapole, magnetic dipole and electric dipole contributions to the expected rate of DM-induced atomic ionisations under the Dirac hypothesis. Our conclusions are based on Monte Carlo simulations and the likelihood ratio test. While the use of asymptotic formulae for the latter is standard in many applications, here it requires a non-trivial extension to the case where one of the hypotheses lies on the boundary of the parameter space. Our results constitute a solid proof of concept about the possibility of using direct detection experiments to reject the Majorana DM hypothesis when the DM interactions are dominated by higher-order electromagnetic moments.

Discrete symmetries played a central role in elucidating the structure of the weak interactions, and they will probably be equally crucial regarding the interactions of the dark matter (DM) particle -- whose nature remains elusive. In this work we show that signals in future direct detection experiments can be used to test, in a model-independent way, for P- and CP-violation in DM-nucleus interactions. The analysis is performed within the most general effective theory for non-relativistic spin-0 DM-nucleus interactions mediated by the exchange of a heavy particle. Assuming an idealised xenon detector, we calculate the expected number of DM signal events required to reject P and CP invariant DM-nucleus interactions. For a DM mass of 30 GeV (or higher), this number lies between about 10 and 300 DM signal events, depending on how P and CP invariance are modeled. Future direct detection experiments, therefore, have the potential to reveal P- and CP-violation in DM interactions, making a decisive step toward the identification of the DM particle.

Despite strong evidence for the existence of large amounts of dark matter (DM) in our Universe, there is no direct indication of its presence in our own solar system. All estimates of the local DM density, crucial for all direct DM searches, rely on extrapolating estimates on much larger scales. We demonstrate for the first time the possibility of measuring the local DM density with a direct detection experiment. It relies on the assumption that incoming DM particles frequently scatter on terrestrial nuclei prior to detection, inducing an additional time-dependence of the signal. We show that for sub-GeV DM, with a spin-independent DM-proton cross section $\sim10^{-32}\,\text{cm}^2$, future direct detection experiments should be able to reconstruct the local DM density with 20% uncertainty.

We compute the projected sensitivity to dark matter (DM) particles in the sub-GeV mass range of future direct detection experiments using germanium and silicon semiconductor targets. We perform this calculation within the dark photon model for DM-electron interactions using the likelihood ratio as a test statistic, Monte Carlo simulations, and background models that we extract from recent experimental data. We present our results in terms of DM-electron scattering cross section values required to reject the background only hypothesis in favour of the background plus DM signal hypothesis with a statistical significance, $\mathcal{Z}$, corresponding to 3 or 5 standard deviations. We also test the stability of our conclusions under changes in the astrophysical parameters governing the local space and velocity distribution of DM in the Milky Way. In the best-case scenario, when a high-voltage germanium detector with an exposure of $50$ kg-year and a CCD silicon detector with an exposure of $1$ kg-year and a dark current rate of $1\times10^{-7}$ counts/pixel/day have simultaneously reported a DM signal, we find that the smallest cross section value compatible with $\mathcal{Z}=3$ ($\mathcal{Z}=5$) is about $8\times10^{-42}$ cm$^2$ ($1\times10^{-41}$ cm$^2$) for contact interactions, and $4\times10^{-41}$ cm$^2$ ($7\times10^{-41}$ cm$^2$) for long-range interactions. Our sensitivity study extends and refine previous works in terms of background models, statistical methods, and treatment of the underlying astrophysical uncertainties.

In the leading paradigm of modern cosmology, about 80% of our Universe's matter content is in the form of hypothetical, as yet undetected particles. These do not emit or absorb radiation at any observable wavelengths, and therefore constitute the so-called Dark Matter (DM) component of the Universe. Detecting the particles forming the Milky Way DM component is one of the main challenges for astroparticle physics and basic science in general. One promising way to achieve this goal is to search for rare DM-electron interactions in low-background deep underground detectors. Key to the interpretation of this search is the response of detectors' materials to elementary DM-electron interactions defined in terms of electron wave functions' overlap integrals. In this work, we compute the response of atomic argon and xenon targets used in operating DM search experiments to general, so far unexplored DM-electron interactions. We find that the rate at which atoms can be ionized via DM-electron scattering can in general be expressed in terms of four independent atomic responses, three of which we identify here for the first time. We find our new atomic responses to be numerically important in a variety of cases, which we identify and investigate thoroughly using effective theory methods. We then use our atomic responses to set 90% confidence level (C.L.) exclusion limits on the strength of a wide range of DM-electron interactions from the null result of DM search experiments using argon and xenon targets.

We investigate the non-relativistic reduction of simplified models for spin 1 dark matter (DM) with the aim of identifying features in the phenomenology of DM-quark interactions which are specific to vector DM. In the case of DM-quark interactions mediated by a spin 1 particle, we find two DM-nucleon interaction operators arising from the non-relativistic reduction of simplified models for spin 1 DM that are specific to spin 1 DM, and which were not considered in previous studies. They are quadratic in the momentum transfer, linear in a symmetric combination of polarisation vectors for the DM particle, and arise from simplified models which do not generate momentum transfer independent operators as leading interactions in the non-relativistic expansion of DM-nucleon scattering amplitudes. Within these simplified models, the new operators cannot be neglected when computing DM signals at direct detection experiments. For example, we find that nuclear recoil energy spectra computed by including or neglecting the new operators can differ by up to one order of magnitude for nuclear recoil energies larger than about 20 keV and DM masses below $50$ GeV. Furthermore, the shape of the expected nuclear recoil spectra depends significantly on whether the new operators are taken into account or not. Finally, neglecting the contribution to DM direct detection signals from the new operators leads to inaccurate conclusions when assessing the compatibility of a future direct detection signal with CMB constraints on the DM relic density, especially when the number of signal events is small, e.g. $\mathcal{O}(1)$.

We study the capture and subsequent annihilation of inelastic dark matter (DM) in the Sun, placing constraints on the DM-nucleon scattering cross section from the null result of the IceCube neutrino telescope. We then compare such constraints with exclusion limits on the same cross section that we derive from XENON1T, PICO and CRESST results. We calculate the cross section for inelastic DM-nucleon scattering within an extension of the effective theory of DM-nucleon interactions which applies to the case of inelastic DM models characterised by a mass splitting between the incoming and outgoing DM particle. We find that for values of the mass splitting parameter larger than about 200 keV, neutrino telescopes place limits on the DM-nucleon scattering cross section which are stronger than the ones from current DM direct detection experiments. The exact mass splitting value for which this occurs depends on whether DM thermalises in the Sun or not. This result applies to all DM-nucleon interactions that generate DM-nucleus scattering cross sections which are independent of the nuclear spin, including the "canonical" spin-independent interaction. We explicitly perform our calculations for a DM candidate with mass of 1 TeV, but our conclusions qualitatively also apply to different masses. Furthermore, we find that exclusion limits from IceCube on the coupling constants of this family of spin-independent interactions are more stringent than the ones from a (hypothetical) reanalysis of XENON1T data based on an extended signal region in nuclear recoil energy. Our results should be taken into account in global analyses of inelastic DM models.

The discovery of dark matter (DM) at XENONnT or LZ would place constraints on DM particle mass and coupling constants. It is interesting to ask when these constraints can be compatible with the DM thermal production mechanism. We address this question within the most general set of renormalisable models that preserve Lorentz and gauge symmetry, and that extend the Standard Model by one DM candidate of mass $m_{\rm DM}$ and one particle of mass $M_{med}$ mediating DM-quark interactions. Our analysis divides into two parts. First, we postulate that XENONnT/LZ has detected $\mu_S\sim\mathcal{O}(100)$ signal events, and use this input to calculate the DM relic density, $\Omega_{DM} h^2$. Then, we identify the regions in the $M_{med} - \Omega_{DM} h^2$ plane which are compatible with the observed signal and with current CMB data. We find that for most of the models considered here, $\mathcal{O}(100)$ signal events at XENONnT/LZ and the DM thermal production are only compatible for resonant DM annihilations, i.e. for $M_{med}\simeq2 m_{DM}$. In this case, XENONnT/LZ would be able to simultaneously measure $m_{DM}$ and $M_{med}$. We also discuss the dependence of our results on $m_{DM}$, $\mu_S$ and the DM spin, and provide analytic expressions for annihilation cross-sections and mediator decay widths for all models considered in this study.

We develop a method to forecast the outcome of the LHC Run 3 based on the hypothetical detection of $\mathcal{O}(100)$ signal events at XENONnT. Our method relies on a systematic classification of renormalisable single-mediator models for dark matter-quark interactions, and is valid for dark matter candidates of spin less than or equal to one. Applying our method to simulated data, we find that at the end of the LHC Run 3 only two mutually exclusive scenarios would be compatible with the detection of $\mathcal{O}(100)$ signal events at XENONnT. In a first scenario, the energy distribution of the signal events is featureless, as for canonical spin-independent interactions. In this case, if a mono-jet signal is detected at the LHC, dark matter must have spin 1/2 and interact with nucleons through a unique velocity-dependent operator. If a mono-jet signal is not detected, dark matter interacts with nucleons through canonical spin-independent interactions. In a second scenario, the spectral distribution of the signal events exhibits a bump at non zero recoil energies. In this second case, a mono-jet signal can be detected at the LHC Run 3, dark matter must have spin 1/2 and interact with nucleons through a unique momentum-dependent operator. We therefore conclude that the observation of $\mathcal{O}(100)$ signal events at XENONnT combined with the detection, or the lack of detection, of a mono-jet signal at the LHC Run 3 would significantly narrow the range of possible dark matter-nucleon interactions. As we argued above, it can also provide key information on the dark matter particle spin.

The theoretical interpretation of dark matter direct detection experiments is hindered by uncertainties of the microphysics governing the dark matter-nucleon interaction, and of the dark matter density and velocity distribution inside the Solar System. These uncertainties are especially relevant when confronting a detection claim to the null results from other experiments, since seemingly conflicting experimental results may be reconciled when relaxing the assumptions about the form of the interaction and/or the velocity distribution. We present in this paper a halo-independent method to calculate the maximum number of events in a direct detection experiment given a set of null search results, allowing for the first time the scattering to be mediated by an arbitrary combination of various interactions (concretely we consider up to 64). We illustrate this method to examine the compatibility of the dark matter interpretation of the three events detected by the silicon detectors in the CDMS-II experiment with the null results from XENON1T and PICO-60.

Dwarf spheroidal galaxies are among the most promising targets for indirect dark matter (DM) searches in $\gamma$-rays. The $\gamma$-ray flux from DM annihilation in a dwarf spheroidal galaxy is proportional to the $J$-factor of the source. The $J$-factor of a dwarf spheroidal galaxy is the line-of-sight integral of the DM mass density squared times $\langle \sigma_{\rm ann} v_{\rm rel} \rangle/(\sigma_{\rm ann} v_{\rm rel})_0$, where $\sigma_{\rm ann} v_{\rm rel}$ is the DM annihilation cross-section times relative velocity $v_{\rm rel}=|{\bf v}_{\rm rel}|$, angle brackets denote average over ${\bf v}_{\rm rel}$, and $(\sigma_{\rm ann} v_{\rm rel})_0$ is the $v_{\rm rel}$-independent part of $\sigma_{\rm ann} v_{\rm rel}$. If $\sigma_{\rm ann} v_{\rm rel}$ is constant in $v_{\rm rel}$, $J$-factors only depend on the DM space distribution in the source. However, if $\sigma_{\rm ann} v_{\rm rel}$ varies with $v_{\rm rel}$, as in the presence of DM self-interactions, $J$-factors also depend on the DM velocity distribution, and on the strength and range of the DM self-interaction. Models for self-interacting DM are increasingly important in the study of the small scale clustering of DM, and are compatible with current cosmological observations. Here we derive the $J$-factor of 20 dwarf spheroidal galaxies from stellar kinematic data under the assumption of Yukawa DM self-interactions. $J$-factors are derived through a profile Likelihood approach, assuming either NFW or cored DM profiles. We also compare our results with $J$-factors derived assuming the same velocity for all DM particles in the target galaxy. We find that this common approximation overestimates the $J$-factors by up to one order of magnitude. $J$-factors for a sample of DM particle masses, self-interaction coupling constants and density profiles are provided electronically, ready to be used in other projects.

If the dark matter particle has spin 0, only two types of WIMP-nucleon interaction can arise from the non-relativistic reduction of renormalisable single-mediator models for dark matter-quark interactions. Based on this crucial observation, we show that about 100 signal events at next generation directional detection experiments can be enough to enable a $2\sigma$ rejection of the spin 0 dark matter hypothesis in favour of alternative hypotheses where the dark matter particle has spin 1/2 or 1. In this context directional sensitivity is crucial, since anisotropy patterns in the sphere of nuclear recoil directions depend on the spin of the dark matter particle. For comparison, about 100 signal events are expected in a CF$_4$ detector operating at a pressure of 30 torr with an exposure of approximately 26,000 cubic-meter-detector days for WIMPs of 100 GeV mass and a WIMP-Fluorine scattering cross-section of 0.25 pb. Comparable exposures are within reach of an array of cubic meter time projection chamber detectors.

We study the process of dark matter particles scattering off $^{3,4}$He with nuclear wave functions computed using an ab initio many-body framework. We employ realistic nuclear interactions from chiral effective field theory at next-to-next-to-leading order (NNLO) and develop an ab initio scheme to compute a general set of different nuclear response functions. In particular, we then perform an accompanying uncertainty quantification on these quantities and study error propagation to physical observables. We find a rich structure of allowed nuclear responses with significant uncertainties for certain spin-dependent interactions. The approach and results that are presented in this Paper establish a new framework for nuclear structure calculations and uncertainty quantification in the context of direct and (certain) indirect searches for dark matter.

Direct detection experiments search for the interactions of Dark Matter (DM) particles with nuclei in terrestrial detectors. But if these interactions are sufficiently strong, DM particles may scatter in the Earth, affecting their distribution in the lab. We present a new analytic calculation of this `Earth-scattering' effect in the regime where DM particles scatter at most once before reaching the detector. We perform the calculation self-consistently, taking into account not only those particles which are scattered away from the detector, but also those particles which are deflected towards the detector. Taking into account a realistic model of the Earth and allowing for a range of DM-nucleon interactions, we present the EarthShadow code, which we make publicly available, for calculating the DM velocity distribution after Earth-scattering. Focusing on low-mass DM, we find that Earth-scattering reduces the direct detection rate at certain detector locations while increasing the rate in others. The Earth's rotation induces a daily modulation in the rate, which we find to be highly sensitive to the detector latitude and to the form of the DM-nucleon interaction. These distinctive signatures would allow us to unambiguously detect DM and perhaps even identify its interactions in regions of the parameter space within the reach of current and future experiments.

In the effective theory of isoscalar and isovector dark matter-nucleon interactions mediated by a heavy spin-1 or spin-0 particle, 8 isotope-dependent nuclear response functions can be generated in the dark matter scattering by nuclei. We compute the 8 nuclear response functions for the 16 most abundant elements in the Sun, i.e. H, $^{3}$He, $^{4}$He, $^{12}$C, $^{14}$N, $^{16}$O, $^{20}$Ne, $^{23}$Na, $^{24}$Mg, $^{27}$Al, $^{28}$Si, $^{32}$S, $^{40}$Ar, $^{40}$Ca, $^{56}$Fe, and $^{58}$Ni, through numerical shell model calculations. We use our response functions to compute the rate of dark matter capture by the Sun for all isoscalar and isovector dark matter-nucleon effective interactions, including several operators previously considered for dark matter direct detection only. We study in detail the dependence of the capture rate on specific dark matter-nucleon interaction operators, and on the different elements in the Sun. We find that a so far neglected momentum dependent dark matter coupling to the nuclear vector charge gives a larger contribution to the capture rate than the constant spin-dependent interaction commonly included in dark matter searches at neutrino telescopes. Our investigation lays the foundations for model independent analyses of dark matter induced neutrino signals from the Sun. The nuclear response functions obtained in this study are listed in analytic form in an appendix, ready to be used in other projects.

We constrain the effective theory of one-body dark matter-nucleon interactions using neutrino telescope observations. We derive exclusion limits on the 28 coupling constants of the theory, exploring interaction operators previously considered in dark matter direct detection only, and using new nuclear response functions recently derived through nuclear structure calculations. We determine for what interactions neutrino telescopes are superior to current direct detection experiments, and show that Hydrogen is not the most important element in the exclusion limit calculation for the majority of the spin-dependent operators.

I calculate the rate of WIMP capture and annihilation in the Earth in the non-relativistic effective theory of dark matter-nucleon interactions. Neglecting operator interference, I consider all Galilean invariant interaction operators that can arise from the exchange of a heavy particle of spin less than or equal to one when WIMPs have spin 0, 1/2 or 1. I compute position and shape of the expected resonances in the mass - capture rate plane and show that Iron is not the most important element in the capture process for many currently ignored interaction operators. I compare these predictions with the recent results of an Earth WIMP analysis of IceCube in the 86-string configuration and set limits on all isoscalar and isovector coupling constants of the effective theory of dark matter-nucleon interactions. For certain interaction operators and for a dark matter particle mass of about 50 GeV, I find that these limits are stronger than those I have previously derived in an analysis of the solar WIMP search performed at IceCube in the 79-string configuration.

We study the capture of WIMP dark matter by the Sun in the non-relativistic effective theory of dark matter self-interactions. The aim is to assess the impact of self-interactions on the expected neutrino flux from the annihilation of WIMPs trapped in the Sun in a model independent manner. We consider all non-relativistic Galilean invariant self-interaction operators that can arise from the exchange of a heavy particle of spin less than or equal to 1 for WIMPs of spin equal to 0, 1/2 and 1. We show that for interaction operators depending at most linearly on the momentum transfer, the WIMP-induced neutrino flux can be enhanced by several orders of magnitude compared to the same flux in absence of self-interactions. This is true even for standard values of the thermally averaged annihilation cross-section. This conclusion impacts the analysis of present and future observations performed at neutrino telescopes.

We study the properties and direct detection prospects of an as of yet neglected population of dark matter (DM) particles moving in orbits gravitationally bound to the Earth. This DM population is expected to form via scattering by nuclei in the Earth's interior. We compute fluxes and nuclear recoil energy spectra expected at direct detection experiments for the new DM population considering detectors with and without directional sensitivity, and different types of target materials and DM-nucleon interactions. DM particles bound to the Earth manifest as a prominent rise in the low-energy part of the observed nuclear recoil energy spectrum. Ultra-low threshold energies of about 1 eV are needed to resolve this effect. Its shape is independent of the DM-nucleus scattering cross-section normalisation.

We examine the dark matter interpretation of the modulation signal reported by the DAMA experiment from the perspective of effective field theories displaying Galilean invariance. We consider the most general effective coupling leading to the elastic scattering of a dark matter particle with spin 0 or 1/2 off a nucleon, and we analyze the compatibility of the DAMA signal with the null results from other direct detection experiments, as well as with the non-observation of a high energy neutrino flux in the direction of the Sun from dark matter annihilation. To this end, we develop a novel semi-analytical approach for comparing experimental results in the high-dimensional parameter space of the non-relativistic effective theory. Assuming the standard halo model, we find a strong tension between the dark matter interpretation of the DAMA modulation signal and the null result experiments. We also list possible ways-out of this conclusion.

We demonstrate that dark matter particles gravitationally bound to the Earth can induce a characteristic nuclear recoil signal at low energies in direct detection experiments. The new spectral feature we predict can provide the ultimate smoking gun for dark matter discovery for experiments with positive signal but unclear background. The new feature is universal, in that the ratio of bound over halo dark matter event rates at detectors is independent of the dark matter-nucleon cross section.

I compare the non-relativistic effective theory of one-body dark matter-nucleon interactions to current dark matter direct detection experiments and neutrino telescope observations, presenting exclusion limits on the coupling constants of the theory. In the analysis of direct detection experiments, I focus on the interference of different dark matter-nucleon interaction operators and on predictions observable at directional detectors. Interpreting neutrino telescope observations, I use new nuclear response functions recently derived through nuclear structure calculations and show that hydrogen is not the most important element in the exclusion limit calculation for the majority of the spin-dependent dark matter-nucleon interaction operators

We revisit astrophysical constraints on models with an ultra-light gravitino, in particular extending the analysis to more general models and the case of R-parity breaking. These constraints allow to restrict the value of the gravitino mass depending on the masses of scalar moduli. We perform both a frequentist and bayesian analysis and we find comparable results, even if the bayesian analysis is in part affected by volume effects. Only a small window of gravitino masses can be excluded by the SN constraints in a model-independent way, while limits obtained with a definite assumption on the scalar masses result more stringent.

The effective theory of isoscalar dark matter-nucleon interactions mediated by heavy spin-one or spin-zero particles depends on 10 coupling constants besides the dark matter particle mass. Here we compare this 11-dimensional effective theory to current observations in a comprehensive statistical analysis of several direct detection experiments, including the recent LUX, SuperCDMS and CDMSlite results. From a multidimensional scan with about 3 million likelihood evaluations, we extract the marginalized posterior probability density functions (a Bayesian approach) and the profile likelihoods (a frequentist approach), as well as the associated credible regions and confidence levels, for each coupling constant vs dark matter mass and for each pair of coupling constants. We compare the Bayesian and frequentist approach in the light of the currently limited amount of data. We find that current direct detection data contain sufficient information to simultaneously constrain not only the familiar spin-independent and spin-dependent interactions, but also the remaining velocity and momentum dependent couplings predicted by the dark matter-nucleon effective theory. For current experiments associated with a null result, we find strong correlations between some pairs of coupling constants. For experiments that claim a signal (i.e., CoGeNT and DAMA), we find that pairs of coupling constants produce degenerate results.

We perform the first comprehensive analysis of the prospects for direct detection of dark matter with future ton-scale detectors in the general 11-dimensional effective theory of isoscalar dark matter-nucleon interactions mediated by a heavy spin-1 or spin-0 particle. The theory includes 8 momentum and velocity dependent dark matter-nucleon interaction operators, besides the familiar spin-independent and spin-dependent operators. From a variegated sample of 27 benchmark points selected in the parameter space of the theory, we simulate independent sets of synthetic data for ton-scale Germanium and Xenon detectors. From the synthetic data, we then extract the marginal posterior probability density functions and the profile likelihoods of the model parameters. The associated Bayesian credible regions and frequentist confidence intervals allow us to assess the prospects for direct detection of dark matter at the 27 benchmark points. First, we analyze the data assuming the knowledge of the correct dark matter nucleon-interaction type, as it is commonly done for the familiar spin-independent and spin-dependent interactions. Then, we analyze the simulations extracting the dark matter-nucleon interaction type from the data directly, in contrast to standard analyses. This second approach requires an extensive exploration of the full 11-dimensional parameter space of the dark matter-nucleon effective theory. Interestingly, we identify 5 scenarios where the dark matter mass and the dark matter-nucleon interaction type can be reconstructed from the data simultaneously. We stress the importance of extracting the dark matter nucleon-interaction type from the data directly, discussing the main challenges found addressing this complex 11-dimensional problem.

We extend the formalism of dark matter directional detection to arbitrary one-body dark matter-nucleon interactions. The new theoretical framework generalizes the one currently used, which is based on 2 types of dark matter-nucleon interaction only. It includes 14 dark matter-nucleon interaction operators, 8 isotope-dependent nuclear response functions, and the Radon transform of the first 2 moments of the dark matter velocity distribution. We calculate the recoil energy spectra at dark matter directional detectors made of CF$_4$, CS$_2$ and $^{3}$He for the 14 dark matter-nucleon interactions, using nuclear response functions recently obtained through numerical nuclear structure calculations. We highlight the new features of the proposed theoretical framework, and present our results for a spherical dark matter halo and for a stream of dark matter particles. This study lays the foundations for model independent analyses of dark matter directional detection experiments.

We compare the general effective theory of one-body dark matter nucleon interactions to current direct detection experiments in a global multidimensional statistical analysis. We derive exclusion limits on the 28 isoscalar and isovector coupling constants of the theory, and show that current data place interesting constraints on dark matter-nucleon interaction operators usually neglected in this context. We characterize the interference patterns that can arise in dark matter direct detection from pairs of dark matter-nucleon interaction operators, or from isoscalar and isovector components of the same operator. We find that commonly neglected destructive interference effects weaken standard direct detection exclusion limits by up to one order of magnitude in the coupling constants.

The Time Renormalization Group (TRG) is an effective method for accurate calculations of the matter power spectrum at the scale of the first baryonic acoustic oscillations. By using a particular variable transformation in the TRG formalism, we can reduce the 2D integral in the source term of the equations of motion for the power spectrum into a series of 1D integrals. The shape of the integrand allows us to pre-compute only thirteen antiderivatives numerically, which can then be reused when evaluating the outer integral. While this introduces a few challenges to keep numerical noise under control, we find that the computation time for nonlinear corrections to the matter power spectrum decreases by a factor of 50. This opens up the possibility to use of TRG for mass production as in Markov Chain Monte Carlo methods. A Fortran code demonstrating this new algorithm has been made publicly available.

Our peculiar motion with respect to the CMB rest frame represents a preferred direction in the observed CMB sky since it induces an apparent deflection of the observed CMB photons (aberration) and a shift in their frequency (Doppler). Both effects distort the multipoles $a_{\ell m}$'s at all $\ell$'s. Such effects are real as it has been recently measured for the first time by Planck according to what was forecast in some recent papers. However, the common lore when estimating a power spectrum from CMB is to consider that Doppler affects only the $\ell=1$ multipole, neglecting any other corrections. In this work we use simulations of the CMB sky in a boosted frame with a peculiar velocity $\beta = v/c = 1.23 \times 10^{-3}$ in order to assess the impact of such effect on power spectrum estimations in different regions of the sky. We show that the boost induces a north-south asymmetry in the power spectrum which is highly significant and non-negligible, of about (0.58 $\pm$ 0.10)% for half-sky cuts when going up to $\ell$ = 2500. We suggest that these effects are relevant and may account for some of the north-south asymmetries seen in the Planck data, being especially important at small scales. Finally we analyze the particular case of the ACT experiment, which observed only a small fraction of the sky and show that it suffers a bias of about 1% on the power spectrum and of similar size on some cosmological parameters: for example the position of the peaks shifts by 0.5% and the overall amplitude of the spectrum is about 0.4% lower than a full-sky case.

We study the gravitational clustering of spherically symmetric overdensities and the statistics of the resulting dark matter halos in the "symmetron model", in which a new long range force is mediated by a $Z_2$ symmetric scalar field. Depending on the initial radius of the overdensity, we identify two distinct regimes: for small initial radii the symmetron mediated force affects the spherical collapse at all redshifts; for initial radii larger than some critical size this force vanishes before collapse because of the symmetron screening mechanism. In both cases overdensities collapse earlier than in the $\Lambda$CDM and statistically tend to form more massive dark matter halos. Regarding the halo-mass function of these objects, we observe order one departures from standard $\Lambda$CDM predictions. The formalism developed here can be easily applied to other models where fifth-forces participate to the dynamics of the gravitational collapse.

Dark matter N-body simulations suggest that the velocity distribution of dark matter is anisotropic. In this work we employ a mass model for the Milky Way whose parameters are determined from a fit to kinematical data. Then we adopt an ansatz for the dark matter phase space distribution which allows to construct self-consistent halo models which feature a degree of anisotropy as a function of the radius such as suggested by the simulations. The resulting velocity distributions are then used for an analysis of current data from dark matter direct detection experiments. We find that velocity distributions which are radially biased at large galactocentric distances (up to the virial radius) lead to an increased high velocity tail of the local dark matter distribution. This affects the interpretation of data from direct detection experiments, especially for dark matter masses around 10 GeV, since in this region the high velocity tail is sampled. We find that the allowed regions in the dark matter mass-cross section plane as indicated by possible hints for a dark matter signal reported by several experiments as well as conflicting exclusion limits from other experiments shift in a similar way when the halo model is varied. Hence, it is not possible to improve the consistency of the data by referring to anisotropic halo models of the type considered in this work.

We review here the status of different dark matter candidates in the context of supersymmetric models, in particular the neutralino as a realization of the WIMP-mechanism and the gravitino. We give a summary of the recent bounds in direct and indirect detection and also of the LHC searches relevant for the dark matter question. We discuss also the implications of the Higgs discovery for the supersymmetric dark matter models and give the prospects for the future years.

The peculiar motion of an observer with respect to the CMB rest frame induces an apparent deflection of the observed CMB photons, {\it i.e.} aberration, and a shift in their frequency, {\it i.e.} Doppler effect. Both effects distort the temperature multipoles $a_{\ell m}$'s via a mixing matrix at {\it any} $\ell$. The common lore when performing a CMB based cosmological parameter estimation is to consider that Doppler affects only the $l=1$ multipole, and neglect any other corrections. In this paper we check the validity of this assumption in parameter estimation for a Planck-like angular resolution, both for a full-sky ideal experiment and also when sky cuts are included to model CMB foreground contaminations with a sky fraction similar to the Planck satellite. Assuming a simple fiducial cosmological model with five parameters, we simulated CMB temperature maps of the sky and added aberration and Doppler effects to the maps. We then analyzed with a MCMC in a Bayesian framework the maps with and without aberration and Doppler effects in order to assess the ability of reconstructing the parameters of the fiducial model. Although the correction on the power spectrum ${C_\ell}$ is larger than the cosmic variance at $\ell>1000$ and potentially important, we find that the bias on the parameters is negligible for Planck.

The complementarity between dark matter searches at colliders and in underground laboratories is an extraordinarily powerful tool in the quest for dark matter. In the vast majority of the analyses conducted so far these dark matter detection strategies have been profitably combined either to perform global fits in the context of certain particle physics models (e.g. the CMSSM) or to estimate the prospects for a direct dark matter detection given the LHC potential of discovering new physics beyond the Standard Model. In this paper we propose an alternative strategy to combine direct and collider dark matter searches: employing the potential of the upcoming generation of 1-ton direct detection experiments, we show that for certain supersymmetric configurations it is possible to translate the information encoded in an hypothetically discovered direct detection signal into classes of expected signals at the LHC. As an illustrative application of our method, we show that for a 60 GeV neutralino thermally produced via resonant annihilations and identified by a 1-ton direct detection experiment, our approach allows to forecast a clearly identifiable prediction for a LHC final state involving three leptons and missing energy. The strategy presented in this paper to systematically translate a direct detection signal into a prediction for the LHC has the potential to significantly strengthen the complementarity between these two dark matter detection strategies.

The peculiar motion of an observer with respect to the CMB rest frame induces a deflection in the arrival direction of the observed photons (also known as CMB aberration) and a Doppler shift in the measured photon frequencies. As a consequence, aberration and Doppler effects induce non trivial correlations between the harmonic coefficients of the observed CMB temperature maps. In this paper we investigate whether these correlations generate a bias on Non-Gaussianity estimators $f_{NL}$. We perform this analysis simulating a large number of temperature maps with Planck-like resolution (lmax $= 2000$) as different realizations of the same cosmological fiducial model (WMAP7yr). We then add to these maps aberration and Doppler effects employing a modified version of the HEALPix code. We finally evaluate a generalization of the Komatsu, Spergel and Wandelt Non-Gaussianity estimator for all the simulated maps, both when peculiar velocity effects have been considered and when these phenomena have been neglected. Using the value $v/c=1.23 \times 10^{-3}$ for our peculiar velocity, we found that the aberration/Doppler induced Non-Gaussian signal is at most of about half of the cosmic variance $\sigma$ for $f_{NL}$ both in a full-sky and in a cut-sky experimental configuration, for local, equilateral and orthogonal estimators. We conclude therefore that when estimating $f_{NL}$ it is safe to ignore aberration and Doppler effects {\it if} the primordial map is already Gaussian. More work is necessary however to assess whether a map which contains Non-Gaussianity can be significantly distorted by a peculiar velocity.

We present a new determination of the local dark matter phase-space density. This result is obtained implementing, in the limit of isotropic velocity distribution and spherical symmetry, Eddington's inversion formula, which links univocally the dark matter distribution function to the density profile, and applying, within a Bayesian framework, a Markov Chain Monte Carlo algorithm to sample mass models for the Milky Way against a broad and variegated sample of dynamical constraints. We consider three possible choices for the dark matter density profile, namely the Einasto, NFW and Burkert profiles, finding that the velocity dispersion, which characterizes the width in the distribution, tends to be larger for the Burkert case, while the escape velocity depends very weakly on the profile, with the mean value we obtain being in very good agreement with estimates from stellar kinematics. The derived dark matter phase-space densities differ significantly--most dramatically in the high velocity tails--from the model usually taken as a reference in dark matter detection studies, a Maxwell-Boltzmann distribution with velocity dispersion fixed in terms of the local circular velocity and with a sharp truncation at a given value of the escape velocity. We discuss the impact of astrophysical uncertainties on dark matter scattering rates and direct detection exclusion limits, considering a few sample cases and showing that the most sensitive ones are those for light dark matter particles and for particles scattering inelastically. As a general trend, when adopting a self-consistent phase-space density, we find that rates are larger, and hence exclusion limits stronger, than with the standard Maxwell-Boltzmann approximation. Tools for applying our result on the local dark matter phase-space density to other dark matter candidates or experimental setups are provided.

Our peculiar velocity with respect to the CMB rest frame is known to induce a large dipole in the CMB. However, the motion of an observer has also the effect of distorting the anisotropies at all scales, as shown by Challinor and Van Leeuwen (2002), due to aberration and Doppler effects. We propose to measure independently our local motion by using off-diagonal two-point correlation functions for high multipoles. We study the observability of the signal for temperature and polarization anisotropies. We point out that Planck can measure the velocity $\beta$ with an error of about 30% and the direction with an error of about 20 degrees. This method constitutes a cross-check, which can be useful to verify that our CMB dipole is due mainly to our velocity or to disentangle the velocity from other possible intrinsic sources. Although in this paper we focus on our peculiar velocity, a similar effect would result also from other intrinsic vectorial distortion of the CMB which would induce a dipolar lensing. Measuring the off-diagonal correlation terms is therefore a test for a preferred direction on the CMB sky.

Alternative cosmologies, based on extensions of General Relativity, predict modified thermal histories in the Early Universe during the pre Big Bang Nucleosynthesis (BBN) era, epoch which is not directly constrained by cosmological observations. When the expansion rate is enhanced with respect to the standard case, thermal relics typically decouple with larger relic abundances. The correct value of the relic abundance is therefore obtained for larger annihilation cross--sections, as compared to standard cosmology. A direct consequence is that indirect detection rates are enhanced. Extending previous analyses of ours, we derive updated astrophysical bounds on the dark matter annihilation cross sections and use them to constrain alternative cosmologies in the pre--BBN era. We also determine the characteristics of these alternative cosmologies in order to provide the correct value of relic abundance for a thermal relic for the (large) annihilation cross--section required to explain the PAMELA results on the positron fraction, therefore providing a "cosmological boost" solution to the dark matter interpretation of the PAMELA data.

We present a novel study on the problem of constructing mass models for the Milky Way, concentrating on features regarding the dark matter halo component. We have considered a variegated sample of dynamical observables for the Galaxy, including several results which have appeared recently, and studied a 7- or 8-dimensional parameter space - defining the Galaxy model - by implementing a Bayesian approach to the parameter estimation based on a Markov Chain Monte Carlo method. The main result of this analysis is a novel determination of the local dark matter halo density which, assuming spherical symmetry and either an Einasto or an NFW density profile is found to be around 0.39 GeV cm$^{-3}$ with a 1-$\sigma$ error bar of about 7%; more precisely we find a $\rho_{DM}(R_0) = 0.385 \pm 0.027 \rm GeV cm^{-3}$ for the Einasto profile and $\rho_{DM}(R_0) = 0.389 \pm 0.025 \rm GeV cm^{-3}$ for the NFW. This is in contrast to the standard assumption that $\rho_{DM}(R_0)$ is about 0.3 GeV cm$^{-3}$ with an uncertainty of a factor of 2 to 3. A very precise determination of the local halo density is very important for interpreting direct dark matter detection experiments. Indeed the results we produced, together with the recent accurate determination of the local circular velocity, should be very useful to considerably narrow astrophysical uncertainties on direct dark matter detection.

This Conference is the fifth of a series of Workshops on High Energy Gamma- ray Experiments, following the Conferences held in Perugia 2003, Bari 2004, Cividale del Friuli 2005, Elba Island 2006. This year the focus was on the use of gamma-ray to study the Dark Matter component of the Universe, the origin and propagation of Cosmic Rays, Extra Large Spatial Dimensions and Tests of Lorentz Invariance.

Scalar-Tensor theories of gravity can be formulated in different frames, most notably, the Einstein and the Jordan one. While some debate still persists in the literature on the physical status of the different frames, a frame transformation in Scalar-Tensor theories amounts to a local redefinition of the metric, and then should not affect physical results. We analyze the issue in a cosmological context. In particular, we define all the relevant observables (redshift, distances, cross-sections, ...) in terms of frame-independent quantities. Then, we give a frame-independent formulation of the Boltzmann equation, and outline its use in relevant examples such as particle freeze-out and the evolution of the CMB photon distribution function. Finally, we derive the gravitational equations for the frame-independent quantities at first order in perturbation theory. From a practical point of view, the present approach allows the simultaneous implementation of the good aspects of the two frames in a clear and straightforward way.

We consider a locally supersymmetric theory where the Planck mass is replaced by a dynamical superfield. This model can be thought of as the Minimal Supersymmetric extension of the Brans-Dicke theory (MSBD). The motivation that underlies this analysis is the research of possible connections between Dark Energy models based on Brans-Dicke-like theories and supersymmetric Dark Matter scenarios. We find that the phenomenology associated with the MSBD model is very different compared to the one of the original Brans-Dicke theory: the gravitational sector does not couple to the matter sector in a universal metric way. This feature could make the minimal supersymmetric extension of the BD idea phenomenologically inconsistent.

A host of dark energy models and non-standard cosmologies predict an enhanced Hubble rate in the early Universe: perfectly viable models, which satisfy Big Bang Nucleosynthesis (BBN), cosmic microwave background and general relativity tests, may nevertheless lead to enhancements of the Hubble rate up to many orders of magnitude. In this paper we show that strong bounds on the pre-BBN evolution of the Universe may be derived, under the assumption that dark matter is a thermal relic, by combining the dark matter relic density bound with constraints coming from the production of cosmic-ray antiprotons by dark matter annihilation in the Galaxy. The limits we derive can be sizable and apply to the Hubble rate around the temperature of dark matter decoupling. For dark matter masses lighter than 100 GeV, the bound on the Hubble-rate enhancement ranges from a factor of a few to a factor of 30, depending on the actual cosmological model, while for a mass of 500 GeV the bound falls in the range 50-500. Uncertainties in the derivation of the bounds and situations where the bounds become looser are discussed. We finally discuss how these limits apply to some specific realizations of non-standard cosmologies: a scalar-tensor gravity model, kination models and a Randall-Sundrum D-brane model.

Using the narrow clustering of the geometrically corrected gamma-ray energies released by gamma-ray bursts, we investigate the possibility of using these sources as standard candles to probe cosmological parameters such as the matter density Omega_m and the cosmological constant energy density Omega_Lambda. By simulating different samples of gamma-ray bursts, we find that Omega_m can be determined with accuracy ~7% with data from 300 sources. We also show that, if Omega = 1 is due to a quintessence field, some of the models proposed in the literature may be discriminated from a Universe with cosmological constant, by a similar-sized sample of gamma-ray bursts.

A scalar field interacting differently with dark matter and baryons may explain why their ratio is of order unity today. We provide three working examples, checking them against the observations of CMB, Large Scale Structure, supernovae Ia, and post-newtonian tests of gravity. Such a scenario could make life much easier for supersymmetric dark matter candidates.

Scalar-tensor theories of gravity provide a consistent framework to accommodate an ultra-light quintessence scalar field. While the equivalence principle is respected by construction, deviations from General Relativity and standard cosmology may show up at nucleosynthesis, CMB, and solar system tests of gravity. After imposing all the bounds coming from these observations, we consider the expansion rate of the universe at WIMP decoupling, showing that it can lead to an enhancement of the dark matter relic density up to few orders of magnitude with respect to the standard case. This effect can have an impact on supersymmetric candidates for dark matter.

In Scalar-Tensor theories of gravity, the expansion rate of the universe in the past may differ considerably from the standard one. After imposing the constraints coming from nucleosynthesis, CMB, type Ia supernovae, and solar system tests of General Relativity, we consider the expansion rate of the universe at WIMP decoupling, showing that it can lead to an enhancement of the dark matter relic density of some orders of magnitude with respect to the standard case. This effect may have a deep impact on most popular candidates for dark matter, as the supersymmetric neutralino.