Riccardo Catena

A novel determination of the local dark matter density

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 8dimensional 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-��

Dark matter relic abundance and scalar-tensor dark energy

Scalar-tensor theories of gravity provide a consistent framework to accommodate an ultralight quintessence scalar field. While the equivalence principle is respected by construction, deviations from general relativity and standard cosmology may show up at nucleosynthesis, cosmic microwave background, and solar system tests of gravity. After imposing all the bounds coming from these observations, we consider the expansion rate of the Universe at weakly interacting massive particle 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.

Einstein and Jordan frames reconciled: A frame-invariant approach to scalar-tensor cosmology

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 cosmic microwave background 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.

The local dark matter phase-space density and impact on WIMP direct detection

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, Eddingtons 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, regardless of the assumed profile, 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.

Form factors for dark matter capture by the Sun in effective theories

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, 3He, 4He, 12C, 14N, 16O, 20Ne, 23Na, 24Mg, 27Al, 28Si, 32S, 40Ar, 40Ca, 56Fe, and 59Ni, 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.

Anisotropic dark matter distribution functions and impact on WIMP direct detection

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.

Global fits of the dark matter-nucleon effective interactions

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.

Prospects for direct detection of dark matter in an effective theory approach

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.

DAMA confronts null searches in the effective theory of dark matter-nucleon interactions

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.

Dark matter signals at neutrino telescopes in effective theories

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.