We present a new method for determining the Galactic gravitational potential based on forward modeling of tidal stellar streams. We use this method to test the performance of smooth and static ...analytic potentials in representing realistic dark matter halos, which have substructure and are continually evolving by accretion. Our FAST-FORWARD method uses a Markov Chain Monte Carlo algorithm to compare, in six-dimensional phase space, an "observed" stream to models created in trial analytic potentials. We analyze a large sample of streams that evolved in the Via Lactea II (VL2) simulation, which represents a realistic Galactic halo potential. The recovered potential parameters are in agreement with the best fit to the global, present-day VL2 potential. However, merely assuming an analytic potential limits the dark matter halo mass measurement to an accuracy of 5%-20%, depending on the choice of analytic parameterization. Collectively, the mass estimates using streams from our sample reach this fundamental limit, but individually they can be highly biased. Individual streams can both under- and overestimate the mass, and the bias is progressively worse for those with smaller perigalacticons, motivating the search for tidal streams at galactocentric distances larger than 70 kpc. We estimate that the assumption of a static and smooth dark matter potential in modeling of the GD-1- and Pal5-like streams introduces an error of up to 50% in the Milky Way mass estimates.
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. Keywords->
We use a particle tagging technique to dynamically populate the N-body Via Lactea II high-resolution simulation with stars. The method is calibrated using the observed luminosity function of Milky ...Way (MW) satellites and the concentration of their stellar populations, and self-consistently follows the accretion and disruption of progenitor dwarfs and the buildup of the stellar halo in a cosmological "live host." Simple prescriptions for assigning stellar populations to collisionless particles are able to reproduce many properties of the observed MW halo and its surviving dwarf satellites, like velocity dispersions, sizes, brightness profiles, metallicities, and spatial distribution. Our model predicts the existence of approximately 1850 subhalos harboring "extremely faint" satellites (with mass-to-light ratios >5 x 10 super(3)) lying beyond the Sloan Digital Sky Survey detection threshold. Of these, about 20 are "first galaxies," i.e., satellites that formed a stellar mass above 10 M sub(middot in circle) before redshift 9. The 10 most luminous satellites (L > 10 super(6) L sub(middot in circle)) in the simulation are hosted by subhalos with peak circular velocities today in the range V sub(max) = 10-40 km s super(-1) that have shed between 80% and 99% of their dark mass after being accreted at redshifts 1.7 < z < 4.6. The satellite maximum circular velocity V sub(max) and stellar line-of-sight velocity dispersion sigma sub(los) today follow the relation V sub(max) = 2.2sigma sub(los). We apply a standard mass estimation algorithm based on Jeans modeling of the line-of-sight velocity dispersion profiles to the simulated dwarf spheroidals and test the accuracy of this technique. The inner (within 300 pc) mass-luminosity relation for currently detectable satellites is nearly flat in our model, in qualitative agreement with the "common mass scale" found in MW dwarfs. We do, however, predict a weak, but significant positive correlation for these objects: M sub(300) is proportional to L super(0.088+ or -0.024).
Bubble chambers and droplet detectors used in dosimetry and dark matter particle search experiments use a superheated metastable liquid in which nuclear recoils trigger bubble nucleation. This ...process is described by the classical heat spike model of F. Seitz Phys. Fluids (1958-1988) 1, 2 (1958)PFLDAS0031-917110.1063/1.1724333, which uses classical nucleation theory to estimate the amount and the localization of the deposited energy required for bubble formation. Here we report on direct molecular dynamics simulations of heat-spike-induced bubble formation. They allow us to test the nanoscale process described in the classical heat spike model. 40 simulations were performed, each containing about 20 million atoms, which interact by a truncated force-shifted Lennard-Jones potential. We find that the energy per length unit needed for bubble nucleation agrees quite well with theoretical predictions, but the allowed spike length and the required total energy are about twice as large as predicted. This could be explained by the rapid energy diffusion measured in the simulation: contrary to the assumption in the classical model, we observe significantly faster heat diffusion than the bubble formation time scale. Finally we examine α-particle tracks, which are much longer than those of neutrons and potential dark matter particles. Empirically, α events were recently found to result in louder acoustic signals than neutron events. This distinction is crucial for the background rejection in dark matter searches. We show that a large number of individual bubbles can form along an α track, which explains the observed larger acoustic amplitudes.
The total mass of dark matter haloes Anderhalden, Donnino; Diemand, Jürg
Monthly notices of the Royal Astronomical Society,
07/2011, Letnik:
414, Številka:
4
Journal Article
Recenzirano
Odprti dostop
The simple, conventional dark matter halo mass definitions commonly used in cosmological simulations ('virial' mass, FoF mass, M
50, 100, 200, ...) only capture part of the collapsed material and are ...therefore inconsistent with the halo mass concept used in analytical treatments of structure formation. Simulations have demonstrated that typical dark matter particle orbits extend out to about 90 per cent of their turnaround radius, which results in apocentre passages outside of the current 'virial' radius on the first and also on the second orbit. Here we describe how the formation history of haloes can be used to identify those particles which took part in the halo collapse, but are missed by conventional group finders because of their remote present location. These particles are added to the part of the halo already identified by FoF. The corrected masses of dark haloes are significantly higher (the median mass increase is 25 per cent) and there is a considerable shift of the halo mass function towards the Press and Schechter form. We conclude that meaningful quantitative comparisons between (semi-)analytic predictions of halo properties (e.g. mass functions, mass accretion rates, merger rates, spatial clustering, etc.) and simulation results will require using the same halo definition in both approaches.
Earth-mass "microhalos" may be the first objects to virialize in the early universe. Their ability to survive the hierarchical clustering process as substructure in the larger halos that form ...subsequently has implications for dark matter detection experiments. We present a large N-body simulation of early substructure in a supersymmetric cold dark matter (SUSY-CDM) scenario characterized by an exponential cutoff in the power spectrum at M sub(c) = 10 super(-6) M sub( ). The simulation resolves a 0.014 M sub( )parent SUSY halo at z = 75 with 14 million particles. On these scales, the effective index of the power spectrum approaches -3, and a range of mass scales collapses almost simultaneously. Compared to a z = 0 galaxy cluster, substructure within our SUSY host is less evident both in phase-space and in physical space, and it is less resistant against tidal disruption. As the universe expands by a factor of 1.3, we find that between 20% and 40% of well-resolved SUSY substructure is destroyed, compared to only 61% in the low-redshift cluster. Nevertheless, SUSY substructure is just as abundant as in z = 0 galaxy clusters; i.e., the normalized mass and circular velocity functions are very similar. The DM self-annihilation g-ray luminosity from bound subhalos and other deviations from a smooth spherical configuration is at least comparable to the spherically averaged signal in the SUSY host, and at least 3 times larger than the spherically averaged signal in the cluster host. Such components must be taken into account when estimating the total cosmological extragalactic g-ray annihilation background. The relative contribution of bound substructure alone to the total annihilation luminosity is about 4 times smaller in the SUSY host than in the z = 0 cluster because of the smaller density contrast of micro-subhalos.
The Milky Way contains several distinct old stellar components that provide a fossil record of its formation. We can understand their spatial distribution and kinematics in a hierarchical formation ...scenario by associating the protogalactic fragments envisaged by Searle & Zinn (1978) with the rare peaks able to cool gas in the cold dark matter density field collapsing at redshift z > 10. We use hierarchical structure formation simulations to explore the kinematics and spatial distribution of these early star-forming structures in galaxy haloes today. Most of the protogalaxies rapidly merge, their stellar contents and dark matter becoming smoothly distributed and forming the inner Galactic halo. The metal-poor globular clusters and old halo stars become tracers of this early evolutionary phase, centrally biased and naturally reproducing the observed steep fall off with radius. The most outlying peaks fall in late and survive to the present day as satellite galaxies. The observed radial velocity dispersion profile and the local radial velocity anisotropy of Milky Way halo stars are successfully reproduced in this model. If this epoch of structure formation coincides with a suppression of further cooling into lower sigma peaks then we can reproduce the rarity, kinematics and spatial distribution of satellite galaxies as suggested by Bullock, Kravtsov & Weinberg (2000). Reionization at z= 12 ± 2 provides a natural solution to the missing satellites problem. Measuring the distribution of globular clusters and halo light on scales from galaxies to clusters could be used to constrain global versus local reionization models. If reionization occurs contemporary, our model predicts a constant frequency of blue globulars relative to the host halo mass, except for dwarf galaxies where the average relative frequencies become smaller.