Profiles of dark matter haloes at high redshift

I study the evolution of halo density profiles as a function of time in the SCDM and ΛCDM cosmologies. Following Del Popolo, I calculate the concentration parameter c = r _v / a and study its time evolution. For a given halo mass, I find that c ( z ) ∝ 1/(1+ z ) in both the ΛCDM and SCDM cosmology, in agreement with the analytic model of Bullock et al. and N -body simulations. In both models, a ( z ) is roughly constant. The present model predicts a stronger evolution of c ( z ) with respect to the Navarro, Frenk & White model. Finally I show some consequences of the results on galaxy modelling.     

Merger histories in warm dark matter structure formation scenarios

Observations on galactic scales seem to be in contradiction with recent high-resolution N -body simulations. This so-called cold dark matter (CDM) crisis has been addressed in several ways, ranging from a change in fundamental physics by introducing self-interacting cold dark matter particles to a tuning of complex astrophysical processes such as global and/or local feedback. All these efforts attempt to soften density profiles and reduce the abundance of satellites in simulated galaxy haloes. In this paper, we explore a different approach that consists of filtering the dark matter power spectrum on small scales, thereby altering the formation history of low-mass objects. The physical motivation for damping these fluctuations lies in the possibility that the dark matter particles have a different nature, i.e. are warm (WDM) rather than cold. We show that this leads to some interesting new results in terms of the merger history and large-scale distribution of low-mass haloes, compared with the standard CDM scenario. However, WDM does not appear to be the ultimate solution, in the sense that it is not able to fully solve the CDM crisis, even though one of the main drawbacks, namely the abundance of satellites, can be remedied. Indeed, the cuspiness of the halo profiles still persists, at all redshifts, and for all haloes and sub-haloes that we investigated. Despite the persistence of the cuspiness problem of DM haloes, WDM seems to be still worth taking seriously, as it alleviates the problems of over-abundant sub-structures in galactic haloes and possibly the lack of angular momentum of simulated disc galaxies. WDM also lessens the need to invoke strong feedback to solve these problems, and may provide a natural explanation of the clustering properties and ages of dwarfs.     

Foregrounds for 21-cm observations of neutral gas at high redshift

We investigate a number of potential foregrounds for an ambitious goal of future radio telescopes such as the Square Kilometer Array (SKA) and the Low Frequency Array (LOFAR): spatial tomography of neutral gas at high redshift in 21-cm emission. While the expected temperature fluctuations due to unresolved radio point sources is highly uncertain, we point out that free–free emission from the ionizing haloes that reionized the Universe should define a minimal bound. This emission is likely to swamp the expected brightness temperature fluctuations, making proposed detections of the angular patchwork of 21-cm emission across the sky unlikely to be viable. Hα observations with JWST could place an upper bound on the contribution of high-redshift sources to the free–free background. An alternative approach is to discern the topology of reionization from spectral features due to 21-cm emission along a pencil-beam slice. This requires tight control of the frequency-dependence of the beam in order to prevent foreground sources from contributing excessive variance. We also investigate potential contamination by galactic and extragalactic radio recombination lines (RRLs). These are unlikely to be show-stoppers, although little is known about the distribution of RRLs away from the Galactic plane. The mini-halo emission signal is always less than that of the intergalactic medium (IGM), making mini-haloes unlikely to be detectable. If they are seen, it will be only in the very earliest stages of structure formation at high redshift, when the spin temperature of the IGM has not yet decoupled from the cosmic microwave background.     

Evidence for merger-driven activity in the clustering of high-redshift quasars

Recently, a very large clustering length has been measured for quasars at a redshift of   z ∼ 4  . In combination with the observed quasar luminosity function, we assess the implications of this clustering for the relationship between quasar luminosity and dark matter halo mass. Our analysis allows for non-linearity and finite scatter in the relation between quasar luminosity and halo mass, as well as a luminosity dependent quasar lifetime. The additional novel ingredient in our modelling is the allowance for an excess in the observed bias over the underlying halo bias owing to the merger driven nature of quasar activity. We find that the observations of clustering and luminosity function can be explained only if both of the following conditions hold: (i) the luminosity to halo mass ratio increases with halo mass; (ii) the observed clustering amplitude is in excess of that expected solely from halo bias. The latter result is statistically significant at the 99 per cent level. Taken together, the observations provide compelling evidence for merger driven quasar activity, with a black-hole growth that is limited by feedback. In difference from previous analyses, we show that there could be scatter in the luminosity–halo mass relation of up to 1 dex, and that quasar clustering cannot be used to estimate the quasar lifetime.     

Fossil H ii regions: self-limiting star formation at high redshift

Recent results from the Wilkinson Microwave Anisotropy Probe ( WMAP ) satellite suggest that the intergalactic medium (IGM) was significantly reionized at redshifts as high as   z ∼ 17  . At this early epoch, the first ionizing sources probably appeared in the shallow potential wells of mini-haloes with virial temperatures   T _vir < 10⁴ K  . Once such an ionizing source turns off, its surrounding H ii region Compton cools and recombines. None the less, we show that the 'fossil' H ii regions left behind remain at high adiabats, prohibiting gas accretion and cooling in subsequent generations of mini-haloes. This greatly amplifies feedback effects explored in previous studies, and early star formation is self-limiting. We quantify this effect to show that star formation in mini-haloes cannot account for the bulk of the electron scattering opacity measured by WMAP , which must be due to more massive objects. We argue that gas entropy, rather than IGM metallicity, regulates the evolution of the global ionizing emissivity and impedes full reionization until lower redshifts. We discuss several important consequences of this early entropy floor for reionization. It reduces gas clumping, curtailing the required photon budget for reionization. An entropy floor also prevents H₂ formation and cooling, due to reduced gas densities: it greatly enhances feedback from ultraviolet photodissociation of H₂. An early X-ray background would also furnish an entropy floor to the entire IGM; thus, X-rays impede rather than enhance H₂ formation. Future 21-cm observations may probe the topology of fossil H ii regions.     

Peaks in the cosmological density field: sensitivity to initial power spectrum, redshift distortions and galaxy halo occupation

We investigate the number density of maxima in the cosmological galaxy density field smoothed with a filter as a probe of clustering. In previous work it has been shown that this statistic is closely related to the slope of the linear power spectrum, even when the directly measured power spectrum is non-linear. In the present paper we investigate the sensitivity of the peak number density to various models with differing power spectra, including rolling index models, cosmologies with massive neutrinos and different baryon densities. We find that rolling index models which have given an improved fit to CMB/LSS (cosmic microwave background/large scale structure) data yield a ∼10 per cent difference in peak density compared to the scale invariant case. Models with 0.3 eV neutrinos have effects of similar magnitude and it should be possible to constrain them with data from current galaxy redshift surveys. Baryon oscillations in the power spectrum also give rise to distinctive features in the peak density. These are preserved without modification when measured from the peak density in fully non-linear N -body simulations. Using the simulations, we also investigate how the peak density is modified in the presence of redshift distortions. Redshift distortions cause a suppression of the number of peaks, largely due to fingers of God overlapping in redshift space. We find that this effect can be modelled by using a modification of the input power spectrum. We also study the results when the simulation density field is traced by galaxies obtained by populating haloes with a halo occupation distribution consistent with observations. The peak number density is consistent with that in the dark matter for filter scales  >4  h ⁻¹ Mpc  , for which we find good agreement with the linear theory predictions. In a companion paper we analyse data from the 2dF Galaxy Redshift Survey.     

Constraints on the clustering, biasing and redshift distribution of radio sources

We discuss how different theoretical predictions for the variance σ ² of the counts-in-cells distribution of radio sources can be matched to measurements from the FIRST survey at different flux limits. The predictions are given by the integration of models for the angular correlation function w ( θ ) for three different functional forms of the redshift distribution N ( z ), different spatial correlation functions that match the observed present-day shape and by different evolutions of the bias b ( z ) with redshift. We also consider the two cases of open and flat universes. Although the predicted w ( θ ) show substantial differences because of differences in the values of N ( z ), these differences are not significant compared to the uncertainties in the current observations. It turns out that, independent of the geometry of the universe and the flux limit, the best fit is provided by models with constant biasing at all times, although the difference between models with epoch-independent bias and models with bias that evolves linearly with redshift is not very large. All models with strong evolution of bias with epoch are ruled out, as they grossly overestimate the amplitude of the variance over the whole range of angular scales sampled by the counts-in-cells analysis. As a further step we directly calculated w _obs( θ ) at 3 mJy from the catalogue and matched it with our models for the angular correlation function, in the hypothesis that the clustering signal comes from two different populations, namely AGN-powered sources and starbursting galaxies. The results are consistent with a scenario for hierarchical clustering where the fainter starbursting galaxies trace the mass at all epochs, while the brighter AGNs are strongly biased, with b ( z ) evolving linearly with redshift, as suggested by some theories of galaxy formation and evolution.     

Non-standard structure formation scenarios

Observations on galactic scales seem to be in contradiction with recent high resolution N-body simulations. This so-called cold dark matter (CDM) crisis has been addressed in several ways, ranging from a change in fundamental physics by introducing self-interacting cold dark matter particles to a tuning of complex astrophysical processes such as global and/or local feedback. All these efforts attempt to soften density profiles and reduce the abundance of satellites in simulated galaxy halos. In this contribution we are exploring the differences between a Warm Dark Matter model and a CDM model where the power on a certain scale is reduced by introducing a narrow negative feature (`dip'). This dip is placed in a way so as to mimic the loss of power in the WDM model: both models have the same integrated power out to the scale where the power of the Dip model rises to the level of the unperturbed CDM spectrum again. Using N-body simulations we show that that the new Dip model appears to be a viable alternative to WDM while being based on different physics: where WDM requires the introduction of a new particle species the Dip stems from anon-standard inflationary period. If we are looking for an alternative to the currently challenged standard ΛCDM structure formation scenario, neither the ΛWDM nor the new Dip model can be ruled out with respect to the analysis presented in this contribution. They both make very similar predictions and the degeneracy between them can only be broken with observations yet to come. This revised version was published online in August 2006 with corrections to the Cover Date.     

Effects of correlation between merging steps on the global halo formation

The excursion set theory of halo formation is modified by adopting the fractional Brownian motion, to account for possible correlation between merging steps. We worked out analytically the conditional mass function, halo merging rate and formation time distribution in the spherical collapse model. We also developed an approximation for the ellipsoidal collapse model and applied it to the calculation of the conditional mass function and the halo formation time distribution. For models in which the steps are positively correlated, the halo merger rate is enhanced when the accreted mass is less than  ∼25 M *  , while for the negatively correlated case this rate is reduced. Compared with the standard model in which the steps are uncorrelated, the models with positively correlated steps produce more aged population in small mass haloes and more younger population in large mass haloes, while for the models with negatively correlated steps the opposite is true. An examination of simulation results shows that a weakly positive correlation between successive merging steps appears to fit best. We have also found a systematic effect in the measured mass function due to the finite volume of simulations. In future work, this will be included in the halo model to accurately predict the three-point correlation function estimated from simulations.     

Tidal alignments as a contaminant of redshift space distortions

We investigate the effect of orientation-dependent selection effects on galaxy clustering in redshift space. It is found that if galaxies are aligned by large-scale tidal fields, then these selection effects give rise to a dependence of the observed galaxy density on the local tidal field, in addition to the well-known dependences on the matter density and radial velocity gradient. This alters the galaxy power spectrum in a way that is different for Fourier modes parallel to and perpendicular to the line of sight. These tidal galaxy alignments can thus mimic redshift space distortions (RSD), and thus result in a bias in the measurement of the velocity power spectrum. If galaxy orientations are affected only by the local tidal field, then the tidal alignment effect has exactly the same scale and angular dependence as the RSDs in the linear regime, so it cannot be projected out or removed by masking small scales in the analysis. We consider several toy models of tidal alignments and orientation-dependent selection, normalize their free parameter (an amplitude) to recent observations, and find that they could bias the velocity amplitude   f ( z ) G ( z )  by 5–10 per cent in some models, although most models give much smaller contamination. We conclude that tidal alignments may be a significant systematic error in RSD measurements that aim to test general relativity via the growth of large-scale structure. We briefly discuss possible mitigation strategies.     

Linear redshift distortions and power in the IRAS Point Source Catalog Redshift Survey

We present a state-of-the-art linear redshift distortion analysis of the recently published IRAS Point Source Catalog Redshift Survey (PSC z ). The procedure involves linear compression into 4096 KarhunenLoève (signal-to-noise) modes culled from a potential pool of 310⁵ modes, followed by quadratic compression into three separate power spectra, the galaxygalaxy, galaxyvelocity and velocityvelocity power spectra. Least squares-fitting to the decorrelated power spectra yields a linear redshift distortion parameter