The Integrated Sachs-Wolfe effect in $f(R)$ gravity

We study the late-time Integrated Sachs-Wolfe (ISW) effect in $f(R)$ gravity using N-body simulations. In the $f(R)$ model under study, the linear growth rate is larger than that in general relativity (GR). This slows down the decay of the cosmic potential and induces a smaller ISW effect on large scales. Therefore, the $dotPhi$ (time derivative of the potential) power spectrum at $k<0.1h$/Mpc is suppressed relative to that in GR. In the non-linear regime, relatively rapid structure formation in $f(R)$ gravity boosts the non-linear ISW effect relative to GR, and the $dotPhi$ power spectrum at $k>0.1h$/Mpc is increased (100$%$ greater on small scales at $z=0$). We explore the detectability of the ISW signal via stacking supercluster and supervoids. The differences in the corresponding ISW cold or hot spots are $sim 20%$ for structures of $sim 100$Mpc/$h$. Such differences are greater for smaller structures, but the amplitude of the signal is lower. The high amplitude of ISW signal detected by Granett et al. can not explained in the $f(R)$ model. We find relatively big differences between $f(R)$ and GR in the transverse bulk motion of matter, and discuss its detectability via the relative frequency shifts of photons from multiple lensed images.

Pan-STARRS1: Galaxy Clustering in the Small Area Survey 2

The Pan-STARRS1 survey is currently obtaining imaging in 5 bands (grizy) for the $3pi$ steradian survey, one of the largest optical surveys ever conducted. The finished survey will have spatially varying depth, due to the survey strategy. This paper presents a method to correct galaxy number counts and galaxy clustering for this potential systematic based on a simplified signal to noise measurement. A star and galaxy separation method calibrated using realistic synthetic images is also presented, along with an approach to mask bright stars. By using our techniques on a ~69 sq. degree region of science verification data this paper shows PS1 measurements of the two point angular correlation function as a function of apparent magnitude agree with measurements from deeper, smaller surveys. Clustering measurements appear reliable down to a magnitude limit of rps<22.5. Additionally, stellar contamination and false detection issues are discussed and quantified. This work is the second of two papers which pave the way for the exploitation of the full $3pi$ survey for studies of large scale structure.

Spatial and luminosity distributions of galactic satellites

We investigate the luminosity functions (LFs) and projected number density profiles of galactic satellites around isolated primaries of different luminosities. We measure these quantities for model satellites placed into the Millennium and Millennium II dark matter simulations by the GALFORM semi-analytic galaxy formation model for different bins of primary galaxy magnitude and we investigate their dependence on satellite luminosity. We compare our model predictions to the data of Guo et al. from the Sloan Digital Sky Survey Data Release 8 (SDSS DR8). First, we use a mock light-cone catalogue to verify that the method we used to count satellites in the SDSS DR8 is unbiased. We find that the radial distributions of model satellites are similar to those around comparable primary galaxies in the SDSS DR8, with only slight differences at low luminosities and small projected radii. However, when splitting the satellites by colour, the model and SDSS satellite systems no longer resemble one another, with many red model satellites, in contrast to the dominant blue fraction at similar luminosity in SDSS. The few model blue satellites are also significantly less centrally concentrated in the halo of their stacked primary than their SDSS counterparts. The implications of this result for the GALFORM model are discussed.

Dark Matter Halo Merger Histories Beyond Cold Dark Matter: I - Methods and Application to Warm Dark Matter

We describe a methodology to accurately compute halo mass functions, progenitor mass functions, merger rates and merger trees in non-cold dark matter universes using a self-consistent treatment of the generalized extended Press-Schechter formalism. Our approach permits rapid exploration of the subhalo population of galactic halos in dark matter models with a variety of different particle properties or universes with rolling, truncated, or more complicated power spectra. We make detailed comparisons of analytically derived mass functions and merger histories with recent warm dark matter cosmological N-body simulations, and find excellent agreement. We show that, once the accretion of smoothly distributed matter is accounted for, coarse-grained statistics such as the mass accretion history of halos can be almost indistinguishable between cold and warm dark matter cases. However, the halo mass function and progenitor mass functions differ significantly, with the warm dark matter cases being strongly suppressed below the free-streaming scale of the dark matter. We demonstrate the importance of using the correct solution for the excursion set barrier first-crossing distribution in warm dark matter - if the solution for a flat barrier is used instead the truncation of the halo mass function is much slower, leading to an overestimate of the number of low mass halos.

Maximum Likelihood Random Galaxy Catalogues and Luminosity Function Estimation

We present a new algorithm to generate a random (unclustered) version of an magnitude limited observational galaxy redshift catalogue. It takes into account both galaxy evolution and the perturbing effects of large scale structure. The key to the algorithm is a maximum likelihood (ML) method for jointly estimating both the luminosity function (LF) and the overdensity as a function of redshift. The random catalogue algorithm then works by cloning each galaxy in the original catalogue, with the number of clones determined by the ML solution. Each of these cloned galaxies is then assigned a random redshift uniformly distributed over the accessible survey volume, taking account of the survey magnitude limit(s) and, optionally, both luminosity and number density evolution. The resulting random catalogues, which can be employed in traditional estimates of galaxy clustering, make fuller use of the information available in the original catalogue and hence are superior to simply fitting a functional form to the observed redshift distribution. They are particularly well suited to studies of the dependence of galaxy clustering on galaxy properties as each galaxy in the random catalogue has the same list of attributes as measured for the galaxies in the genuine catalogue. The derivation of the joint overdensity and LF estimator reveals the limit in which the ML estimate reduces to the standard 1/Vmax LF estimate, namely when one makes the prior assumption that the are no fluctuations in the radial overdensity. The new ML estimator can be viewed as a generalization of the 1/Vmax estimate in which Vmax is replaced by a density corrected Vdc,max.

Which halos host Herschel-ATLAS galaxies in the local Universe?

We measure the projected cross-correlation between low redshift (z < 0.5) far-IR selected galaxies in the SDP field of the Herschel-ATLAS (H-ATLAS) survey and optically selected galaxies from the Galaxy and Mass Assembly (GAMA) redshift survey. In order to obtain robust correlation functions, we restrict the analysis to a subset of 969 out of 6900 H-ATLAS galaxies, which have reliable optical counterparts with r<19.4 mag and well-determined spectroscopic redshifts. The overlap region between the two surveys is 12.6 sq. deg; the matched sample has a median redshift of z ~ 0.2. The cross-correlation of GAMA and H-ATLAS galaxies within this region can be fitted by a power law, with correlation length r_0 ~ 4.63 +/- 0.51 Mpc. Comparing with the corresponding auto-correlation function of GAMA galaxies within the SDP field yields a relative bias (averaged over 2-8 Mpc) of H-ATLAS and GAMA galaxies of b_H/b_G ~ 0.6. Combined with clustering measurements from previous optical studies, this indicates that most of the low redshift H-ATLAS sources are hosted by halos with masses comparable to that of the Milky Way. The correlation function appears to depend on the 250 um luminosity, L_250, with bright (median luminosity u L_250 ~ 1.6 x 10^10 L_sun) objects being somewhat more strongly clustered than faint ( u L_250 ~ 4.0 x 10^9 L_sun) objects. This implies that galaxies with higher dust-obscured star formation rates are hosted by more massive halos.

The galaxy power spectrum: precision cosmology from large scale structure?

Published galaxy power spectra from the 2dFGRS and SDSS are not in good agreement. We revisit this issue by analyzing both the 2dFGRS and SDSS DR5 catalogues using essentially identical techniques. We confirm that the 2dFGRS exhibits relatively more large scale power than the SDSS, or, equivalently, SDSS has more small scale power. We demonstrate that this difference is due to the r-band selected SDSS catalogue being dominated by more strongly clustered red galaxies, which have a stronger scale dependent bias. The power spectra of galaxies of the same rest frame colours from the two surveys match well. If not accounted for, the difference between the SDSS and 2dFGRS power spectra causes a bias in the obtained constraints on cosmological parameters which is larger than the uncertainty with which they are determined. We also found that the correction developed by Cole et al.(2005) to model the distortion in the shape of the power spectrum due to non-linear evolution and scale dependent bias is not able to reconcile the constraints obtained from the 2dFGRS and SDSS power spectra. Intriguingly, the model is able to describe the differences between the 2dFGRS and the much more strongly clustered LRG sample, which exhibits greater nonlinearities. This shows that more work is needed to understand the relation between the galaxy power spectrum and the linear perturbation theory prediction for the power spectrum of matter fluctuations. It is therefore important to accurately model these effects to get precise estimates of cosmological parameters from these power spectra and from future galaxy surveys like Pan-STARRS, or the Dark Energy Survey, which will use selection criteria similar to the one of SDSS.

Generating Dark Matter Halo Merger Trees

We present a new Monte-Carlo algorithm to generate merger trees describing the formation history of dark matter halos. The algorithm is a modification of the algorithm of Cole et al (2000) used in the GALFORM semi-analytic galaxy formation model. As such, it is based on the Extended Press-Schechter theory and so should be applicable to hierarchical models with a wide range of power spectra and cosmological models. It is tuned to be in accurate agreement with the conditional mass functions found in the analysis of merger trees extracted from the LCDM Millennium N-body simulation. We present a comparison of its predictions not only with these conditional mass functions, but also with additional statistics of the Millennium Simulation halo merger histories. In all cases we find it to be in good agreement with the Millennium Simulation and thus it should prove to be a very useful tool for semi-analytic models of galaxy formation and for modelling hierarchical structure formation in general. We have made our merger tree generation code and code to navigate the trees available at http://star-www.dur.ac.uk/~cole/merger_trees .

Measuring the Baryon Acoustic Oscillation scale using the SDSS and 2dFGRS

We introduce a method to constrain general cosmological models using Baryon Acoustic Oscillation (BAO) distance measurements from galaxy samples covering different redshift ranges, and apply this method to analyse samples drawn from the SDSS and 2dFGRS. BAO are detected in the clustering of the combined 2dFGRS and SDSS main galaxy samples, and measure the distance--redshift relation at z=0.2. BAO in the clustering of the SDSS luminous red galaxies measure the distance--redshift relation at z=0.35. The observed scale of the BAO calculated from these samples and from the combined sample are jointly analysed using estimates of the correlated errors, to constrain the form of the distance measure D_V(z)=[(1+z)^2D_A^2cz/H(z)]^(1/3). Here D_A is the angular diameter distance, and H(z) is the Hubble parameter. This gives r_s/D_V(0.2)=0.1980+/-0.0058 and r_s/D_V(0.35)=0.1094+/-0.0033 (1sigma errors), with correlation coefficient of 0.39, where r_s is the comoving sound horizon scale at recombination. Matching the BAO to have the same measured scale at all redshifts then gives D_V(0.35)/D_V(0.2)=1.812+/-0.060. The recovered ratio is roughly consistent with that predicted by the higher redshift SNLS supernovae data for Lambda cosmologies, but does require slightly stronger cosmological acceleration at low redshift. If we force the cosmological model to be flat with constant w, then we find Om_m=0.249+/-0.018 and w=-1.004+/-0.089 after combining with the SNLS data, and including the WMAP measurement of the apparent acoustic horizon angle in the CMB.

The statistical properties of LCDM halo formation

We present a comparison of the statistical properties of dark matter halo merger trees extracted from the Millennium Simulation with Extended Press-Schechter (EPS) formalism and the related GALFORM Monte-Carlo method for generating ensembles of merger trees. The volume, mass resolution and output frequency make the Millennium Simulation a unique resource for the study of the hierarchical growth of structure. We construct the merger trees of present day friends-of-friends groups and calculate a variety of statistics that quantify the masses of their progenitors as a function of redshift; accretion rates; and the redshift distribution of their most recent major merger. We also look in the forward direction and quantify the present day mass distribution of halos into which high redshift progenitors of a specific mass become incorporated. We find that EPS formalism and its Monte-Carlo extension capture the qualitative behaviour of all these statistics but, as redshift increases they systematically underestimate the masses of the most massive progenitors. This shortcoming is worst for the Monte-Carlo algorithm. We present a fitting function to a scaled version of the progenitor mass distribution and show how it can be used to make more accurate predictions of both progenitor and final halo mass distributions.