Through a Smoother Lens: An expected absence of LCDM substructure detections from hydrodynamic and dark matter only simulations

A fundamental prediction of the cold dark matter cosmology is the existence of a large number of dark subhalos around galaxies, most of which should be entirely devoid of stars. Confirming the existence of dark substructures stands among the most important empirical challenges in modern cosmology: if they are found and quantified with the mass spectrum expected, then this would close the door on a vast array of competing theories. But in order for observational programs of this kind to reach fruition, we need robust predictions. Here we explore substructure predictions for lensing using galaxy lens-like hosts at z=0.2 from the Illustris simulations both in full hydrodynamics and dark matter only. We quantify substructures more massive than ~ 10^9 M_sun, comparable to current lensing detections derived from HST, Keck, and ALMA. The addition of full hydrodynamics reduces the overall subhalo mass function by about a factor of two. Even for the dark matter only runs, most (~ 85%) lines of sight through projected cylinders of size close to an Einstein radius contain no substructures larger than 10^9 M_sun. The fraction of empty sight lines rises to ~ 95% in full physics simulations. This suggests we will likely need hundreds of strong lensing systems suitable for substructure studies, as well as predictions that include the effects of baryon physics on substructure, to properly constrain cosmological models. Fortunately, the field is poised to fulfill these requirements.Comment: 11 pages, 9 figure

Direct Detection of Dark Matter Debris Flows

Tidal stripping of dark matter from subhalos falling into the Milky Way produces narrow, cold tidal streams as well as more spatially extended "debris flows" in the form of shells, sheets, and plumes. Here we focus on the debris flow in the Via Lactea II simulation, and show that this incompletely phase-mixed material exhibits distinctive high velocity behavior. Unlike tidal streams, which may not necessarily intersect the Earth's location, debris flow is spatially uniform at 8 kpc and thus guaranteed to be present in the dark matter flux incident on direct detection experiments. At Earth-frame speeds greater than 450 km/s, debris flow comprises more than half of the dark matter at the Sun's location, and up to 80% at even higher speeds. Therefore, debris flow is most important for experiments that are particularly sensitive to the high speed tail of the dark matter distribution, such as searches for light or inelastic dark matter or experiments with directional sensitivity. We show that debris flow yields a distinctive recoil energy spectrum and a broadening of the distribution of incidence direction.Comment: 22 pages, 7 figures; accepted for publication in PR

Galactic accretion and the outer structure of galaxies in the CDM model

We have combined the semi-analytic galaxy formation model of Guo et al. (2011) with the particle-tagging technique of Cooper et al. (2010) to predict galaxy surface brightness profiles in a representative sample of ~1900 massive dark matter haloes (10^12--10^14 M_sol) from the Millennium II Lambda-CDM N-body simulation. Here we present our method and basic results focusing on the outer regions of galaxies, consisting of stars accreted in mergers. These simulations cover scales from the stellar haloes of Milky Way-like galaxies to the 'cD envelopes' of groups and clusters, and resolve low surface brightness substructure such as tidal streams. We find that the surface density of accreted stellar mass around the central galaxies of dark matter haloes is well described by a Sersic profile, the radial scale and amplitude of which vary systematically with halo mass (M_200). The total stellar mass surface density profile breaks at the radius where accreted stars start to dominate over stars formed in the galaxy itself. This break disappears with increasing M_200 because accreted stars contribute more of the total mass of galaxies, and is less distinct when the same galaxies are averaged in bins of stellar mass, because of scatter in the relation between M_star and M_200. To test our model we have derived average stellar mass surface density profiles for massive galaxies at z~0.08 by stacking SDSS images. Our model agrees well with these stacked profiles and with other data from the literature, and makes predictions that can be more rigorously tested by future surveys that extend the analysis of the outer structure of galaxies to fainter isophotes. We conclude that it is likely that the outer structure of the spheroidal components of galaxies is largely determined by collisionless merging during their hierarchical assemblyComment: Accepted by MNRAS. Shortened following referee's report, conclusions unchanged. 21 pages, 15 figure

The Mass Profile and Accretion History of Cold Dark Matter Halos

We use the Millennium Simulation series to study the relation between the accretion history (MAH) and mass profile of cold dark matter halos. We find that the mean density within the scale radius, r_{-2} (where the halo density profile has isothermal slope), is directly proportional to the critical density of the Universe at the time when the main progenitor's virial mass equals the mass enclosed within r_{-2}. Scaled to these characteristic values of mass and density, the mean MAH, expressed in terms of the critical density of the Universe, M(\rho_{crit}(z)), resembles that of the enclosed density profile, M(), at z=0. Both follow closely the NFW profile, suggesting that the similarity of halo mass profiles originates from the mass-independence of halo MAHs. Support for this interpretation is provided by outlier halos whose accretion histories deviate from the NFW shape; their mass profiles show correlated deviations from NFW and are better approximated by Einasto profiles. Fitting both M() and M(\rho_{crit}) with either NFW or Einasto profiles yield concentration and shape parameters that are correlated, confirming and extending earlier work linking the concentration of a halo with its accretion history. These correlations also confirm that halo structure is insensitive to initial conditions: only halos whose accretion histories differ greatly from the NFW shape show noticeable deviations from NFW in their mass profiles. As a result, the NFW profile provides acceptable fits to hot dark matter halos, which do not form hierarchically, and for fluctuation power spectra other than CDM. Our findings, however, predict a subtle but systematic dependence of mass profile shape on accretion history which, if confirmed, would provide strong support for the link between accretion history and halo structure we propose here.Comment: 12 pages, 8 figures, MNRAS 432 1103L (2013

Satellite abundances around bright isolated galaxies

We study satellite galaxy abundances in SDSS by counting photometric galaxies around isolated bright primaries. We present results as a function of the luminosity, stellar mass and colour of the satellites, and of the stellar mass and colour of the primaries. For massive primaries the luminosity and stellar mass functions of satellites are similar in shape to those of field galaxies, but for lower mass primaries they are significantly steeper. The steepening is particularly marked for the stellar mass function. Satellite abundance increases strongly with primary stellar mass, approximately in proportion to expected dark halo mass. Massive red primaries have up to a factor of 2 more satellites than blue ones of the same stellar mass. Satellite galaxies are systematically redder than field galaxies of the same stellar mass. Satellites are also systematically redder around more massive primaries. At fixed primary mass, they are redder around red primaries. We select similarly isolated galaxies from mock catalogues based on the simulations of Guo et al.(2011) and analyze them in parallel with the SDSS data. The simulation reproduces all the above trends qualitatively, except for the steepening of the satellite luminosity and stellar mass functions. Model satellites, however, are systematically redder than in the SDSS, particularly at low mass and around low-mass primaries. Simulated haloes of a given mass have satellite abundances that are independent of central galaxy colour, but red centrals tend to have lower stellar masses, reflecting earlier quenching of their star formation by feedback. This explains the correlation between satellite abundance and primary colour in the simulation. The correlation between satellite colour and primary colour arises because red centrals live in haloes which are more massive, older and more gas-rich, so that satellite quenching is more efficient.Comment: 29 pages, 24 figure

A fitting formula for the merger timescale of galaxies in hierarchical clustering

We study galaxy mergers using a high-resolution cosmological hydro/N-body simulation with star formation, and compare the measured merger timescales with theoretical predictions based on the Chandrasekhar formula. In contrast to Navarro et al., our numerical results indicate, that the commonly used equation for the merger timescale given by Lacey and Cole, systematically underestimates the merger timescales for minor mergers and overestimates those for major mergers. This behavior is partly explained by the poor performance of their expression for the Coulomb logarithm, \ln (m_pri/m_sat). The two alternative forms \ln (1+m_pri/m_sat) and 1/2\ln [1+(m_pri/m_sat)^2] for the Coulomb logarithm can account for the mass dependence of merger timescale successfully, but both of them underestimate the merger time scale by a factor 2. Since \ln (1+m_pri/m_sat) represents the mass dependence slightly better we adopt this expression for the Coulomb logarithm. Furthermore, we find that the dependence of the merger timescale on the circularity parameter \epsilon is much weaker than the widely adopted power-law \epsilon^{0.78}, whereas 0.94*{\epsilon}^{0.60}+0.60 provides a good match to the data. Based on these findings, we present an accurate and convenient fitting formula for the merger timescale of galaxies in cold dark matter models.Comment: 16 pages, 14 figures, accepted for publication in ApJ, minor changes in the last few sentences of the discussio

Bound and unbound substructures in Galaxy-scale Dark Matter haloes

We analyse the coarse-grained phase-space structure of the six Galaxy-scale dark matter haloes of the Aquarius Project using a state-of-the-art 6D substructure finder. Within r_50, we find that about 35% of the mass is in identifiable substructures, predominantly tidal streams, but including about 14% in self-bound subhaloes. The slope of the differential substructure mass function is close to -2, which should be compared to around -1.9 for the population of self-bound subhaloes. Near r_50 about 60% of the mass is in substructures, with about 30% in self-bound subhaloes. The inner 35 kpc of the highest resolution simulation has only 0.5% of its mass in self-bound subhaloes, but 3.3% in detected substructure, again primarily tidal streams. The densest tidal streams near the solar position have a 3-D mass density about 1% of the local mean, and populate the high velocity tail of the velocity distribution.Comment: Submitted to MNRAS on 12/10/2010, 11 pages, 10 figure

Red Galaxy Growth and the Halo Occupation Distribution

We have traced the past 7 Gyr of red galaxy stellar mass growth within dark matter halos. We have determined the halo occupation distribution, which describes how galaxies reside within dark matter halos, using the observed luminosity function and clustering of 40,696 0.2<z<1.0 red galaxies in Bootes. Half of 10^{11.9} Msun/h halos host a red central galaxy, and this fraction increases with increasing halo mass. We do not observe any evolution of the relationship between red galaxy stellar mass and host halo mass, although we expect both galaxy stellar masses and halo masses to evolve over cosmic time. We find that the stellar mass contained within the red population has doubled since z=1, with the stellar mass within red satellite galaxies tripling over this redshift range. In cluster mass halos most of the stellar mass resides within satellite galaxies and the intra-cluster light, with a minority of the stellar mass residing within central galaxies. The stellar masses of the most luminous red central galaxies are proportional to halo mass to the power of a third. We thus conclude that halo mergers do not always lead to rapid growth of central galaxies. While very massive halos often double in mass over the past 7 Gyr, the stellar masses of their central galaxies typically grow by only 30%.Comment: Accepted for publication in the ApJ. 34 pages, 22 Figures, 5 Table

Gravitational detection of a low-mass dark satellite at cosmological distance

The mass-function of dwarf satellite galaxies that are observed around Local Group galaxies substantially differs from simulations based on cold dark matter: the simulations predict many more dwarf galaxies than are seen. The Local Group, however, may be anomalous in this regard. A massive dark satellite in an early-type lens galaxy at z = 0.222 was recently found using a new method based on gravitational lensing, suggesting that the mass fraction contained in substructure could be higher than is predicted from simulations. The lack of very low mass detections, however, prohibited any constraint on their mass function. Here we report the presence of a 1.9 +/- 0.1 x 10^8 M_sun dark satellite in the Einstein-ring system JVAS B1938+666 at z = 0.881, where M_sun denotes solar mass. This satellite galaxy has a mass similar to the Sagittarius galaxy, which is a satellite of the Milky Way. We determine the logarithmic slope of the mass function for substructure beyond the local Universe to be alpha = 1.1^+0.6_-0.4, with an average mass-fraction of f = 3.3^+3.6_-1.8 %, by combining data on both of these recently discovered galaxies. Our results are consistent with the predictions from cold dark matter simulations at the 95 per cent confidence level, and therefore agree with the view that galaxies formed hierarchically in a Universe composed of cold dark matter.Comment: 25 pages, 7 figures, accepted for publication in Nature (19 January 2012