Cold Dark Matter Simulations Research Articles

Weighing the Local Group in the presence of dark energy

Abstract We revise the mass estimate of the Local Group (LG) when dark energy (in the form of the cosmological constant) is incorporated into the timing argument (TA) mass estimator for the LG. Assuming the age of the Universe and the cosmological constant according to the recent values from the Planck cosmic microwave background experiment, we find the mass of the LG to be MTAΛ = (4.73 ± 1.03) × 1012 M⊙ which is 13 per cent higher than the classical TA mass estimate. This partly explains the discrepancy between earlier results from Lambda cold dark matter simulations and the classical TA. When a similar analysis is performed on 16 LG-like galaxy pairs from the Constrained Local Universe Simulations, we find that the scatter in the ratio of the virial to the TA estimated mass is given by Mvir/MTAΛ = 1.04 ± 0.16. Applying it to the LG mass estimation we find a calibrated Mvir = (4.92 ± 1.08(obs.) ± 0.79(sys.)) × 1012 M⊙.

The warm dark matter halo mass function below the cut-off scale

Warm Dark Matter (WDM) cosmologies are a viable alternative to the Cold Dark Matter (CDM) scenario. Unfortunately, an accurate scrutiny of the WDM predictions with N-body simulations has proven difficult due to numerical artefacts. Here, we report on cosmological simulations that, for the first time, are devoid of those problems, and thus, are able to accurately resolve the WDM halo mass function well below the cut-off. We discover a complex picture, with perturbations at different evolutionary stages populating different ranges in the halo mass function. On the smallest mass scales we can resolve, identified objects are typically centres of filaments that are starting to collapse. On intermediate mass scales, objects typically correspond to fluctuations that have collapsed and are in the process of relaxation, whereas the high mass end is dominated by objects similar to haloes identified in CDM simulations. We then explicitly show how the formation of low-mass haloes is suppressed, which translates into a strong cut-off in the halo mass function. This disfavours some analytic formulations that predict a halo mass function that would extend well below the free streaming mass. We argue for a more detailed exploration of the formation of the smallest structures expected to form in a given cosmology, which, we foresee, will advance our overall understanding of structure formation.

LUMINOUS SATELLITES VERSUS DARK SUBHALOS: CLUSTERING IN THE MILKY WAY

The observed population of the Milky Way satellite galaxies offer a unique testing ground for galaxy formation theory on small-scales. Our novel approach was to investigate the clustering of the known Milky Way satellite galaxies and to quantify the amount of substructure within their distribution using a two-point correlation function statistic in each of three spaces: configuration space, line-of-sight velocity space, and four-dimensional phase-space. These results were compared to those for three sets of subhaloes in the Via Lactea II Cold Dark Matter simulation defined to represent the luminous dwarfs. We found no evidence at a significance level above 2-sigma of substructure within the distribution of the Milky Way satellite galaxies in any of the three spaces. The "luminous" subhalo sets are more strongly clustered than are the Milky Way satellites in all three spaces and over a broader range of scales in four-dimensional phase-space. Each of the "luminous" subhalo sets are clustered as a result of substructure within their line-of-sight velocity space distributions at greater than 3-sigma significance, whereas the Milky Way satellite galaxies are randomly distributed in line-of-sight velocity space. While our comparison is with only one Cold Dark Matter simulation, the inconsistencies between the Milky Way satellite galaxies and the Via Lactea II subhalo sets for all clustering methods suggest a potential new 'small-scale' tension between Cold Dark Matter theory and the observed Milky Way satellites. Future work will obtain a more robust comparison between the observed Milky Way satellites and Cold Dark Matter theory by studying additional simulations.

Indirect detection of self-interacting asymmetric dark matter

Self-interacting dark matter resolves the issue of cuspy profiles that appear in noninteracting cold dark matter simulations; it may additionally resolve the so-called ``too big to fail'' problem in structure formation. Asymmetric dark matter provides a natural explanation of the comparable densities of baryonic matter and dark matter. In this paper, we discuss unique indirect detection signals produced by a minimal model of self-interacting asymmetric scalar dark matter. Through the formation of dark matter bound states, a dark force mediator particle may be emitted; the decay of this particle may produce an observable signal. We estimate the produced signal and explicitly demonstrate parameters for which the signal exceeds current observations.

A BARYONIC SOLUTION TO THE MISSING SATELLITES PROBLEM

It has been demonstrated that the inclusion of baryonic physics can alter the dark matter densities in the centers of low-mass galaxies, making the central dark matter slope more shallow than predicted in pure cold dark matter simulations. This flattening of the dark matter profile can occur in the most luminous subhalos around Milky Way-mass galaxies. Zolotov et al. (2012) have suggested a correction to be applied to the central masses of dark matter-only satellites in order to mimic the affect of (1) the flattening of the dark matter cusp due to supernova feedback in luminous satellites, and (2) enhanced tidal stripping due to the presence of a baryonic disk. In this paper, we apply this correction to the z=0 subhalo masses from the high resolution, dark matter-only Via Lactea II (VL2) simulation, and find that the number of massive subhalos is dramatically reduced. After adopting a stellar mass to halo mass relationship for the VL2 halos, and identifying subhalos that are (1) likely to be destroyed by stripping and (2) likely to have star formation suppressed by photo-heating, we find that the number of massive, luminous satellites around a Milky Way-mass galaxy is in agreement with the number of observed satellites around the Milky Way or M31. We conclude that baryonic processes have the potential to solve the missing satellites problem.

MEASURING DARK MATTER PROFILES NON-PARAMETRICALLY IN DWARF SPHEROIDALS: AN APPLICATION TO DRACO

We introduce a novel implementation of orbit-based (or Schwarzschild) modeling that allows dark matter density profiles to be calculated non-parametrically in nearby galaxies. Our models require no assumptions to be made about velocity anisotropy or the dark matter profile. The technique can be applied to any dispersion-supported stellar system, and we demonstrate its use by studying the Local Group dwarf spheroidal (dSph) galaxy Draco. We use existing kinematic data at larger radii and also present 12 new radial velocities within the central 13 pc obtained with the VIRUS-W integral field spectrograph on the 2.7m telescope at McDonald Observatory. Our non-parametric Schwarzschild models find strong evidence that the dark matter profile in Draco is cuspy for 20 < r < 700 pc. The profile for r > 20 pc is well-fit by a power law with slope \alpha=-1.0 +/- 0.2, consistent with predictions from Cold Dark Matter (CDM) simulations. Our models confirm that, despite its low baryon content relative to other dSphs, Draco lives in a massive halo.

Not too big, not too small: the dark haloes of the dwarf spheroidals in the Milky Way

We present a new analysis of the Aquarius simulations done in combination with a semi-analytic galaxy formation model. Our goal is to establish whether the subhaloes present in Λ cold dark matter simulations of Milky Way (MW) like systems could host the dwarf spheroidal (dSph) satellites of our Galaxy. Our analysis shows that, contrary to what has been assumed in most previous work, the mass profiles of subhaloes are generally not well fitted by Navarro–Frenk–White models but that Einasto profiles are preferred. We find that for shape parameters α = 0.2–0.5 and vmax = 10–30 km s−1 there is very good correspondence with the observational constraints obtained for the nine brightest dSphs of the MW. However, to explain the internal dynamics of these systems as well as the number of objects of a given circular velocity the total mass of the MW should be ∼8 × 1011 M⊙, a value that is in agreement with many recent determinations, and at the low-mass end of the range explored by the Aquarius simulations. Our simulations show important scatter in the number of bright satellites, even when the Aquarius MW-like hosts are scaled to a common mass, and we find no evidence for a missing population of massive subhaloes in the Galaxy. This conclusion is also supported when we examine the dynamics of the satellites of M31.

Stochastic star formation and early galactic nucleosynthesis

AbstractWe discuss calculations of the star formation, nucleosynthesis, and stochastic evolution of proto-galactic clouds in a galactic chemical evolution model which is motivated by cold dark matter simulations of hierarchical galaxy formation (Saleh et al. 2006; Lan et al. 2010). We utilize SN-induced and dark matter halo formation-induced star formation within a model that follows the evolution of chemical enrichment and energy input to the clouds via Type II, Ia supernovae and stellar winds.

CLASH: MASS DISTRIBUTION IN AND AROUND MACS J1206.2-0847 FROM A FULL CLUSTER LENSING ANALYSIS

We derive an accurate mass distribution of the galaxy cluster MACS J1206.2-0847 (z=0.439) from a combined weak-lensing distortion, magnification, and strong-lensing analysis of wide-field Subaru BVRIz' imaging and our recent 16-band Hubble Space Telescope observations taken as part of the Cluster Lensing And Supernova survey with Hubble (CLASH) program. We find good agreement in the regions of overlap between several weak and strong lensing mass reconstructions using a wide variety of modeling methods, ensuring consistency. The Subaru data reveal the presence of a surrounding large scale structure with the major axis running approximately north-west south-east (NW-SE), aligned with the cluster and its brightest galaxy shapes, showing elongation with a \sim 2:1 axis ratio in the plane of the sky. Our full-lensing mass profile exhibits a shallow profile slope dln\Sigma/dlnR\sim -1 at cluster outskirts (R>1Mpc/h), whereas the mass distribution excluding the NW-SE excess regions steepens further out, well described by the Navarro-Frenk-White form. Assuming a spherical halo, we obtain a virial mass M_{vir}=(1.1\pm 0.2\pm 0.1)\times 10^{15} M_{sun}/h and a halo concentration c_{vir} = 6.9\pm 1.0\pm 1.2 (\sim 5.7 when the central 50kpc/h is excluded), which falls in the range 4< <c> <7 of average c(M,z) predictions for relaxed clusters from recent Lambda cold dark matter simulations. Our full lensing results are found to be in agreement with X-ray mass measurements where the data overlap, and when combined with Chandra gas mass measurements, yield a cumulative gas mass fraction of 13.7^{+4.5}_{-3.0}% at 0.7Mpc/h (\approx 1.7r_{2500}), a typical value observed for high mass clusters.

Subhaloes in self-interacting galactic dark matter haloes

We present N-body simulations of a new class of self-interacting dark matter models, which do not violate any astrophysical constraints due to a non-power-law velocity dependence of the transfer cross section which is motivated by a Yukawa-like new gauge boson interaction. Specifically, we focus on the formation of a Milky Way-like dark matter halo taken from the Aquarius project and re-simulate it for a couple of representative cases in the allowed parameter space of this new model. We find that for these cases, the main halo only develops a small core (~1 kpc) followed by a density profile identical to that of the standard cold dark matter scenario outside of that radius. Neither the subhalo mass function nor the radial number density of subhaloes are altered in these models but there is a significant change in the inner density structure of subhaloes resulting in the formation of a large density core. As a consequence, the inner circular velocity profiles of the most massive subhaloes differ significantly from the cold dark matter predictions and we demonstrate that they are compatible with the observational data of the brightest Milky Way dSphs in such a velocity-dependent self-interacting dark matter scenario. Specifically, and contrary to the cold dark matter case, there are no subhaloes that are more concentrated than what is inferred from the kinematics of the Milky Way dSphs. We conclude that these models offer an interesting alternative to the cold dark matter model that can reduce the recently reported tension between the brightest Milky Way satellites and the dense subhaloes found in cold dark matter simulations.

THE COSMOLOGICAL IMPACT OF LUMINOUS TeV BLAZARS. III. IMPLICATIONS FOR GALAXY CLUSTERS AND THE FORMATION OF DWARF GALAXIES

A subset of blazars emit TeV gamma rays which annihilate and pair produce on the extragalactic background light. We have argued in Broderick et al. (2011, Paper I) that plasma beam instabilities can dissipate the pairs' energy locally. This heats the intergalactic medium and dramatically increases its entropy after redshift z~2, with important implications for structure formation: (1) This suggests a scenario for the origin of the cool core (CC)/non-cool core (NCC) bimodality in galaxy clusters and groups. Early forming galaxy groups are unaffected because they can efficiently radiate the additional entropy, developing a CC. However, late forming groups do not have sufficient time to cool before the entropy is gravitationally reprocessed through successive mergers - counteracting cooling and raising the core entropy further. Hence blazar heating works different than feedback by active galactic nuclei, which balances radiative cooling but is unable to transform CC into NCC clusters due to the weak coupling to the cluster gas. (2) We predict a suppression of the Sunyaev-Zel'dovich power spectrum on angular scales smaller than 5' due to the globally reduced central pressure of groups and clusters forming after z~1. (3) Our redshift dependent entropy floor increases the characteristic halo mass below which dwarf galaxies cannot form by a factor of ~10 (50) at mean density (in voids) over that found in models that include photoionization alone. This prevents the formation of late forming dwarf galaxies (z<2) with masses ranging from 10^{10} to 10^{11} M_sun for redshifts z~2 to 0, respectively. This may help resolve the "missing satellite problem" in the Milky Way and the "void phenomenon" of the low observed abundances of dwarf satellites compared to cold dark matter simulations and may bring the observed early star formation histories into agreement with galaxy formation models. (abridged)

THE DARK MATTER DENSITY PROFILE OF THE FORNAX DWARF

We construct axisymmetric Schwarzschild models to measure the mass profile of the local group dwarf galaxy Fornax. These models require no assumptions to be made about the orbital anisotropy of the stars, as is the case for commonly used Jeans models. We test a variety of parameterizations of dark matter density profiles and find cored models with uniform density rho_c = (1.6 +/- 0.1) x 10^-2 M_sun pc^-3 fit significantly better than the cuspy halos predicted by cold dark matter simulations. We also construct models with an intermediate-mass black hole, but are unable to make a detection. We place a 1-sigma upper limit on the mass of a potential intermediate-mass black hole at M_BH < 3.2 x 10^4 M_sun.

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

The mass function of dwarf satellite galaxies that are observed around Local Group galaxies differs substantially 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 a redshift of 0.222 was recently found using a 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) × 10(8) M dark satellite galaxy in the Einstein ring system JVAS B1938+666 (ref. 11) at a redshift of 0.881, where M denotes the solar mass. This satellite galaxy has a mass similar to that of 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 1.1(+0.6)(-0.4), with an average mass fraction of 3.3(+3.6)(-1.8) per cent, 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.

WHERE ARE THE FOSSILS OF THE FIRST GALAXIES? II. TRUE FOSSILS, GHOST HALOS, AND THE MISSING BRIGHT SATELLITES

We use a new set of cold dark matter simulations of the local universe to investigate the distribution of fossils of primordial dwarf galaxies within, and around the Milky Way. Throughout, we build upon previous results showing agreement between the observed stellar properties of a subset of the ultra-faint dwarfs and our simulated fossils. Here, we show that fossils of the first galaxies have galactocentric distributions and cumulative luminosity functions consistent with observations. In our model there are ~ 300 luminous satellites orbiting the Milky Way, ~50-70% of which are well preserved fossils, with this fraction decreasing with galactocentric distance. Within the Milky Way virial radius, the majority of these fossils have luminosities L_V<10^5 L_solar. This work produces an overabundance of bright dwarf satellites (L_V > 10^4 L_solar) with respect to observations where observations are nearly complete. The "bright satellite problem" is most evident in the outer parts of the Milky Way. We estimate that, although relatively bright, the primordial stellar populations are very diffuse, producing a population with surface brightnesses below surveys` detection limits and are easily stripped by tidal forces. Although we cannot yet present unmistakable evidence for the existence of the fossils of first galaxies in the Local Group, the results of our studies suggest observational strategies that may demonstrate their existence. Primarily, the detection of "ghost halos" of primordial stars around isolated dwarfs would prove that stars formed in minihalos (M<10^8 M_solar) before reionization, and strongly suggest that at least a fraction of the ultra-faint dwarfs are fossils of the first galaxies.

THE VELOCITY WIDTH FUNCTION OF GALAXIES FROM THE 40% ALFALFA SURVEY: SHEDDING LIGHT ON THE COLD DARK MATTER OVERABUNDANCE PROBLEM

The ongoing Arecibo Legacy Fast ALFA (ALFALFA) survey is a wide-area, extragalactic HI-line survey conducted at the Arecibo Observatory. Sources have so far been extracted over ~3,000 sq.deg of sky (40% of its final area), resulting in the largest HI-selected sample to date. We measure the space density of HI-bearing galaxies as a function of their observed velocity width (uncorrected for inclination) down to w = 20 km/s, a factor of 2 lower than the previous generation HI Parkes All-Sky Survey. We confirm previous results that indicate a substantial discrepancy between the observational distribution and the theoretical one expected in a cold dark matter (CDM) universe, at low widths. In particular, a comparison with synthetic galaxy samples populating state-of-the-art CDM simulations imply a factor of ~8 difference in the abundance of galaxies with w = 50 km/s (increasing to a factor of ~100 when extrapolated to the ALFALFA limit of w = 20 km/s). We furthermore identify possible solutions, including a keV warm dark matter scenario and the fact that HI disks in low mass galaxies are usually not extended enough to probe the full amplitude of the galactic rotation curve. In this latter case, we can statistically infer the relationship between the measured HI rotational velocity of a galaxy and the mass of its host CDM halo. Observational verification of the presented relationship at low velocities would provide an important test of the validity of the established dark matter model.

Too big to fail? The puzzling darkness of massive Milky Way subhaloes

Abstract We show that dissipationless Λ cold dark matter simulations predict that the majority of the most massive subhaloes of the Milky Way are too dense to host any of its bright satellites (LV &amp;gt; 105 L⊙). These dark subhaloes have peak circular velocities at infall of Vinfall= 30–70 km s-1 and infall masses of (0.2–4) × 1010 M⊙. Unless the Milky Way is a statistical anomaly, this implies that galaxy formation becomes effectively stochastic at these masses. This is in marked contrast to the well-established monotonic relation between galaxy luminosity and halo circular velocity (or halo mass) for more massive haloes. We show that at least two (and typically four) of these massive dark subhaloes are expected to produce a larger dark matter annihilation flux than Draco. It may be possible to circumvent these conclusions if baryonic feedback in dwarf satellites or different dark matter physics can reduce the central densities of massive subhaloes by order unity on a scale of 0.3–1 kpc.

DARK AND LUMINOUS MATTER IN THINGS DWARF GALAXIES

We present mass models for the dark matter component of seven dwarf galaxies taken from "The HI Nearby Galaxy Survey" (THINGS) and compare these with those from numerical Lambda Cold Dark Matter (LCDM) simulations. The THINGS high-resolution data significantly reduce observational uncertainties and thus allow us to derive accurate dark matter distributions in these systems. We here use the bulk velocity fields when deriving the rotation curves of the galaxies. Compared to other types of velocity fields, the bulk velocity field minimizes the effect of small-scale random motions more effectively and traces the underlying kinematics of a galaxy more properly. The "Spitzer Infrared Nearby Galaxies Survey" (SINGS) 3.6 micron and ancillary optical data are used for separating the baryons from their total matter content in the galaxies. The sample dwarf galaxies are found to be dark matter dominated over most radii. We find discrepancies between the derived dark matter distributions of the galaxies and those of LCDM simulations, even after corrections for non-circular motions have been applied. The observed solid body-like rotation curves of the galaxies rise too slowly to reflect the cusp-like dark matter distribution in CDM halos. Instead, they are better described by core-like models such as pseudo-isothermal halo models dominated by a central constant-density core. The mean value of the logarithmic inner slopes of the mass density profiles is alpha = -0.29 +- 0.07. They are significantly different from the steep slope of ~ -1.0 inferred from previous dark-matter-only simulations, and are more consistent with shallower slopes found in recent LCDM simulations of dwarf galaxies in which the effects of baryonic feedback processes are included.

The diversity and similarity of simulated cold dark matter haloes

We study the mass, velocity dispersion, and anisotropy profiles of $\Lambda$CDM halos using a suite of N-body simulations of unprecedented numerical resolution (the {\it Aquarius Project}). Our analysis confirms a number of results claimed by earlier work, and clarifies a few issues where conflicting claims may be found in the recent literature. The spherically-averaged density profile becomes progressively shallower inwards and, at the innermost resolved radius, the logarithmic slope is $\gamma \equiv -$d$\ln\rho/$d$\ln r \simlt 1$. Asymptotic inner slopes as steep as the recently claimed $\rho \propto r^{-1.2}$ are clearly ruled out. The radial dependence of $\gamma$ is well approximated by a power-law, $\gamma \propto r^{\alpha}$ (the Einasto profile). The shape parameter, $\alpha$, varies slightly but significantly from halo to halo, implying that the mass profiles of $\Lambda$CDM halos are not strictly universal: different halos cannot, in general, be rescaled to look identical. Departures from similarity are also seen in velocity dispersion profiles and correlate with those in density profiles so as to preserve a power-law form for the spherically averaged pseudo-phase-space density, $\rho/\sigma^3\propto r^{-1.875}$. Our conclusions are reliable down to radii below 0.4% of the virial radius, providing well-defined predictions for halo structure when baryonic effects are neglected, and thus an instructive theoretical template against which the modifications induced by the baryonic components of real galaxies can be judged.

Simulations of quintessential cold dark matter: beyond the cosmological constant

We study the non-linear growth of cosmic structure in different dark energy models, using large volume N-body simulations. We consider a range of quintessence models which feature both rapidly and slowly varying dark energy equations of state, and compare the growth of structure to that in a universe with a cosmological constant. We use a four-parameter equation of state for the dark energy which accurately reproduces the quintessence dynamics over a wide range of redshifts. The adoption of a quintessence model changes the expansion history of the universe, the form of the linear theory power spectrum and can alter key observables, such as the horizon scale and the distance to last scattering. We incorporate these effects into our simulations in stages to isolate the impact of each on the growth of structure. The difference in structure formation can be explained to first order by the difference in growth factor at a given epoch; this scaling also accounts for the non-linear growth at the 15 per cent level. We find that quintessence models that are different from Λ cold dark matter (ΛCDM) both today and at high redshifts (z∼ 1000), and which feature late (z < 2), rapid transitions in the equation of state, can have identical baryonic acoustic oscillation (BAO) peak positions to those in ΛCDM. We find that these models have higher abundances of dark matter haloes at z > 0 compared to ΛCDM and so measurements of the mass function should allow us to distinguish these quintessence models from a cosmological constant. However, we find that a second class of quintessence models, whose equation of state makes an early (z > 2) rapid transition to w=−1, cannot be distinguished from ΛCDM using measurements of the mass function or the BAO, even if these models have non-negligible amounts of dark energy at early times.

Non-circular motions and the cusp-core discrepancy in dwarf galaxies

Context. The cusp-core discrepancy is one of the major problems in astrophysics. It results from comparing the observed mass distribution of galaxies with the predictions of Cold Dark Matter simulations. The latter predict a cuspy density profile in the inner parts of galaxies, whereas observations of dwarf and low surface brightness galaxies show a constant density core. Aims. We want to determine the shape of the dark matter potential in the nuclear regions of a sample of six nearby irregular dwarf galaxies. Methods. In order to quantify the amount of non-circular motions which could potentially affect a mass decomposition, we first perform a harmonic decomposition of the HI Hermite velocity fields of all sample galaxies. We then decompose the HI rotation curves into different mass components by fitting NFW and pseudo-isothermal halo models to the HI rotation curves using a chi^2 minimisation. We model the minimum-disc, the minimum-disc+gas, and the maximum-disc cases. Results. The non-circular motions are in all cases studied here of the order of only a few km/s (generally corresponding to less than 25% of the local rotation velocity), which means that they do not significantly affect the rotation curves. The observed rotation curves can better be described by the cored pseudo-isothermal halo than by the NFW halo. The slopes of the dark matter density profiles confirm this result and are in good agreement with previous studies. The quality of the fits can often be improved when including the baryons, which suggests that they contribute significantly to the inner part of the density profile of dwarf galaxies.