Lambda CDM Simulations Research Articles

Is The Vast Polar Structure Of Dwarf Galaxies A Serious Problem For ΛCDM?

The dwarf galaxies around the Milky Way are distributed in a so-called vast polar structure (VPOS) that may be in conflict with Lambda CDM simulations. Here, we seek to determine if the VPOS poses a serious challenge to the Lambda cold dark matter paradigm on galactic scales. Specifically, we investigate if the VPOS remains coherent as a function of time. Using the measured Hubble Space Telescope (HST) proper motions and associated uncertainties, we integrate the orbits of the classical Milky Way satellites backwards in time and find that the structure disperses well before a dynamical time. We also examine in particular Leo I and Leo II using their most recent proper motion data, both of which have extreme kinematic properties, but these satellites do not appear to drive the polar fit that is seen at the present day. We have studied the effect of the uncertainties on the HST proper motions on the coherence of the VPOS as a function of time. We find that 8 of the 11 classical dwarfs have reliable proper motions; for these 8, the VPOS also loses significance in less than a dynamical time, indicating that the VPOS is not a dynamically stable structure. Obtaining more accurate proper motion measurements of Ursa Minor, Sculptor, and Carina would bolster these conclusions.

AN INCREASING STELLAR BARYON FRACTION IN BRIGHT GALAXIES AT HIGH REDSHIFT

Recent observations have shown that the characteristic luminosity of the rest-frame ultraviolet (UV) luminosity function does not significantly evolve at 4 < z < 7 and is approximately M*_UV ~ -21. We investigate this apparent non-evolution by examining a sample of 178 bright, M_UV < -21 galaxies at z=4 to 7, analyzing their stellar populations and host halo masses. Including deep Spitzer/IRAC imaging to constrain the rest-frame optical light, we find that M*_UV galaxies at z=4-7 have similar stellar masses of log(M/Msol)=9.6-9.9 and are thus relatively massive for these high redshifts. However, bright galaxies at z=4-7 are less massive and have younger inferred ages than similarly bright galaxies at z=2-3, even though the two populations have similar star formation rates and levels of dust attenuation. We match the abundances of these bright z=4-7 galaxies to halo mass functions from the Bolshoi Lambda-CDM simulation to estimate the halo masses. We find that the typical halo masses in ~M*_UV galaxies decrease from log(M_h/Msol)=11.9 at z=4 to log(M_h/Msol)=11.4 at z=7. Thus, although we are studying galaxies at a similar mass across multiple redshifts, these galaxies live in lower mass halos at higher redshift. The stellar baryon fraction in units of the cosmic mean Omega_b/Omega_m rises from 5.1% at z=4 to 11.7% at z=7; this evolution is significant at the ~3-sigma level. This rise does not agree with simple expectations of how galaxies grow, and implies that some effect, perhaps a diminishing efficiency of feedback, is allowing a higher fraction of available baryons to be converted into stars at high redshifts.

A robust determination of Milky Way satellite properties using hierarchical mass modelling

We introduce a new methodology to robustly determine the mass profile, as well as the overall distribution, of Local Group satellite galaxies. Specifically we employ a statistical multilevel modelling technique, Bayesian hierarchical modelling, to simultaneously constrain the properties of individual Local Group Milky Way satellite galaxies and the characteristics of the Milky Way satellite population. We show that this methodology reduces the uncertainty in individual dwarf galaxy mass measurements up to a factor of a few for the faintest galaxies. We find that the distribution of Milky Way satellites inferred by this analysis, with the exception of the apparent lack of high-mass haloes, is consistent with the Lambda cold dark matter (Lambda-CDM) paradigm. In particular we find that both the measured relationship between the maximum circular velocity and the radius at this velocity, as well as the inferred relationship between the mass within 300 pc and luminosity, match the values predicted by Lambda-CDM simulations for halos with maximum circular velocities below 20 km/sec. Perhaps more striking is that this analysis seems to suggest a more cusped "average" halo shape that is shared by these galaxies. While this study reconciles many of the observed properties of the Milky Way satellite distribution with that of Lambda-CDM simulations, we find that there is still a deficit of satellites with maximum circular velocities of 20-40 km/sec.

HIGH-RESOLUTION MASS MODELS OF DWARF GALAXIES FROM LITTLE THINGS

We present high-resolution rotation curves and mass models of 26 dwarf galaxies from LITTLE THINGS. LITTLE THINGS is a high-resolution Very Large Array HI survey for nearby dwarf galaxies in the local volume within 11 Mpc. The rotation curves of the sample galaxies derived in a homogeneous and consistent manner are combined with Spitzer archival 3.6 micron and ancillary optical U, B, and V images to construct mass models of the galaxies. We decompose the rotation curves in terms of the dynamical contributions by baryons and dark matter halos, and compare the latter with those of dwarf galaxies from THINGS as well as Lambda CDM SPH simulations in which the effect of baryonic feedback processes is included. Being generally consistent with THINGS and simulated dwarf galaxies, most of the LITTLE THINGS sample galaxies show a linear increase of the rotation curve in their inner regions, which gives shallower logarithmic inner slopes alpha of their dark matter density profiles. The mean value of the slopes of the 26 LITTLE THINGS dwarf galaxies is alpha =-0.32 +/- 0.24 which is in accordance with the previous results found for low surface brightness galaxies (alpha = -0.2 +/- 0.2) as well as the seven THINGS dwarf galaxies (alpha =-0.29 +/- 0.07). However, this significantly deviates from the cusp-like dark matter distribution predicted by dark-matter-only Lambda CDM simulations. Instead our results are more in line with the shallower slopes found in the Lambda CDM SPH simulations of dwarf galaxies in which the effect of baryonic feedback processes is included. In addition, we discuss the central dark matter distribution of DDO 210 whose stellar mass is relatively low in our sample to examine the scenario of inefficient supernova feedback in low mass dwarf galaxies predicted from recent Lambda SPH simulations of dwarf galaxies where central cusps still remain.

Evolution of star clusters in a cosmological tidal field

We present a method to couple N-body star cluster simulations to a cosmological tidal field, using the Astrophysical Multipurpose Software Environment. We apply this method to star clusters embedded in the CosmoGrid dark matter-only LambdaCDM simulation. Our star clusters are born at z = 10 (corresponding to an age of the Universe of about 500 Myr) by selecting a dark matter particle and initializing a star cluster with 32,000 stars on its location. We then follow the dynamical evolution of the star cluster within the cosmological environment. We compare the evolution of star clusters in two Milky-Way size haloes with a different accretion history. The mass loss of the star clusters is continuous irrespective of the tidal history of the host halo, but major merger events tend to increase the rate of mass loss. From the selected two dark matter haloes, the halo that experienced the larger number of mergers tends to drive a smaller mass loss rate from the embedded star clusters, even though the final masses of both haloes are similar. We identify two families of star clusters: native clusters, which become part of the main halo before its final major merger event, and the immigrant clusters, which are accreted upon or after this event; native clusters tend to evaporate more quickly than immigrant clusters. Accounting for the evolution of the dark matter halo causes immigrant star clusters to retain more mass than when the z=0 tidal field is taken as a static potential. The reason for this is the weaker tidal field experienced by immigrant star clusters before merging with the larger dark matter halo.

Simulations of Galileon modified gravity: Clustering statistics in real and redshift space

We use N-body simulations to study the statistics of massive halos and redshift space distortions for theories with a standard \Lambda CDM expansion history and a galileon-type scalar field. The extra scalar field increases the gravitational force, leading to enhanced structure formation. We compare our measurements of the real space matter power spectrum and halo properties with fitting formula for estimating these quantities analytically. We find that a model for power spectrum, halo mass-function and halo bias, derived from \Lambda CDM simulations can fit the results from our simulations of modified gravity when \sigma_8 is appropriately adjusted. We also study the redshift space distortions in the two point correlation function measured from these simulations, finding a difference in the ratio of the redshift space to real space clustering amplitude relative to standard gravity on all scales. We find enhanced clustering on scales r>10 Mpc/h and increased damping of the correlation function for scales r<9 Mpc/h. The boost in the clustering on large scales due to the enhanced gravitational forces cannot be mimicked in a standard gravity model by simply changing \sigma_8. This result illustrates the usefulness of redshift space distortion measurements as a probe of modifications to General Relativity.

THE CLUSTERING OF ALFALFA GALAXIES: DEPENDENCE ON H I MASS, RELATIONSHIP WITH OPTICAL SAMPLES, AND CLUES OF HOST HALO PROPERTIES

We use a sample of ~6000 galaxies detected by the Arecibo Legacy Fast ALFA (ALFALFA) 21cm survey, to measure the clustering properties of HI-selected galaxies. We find no convincing evidence for a dependence of clustering on the galactic atomic hydrogen (HI) mass, over the range M_HI ~ 10^{8.5} - 10^{10.5} M_sun. We show that previously reported results of weaker clustering for low-HI mass galaxies are probably due to finite-volume effects. In addition, we compare the clustering of ALFALFA galaxies with optically selected samples drawn from the Sloan Digital Sky Survey (SDSS). We find that HI-selected galaxies cluster more weakly than even relatively optically faint galaxies, when no color selection is applied. Conversely, when SDSS galaxies are split based on their color, we find that the correlation function of blue optical galaxies is practically indistinguishable from that of HI-selected galaxies. At the same time, SDSS galaxies with red colors are found to cluster significantly more than HI-selected galaxies, a fact that is evident in both the projected as well as the full two-dimensional correlation function. A cross-correlation analysis further reveals that gas-rich galaxies "avoid" being located within ~3 Mpc of optical galaxies with red colors. Next, we consider the clustering properties of halo samples selected from the Bolshoi LambdaCDM simulation. A comparison with the clustering of ALFALFA galaxies suggests that galactic HI mass is not tightly related to host halo mass, and that a sizable fraction of subhalos do not host HI galaxies. Lastly, we find that we can recover fairly well the correlation function of HI galaxies by just excluding halos with low spin parameter. This finding lends support to the hypothesis that halo spin plays a key role in determining the gas content of galaxies.

Assembly of filamentary void galaxy configurations

We study the formation and evolution of filamentary configurations of dark matter haloes in voids. Our investigation uses the high-resolution LambdaCDM simulation CosmoGrid to look for void systems resembling the VGS_31 elongated system of three interacting galaxies that was recently discovered by the Void Galaxy Survey (VGS) inside a large void in the SDSS galaxy redshift survey. HI data revealed these galaxies to be embedded in a common elongated envelope, possibly embedded in intravoid filament. In the CosmoGrid simulation we look for systems similar to VGS_31 in mass, size and environment. We find a total of eight such systems. For these systems, we study the distribution of neighbour haloes, the assembly and evolution of the main haloes and the dynamical evolution of the haloes, as well as the evolution of the large-scale structure in which the systems are embedded. The spatial distribution of the haloes follows that of the dark matter environment. We find that VGS_31-like systems have a large variation in formation time, having formed between 10 Gyr ago and the present epoch. However, the environments in which the systems are embedded evolved resemble each other substantially. Each of the VGS_31-like systems is embedded in an intra-void wall, that no later than z = 0.5 became the only prominent feature in its environment. While part of the void walls retain a rather featureless character, we find that around half of them are marked by a pronounced and rapidly evolving substructure. Five haloes find themselves in a tenuous filament of a few Mpc/h long inside the intra-void wall. Finally, we compare the results to observed data from VGS_31. Our study implies that the VGS_31 galaxies formed in the same (proto)filament, and did not meet just recently. The diversity amongst the simulated halo systems indicates that VGS_31 may not be typical for groups of galaxies in voids.

ARGOS – IV. The kinematics of the Milky Way bulge

We present the kinematic results from our ARGOS spectroscopic survey of the Galactic bulge of the Milky Way. Our aim is to understand the formation of the Galactic bulge. We examine the kinematics of about 17,400 stars in the bulge located within 3.5 kpc of the Galactic centre, identified from the 28,000 star ARGOS survey. We aim to determine if the formation of the bulge has been internally driven from disk instabilities as suggested by its boxy shape, or if mergers have played a significant role as expected from Lambda CDM simulations. From our velocity measurements across latitudes b = -5 deg, -7.5 deg and -10 deg we find the bulge to be a cylindrically rotating system that transitions smoothly out into the disk. Within the bulge, we find a kinematically distinct metal-poor population ([Fe/H] < -1.0) that is not rotating cylindrically. The 5% of our stars with [Fe/H] < -1.0 are a slowly rotating spheroidal population, which we believe are stars of the metal weak thick disk and halo which presently lie in the inner Galaxy. The kinematics of the two bulge components that we identified in ARGOS paper III (mean [Fe/H] = -0.25 and [Fe/H] = +0.15, respectively) demonstrate that they are likely to share a common formation origin and are distinct from the more metal poor populations of the thick disk and halo which are colocated inside the bulge. We do not exclude an underlying merger generated bulge component but our results favour bulge formation from instabilities in the early thin disk.

Weak lensing from space: first cosmological constraints from three-point shear statistics★

We use weak lensing data from the Hubble Space Telescope COSMOS survey to measure the second- and third-moments of the cosmic shear field, estimated from about 450,000 galaxies with average redshift <z> ~ 1.3. We measure two- and three-point shear statistics using a tree-code, dividing the signal in E, B and mixed components. We present a detection of the third-order moment of the aperture mass statistic and verify that the measurement is robust against systematic errors caused by point spread function (PSF) residuals and by the intrinsic alignments between galaxies. The amplitude of the measured three-point cosmic shear signal is in very good agreement with the predictions for a WMAP7 best-fit model, whereas the amplitudes of potential systematics are consistent with zero. We make use of three sets of large Lambda CDM simulations to test the accuracy of the cosmological predictions and to estimate the influence of the cosmology-dependent covariance. We perform a likelihood analysis using the measurement and find that the Omega_m-sigma_8 degeneracy direction is well fitted by the relation: sigma_8 (Omega_m/0.30)^(0.49)=0.78+0.11/-0.26. We present the first measurement of a more generalised three-point shear statistic and find a very good agreement with the WMAP7 best-fit cosmology. The cosmological interpretation of this measurement gives sigma_8 (Omega_m/0.30)^(0.46)=0.69 +0.08/-0.14. Furthermore, the combined likelihood analysis of this measurement with the measurement of the second order moment of the aperture mass improves the accuracy of the cosmological constraints, showing the high potential of this combination of measurements to infer cosmological constraints.

Clustering and descendants of MUSYC galaxies at z &lt; 1.5

We measure the evolution of galaxy clustering out to a redshift of z~1.5 using data from two MUSYC fields, the Extended Hubble Deep Field South (EHDF-S) and the Extended Chandra Deep Field South (ECDF-S). We use photometric redshift information to calculate the projected-angular correlation function, omega(sigma), from which we infer the projected correlation function Xi(sigma). We demonstrate that this technique delivers accurate measurements of clustering even when large redshift measurement errors affect the data. To this aim we use two mock MUSYC fields extracted from a LambdaCDM simulation populated with GALFORM semi-analytic galaxies which allow us to assess the degree of accuracy of our estimates of Xi(sigma) and to identify and correct for systematic effects in our measurements. We study the evolution of clustering for volume limited subsamples of galaxies selected using their photometric redshifts and rest-frame r-band absolute magnitudes. We find that the real-space correlation length r_0 of bright galaxies, M_r<-21 (rest-frame) can be accurately recovered out to z~1.5, particularly for ECDF-S given its near-infrared photometric coverage. There is mild evidence for a luminosity dependent clustering in both fields at the low redshift samples (up to <z>=0.57), where the correlation length is higher for brighter galaxies by up to 1Mpc/h between median rest-frame r-band absolute magnitudes of -18 to -21.5. As a result of the photometric redshift measurement, each galaxy is assigned a best-fit template; we restrict to E and E+20%Sbc types to construct subsamples of early type galaxies (ETGs). Our ETG samples show a strong increase in r_0 as the redshift increases, making it unlikely (95% level) that ETGs at median redshift z_med=1.15 are the direct progenitors of ETGs at z_med=0.37 with equivalent passively evolved luminosities. (ABRIDGED)

Offset between dark matter and ordinary matter: evidence from a sample of 38 lensing clusters of galaxies

We compile a sample of 38 galaxy clusters which have both X-ray and strong lensing observations, and study for each cluster the projected offset between the dominant component of baryonic matter center (measured by X-rays) and the gravitational center (measured by strong lensing). Among the total sample, 45% clusters have offsets >10". The >10" separations are significant, considering the arcsecond precision in the measurement of the lensing/X-ray centers. This suggests that it might be a common phenomenon in unrelaxed galaxy clusters that gravitational field is separated spatially from the dominant component of baryonic matter. It also has consequences for lensing models of unrelaxed clusters since the gas mass distribution may differ from the dark matter distribution and give perturbations to the modeling. Such offsets can be used as a statistical tool for comparison with the results of Lambda-CDM simulations and to test the modified dynamics.

Cosmological simulations of normal-branch braneworld gravity

We introduce a cosmological model based on the normal branch of DGP braneworld gravity with a smooth dark energy component on the brane. The expansion history in this model is identical to LambdaCDM, thus evading all geometric constraints on the DGP cross-over scale r_c. This model can serve as a first approximation to more general braneworld models whose cosmological solutions have not been obtained yet. We study the formation of large scale structure in this model in the linear and non-linear regime using N-body simulations for different values of r_c. The simulations use the code presented in (F.S., arXiv:0905.0858) and solve the full non-linear equation for the brane-bending mode in conjunction with the usual gravitational dynamics. The brane-bending mode is attractive rather than repulsive in the DGP normal branch, hence the sign of the modified gravity effects is reversed compared to those presented in arXiv:0905.0858. We compare the simulation results with those of ordinary LambdaCDM simulations run using the same code and initial conditions. We find that the matter power spectrum in this model shows a characteristic enhancement peaking at k ~ 0.7 h/Mpc. We also find that the abundance of massive halos is significantly enhanced. Other results presented here include the density profiles of dark matter halos, and signatures of the brane-bending mode self-interactions (Vainshtein mechanism) in the simulations. Independently of the expansion history, these results can be used to place constraints on the DGP model and future generalizations through their effects on the growth of cosmological structure.

Dark matter halo growth - II. Diffuse accretion and its environmental dependence

Dark matter haloes in Lambda CDM simulations grow by mergers with other haloes as well as accretion of "diffuse" non-halo material. We quantify the mass growth rates via these two processes, dM_mer/dt and dM_dif/dt, and their dependence on halo environment using the ~500,000 haloes in the Millennium simulation. Adopting a local mass density parameter as a measure of halo environment, we find the two rates show strong but opposite environmental dependence, with mergers playing an increasingly important role for halo growths in overdense regions and diffuse accretion dominating growth in voids. This behaviour is independent of the mass cuts used to define haloes vs non-haloes. For galaxy-scale haloes, these two opposite correlations largely cancel out, but a weak environmental dependence remains that results in a slightly lower mean total growth rate, and hence an earlier mean formation redshift, for haloes in denser regions. The mean formation redshift of the ~5000 cluster-mass haloes, on the other hand, appears to have no correlation with halo environment. The origin of the positive correlation of dM_mer/dt with local density can be traced to the surrounding mass reservoir outside the haloes, where more progenitor haloes are available in denser regions. The negative correlation of dM_dif/dt with density, however, is not explained by the available diffuse mass in the reservoir, which is in fact larger in denser regions. The non-halo component may therefore be partially comprised of truly diffuse dark matter particles that are dynamically hotter and are accreted at a suppressed rate in denser regions. We also discuss the implications of these results for how to modify the Extended Press-Schechter model of halo growth, which in its original form does not predict environmental dependence.

The graininess of dark matter haloes

We use the recently completed one billion particle Via Lactea II LambdaCDM simulation to investigate local properties like density, mean velocity, velocity dispersion, anisotropy, orientation and shape of the velocity dispersion ellipsoid, as well as 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 can not 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.

Discreteness Effects in ΛCDM Simulations: A Wavelet‐Statistical View

The effects of particle discreteness in N-body LambdaCDM simulations are still an intensively debated issue. In this paper we explore such effects, taking into account the scatter caused by the randomness of the initial conditions and focusing on the statistical properties of the cosmological density field. For this purpose, we run large sets of LambdaCDM simulations and analyze them using a wide variety of diagnostics, including new and powerful wavelet statistics. Among other facts, we point out (1) that dynamical evolution does not propagate discreteness noise up from the small scales at which it is introduced and (2) that one should aim to satisfy the condition e~2d, where e is the force resolution and d is the interparticle distance. We clarify what such a condition means and how to implement it in modern cosmological codes.

Two‐dimensional Topology of the Sloan Digital Sky Survey

We present the topology of a volume-limited sample of 11,884 galaxies, selected from an apparent-magnitude limited sample of over 100,000 galaxies observed as part of the Sloan Digital Sky Survey (SDSS). The data currently cover three main regions on the sky: one in the Galactic north and one in the south, both at zero degrees declination, and one area in the north at higher declination. Each of these areas covers a wide range of survey longitude but a narrow range of survey latitude, allowing the two dimensional genus to be measured. The genus curves of the SDSS sub-samples are similar, after appropriately normalizing these measurements for the different areas. We sum the genus curves from the three areas to obtain the total genus curve of the SDSS. The total curve has a shape similar to the genus curve derived from mock catalogs drawn from the Hubble Volume Lambda CDM simulation and is similar to that of a Gaussian random field. Likewise, comparison with the genus of the 2dFGRS, after normalization for the difference in area, reveals remarkable similarity in the topology of these samples. We test for the effects of galaxy type segregation by splitting the SDSS data into thirds, based on the u^*- r^* colors of the galaxies, and measure the genus of the reddest and bluest sub-samples. This red/blue split in u^*- r^* is essentially a split by morphology (Strateva et al. 2001). We find that the genus curve for the reddest galaxies exhibits a ``meatball'' shift of the topology -- reflecting the concentration of red galaxies in high density regions -- compared to the bluest galaxies and the full sample, in agreement with predictions from simulations.

Galaxy Halo Masses from Galaxy‐Galaxy Lensing

We present measurements of the extended dark halo profiles of bright early-type galaxies at redshifts 0.1 < z < 0.9 obtained via galaxy-galaxy lensing analysis of images taken at the Canada-France-Hawaii Telescope using the UH8K CCD mosaic camera. Six 05 × 05 fields were observed for a total of 2 hr each in I and V, resulting in catalogs containing ~20,000 galaxies per field. We used V-I color and I magnitude to select bright early-type galaxies as the lens galaxies, yielding a sample of massive lenses with fairly well-determined redshifts and absolute magnitudes M ~ M* ± 1. We paired these with faint galaxies lying at angular distances 20'' < θ < 60'', corresponding to physical radii of 26 < r < 77 h-1 kpc (z = 0.1) and 105 < r < 315 h-1 kpc (z = 0.9), and computed the mean tangential shear γT(θ) of the faint galaxies. The shear falls off with radius roughly as γT ∝ 1/θ as expected for flat rotation curve halos. The shear values were weighted in proportion to the square root of the luminosity of the lens galaxy. This is optimal if the halo mass at a given radius varies as M ∝ , as is the case at smaller radii, and in this context our results give a value for the average mean rotation velocity of an L* galaxy halo at r ~ 50-200 h-1 kpc of v* = 238 km s-1 for a flat lambda (Ωm0 = 0.3,Ωλ0 = 0.7) cosmology (v* = 269 km s-1 for Einstein-de Sitter) and with little evidence for evolution with redshift. These halo masses are somewhat (2-3 times) lower than a simple perfectly flat rotation curve extrapolation from smaller scale dynamical measurements. They are also considerably lower than the masses of halos found from the best-studied X-ray halos. They do, however, agree extremely well with the masses of halos of the same abundance in lambda-CDM simulations. We find a mass-to-light ratio of M/LB ≃ 121 ± 28h(r/100 h-1 kpc) (for L* galaxies) and these halos constitute Ω ≃ 0.04 ± 0.01(r/100 h-1 kpc) of closure density.