The Mass Function of Dark Halos in Superclusters and Voids

A modification of the Press-Schechter theory allowing for presence of a background large-scale structure (LSS) - a supercluster or a void, is proposed. The LSS is accounted as the statistical constraints in form of linear functionals of the random overdensity field. The deviation of the background density within the LSS is interpreted in a pseudo-cosmological sense. Using the constraints formalism may help us to probe non-trivial spatial statistics of haloes, e.g. edge and shape effects on boundaries of the superclusters and voids. Parameters of the constraints are connected to features of the LSS: its mean overdensity, a spatial scale and a shape, and spatial momenta of higher orders. It is shown that presence of a non-virialized LSS can lead to an observable deviation of the mass function. This effect is exploited to build a procedure to recover parameters of the background perturbation from the observationally estimated mass function.Comment: 23 pages, 6 figures; to be appeared in Astronomy Reports, 2014, Vol. 58, No. 6, pp. 386-39

The Small Scale Velocity Dispersion of Galaxies: A Comparison of Cosmological Simulations

The velocity dispersion of galaxies on small scales ($r\sim1h^{-1}$ Mpc), $\sigma_{12}(r)$, can be estimated from the anisotropy of the galaxy-galaxy correlation function in redshift space. We apply this technique to ``mock-catalogs'' extracted from N-body simulations of several different variants of Cold Dark Matter dominated cosmological models to obtain results which may be consistently compared to similar results from observations. We find a large variation in the value of $\sigma_{12}(1 h^{-1} Mpc)$ in different regions of the same simulation. We conclude that this statistic should not be considered to conclusively rule out any of the cosmological models we have studied. We attempt to make the statistic more robust by removing clusters from the simulations using an automated cluster-removing routine, but this appears to reduce the discriminatory power of the statistic. However, studying $\sigma_{12}$ as clusters with different internal velocity dispersions are removed leads to interesting information about the amount of power on cluster and subcluster scales. We also compute the pairwise velocity dispersion directly and compare this to the values obtained using the Davis-Peebles method, and find that the agreement is fairly good. We evaluate the models used for the mean streaming velocity and the pairwise peculiar velocity distribution in the original Davis-Peebles method by comparing the models with the results from the simulations.Comment: 20 pages, uuencoded (Latex file + 8 Postscript figures), uses AAS macro

Conditional Mass Functions and Merger Rates of Dark Matter Halos in the Ellipsoidal Collapse Model

Analytic models based on spherical and ellipsoidal gravitational collapse have been used to derive the mass functions of dark matter halos and their progenitors (the conditional mass function). The ellipsoidal model generally provides a better match to simulation results, but there has been no simple analytic expression in this model for the conditional mass function that is accurate for small time steps, a limit that is important for generating halo merger trees and computing halo merger rates. We remedy the situation by deriving accurate analytic formulae for the first-crossing distribution, the conditional mass function, and the halo merger rate in the ellipsoidal collapse model in the limit of small look-back times. We show that our formulae provide a closer match to the Millennium simulation results than those in the spherical collapse model and the ellipsoidal model of Sheth & Tormen (2002).Comment: 5 pages, 3 figures, accepted by MNRAS letter

The DEEP2 Galaxy Redshift Survey: The Evolution of Void Statistics from z~1 to z~0

We present measurements of the void probability function (VPF) at z~1 using data from the DEEP2 Redshift Survey and its evolution to z~0 using data from the Sloan Digital Sky Survey (SDSS). We measure the VPF as a function of galaxy color and luminosity in both surveys and find that it mimics trends displayed in the two-point correlation function, $\xi$; namely that samples of brighter, red galaxies have larger voids (i.e. are more strongly clustered) than fainter, blue galaxies. We also clearly detect evolution in the VPF with cosmic time, with voids being larger in comoving units at z~0. We find that the reduced VPF matches the predictions of a `negative binomial' model for galaxies of all colors, luminosities, and redshifts studied. This model lacks a physical motivation, but produces a simple analytic prediction for sources of any number density and integrated two-point correlation function, \bar{\xi}. This implies that differences in the VPF across different galaxy populations are consistent with being due entirely to differences in the population number density and \bar{\xi}. The robust result that all galaxy populations follow the negative binomial model appears to be due to primarily to the clustering of dark matter halos. The reduced VPF is insensitive to changes in the parameters of the halo occupation distribution, in the sense that halo models with the same \bar{\xi} will produce the same VPF. For the wide range of galaxies studied, the VPF therefore does not appear to provide useful constraints on galaxy evolution models that cannot be gleaned from studies of \bar{\xi} alone. (abridged)Comment: 17 pages, 15 figures, ApJ accepte

Precision Determination of the Mass Function of Dark Matter Halos

The predicted mass function of dark matter halos is essential in connecting observed galaxy cluster counts and models of galaxy clustering to the properties of the primordial density field. We determine the mass function in the concordance $\Lambda$CDM cosmology, as well as its uncertainty, using sixteen $1024^3$-particle nested-volume dark-matter simulations, spanning a mass range of over five orders of magnitude. Using the nested volumes and single-halo tests, we find and correct for a systematic error in the friends-of-friends halo-finding algorithm. We find a fitting form and full error covariance for the mass function that successfully describes the simulations' mass function and is well-behaved outside the simulations' resolutions. Estimated forecasts of uncertainty in cosmological parameters from future cluster count surveys have negligible contribution from remaining statistical uncertainties in the central cosmology multiplicity function. There exists a potentially non-negligible cosmological dependence (non-universality) of the halo multiplicity function.Comment: 4 pages, 3 figures, submitted to ApJ

Efficient Simulations of Early Structure Formation and Reionization

We present a method to construct semi-numerical ``simulations'', which can efficiently generate realizations of halo distributions and ionization maps at high redshifts. Our procedure combines an excursion-set approach with first-order Lagrangian perturbation theory and operates directly on the linear density and velocity fields. As such, the achievable dynamic range with our algorithm surpasses the current practical limit of N-body codes by orders of magnitude. This is particularly significant in studies of reionization, where the dynamic range is the principal limiting factor. We test our halo-finding and HII bubble-finding algorithms independently against N-body simulations with radiative transfer and obtain excellent agreement. We compute the size distributions of ionized and neutral regions in our maps. We find even larger ionized bubbles than do purely analytic models at the same volume-weighted mean hydrogen neutral fraction. We also generate maps and power spectra of 21-cm brightness temperature fluctuations, which for the first time include corrections due to gas bulk velocities. We find that velocities widen the tails of the temperature distributions and increase small-scale power, though these effects quickly diminish as reionization progresses. We also include some preliminary results from a simulation run with the largest dynamic range to date: a 250 Mpc box that resolves halos with masses M >~ 2.2 x10^8 M_sun. We show that accurately modeling the late stages of reionization requires such large scales. The speed and dynamic range provided by our semi-numerical approach will be extremely useful in the modeling of early structure formation and reionization.Comment: 13 pages, 10 figures; ApJ submitte

Close Pairs as Proxies for Galaxy Cluster Mergers

Galaxy cluster merger statistics are an important component in understanding the formation of large-scale structure. Unfortunately, it is difficult to study merger properties and evolution directly because the identification of cluster mergers in observations is problematic. We use large N-body simulations to study the statistical properties of massive halo mergers, specifically investigating the utility of close halo pairs as proxies for mergers. We examine the relationship between pairs and mergers for a wide range of merger timescales, halo masses, and redshifts (0<z<1). We also quantify the utility of pairs in measuring merger bias. While pairs at very small separations will reliably merge, these constitute a small fraction of the total merger population. Thus, pairs do not provide a reliable direct proxy to the total merger population. We do find an intriguing universality in the relation between close pairs and mergers, which in principle could allow for an estimate of the statistical merger rate from the pair fraction within a scaled separation, but including the effects of redshift space distortions strongly degrades this relation. We find similar behavior for galaxy-mass halos, making our results applicable to field galaxy mergers at high redshift. We investigate how the halo merger rate can be statistically described by the halo mass function via the merger kernel (coagulation), finding an interesting environmental dependence of merging: halos within the mass resolution of our simulations merge less efficiently in overdense environments. Specifically, halo pairs with separations less than a few Mpc/h are more likely to merge in underdense environments; at larger separations, pairs are more likely to merge in overdense environments.Comment: 12 pages, 9 figures; Accepted for publication in ApJ. Significant additions to text and two figures changed. Added new findings on the universality of pair mergers and added analysis of the effect of FoF linking length on halo merger

Scale-dependent Galaxy Bias

We present a simple heuristic model to demonstrate how feedback related to the galaxy formation process can result in a scale-dependent bias of mass versus light, even on very large scales. The model invokes the idea that galaxies form initially in locations determined by the local density field, but the subsequent formation of galaxies is also influenced by the presence of nearby galaxies that have already formed. The form of bias that results possesses some features that are usually described in terms of stochastic effects, but our model is entirely deterministic once the density field is specified. Features in the large-scale galaxy power spectrum (such as wiggles that might in an extreme case mimic the effect of baryons on the primordial transfer function) could, at least in principle, arise from spatial modulations of the galaxy formation process that arise naturally in our model. We also show how this fully deterministic model gives rise to apparently stochasticity in the galaxy distribution.Comment: 14 pages, 2 figures, typos corrected, discussion added and references corrected; matches version accepted by JCA

Gravity and Large-Scale Non-local Bias

The relationship between galaxy and matter overdensities, bias, is most often assumed to be local. This is however unstable under time evolution, we provide proofs under several sets of assumptions. In the simplest model galaxies are created locally and linearly biased at a single time, and subsequently move with the matter (no velocity bias) conserving their comoving number density (no merging). We show that, after this formation time, the bias becomes unavoidably non-local and non-linear at large scales. We identify the non-local gravitationally induced fields in which the galaxy overdensity can be expanded, showing that they can be constructed out of the invariants of the deformation tensor (Galileons). In addition, we show that this result persists if we include an arbitrary evolution of the comoving number density of tracers. We then include velocity bias, and show that new contributions appear, a dipole field being the signature at second order. We test these predictions by studying the dependence of halo overdensities in cells of fixed matter density: measurements in simulations show that departures from the mean bias relation are strongly correlated with the non-local gravitationally induced fields identified by our formalism. The effects on non-local bias seen in the simulations are most important for the most biased halos, as expected from our predictions. The non-locality seen in the simulations is not fully captured by assuming local bias in Lagrangian space. Accounting for these effects when modeling galaxy bias is essential for correctly describing the dependence on triangle shape of the galaxy bispectrum, and hence constraining cosmological parameters and primordial non-Gaussianity. We show that using our formalism we remove an important systematic in the determination of bias parameters from the galaxy bispectrum, particularly for luminous galaxies. (abridged)Comment: 26 pages, 9 figures. v2: improved appendix

Visualizing climate change impact with ubiquitous spatial technologies

Copyright confirmation in progress. Any queries to [email protected] Ma