Cosmological perturbation theory and the spherical collapse model — III. The velocity divergence field and the Ω dependence

Cosmological perturbation theory (PT) is a useful tool to study the cumulants of the density and velocity fields in the large-scale structure of the Universe. In Papers I and II of this series we saw that the spherical collapse (SC) model provides the exact solution to PT at tree-level and gives a good approximation to the loop corrections (next-to-leading orders), indicating negligible tidal effects. Here, we derive predictions for the (smoothed) cumulants of the velocity divergence field θ ≡ ▽ ⊙  v for an irrotational fluid in the SC model. By comparing these with the exact analytic results of Scoccimarro &38; Frieman, it is shown that, at least for the unsmoothed case, the loop corrections to the cumulants of θ are dominated by tidal effects. However, most of the tidal contribution seems to cancel out when computing the hierarchical ratios, T _J  = 〈θ ^J 〉 / 〈θ²〉  ^{J −1}. We also extend the work presented in Papers I and II to give predictions for the cumulants of the density and velocity divergence fields in non-flat spaces. In particular, we show the equivalence between the spherically symmetric solution to the equations of motion in the SC model (given in terms of the density) and that of the Lagrangian PT approach (given in terms of the displacement field). It is shown that the Ω dependence is very weak for both cosmic fields even at one loop (a 10 per cent effect at most), except for the overall factor f (Ω) that couples to the velocity divergence.     

The non-linear redshift-space power spectrum of galaxies

We study the power spectrum of galaxies in redshift space, with third-order perturbation theory to include corrections that are absent in linear theory. We assume a local bias for the galaxies: i.e., the galaxy density is sampled from some local function of the underlying mass distribution. We find that the effect of the non-linear bias in real space is to introduce two new features: first, there is a contribution to the power which is constant with wavenumber, whose nature we reveal as essentially a shot-noise term. In principle this contribution can mask the primordial power spectrum, and could limit the accuracy with which the latter might be measured on very large scales. Secondly, the effect of second- and third-order bias is to modify the effective bias (defined as the square root of the ratio of galaxy power spectrum to matter power spectrum). The effective bias is almost scale-independent over a wide range of scales. These general conclusions also hold in redshift space. In addition, we have investigated the distortion of the power spectrum by peculiar velocities, which may be used to constrain the density of the Universe. We look at the quadrupole-to-monopole ratio, and find that higher order terms can mimic linear theory bias, but the bias implied is neither the linear bias, nor the effective bias referred to above. We test the theory with biased N -body simulations, and find excellent agreement in both real and redshift space, providing the local biasing is applied on a scale whose fractional rms density fluctuations are < 0.5.     

An excursion set model for the distribution of dark matter and dark matter haloes

A model of the gravitationally evolved dark matter distribution, in the Eulerian space, is developed. It is a simple extension of the excursion set model that is commonly used to estimate the mass function of collapsed dark matter haloes. In addition to describing the evolution of the Eulerian space distribution of the haloes, the model allows one to describe the evolution of the dark matter itself. It can also be used to describe density profiles, on scales larger than the virial radius of these haloes, and to quantify the way in which matter flows in and out of Eulerian cells. When the initial Lagrangian space distribution is white noise Gaussian, the model suggests that the Inverse Gaussian distribution should provide a reasonably good approximation to the evolved Eulerian density field, in agreement with numerical simulations. Application of this model to clustering from more general Gaussian initial conditions is discussed at the end.     

Gamma-ray bursts from stellar remnants: probing the Universe at high redshift

A gamma-ray burst (GRB) releases an amount of energy similar to that of a supernova explosion, which combined with its rapid variability suggests an origin related to neutron stars or black holes. Since these compact stellar remnants form from the most massive stars not long after their birth, GRBs should trace the star formation rate in the Universe; we show that the GRB flux distribution is consistent with this. Because of the strong evolution of the star formation rate with redshift, it follows that the dimmest known bursts have z  ∼ 6, much above the value usually quoted and beyond the most distant quasars. This explains the absence of bright galaxies in well-studied GRB error boxes. The increased distances imply a peak luminosity of 8.3 × 10⁵¹ erg s⁻¹ and a rate density of 0.025 per million years per galaxy. These values are 20 times higher and 150 times lower, respectively, than are implied by fits with non-evolving GRB rates. This means either that GRBs are caused by a much rarer phenomenon than mergers of binary neutron stars, or that their gamma-ray emission is often invisible to us due to beaming. Precise burst locations from optical transients will discriminate between the various models for GRBs from stellar deaths, because the distance between progenitor birth place and burst varies greatly among them. The dimmest GRBs are then the most distant known objects, and may probe the Universe at an age when the first stars were forming.     

The UVX quasar optical luminosity function and its evolution

The recently finished Edinburgh UVX quasar survey at B <18 is used together with other complete samples to estimate the shape and evolution of the optical luminosity function in the redshift range 0.3< z <2.2. There is a significantly higher space density of quasars at high luminosity and low redshift than previously found in the PG sample of Schmidt & Green, with the result that the shape of the luminosity function at low redshifts ( z <1) is seen to be consistent with a single power law. At higher redshifts the slope of the power law at high luminosities appears to steepen significantly. There does not appear to be any consistent break feature which could be used as a tracer of luminosity evolution in the population.     

Red quasars not so dusty

Webster et al. claimed that up to 80 per cent of QSOs may be obscured by dust. They inferred the presence of this dust from the remarkably broad range of B − K optical–infrared colours of a sample of flat-spectrum PKS radio QSOs. If such dust is typical of QSOs, it will have rendered invisible most of those which would otherwise have been detected by optical surveys. We used the William Herschel Telescope on La Palma to obtain K infrared images of 54 B3 radio QSOs selected at low frequency (mainly steep-spectrum), and we find that although several have very red optical–infrared colours, most of these can be attributed to an excess of light in K rather than a dust-induced deficit in B. We present evidence that some of the infrared excess comes from the light of stars in the host galaxy (some, as previously suggested, comes from synchrotron radiation associated with flat-spectrum radio sources). The B − K colours of the B3 QSOs provide no evidence for a large reddened population. Either the Webster et al. QSOs are atypical in having such large extinctions, or their reddening is not due to dust; either way, the broad range of their B − K colours does not provide evidence that a large fraction of QSOs has been missed from optical surveys.