BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

C. M. Baugh

C. M. Baugh appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.COastro-phastro-ph.GAhep-phhep-th

Catalog footprint

What is connected

81works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

81 published item(s)

preprint2026arXiv

The PAU Survey: Measuring intrinsic galaxy alignments in deep wide fields as a function of colour, luminosity, stellar mass and redshift

We present the measurements and constraints of intrinsic alignments (IA) in the Physics of the Accelerating Universe Survey (PAUS) deep wide fields, which include the W1 and W3 fields from the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) and the G09 field from the Kilo-Degree Survey (KiDS). Our analyses cover 51deg$^{2}$, in the photometric redshift (photo-$z$) range $0.1 < z_{\mathrm{b}} < 1$ and a magnitude limit $i_{\mathrm{AB}}<22$. The precise photo-$z$s and the luminosity coverage of PAUS enable robust IA measurements, which are key for setting informative priors for upcoming stage-IV surveys. For red galaxies, we detect an increase in IA amplitude with both luminosity and stellar mass, extending previous results towards fainter and less massive regimes. As a function of redshift, we observe strong IA signals at intermediate ($z_{\mathrm{b}}\sim0.55$) and high ($z_{\mathrm{b}}\sim0.75$) redshift bins. However, we find no significant trend of IA evolution with redshift after accounting for the varying luminosities across redshift bins, consistent with the literature. For blue galaxies, no significant IA signal is detected, with $A_{1}=0.68_{-0.51}^{+0.53}$ when splitting only by galaxy colour, yielding some of the tightest constraints to date for the blue population and constraining a regime of very faint and low-mass galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Modelling emission lines in star forming galaxies

We present a new model to compute the luminosity of emission lines in star forming galaxies and apply this in the semi-analytical galaxy formation code GALFORM. The model combines a pre-computed grid of HII region models with an empirical determination of how the properties of HII regions depend on the macroscopic properties of galaxies based on observations of local galaxies. The new model gives a very good reproduction of the locus of star-forming galaxies on standard line ratio diagnostic diagrams. The new model shows evolution in the locus of star forming galaxies with redshift on this line ratio diagram, with a good match to the observed line ratios at $z=1.6$. The model galaxies at high redshift have gas densities and ionisation parameters that are predicted to be $\approx 2-3$ times higher than in local star forming galaxies, which is partly driven by the changing selection with redshift to mimic the observational selection. Our results suggest that the observed evolution in emission line ratios requires other HII region properties to evolve with redshift, such as the gas density, and cannot be reproduced by HII model grids that only allow the gas metallicity and ionisation parameter to vary.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Modelling the quenching of star formation activity from the evolution of the colour-magnitude relation in VIPERS

We study the evolution of the colour-magnitude relation for galaxies in the VIMOS Public Extragalactic Redshift Survey (VIPERS) by introducing the concept of the bright edge, and use this to derive constraints on the quenching of star formation activity in galaxies over the redshift range $0.5 < z < 1.1$. The bright-edge of the colour-magnitude diagram evolves with little dependence on galaxy colour, and therefore on the amount of star formation taking place in bright galaxies. We modelled this evolution with delayed exponential star formation histories (SFHs), to better understand the time-scale of the turn-off in star formation activity. We show that using SFHs without quenching, the transition from the blue cloud to the red sequence is too slow. This indicates that a scenario purely driven by the consumption of the gas inside each galaxy does not reproduce the observed evolution of the colour-magnitude bright edge. Among the quenching scenarios explored, the one that best matches the observations assumes that galaxies stop their star formation at a randomly selected time with a uniform distribution up to $2.5$ Gyr. We argue that quenching is required over a wide range of stellar masses. Qualitatively similar evolution of the bright edge is found in the predictions of a semi-analytical galaxy formation model, but quantitatively there are marked differences with the observations. This illustrates the utility of the bright edge as a test of galaxy formation models. The evolution changes and no longer matches the observed trend if feedback from heating by active galactic nuclei is turned off.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dark matter-radiation interactions: the structure of Milky Way satellite galaxies

In the thermal dark matter (DM) paradigm, primordial interactions between DM and Standard Model particles are responsible for the observed DM relic density. In Boehm et al. (2014), we showed that weak-strength interactions between DM and radiation (photons or neutrinos) can erase small-scale density fluctuations, leading to a suppression of the matter power spectrum compared to the collisionless cold DM (CDM) model. This results in fewer DM subhaloes within Milky Way-like DM haloes, implying a reduction in the abundance of satellite galaxies. Here we use very high resolution N-body simulations to measure the dynamics of these subhaloes. We find that when interactions are included, the largest subhaloes are less concentrated than their counterparts in the collisionless CDM model and have rotation curves that match observational data, providing a new solution to the "too big to fail" problem.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Measuring galaxy environment with the synergy of future photometric and spectroscopic surveys

We exploit the synergy between low-resolution spectroscopy and photometric redshifts to study environmental effects on galaxy evolution in slitless spectroscopic surveys from space. As a test case, we consider the future Euclid Deep survey (~40deg$^2$), which combines a slitless spectroscopic survey limited at H$α$ flux $\geq5\times 10^{-17}$ erg cm$^{-2}$ s$^{-1}$ and a photometric survey limited in H-band ($H\leq26$). We use Euclid-like galaxy mock catalogues, in which we anchor the photometric redshifts to the 3D galaxy distribution of the available spectroscopic redshifts. We then estimate the local density contrast by counting objects in cylindrical cells with radius from 1 to 10 h$^{-1}$Mpc over the redshift range 0.9<z<1.8. We compare this density field with the one computed in a mock catalogue with the same depth as the Euclid Deep survey (H=26) but without redshift measurement errors. We find that our method successfully separates high from low density environments (the last from the first quintile of the density distribution), with higher efficiency at low redshift and large cell: the fraction of low density regions mistaken by high density peaks is <1% for all scales and redshifts explored, but for scales of 1 h$^{-1}$Mpc for which is a few percent. These results show that we can efficiently study environment in photometric samples if spectroscopic information is available for a smaller sample of objects that sparsely samples the same volume. We demonstrate that these studies are possible in the Euclid Deep survey, i.e. in a redshift range in which environmental effects are different from those observed in the local universe, hence providing new constraints for galaxy evolution models.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The clustering and halo occupation distribution of Lyman-break galaxies at $z\sim4$

We investigate the clustering of Lyman-break galaxies (LBGs) at $z\sim4$. Using the hierarchical galaxy formation model GALFORM, we predict, for the first time using a semi-analytical model with feedback from active galactic nuclei (AGN), the angular correlation function (ACF) of LBGs and find agreement within $3\,σ$ with new measurements of the ACF from surveys including the Hubble eXtreme Deep Field (XDF) and CANDELS field. Our simulations confirm the conclusion reached using independent models that although the predicted ACFs reproduce the trend of increased clustering with luminosity, the dependence is less strong than observed. We find that for the detection limits of the XDF field central LBGs at $z\sim 4$ predominantly reside in haloes of mass $\sim 10^{11}-10^{12}h^{-1}M_{\rm \odot}$ and that satellites reside in larger haloes of mass $\sim 10^{12}-10^{13}h^{-1}M_{\rm \odot}$. The model predicts fewer bright satellite LBGs at $z\sim4$ than is inferred from measurements of the ACF at small scales. By analysing the halo occupation distribution (HOD) predicted by the model, we find evidence that AGN feedback affects the HOD of central LBGs in massive haloes. This is a new high-redshift test of this important feedback mechanism. We investigate the effect of photometric errors in the observations on the ACF predictions. We find that the observational uncertainty in the galaxy luminosity reduces the clustering amplitude and that this effect increases towards faint galaxies, particularly on small scales. To compare properties of model with observed LBGs this uncertainty must be considered.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The evolution of the galaxy content of dark matter haloes

We use the halo occupation distribution (HOD) framework to characterise the predictions from two independent galaxy formation models for the galactic content of dark matter haloes and its evolution with redshift. Our galaxy samples correspond to a range of fixed number densities defined by stellar mass and span $0 \le z \le 3$. We find remarkable similarities between the model predictions. Differences arise at low galaxy number densities which are sensitive to the treatment of heating of the hot halo by active galactic nuclei. The evolution of the form of the HOD can be described in a relatively simple way, and we model each HOD parameter using its value at $z=0$ and an additional evolutionary parameter. In particular, we find that the ratio between the characteristic halo masses for hosting central and satellite galaxies can serve as a sensitive diagnostic for galaxy evolution models. Our results can be used to test and develop empirical studies of galaxy evolution and can facilitate the construction of mock galaxy catalogues for future surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Extraordinary Amount of Substructure in the Hubble Frontier Fields Cluster Abell 2744

We present a joint optical/X-ray analysis of the massive galaxy cluster Abell 2744 (z=0.308). Our strong- and weak-lensing analysis within the central region of the cluster, i.e., at R<1Mpc from the brightest cluster galaxy, reveals eight substructures, including the main core. All of these dark-matter halos are detected with a significance of at least 5sigma and feature masses ranging from 0.5 to 1.4x10^{14}Msun within R<150kpc. Merten et al. (2011) and Medezinski et al. (2016) substructures are also detected by us. We measure a slightly higher mass for the main core component than reported previously and attribute the discrepancy to the inclusion of our tightly constrained strong-lensing mass model built on Hubble Frontier Fields data. X-ray data obtained by XMM-Newton reveal four remnant cores, one of them a new detection, and three shocks. Unlike Merten et al. (2011), we find all cores to have both dark and luminous counterparts. A comparison with clusters of similar mass in the MXXL simulations yields no objects with as many massive substructures as observed in Abell 2744, confirming that Abell 2744 is an extreme system. We stress that these properties still do not constitute a challenge to $Λ$CDM, as caveats apply to both the simulation and the observations: for instance, the projected mass measurements from gravitational lensing and the limited resolution of the sub-haloes finders. We discuss implications of Abell 2744 for the plausibility of different dark-matter candidates and, finally, measure a new upper limit on the self-interaction cross-section of dark matter of sigma_{DM}<1.28cm2/g(68\% CL), in good agreement with previous results from Harvey et al. (2015).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The spatial distribution of neutral hydrogen as traced by low HI mass galaxies

The formation and evolution of galaxies with low neutral atomic hydrogen (HI) masses, M$_{\rm HI}$$<$10$^{8}h^{-2}$M$_{\odot}$, are affected by host dark matter halo mass and photoionisation feedback from the UV background after the end of reionization. We study how the physical processes governing the formation of galaxies with low HI mass are imprinted on the distribution of neutral hydrogen in the Universe using the hierarchical galaxy formation model, GALFORM. We calculate the effect on the correlation function of changing the HI mass detection threshold at redshifts $0 \le z \le 0.5$. We parameterize the clustering as $ξ(r)=(r/r_{0})^{-γ}$ and we find that including galaxies with M$_{\rm HI}$$<$10$^{8}h^{-2}$M$_{\odot}$ increases the clustering amplitude $r_{0}$ and slope $γ$ compared to samples of higher HI masses. This is due to these galaxies with low HI masses typically being hosted by haloes with masses greater than 10$^{12}{h}^{-1}$M$_{\odot}$, and is in contrast to optically selected surveys for which the inclusion of faint, blue galaxies lowers the clustering amplitude. We show the HI mass function for different host dark matter halo masses and galaxy types (central or satellite) to interpret the values of $r_{0}$ and $γ$ of the clustering of HI-selected galaxies. We also predict the contribution of low HI mass galaxies to the 21cm intensity mapping signal. We calculate that a dark matter halo mass resolution better than $\sim$10$^{10}{h}^{-1}$M$_{\odot}$ at redshifts higher than 0.5 is required in order to predict converged 21cm brightness temperature fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A hybrid multi resolution scheme to efficiently model the structure of reionization on the largest scales

Redshifted 21cm measurements of the structure of ionised regions that grow during reionization promise to provide a new probe of early galaxy and structure formation. One of the challenges of modelling reionization is to account both for the sub-halo scale physics of galaxy formation and the regions of ionization on scales that are many orders of magnitude larger. To bridge this gap we first calculate the statistical relationship between ionizing luminosity and Mpc-scale overdensity using detailed models of galaxy formation computed using relatively small volume - ($\sim$100Mpc/$h$)$^{3}$, high resolution dark matter simulations. We then use a Monte-Carlo technique to apply this relationship to reionization of the intergalactic medium within large volume dark matter simulations - ($>$1Gpc/$h$)$^{3}$. The resulting simulations can be used to address the contribution of very large scale clustering of galaxies to the structure of reionization, and show that volumes larger than 500Mpc/$h$ are required to probe the largest reionization features mid-way through reionization. As an example application of our technique, we demonstrate that the predicted 21cm power spectrum amplitude and gradient could be used to determine the importance of supernovae feedback for early galaxy formation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dark matter-radiation interactions: the impact on dark matter haloes

Interactions between dark matter (DM) and radiation (photons or neutrinos) in the early Universe suppress density fluctuations on small mass scales. Here we perform a thorough analysis of structure formation in the fully non-linear regime using N-body simulations for models with DM-radiation interactions and compare the results to a traditional calculation in which DM only interacts gravitationally. Significant differences arise due to the presence of interactions, in terms of the number of low-mass DM haloes and their properties, such as their spin and density profile. These differences are clearly seen even for haloes more massive than the scale on which density fluctuations are suppressed. We also show that semi-analytical descriptions of the matter distribution in the non-linear regime fail to reproduce our numerical results, emphasizing the challenge of predicting structure formation in models with physics beyond collisionless DM.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Galaxy Formation & Dark Matter Modelling in the Era of the Square Kilometre Array

Theoretical galaxy formation models are an established and powerful tool for interpreting the astrophysical significance of observational data, particularly galaxy surveys. Such models have been utilised with great success by optical surveys such as 2dFGRS and SDSS, but their application to radio surveys of cold gas in galaxies has been limited. In this chapter we describe recent developments in the modelling of the cold gas properties in the models, and how these developments are essential if they are to be applied to cold gas surveys of the kind that will be carried out with the SKA. By linking explicitly a galaxy's star formation rate to the abundance of molecular hydrogen in the galaxy rather than cold gas abundance, as was assumed previously, the latest models reproduce naturally many of the global atomic and molecular hydrogen properties of observed galaxies. We review some of the key results of the latest models and highlight areas where further developments are necessary. We discuss also how model predictions can be most accurately compared with observational data, what challenges we expect when creating synthetic galaxy surveys in the SKA era, and how the SKA can be used to test models of dark matter.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The galaxy - dark matter halo connection: which galaxy properties are correlated with the host halo mass?

We demonstrate how the properties of a galaxy depend on the mass of its host dark matter subhalo, using two independent models of galaxy formation. For the cases of stellar mass and black hole mass, the median property value displays a monotonic dependence on subhalo mass. The slope of the relation changes for subhalo masses for which heating by active galactic nuclei becomes important. The median property values are predicted to be remarkably similar for central and satellite galaxies. The two models predict considerable scatter around the median property value, though the size of the scatter is model dependent. There is only modest evolution with redshift in the median galaxy property at a fixed subhalo mass. Properties such as cold gas mass and star formation rate, however, are predicted to have a complex dependence on subhalo mass. In these cases subhalo mass is not a good indicator of the value of the galaxy property. We illustrate how the predictions in the galaxy property - subhalo mass plane differ from the assumptions made in empirical models of galaxy clustering by reconstructing the model output using a subhalo abundance matching scheme. In its simplest form, abundance matching generally does not reproduce the clustering predicted by the models, typically resulting in an overprediction of the clustering signal. We show how the basic abundance matching scheme can be extended to reproduce the model predictions more faithfully, which has implications for the analysis of galaxy clustering, particularly for low abundance samples.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The HI mass function as a probe of photoionisation feedback on low mass galaxy formation

We explore the galaxy formation physics governing the low mass end of the HI mass function in the local Universe. Specifically, we predict the effects on the HI mass function of varying i) the strength of photoionisation feedback and the redshift of the end of the epoch of reionization, ii) the cosmology, iii) the supernovae feedback prescription, and iv) the efficiency of star formation. We find that the shape of the low-mass end of the HI mass function is most affected by the critical halo mass below which galaxy formation is suppressed by photoionisation heating of the intergalactic medium. We model the redshift dependence of this critical dark matter halo mass by requiring a match to the low-mass end of the HI mass function. The best fitting critical dark matter halo mass decreases as redshift increases in this model, corresponding to a circular velocity of $\sim 50 \, {\rm km \,s}^{-1}$ at $z=0$, $\sim 30 \, {\rm km\, s}^{-1}$ at $z \sim 1$ and $\sim 12 \, {\rm km \, s}^{-1}$ at $z=6$. We find that an evolving critical halo mass is required to explain both the shape and abundance of galaxies in the HI mass function below $M_{\rm HI} \sim 10^{8} h^{-2} {\rm M_{\odot}}$. The model makes specific predictions for the clustering strength of HI-selected galaxies with HI masses > $10^{6} h^{-2} {\rm M_{\odot}}$ and $> 10^{7} h^{-2} {\rm M_{\odot}}$ and for the relation between the HI and stellar mass contents of galaxies which will be testable with upcoming surveys with the Square Kilometre Array and its pathfinders. We conclude that measurements of the HI mass function at $z \ge 0$ will lead to an improvement in our understanding of the net effect of photoionisation feedback on galaxy formation and evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The origin of the atomic and molecular gas contents of early-type galaxies. II. Misaligned gas accretion

We study the origin of the wide distribution of angles between the angular momenta of the stellar and gas components, $α_{\rm G,S}$, in early-type galaxies (ETGs). We use the GALFORM model of galaxy formation, set in the $Λ$ cold dark matter framework, and coupled it with a Monte-Carlo simulation to follow the angular momenta flips driven by matter accretion onto haloes and galaxies. We consider a gas disk to be misaligned with respect to the stellar body if $α_{\rm G,S}>30$~degrees. By assuming that the only sources of misaligments in galaxies are galaxy mergers, we place a lower limit of $2-5$ per cent on the fraction of ETGs with misaligned gas/stellar components. These low fractions are inconsistent with the observed value of $\approx 42\pm 6$ per cent in ATLAS$^{\rm 3D}$. In the more general case, in which smooth gas accretion in addition to galaxy mergers can drive misalignments, our calculation predicts that $\approx 46$ per cent of ETGs have $α_{\rm G,S}>30$~degrees. In this calculation, we find correlations between $α_{\rm G,S}$ and stellar mass, cold gas fraction and star formation rate, such that ETGs with high masses, low cold gas fractions and low star formation rates are more likely to display aligned cold gas and stellar components. We confirm these trends observationally for the first time using ATLAS$^{\rm 3D}$ data. We argue that the high fraction of misaligned gas discs observed in ETGs is mostly due to smooth gas accretion (e.g. cooling from the hot halo of galaxies) which takes place after most of the stellar mass of the galaxy is in place and comes misaligned with respect to the stellar component. Galaxies that have accreted most of their cold gas content prior to the time where most of the stellar mass was in place show aligned components.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A new model of galaxy formation: How sensitive are predicted galaxy luminosities to the choice of SPS model?

We present a new release of the GALFORM semi-analytical model of galaxy formation and evolution, which exploits a Millennium Simulation-class N-body run performed with the WMAP7 cosmology. We use this new model to study the impact of the choice of stellar population synthesis (SPS) model on the predicted evolution of the galaxy luminosity function. The semi-analytical model is run using seven different SPS models. In each case we obtain the rest-frame luminosity function in the far-ultra-violet, optical and near-infrared (NIR) wavelength ranges. We find that both the predicted rest-frame ultra-violet and optical luminosity function are insensitive to the choice of SPS model. However, we find that the predicted evolution of the rest-frame NIR luminosity function depends strongly on the treatment of the thermally pulsating asymptotic giant branch (TP-AGB) stellar phase in the SPS models, with differences larger than a factor of 2 for model galaxies brighter than $M_{\rm AB}(K)-5$log$h<-22$ ($\sim$L$_*$ for $0\leq z\leq 1.5$). We have also explored the predicted number counts of galaxies, finding remarkable agreement between the results with different choices of SPS model, except when selecting galaxies with very red optical-NIR colours. The predicted number counts of these extremely red galaxies appear to be more affected by the treatment of star formation in disks than by the treatment of TP-AGB stars in the SPS models.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The origin of the atomic and molecular gas contents of early-type galaxies. I. A new test of galaxy formation physics

We study the atomic (HI) and molecular hydrogen (H2) contents of early-type galaxies (ETGs) and their gas sources using the GALFORM model of galaxy formation. This model uses a self-consistent calculation of the star formation rate (SFR), which depends on the H2 content of galaxies. We first present a new analysis of HIPASS and ATLAS3D surveys, with special emphasis on ETGs. The model predicts HI and H2 contents of ETGs in agreement with the observations from these surveys only if partial ram pressure stripping of the hot gas is included, showing that observations of neutral gas in `quenched' galaxies place stringent constraints on the treatment of the hot gas in satellites. We find that ~90% of ETGs at z=0 have neutral gas contents supplied by radiative cooling from their hot halos, 8% were supplied by gas accretion from minor mergers that took place in the last 1Gyr, while 2% were supplied by mass loss from old stars. The model predicts neutral gas fractions strongly decreasing with increasing bulge fraction. This is due to the impeded disk regeneration in ETGs, resulting from both active galactic nuclei feedback and environmental quenching by partial ram pressure stripping of the hot gas.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Using the Milky Way satellites to study interactions between cold dark matter and radiation

The cold dark matter (CDM) model faces persistent challenges on small scales. In particular, taken at face value, the model significantly overestimates the number of satellite galaxies around the Milky Way. Attempts to solve this problem remain open to debate and have even led some to abandon CDM altogether. However, current simulations are limited by the assumption that dark matter feels only gravity. Here, we show that including interactions between CDM and radiation (photons or neutrinos) leads to a dramatic reduction in the number of satellite galaxies, alleviating the Milky Way satellite problem and indicating that physics beyond gravity may be essential to make accurate predictions of structure formation on small scales. The methodology introduced here gives constraints on dark matter interactions that are significantly improved over those from the cosmic microwave background.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Which galaxies dominate the neutral gas content of the Universe?

We study the contribution of galaxies with different properties to the global densities of star formation rate (SFR), atomic (HI) and molecular hydrogen (H2) as a function of redshift. We use the GALFORM model of galaxy formation, which is set in the LCDM framework. This model includes a self-consistent calculation of the SFR, which depends on the H2 content of galaxies. The predicted SFR density and how much of this is contributed by galaxies with different stellar masses and infrared luminosities are in agreement with observations. The model predicts a modest evolution of the HI density at z<3, which is also in agreement with the observations. The HI density is predicted to be always dominated by galaxies with SFR<1Msun/yr. This contrasts with the H2 density, which is predicted to be dominated by galaxies with SFR>10Msun/yr. Current high-redshift galaxy surveys are limited to detect carbon monoxide in galaxies with SFR>30Msun/yr, which in our model make up, at most, 20% of the H2 in the universe. In terms of stellar mass, the predicted H2 density is dominated by massive galaxies, Mstellar>10^10Msun, while the HI density is dominated by low mass galaxies, Mstellar<10^9Msun. In the context of upcoming neutral gas surveys, we suggest that the faint nature of the galaxies dominating the HI content of the Universe will hamper the identification of optical counterparts, while for H2, we expect follow up observations of molecular emission lines of already existing galaxy catalogues to be able to uncover the H2 density of the Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Constraints on black hole fuelling modes from the clustering of X-ray AGN

We present a clustering analysis of X-ray selected AGN by compiling X-ray samples from the literature and re-estimating the dark matter (DM) halo masses of AGN in a uniform manner. We find that moderate luminosity AGN (Lx(2-10 keV)=10^42-10^44 erg/sec) in the z=0-1.3 Universe are typically found in DM haloes with masses of ~10^13 Msun. We then compare our findings to the theoretical predictions of the coupled galaxy and black hole formation model GALFORM. We find good agreement when our calculation includes the hot-halo mode of accretion onto the central black hole. This type of accretion, which is additional to the common cold accretion during disk instabilities and galaxy mergers, is tightly coupled to the AGN feedback in the model. The hot-halo mode becomes prominent in DM haloes with masses greater than ~10^12.5 Msun, where AGN feedback typically operates, giving rise to a distinct class of moderate luminosity AGN that inhabit rich clusters and superclusters. Cold gas fuelling of the black hole cannot produce the observationally inferred DM halo masses of X-ray AGN. Switching off AGN feedback in the model results in a large population of luminous quasars (Lx(2-10 keV) > 10^44 erg/sec) in DM haloes with masses up to ~10^14 Msun, which is inconsistent with the observed clustering of quasars. The abundance of hot-halo AGN decreases significantly in the z~3-4 universe. At such high redshifts, the cold accretion mode is solely responsible for shaping the environment of moderate luminosity AGN. Our analysis supports two accretion modes (cold and hot) for the fuelling of supermassive black holes and strongly underlines the importance of AGN feedback in cosmological models both of galaxy formation and black hole growth.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

How well can we really estimate the stellar masses of galaxies from broadband photometry?

The estimated stellar masses of galaxies are widely used to characterize how the galaxy population evolves over cosmic time. If stellar masses can be estimated in a robust manner, free from any bias, global diagnostics such as the stellar mass function can be used to constrain the physics of galaxy formation. We explore how galaxy stellar masses, estimated by fitting broad-band spectral energy distributions (SEDs) with stellar population models, can be biased as a result of commonly adopted assumptions for the star-formation and chemical enrichment histories, recycled fractions and dust attenuation curves of galaxies. We apply the observational technique of broad-band SED fitting to model galaxy SEDs calculated by the theoretical galaxy formation model GALFORM, isolating the effect of each of these assumptions. We find that, averaged over the entire galaxy population, the common assumption of exponentially declining star-formation histories does not adversely affect stellar mass estimation. We show that fixing the metallicity in SED fitting or using sparsely sampled metallicity grids can introduce mass dependent systematics into stellar mass estimates. We find that the common assumption of a star-dust geometry corresponding to a uniform foreground dust screen can cause the stellar masses of dusty model galaxies to be significantly underestimated. Finally, we show that stellar mass functions recovered by applying SED fitting to model galaxies at high redshift can differ significantly in both shape and normalization from the intrinsic mass functions predicted by a given model. Given these differences, our methodology of using stellar masses estimated from model galaxy SEDs offers a new, self-consistent way to compare model predictions with observations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Luminosity bias: from haloes to galaxies

Large surveys of the local Universe have shown that galaxies with different intrinsic properties, such as colour, luminosity and morphological type display a range of clustering amplitudes. Galaxies are therefore not faithful tracers of the underlying matter distribution. This modulation of galaxy clustering, called bias, contains information about the physics behind galaxy formation. It is also a systematic to be overcome before the large-scale structure of the Universe can be used as a cosmological probe. Two types of approaches have been developed to model the clustering of galaxies. The first class is empirical and filters or weights the distribution of dark matter to reproduce the measured clustering. In the second approach an attempt is made to model the physics which governs fate of baryons in order to predict the number of galaxies in dark matter haloes. I will review the development of both approaches and summarize what we have learnt about galaxy bias.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Scaling relations for galaxy clusters in the Millennium-XXL simulation

We present a very large high-resolution cosmological N-body simulation, the Millennium-XXL or MXXL, which uses 303 billion particles to represent the formation of dark matter structures throughout a 4.1Gpc box in a LambdaCDM cosmology. We create sky maps and identify large samples of galaxy clusters using surrogates for four different observables: richness estimated from galaxy surveys, X-ray luminosity, integrated Sunyaev-Zeldovich signal, and lensing mass. The unprecedented combination of volume and resolution allows us to explore in detail how these observables scale with each other and with cluster mass. The scatter correlates between different mass-observable relations because of common sensitivities to the internal structure, orientation and environment of clusters, as well as to line-of-sight superposition of uncorrelated structure. We show that this can account for the apparent discrepancies uncovered recently between the mean thermal SZ signals measured for optically and X-ray selected clusters by stacking data from the Planck satellite. Related systematics can also affect inferences from extreme clusters detected at high redshift. Our results illustrate that cosmological conclusions from galaxy cluster surveys depend critically on proper modelling, not only of the relevant physics, but also of the full distribution of the observables and of the selection biases induced by cluster identification procedures.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The journey of QSO haloes from z=6 to the present

We apply a recently developed scaling technique to the Millennium-XXL, one of the largest cosmological N-body simulations carried out to date 3x10^11 particles within a cube of volume ~70Gpc^3). This allows us to investigate the cosmological parameter dependence of the mass and evolution of haloes in the extreme high-mass tail of the z=6 distribution. We assume these objects to be likely hosts for the population of rare but ultraluminous high-redshift quasars discovered by the Sloan Digital Sky Survey. Haloes with a similar abundance to these quasars have a median mass of 9x10^12 Msun in the currently preferred cosmology, but do not evolve into equally extreme objects at z=0. Rather, their descendants span the full range conventionally assigned to present-day clusters, 6x10^13 to 2.5x10^15 Msun for this same cosmology. The masses both at z=6 and at z=0 shift up or down by factors exceeding two if cosmological parameters are pushed to the boundaries of the range discussed in published interpretations of data from the WMAP satellite. The main factor determining the future growth of a high-mass z=6 halo is the mean overdensity of its environment on scales of 7 to 14 Mpc, and descendant masses can be predicted 6 to 8 times more accurately if this density is known than if it is not. All these features are not unique to extreme high-z haloes, but are generic to hierarchical growth. Finally, we find that extreme haloes at z=6 typically acquired about half of their total mass in the preceding 100 Myr, implying very large recent accretion rates which may be related to the large black hole masses and high luminosities of the SDSS quasars.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The most luminous quasars do not live in the most massive dark matter haloes at any redshift

Quasars represent the brightest Active Galactic Nuclei (AGN) in the Universe and are thought to indicate the location of prodigiously growing Black Holes (BHs), with luminosities as high as 10^48 erg/sec. It is often expected though that such an extremely energetic process will take place in the most massive bound structures in the dark matter (DM) distribution. We show that in contrast to this expectation, in a galaxy formation model which includes AGN feedback, quasars are predicted to live in average DM halo environments with typical masses of a few times 10^12 Msun. This fundamental prediction arises from the fact that quasar activity (i.e., BH accretion with luminosity greater than 10^46 erg/sec) is inhibited in DM haloes where AGN feedback operates. The galaxy hosts of quasars in our simulations are identified with over massive (in gas and stars) spheroidal galaxies, in which BH accretion is triggered via a galaxy merger or secular processes. We further show that the z=0 descendants of high redshift (z~6) QSOs span a wide range of morphologies, galaxy and halo masses. The z~6 BHs typically grow only by a modest factor by the present day. Remarkably, high redshift QSOs never inhabit the largest DM haloes at that time and their descendants are very seldom found in the most massive haloes at z=0. We also show that observationally it is very likely to find an enhancement in the abundance of galaxies around quasars at z~5. However, these enhancements are considerably weaker compared to the overdensities expected at the extreme peaks of the DM distribution. Thus, it is very unlikely that a quasar detected in the $z\gtrsim5$ Universe pinpoints the location of the progenitors of superclusters in the local Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Clustering properties of BzK-selected galaxies in GOODS-N: environmental quenching and triggering of star formation at z ~ 2

Using a sample of BzK-selected galaxies at z~2 identified from the CFHT/WIRCAM near-infrared survey of GOODS-North, we discuss the relation between star formation rate (SFR), specific star formation rate (SSFR), and stellar mass (M_{*}), and the clustering of galaxies as a function of these parameters. For star-forming galaxies (sBzKs), the UV-based SFR, corrected for extinction, scales with the stellar mass as SFR ~ M_{*}^{alpha} with alpha = 0.74+/-0.20 down to M_{*} ~ 10^{9} M_{solar}, indicating a weak dependence on the stellar mass of the SSFR. We also measure the angular correlation function and hence infer the correlation length for sBzK galaxies as a function of M_{*}, SFR, and SSFR, as well as K-band apparent magnitude. We show that passive galaxies (pBzKs) are more strongly clustered than sBzK galaxies at a given stellar mass, mirroring the color-density relation seen at lower redshifts. We also find that the correlation length of sBzK galaxies ranges from 4 to 20 h^{-1}Mpc, being a strong function of M_{K}, M_{*}, and SFR. On the other hand, the clustering dependence on SSFR changes abruptly at 2x10^{-9} yr^{-1}, which is the typical value for "main sequence" star-forming galaxies at z~2. We show that the correlation length reaches a minimum at this characteristic value, and is larger for galaxies with both smaller and larger SSFRs; a dichotomy that is only marginally implied from the predictions of the semi-analytical models. Our results suggest that there are two types of environmental effects at work at z~2. Stronger clustering for relatively quiescent galaxies implies that the environment has started to play a role in quenching star formation. At the same time, stronger clustering for galaxies with elevated SSFRs ("starbursts") might be attributed to an increased efficiency for galaxy interactions and mergers in dense environments.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Massive, red galaxies in a hierarchical universe-II Clustering of Extremely Red Objects

We present predictions for the clustering of Extremely Red Objects (EROs) in a Λcold dark matter Universe, using a semi-analytical galaxy formation model in combination with a cosmological N-body simulation. EROs are red, massive galaxies observed at 0.7< z < 3, and their numbers and properties have posed a challenge to hierarchical galaxy formation models. We analyse the halo occupation distribution and two-point correlation function of EROs, exploring how these quantities change with apparent magnitude, colour cut and redshift. Our model predicts a halo occupation distribution that is significantly different from that typically assumed. This is due to the inclusion of AGN feedback, which changes the slope and scatter of the luminosity-host halo mass relation above the mass where AGN feedback first becomes important. We predict that, on average, dark matter haloes with masses above 10^{13}h^{-1}Msun host at least one ERO at 1.5 < z < 2.5. Taking into account sample variance in observational estimates, the predicted angular clustering for EROs with either (R-K)>5 or (i-K)>4.5 is in reasonable agreement with current observations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the role of feedback in shaping the cosmic abundance and clustering of neutral atomic hydrogen in galaxies

We investigate the impact of feedback - from supernovae (SNe), active galactic nuclei (AGN) and a photo-ionizing background at high redshifts - on the neutral atomic hydrogen (HI) mass function, the $b_{\rm J}$ band luminosity function, and the spatial clustering of these galaxies at $z$=0. We use a version of the semi-analytical galaxy formation model GALFORM that calculates self-consistently the amount of HI in a galaxy as a function of cosmic time and links its star formation rate to its mass of molecular hydrogen (H$_2$). We find that a systematic increase or decrease in the strength of SNe feedback leads to a systematic decrease or increase in the amplitudes of the luminosity and HI mass functions, but has little influence on their overall shapes. Varying the strength of AGN feedback influences only the numbers of the brightest or most HI massive galaxies, while the impact of varying the strength of photo-ionization feedback is restricted to changing the numbers of the faintest or least HI massive galaxies.Our results suggest that the HI mass function is a more sensitive probe of the consequences of cosmological reionization for galaxy formation than the luminosity function. We find that increasing the strength of any of the modes of feedback acts to weaken the clustering strength of galaxies, regardless of their HI-richness. In contrast, weaker AGN feedback has little effect on the clustering strength whereas weaker SNe feedback increases the clustering strength of HI-poor galaxies more strongly than HI-rich galaxies. These results indicate that forthcoming HI surveys on next generation radio telescopes such as the Square Kilometre Array and its pathfinders will be exploited most fruitfully as part of multiwavelength survey campaigns.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The clustering of H-alpha emitters at z=2.23 from HiZELS

(abridged) We present a clustering analysis of 370 high-confidence H-alpha emitters (HAEs) at z=2.23. The HAEs are detected in the Hi-Z Emission Line Survey (HiZELS), a large-area blank field 2.121um narrowband survey using the United Kingdom Infrared Telescope (UKIRT) Wide Field Camera (WFCAM). Averaging the two-point correlation function of HAEs in two ~1 degree scale fields (United Kingdom Infrared Deep Sky Survey/Ultra Deep Survey [UDS] and Cosmological Evolution Survey [COSMOS] fields) we find a clustering amplitude equivalent to a correlation length of r_0=3.7+/-0.3 Mpc/h for galaxies with star formation rates of >7 M_sun/yr. The data are well-fitted by the expected correlation function of Cold Dark Matter, scaled by a bias factor: omega_HAE=b^2 omega_DM where b=2.4^{+0.1}_{-0.2}. The corresponding 'characteristic' mass for the halos hosting HAEs is log(M_h/[M_sun/h])=11.7+/-0.1. Comparing to the latest semi-analytic GALFORM predictions for the evolution of HAEs in a LCDM cosmology, we find broad agreement with the observations, with GALFORM predicting a HAE correlation length of ~4 Mpc/h. Motivated by this agreement, we exploit the simulations to construct a parametric model of the halo occupation distribution (HOD) of HAEs, and use this to fit the observed clustering. Our best-fitting HOD can adequately reproduce the observed angular clustering of HAEs, yielding an effective halo mass and bias in agreement with that derived from the scaled omega_DM fit, but with the relatively small sample size the current data provide a poor constraint on the HOD. Our results support the broad picture that 'typical' (~L*) star-forming galaxies have been hosted by dark matter haloes with M_h<10^12 M_sun/h since z~2, but with a broad occupation distribution and clustering that is likely to be a strong function of luminosity.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The contribution of star-forming galaxies to fluctuations in the cosmic background light

Star-forming galaxies which are too faint to be detected individually produce intensity fluctuations in the cosmic background light. This contribution needs to be taken into account as a foreground when using the primordial signal to constrain cosmological parameters. The extragalactic fluctuations are also interesting in their own right as they depend on the star formation history of the Universe and the way in which this connects with the formation of cosmic structure. We present a new framework which allows us to predict the occupation of dark matter haloes by star-forming galaxies and uses this information, in conjunction with an N-body simulation of structure formation, to predict the power spectrum of intensity fluctuations in the infrared background. We compute the emission from galaxies at far-infrared, millimetre and radio wavelengths. Our method gives accurate predictions for the clustering of galaxies both for the one halo and two halo terms. We illustrate our new framework using a previously published model which reproduces the number counts and redshift distribution of galaxies selected by their emission at $850\,μ$m. Without adjusting any of the model parameters, the predictions show encouraging agreement at high frequencies and on small angular scales with recent estimates of the extragalactic fluctuations in the background made from early data analysed by the Planck Collaboration. There are, however, substantial discrepancies between the model predictions and observations on large angular scales and at low frequencies, which illustrates the usefulness of the intensity fluctuations as a constraint on galaxy formation models.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The UV colours and dust attenuation of Lyman-break galaxies

Using GALFORM, a semi-analytical model of galaxy formation in the Lambda cold dark matter cosmology, we study the rest-frame ultraviolet (UV) colours of Lyman-break galaxies (LBGs) in the redshift range 2.5 < z < 10. As the impact of dust on UV luminosity can be dramatic, our model includes a self-consistent computation of dust attenuation based on a radiative transfer model. We find that intrinsically brighter galaxies suffer stronger dust attenuation than fainter ones, though the relation has a large scatter. The model predicts galaxies with UV colours consistent with the colour selection regions designed to select LBGs in observational surveys. We find that the drop-out technique that selects LBGs based on two rest-frame UV colours is robust and effective, selecting more than 70 per cent of UV bright galaxies at a given redshift. We investigate the impact on the predicted UV colours of varying selected model parameters. We find that the UV colours are most sensitive to the modelling of dust attenuation and in particular, to the extinction curve used in the radiative transfer calculation. If we assume a Milky Way dust extinction curve, the predicted UV continuum slopes are, in general, bluer than observed. However, we find that the opposite is true when using the Small Magellanic Cloud dust extinction curve. This demonstrates the strong dependence of UV colours on dust properties and highlights the inadequacy of using the UV continuum slope as a tracer of dust attenuation without any further knowledge of the galaxy inclination or dust characteristics in high redshift galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Are the 2dFGRS superstructures a problem for hierarchical models?

We introduce an objective method to assess the probability of finding extreme events in the distribution of cold dark matter such as voids, overdensities or very high mass haloes. Our approach uses an ensemble of N-body simulations of the hierar- chical clustering of dark matter to find extreme structures. The frequency of extreme events, in our case the cell or smoothing volume with the highest count of cluster-mass dark matter haloes, is well described by a Gumbel distribution. This distribution can then be used to forecast the probability of finding even more extreme events, which would otherwise require a much larger ensemble of simulations to quantify. We use our technique to assess the chance of finding concentrations of massive clusters or super- clusters, like the two found in the two-degree field galaxy redshift survey (2dFGRS), using a counts-in-cells analysis. The Gumbel distribution gives an excellent descrip- tion of the distribution of extreme cell counts across two large ensembles of simulations covering different cosmologies, and also when measuring the clustering in both real and redshift space. We find examples of structures like those found in the 2dFGRS in the simulations. The chance of finding such structures in a volume equal to that of the 2dFGRS is around 2%.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Testing dark energy using pairs of galaxies in redshift space

The distribution of angles subtended between pairs of galaxies and the line of sight,which is uniform in real space, is distorted by their peculiar motions, and has been proposed as a probe of cosmic expansion. We test this idea using N-body simulations of structure formation in a cold dark matter universe with a cosmological constant and in two variant cosmologies with different dark energy models. We find that the distortion of the distribution of angles is sensitive to the nature of dark energy. However, for the first time, our simulations also reveal dependences of the normalization of the distribution on both redshift and cosmology that have been neglected in previous work. This introduces systematics that severely limit the usefulness of the original method. Guided by our simulations, we devise a new, improved test of the nature of dark energy. We demonstrate that this test does not require prior knowledge of the background cosmology and that it can even distinguish between models that have the same baryonic acoustic oscillations and dark matter halo mass functions. Our technique could be applied to the completed BOSS galaxy redshift survey to constrain the expansion history of the Universe to better than 2%. The method will also produce different signals for dark energy and modified gravity cosmologies even when they have identical expansion histories, through the different peculiar velocities induced in these cases.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The nature and descendants of Lyman-break galaxies in the LambdaCDM cosmology

We predict the formation histories, properties and descendants of Lyman-break galaxies (LBGs) in the Lambda-CDM cosmology. In our model, which incorporates a top-heavy IMF in starbursts, we find that most LBGs are starbursts triggered by minor mergers of galaxies. The duration of the LBG phase is predicted to be quite short, ~20-60 Myr. We investigate the distributions of stellar and halo masses and morphologies for bright (L_UV > L*_UV) and faint (L_UV > 0.1 L*_UV) LBGs at z=3, z=6 and z=10 (where we classify LBGs according to their rest-frame UV luminosities relative the observed characteristic luminosity L*_UV at z \approx 3). Bright LBGs at z=3 are predicted to have median stellar masses ~ 1x10^9 Msun/h and host halo masses ~ 3x10^{11} Msun/h, and to be typically mildly disk-dominated in stellar mass. On the other hand, faint LBGs at z=10 are predicted to have median stellar masses of only ~ 1x10^7 Msun/h and host halo masses 2x10^{10} Msun/h, and to be generally bulge-dominated. Bright LBGs at z=3 evolve to present-day galaxies with median stellar mass ~ 5x^{10} Msun/h (comparable to the Milky Way), consisting of roughly equal numbers of disk- and bulge-dominated systems, and hosted by halos with median mass ~2x10^{13} Msun/h (corresponding to medium-size galaxy groups). The model predicts that 40% of Milky Way mass galaxies at the present-day have a bright LBG progenitor in the redshift range 3<z<4, while 95% have a faint LBG progenitor in the same redshift range, and 7% have a faint LBG progenitor at 10<z<11. With our multiwavelength model, we also investigate the overlap between the LBG population and that of submillimetre selected galaxies (SMGs); at z=3, only ~1% of bright LBGs are also predicted to also be bright SMGs (with an 850 mum flux in excess of 5 mJy).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Extremely Red Objects in a hierarchical universe

We analyse whether hierarchical formation models based on Lambda cold dark matter cosmology can produce enough massive red galaxies to match observations. For this purpose, we compare with observations the predictions from two published models for the abundance and redshift distribution of Extremely Red Objects (EROs), which are red, massive galaxies observed at z >= 1. One of the models invokes a "superwind" to regulate star formation in massive haloes and the other suppresses cooling through "radio-mode" AGN feedback. The first one underestimates the number counts of EROs by an order of magnitude, whereas the radio-mode AGN feedback model gives excellent agreement with the number counts of EROs and redshift distribution of K-selected galaxies. This study highlights the need to consider AGN feedback in order to understand the formation and evolution of massive galaxies at z >= 1.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Grand unification of AGN activity in the LambdaCDM cosmology

We track the co-evolution of supermassive black holes (SMBHs) and their host galaxies. The calculation is embedded in the GALFORM semi-analytical model which simulates the formation and evolution of galaxies in a cold dark matter (CDM) universe. During the evolution of the host galaxy, hot and cold gas are added to the SMBH by flows triggered by halo gas cooling, disc instabilities and galaxy mergers. This builds up the mass and spin of the BH, and the resulting accretion power regulates the gas cooling and subsequent star formation. The accretion flow is assumed to form a geometrically thin cool disc when the accretion rate exceeds 0.01\dot{M}_Edd, and a geometrically thick, radiatively inefficient hot flow when the accretion rate falls below this value. The resulting quasar optical luminosity function matches observations very well, and the mass of the SMBH correlates with the mass of the galaxy bulge as observed. The BH spin distribution depends strongly on whether the gas in any given accretion episode remains in the same plane (prolonged accretion) or whether, due to self-gravity, it fragments into multiple, randomly aligned accretion episodes (chaotic accretion). In the chaotic accretion model there is a clear correlation of spin with SMBH mass. Massive BHs (M>5\times10^8\Msun) are hosted by giant elliptical galaxies and are rapidly spinning, while lower mass BHs are hosted in spiral galaxies and have much lower spin. Using the Blandford-Znajek mechanism for jet production to calculate the jet power, our model is able to reproduce the radio loudness of radio galaxies, LINERS and Seyferts. This is the first confirmation that a CDM galaxy formation model can reproduce the observed phenomenology of AGN.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

How BAO measurements can fail to detect quintessence

We model the nonlinear growth of cosmic structure in different dark energy models, using large volume N-body simulations. We consider a range of quintessence models which feature both rapidly and slowly varying dark energy equations of state, and compare the growth of structure to that in a universe with a cosmological constant. The adoption of a quintessence model changes the expansion history of the universe, the form of the linear theory power spectrum and can alter key observables, such as the horizon scale and the distance to last scattering. The difference in structure formation can be explained to first order by the difference in growth factor at a given epoch; this scaling also accounts for the nonlinear growth at the 15% level. We find that quintessence models which feature late $(z<2)$, rapid transitions towards $w=-1$ in the equation of state, can have identical baryonic acoustic oscillation (BAO) peak positions to those in $Λ$CDM, despite being very different from $Λ$CDM both today and at high redshifts $(z \sim 1000)$. We find that a second class of models which feature non-negligible amounts of dark energy at early times cannot be distinguished from $Λ$CDM using measurements of the mass function or the BAO. These results highlight the need to accurately model quintessence dark energy in N-body simulations when testing cosmological probes of dynamical dark energy.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Massive, red galaxies in a hierarchical universe I. Counts of Extremely Red Objects and basic properties

We present predictions for the abundance and nature of Extremely Red Objects (EROs) in the Lambda cold dark matter model. EROs are red, massive galaxies observed at z>= 1 and their numbers and properties pose a challenge to hierarchical galaxy formation models. We compare the predictions from two published models, one of which invokes a "superwind" to regulate star formation in massive haloes and the other which suppresses gas cooling in haloes through "radio-mode" AGN feedback. The superwind model underestimates the number counts of EROs by an order of magnitude, whereas the radio-mode AGN feedback model gives excellent agreement with the number counts and redshift distribution of EROs. In the AGN feedback model the ERO population is dominated by old, passively evolving galaxies, whereas observations favour an equal split between old galaxies and dusty starbursts. Also, the model predicts a more extended redshift distribution of passive galaxies than is observed. These comparisons suggest that star formation may be quenched too efficiently in this model.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Modelling the dusty universe II: The clustering of submillimetre-selected galaxies

We combine the GALFORM semi-analytical model of galaxy formation, which predicts the star formation and merger histories of galaxies, the GRASIL spectro-photometric code, which calculates the spectral energy distributions (SEDs) of galaxies self-consistently including reprocessing of radiation by dust, and artificial neural networks (ANN), to investigate the clustering properties of galaxies selected by their emission at submillimetre wavelengths (SMGs). We use the Millennium Simulation to predict the spatial and angular distribution of SMGs. At redshift z = 2, we find that these galaxies are strongly clustered, with a comoving correlation length of r0 = 5.6 \pm 0.9 Mpc/h for galaxies with 850μm flux densities brighter than 5 mJy, in agreement with observations. We predict that at higher redshifts these galaxies trace denser and increasingly rarer regions of the universe. We present the predicted dependence of the clustering on luminosity, submillimetre colour, halo and total stellar masses. Interestingly, we predict tight relations between correlation length and halo and stellar masses, independent of sub-mm luminosity.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Modelling the spectral energy distribution of galaxies: introducing the artificial neural network

The spectral energy distribution of galaxies is a complex function of the star formation history and geometrical arrangement of stars and gas in galaxies. The computation of the radiative transfer of stellar radiation through the dust distribution is time-consuming. This aspect becomes unacceptable in particular when dealing with the predictions by semi-analytical galaxy formation models populating cosmological volumes, to be then compared with multi-wavelength surveys. Mainly for this aim, we have implemented an artificial neural network algorithm into the spectro-photometric and radiative transfer code GRASIL in order to compute the spectral energy distribution of galaxies in a short computing time. This allows to avoid the adoption of empirical templates that may have nothing to do with the mock galaxies output by models. The ANN has been implemented to compute the dust emission spectrum (the bottleneck of the computation), and separately for the star-forming molecular clouds and the diffuse dust (due to their different properties and dependencies). We have defined the input neurons effectively determining their emission, which means this implementation has a general applicability and is not linked to a particular galaxy formation model. We have trained the net for the disc and spherical geometries, and tested its performance to reproduce the SED of disc and starburst galaxies, as well as for a semi-analytical model for spheroidal galaxies. We have checked that for this model both the SEDs and the galaxy counts in the Herschel bands obtained with the ANN approximation are almost superimposed to the same quantities obtained with the full GRASIL. We conclude that this method appears robust and advantageous, and will present the application to a more complex SAM in another paper.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Predictions for Herschel from LambdaCDM: unveiling the cosmic star formation history

We use a model for the evolution of galaxies in the far-IR based on the LambdaCDM cosmology to make detailed predictions for upcoming cosmological surveys with the Herschel Space Observatory. We use the combined GALFORM semi-analytical galaxy formation model and GRASIL spectrophotometric code to compute galaxy SEDs including the reprocessing of radiation by dust. The model, which is the same as that in Baugh et al. (2005), assumes two different IMFs: a normal solar neighbourhood IMF for quiescent star formation in disks, and a very top-heavy IMF in starbursts triggered by galaxy mergers. We have shown previously that the top-heavy IMF appears necessary to explain the number counts and redshifts of faint sub-mm galaxies. In this paper, we present predictions for galaxy luminosity functions, number counts and redshift distributions in the Herschel imaging bands. We find that source confusion will be a serious problem in the deepest planned surveys. We also show predictions for physical properties such as star formation rates and stellar, gas and halo masses, together with fluxes at other wavelengths (from the far-UV to the radio) relevant for multi-wavelength follow-up observations. We investigate what fraction of the total IR emission from dust and of the high-mass star formation over the history of the Universe should be resolved by planned surveys with Herschel, and find a fraction ~30-50%, depending on confusion. Finally, we show that galaxies in Herschel surveys should be significantly clustered.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The evolution of AGN across cosmic time: what is downsizing?

We use a coupled model of the formation and evolution of galaxies and black holes (BH) to study the evolution of active galactic nuclei (AGN) in a cold dark matter universe. The model predicts the BH mass, spin and mass accretion history. BH mass grows via accretion triggered by discs becoming dynamically unstable or galaxy mergers (called the starburst mode) and accretion from quasi-hydrostatic hot gas haloes (called the hot-halo mode). By taking into account AGN obscuration, we obtain a very good fit to the observed luminosity functions (LF) of AGN (optical, soft and hard X-ray, and bolometric) for a wide range of redshifts (0<z<6). The model predicts a hierarchical build up of BH mass, with the typical mass of actively growing BHs increasing with decreasing redshift. Remarkably, despite this, we find downsizing in the AGN population, in terms of the differential growth with redshift of the space density of faint and bright AGN. This arises naturally from the interplay between the starburst and hot-halo accretion modes. The faint end of the LF is dominated by massive BHs experiencing quiescent accretion via a thick disc, primarily during the hot-halo mode. The bright end of the LF, on the other hand, is dominated by AGN which host BHs accreting close to or in excess of the Eddington limit during the starburst mode. The model predicts that the comoving space density of AGN peaks at z~3, similar to the star formation history. However, when taking into account obscuration, the space density of faint AGN peaks at lower redshift (z<2) than that of bright AGN (z~2-3). This implies that the cosmic evolution of AGN is shaped in part by obscuration.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The evolution of Lyman-break galaxies in CDM

We make a detailed investigation of the properties of Lyman-break galaxies (LBGs) in the LambdaCDM model. We present predictions for two published variants of the GALFORM semi-analytical model: the Baugh et al. (2005) model, which has star formation at high redshifts dominated by merger-driven starbursts with a top-heavy IMF, and the Bower et al. (2006) model, which has AGN feedback and a standard Solar neighbourhood IMF throughout. We show predictions for the evolution of the rest-frame far-UV luminosity function in the redshift range z=3-20, and compare with the observed luminosity functions of LBGs at z=3-10. We find that the Baugh et al. model is in excellent agreement with these observations, while the Bower et al. model predicts too many high-luminosity LBGs. Dust extinction, which is predicted self-consistently based on galaxy gas contents, metallicities and sizes, is found to have a large effect on LBG luminosities. We compare predictions for the size evolution of LBGs at different luminosities with observational data for 2<z<7, and find the Baugh et al. model to be in good agreement. We present predictions for stellar, halo and gas masses, star formation rates, circular velocities, bulge-to-disk ratios, gas and stellar metallicities and clustering bias, as functions of far-UV luminosity and redshift. We find broad consistency with current observational constraints. We then present predictions for the abundance and angular sizes of LBGs out to very high redshift (z<20), finding that planned deep surveys with JWST should detect objects out to z<15. The typical UV luminosities of galaxies are predicted to be very low at high redshifts, which has implications for detecting the galaxies responsible for reionizing the IGM; for example, at z=10, 50% of the ionizing photons are expected to be produced by galaxies fainter than M_AB(1500A)-5logh ~ -15.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Parameter Space of Galaxy Formation

Semi-analytic models are a powerful tool for studying the formation of galaxies. However, these models inevitably involve a significant number of poorly constrained parameters that must be adjusted to provide an acceptable match to the observed universe. In this paper, we set out to quantify the degree to which observational data-sets can constrain the model parameters. By revealing degeneracies in the parameter space we can hope to better understand the key physical processes probed by the data. We use novel mathematical techniques to explore the parameter space of the GALFORM semi-analytic model. We base our investigation on the Bower et al. 2006 version of GALFORM, adopting the same methodology of selecting model parameters based on an acceptable match to the local bJ and K luminosity functions. The model contains 16 parameters that are poorly constrained, and we investigate this parameter space using the Model Emulator technique, constructing a Bayesian approximation to the GALFORM model that can be rapidly evaluated at any point in parameter space. By combining successive waves of emulation, we show that only 0.26% of the initial volume is of interest for further exploration. However, within this region we show that the Bower et al. 2006 model is only one choice from an extended sub-space of model parameters that can provide equally acceptable fits. We explore the geometry of this region and begin to explore the physical connections between parameters that are exposed by this analysis. We also consider the impact of adding additional observational data to further constrain the parameter space.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Redshift Evolution of the Mass Function of Cold Gas in Hierarchical Galaxy Formation Models

(Abridged) Accurately predicting how the cosmic abundance of neutral hydrogen evolves with redshift is a challenging problem facing modellers of galaxy formation. We investigate the predictions of four currently favoured semi-analytical galaxy formation models applied to the Millennium simulation for the mass function of cold neutral gas (atomic and molecular) in galaxies as a function of redshift, and we use these predictions to construct number counts for the next generation of all-sky neutral atomic hydrogen (HI) surveys. Despite the different implementations of the physical ingredients of galaxy formation, we find that the model predictions are broadly consistent with one another; the key differences reflect how the models treat AGN feedback and how the timescale for star formation evolves with redshift. The models produce mass functions of cold gas in galaxies that are generally in good agreement with HI surveys at $z$=0. Interestingly we find that these mass functions do not evolve significantly with redshift. Adopting a simple conversion factor for cold gas mass to HI mass that we apply to all galaxies at all redshifts, we derive mass functions of HI in galaxies from the predicted mass functions of cold gas, which we use to predict the number counts of sources likely to be detected by HI surveys on next generation radio telescopes such as the Square Kilometre Array and its pathfinders. We show also how adopting a conversion factor for cold gas mass to HI mass that varies from galaxy to galaxy impacts on number counts. In addition, we examine how the typical angular sizes of galaxies vary with redshift. Taken together, these results make clear that forthcoming HI surveys will provide important and powerful tests of theoretical galaxy formation models.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The role of submillimetre galaxies in hierarchical galaxy formation

We study the role of submillimetre galaxies (SMGs) in the galaxy formation process in the Lambda Cold Dark Matter cosmology. We use the Baugh et al. (2005) semi-analytical model, which matches the observed SMG number counts and redshift distribution by assuming a top-heavy initial mass function (IMF) in bursts triggered by galaxy mergers. We build galaxy merger trees and follow the evolution and properties of SMGs and their descendants. Our primary sample of model SMGs consists of galaxies which had 850 mu fluxes brighter than 5 mJy at some redshift z>1. Our model predicts that the present-day descendants of such SMGs cover a wide range of stellar masses ~ 10^{10} - 10^{12} Msun/h, with a median ~ 10^{11} Msun/h, and that more than 70% of these descendants are bulge-dominated. More than 50% of present day galaxies with stellar masses larger than 7 x 10^{11} Msun/h are predicted to be descendants of such SMGs. We find that although SMGs make an important contribution to the total star formation rate at z~2, the final stellar mass produced in the submillimetre phase contributes only 0.2% of the total present-day stellar mass, and 2% of the stellar mass of SMG descendants, in stark contrast to the popular picture in which the SMG phase marks the production of the bulk of the mass of present day massive ellipticals.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The spatial distribution of cold gas in hierarchical galaxy formation models

The distribution of cold gas in dark matter haloes is driven by key processes in galaxy formation: gas cooling, galaxy mergers, star formation and reheating of gas by supernovae. We compare the predictions of four different galaxy formation models for the spatial distribution of cold gas. We find that satellite galaxies make little contribution to the abundance or clustering strength of cold gas selected samples, and are far less important than they are in optically selected samples. The halo occupation distribution function of present-day central galaxies with cold gas mass > 10^9 h^-1 Msun is peaked around a halo mass of ~ 10^11 h^-1 Msun, a scale that is set by the AGN suppression of gas cooling. The model predictions for the projected correlation function are in good agreement with measurements from the HI Parkes All-Sky Survey. We compare the effective volume of possible surveys with the Square Kilometre Array with those expected for a redshift survey in the near-infrared. Future redshift surveys using neutral hydrogen emission will be competitive with the most ambitious spectroscopic surveys planned in the near-infrared.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Cosmological parameter constraints from SDSS luminous red galaxies: a new treatment of large-scale clustering

We apply a new model for the spherically averaged correlation function at large pair separations to the measurement of the clustering of luminous red galaxies (LRGs) made from the SDSS by Cabre and Gaztanaga(2009). Our model takes into account the form of the BAO peak and the large scale shape of the correlation function. We perform a Monte Carlo Markov chain analysis for different combinations of datasets and for different parameter sets. When used in combination with a compilation of the latest CMB measurements, the LRG clustering and the latest supernovae results give constraints on cosmological parameters which are comparable and in remarkably good agreement, resolving the tension reported in some studies. The best fitting model in the context of a flat, Lambda-CDM cosmology is specified by Omega_m=0.261+-0.013, Omega_b=0.044+-0.001, n_s=0.96+-0.01, H_0=71.6+-1.2 km/s/Mpc and sigma_8=0.80+-0.02. If we allow the time-independent dark energy equation of state parameter to vary, we find results consistent with a cosmological constant at the 5% level using all data sets: w_DE=-0.97+-0.05. The large scale structure measurements by themselves can constrain the dark energy equation of state parameter to w_DE=-1.05+-0.15, independently of CMB or supernovae data. We do not find convincing evidence for an evolving equation of state. We provide a set of "extended distance priors" that contain the most relevant information from the CMB power spectrum and the shape of the LRG correlation function which can be used to constrain dark energy models and spatial curvature. Our model should provide an accurate description of the clustering even in much larger, forthcoming surveys, such as those planned with NASA's JDEM or ESA's Euclid mission.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Empirical H-alpha emitter count predictions for dark energy surveys

Future galaxy redshift surveys aim to measure cosmological quantities from the galaxy power spectrum. A prime example is the detection of baryonic acoustic oscillations (BAOs), providing a standard ruler to measure the dark energy equation of state, w(z), to high precision. The strongest practical limitation for these experiments is how quickly accurate redshifts can be measured for sufficient galaxies to map the large-scale structure. A promising strategy is to target emission-line (i.e. star-forming) galaxies at high-redshift (z~0.5-2); not only is the space density of this population increasing out to z~2, but also emission-lines provide an efficient method of redshift determination. Motivated by the prospect of future dark energy surveys targeting H-alpha emitters at near-infrared wavelengths (i.e. z>0.5), we use the latest empirical data to model the evolution of the H-alpha luminosity function out to z~2, and thus provide predictions for the abundance of H-alpha emitters for practical limiting fluxes. We caution that the estimates presented in this work must be tempered by an efficiency factor, epsilon, giving the redshift success rate from these potential targets. For a range of practical efficiencies and limiting fluxes, we provide an estimate of nP_{0.2}, where n is the 3D galaxy number density and P_{0.2} is the galaxy power spectrum evaluated at k=0.2h/Mpc. Ideal surveys must provide nP_{0.2}>1 in order to balance shot-noise and cosmic variance errors. We show that a realistic emission-line survey (epsilon=0.5) could achieve nP_{0.2}=1 out to z~1.5 with a limiting flux of 10^{-16} erg/s/cm^{-2}. If the limiting flux is a factor 5 brighter, then this goal can only be achieved out to z~0.5, highlighting the importance of survey depth and efficiency in cosmological redshift surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Modelling galaxy clustering: Is new physics needed in galaxy formation models?

The clustering amplitude of galaxies depends on their intrinsic luminosity. We compare the properties of publicly available galaxy formation models with clustering measurements from the two-degree field galaxy redshift survey. The model predictions show the same qualitative behaviour as the data but fail to match the observations at the level of accuracy at which current measurements can be made. We demonstrate that this is due to the model producing too many satellite galaxies in massive haloes. We implement simple models to describe two new processes, satellite-satellite mergers and the tidal dissolution of satellites to investigate their impact on the predicted clustering. We find that both processes need to be included in order to produce a model which matches the observations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Modelling the dusty universe I: Introducing the artificial neural network and first applications to luminosity and colour distributions

We introduce a new technique based on artificial neural networks which allows us to make accurate predictions for the spectral energy distributions (SEDs) of large samples of galaxies, at wavelengths ranging from the far-ultra-violet to the sub-millimetre and radio. The neural net is trained to reproduce the SEDs predicted by a hybrid code comprised of the GALFORM semi-analytical model of galaxy formation, which predicts the full star formation and galaxy merger histories, and the GRASIL spectro-photometric code, which carries out a self-consistent calculation of the SED, including absorption and emission of radiation by dust. Using a small number of galaxy properties predicted by GALFORM, the method reproduces the luminosities of galaxies in the majority of cases to within 10% of those computed directly using GRASIL. The method performs best in the sub-mm and reasonably well in the mid-infrared and the far-ultra-violet. The luminosity error introduced by the method has negligible impact on predicted statistical distributions, such as luminosity functions or colour distributions of galaxies. We use the neural net to predict the overlap between galaxies selected in the rest-frame UV and in the observer-frame sub-mm at z=2. We find that around half of the galaxies with a 850um flux above 5 mJy should have optical magnitudes brighter than R_AB < 25 mag. However, only 1% of the galaxies selected in the rest-frame UV down to R_AB < 25 mag should have 850um fluxes brighter than 5 mJy. Our technique will allow the generation of wide-angle mock catalogues of galaxies selected at rest-frame UV or mid- and far-infrared wavelengths.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Probing dark energy with future redshift surveys: A comparison of emission line and broad band selection in the near infrared

Future galaxy surveys will map the galaxy distribution in the redshift interval $0.5<z<2$ using near-infrared cameras and spectrographs. The primary science goal of such surveys is to constrain the nature of the dark energy by measuring the large-scale structure of the Universe. This requires a tracer of the underlying dark matter which maximizes the useful volume of the survey. We investigate two potential survey selection methods: an emission line sample based on the \ha line and a sample selected in the H-band. We present predictions for the abundance and clustering of such galaxies, using two published versions of the \galform galaxy formation model. Our models predict that \ha selected galaxies tend to avoid massive dark matter haloes and instead trace the surrounding filamentary structure; H-band selected galaxies, on the other hand, are found in the highest mass haloes. This has implications for the measurement of the rate at which fluctuations grow due to gravitational instability. We use mock catalogues to compare the effective volumes sampled by a range of survey configurations. To give just two examples: a redshift survey down to $H_{\rm AB}=22$ samples an effective volume that is $\sim 5-10$ times larger than that probed by an \ha survey with $\logfha > -15.4$; a flux limit of at least $\logfha = -16$ is required for an \ha sample to become competitive in effective volume.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

The Fate of Substructures in Cold Dark Matter Haloes

We use the Millennium Simulation, a large, high resolution N-body simulation of the evolution of structure in a LambdaCDM cosmology, to study the properties and fate of substructures within a large sample of dark matter haloes. We find that the subhalo mass function departs significantly from a power law at the high mass end. We also find that the radial and angular distributions of substructures depend on subhalo mass. In particular, high mass subhaloes tend to be less radially concentrated and to have angular distributions closer to the direction perpendicular to the spin of the host halo than their less massive counterparts. We find that mergers between subhaloes occur. These tend to be between substructures that were already dynamically associated before accretion into the main halo. For subhaloes larger than 0.001 times the mass of the host halo, it is more likely that the subhalo will merge with the central or main subhalo than with another subhalo larger than itself. For lower masses, subhalo-subhalo mergers become equally likely to mergers with the main subhalo. Our results have implications for the variation of galaxy properties with environment and for the treatment of mergers in galaxy formation models.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Can the faint sub-mm galaxies be explained in the Lambda-CDM model?

We present predictions for the abundance of sub-mm galaxies (SMGs) and Lyman-break galaxies (LBGs) in the $Λ$CDM cosmology. A key feature of our model is the self-consistent calculation of the absorption and emission of radiation by dust. The new model successfully matches the LBG luminosity function, as well reproducing the properties of the local galaxy population in the optical and IR. The model can also explain the observed galaxy number counts at $850\mum$, but only if we assume a top-heavy IMF for the stars formed in bursts. The predicted redshift distribution of SMGs depends relatively little on their flux over the range 1-$10\mjy$, with a median value of $z\approx 2.0$ at a flux of $5\mjy$, in very good agreement with the recent measurement by Chapman et al The counts of SMGs are predicted to be dominated by ongoing starbursts. However, in the model these bursts are responsible for making only a few per cent of the stellar mass locked up in massive ellipticals at the present day.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

The 2dF Galaxy Redshift Survey: luminosity functions by density environment and galaxy type

We use the 2dF Galaxy Redshift Survey to measure the dependence of the bJ-band galaxy luminosity function on large-scale environment, defined by density contrast in spheres of radius 8h-1Mpc, and on spectral type, determined from principal component analysis. We find that the galaxy populations at both extremes of density differ significantly from that at the mean density. The population in voids is dominated by late types and shows, relative to the mean, a deficit of galaxies that becomes increasingly pronounced at magnitudes brighter than M_bJ-5log10h <-18.5. In contrast, cluster regions have a relative excess of very bright early-type galaxies with M_bJ-5log10h < -21. Differences in the mid to faint-end population between environments are significant: at M_bJ-5log10h=-18 early and late-type cluster galaxies show comparable abundances, whereas in voids the late types dominate by almost an order of magnitude. We find that the luminosity functions measured in all density environments, from voids to clusters, can be approximated by Schechter functions with parameters that vary smoothly with local density, but in a fashion which differs strikingly for early and late-type galaxies. These observed variations, combined with our finding that the faint-end slope of the overall luminosity function depends at most weakly on density environment, may prove to be a significant challenge for models of galaxy formation.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Galaxy Evolution and Cosmology with the Square Kilometre Array

The present-day Universe is seemingly dominated by dark energy and dark matter, but mapping the normal (baryonic) content remains vital for both astrophysics - understanding how galaxies form - and astro-particle physics - inferring properties of the dark components. The Square Kilometre Array (SKA) will provide the only means of studying the cosmic evolution of neutral Hydrogen (HI) which, alongside information on star formation from the radio continuum, is needed to understand how stars formed from gas within dark-matter over-densities and the roles of gas accretion and galaxy merging. `All hemisphere' HI redshift surveys to redshift 1.5 are feasible with wide-field-of-view realizations of the SKA and, by measuring the galaxy power spectrum in exquisite detail, will allow the first precise studies of the equation-of-state of dark energy. The SKA will be capable of other uniquely powerful cosmological studies including the measurement of the dark-matter power spectrum using weak gravitational lensing, and the precise measurement of H0 using extragalactic water masers. The SKA is likely to become the premier dark-energy-measuring machine, bringing breakthroughs in cosmology beyond those likely to be made possible by combining CMB (e.g. Planck), optical (e.g. LSST, SNAP) and other early-21st-century datasets.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Galaxy groups in the 2dFGRS: the group-finding algorithm and the 2PIGG catalogue

The construction of a catalogue of galaxy groups from the 2-degree Field Galaxy Redshift Survey (2dFGRS) is described. Groups are identified by means of a friends-of-friends percolation algorithm which has been thoroughly tested on mock versions of the 2dFGRS generated from cosmological N-body simulations. The tests suggest that the algorithm groups all galaxies that it should be grouping, with an additional 40% of interlopers. About 55% of the ~190000 galaxies considered are placed into groups containing at least two members of which ~29000 are found. Of these, ~7000 contain at least four galaxies, and these groups have a median redshift of 0.11 and a median velocity dispersion of 260km/s. This 2dFGRS Percolation-Inferred Galaxy Group (2PIGG) catalogue represents the largest available homogeneous sample of galaxy groups. It is publicly available on the WWW.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Galaxy groups in the 2dFGRS: the luminous content of the groups

The 2dFGRS Percolation-Inferred Galaxy Group (2PIGG) catalogue of ~29000 objects is used to study the luminous content of galaxy systems of various sizes. Mock galaxy catalogues constructed from cosmological simulations are used to gauge the accuracy with which intrinsic group properties can be recovered. A Schechter function is found to be a reasonable fit to the galaxy luminosity functions in groups of different mass in the real data. The characteristic luminosity L* is larger for more massive groups. However, the mock data show that the shape of the recovered luminosity function is expected to differ from the true shape, and this must be allowed for when interpreting the data. The variation of halo mass-to-light ratio with group size is studied in both these wavebands. A robust trend of increasing M/L with increasing group luminosity is found in the 2PIGG data. From groups with L_bj=10^{10}Lsol to those 100 times more luminous, the typical bj-band M/L increases by a factor of 5, whereas the rf-band M/L grows by a factor of 3.5. These trends agree well with the simulations, which also predict a minimum M/L on a scale corresponding to the Local Group. Our data indicate that if such a minimum exists, then it must occur at L<~10^{10}Lsol, below the range accurately probed by the 2PIGG catalogue. According to the mock data, the bj M/Ls of the largest groups are expected to be approximately 1.1 times the global value. Assuming that this correction applies to the real data yields an estimate of Omega_m=0.26+/-0.03 (statistical).

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

The 2dF Galaxy Redshift Survey: Hierarchical galaxy clustering

We use the two-degree field Galaxy Redshift Survey (2dFGRS) to test the hierarchical scaling hypothesis: namely, that the $p$-point galaxy correlation functions can be written in terms of the two point correlation function or variance. This scaling is expected if an initially Gaussian distribution of density fluctuations evolves under the action of gravitational instability. We measure the volume averaged $p$-point correlation functions using a counts in cells technique applied to a volume limited sample of 44,931 $L_*$ galaxies. We demonstrate that $L_{*}$ galaxies display hierarchical clustering up to order $p=6$ in redshift space. The variance measured for $L_{*}$ galaxies is in excellent agreement with the predictions from a $Λ$-cold dark matter N-body simulation. This applies to all cell radii considered, $0.3<(R/h^{-1}{\rm Mpc})<30$. However, the higher order correlation functions of $L_*$ galaxies have a significantly smaller amplitude than is predicted for the dark matter for $R<10h^{-1}$Mpc. This disagreement implies that a non-linear bias exists between the dark matter and $L_*$ galaxies on these scales. We also show that the presence of two rare, massive superclusters in the 2dFGRS has an impact on the higher-order clustering moments measured on large scales.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

The 2dF Galaxy Redshift Survey: higher order galaxy correlation functions

We measure moments of the galaxy count probability distribution function in the two-degree field galaxy redshift survey (2dFGRS). The survey is divided into volume limited subsamples in order to examine the dependence of the higher order clustering on galaxy luminosity. We demonstrate the hierarchical scaling of the averaged p-point galaxy correlation functions, xibar_p, up to p=6. The hierarchical amplitudes, S_p = xibar_p/xibar_2^{p-1}, are approximately independent of the cell radius used to smooth the galaxy distribution on small to medium scales. On larger scales we find the higher order moments can be strongly affected by the presence of rare, massive superstructures in the galaxy distribution. The skewness S_3 has a weak dependence on luminosity, approximated by a linear dependence on log luminosity. We discuss the implications of our results for simple models of linear and non-linear bias that relate the galaxy distribution to the underlying mass.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

The 2dF Galaxy Redshift Survey: stochastic relative biasing between galaxy populations

It is well known that the clustering of galaxies depends on galaxy type.Such relative bias complicates the inference of cosmological parameters from galaxy redshift surveys, and is a challenge to theories of galaxy formation and evolution. In this paper we perform a joint counts-in-cells analysis on galaxies in the 2dF Galaxy Redshift Survey, classified by both colour and spectral type, eta, as early or late type galaxies. We fit three different models of relative bias to the joint probability distribution of the cell counts, assuming Poisson sampling of the galaxy density field. We investigate the nonlinearity and stochasticity of the relative bias, with cubical cells of side 10Mpc \leq L \leq 45Mpc (h=0.7). Exact linear bias is ruled out with high significance on all scales. Power law bias gives a better fit, but likelihood ratios prefer a bivariate lognormal distribution, with a non-zero `stochasticity' - i.e. scatter that may result from physical effects on galaxy formation other than those from the local density field. Using this model, we measure a correlation coefficient in log-density space (r_LN) of 0.958 for cells of length L=10Mpc, increasing to 0.970 by L=45Mpc. This corresponds to a stochasticity sigma_b/bhat of 0.44\pm0.02 and 0.27\pm0.05 respectively. For smaller cells, the Poisson sampled lognormal distribution presents an increasingly poor fit to the data, especially with regard to the fraction of completely empty cells. We compare these trends with the predictions of semianalytic galaxy formation models: these match the data well in terms of overall level of stochasticity, variation with scale, and fraction of empty cells.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

The 2dF Galaxy Redshift Survey: The clustering of galaxy groups

We measure the clustering of galaxy groups in the 2dFGRS Percolation-Inferred Galaxy Group (2PIGG) catalogue. The 2PIGG sample has 29,000 groups with at least two members. The clustering amplitude of the full 2PIGG catalogue is weaker than that of 2dFGRS galaxies, in agreement with theoretical predictions. We have subdivided the 2PIGG catalogue into samples that span a factor of 25 in median total luminosity. Our correlation function measurements span an unprecedented range of clustering strengths, connecting the regimes probed by groups fainter than L* galaxies and rich clusters. There is a steady increase in clustering strength with group luminosity; the most luminous groups are ten times more strongly clustered than the full 2PIGG catalogue. We demonstrate that the 2PIGG results are in very good agreement with the clustering of groups expected in the LCDM model.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

The 2dF Galaxy Redshift Survey: voids and hierarchical scaling models

We measure the redshift space reduced void probability function (VPF) for 2dFGRS volume limited galaxy samples covering the absolute magnitude range M_bJ-5logh=-18 to -22. Theoretically, the VPF connects the distribution of voids to the moments of galaxy clustering of all orders, and can be used to discriminate clustering models in the weakly non-linear regime. The reduced VPF measured from the 2dFGRS is in excellent agreement with the paradigm of hierarchical scaling of the galaxy clustering moments. The accuracy of our measurement is such that we can rule out, at a very high significance, popular models for galaxy clustering, including the lognormal distribution. We demonstrate that the negative binomial model gives a very good approximation to the 2dFGRS data over a wide range of scales, out to at least 20h-1Mpc. Conversely, the reduced VPF for dark matter in a LambdaCDM universe does appear to be lognormal on small scales but deviates significantly beyond \approx 4h-1Mpc. We find little dependence of the 2dFGRS reduced VPF on galaxy luminosity. Our results hold independently in both the north and south Galactic pole survey regions.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Where are the stars?

The 2dFGRS is used in conjunction with the 2MASSXSC to study the near-IR light and stellar mass content of the local Universe. Mock galaxy catalogues, constructed from cosmological N-body simulations and semi-analytical galaxy formation models, are used to gauge the accuracy with which quantities can be recovered. The mean luminosity densities of the Universe are found to be rho_J=(3.57+/-0.11)*10^8 h Lsol/Mpc^3 and rho_KS=(7.04+/-0.23)*10^8 h Lsol/Mpc^3 (statistical uncertainty). Using the 2PIGG catalogue, the group dynamical mass-to-light ratio in the K_S band is found to increase by a factor of ~3 when going from groups with total bJ-band luminosities of 3*10^10 h^-2 Lsol to rich clusters, which have typical values of Upsilon_K~80 h Upsilon_sol. Taking into account the bias introduced by uncertainties in estimating galaxy stellar masses from luminosities, a value of Omega_* h=(0.99+/-0.03)*10^-3 is measured, assuming that a Kennicutt stellar IMF is applicable to all galaxies. The 2PIGGs are then used to study the distribution of the stellar content of the local Universe. The three main conclusions are: (1) a slowly rising stellar M/L_KS is found with the clusters having the largest value of ~0.6 Upsilon_sol, (2) in contrast, the fraction of mass in stars decreases with increasing group size, reaching ~5*10^-3 h for the rich clusters, and (3) in answer to the question posed in the title, most stellar mass is contained in Local Group-sized objects (M~2*10^12 h^-1 Msol) with only ~2% in clusters with M>5*10^14 h^-1 Msol.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

The 2dF Galaxy Redshift Survey: galaxy clustering per spectral type

We have calculated the two-point correlation functions in redshift space, xi(sigma,pi), for galaxies of different spectral types in the 2dF Galaxy Redshift Survey. Using these correlation functions we are able to estimate values of the linear redshift-space distortion parameter, beta = Omega_m^0.6/b, the pairwise velocity dispersion, a, and the real-space correlation function, xi(r), for galaxies with both relatively low star-formation rates (for which the present rate of star formation is less than 10% of its past averaged value) and galaxies with higher current star-formation activity. At small separations, the real-space clustering of passive galaxies is very much stronger than that of the more actively star-forming galaxies; the correlation-function slopes are respectively 1.93 and 1.50, and the relative bias between the two classes is a declining function of radius. On scales larger than 10 h^-1 Mpc there is evidence that the relative bias tends to a constant, b(passive)/b(active) ~ 1. This result is consistent with the similar degrees of redshift-space distortions seen in the correlation functions of the two classes -- the contours of xi(sigma,pi) require beta(active)=0.49+/-0.13, and beta(passive)=0.48+/-0.14. The pairwise velocity dispersion is highly correlated with beta. However, despite this a significant difference is seen between the two classes. Over the range 8-20 h^-1 Mpc, the pairwise velocity dispersion has mean values 416+/-76 km/s and 612+/-92 km/s for the active and passive galaxy samples respectively. This is consistent with the expectation from morphological segregation, in which passively evolving galaxies preferentially inhabit the cores of high-mass virialised regions.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

The 2dF Galaxy Redshift Survey: constraints on cosmic star-formation history from the cosmic spectrum

We present the first results on the history of star formation in the Universe based on the `cosmic spectrum', in particular, the volume-averaged, luminosity-weighted, stellar absorption line spectrum of present day galaxies from the 2dFGRS. This method is novel in that unlike previous studies it is not an estimator based on total luminosity density. The cosmic spectrum is fitted with models of population synthesis, tracing the history of star formation prior to the epoch of the observed galaxies, using a method we have developed which decouples continuum and spectral-line variations and is robust against spectrophotometric uncertainties. The cosmic spectrum can only be fitted with models incorporating chemical evolution and indicates there was a peak of star-formation rate in the past of at least three times the current value and that the increase back to z=1, assuming it scales as (1+z)^beta, has a strong upper limit of beta<5. We find in the general case there is some model degeneracy between star formation at low and high redshift. However, if we incorporate previous work on star formation at z<1 we can put strong upper limits on the SFR at z>1: e.g., if beta>2 then the SFR for 1<z<5 scales as (1+z)^alpha with alpha<2. This is equivalent to stating that no more than 80% of stars in the Universe formed at z>1. Our results are consistent with the best-fit results from compilations of cosmic SFR estimates based on UV luminosity density, which give 1.8<beta<2.9 and -1.0<alpha<0.7, and are also consistent with estimates of Omega_stars based on the K-band luminosity density.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

The 2dF Galaxy Redshift Survey: Galaxy luminosity functions per spectral type

We calculate the optical bj luminosity function of the 2dF Galaxy Redshift Survey (2dFGRS) for different subsets defined by their spectral properties. These spectrally selected subsets are defined using a new parameter, eta, which is a linear combination of the first two projections derived from a Principal Component Analysis. This parameter eta identifies the average emission and absorption line strength in the galaxy rest-frame spectrum and hence is a useful indicator of the present star formation. We use a total of 75,000 galaxies in our calculations, chosen from a sample of high signal-to-noise ratio, low redshift galaxies observed before January 2001. We find that there is a systematic steepening of the faint end slope (alpha) as one moves from passive (alpha = -0.54) to active (alpha = -1.50) star-forming galaxies, and that there is also a corresponding faintening of the rest-frame characteristic magnitude M* - 5\log_10(h) (from -19.6 to -19.2). We also show that the Schechter function provides a poor fit to the quiescent (Type 1) LF for very faint galaxies (M - 5log_10(h) fainter than -16.0), perhaps suggesting the presence of a significant dwarf population. The luminosity functions presented here give a precise confirmation of the trends seen previously in a much smaller preliminary 2dFGRS sample, and in other surveys. We also present a new procedure for determining self-consistent K-corrections and investigate possible fibre-aperture biases.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

The 2dF Galaxy Redshift Survey: The dependence of galaxy clustering on luminosity and spectral type

We investigate the dependence of galaxy clustering on luminosity and spectral type using the 2dF Galaxy Redshift Survey (2dFGRS). Spectral types are assigned using the principal component analysis of Madgwick et al. We divide the sample into two broad spectral classes: galaxies with strong emission lines (`late-types'), and more quiescent galaxies (`early-types'). We measure the clustering in real space, free from any distortion of the clustering pattern due to peculiar velocities, for a series of volume-limited samples. The projected correlation functions of both spectral types are well described by a power law for transverse separations in the range 2 < (sigma/Mpc/h) < 15, with a marginally steeper slope for early-types than late-types. Both early and late types have approximately the same dependence of clustering strength on luminosity, with the clustering amplitude increasing by a factor of ~2.5 between L* and 4 L*. At all luminosities, however, the correlation function amplitude for the early-types is ~50% higher than that of the late-types. These results support the view that luminosity, and not type, is the dominant factor in determining how the clustering strength of the whole galaxy population varies with luminosity.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

The Effects of Photoionization on Galaxy Formation -- III: Environmental Dependence in the Luminosity Function

Using semi-analytic modeling techniques, we calculate the luminosity function of galaxy populations residing in cold dark matter halos of different mass. We pay particular attention to the influence of the reionization of the Universe on the number of faint galaxies and to the effects of dynamical friction and tidal limitation of satellites on the number of bright galaxies. We find substantial differences in the shapes of the galaxy luminosity functions in halos of different mass which reflect generic features of the cold dark matter model of galaxy formation and thus offer the opportunity to test it. We then consider how the individual halo luminosity functions combine together to produce the global luminosity function. Surprisingly, the global function ends up having a shallower faint end slope than those of the constituent halo luminosity functions. We compare our model predictions with the limited datasets compiled by Trentham & Hodgkin. We find good agreement with the luminosity functions measured in the Virgo and Coma clusters but significant disagreement with the luminosity functions measured in the Local Group and Ursa Minor cluster. We speculate on possible inadequacies in our modeling and in the existing observational samples. The luminosity functions of galaxies in groups and clusters identified in the 2dF and SDSS galaxy redshift surveys offer the prospect of testing galaxy formation models in detail.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

A measurement of the cosmological mass density from clustering in the 2dF Galaxy Redshift Survey

The large-scale structure in the distribution of galaxies is thought to arise from the gravitational instability of small fluctuations in the initial density field of the universe. A key test of this hypothesis is that superclusters of galaxies in the process of formation should generate systematic infall of other galaxies. This would be evident in the pattern of recessional velocities, causing an anisotropy in the inferred spatial clustering of galaxies. Here we report a precise measurement of this clustering, using the redshifts of more than 141,000 galaxies from the two-degree-field galaxy redshift survey. We determine the parameter beta = Omega^{0.6}/b = 0.43 +- 0.07, where Omega is the total mass-density parameter and b is a measure of the `bias' of the luminous galaxies in the survey. Combined with the anisotropy of the cosmic microwave background, our results favour a low-density universe with Omega approximately 0.3.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

The 2dF Galaxy Redshift Survey: Luminosity dependence of galaxy clustering

We investigate the dependence of the strength of galaxy clustering on intrinsic luminosity using the Anglo-Australian two degree field galaxy redshift survey (2dFGRS). The 2dFGRS is over an order of magnitude larger than previous redshift surveys used to address this issue. We measure the projected two-point correlation function of galaxies in a series of volume-limited samples. The projected correlation function is free from any distortion of the clustering pattern induced by peculiar motions and is well described by a power-law in pair separation over the range 0.1 < r /h Mpc < 10. The clustering of L* galaxies in real space is well fit by a correlation length r0 = 4.9 +/- 0.3 /h Mpc and power-law slope gamma = 1.71 +/- 0.06. The clustering amplitude increases slowly with absolute magnitude for galaxies fainter than M*, but rises more strongly at higher luminosities. At low luminosities, our results agree with measurements from the SSRS2 by Benoist et al. However, we find a weaker dependence of clustering strength on luminosity at the highest luminosities. The correlation function amplitude increases by a factor of 4.0 between $M_{b_{J}} -5\log_{10}h = -18$ and -22.5, and the most luminous galaxies are 3.0 times more strongly clustered than L* galaxies. The power-law slope of the correlation function shows remarkably little variation for samples spanning a factor of 20 in luminosity. Our measurements are in very good agreement with the predictions of the hierarchical galaxy formation models of Benson et al.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

The seeds of rich galaxy clusters in the Universe

The discovery of a population of young galaxies at an epoch when the universe was about one tenth of its current age has shed new light on the question of when and how galaxies formed. Within the context of popular models this is the population of primeval galaxies that built themselves up to the size of present--day galaxies through the process of repeated mergers. But the recent detection of a large concentration of these primeval galaxies (Steidel et al. 1997) appears to be incompatible with hierarchical clustering models, which generally predict that clusters of this size are fully formed later in time. Here we use a combination of two powerful theoretical techniques --semi-analytic modelling and N-body simulations-- to show that such large concentrations should be quite common in a universe dominated by cold dark matter, and that they are the progenitors of the rich galaxy clusters seen today. We predict the clustering properties of primeval galaxies which should, when compared with data that will be collected in the near future, test our understanding of galaxy formation within the framework of a universe dominated by cold dark matter.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

A wide-field K-band survey - II. Galaxy clustering

We present the first measurement of the angular correlation function in a $K$-selected galaxy survey, from two fields covering almost 10 square degrees. The angular correlation function at small angles is well described by a $θ^{-0.8}$ power law, as for optically selected samples. The clustering amplitude is reduced by a factor of $\sim 4$ between galaxies with $K<15$ and a fainter magnitude slice with $15<K<16$. This allows us to place constraints upon the redshift distribution of the galaxies and their spatial correlation function. We find no clear evidence for a change in clustering amplitude when galaxies are selected by their observed $B-K$ colours.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Evolution of the Hubble Sequence in Hierarchical Models for Galaxy Formation

We present a model for the broad morphological distinction between the disk and spheroidal components of galaxies. Elaborating on the hierarchical clustering scheme of galaxy formation proposed by Cole et al., we assume that galaxies form stars quiescently in a disk until they are disrupted into a spheroidal configuration by mergers. We calculate formation and merging histories, and the evolution in colour, luminosity and morphology of the galaxy populations in different environments. Roughly $ 50\ %$ of ellipticals in our model have undergone a major merger since $z = 0.5$, yet in spite of this we find that cluster ellipticals have colour-magnitude diagrams with remarkably small scatter. The morphological mix of galaxies that become rich cluster members at high redshift is dominated by spiral galaxies, due to the long timescale for galaxy mergers compared with the timescale for cluster assembly at high redshift. The assembly of low redshift clusters is slower, allowing more galaxy mergers to occur in the progenitor halos. As a result $z=0$ rich clusters become E/S0 dominated and we find a ``Butcher-Oemler'' effect that becomes weaker for poorer groups at high redshift. The field luminosity function of red galaxies shows little evolution out to $z\simeq 1$ and the reddest galaxies at these redshifts are as bright as their local counterparts. The blue luminosity function, on the other hand, evolves rapidly with redshift, increasing its characteristic luminosity and becoming steeper at the faint end. These trends are similar to those recently observed in the Canada-France Redshift Survey.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Faint counts as a function of morphological type in a hierarchical merger model

The unprecedented resolution of the refurbished Wide Field and Planetary Camera 2 (WFPC2) on the Hubble Space Telescope (HST) has led to major advances in our understanding of galaxy formation. The high image quality in the Medium Deep Survey and Hubble Deep Field has made it possible, for the first time, to classify faint distant galaxies according to morphological type. These observations have revealed a large population of galaxies classed as irregulars or which show signs of recent merger activity. Their abundance rises steeply with apparent magnitude, providing a likely explanation for the large number of blue galaxies seen at faint magnitudes. We demonstrate that such a population arises naturally in a model in which structure forms hierarchically and which is dynamically dominated by cold dark matter. The number counts of irregular, spiral and elliptical galaxies as a function of magnitude seen in the HST data are well reproduced in this model.We present detailed predictions for the outcome of spectroscopic follow-up observations of the HST surveys. By measuring the redshift distributions of faint galaxies of different morphological types, these programmes will provide a test of the hierarchical galaxy formation paradigm and might distinguish between models with different cosmological parameters.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Large scale fluctuations in the distribution of galaxies

The null hypothesis that the three dimensional power spectrum measured from the APM Survey is consistent with the one dimensional power spectrum measured from the pencil beams surveys of Broadhurst et al and of Szalay et al is tested. The external estimates of the mean power that we make are sensitive to details of the model for the survey geometry and to the assumed level of the distortion of the pattern of galaxy clustering caused by peculiar motions. We find that the measured 3D clustering of galaxies can account for the presence of peaks in the one dimensional power spectrum, but is less successful in recovering the detailed appearence of the observations. We find no strong evidence for any additional large scale structure in the deep pencil beams beyond that recovered from the APM Survey. We conclude therefore that it is unlikey that large scale structure can be responsible for the steep local number counts of bright galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Testing ansatze for Quasi-nonlinear Clustering: the linear APM power spectrum

We compare the accuracy of published formulae that transform the linear perturbation theory power spectrum into the nonlinear regime against the results of an ensemble of large N-body simulations. We find that the modified transformation given by Jain, Mo and White (1995) performs well for initial fluctuations that have a power law power spectrum, but is less succesful for scale dependent spectra. The inverse transformation to obtain the linear power spectrum is applied to the APM Galaxy Survey power spectrum measured by Baugh \& Efstathiou (1993). The resulting linear spectrum is used to generate the initial density fluctuations in an N-body simulation, which is evolved to match the measured APM amplitude on large scales. We find very good agreement between the final power spectrum of the simulation and the measured APM power spectrum. However, the higher moments for the particle distribution only match the ones recovered from the APM Survey on large scales, $R \simgt 10 \mpc$. On small scales, $R \simlt 10 \mpc$, the APM estimations give smaller amplitudes, indicating that non-gravitational effects, such as biasing, are important on those scales. Our approach can be used to constrain a model of how light from galaxies traces the underlying mass distribution.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Galaxy formation and evolution: what to expect from hierarchical clustering models

We give a brief review of current theoretical work in galaxy formation. Recent results from N-body and N-body/hydrodynamic simulations, and from semianalytic modelling are discussed. We present updated versions of some figures from Cole et al (1994). In particular, we show the effect of using the revised stellar population synthesis model of Bruzual and Charlot, which results in a much better match to the observed colour distribution of galaxies than before. We also compare the model output with recently available data on the galaxy luminosity function and the redshift distribution of galaxies in the B and K bands. The form of the Tully-Fisher relation at high redshift predicted by our semi-analytic scheme for galaxy formation is given.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

The Real Space Correlation Function Measured from the APM Galaxy Survey

We present a determination of the real space correlation function for galaxies in the APM Survey, with magnitudes in the range $17 < b_J < 20$. We recover a power law form for the correlation function $ξ(r) = (r/4.5)^{-1.7}$ on scales $r < 4 Mpc/h$, for a flat universe and clustering that is fixed in comoving coordinates.If we assume that clustering evolves according to linear perturbation theory, we obtain $ξ(r) = (r/5.25)^{-1.7}$. We recover a shoulder feature in $ξ(r)$ for scales $4 < r < 25 Mpc/h$, with the correlation function rising above the power law, before becoming consistent with zero around $r > 40 Mpc/h$. The form we obtain for $ξ(r)$ is smooth on large scales, allowing a determination of the distortions in the correlation function in redshift space, caused by the peculiar motions of galaxies to be made. Comparing the real space correlation function with the measurement of the redshift space correlation function by Loveday et al (1995), we find that a value for $omega^{0.6}/b$ (where $ω$ is the density parameter of the universe and b is the bias factor between the density of light and mass) of unity is ruled out at the 2 sigma level.

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

A comparison of the evolution of the density field in perturbation theory and numerical simulations - II Counts in cells analysis

We present a detailed comparison of the predictions of perturbation theory for the averaged J-point correlation functions, $\xibar_J$, with the results of numerical simulations of gravitational clustering. We have carried out a systematic analysis of this method using ensembles of simulations with different numbers of particles, different box sizes and using different particle arrangements and clustering amplitudes in the initial conditions. We estimate $\xibar_J$, for $J=2-10$, from moments of counts-in-cells. We find significant non-linear effects in the variance, $J=2$, even at scales as large as $R \sim 30 \Mpc$. Perturbation theory gives remarkable agreement at large scales, where $\xibar_2 \simlt 1$, with the measured hierarchical amplitudes $S_J=\xibar_J/\xibar_2^{J-1}$. We have followed the evolution of $\xibar_J$ in time and find that for a change in the effective redshift of at least $Δz \simeq 2$ the amplitudes $S_J$ remain unchanged, despite the fact that the $\xibar_J$ have evolved by large factors, $\simeq 10^{J-1}$. We illustrate how these results can be applied to interpret the clustering in galaxy surveys and conclude that the observed hierarchical pattern in the APM is compatible with gravitational evolution in unbiased, initially Gaussian, models.

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

Hierarchical correlations in models of galaxy clustering

We present a comparison of the predictions of perturbation theory for the hierarchical J-order correlation amplitudes $S_J=\xibar_J/\xibar_2^{J-1}$, $J=3-10$, with the results of large numerical simulations of gravitational clustering. We consider two different initial power spectra of density fluctuations, one being flatter than the other. The clustering amplitudes $S_J$ measured in these models are different at all scales. In each case, the perturbation theory predictions give an excellent agreement with the simulations on scales for which the variance is approximately linear, i.e. for scales on which $\xibar_2 \simlt 1$. For cells of radius $R \simgt 30 \Mpc$, the sampling variance arising from the finite size of our simulations dominates the results. We also find that $S_J$ are roughly independent of time, $Ω$, or $λ$, as expected. On comparing with the observations, one can use these results to discriminate between models for structure formation.

0 citations0 reviewsNo code linkedOpen access