Doug Potter

Clumps and streams in the local dark matter distribution

In cold dark matter cosmological models, structures form and grow through the merging of smaller units. Numerical simulations have shown that such merging is incomplete; the inner cores of haloes survive and orbit as ‘subhaloes’ within their hosts. Here we report a simulation that resolves such substructure even in the very inner regions of the Galactic halo. We find hundreds of very concentrated dark matter clumps surviving near the solar circle, as well as numerous cold streams. The simulation also reveals the fractal nature of dark matter clustering: isolated haloes and subhaloes contain the same relative amount of substructure and both have cusped inner density profiles. The inner mass and phase-space densities of subhaloes match those of recently discovered faint, dark-matter-dominated dwarf satellite galaxies, and the overall amount of substructure can explain the anomalous flux ratios seen in strong gravitational lenses. Subhaloes boost γ-ray production from dark matter annihilation by factors of 4 to 15 relative to smooth galactic models. Local cosmic ray production is also enhanced, typically by a factor of 1.4 but by a factor of more than 10 in one per cent of locations lying sufficiently close to a large subhalo. (These estimates assume that the gravitational effects of baryons on dark matter substructure are small.)

Fundamental differences between SPH and grid methods

We have carried out a comparison study of hydrodynamical codes by investigating their performance in modelling interacting multiphase fluids. The two commonly used techniques of grid and smoothed particle hydrodynamics (SPH) show striking differences in their ability to model processes that are fundamentally important across many areas of astrophysics. Whilst Eulerian grid based methods are able to resolve and treat important dynamical instabilities, such as Kelvin-Helmholtz or Rayleigh-Taylor, these processes are poorly or not at all resolved by existing SPH techniques. We show that the reason for this is that SPH, at least in its standard implementation, introduces spurious pressure forces on particles in regions where there are steep density gradients. This results in a boundary gap of the size of an SPH smoothing kernel radius over which interactions are severely damped.

Concentration, spin and shape of dark matter haloes: scatter and the dependence on mass and environment

We use a series of cosmological N-body simulations for a flat A cold dark matter (ACDM) cosmology to investigate the structural properties of dark matter haloes, at redshift zero, in the mass range 3 x 10 9 h -1 ≤ M vir ≤ 3 x 10 13 h -1 M ⊙ . These properties include the concentration parameter, c, the spin parameter, λ, and the mean axis ratio, q. For the concentration-mass relation we find c oc M -0.11 vir agreement with the model proposed by Bullock et al., but inconsistent with the alternative model of Eke et al. The normalization of the concentration-mass relation, however, is 15 per cent lower than suggested by Bullock et al. The results for X and q are in good agreement with previous studies, when extrapolated to the lower halo masses probed here, while c and λ are anticorrelated, in that high-spin haloes have, on average, lower concentrations. In an attempt to remove unrelaxed haloes from the sample, we compute for each halo the offset parameter, x off , defined as the distance between the most bound particle and the centre of mass, in units of the virial radius. Removing haloes with large x off increases the mean concentration by ∼ 10 per cent, lowers the mean spin parameter by ∼ 15 per cent, and removes the most prolate haloes. In addition, it largely removes the anticorrelation between c and λ, though not entirely. We also investigate the relation between halo properties and their large-scale environment density. For low-mass haloes we find that more concentrated haloes live in denser environments than their less concentrated counterparts of the same mass, consistent with recent correlation function analyses. Note, however, that the trend is weak compared to the scatter. For the halo spin parameters we find no environment dependence, while there is a weak indication that the most spherical haloes reside in slightly denser environments. Finally, using a simple model for disc galaxy formation we show that haloes that host low surface brightness galaxies are expected to be hosted by a biased subset of haloes. Not only do these haloes have spin parameters that are larger than average, they also have concentration parameters that are ∼15 per cent lower than the average at a given halo mass. We discuss the implications of all these findings for the claimed disagreement between halo concentrations inferred from low surface brightness rotation curves, and those expected for a ACDM cosmology.

Haloes gone MAD: The Halo-Finder Comparison Project

We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends, spherical-overdensity and phase-space-based algorithms. We

Dark matter direct detection with non-Maxwellian velocity structure

The velocity distribution function of dark matter particles is expected to show significant departures from a Maxwell-Boltzmann distribution. This can have profound effects on the predicted dark matter - nucleon scattering rates in direct detection experiments, especially for dark matter models in which the scattering is sensitive to the high velocity tail of the distribution, such as inelastic dark matter (iDM) or light (few GeV) dark matter (LDM), and for experiments that require high energy recoil events, such as many directionally sensitive experiments. Here we determine the velocity distribution functions from two of the highest resolution numerical simulations of Galactic dark matter structure (Via Lactea II and GHALO), and study the effects for these scenarios. For directional detection, we find that the observed departures from Maxwell-Boltzmann increase the contrast of the signal and change the typical direction of incoming DM particles. For iDM, the expected signals at direct detection experiments are changed dramatically: the annual modulation can be enhanced by more than a factor two, and the relative rates of DAMA compared to CDMS can change by an order of magnitude, while those compared to CRESST can change by a factor of two. The spectrum of the signal can also change dramatically, with many features arising due to substructure. For LDM the spectral effects are smaller, but changes do arise that improve the compatibility with existing experiments. We find that the phase of the modulation can depend upon energy, which would help discriminate against background should it be found.

The graininess of dark matter haloes

We use the recently completed one billion particle Via Lactea II A cold dark matter simulation to investigate local properties like density, mean velocity, velocity dispersion, anisotropy, orientation and shape of the velocity dispersion ellipsoid, as well as the structure in velocity space of dark matter haloes. We show that at the same radial distance from the halo centre, these properties can deviate by orders of magnitude from the canonical, spherically averaged values, a variation that can only be partly explained by triaxiality and the presence of subhaloes. The mass density appears smooth in the central relaxed regions but spans four orders of magnitude in the outskirts, both because of the presence of subhaloes as well as of underdense regions and holes in the matter distribution. In the inner regions, the local velocity dispersion ellipsoid is aligned with the shape ellipsoid of the halo. This is not true in the outer parts where the orientation becomes more isotropic. The clumpy structure in local velocity space of the outer halo cannot be well described by a smooth multivariate normal distribution. Via Lactea II also shows the presence of cold streams made visible by their high 6D phase space density. Generally, the structure of dark matter haloes shows a high degree of graininess in phase space that cannot be described by a smooth distribution function.

Fossil Remnants of Reionization in the Halo of the Milky Way

Our recently completed one billion particle Via Lactea II simulation of a Milky Way-sized dark matter halo resolves over 50,000 gravitationally bound clumps orbiting today within the virialized region of the main host. About 2300 of these subhalos have one or more progenitors with -->M > 106 M? at redshift -->z = 11, i.e., massive enough for their gas to have cooled via excitation of H2 and fragmented prior to the epoch of cosmic reionization. We count 4500 such progenitors: if these were able to convert a fraction of their gas content into very metal-poor stars with a Salpeter initial mass function (IMF), they would be shining today with a visual magnitude -->MV = 6.7 per solar mass in stars. Assuming a universal baryon fraction, we show that mean star formation efficiencies as low as 0.1% in progenitors 108 M? would overproduce the abundance of the faint Galactic dwarf spheroidals observed by the Sloan Digital Sky Survey. Star formation at first light must either have occurred with an IMF lacking stars below 0.9 M?, or was intrinsically very inefficient in small dark matter halos. If the latter, our results may be viewed as another hint that there is a minimum scale in galaxy formation.

Towards an accurate mass function for precision cosmology

Cosmological surveys aim to use the evolution of the abundance of galaxy clusters to accurately constrain the cosmological model. In the context of CDM, we show that it is possible to achieve the required percent level accuracy in the halo mass function with gravity-only cosmological simulations, and we provide simulation start and run parameter guidelines for doing so. Some previous works have had sucient statistical precision, but lacked robust verication of absolute accuracy. Convergence tests of the mass function with, for example, simulation start redshift can exhibit false convergence of the mass function due to counteracting errors, potentially misleading one to infer overly optimistic estimations of simulation accuracy. Percent level accuracy is possible if initial condition particle mapping uses second order Lagrangian Perturbation Theory, and if the start epoch is between 10 and 50 expansion factors before the epoch of halo formation of interest. The mass function for halos with fewer than 1000 particles is highly sensitive to simulation parameters and start redshift, implying a practical minimum mass resolution limit due to mass discreteness. The narrow range in converged start redshift suggests that it is not presently possible for a single simulation to capture accurately the cluster mass function while also starting early enough to model accurately the numbers of reionisation era galaxies, whose baryon feedback processes may aect later cluster properties. Ultimately, to fully exploit current and future cosmological surveys will require accurate modeling of baryon physics and observable properties, a formidable challenge for which accurate gravity-only simulations are just an initial step.

Matter power spectrum and the challenge of percent accuracy

Future galaxy surveys require one percent precision in the theoretical knowledge of the power spectrum over a large range including very nonlinear scales. While this level of accuracy is easily obtained in the linear regime with perturbation theory, it represents a serious challenge for small scales where numerical simulations are required. In this paper we quantify the precision of present-day

Subhaloes gone Notts: spin across subhaloes and finders

N