Voids in a ΛCDM universe

We study the formation and evolution of voids in the dark matter distribution using various simulations of the popular Λ cold dark matter cosmogony. We identify voids by requiring them to be regions of space with a mean overdensity of −0.8 or less – roughly the equivalent of using a spherical overdensity group finder for haloes. Each of the simulations contains thousands of voids. The distribution of void sizes in the different simulations shows good agreement when differences in particle and grid resolution are accounted for. Voids very clearly correspond to minima in the smoothed initial density field. Apart from a very weak dependence on the mass resolution, the rescaled mass profiles of voids in the different simulations agree remarkably well. We find a universal void mass profile of the form  ρ(< r )/ρ( r _eff) ∝ exp[( r / r _eff)^α]  , where r _eff is the effective radius of a void and  α∼ 2  . The mass function of haloes in voids is steeper than that of haloes that populate denser regions. In addition, the abundances of void haloes seem to evolve somewhat more strongly between redshifts ∼1 and 0 than the global abundances of haloes.     

Environmental dependence of dark matter halo growth – I. Halo merger rates

In an earlier paper, we quantified the mean merger rate of dark matter haloes as a function of redshift z , descendant halo mass M ₀, and progenitor halo mass ratio ξ using the Millennium simulation of the Λ cold dark matter cosmology. Here, we broaden that study and investigate the dependence of the merger rate of haloes on their surrounding environment. A number of local mass overdensity variables, both including and excluding the halo mass itself, are tested as measures of a halo's environment. The simple functional dependence on   z , M ₀  , and ξ of the merger rate found in our earlier work, is largely preserved in different environments, but we find that the overall amplitude of the merger rate has a strong positive correlation with the environmental densities. For galaxy-mass haloes, we find mergers to occur ∼2.5 times more frequently in the densest regions than in voids at both   z = 0  and higher redshifts. Higher mass haloes show similar trends. We present a fitting form for this environmental dependence that is a function of both mass and local density and valid out to   z = 2  . The amplitude of the progenitor (or conditional) mass function shows a similar correlation with local overdensity, suggesting that the extended Press–Schechter model for halo growth needs to be modified to incorporate environmental effects.     

Testing tidal-torque theory – I. Spin amplitude and direction

We evaluate the success of linear tidal-torque theory (TTT) in predicting galactic-halo spin using a cosmological N -body simulation with thousands of well-resolved haloes. The protohaloes are identified by tracing today's haloes back to the initial conditions. The TTT predictions for the protohaloes match, on average, the spin amplitudes of the virialized haloes of today, if linear growth is assumed until ∼ t ₀/3, or  55–70  per cent of the halo effective turn-around time. This makes it a useful qualitative tool for understanding certain average properties of galaxies, such as total spin and angular momentum distribution within haloes, but with a random scatter of the order of the signal itself. Non-linear changes in spin direction cause a mean error of ∼50° in the TTT prediction at t ₀, such that the linear spatial correlations of spins on scales ≥1  h ⁻¹ Mpc are significantly weakened by non-linear effects. This questions the usefulness of TTT for predicting intrinsic alignments in the context of gravitational lensing. We find that the standard approximations made in TTT, including a second-order expansion of the Zel'dovich potential and a smoothing of the tidal field, provide close-to-optimal results.     

How big were the first cosmological objects?

We calculate the cooling times at constant density for haloes with virial temperatures from 100 K to  1×10⁵ K  that originate from a 3 σ fluctuation of a CDM power spectrum in three different cosmologies. Our intention is to determine the first objects that can cool to low temperatures, but not to follow their dynamical evolution. We identify two generations of haloes: those with low virial temperatures,   T _vir≲9000 K  that remain largely neutral, and those with larger virial temperatures that become ionized. The lower temperature, lower mass haloes are the first to cool to 75 per cent of their virial temperature. The precise temperature and mass of the first objects are dependent upon the molecular hydrogen (H₂) cooling function and the cosmological model. The higher mass haloes collapse later but, in this paradigm, cool much more efficiently once they have done so, first via electronic transitions and then via molecular cooling: in fact, a greater residual ionization once the haloes cool below 9000 K results in an enhanced H₂ production and hence a higher cooling rate at low temperatures than for the lower mass haloes, so that within our constant-density model it is the former that are the first to cool to really low temperatures. We discuss the possible significance of this result in the context of CDM models in which the shallow slope of the initial fluctuation spectrum on small scales leads to a wide range of halo masses (of differing overdensities) collapsing over a small redshift interval. This 'crosstalk' is sufficiently important that both high- and low-mass haloes collapse during the lifetimes of the massive stars which may be formed at these epochs. Further investigation is thus required to determine which generation of haloes plays the dominant role in early structure formation.     

The orientation of galaxy dark matter haloes around cosmic voids

Using the Millennium N -body Simulation we explore how the shape and angular momentum of galaxy dark matter haloes surrounding the largest cosmological voids are oriented. We find that the major and intermediate axes of the haloes tend to lie parallel to the surface of the voids, whereas the minor axis points preferentially in the radial direction. We have quantified the strength of these alignments at different radial distances from the void centres. The effect of these orientations is still detected at distances as large as 2.2 R _void from the void centre. Taking a subsample of haloes expected to contain disc-dominated galaxies at their centres we detect, at the 99.9 per cent confidence level, a signal that the angular momentum of those haloes tends to lie parallel to the surface of the voids. Contrary to the alignments of the inertia axes, this signal is only detected in shells at the void surface  (1 < R < 1.07  R _void)  and disappears at larger distances. This signal, together with the similar alignment observed using real spiral galaxies, strongly supports the prediction of the Tidal Torque theory that both dark matter haloes and baryonic matter have acquired, conjointly, their angular momentum before the moment of turnaround.     

On the distribution of haloes, galaxies and mass

The stochasticity in the distribution of dark haloes in the cosmic density field is reflected in the distribution function   P _V ( N _h| δ _m)  , which gives the probability of finding N _h haloes in a volume V with mass density contrast δ _m. We study the properties of this function using high-resolution N -body simulations, and find that   P _V ( N _h| δ _m)  is significantly non-Poisson. The ratio between the variance and the mean goes from ∼1 (Poisson) at  1+ δ _m≪1  to <1 (sub-Poisson) at  1+ δ _m∼1  to >1 (super-Poisson) at  1+ δ _m≫1  . The mean bias relation is found to be well described by halo bias models based on the Press–Schechter formalism. The sub-Poisson variance can be explained as a result of halo exclusion, while the super-Poisson variance at high δ _m may be explained as a result of halo clustering. A simple phenomenological model is proposed to describe the behaviour of the variance as a function of δ _m. Galaxy distribution in the cosmic density field predicted by semi-analytic models of galaxy formation shows similar stochastic behaviour. We discuss the implications of the stochasticity in halo bias to the modelling of higher order moments of dark haloes and of galaxies.     

The Aspen–Amsterdam void finder comparison project

Despite a history that dates back at least a quarter of a century, studies of voids in the large-scale structure of the Universe are bedevilled by a major problem: there exist a large number of quite different void-finding algorithms, a fact that has so far got in the way of groups comparing their results without worrying about whether such a comparison in fact makes sense. Because of the recent increased interest in voids, both in very large galaxy surveys and in detailed simulations of cosmic structure formation, this situation is very unfortunate. We here present the first systematic comparison study of 13 different void finders constructed using particles, haloes, and semi-analytical model galaxies extracted from a subvolume of the Millennium simulation. This study includes many groups that have studied voids over the past decade. We show their results and discuss their differences and agreements. As it turns out, the basic results of the various methods agree very well with each other in that they all locate a major void near the centre of our volume. Voids have very underdense centres, reaching below 10 per cent of the mean cosmic density. In addition, those void finders that allow for void galaxies show that those galaxies follow similar trends. For example, the overdensity of void galaxies brighter than   m _B =−20  is found to be smaller than about −0.8 by all our void finding algorithms.     

Evolution of the mass function of dark matter haloes

We use a high-resolution ΛCDM numerical simulation to calculate the mass function of dark matter haloes down to the scale of dwarf galaxies, back to a redshift of 15, in a  50 h ⁻¹ Mpc  volume containing 80 million particles. Our low-redshift results allow us to probe low-σ density fluctuations significantly beyond the range of previous cosmological simulations. The Sheth & Tormen mass function provides an excellent match to all of our data except for redshifts of 10 and higher, where it overpredicts halo numbers increasingly with redshift, reaching roughly 50 per cent for the  10¹⁰–10¹¹ M_⊙  haloes sampled at redshift 15. Our results confirm previous findings that the simulated halo mass function can be described solely by the variance of the mass distribution, and thus has no explicit redshift dependence. We provide an empirical fit to our data that corrects for the overprediction of extremely rare objects by the Sheth & Tormen mass function. This overprediction has implications for studies that use the number densities of similarly rare objects as cosmological probes. For example, the number density of high-redshift  ( z ≃ 6) QSOs  , which are thought to be hosted by haloes at 5σ peaks in the fluctuation field, are likely to be overpredicted by at least a factor of 50 per cent. We test the sensitivity of our results to force accuracy, starting redshift and halo-finding algorithm.     

The assembly bias of dark matter haloes to higher orders

We use an extremely large volume  (2.4  h ⁻³ Gpc³)  , high-resolution N -body simulation to measure the higher order clustering of dark matter haloes as a function of mass and internal structure. As a result of the large simulation volume and the use of a novel 'cross-moment' counts-in-cells technique which suppresses discreteness noise, we are able to measure the clustering of haloes corresponding to rarer peaks than was possible in previous studies; the rarest haloes for which we measure the variance are 100 times more clustered than the dark matter. We are able to extract, for the first time, halo bias parameters from linear up to fourth order. For all orders measured, we find that the bias parameters are a strong function of mass for haloes more massive than the characteristic mass   M _*  . Currently, no theoretical model is able to reproduce this mass dependence closely. We find that the bias parameters also depend on the internal structure of the halo up to fourth order. For haloes more massive than   M _*  , we find that the more concentrated haloes are more weakly clustered than the less concentrated ones. We see no dependence of clustering on concentration for haloes with masses   M < M _*  ; this is contrary to the trend reported in the literature when segregating haloes by their formation time. Our results are insensitive to whether haloes are labelled by the total mass returned by the friends-of-friends group finder or by the mass of the most massive substructure. This implies that our conclusions are not an artefact of the particular choice of group finding algorithm. Our results will provide important input to theoretical models of galaxy clustering.     

AzTEC Millimetre Survey of the COSMOS field – II. Source count overdensity and correlations with large-scale structure

We report an overdensity of bright submillimetre galaxies (SMGs) in the 0.15 deg² AzTEC/COSMOS survey and a spatial correlation between the SMGs and the optical-IR galaxy density at   z ≲ 1.1  . This portion of the COSMOS field shows a  ∼3σ  overdensity of robust SMG detections when compared to a background, or 'blank-field', population model that is consistent with SMG surveys of fields with no extragalactic bias. The SMG overdensity is most significant in the number of very bright detections (14 sources with measured fluxes   S _1.1 mm > 6  mJy), which is entirely incompatible with sample variance within our adopted blank-field number densities and infers an overdensity significance of  ≫ 4σ  . We find that the overdensity and spatial correlation to optical-IR galaxy density are most consistent with lensing of a background SMG population by foreground mass structures along the line of sight, rather than physical association of the SMGs with the   z ≲ 1.1  galaxies/clusters. The SMG positions are only weakly correlated with weak-lensing maps, suggesting that the dominant sources of correlation are individual galaxies and the more tenuous structures in the survey region, and not the massive and compact clusters. These results highlight the important roles cosmic variance and large-scale structure can play in the study of SMGs.     

Deep galaxy counts, extragalactic background light and the stellar baryon budget

We assess the constraints imposed by the observed extragalactic background light (EBL) on the cosmic history of star formation and the stellar-mass density today. The logarithmic slope of the galaxy number–magnitude relation from the Southern Hubble Deep Field imaging survey is flatter than 0.4 in all seven UBVIJHK optical bandpasses, i.e. the light from resolved galaxies has converged from the UV to the near-IR. We find a lower limit to the surface brightness of the optical extragalactic sky of about 15 nW m⁻² sr⁻¹, comparable to the intensity of the far-IR background from COBE data. Assuming a Salpeter initial mass function with a lower cut-off consistent with observations of M subdwarf disc stars, we set a lower limit of Ω_g+s h ²>0.0013  I ₅₀ to the visible (processed gas + stars) mass density required to generate an EBL at a level of 50  I ₅₀ nW m⁻² sr⁻¹; our 'best-guess' value is Ω_g+s h ²≈0.0031  I ₅₀. Motivated by the recent microlensing results of the MACHO collaboration, we consider the possibility that massive dark haloes around spiral galaxies are composed of faint white dwarfs, and show that only a small fraction (≲5 per cent) of the nucleosynthetic baryons can be locked in the remnants of intermediate-mass stars forming at z _F≲5, as the bright early phases of such haloes would otherwise overproduce the observed EBL.     

The distribution function of dark matter in massive haloes

We study the distribution function (DF) of dark matter particles in haloes of mass range  10¹⁴–10¹⁵ M_⊙  . In the numerical part of this work we measure the DF for a sample of relaxed haloes formed in the simulation of a standard Λ cold dark matter (ΛCDM) model. The DF is expressed as a function of energy E and the absolute value of the angular momentum L , a form suitable for comparison with theoretical models. By proper scaling we obtain the results that do not depend on the virial mass of the haloes. We demonstrate that the DF can be separated into energy and angular momentum components and propose a phenomenological model of the DF in the form     . This formulation involves three parameters describing the anisotropy profile in terms of its asymptotic values (β₀ and  β_∞  ) and the scale of transition between them ( L ₀). The energy part   f _E ( E )  is obtained via inversion of the integral for spatial density. We provide a straightforward numerical scheme for this procedure as well as a simple analytical approximation for a typical halo formed in the simulation. The DF model is extensively compared with the simulations: using the model parameters obtained from fitting the anisotropy profile, we recover the DF from the simulation as well as the profiles of the dispersion and kurtosis of radial and tangential velocities. Finally, we show that our DF model reproduces the power-law behaviour of phase-space density   Q =ρ( r )/σ³( r )  .     

Figure rotation of dark haloes in cold dark matter simulations

We investigate the figure rotation of dark matter haloes identified in Λ cold dark matter (CDM) simulations. We find that when strict criteria are used to select suitable haloes for study, five of the 222 haloes identified in our   z = 0  simulation output undergo coherent figure rotation over a  5 h ⁻¹ Gyr  period. We discuss the effects of varying the selection criteria and find that pattern speeds for a much larger fraction of the haloes can be measured when the criteria are relaxed. Pattern speeds measured over a  1 h ⁻¹ Gyr  period follow a lognormal distribution, centred at  Ω_p= 0.2 h rad Gyr⁻¹  with a maximum value of 0.94 h rad Gyr⁻¹. Over a  5 h ⁻¹ Gyr  period, the average pattern speed of a halo is about  0.1 h rad Gyr⁻¹  and the largest pattern speed found is  0.24 h rad Gyr⁻¹  . Less than half of the selected haloes showed alignment between their figure rotation axis and minor axis, the exact fraction being somewhat dependent on how one defines a halo. While the pattern speeds observed are lower than those generally thought capable of causing spiral structure, we note that coherent figure rotation is found over very long periods and argue that further simulations would be required before strong conclusions about spiral structure in all galaxies could be drawn. We find no correlation between halo properties such as total mass and the pattern speed.     

Mass models of the Milky Way

A parametrized model of the mass distribution within the Milky Way is fitted to the available observational constraints. The most important single parameter is the ratio of the scalelength R _d* of the stellar disc to R ₀. The disc and bulge dominate v _c( R ) at R ≲ R ₀ only for R _d,*/ R ₀≲0.3. Since the only knowledge we have of the halo derives from studies like the present one, we allow it to contribute to the density at all radii. When allowed this freedom, however, the halo causes changes in assumptions relating to R  ≪  R ₀ to affect profoundly the structure of the best-fitting model at R  ≫  R ₀. For example, changing the disc slightly from an exponential surface-density profile significantly changes the form of v _c( R ) at R  ≫  R ₀, where the disc makes a negligible contribution to v _c. Moreover, minor changes in the constraints can cause the halo to develop a deep hole at its centre that is not physically plausible. These problems call into question the proposition that flat rotation curves arise because galaxies have physically distinct haloes rather than outwards-increasing mass-to-light ratios.   The mass distribution of the Galaxy and the relative importance of its various components will remain very uncertain until more observational data can be used to constrain mass models. Data that constrain the Galactic force field at z ≳ R and at R  >  R ₀ are especially important.     

Scale-dependent bias induced by local non-Gaussianity: a comparison to N-body simulations

We investigate the effect of primordial non-Gaussianity of the local f _NL type on the auto- and cross-power spectra of dark matter haloes using simulations of the Λ cold dark matter cosmology. We perform a series of large N -body simulations of both positive and negative f _NL, spanning the range between 10 and 100. Theoretical models predict a scale-dependent bias correction  Δ b ( k , f _NL)  that depends on the linear halo bias   b ( M )  . We measure the power spectra for a range of halo mass and redshifts covering the relevant range of existing galaxy and quasar populations. We show that auto- and cross-correlation analyses of bias are consistent with each other. We find that for low wavenumbers with   k < 0.03  h  Mpc⁻¹  the theory and the simulations agree well with each other for biased haloes with   b ( M ) > 1.5  . We show that a scale-independent bias correction improves the comparison between theory and simulations on smaller scales, where the scale-dependent effect rapidly becomes negligible. The current limits on f _NL from Slosar et al. come mostly from very large scales   k < 0.01  h  Mpc⁻¹  and, therefore, remain valid. For the halo samples with   b ( M ) < 1.5 − 2  , we find that the scale-dependent bias from non-Gaussianity actually exceeds the theoretical predictions. Our results are consistent with the bias correction scaling linearly with f _NL.     

An excess of damped Lyman α galaxies near quasi-stellar objects

We present a sample of 33 damped Lyman α systems (DLAs) discovered in the Sloan Digital Sky Survey (SDSS) whose absorption redshifts ( z _abs) are within 6000 km s⁻¹ of the quasi-stellar object's (QSO) systemic redshift ( z _sys). Our sample is based on  731 2.5 < z _sys < 4.5  non-broad absorption line (non-BAL) QSOs from Data Release 3 (DR3) of the SDSS. We estimate that our search is ≈100 per cent complete for absorbers with N (H  i )  ≥ 2 × 10²⁰ cm⁻²  . The derived number density of DLAs per unit redshift, n ( z ), within  Δ v < 6000 km s⁻¹  is higher (3.5σ significance) by almost a factor of 2 than that of intervening absorbers observed in the SDSS DR3, i.e. there is evidence for an overdensity of galaxies near the QSOs. This provides a physical motivation for excluding DLAs at small velocity separations in surveys of intervening 'field' DLAs. In addition, we find that the overdensity of proximate DLAs is independent of the radio-loudness of the QSO, consistent with the environments of radio-loud and radio-quiet QSOs being similar.     

Molecular hydrogen in damped Lyα systems: clues to interstellar physics at high redshift

In order to interpret H₂ quasar absorption-line observations of damped Lyα systems (DLAs) and subDLAs, we model their H₂ abundance as a function of dust-to-gas ratio, including H₂ self-shielding and dust extinction against dissociating photons. Then, we constrain the physical state of the gas by using H₂ data. Using H₂ excitation data for DLAs with H₂ detections, we derive a gas density  1.5 ≲ log n (cm⁻³) ≲ 2.5  , temperature  1.5 ≲ log T (K) ≲ 3  , and an internal ultraviolet (UV) radiation field (in units of the Galactic value)  0.5 ≲ log χ≲ 1.5  . We then find that the observed relation between the molecular fraction and the dust-to-gas ratio of the sample is naturally explained by the above conditions. However, it is still possible that H₂ deficient DLAs and subDLAs with H₂ fractions less than  ∼10⁻⁶  are in a more diffuse and warmer state. The efficient photodissociation by the internal UV radiation field explains the extremely small H₂ fraction  (≲10⁻⁶)  observed for  κ≲ 1/30  (κ is the dust-to-gas ratio in units of the Galactic value); H₂ self-shielding causes a rapid increase in, and large variations of, H₂ abundance for  κ≳ 1/30  . We finally propose an independent method to estimate the star formation rates of DLAs from H₂ abundances; such rates are then critically compared with those derived from other proposed methods. The implications for the contribution of DLAs to the cosmic star formation history are briefly discussed.     

A constant dark matter halo surface density in galaxies

We confirm and extend the recent finding that the central surface density  μ_0D≡ r ₀ρ₀  of galaxy dark matter haloes, where r ₀ and  ρ₀  are the halo core radius and central density, is nearly constant and independent of galaxy luminosity. Based on the co-added rotation curves (RCs) of ∼1000 spiral galaxies, the mass models of individual dwarf irregular and spiral galaxies of late and early types with high-quality RCs, and the galaxy–galaxy weak-lensing signals from a sample of spiral and elliptical galaxies, we find that  log μ_0D= 2.15 ± 0.2  in units of  log(M_⊙ pc⁻²)  . We also show that the observed kinematics of Local Group dwarf spheroidal galaxies are consistent with this value. Our results are obtained for galactic systems spanning over 14 mag, belonging to different Hubble types and whose mass profiles have been determined by several independent methods. In the same objects, the approximate constancy of  μ_0D  is in sharp contrast to the systematical variations, by several orders of magnitude, of galaxy properties, including  ρ₀  and central stellar surface density.     

Peculiar velocities of galaxies and clusters

We present a simple model for the shape of the distribution function of galaxy peculiar velocities. We show how both non-linear and linear theory terms combine to produce a distribution which has an approximately Gaussian core with exponential wings. The model is easily extended to study how the statistic depends on the type of particle used to trace the velocity field (dark matter particles, dark matter haloes, galaxies), and on the density of the environment in which the test particles are located. Comparisons with simulations suggest that our model is accurate. We also show that the evolution of the peculiar velocities depends on the local, rather than the global, density. Since clusters populate denser regions on average, using cluster velocities with the linear theory scaling may lead to an overestimate of the global value of Ω₀. Conversely, using linear theory with the global value of Ω₀ to scale cluster velocities from the initial to the present time results in an underestimate of their true velocities. In general, however, the directions of motions of haloes are rather well described by linear theory. Our results help to simplify models of redshift-space distortions considerably.