Search

Models in which scalar field dark energy interacts with dark matter via a pure momentum coupling have previously been found to potentially ease the structure formation tension between early and late-Universe observations. In this paper, we explore the physical mechanism underlying this feature. We argue analytically that the perturbation growth equations imply the suppression of structure growth, illustrating our discussion with numerical calculations. Then we generalize the previously studied quadratic coupling between the dark energy and dark matter to a more general power-law case, also allowing for the slope of the dark energy exponential potential to vary. We find that the structure growth suppression is a generic feature of power-law couplings and it can, for a range of parameter values, be larger than previously found.

A small fraction, $f$, of cosmic string loops can collapse to form Primordial Black Holes (PBHs). Constraints on the abundance of PBHs can therefore be used to constrain $f$. We update these calculations, taking into account the PBH extended mass function, and find $f < 10^{-31} (G ��/ c^2)^{-3/2}$. This is roughly two orders of magnitude tighter than previous constraints. The improvement from the tighter constraints on the abundance of PBHs is partly off-set by refinements to the theoretical calculation of the cosmic string loop formation rate., 5 pages, v2: takes into account only a fraction of string network energy goes into long loops, constraints consequently weakened by one order of magnitude (thank you to referee for pointing this out) + other minor changes

We investigate cosmological models in which dynamical dark energy consists of a scalar field whose present-day value is controlled by a coupling to the neutrino sector. The behaviour of the scalar field depends on three functions: a kinetic function, the scalar field potential, and the scalar field-neutrino coupling function. We present an analytic treatment of the background evolution during radiation- and matter-domination for exponential and inverse power law potentials, and find a relaxation of constraints compared to previous work on the amount of early dark energy in the exponential case. We then carry out a numerical analysis of the background cosmology for both types of potential and various illustrative choices of the kinetic and coupling functions. By applying bounds from Planck on the amount of early dark energy, we are able to constrain the magnitude of the kinetic function at early times., 10 pages, 6 figures; added references, corrected typos, added intermediate steps to calculation in appendix

Cosmological observations indicate that most of the matter in the Universe is Dark Matter. Dark Matter in the form of Weakly Interacting Massive Particles (WIMPs) can be detected directly, via its elastic scattering off target nuclei. Most current direct detection experiments only measure the energy of the recoiling nuclei. However, directional detection experiments are sensitive to the direction of the nuclear recoil as well. Due to the Sun's motion with respect to the Galactic rest frame, the directional recoil rate has a dipole feature, peaking around the direction of the Solar motion. This provides a powerful tool for demonstrating the Galactic origin of nuclear recoils and hence unambiguously detecting Dark Matter. Furthermore, the directional recoil distribution depends on the WIMP mass, scattering cross section and local velocity distribution. Therefore, with a large number of recoil events it will be possible to study the physics of Dark Matter in terms of particle and astrophysical properties. We review the potential of directional detectors for detecting and characterizing WIMPs., 57 pages, 23 figures, to appear in Physics Reports

Primordial Black Holes (PBHs) can form in the radiation dominated early Universe from the collapse of large density perturbations produced by inflation. A power-law parameterisation of the primordial power spectrum is often used to extrapolate from cosmological scales, where the amplitude of the perturbations is well-measured by Cosmic Microwave Background and Large Scale Structure observations, down to the small scales on which PBHs may form. We show that this typically leads to large errors in the amplitude of the fluctuations on small scales, and hence extremely inaccurate calculations of the abundance of PBHs formed., 7 pages, 4 figures, v2: version to appear in Phys. Rev. D with minor changes to text

We study CMB constraints on a model with a cosmological constant and a fraction of dark matter non-minimally coupled to a massless scalar field. In this scenario, there is an extra gravity-like fifth force which can affect the evolution of the Universe enough to have a discernible effect on measurements of cosmological parameters. Using Planck and WMAP polarisation data, we find that up to half of the dark matter can be coupled. The coupling can also be several times larger than in models with a single species of cold dark matter coupled to a quintessence scalar field, as the scalar field does not play the role of dark energy and is therefore less constrained by the data., Comment: 5 pages, 4 figures

The analysis of signals in directional dark matter (DM) detectors typically assumes that the directions of nuclear recoils can be measured in the Galactic rest frame. However, this is not possible with all directional detection technologies. In nuclear emulsions, for example, the recoil events must be detected and measured after the exposure time of the experiment. Unless the entire detector is mounted and rotated with the sidereal day, the recoils cannot be reoriented in the Galactic rest frame. We examine the effect of this `time integration' on the primary goals of directional detection, namely: (1) confirming that the recoils are anisotropic; (2) measuring the median recoil direction to confirm their Galactic origin; and (3) probing below the neutrino floor. We show that after time integration the DM recoil distribution retains a preferred direction and is distinct from that of Solar neutrino-induced recoils. Many of the advantages of directional detection are therefore preserved and it is not crucial to mount and rotate the detector. Rejecting isotropic backgrounds requires a factor of 2 more signal events compared with an experiment with event time information, whereas a factor of 1.5-3 more events are needed to measure a median direction in agreement with the expectation for DM. We also find that there is still effectively no neutrino floor in a time-integrated directional experiment. However to reach a cross section an order of magnitude below the floor, a factor of 8 larger exposure is required than with a conventional directional experiment. We also examine how the sensitivity is affected for detectors with only 2D recoil track readout, and/or no head-tail measurement. As for non-time-integrated experiments, 2D readout is not a major disadvantage, though a lack of head-tail sensitivity is., 15 pages, 11 figures. Version published in PRD

Terrestrial searches for the conversion of dark matter axions or axion-like particles into photons inside magnetic fields are sensitive to the phase space structure of the local Milky Way halo. We simulate signals in a hypothetical future experiment based on the Axion Dark Matter eXperiment (ADMX) that could be performed once the axion has been detected and a frequency range containing the axion mass has been identified. We develop a statistical analysis to extract astrophysical parameters, such as the halo velocity dispersion and laboratory velocity, from such data and find that with only a few days integration time a level of precision can be reached matching that of astronomical observations. For longer experiments lasting up to a year in duration we find that exploiting the modulation of the power spectrum in time allows accurate measurements of the Solar peculiar velocity with an accuracy that would improve upon astronomical observations. We also simulate signals based on results from N-body simulations and find that finer substructure in the form of tidal streams would show up prominently in future data, even if only a subdominant contribution to the local dark matter distribution. In these cases it would be possible to reconstruct all the properties of a dark matter stream using the time and frequency dependence of the signal. Finally we consider the detection prospects for a network of streams from tidally disrupted axion miniclusters. These features appear much more prominently in the resolved spectrum than suggested by calculations based on a scan over a range of resonant frequencies, making the detection of axion minicluster streams more viable than previously thought. These results confirm that haloscope experiments in a post-discovery era are able to perform "axion astronomy"., Comment: 17 pages, 8 figures, references and minor detail added. Matches published version

There has recently been interest in multi-Solar mass Primordial Black Holes (PBHs) as a dark matter (DM) candidate. There are various microlensing, dynamical and accretion constraints on the abundance of PBHs in this mass range. Taken at face value these constraints exclude multi-Solar mass PBHs making up all of the DM for both delta-function and extended mass functions. However the stellar microlensing event rate depends on the density and velocity distribution of the compact objects along the line of sight to the Magellanic Clouds. We study the dependence of the constraints on the local dark matter density and circular speed and also consider models where the velocity distribution varies with radius. We find that the largest mass constrained by stellar microlensing can vary by an order of magnitude. In particular the constraints are significantly weakened if the velocity dispersion of the compact objects is reduced. The change is not sufficiently large to remove the tension between the stellar microlensing and dynamical constraints. However this demonstrates that it is crucial to take into account astrophysical uncertainties when calculating and comparing constraints. We also confirm the recent finding that the tension between the constraints is in fact increased for realistic, finite width mass functions., Comment: 6 pages, 5 figures, v2: version to appear in PRD with cosmetic changes to figs. 1 & 4

The differential event rate in Weakly Interacting Massive Particle (WIMP) direct detection experiments depends on the local dark matter density and velocity distribution. Accurate modelling of the local dark matter distribution is therefore required to obtain reliable constraints on the WIMP particle physics properties. Data analyses typically use a simple Standard Halo Model which might not be a good approximation to the real Milky Way (MW) halo. We review observational determinations of the local dark matter density, circular speed and escape speed and also studies of the local dark matter distribution in simulated MW-like galaxies. We discuss the effects of the uncertainties in these quantities on the energy spectrum and its time and direction dependence. Finally we conclude with an overview of various methods for handling these astrophysical uncertainties., Comment: 10 pages, v2: to appear in JPhysG focus issue on Dark Matter, minor changes

The recent discovery of gravitational waves from mergers of $\sim 10 \, M_{\odot}$ black hole binaries has stimulated interested in Primordial Black Hole dark matter in this mass range. Microlensing and dynamical constraints exclude all of the dark matter being in compact objects with a delta function mass function in the range $10^{-7} \lesssim M/ M_{\odot} \lesssim 10^{5}$. However it has been argued that all of the dark matter could be composed of compact objects in this range with an extended mass function. We explicitly recalculate the microlensing and dynamical constraints for compact objects with an extended mass function which replicates the PBH mass function produced by inflation models. We find that the microlensing and dynamical constraints place conflicting constraints on the width of the mass function, and do not find a mass function which satisfies both constraints., Comment: 6 pages, 4 figures, v2: version published in Phys. Rev. D with additional references

When comparing constraints on the weakly interacting massive particle (WIMP) properties from direct and indirect detection experiments it is crucial that the assumptions made about the dark matter (DM) distribution are realistic and consistent. For instance, if the Fermi-LAT Galactic center GeV gamma-ray excess was due to WIMP annihilation, its morphology would be incompatible with the standard halo model that is usually used to interpret data from direct detection experiments. In this article, we calculate exclusion limits from direct detection experiments using self-consistent velocity distributions, derived from mass models of the Milky Way where the DM halo has a generalized Navarro-Frenk-White profile. We use two different methods to make the mass model compatible with a DM interpretation of the Galactic center gamma-ray excess. First, we fix the inner slope of the DM density profile to the value that best fits the morphology of the excess. Second, we allow the inner slope to vary and include the morphology of the excess in the data sets used to constrain the gravitational potential of the Milky Way. The resulting direct detection limits differ significantly from those derived using the standard halo model, in particular for light WIMPs, due to the differences in both the local DM density and velocity distribution, D. G. C. is funded by the STFC. M. F. gratefully acknowledges support from the Netherlands Organization for Scientific Research (NWO) through a Vidi grant (P. I.: Dr. Shin’ichiro Ando). A. M. G. acknowledges support from STFC Grant No. ST/L000393/1. M. P. is supported under the ERC Advanced Grant SPLE under Contract No. ERC-2012-ADG-20120216-320421. We acknowledge support of the Consolider-Ingenio 2010 program under Grant No. MULTIDARK CSD2009- 00064, the Spanish MICINN under Grant No. FPA2013- 44773-P and the Spanish MINECO Centro de Excelencia Severo Ochoa Program under Grant No. SEV-2012-0249

The direct detection of dark matter on Earth depends crucially on its density and its velocity distribution on a milliparsec scale. Conventional N-body simulations are unable to access this scale, making the development of other approaches necessary. In this paper, we apply the method developed in Fantin et al. 2008 to a cosmologically-based merger tree, transforming it into a useful instrument to reproduce and analyse the merger history of a Milky Way-like system. The aim of the model is to investigate the implications of any ultra-fine structure for the current and next generation of directional dark matter detectors. We find that the velocity distribution of a Milky Way-like Galaxy is almost smooth, due to the overlap of many streams of particles generated by multiple mergers. Only the merger of a 10^10 Msun analyse can generate significant features in the ultra-local velocity distribution, detectable at the resolution attainable by current experiments.

The search for weakly interacting massive particles (WIMPs) by direct detection faces an encroaching background due to coherent neutrino-nucleus scattering. As the sensitivity of these experiments improves, the question of how to best distinguish a dark matter signal from neutrinos will become increasingly important. A proposed method of overcoming this so-called 'neutrino floor' is to utilize the directional signature that both neutrino and dark matter induced recoils possess. We show that directional experiments can indeed probe WIMP-nucleon cross-sections below the neutrino floor with little loss in sensitivity due to the neutrino background. In particular we find at low WIMP masses (around 6 GeV) the discovery limits for directional detectors penetrate below the non-directional limit by several orders of magnitude. For high WIMP masses (around 100 GeV), the non-directional limit is overcome by a factor of a few. Furthermore we show that even for directional detectors which can only measure 1- or 2-dimensional projections of the 3-dimensional recoil track, the discovery potential is only reduced by a factor of 3 at most. We also demonstrate that while the experimental limitations of directional detectors, such as sense recognition and finite angular resolution, have a detrimental effect on the discovery limits, it is still possible to overcome the ultimate neutrino background faced by non-directional detectors., 16 pages, 12 figures, published version

With positive signals from multiple direct detection experiments it will, in principle, be possible to measure the mass and cross sections of weakly-interacting massive particle (WIMP) dark matter. Recent work has shown that, with a polynomial parameterisation of the WIMP speed distribution, it is possible to make an unbiased measurement of the WIMP mass, without making any astrophysical assumptions. However, direct detection experiments are not sensitive to low-speed WIMPs and, therefore, any model-independent approach will lead to a bias in the cross section. This problem can be solved with the addition of measurements of the flux of neutrinos from the Sun. This is because the flux of neutrinos produced from the annihilation of WIMPs which have been gravitationally captured in the Sun is sensitive to low-speed WIMPs. Using mock data from next-generation direct detection experiments and from the IceCube neutrino telescope, we show that the complementary information from IceCube on low-speed WIMPs breaks the degeneracy between the cross section and the speed distribution. This allows unbiased determinations of the WIMP mass and spin-independent and spin-dependent cross sections to be made, and the speed distribution to be reconstructed. We use two parameterisations of the speed distribution: binned and polynomial. While the polynomial parameterisation can encompass a wider range of speed distributions, this leads to larger uncertainties in the particle physics parameters., Comment: 24 pages, 13 figures. Corrected error in likelihood; changes to discussion of Benchmark C. Matches version accepted in PRD

The WIMP direct detection event rate (and in particular its annual modulation) depends on the local dark matter distribution and also the Earth’s velocity with respect to the Galactic rest frame. I discuss the eect of these astrophysical inputs on the expected WIMP signals, focussing on exclusions limits and the size and phase of the annual modulation signal.

The formation of large scale structure is independent of the nature of the cold dark matter (CDM), however the fate of very small scale inhomogeneities depends on the micro-physics of the CDM particles. We investigate the matter power spectrum for scales that enter the Hubble radius well before matter-radiation equality, and follow its evolution until the time when the first inhomogeneities become non-linear. Our focus lies on weakly interacting massive particles (WIMPs), and as a concrete example we analyze the case when the lightest supersymmetric particle is a bino. We show that collisional damping and free-streaming of WIMPs lead to a matter power spectrum with a sharp cut-off at about 10^(-6) M_odot and a maximum close to that cut-off. We also calculate the transfer function for the growth of the inhomogeneities in the linear regime. These three effects (collisional damping, free-streaming and gravitational growth) are combined to provide a WMAP normalized primordial CDM power spectrum, which could serve as an input for high resolution CDM simulations. The smallest inhomogeneities typically enter the non-linear regime at a redshift of about 60.

Weakly Interacting Massive Particle (WIMP) direct detection experiments are just reaching the sensitivity required to detect Galactic dark matter in the form of neutralinos (or indeed any stable weakly interacting particle). Detection strategies and data analyses are often based on the simplifying assumption of a standard spherical, isothermal halo model, but observations and numerical simulations indicate that galaxy halos are in fact triaxial and anisotropic, and contain substructure. The annual modulation and direction dependence of the event rate (due to the motion of the Earth) provide the best prospects of distinguishing WIMP scattering from background events, however these signals depend sensitively on the local WIMP velocity distribution. I briefly review the status of WIMP direct detection experiments before discussing the dependence of the annual modulation signal on astrophysical input, in particular the structure of the Milky Way halo, and the possibility that the local WIMP distribution is not smooth., 6 pages, 2 figures, to appear in the proceedings of IAU Symposium 220 "Dark matter in galaxies", ASP, Eds: S. Ryder, D.J. Pisano, M. Walker, K. Freeman

Primordial Black Holes (PBHs) are, typically light, black holes which can form in the early Universe. There are a number of formation mechanisms, including the collapse of large density perturbations, cosmic string loops and bubble collisions. The number of PBHs formed is tightly constrained by the consequences of their evaporation and their lensing and dynamical effects. Therefore PBHs are a powerful probe of the physics of the early Universe, in particular models of inflation. They are also a potential cold dark matter candidate.

Directional detection of WIMPs, in which the energies and directions of the recoiling nuclei are measured, currently presents the only prospect for probing the local velocity distribution of Galactic dark matter. We investigate the extent to which future directional detectors would be capable of probing dark matter substructure in the form of streams. We analyse the signal expected from a Sagittarius-like stream and also explore the full parameter space of stream speed, direction, dispersion and density. Using a combination of non-parametric directional statistics, a profile likelihood ratio test and Bayesian parameter inference we find that within acceptable exposure times (O(10) kg yr for cross sections just below the current exclusion limits) future directional detectors will be sensitive to a wide range of stream velocities and densities. We also examine and discuss the importance of the energy window of the detector., 18 pages, 12 figures, 1st revision, submitted to Phys. Rev. D

A novel method for characterising the full spectrum of deuteron ions emitted by laser driven multi-species ion sources is discussed. The procedure is based on using differential filtering over the detector of a Thompson parabola ion spectrometer, which enables discrimination of deuterium ions from heavier ion species with the same charge-to-mass ratio (such as C6+, O8+, etc.). Commonly used Fuji Image plates were used as detectors in the spectrometer, whose absolute response to deuterium ions over a wide range of energies was calibrated by using slotted CR-39 nuclear track detectors. A typical deuterium ion spectrum diagnosed in a recent experimental campaign is presented., 7 pages, 7 figures

Dark Matter (DM) direct detection experiments usually assume the simplest possible 'Standard Halo Model' for the Milky Way (MW) halo in which the velocity distribution is Maxwellian. This model assumes that the MW halo is an isotropic, isothermal sphere, hypotheses that are unlikely to be valid in reality. An alternative approach is to derive a self-consistent solution for a particular mass model of the MW (i.e. obtained from its gravitational potential) using the Eddington formalism, which assumes isotropy. In this paper we extend this approach to incorporate an anisotropic phase-space distribution function. We perform Bayesian scans over the parameters defining the mass model of the MW and parameterising the phase-space density, implementing constraints from a wide range of astronomical observations. The scans allow us to estimate the precision reached in the reconstruction of the velocity distribution (for different DM halo profiles). As expected, allowing for an anisotropic velocity tensor increases the uncertainty in the reconstruction of f(v) but the distribution can still be determined with a precision of a factor of 4-5. The mean velocity distribution resembles the isotropic case, however the amplitude of the high-velocity tail is up to a factor of 2 larger. Our results agree with the phenomenological parametrization proposed in Mao et al. (2013) as a good fit to N-body simulations (with or without baryons), since their velocity distribution is contained in our 68% credible interval., 17 pages, 6 figures

Primordial black holes may form in the early Universe, for example from the collapse of large amplitude density perturbations predicted in some inflationary models. Light black holes undergo Hawking evaporation, the energy injection from which is constrained both at the epoch of nucleosynthesis and at the present. The failure as yet to unambiguously detect primordial black holes places important constraints. In this article, we are particularly concerned with the dependence of these constraints on the model for the complete cosmological history, from the time of formation to the present. Black holes presently give the strongest constraint on the spectral index $n$ of density perturbations, though this constraint does require $n$ to be constant over a very wide range of scales., 8 pages LaTeX file, using elsart.sty, with three figures incorporated using epsf. To appear, proceedings of DM98, Los Angeles (ed D Cline, Elsevier)

Many models have been studied that contain more than one species of dark matter and some of these couple the Cold Dark Matter (CDM) to a light scalar field. In doing this we introduce additional long range forces, which in turn can significantly affect our estimates of cosmological parameters if not properly accounted for. It is, therefore, important to study these models and their resulting cosmological implications. We present a model in which a fraction of the total cold dark matter density is coupled to a scalar field. We study the background and perturbation evolution and calculate the resulting Cosmic Microwave Background anisotropy spectra. The greater the fraction of dark matter coupled to the scalar field and the stronger the coupling strength, the greater the deviation of the background evolution from LCDM. Previous work, with a single coupled dark matter species, has found an upper limit on the coupling strength of order O(0.1). We find that with a coupling of this magnitude more than half the dark matter can be coupled to a scalar field without producing any significant deviations from LCDM., 10 pages, 4 figures, minor changes

Determining the dark matter (DM) mass is of paramount importance for understanding dark matter. We present a novel parametrization of the DM speed distribution which will allow the DM mass to be accurately measured using data from Weakly Interacting Massive Particle (WIMP) direct detection experiments. Specifically, we parametrize the natural logarithm of the speed distribution as a polynomial in the speed v. We demonstrate, using mock data from upcoming experiments, that by fitting the WIMP mass and interaction cross-section, along with the polynomial coefficients, we can accurately reconstruct both the WIMP mass and speed distribution. This new method is the first demonstration that an accurate, unbiased reconstruction of the WIMP mass is possible without prior assumptions about the distribution function. We anticipate that this technique will be invaluable in the analysis of future experimental data., 5 pages, 3 figures. Minor changes. Matches version published in PRL

The search for weakly-interacting massive particle (WIMP) dark matter is multi-pronged. Ultimately, the WIMP-dark-matter picture will only be confirmed if different classes of experiments see consistent signals and infer the same WIMP properties. In this work, we review the ideas, methods, and status of direct-detection searches. We focus in particular on extracting WIMP physics (WIMP interactions and phase-space distribution) from direct-detection data in the early discovery days when multiple experiments see of order dozens to hundreds of events. To demonstrate the essential complementarity of different direct-detection experiments in this context, we create mock data intended to represent the data from the near-future Generation 2 experiments. We consider both conventional supersymmetry-inspired benchmark points (with spin-independent and -dependent elastic cross sections just below current limits), as well as benchmark points for other classes of models (inelastic and effective-operator paradigms). We also investigate the effect on parameter estimation of loosening or dropping the assumptions about the local WIMP phase-space distribution. We arrive at two main conclusions. Firstly, teasing out WIMP physics with experiments depends critically on having a wide set of detector target materials, spanning a large range of target nuclear masses and spin-dependent sensitivity. It is also highly desirable to obtain data from low-threshold experiments. Secondly, a general reconstruction of the local WIMP velocity distribution, which will only be achieved if there are multiple experiments using different target materials, is critical to obtaining a robust and unbiased estimate of the WIMP mass., Comment: 39 pages, 26 figures, submitted to Physics of the Dark Universe

Primordial Black Holes (PBH's) can form in the early Universe from the collapse of large density fluctuations. Tight observational limits on their abundance constrain the amplitude of the primordial fluctuations on very small scales which can not otherwise be constrained, with PBH's only forming from the extremely rare large fluctuations. The number of PBH's formed is therefore sensitive to small changes in the shape of the tail of the fluctuation distribution, which itself depends on the amount of non-Gaussianity present. We study, for the first time, how quadratic and cubic local non-Gaussianity of arbitrary size (parameterised by f_nl and g_nl respectively) affects the PBH abundance and the resulting constraints on the amplitude of the fluctuations on very small scales. Intriguingly we find that even non-linearity parameters of order unity have a significant impact on the PBH abundance. The sign of the non-Gaussianity is particularly important, with the constraint on the allowed fluctuation amplitude tightening by an order of magnitude as f_nl changes from just -0.5 to 0.5. We find that if PBH's are observed in the future, then regardless of the amplitude of the fluctuations, non-negligible negative f_nl would be ruled out. Finally we show that g_nl can have an even larger effect on the number of PBH's formed than f_nl., 9 pages, 5 figures, v2: version to appear in Phys. Rev. D with minor changes, v3: typos corrected (including factor of 1/2 in erfc prefactor), no changes to results

We study the effect of curvaton decay on the primordial power spectrum. Using analytical approximations and also numerical calculations, we find that the power spectrum is enhanced during the radiation dominated era after the curvaton decay. The amplitude of the Bardeen potential is controlled by the fraction of the energy density in the curvaton at the time of curvaton decay. We show that the enhancement in the amplitude of the primordial curvature perturbation is, however, not large enough to lead to primordial black hole overproduction on scales which re-enter the horizon after the time of curvaton decay., Comment: 13 pages, 6 figures

The CoGeNT and CRESST WIMP direct detection experiments have recently observed excesses of nuclear recoil events, while the DAMA/LIBRA experiment has a long standing annual modulation signal. It has been suggested that these excesses may be due to light mass, m_chi ~ 5-10 GeV, WIMPs. The Earth's motion with respect to the Galactic rest frame leads to a directional dependence in the WIMP scattering rate, providing a powerful signal of the Galactic origin of any recoil excess. We investigate whether direct detection experiments with directional sensitivity have the potential to observe this anisotropic scattering rate with the elastically scattering light WIMPs proposed to explain the observed excesses. We find that the number of recoils required to detect an anisotropic signal from light WIMPs at 5 sigma significance varies from 7 to more than 190 over the set of target nuclei and energy thresholds expected for directional detectors. Smaller numbers arise from configurations where the detector is only sensitive to recoils from the highest speed, and hence most anisotropic, WIMPs. However, the event rate above threshold is very small in these cases, leading to the need for large experimental exposures to accumulate even a small number of events. To account for this sensitivity to the tail of the WIMP velocity distribution, whose shape is not well known, we consider two exemplar halo models spanning the range of possibilities. We also note that for an accurate calculation the Earth's orbital speed must be averaged over. We find that the exposures required to detect 10 GeV WIMPs at a WIMP-proton cross-section of 10^-4 pb are of order 10^3 kg day for a 20 keV energy threshold, within reach of planned directional detectors. Lower WIMP masses require higher exposures and/or lower energy thresholds for detection., Comment: 9 pages, 6 figures, v2: version to appear in Phys. Rev. D with additional discussion

A sufficiently light scalar field slowly evolving in a potential can account for the dark energy that presently dominates the universe. This quintessence field is expected to couple directly to matter components, unless some symmetry of a more fundamental theory protects or suppresses it. Such a coupling would leave distinctive signatures in the background expansion history of the universe and on cosmic structure formation, particularly at galaxy cluster scales. Using semi--analytic expressions for the CDM halo mass function, we make predictions for halo abundance in models where the quintessence scalar field is coupled to cold dark matter, for a variety of quintessence potentials. We evaluate the linearly extrapolated density contrast at the redshift of collapse using the spherical collapse model and we compare this result to the corresponding prediction obtained from the non--linear perturbation equations in the Newtonian limit. For all the models considered in this work, if there is a continuous flow of energy from the quintessence scalar field to the CDM component, then the predicted number of CDM haloes can only lie below that of $\Lambda$CDM, when each model shares the same cosmological parameters today. In the last stage of our analysis we perform a global MCMC fit to data to find the best fit values for the cosmological model parameters. We find that for some forms of the quintessence potential, coupled dark energy models can offer a viable alternative to $\Lambda$CDM in light of the recent detections of massive high--$z$ galaxy clusters, while other models of coupled quintessence predict a smaller number of massive clusters at high redshift compared to $\Lambda$CDM., Comment: 18 pages, 10 figures. Minor typos corrected, discussion extended. Version to appear in PRD

The signals expected in WIMP direct detection experiments depend on the ultra-local dark matter distribution. Observations probe the local density, circular speed and escape speed, while simulations find velocity distributions that deviate significantly from the standard Maxwellian distribution. We calculate the energy, time and direction dependence of the event rate for a range of velocity distributions motivated by recent observations and simulations, and also investigate the uncertainty in the determination of WIMP parameters. The dominant uncertainties are the systematic error in the local circular speed and whether or not the MW has a high density dark disc. In both cases there are substantial changes in the mean differential event rate and the annual modulation signal, and hence exclusion limits and determinations of the WIMP mass. The uncertainty in the shape of the halo velocity distribution is less important, however it leads to a 5% systematic error in the WIMP mass. The detailed direction dependence of the event rate is sensitive to the velocity distribution. However the numbers of events required to detect anisotropy and confirm the median recoil direction do not change substantially., 21 pages, 7 figures, v2 version to appear in JCAP, minor changes

Direct detection experiments have reached the sensitivity required to detect dark matter WIMPs. Demonstrating that a putative signal is due to WIMPs, and not backgrounds, is a major challenge however. The direction dependence of the WIMP scattering rate provides a potential WIMP `smoking gun'. If the WIMP distribution is predominantly smooth, the Galactic recoil distribution is peaked in the direction opposite to the direction of Solar motion. Previous studies have found that, for an ideal detector, of order 10 WIMP events would be sufficient to reject isotropy, and rule out an isotropic background. We examine how the median recoil direction could be used to confirm the WIMP origin of an anisotropic recoil signal. Specifically we determine the number of events required to reject the null hypothesis that the median direction is random (corresponding to an isotropic Galactic recoil distribution) at 95% confidence. We find that for zero background 31 events are required, a factor of roughly 2 more than are required to simply reject isotropy. We also investigate the effect of a non-zero isotropic background. As the background rate is increased the number of events required increases, initially fairly gradually and then more rapidly, once the signal becomes subdominant. We also discuss the effect of features in the speed distribution at large speeds, as found in recent high resolution simulations, on the median recoil direction., Comment: 4 pages, 2 figures, version to appear in Phys. Rev. D as rapid communication, minor changes

Dark matter plays a fundamental role in theories of formation and evolution of galaxies, so every attempt to model galaxy formation and evolution has to consider the presence of dark halos. Moreover, mergers and accretions appear to be fundamental driving mechanisms in determining the present day properties of galaxies. Our project studies the ultra‐fine distribution of dark matter in the Milky Way and the implications for the next generation of dark matter particle detectors. We have developed a model for halo mergers in a Milky Way‐like galaxy. The signals expected in the lab‐based dark matter detection experiments depend on the phase‐space distribution on submilliparsec scales: with our numerical technique it is possible to resolve structures, produced by mergers of subhalos with a larger parent halo, inaccessible to conventional N‐body simulations. When applied in a cosmological context, this method becomes a powerful instrument to reproduce and analyse the complete multiple merger history of a Milky ...

Ultracompact minihalos (UCMHs) are dense dark matter structures which can form from large density perturbations shortly after matter-radiation equality. If dark matter is in the form of Weakly Interacting Massive Particles (WIMPs), then UCMHs may be detected via their gamma-ray emission. We investigate how the {\em{Fermi}} satellite could constrain the abundance of UCMHs and place limits on the power spectrum of the primordial curvature perturbation. Detection by {\em Fermi} would put a lower limit on the UCMH halo fraction. The smallest detectable halo fraction, $f_{\rm UCMH} \gtrsim 10^{-7}$, is for $M_{\rm UCMH} \sim 10^{3} M_{\odot}$. If gamma-ray emission from UCMHs is not detected, an upper limit can be placed on the halo fraction. The bound is tightest, $f_{\rm UCMH} \lesssim 10^{-5}$, for $M_{\rm UCMH} \sim 10^{5} M_{\odot}$. The resulting upper limit on the power spectrum of the primordial curvature perturbation in the event of non-detection is in the range $\mathcal{P_R} \lesssim 10^{-6.5}- 10^{-6}$ on scales $k \sim 10^{1}-10^{6} \, {\rm Mpc}^{-1}$. This is substantially tighter than the existing constraints from primordial black hole formation on these scales, however it assumes that dark matter is in the form of WIMPs and UCMHs are not disrupted during the formation of the Milky Way halo., Comment: 5 pages, 2 figures, version to appear in Phys. Rev. D, minor changes

Primordial black hole (PBH) abundance limits constrain the primordial power spectrum, and hence models of inflation, on scales far smaller than those probed by cosmological observations. Single field inflation models which are compatible with all cosmological data can have large enough perturbations on small scales to overproduce PBHs, and hence be excluded. The standard formulae for the amplitude of perturbations do not hold for modes that exit the horizon close to the end of inflation however. We use a modified flow analysis to identify models of inflation where the amplitude of perturbations on small scales is large. For these models we then carry out a numerical evolution of the perturbations and use the PBH constraints on the power spectrum to eliminate models which overproduce PBHs. Significant PBH formation can occur in models in which inflation can continue indefinitely and is ended via a secondary mechanism. We demonstrate that PBHs constrain these types of inflation models and show that a numerical evaluation of the power spectrum decreases the number of otherwise viable models of inflation., Comment: 5 pages, 5 figures v3: version to appear in Phys. Rev. D with minor changes

We present the case for a dark matter detector with directional sensitivity. This document was developed at the 2009 CYGNUS workshop on directional dark matter detection, and contains contributions from theorists and experimental groups in the field. We describe the need for a dark matter detector with directional sensitivity; each directional dark matter experiment presents their project's status; and we close with a feasibility study for scaling up to a one ton directional detector, which would cost around $150M., 48 pages, 37 figures, whitepaper on direct dark matter detection with directional sensitivity

Primordial black holes (PBHs) can form in the early Universe via the collapse of large density perturbations. There are tight constraints on the abundance of PBHs formed due to their gravitational effects and the consequences of their evaporation. These abundance constraints can be used to constrain the primordial power spectrum, and hence models of inflation, on scales far smaller than those probed by cosmological observations. We compile, and where relevant update, the constraints on the abundance of PBHs before calculating the constraints on the curvature perturbation, taking into account the growth of density perturbations prior to horizon entry. We consider two simple parameterizations of the curvature perturbation spectrum on the scale of interest: constant and power-law. The constraints from PBHs on the amplitude of the power spectrum are typically in the range 10^{-2}-10^{-1} with some scale dependence., 10 pages, 2 figures, version to appear in Phys. Rev. D with minor change to calculation of constraints for spectral index not equal to one

Various laboratory-based experiments are underway attempting to detect dark matter directly. The event rates and detailed signals expected in these experiments depend on the dark matter phase space distribution on sub-milliparsec scales. These scales are many orders of magnitude smaller than those that can be resolved by conventional N-body simulations, so one cannot hope to use such tools to investigate the effect of mergers in the history of the Milky Way on the detailed phase-space structure probed by the current experiments. In this paper we present an alternative approach to investigating the results of such mergers, by studying a simplified model for a merger of a sub-halo with a larger parent halo. With an appropriate choice of parent halo potential, the evolution of material from the sub-halo can be expressed analytically in action-angle variables, so it is possible to obtain its entire orbit history very rapidly without numerical integration. Furthermore by evolving backwards in time, we can obtain arbitrarily-high spatial resolution for the current velocity distribution at a fixed point. Although this model cannot provide a detailed quantitative comparison with the Milky Way, its properties are sufficiently generic that it offers qualitative insight into the expected structure arising from a merger at a resolution that cannot be approached with full numerical simulations. Preliminary results indicate that the velocity-space distribution of dark matter particles remains characterized by discrete and well-defined peaks over an extended period of time, both for single and multi-merging systems, in contrast to the simple smooth velocity distributions sometimes assumed in predicting laboratory experiment detection rates., 6 pages, 6 figures, accepted for publication in MNRAS

We study the consequences of limited recoil sense reconstruction on the number of events required to reject isotropy and detect a WIMP signal using a directional detector. For a constant probability of determining the sense correctly, 3-d read-out and zero background, we find that as the probability is decreased from 1.0 to 0.75 the number of events required to reject isotropy using the mean angle statistic is increased by a factor of a few. As the probability is decreased further the number of events required using this statistic increases sharply, and in fact isotropy can be rejected more easily by discarding the sense information and using axial statistics. This however requires an order of magnitude more events than vectorial data with perfect sense determination. We also consider energy dependent probabilities of correctly measuring the sense, 2-d read-out and non-zero background. Our main conclusion regarding the sense determination is that correctly determining the sense of the abundant, but less anisotropic, low energy recoils is most important for minimising the number of events required., 5 pages, 3 figures, version to appear in PRD, minor changes

The energy spectrum of nuclear recoils in Weakly Interacting Massive Particle (WIMP) direct detection experiments depends on the underlying WIMP mass (strongly for light WIMPs, weakly for heavy WIMPs). We discuss how the accuracy with which the WIMP mass could be determined by a single direct detection experiment depends on the detector configuration and the WIMP properties. In particular we examine the effects of varying the underlying WIMP mass, the detector target nucleus, exposure, energy threshold and maximum energy, the local velocity distribution and the background event rate and spectrum., Comment: 4 pages, 1 figure, to appear in the proceedings of Identification of Dark Matter (idm2008), Stockholm, 18-22 August 2008

We study the dependence of the exposure required to directly detect a WIMP directional recoil signal on the capabilities of a directional detector. Specifically we consider variations in the nuclear recoil energy threshold, the background rate, whether the detector measures the recoil momentum vector in 2 or 3 dimensions and whether or not the sense of the momentum vector can be determined. We find that the property with the biggest effect on the required exposure is the measurement of the momentum vector sense. If the detector cannot determine the recoil sense, the exposure required is increased by an order of magnitude for 3-d read-out and two orders of magnitude for 2-d read-out. For 2-d read-out the required exposure, in particular if the senses can not be measured, can be significantly reduced by analyzing the reduced angles with the, time dependent, projected direction of solar motion subtracted. The background rate effectively places a lower limit on the WIMP cross-section to which the detector is sensitive; it will be very difficult to detect WIMPs with a signal rate more than an order of magnitude below the background rate. Lowering the energy threshold also reduces the required exposure, but only for thresholds above 20 keV., 18 pages, 7 figures, version to appear in Astroparticle Physics with minor changes

We calculate the constraints on the primordial curvature perturbation at the end of inflation from the present day abundance of Primordial Black Holes (PBHs), as a function of the reheat temperature T_{\rm RH}. We first extend recent work on the formation of PBHs on scales which remain within the horizon during inflation and calculate the resulting constraints on the curvature perturbation. We then evaluate the constraint from PBHs that form, more conventionally, from super-horizon perturbations. The constraints apply for T_{\rm RH} < 10^{8} GeV and the inclusion of sub-horizon PBHs leads to a limit which is roughly three times tighter than the bound from super-horizon PBHs., 8 pages, 2 figures, version to appear in JCAP, minor changes

We study, analytically and numerically, the energy input into dark matter mini-haloes by interactions with stars. We find that the fractional energy input in simulations of Plummer spheres agrees well with the impulse approximation for small and large impact parameters, with a rapid transition between these two regimes. Using the impulse approximation the fractional energy input at large impact parameters is fairly independent of the mass and density profile of the mini-halo, however low-mass mini-haloes experience a greater fractional energy input in close encounters. We formulate a fitting function which encodes these results and use it to estimate the disruption timescales of mini-haloes, taking into account the stellar velocity dispersion and mass distribution. For mini-haloes with mass M< {\cal O}(10^{-7} M_{\odot} on typical orbits which pass through the disc, we find that the estimated disruption timescales are independent of mini-halo mass, and are of order the age of the Milky Way. For more massive mini-haloes the estimated disruption timescales increase rapidly with increasing mass., 10 pages, 6 figures. Version to appear in MNRAS. Significant changes to presentation, but conclusions essentially unchanged

The dark matter distribution on small scales may depend on the properties of the first generation of dark matter halos to form, which is in turn determined by the microphysics of the dark matter particles. We overview the microphysics of WIMPs and the calculation of the collisional damping and free streaming scales. We then present the resulting density perturbation power spectrum and the red‐shift at which typical halos form, taking into account the effect of uncertainties in the WIMP properties (mass and interaction channel) and the primordial power spectrum. Finally we review simulations of the formation of the first WIMPy halos and discuss their subsequent evolution, in particular disruption due to encounters with stars.

The direction dependence of the WIMP direct detection rate provides a powerful tool for distinguishing a WIMP signal from possible backgrounds. We study the number of events required to discriminate a WIMP signal from an isotropic background for a detector with 2-d readout using nonparametric circular statistics. We also examine the number of events needed to (i) detect a deviation from rotational symmetry, due to flattening of the Milky Way halo and (ii) detect a deviation in the mean direction due to a tidal stream. If the senses of the recoils are measured then of order 20---70 events (depending on the plane of the 2-d readout and the detector location) will be sufficient to reject isotropy of the raw recoil angles at $90%$ confidence. If the senses can not be measured these number increase by roughly 2 orders of magnitude (compared with an increase of 1 order of magnitude for the case of full 3-d readout). The distributions of the reduced angles, with the (time-dependent) direction of solar motion subtracted, are far more anisotropic, however, and if the isotropy tests are applied to these angles then the numbers of events required are similar to the case of 3-d readout. A deviation from rotational symmetry will only be detectable if the Milky Way halo is significantly flattened. The deviation in the mean direction due to a tidal stream is potentially detectable, however, depending on the density and direction of the stream. The meridian plane (which contains the Earth's spin axis) is, for all detector locations, the optimum readout plane for rejecting isotropy. However readout in this plane can not be used for detecting flattening of the Milky Way halo or a stream with direction perpendicular to the galactic plane. In these cases the optimum readout plane depends on the detector location.

Chameleon fields are scalar fields whose mass depends on the ambient matter density. We investigate the effects of these fields on the growth of density perturbations on sub-galactic scales and the formation of the first dark matter halos. Density perturbations on comoving scales $R < 1 {\rm pc}$ go non--linear and collapse to form structure much earlier than in standard $\Lambda$CDM cosmology. The resulting mini-halos are hence more dense and resilient to disruption. We therefore expect (provided that the density perturbations on these scales have not been erased by damping processes) that the dark matter distribution on small scales would be more clumpy in chameleon cosmology than in the $\Lambda$CDM model., Comment: 13 pages, 4 figures

The direction dependence of the event rate in WIMP direct detection experiments provides a powerful tool for distinguishing WIMP events from potential backgrounds. We use a variety of (nonparametric) statistical tests to examine the number of events required to distinguish a WIMP signal from an isotropic background when the uncertainty in the reconstruction of the nuclear recoil direction is included in the calculation of the expected signal. We consider a range of models for the Milky Way halo, and also study rotational symmetry tests aimed at detecting nonsphericity/isotropy of the Milky Way halo. Finally we examine ways of detecting tidal streams of WIMPs. We find that if the senses of the recoils are known then of order ten events will be sufficient to distinguish a WIMP signal from an isotropic background for all of the halo models considered, with the uncertainties in reconstructing the recoil direction only mildly increasing the required number of events. If the senses of the recoils are not known the number of events required is an order of magnitude larger, with a large variation between halo models, and the recoil resolution is now an important factor. The rotational symmetry tests require of order a thousand events to distinguish between spherical and significantly triaxial halos, however a deviation of the peak recoil direction from the direction of the solar motion due to a tidal stream could be detected with of order a hundred events, regardless of whether the sense of the recoils is known.