Source author record

C. S. Frenk

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

ResearcherUnclaimed source record

astro-ph.CO astro-ph astro-ph.GA hep-ph astro-ph.IM hep-th

Catalog footprint

What is connected

81works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2026arXiv

DESI Data Release 1: Stellar Catalogue

In this paper we present the stellar Value-Added Catalogue (VAC) based on the DESI Data Release 1. This VAC contains stellar parameter, abundance and radial velocity measurements for more than 4 million stars. It also contains, for the first time, measurements from individual epochs for more than a million stars with at least two observations. The main contribution to the catalogue comes from the bright program of the main survey, which includes $\sim $2.5 million stars, and the backup program, which includes $\sim $ 1 million stars. The combined magnitude range for the stars in the catalogue extends from ${\it Gaia}$ G $\sim 12$ to G $\sim 21$. For the magnitude range $17.5<G<21$ this catalogue represents a factor of 10 increase in the number of stars with radial velocity and abundance measurements compared to existing surveys. Despite DESI's resolution (R $\sim 2500-5000$), the median radial velocity uncertainty for stars in the catalogue is better than 1 km s$^{-1}$. The stellar parameters and abundances of stars in DESI are measured by two independent pipelines, and after applying a temperature-dependent calibration, [Fe/H] abundances of high signal-to-noise stars are accurate to better than $\sim$ 0.1 dex when compared to high-resolution surveys. The catalogue probes different Galactic components including a particularly large number of distant stars: tens of thousands of stars further than 10 kpc, and thousands further than 50 kpc. The catalogue also contains several thousand extremely metal-poor stars with ${\rm [Fe/H]}<-3$. The released sample of stars includes measurements for thousands of stars that are members of dwarf galaxies, open and globular clusters as well as members of several dozen stellar streams. The next public DESI data release is expected in less than two years and will contain three times as many stars as DR1.

preprint2016arXiv

Simulated Milky Way analogues: implications for dark matter indirect searches

We study high-resolution hydrodynamic simulations of Milky Way type galaxies obtained within the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) project, and identify the those that best satisfy observational constraints on the Milky Way total stellar mass, rotation curve, and galaxy shape. Contrary to mock galaxies selected on the basis of their total virial mass, the Milky Way analogues so identified consistently exhibit very similar dark matter profiles inside the solar circle, therefore enabling more accurate predictions for indirect dark matter searches. We find in particular that high resolution simulated haloes satisfying observational constraints exhibit, within the inner few kiloparsecs, dark matter profiles shallower than those required to explain the so-called Fermi GeV excess via dark matter annihilation.

preprint2015arXiv

Evolution of galaxy stellar masses and star formation rates in the EAGLE simulations

We investigate the evolution of galaxy masses and star formation rates in the Evolution and Assembly of Galaxies and their Environment (EAGLE) simulations. These comprise a suite of hydrodynamical simulations in a $Λ$CDM cosmogony with subgrid models for radiative cooling, star formation, stellar mass loss, and feedback from stars and accreting black holes. The subgrid feedback was calibrated to reproduce the observed present-day galaxy stellar mass function and galaxy sizes. Here we demonstrate that the simulations reproduce the observed growth of the stellar mass density to within 20 per cent. The simulation also tracks the observed evolution of the galaxy stellar mass function out to redshift z = 7, with differences comparable to the plausible uncertainties in the interpretation of the data. Just as with observed galaxies, the specific star formation rates of simulated galaxies are bimodal, with distinct star forming and passive sequences. The specific star formation rates of star forming galaxies are typically 0.2 to 0.4 dex lower than observed, but the evolution of the rates track the observations closely. The unprecedented level of agreement between simulation and data makes EAGLE a powerful resource to understand the physical processes that govern galaxy formation.

preprint2015arXiv

Galaxy And Mass Assembly (GAMA): end of survey report and data release 2

The Galaxy And Mass Assembly (GAMA) survey is one of the largest contemporary spectroscopic surveys of low-redshift galaxies. Covering an area of ~286 deg^2 (split among five survey regions) down to a limiting magnitude of r < 19.8 mag, we have collected spectra and reliable redshifts for 238,000 objects using the AAOmega spectrograph on the Anglo-Australian Telescope. In addition, we have assembled imaging data from a number of independent surveys in order to generate photometry spanning the wavelength range 1 nm - 1 m. Here we report on the recently completed spectroscopic survey and present a series of diagnostics to assess its final state and the quality of the redshift data. We also describe a number of survey aspects and procedures, or updates thereof, including changes to the input catalogue, redshifting and re-redshifting, and the derivation of ultraviolet, optical and near-infrared photometry. Finally, we present the second public release of GAMA data. In this release we provide input catalogue and targeting information, spectra, redshifts, ultraviolet, optical and near-infrared photometry, single-component Sérsic fits, stellar masses, H$α$-derived star formation rates, environment information, and group properties for all galaxies with r < 19.0 mag in two of our survey regions, and for all galaxies with r < 19.4 mag in a third region (72,225 objects in total). The database serving these data is available at http://www.gama-survey.org/.

preprint2015arXiv

The impact of angular momentum on black hole accretion rates in simulations of galaxy formation

Feedback from energy liberated by gas accretion onto black holes (BHs) is an attractive mechanism to explain the exponential cut-off at the massive end of the galaxy stellar mass function (SMF). Semi-analytic models of galaxy formation in which this form of feedback is assumed to suppress cooling in haloes where the gas cooling time is large compared to the dynamical time do indeed achieve a good match to the observed SMF. Furthermore, hydrodynamic simulations of individual halos in which gas is assumed to accrete onto the central BH at the Bondi rate have shown that a self-regulating regime is established in which the BH grows just enough to liberate an amount of energy comparable to the thermal energy of the halo. However, this process is efficient at suppressing the growth not only of massive galaxies but also of galaxies like the Milky Way, leading to disagreement with the observed SMF. The Bondi accretion rate, however, is inappropriate when the accreting material has angular momentum. We present an improved accretion model that takes into account the circularisation and subsequent viscous transport of infalling material and include it as a "subgrid" model in hydrodynamic simulations of the evolution of halos with a wide range of masses. The resulting accretion rates are generally low in low mass ($\lsim 10^{11.5} \msun$) halos, but show outbursts of Eddington-limited accretion during galaxy mergers. During outbursts these objects strongly resemble quasars. In higher mass haloes, gas accretion occurs continuously, typically at $~10$ % of the Eddington rate, which is conducive to the formation of radio jets. The resulting dependence of the accretion behaviour on halo mass induces a break in the relation between galaxy stellar mass and halo mass in these simulations that matches observations.

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.

preprint2013arXiv

Galaxy And Mass Assembly (GAMA): Spectroscopic analysis

The Galaxy And Mass Assembly (GAMA) survey is a multiwavelength photometric and spectroscopic survey, using the AAOmega spectrograph on the Anglo-Australian Telescope to obtain spectra for up to ~300000 galaxies over 280 square degrees, to a limiting magnitude of r_pet < 19.8 mag. The target galaxies are distributed over 0<z<0.5 with a median redshift of z~0.2, although the redshift distribution includes a small number of systems, primarily quasars, at higher redshifts, up to and beyond z=1. The redshift accuracy ranges from sigma_v~50km/s to sigma_v~100km/s depending on the signal-to-noise of the spectrum. Here we describe the GAMA spectroscopic reduction and analysis pipeline. We present the steps involved in taking the raw two-dimensional spectroscopic images through to flux-calibrated one-dimensional spectra. The resulting GAMA spectra cover an observed wavelength range of 3750<lambda<8850 A at a resolution of R~1300. The final flux calibration is typically accurate to 10-20%, although the reliability is worse at the extreme wavelength ends, and poorer in the blue than the red. We present details of the measurement of emission and absorption features in the GAMA spectra. These measurements are characterised through a variety of quality control analyses detailing the robustness and reliability of the measurements. We illustrate the quality of the measurements with a brief exploration of elementary emission line properties of the galaxies in the GAMA sample. We demonstrate the luminosity dependence of the Balmer decrement, consistent with previously published results, and explore further how Balmer decrement varies with galaxy mass and redshift. We also investigate the mass and redshift dependencies of the [NII]/Halpha vs [OIII]/Hbeta spectral diagnostic diagram, commonly used to discriminate between star forming and nuclear activity in galaxies.

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.

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.

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.

preprint2012arXiv

Galaxy and Mass Assembly (GAMA): Colour and luminosity dependent clustering from calibrated photometric redshifts

We measure the two-point angular correlation function of a sample of 4,289,223 galaxies with r < 19.4 mag from the Sloan Digital Sky Survey as a function of photometric redshift, absolute magnitude and colour down to M_r - 5log h = -14 mag. Photometric redshifts are estimated from ugriz model magnitudes and two Petrosian radii using the artificial neural network package ANNz, taking advantage of the Galaxy and Mass Assembly (GAMA) spectroscopic sample as our training set. The photometric redshifts are then used to determine absolute magnitudes and colours. For all our samples, we estimate the underlying redshift and absolute magnitude distributions using Monte-Carlo resampling. These redshift distributions are used in Limber's equation to obtain spatial correlation function parameters from power law fits to the angular correlation function. We confirm an increase in clustering strength for sub-L* red galaxies compared with ~L* red galaxies at small scales in all redshift bins, whereas for the blue population the correlation length is almost independent of luminosity for ~L* galaxies and fainter. A linear relation between relative bias and log luminosity is found to hold down to luminosities L~0.03L*. We find that the redshift dependence of the bias of the L* population can be described by the passive evolution model of Tegmark & Peebles (1998). A visual inspection of a random sample of our r < 19.4 sample of SDSS galaxies reveals that about 10 per cent are spurious, with a higher contamination rate towards very faint absolute magnitudes due to over-deblended nearby galaxies. We correct for this contamination in our clustering analysis.

preprint2012arXiv

Galaxy And Mass Assembly (GAMA): The 0.013 < z < 0.1 cosmic spectral energy distribution from 0.1 micron to 1mm

We use the GAMA I dataset combined with GALEX, SDSS and UKIDSS imaging to construct the low-redshift (z<0.1) galaxy luminosity functions in FUV, NUV, ugriz, and YJHK bands from within a single well constrained volume of 3.4 x 10^5 (Mpc/h)^{3}. The derived luminosity distributions are normalised to the SDSS DR7 main survey to reduce the estimated cosmic variance to the 5 per cent level. The data are used to construct the cosmic spectral energy distribution (CSED) from 0.1 to 2.1 \mum free from any wavelength dependent cosmic variance for both the elliptical and non-elliptical populations. The two populations exhibit dramatically different CSEDs as expected for a predominantly old and young population respectively. Using the Driver et al. (2008) prescription for the azimuthally averaged photon escape fraction, the non-ellipticals are corrected for the impact of dust attenuation and the combined CSED constructed. The final results show that the Universe is currently generating (1.8 +/- 0.3) x 10^{35} h W Mpc^{-3} of which (1.2 +/- 0.1) x 10^{35} h W Mpc^{-3} is directly released into the inter-galactic medium and (0.6 +/- 0.1) x 10^{35} h W Mpc^{-3} is reprocessed and reradiated by dust in the far-IR. Using the GAMA data and our dust model we predict the mid and far-IR emission which agrees remarkably well with available data. We therefore provide a robust description of the pre- and post dust attenuated energy output of the nearby Universe from 0.1micron to 0.6mm. The largest uncertainty in this measurement lies in the mid and far-IR bands stemming from the dust attenuation correction and its currently poorly constrained dependence on environment, stellar mass, and morphology.

preprint2012arXiv

Galaxy And Mass Assembly (GAMA): The mass-metallicity relationship

Context: The mass-metallicity relationship (MMR) of star-forming galaxies is well-established, however there is still some disagreement with respect to its exact shape and its possible dependence on other observables. Aims: We measure the MMR in the Galaxy And Mass Assembly (GAMA) survey. We compare our measured MMR to that measured in the Sloan Digital Sky Survey (SDSS) and study the dependence of the MMR on various selection criteria to identify potential causes for disparities seen in the literature. Methods: We use strong emission line ratio diagnostics to derive oxygen abundances. We then apply a range of selection criteria for the minimum signal-to-noise in various emission lines, as well as the apparent and absolute magnitude to study variations in the inferred MMR. Results: The shape and position of the MMR can differ significantly depending on the metallicity calibration and selection used. After selecting a robust metallicity calibration amongst those tested, we find that the mass-metallicity relation for redshifts 0.061< z<0.35 in GAMA is in reasonable agreement with that found in the SDSS despite the difference in the luminosity range probed. Conclusions: In view of the significant variations of the MMR brought about by reasonable changes in the sample selection criteria and method, we recommend that care be taken when comparing the MMR from different surveys and studies directly. We also conclude that there could be a modest level of evolution over 0.06<z<0.35 within the GAMA sample.

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.

preprint2012arXiv

Principles of supernova-driven winds

The formation of galaxies is regulated by a balance between the supply of gas and the rate at which it is ejected. Traditional explanations of gas ejection equate the energy required to escape the galaxy or host halo to an estimate for the energy yield from supernovae. This yield is usually assumed to be a constant fraction of the total available from the supernova, or is derived from the assumption of a consistent momentum yield. By applying these ideas in the context of a cold dark matter cosmogony, we derive a 1st-order analytic connection between these working assumptions and the expected relationship between baryon content and galaxy circular velocity, and find that these quick predictions straddle recent observational estimates. To examine the premises behind these theories in more detail, we then explore their applicability to a set of gasdynamical simulations of idealised galaxies. We show that different premises dominate to differing degrees in the simulated outflow, depending on the mass of the system and the resolution with which it is simulated. Using this study to anticipate the emergent behaviour at arbitrarily high resolution, we motivate more comprehensive analytic model which allows for the range of velocities with which the gas may exit the system, and incorporates both momentum and energy-based constraints on the outflow. Using a trial exit velocity distribution, this is shown to be compatible with the observed baryon fractions in intermediate-mass systems, but implies that current estimates for low-mass systems can not be solely accounted for by supernova winds under commonly-held assumptions.

preprint2012arXiv

The Aquila comparison Project: The Effects of Feedback and Numerical Methods on Simulations of Galaxy Formation

We compare the results of thirteen cosmological gasdynamical codes used to simulate the formation of a galaxy in the LCDM structure formation paradigm. The various runs differ in their hydrodynamical treatment (SPH, moving-mesh and AMR) but share the same initial conditions and adopt their latest published model of cooling, star formation and feedback. Despite the common halo assembly history, we find large code-to-code variations in the stellar mass, size, morphology and gas content of the galaxy at z=0, due mainly to the different implementations of feedback. Compared with observation, most codes tend to produce an overly massive galaxy, smaller and less gas-rich than typical spirals, with a massive bulge and a declining rotation curve. A stellar disk is discernible in most simulations, though its prominence varies widely from code to code. There is a well-defined trend between the effects of feedback and the severity of the disagreement with observation. Models that are more effective at limiting the baryonic mass of the galaxy come closer to matching observed galaxy scaling laws, but often to the detriment of the disk component. Our conclusions hold at two different numerical resolutions. Some differences can also be traced to the numerical techniques: more gas seems able to cool and become available for star formation in grid-based codes than in SPH. However, this effect is small compared to the variations induced by different feedback prescriptions. We conclude that state-of-the-art simulations cannot yet uniquely predict the properties of the baryonic component of a galaxy, even when the assembly history of its host halo is fully specified. Developing feedback algorithms that can effectively regulate the mass of a galaxy without hindering the formation of high-angular momentum stellar disks remains a challenge.

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.

preprint2012arXiv

The Phoenix Project: the Dark Side of Rich Galaxy Clusters

[abridged] We introduce the Phoenix Project, a set of $Λ$CDM simulations of the dark matter component of nine rich galaxy clusters. Each cluster is simulated at least at two different numerical resolutions. For eight of them, the highest resolution corresponds to $\sim 130$ million particles within the virial radius, while for one this number is over one billion. We study the structure and substructure of these systems and contrast them with six galaxy-sized dark matter haloes from the Aquarius Project, simulated at comparable resolution. This comparison highlights the approximate mass invariance of CDM halo structure and substructure. We find little difference in the spherically-averaged mass, pseudo-phase-space density, and velocity anisotropy profiles of Aquarius and Phoenix haloes. When scaled to the virial properties of the host halo, the abundance and radial distribution of subhaloes are also very similar, despite the fact that Aquarius and Phoenix haloes differ by roughly three decades in virial mass. The most notable difference is that cluster haloes have been assembled more recently and are thus significantly less relaxed than galaxy haloes, which leads to decreased regularity, increased halo-to-halo scatter and sizable deviations from the mean trends. This accentuates the effects of the strong asphericity of individual clusters on surface density profiles, which may vary by up to a factor of three at a given radius, depending on projection. The high apparent concentration reported for some strong-lensing clusters might very well reflect these effects. A more recent assembly also explains why substructure in some Phoenix haloes is slightly more abundant than in Aquarius, especially in the inner regions. Resolved subhaloes nevertheless contribute only $11 \pm 3%$ of the virial mass in Phoenix clusters. .

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.

preprint2011arXiv

A Two-Point Correlation Function For Galactic Halo Stars

We describe a correlation function statistic that quantifies the amount of spatial and kinematic substructure in the stellar halo. We test this statistic using model stellar halo realizations constructed from the Aquarius suite of six high-resolution N-body simulations in combination with the Galform semi-analytic galaxy formation model. These simulations show considerable scatter in the properties of stellar haloes. We find that our statistic can distinguish between these plausible alternatives for the global structure of the Milky Way stellar halo. We compare with observational data and show that pencil beam surveys of ~100 tracer stars (such as the Spaghetti Survey) are not sufficient to constrain the degree of structure in the Milky Way halo with this statistic. Larger area surveys with >1000 tracer stars (such as BHB stars in the Sloan Digital Sky Survey) provide much tighter constraints on comparisons between models and data. In our simulations, we find examples of haloes with spatial and kinematic substructure consistent with the available Milky Way data.

preprint2011arXiv

Galaxy and Mass Assembly (GAMA): The GAMA Galaxy Group Catalogue (G3Cv1)

Using the complete GAMA-I survey covering ~142 sq. deg. to r=19.4, of which ~47 sq. deg. is to r=19.8, we create the GAMA-I galaxy group catalogue (G3Cv1), generated using a friends-of-friends (FoF) based grouping algorithm. Our algorithm has been tested extensively on one family of mock GAMA lightcones, constructed from Lambda-CDM N-body simulations populated with semi-analytic galaxies. Recovered group properties are robust to the effects of interlopers and are median unbiased in the most important respects. G3Cv1 contains 14,388 galaxy groups (with multiplicity >= 2$), including 44,186 galaxies out of a possible 110,192 galaxies, implying ~40% of all galaxies are assigned to a group. The similarities of the mock group catalogues and G3Cv1 are multiple: global characteristics are in general well recovered. However, we do find a noticeable deficit in the number of high multiplicity groups in GAMA compared to the mocks. Additionally, despite exceptionally good local spatial completeness, G3Cv1 contains significantly fewer compact groups with 5 or more members, this effect becoming most evident for high multiplicity systems. These two differences are most likely due to limitations in the physics included of the current GAMA lightcone mock. Further studies using a variety of galaxy formation models are required to confirm their exact origin.

preprint2011arXiv

Galaxy And Mass Assembly (GAMA): the red fraction and radial distribution of satellite galaxies

We investigate the properties of satellite galaxies that surround isolated hosts within the redshift range 0.01 < z < 0.15, using data taken as part of the Galaxy And Mass Assembly survey. Making use of isolation and satellite criteria that take into account stellar mass estimates, we find 3514 isolated galaxies of which 1426 host a total of 2998 satellites. Separating the red and blue populations of satellites and hosts, using colour-mass diagrams, we investigate the radial distribution of satellite galaxies and determine how the red fraction of satellites varies as a function of satellite mass, host mass and the projected distance from their host. Comparing the red fraction of satellites to a control sample of small neighbours at greater projected radii, we show that the increase in red fraction is primarily a function of host mass. The satellite red fraction is about 0.2 higher than the control sample for hosts with 11.0 < log M_* < 11.5, while the red fractions show no difference for hosts with 10.0 < log M_* < 10.5. For the satellites of more massive hosts the red fraction also increases as a function of decreasing projected distance. Our results suggest that the likely main mechanism for the quenching of star formation in satellites hosted by isolated galaxies is strangulation.

preprint2011arXiv

Galaxy And Mass Assembly (GAMA): The star formation rate dependence of the stellar initial mass function

The stellar initial mass function (IMF) describes the distribution in stellar masses produced from a burst of star formation. For more than fifty years, the implicit assumption underpinning most areas of research involving the IMF has been that it is universal, regardless of time and environment. We measure the high-mass IMF slope for a sample of low-to-moderate redshift galaxies from the Galaxy And Mass Assembly survey. The large range in luminosities and galaxy masses of the sample permits the exploration of underlying IMF dependencies. A strong IMF-star formation rate dependency is discovered, which shows that highly star forming galaxies form proportionally more massive stars (they have IMFs with flatter power-law slopes) than galaxies with low star formation rates. This has a significant impact on a wide variety of galaxy evolution studies, all of which rely on assumptions about the slope of the IMF. Our result is supported by, and provides an explanation for, the results of numerous recent explorations suggesting a variation of or evolution in the IMF.

preprint2011arXiv

Galaxy and Mass Assembly (GAMA): ugriz galaxy luminosity functions

Galaxy and Mass Assembly (GAMA) is a project to study galaxy formation and evolution, combining imaging data from ultraviolet to radio with spectroscopic data from the AAOmega spectrograph on the Anglo-Australian Telescope. Using data from phase 1 of GAMA, taken over three observing seasons, and correcting for various minor sources of incompleteness, we calculate galaxy luminosity functions (LFs) and their evolution in the ugriz passbands. At low redshift, z < 0.1, we find that blue galaxies, defined according to a magnitude-dependent but non-evolving colour cut, are reasonably well fit over a range of more than ten magnitudes by simple Schechter functions in all bands. Red galaxies, and the combined blue-plus-red sample, require double power-law Schechter functions to fit a dip in their LF faintward of the characteristic magnitude M* before a steepening faint end. This upturn is at least partly due to dust-reddened disk galaxies. We measure evolution of the galaxy LF over the redshift range 0.002 < z < 0.5 using both a parametric fit and by measuring binned LFs in redshift slices. The characteristic luminosity L* is found to increase with redshift in all bands, with red galaxies showing stronger luminosity evolution than blue galaxies. The comoving number density of blue galaxies increases with redshift, while that of red galaxies decreases, consistent with prevailing movement from blue cloud to red sequence. As well as being more numerous at higher redshift, blue galaxies also dominate the overall luminosity density beyond redshifts z = 0.2. At lower redshifts, the luminosity density is dominated by red galaxies in the riz bands, by blue galaxies in u and g.

preprint2011arXiv

Mismatch and Misalignment: Dark Haloes and Satellites of Disc Galaxies

(Abridged) We study the phase-space distribution of satellite galaxies associated with late-type galaxies in the GIMIC suite of simulations. GIMIC consists of re-simulations of 5 cosmologically representative regions from the Millennium simulation, which have higher resolution and incorporate baryonic physics. Whilst the disc of the galaxy is well aligned with the inner regions (r ~ 0.1r_200) of the dark matter halo, both in shape and angular momentum, there can be substantial misalignments at larger radii (r ~ r_200). Misalignments of > 45 deg are seen in ~ 30% of our sample. We find that the satellite population aligns with the shape (and angular momentum) of the outer dark matter halo. However, the alignment with the galaxy is weak owing to the mismatch between the disc and dark matter halo. Roughly 20% of the satellite systems with ten bright galaxies within r_200 exhibit a polar spatial alignment with respect to the galaxy --- an orientation reminiscent of the classical satellites of the Milky Way. We find that a small fraction (~10%) of satellite systems show evidence for rotational support which we attribute to group infall. There is a bias towards satellites on prograde orbits relative to the spin of the dark matter halo (and to a lesser extent with the angular momentum of the disc). This preference towards co-rotation is stronger in the inner regions of the halo where the most massive satellites accreted at relatively early times are located. We attribute the anisotropic spatial distribution and angular momentum bias of the satellites at z=0 to their directional accretion along the major axes of the dark matter halo. The satellite galaxies have been accreted relatively recently compared to the dark matter mass and have experienced less phase-mixing and relaxation --- the memory of their accretion history can remain intact to z=0.

preprint2011arXiv

Substructure in the stellar halos of the Aquarius simulations

We characterize substructure in the simulated stellar halos of Cooper et al. (2010) which were formed by the disruption of satellite galaxies within the cosmological N-body simulations of galactic halos of the Aquarius Project. These stellar halos exhibit a wealth of tidal features: broad overdensities and very narrow faint streams akin to those observed around the Milky Way. The substructures are distributed anisotropically on the sky, a characteristic that should become apparent in the next generation of photometric surveys. The normalized RMS of the density of stars on the sky appears to be systematically larger for our halos compared to the value estimated for the Milky Way from main sequence turn-off stars in the Sloan Digital Sky Survey. We show that this is likely to be due in part to contamination by faint QSOs and redder main sequence stars, and might suggest that ~10% of the Milky Way halo stars have formed in-situ.

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).

preprint2011arXiv

Where will supersymmetric dark matter first be seen?

If the dark matter consists of supersymmetric particles, Gamma-ray observatories such as the Large Area Telescope aboard the Fermi satellite may detect annihilation radiation from the haloes of galaxies and galaxy clusters. Much recent effort has been devoted to searching for this signal around the Milky Way's dwarf satellites. Using a new suite of high-resolution simulations of galaxy cluster haloes (the Phoenix Project), together with the Aquarius simulations of Milky-Way-like galaxy haloes, we show that higher signal-to-noise and equally clean signals are, in fact, predicted to come from nearby rich galaxy clusters. Most of the cluster emission is produced by small subhaloes with masses less than that of the Sun. The large range of mass scales covered by our two sets of simulations allows us to deduce a physically motivated extrapolation to these small (and unresolved) masses. Since tidal effects destroy subhaloes in the dense inner regions of haloes, most cluster emission is then predicted to come from large radii, implying that the nearest and brightest systems should be much more extended than Fermi's angular resolution limit. The most promising targets for detection are clusters such as Coma and Fornax, but detection algorithms must be tuned to the predicted profile of the emission if they are to maximize the chance of finding this weak signal.

preprint2010arXiv

Assembly History and Structure of Galactic Cold Dark Matter Halos

We use the Aquarius simulation series to study the imprint of assembly history on the structure of Galaxy-mass cold dark matter halos. Our results confirm earlier work regarding the influence of mergers on the mass density profile and the inside-out growth of halos. The inner regions that contain the visible galaxies are stable since early times and are significantly affected only by major mergers. Particles accreted diffusely or in minor mergers are found predominantly in the outskirts of halos. Our analysis reveals trends that run counter to current perceptions of hierarchical halo assembly. For example, major mergers (i.e. those with progenitor mass ratios greater than 1:10) contribute little to the total mass growth of a halo, on average less than 20 per cent for our six Aquarius halos. The bulk is contributed roughly equally by minor mergers and by "diffuse" material which is not resolved into individual objects. This is consistent with modeling based on excursion-set theory which suggests that about half of this diffuse material should not be part of a halo of any scale. Interestingly, the simulations themselves suggest that a significantly fraction is not truly diffuse, since it was ejected from earlier halos by mergers prior to their joining the main system. The Aquarius simulations resolve halos to much lower mass scales than are expected to retain gas or form stars. These results thus confirm that most of the baryons from which visible galaxies form are accreted diffusely, rather than through mergers, and they suggest that only relatively rare major mergers will affect galaxy structure at later times.

preprint2010arXiv

Galactic stellar haloes in the CDM model

We present six simulations of galactic stellar haloes formed by the tidal disruption of accreted dwarf galaxies in a fully cosmological setting. Our model is based on the Aquarius project, a suite of high resolution N-body simulations of individual dark matter haloes. We tag subsets of particles in these simulations with stellar populations predicted by the Galform semi-analytic model. Our method self-consistently tracks the dynamical evolution and disruption of satellites from high redshift. The luminosity function and structural properties of surviving satellites, which agree well with observations, suggest that this technique is appropriate. We find that accreted stellar haloes are assembled between 1<z<7 from less than 5 significant progenitors. These progenitors are old, metal-rich satellites with stellar masses similar to the brightest Milky Way dwarf spheroidals. In contrast to previous stellar halo simulations, we find that several of these major contributors survive as self-bound systems to the present day. Both the number of these significant progenitors and their infall times are inherently stochastic. This results in great diversity among our stellar haloes, which amplifies small differences between the formation histories of their dark halo hosts. The masses and density/surface-brightness profiles of the stellar haloes are consistent with expectations from the Milky Way and M31. Each halo has a complex structure, consisting of well-mixed components, tidal streams, shells and other subcomponents. This structure is not adequately described by smooth models. We find one example of an accreted thick disk. Contrasts in age and metallicity between halo stars and those in surviving satellites are in broad agreement with recent observations. [Abridged]

preprint2010arXiv

Galaxy And Mass Assembly (GAMA): Galaxies at the faint end of the Halpha luminosity function

We present an analysis of the properties of the lowest Halpha-luminosity galaxies (L_Halpha<4x10^32 W; SFR<0.02 Msun/yr) in the Galaxy And Mass Assembly (GAMA) survey. These galaxies make up the the rise above a Schechter function in the number density of systems seen at the faint end of the Halpha luminosity function. Above our flux limit we find that these galaxies are principally composed of intrinsically low stellar mass systems (median stellar mass =2.5x10^8 Msun) with only 5/90 having stellar masses M>10^10 Msun. The low SFR systems are found to exist predominantly in the lowest density environments (median density ~0.02 galaxy Mpc^-2 with none in environments more dense than ~1.5 galaxy Mpc^-2). Their current specific star formation rates (SSFR; -8.5 < log(SSFR[yr^-1])<-12.) are consistent with their having had a variety of star formation histories. The low density environments of these galaxies demonstrates that such low-mass, star-forming systems can only remain as low-mass and forming stars if they reside sufficiently far from other galaxies to avoid being accreted, dispersed through tidal effects or having their gas reservoirs rendered ineffective through external processes.

preprint2010arXiv

Galaxy and Mass Assembly (GAMA): survey diagnostics and core data release

The Galaxy And Mass Assembly (GAMA) survey has been operating since February 2008 on the 3.9-m Anglo-Australian Telescope using the AAOmega fibre-fed spectrograph facility to acquire spectra with a resolution of R~1300 for 120,862 SDSS selected galaxies. The target catalogue constitutes three contiguous equatorial regions centred at 9h (G09), 12h (G12) and 14.5h (G15) each of 12 x 4 sq.deg to limiting fluxes of r < 19.4, r < 19.8, and r < 19.4 mag respectively (and additional limits at other wavelengths). Spectra and reliable redshifts have been acquired for over 98 per cent of the galaxies within these limits. Here we present the survey footprint, progression, data reduction, redshifting, re-redshifting, an assessment of data quality after 3 years, additional image analysis products (including ugrizYJHK photometry, Sersic profiles and photometric redshifts), observing mask, and construction of our core survey catalogue (GamaCore). From this we create three science ready catalogues: GamaCoreDR1 for public release, which includes data acquired during year 1 of operations within specified magnitude limits (February 2008 to April 2008); GamaCoreMainSurvey containing all data above our survey limits for use by the GAMA team and collaborators; and GamaCoreAtlasSv containing year 1, 2 and 3 data matched to Herschel-ATLAS Science Demonstration data. These catalogues along with the associated spectra, stamps and profiles can be accessed via the GAMA website: http://www.gama-survey.org/

preprint2010arXiv

Galaxy And Mass Assembly (GAMA): The input catalogue and star-galaxy separation

We describe the spectroscopic target selection for the Galaxy And Mass Assembly (GAMA) survey. The input catalogue is drawn from the Sloan Digital Sky Survey (SDSS) and UKIRT Infrared Deep Sky Survey (UKIDSS). The aim is to measure redshifts for galaxies in three 4x12 degree regions at 9h, 12h and 14.5h, on the celestial equator, with magnitude selections r < 19.4, z < 18.2 and K(AB) < 17.6 over all three regions, and r < 19.8 in the 12-h region. The target density is 1080 per sq. deg. in the 12-h region and 720 per sq. deg. in the other regions. The average GAMA target density and area are compared with completed and ongoing galaxy redshift surveys. The GAMA survey implements a highly complete star-galaxy separation that jointly uses an intensity-profile separator (delta(sg) = r-band psf mag - model mag) as per the SDSS) and a colour separator. The colour separator is defined as delta(sg,jk) = J - K - f(g-i), where f(g-i) is a quadratic fit to the J-K colour of the stellar locus over the range 0.3 < g-i < 2.3. All galaxy populations investigated are well separated with delta(sg,jk) > 0.2. From two years out of a three-year AAOmega program on the Anglo-Australian Telescope, we have obtained 79599 unique galaxy redshifts. Previously known redshifts in the GAMA region bring the total up to 98497. The median galaxy redshift is 0.2 with 99% at z < 0.5. We present some of the global statistical properties of the survey, including K-band galaxy counts, colour-redshift relations and preliminary n(z).

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.

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.

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.

preprint2010arXiv

Substructure Lensing: Effects of Galaxies, Globular Clusters & Satellite Streams

Lensing flux-ratio anomalies have been frequently observed and taken as evidence for the presence of abundant dark matter substructures in lensing galaxies, as predicted by the cold dark matter (CDM) model of cosmogony. In previous work, we examined the cusp-caustic relations of the multiple images of background quasars lensed by galaxy-scale dark matter haloes, using a suite of high-resolution N-body simulations (the Aquarius simulations). In this work, we extend our previous calculations to incorporate both the baryonic and diffuse dark components in lensing haloes. We include in each lensing simulation: (1) a satellite galaxy population derived from a semi-analytic model applied to the Aquarius haloes, (2) an empirical Milky-Way globular cluster population and (3) satellite streams (diffuse dark component) identified in the simulations. Accounting for these extra components, we confirm our earlier conclusion that the abundance of intrinsic substructures (dark or bright, bound or diffuse) is not sufficient to explain the observed frequency of cusp-caustic violations in the CLASS survey. We conclude that the observed effect could be the result of the small number statistics of CLASS, or intergalactic haloes along the line of sight acting as additional sources of lensing flux anomalies. Another possibility is that this discrepancy signals a failure of the CDM model.

preprint2010arXiv

The earliest stars and their relics in the Milky Way

We have implemented a simple model to identify the likely sites of the first stars and galaxies in the high-resolution simulations of the formation of galactic dark matter halos of the Aquarius Project. The first star in a galaxy like the Milky Way formed around redshift z=35; by z=10, the young galaxy contained up to 30000 dark matter haloes capable of forming stars by molecular hydrogen cooling. These minihaloes were strongly clustered and feedback may have severely limited the actual number of Population III stars that formed. By the present day, the remnants of the first stars would be strongly concentrated to the centre of the main halo. If a second generation of long-lived stars formed near the first (the first star relics), we would expect to find half of them within 30kpc/h of the galactic centre and a significant fraction in satellites where they may be singled out by their anomalous metallicity patterns. The first halo in which gas could cool by atomic hydrogen line radiation formed at z=25; by z=10, the number of such first galaxies had increased to ~300. Feedback might have decreased the number of first galaxies at the time they undergo significant star formation, but not the number that survive to the present because near neighbours merge. Half of all the "first galaxies" that form before z=10 merge with the main halo before z~3 and most lose a significant fraction of their mass. However, today there should still be more than 20 remnants orbiting within the central ~30kpc/h of the Milky Way. These satellites have circular velocities of a few kilometers per second or more, comparable to those of known Milky Way dwarfs. They are a promising hunting ground for the remnants of the earliest epoch of star formation.

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.

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.

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.

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.

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.

preprint2010arXiv

The statistics of the subhalo abundance of dark matter haloes

We study the population statistics of the surviving subhaloes of LCDM dark matter haloes using a set of very high resolution N-body simulations. These include both simulations of representative regions of the Universe and ultra-high resolution resimulations of individual dark matter haloes. We find that more massive haloes tend to have a larger mass fraction in subhaloes. For example, cluster size haloes typically have 7.5 percent of their mass within R200 in substructures of fractional mass larger than 1e-5, which is 25 percent higher than galactic haloes. There is, however, a large variance in the subhalo mass fraction from halo to halo, whereas the subhalo abundance shows much higher regularity. For dark matter haloes of fixed mass, the subhalo abundance decreases by 30 percent between redshift 2 and 0. The subhalo abundance function correlates with the host halo concentration parameter and formation redshift. However, the intrinsic scatter is not significantly reduced for narrow ranges of concentration parameter or formation redshift, showing that they are not the dominant parameters that determine the subhalo abundance in a halo.

preprint2009arXiv

Effects of Dark Matter Substructures on Gravitational Lensing: Results from the Aquarius Simulations

We use high-resolution Aquarius simulations of Milky Way-sized haloes in the LCDM cosmology to study the effects of dark matter substructures on gravitational lensing. Each halo is resolved with ~ 10^8 particles (at a mass resolution ~ 10^3-4 M_sun/h) within its virial radius. Subhaloes with masses larger than 10^5 M_sun/h are well resolved, an improvement of at least two orders of magnitude over previous lensing studies. We incorporate a baryonic component modelled as a Hernquist profile and account for the response of the dark matter via adiabatic contraction. We focus on the "anomalous" flux ratio problem, in particular on the violation of the cusp-caustic relation due to substructures. We find that subhaloes with masses less than ~ 10^8 M_sun/h play an important role in causing flux anomalies; such low mass subhaloes have been unresolved in previous studies. There is large scatter in the predicted flux ratios between different haloes and between different projections of the same halo. In some cases, the frequency of predicted anomalous flux ratios is comparable to that observed for the radio lenses, although in most cases it is not. The probability for the simulations to reproduce the observed violations of the cusp lenses is about 0.001. We therefore conclude that the amount of substructure in the central regions of the Aquarius haloes is insufficient to explain the observed frequency of violations of the cusp-caustic relation. These conclusions are based purely on our dark matter simulations which ignore the effect of baryons on subhalo survivability.

preprint2009arXiv

Galaxy and Mass Assembly (GAMA): Optimal Tiling of Dense Surveys with a Multi-Object Spectrograph

A heuristic greedy algorithm is developed for efficiently tiling spatially dense redshift surveys. In its first application to the Galaxy and Mass Assembly (GAMA) redshift survey we find it rapidly improves the spatial uniformity of our data, and naturally corrects for any spatial bias introduced by the 2dF multi object spectrograph. We make conservative predictions for the final state of the GAMA redshift survey after our final allocation of time, and can be confident that even if worse than typical weather affects our observations, all of our main survey requirements will be met.

preprint2009arXiv

Galaxy morphology in the LambdaCDM cosmology

We investigate the origins of galaxy morphology (defined by bulge-to-total K-band luminosity) in the LambdaCDM cosmology using two galaxy formation models ,based on the Millennium simulation, by Bower et al. (the Durham model) and De Lucia & Blaizot (the MPA model). There are many similarities, but also fundamental disagreements in the predictions of the two models for galaxy morphology. For example, taking into account uncertainties in the available observational data, both produce a realistic morphological mix today, but its evolution is very different. A main cause of this and other differences is the treatment of disk instabilities which play a more prominent role in the Durham model. Our analysis confirms previous theoretical predictions that elliptical galaxies form most of their stars before the bulk of the galaxy is assembled. Spirals tend to have later `assembly' times as a consequence of in-situ star formation. With the exception of the brightest ellipticals (stellar mass > 2.5e11 M_sun/h), we find that major mergers are not the primary mechanism by which most spheroids (ellipticals and spiral bulges) assemble their mass. In fact, the majority of ellipticals (and the overwhelming majority of spirals) never experience a major merger (above the resolution limit of our simulation.) Most ellipticals and spiral bulges acquire their stellar mass through minor mergers or disk instabilities. These conclusions are common to both models. The rotation properties of spheroids may help to constrain the importance of disk instabilities in these models.

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.

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.

preprint2009arXiv

The role of black holes in galaxy formation and evolution

Virtually all massive galaxies, including our own, host central black holes ranging in mass from millions to billions of solar masses. The growth of these black holes releases vast amounts of energy that powers quasars and other weaker active galactic nuclei. A tiny fraction of this energy, if absorbed by the host galaxy, could halt star formation by heating and ejecting ambient gas. A central question in galaxy evolution is the degree to which this process has caused the decline of star formation in large elliptical galaxies, which typically have little cold gas and few young stars, unlike spiral galaxies.

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.

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.

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.

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).

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.

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.

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.

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.

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.

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.

preprint2003arXiv

Stable clustering, the halo model and nonlinear cosmological power spectra

We present the results of a large library of cosmological N-body simulations, using power-law initial spectra. The nonlinear evolution of the matter power spectra is compared with the predictions of existing analytic scaling formulae based on the work of Hamilton et al. The scaling approach has assumed that highly nonlinear structures obey `stable clustering' and are frozen in proper coordinates. Our results show that, when transformed under the self-similarity scaling, the scale-free spectra define a nonlinear locus that is clearly shallower than would be required under stable clustering. Furthermore, the small-scale nonlinear power increases as both the power-spectrum index n and the density parameter Omega decrease, and this evolution is not well accounted for by the previous scaling formulae. This breakdown of stable clustering can be understood as resulting from the modification of dark-matter haloes by continuing mergers. These effects are naturally included in the analytic `halo model' for nonlinear structure; using this approach we are able to fit both our scale-free results and also our previous CDM data. This approach is more accurate than the commonly-used Peacock--Dodds formula and should be applicable to more general power spectra. Code to evaluate nonlinear power spectra using this method is available from http://as1.chem.nottingham.ac.uk/~res/software.html Following publication, we will make the power-law simulation data available through the Virgo website http://www.mpa-garching.mpg.de/Virgo

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.

preprint2002arXiv

Radio sources in the 2dF Galaxy Redshift Survey - II. Local radio luminosity functions for AGN and star-forming galaxies at 1.4 GHz

We have cross-matched the 1.4 GHz NRAO VLA Sky Survey (NVSS) with the first 210 fields observed in the 2dF Galaxy Redshift Survey (2dFGRS), covering an effective area of 325 square degrees (about 20% of the final 2dFGRS area). This yields a set of optical spectra of 912 candidate NVSS counterparts, of which we identify 757 as genuine radio IDs - the largest and most homogeneous set of radio-source spectra ever obtained. The 2dFGRS radio sources span the redshift range z=0.005 to 0.438, and are a mixture of active galaxies (60%) and star-forming galaxies (40%). About 25% of the 2dFGRS radio sources are spatially resolved by NVSS, and the sample includes three giant radio galaxies with projected linear size greater than 1 Mpc. The high quality of the 2dF spectra means we can usually distinguish unambiguously between AGN and star-forming galaxies. We have made a new determination of the local radio luminosity function at 1.4 GHz for both active and star-forming galaxies, and derive a local star-formation density of 0.022+/-0.004 solar masses per year per cubic Mpc. (Ho=50 km/s/Mpc).

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.

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.

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.

preprint2002arXiv

The 2dF Galaxy Redshift Survey: the luminosity function of cluster galaxies

We have determined the composite luminosity function (LF) for galaxies in 60 clusters from the 2dF Galaxy Redshift Survey. The LF spans the range $-22.5<M_{b_{\rm J}}<-15$, and is well-fitted by a Schechter function with ${M_{b_{\rm J}}}^{*}=-20.07\pm0.07$ and $α=-1.28\pm0.03$ ($H_0$=100 km s$^{-1}$ Mpc$^{-1}$, $Ω_M$=0.3, $Ω_Λ$=0.7). It differs significantly from the field LF of \cite{mad02}, having a characteristic magnitude that is approximately 0.3 mag brighter and a faint-end slope that is approximately 0.1 steeper. There is no evidence for variations in the LF across a wide range of cluster properties. However the LF of early-type galaxies in clusters is both brighter and steeper than its field counterpart. The differences between the field and cluster LFs for the various spectral types can be qualitatively explained by the suppression of star formation in the dense cluster environment, together with mergers to produce the brightest early-type galaxies.

preprint2002arXiv

The 2dF Galaxy Redshift Survey: The population of nearby radio galaxies at the 1 mJy level

We use redshift determinations and spectral analysis of galaxies in the 2dF Galaxy Redshift Survey to study the properties of local radio sources with $S\ge 1$ mJy. 557 objects drawn from the FIRST survey, corresponding to 2.3 per cent of the total radio sample, are found in the 2dFGRS cataloge within the area $9^h 48^m \simlt {\rm RA}({\rm 2000}) \simlt 14^h 32^m$ and $-2.77^\circ \simlt {\rm dec}({\rm 2000}) \simlt 2.25^\circ$, down to a magnitude limit $b_J= 19.45$. The excellent quality of 2dF spectra allows us to divide these sources into classes, according to their optical spectra. Radio galaxies (i.e. absorption systems, mainly belonging to the class of FR I sources, whose spectra may or may not show emission lines due to AGN activity) make up 63 per cent of the spectroscopic sample; they are characterized by relatively high radio-to-optical ratios, red colours and luminosities $10^{21}\simlt P_{1.4{\rm GHz}} W Hz^${-1}$ sr$^{-1}$ \simlt 10^{24}$. Actively star-forming galaxies instead contribute with 30 per cent to the total sample. These objects -- principally found at low redshifts ($z \simlt 0.1$) -- in general show low radio-to-optical ratios, blue colours and low radio luminosities. We also found 18 Seyfert 2 galaxies and 4 Seyfert 1's. Analysis of the local radio luminosity function shows that radio galaxies are well described by models that assume pure luminosity evolution, at least down to radio powers $P_{1.4{\rm GHz}}\simlt 10^{20.5} W Hz$^{-1}$ sr$^{-1}$ Late-type galaxies, whose relative contribution to the RLF is found to be lower than was predicted by previous work, present a luminosity function which is compatible with the IRAS galaxy LF, and therefore plausibly constitute the radio counterpart of the dusty spirals and starbursts that dominate the counts at 60 $μ$m.

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.

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.

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.

preprint1998arXiv

Evolution of structure in cold dark matter universes

We present an analysis of the clustering evolution of dark matter in four cold dark matter (CDM) cosmologies. We use a suite of high resolution, 17-million particle, N-body simulations which sample volumes large enough to give clustering statistics with unprecedented accuracy. We investigate both a flat and an open model with Omega_0=0.3, and two models with Omega=1, one with the standard CDM power spectrum and the other with the same power spectrum as the Omega_0=0.3 models. The amplitude of primordial fluctuations is set so that the models reproduce the observed abundance of rich galaxy clusters by the present day. The mass 2-point correlation function and power spectrum of all the simulations differ significantly from those of the observed galaxy distribution, in both shape and amplitude. Thus, for any of these models to provide an acceptable representation of reality, the distribution of galaxies must be biased relative to the mass in a non-trivial, scale-dependent, fashion. In the Omega=1 models the required bias is always greater than unity, but in the Omega_0=0.3 models an "antibias" is required on scales smaller than \sim 5\hmpc. The mass correlation functions in the simulations are well fit by recently published analytic models. The velocity fields are remarkably similar in all the models, whether they be characterised as bulk flows, 1-particle or pairwise velocity dispersions. This similarity is a direct consequence of our adopted normalisation. The small-scale pairwise velocity dispersion of the dark matter is somewhat larger than recent determinations from galaxy redshift surveys, but the bulk-flows predicted by our models are broadly in agreement with most available data.

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.

preprint1997arXiv

Projection effects in cluster catalogues

We investigate the importance of projection effects in the identification of galaxy clusters in 2D galaxy maps and their effect on the estimation of cluster velocity dispersions. A volume limited galaxy catalogue that was derived from a Standard CDM N-body simulation was used. We select clusters using criteria that match those employed in the construction of real cluster catalogues and find that our mock Abell cluster catalogues are heavily contaminated and incomplete. Over one third (34 per cent) of clusters of richness class R>=1 are miclassifications arising from the projection of one or more clumps onto an intrinsically poor cluster. Conversely, 32 per cent of intrinsically rich clusters are missed from the R>=1 catalogues, mostly because of statistical fluctuations in the background count. Selection by X-ray luminosity rather than optical richness reduces, but does not completely eliminate, these problems. Contamination by unvirialised sub-clumps near a cluster leads to a considerable overestimation of the cluster velocity dispersion, even if the analysis is restricted only to galaxies close to the cluster centre. Thus, the distribution of cluster masses - often used to test cosmological models - is a highly unreliable statistic. However, the median value of the distribution is considerably more robust. Improved estimates of the cluster velocity dispersion distribution require constructing new cluster catalogues in which clusters are defined according to the number of galaxies within a radius about three times smaller than the Abell radius.

preprint1997arXiv

The structure of galaxy clusters in different cosmologies

We investigate the internal structure of clusters of galaxies in high-resolution N-body simulations of 4 different cosmologies. There is a higher proportion of disordered clusters in critical-density than in low-density universes, although the structure of relaxed clusters is very similar in each. Crude measures of substructure, such as the shift in the position of the centre-of-mass as the density threshold is varied, can distinguish the two in a sample of just 20 or so clusters; it is harder to differentiate between clusters in open and flat models with the same density parameter. Most clusters are in a quasi-steady state within the virial radius and are well-described by the density profile of Navarro, Frenk & White (1995).

preprint1996arXiv

A wide-field K-band survey - I. Galaxy counts in B, V, I and K

We present bright galaxy number counts measured with linear detectors in the B, V, I, and K bands in two fields covering nearly 10 square degrees. All of our measurements are consistent with passive evolution models, and do not confirm the steep slope measured in other surveys at bright magnitudes. Throughout the range 16 < B < 19, our B-band counts are consistent with the "high normalization" models proposed to reduce the faint blue galaxy problem. Our K-band counts agree with previous measurements, and have reached a fair sample of the universe in the magnitude range where evolution and K-corrections are well understood.

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.

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.

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.

preprint1996arXiv

The Virgo consortium: simulations of dark matter and galaxy clustering

We report on work in progress by the Virgo consortium, a collaboration set up to carry out large simulations of the formation of galaxies and large-scale structure exploiting the latest generation of parallel supercomputers. We show results of $256^3$ particle N-body simulations of the clustering evolution of dark matter in four cold dark matter models with different cosmological parameters. The high resolution and large volume of these simulations allows us to determine reliably the mass autocorrelation function for pair separations in the range $40\hkpc$ to $20\hmpc$. Comparison of these with the observed galaxy correlation function shows that for any of these models to be viable, the distribution of galaxies must be biased relative to the distribution of mass in a non-trivial, scale-dependent fashion. In particular, low $Ω_0$ models require the galaxies to be more ``weakly'' clustered than the mass at small and intermediate pair separations. Simulations which include the evolution of gas show that cold gas knots form with approximately the abundance expected on theoretical grounds, although a few excessively massive objects grow near the centres of rich clusters. The locations where these cold gas knots form are, in general, biased relative to the distribution of mass in a scale-dependent way. Some of these biases have the required sign but they are, for the most part, weaker than is necessary for agreement with observations. The antibias present in our low $Ω_0$ N-body/SPH simulation appears to be related to the merging and disruption of galaxies in rich clusters.

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.

C. S. Frenk

What is connected

Connect this record

See the researcher in context

Building this map preview

81 published item(s)

DESI Data Release 1: Stellar Catalogue

Simulated Milky Way analogues: implications for dark matter indirect searches

Evolution of galaxy stellar masses and star formation rates in the EAGLE simulations

Galaxy And Mass Assembly (GAMA): end of survey report and data release 2

The impact of angular momentum on black hole accretion rates in simulations of galaxy formation

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

Galaxy And Mass Assembly (GAMA): Spectroscopic analysis

Scaling relations for galaxy clusters in the Millennium-XXL simulation

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

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

Galaxy and Mass Assembly (GAMA): Colour and luminosity dependent clustering from calibrated photometric redshifts

Galaxy And Mass Assembly (GAMA): The 0.013 < z < 0.1 cosmic spectral energy distribution from 0.1 micron to 1mm

Galaxy And Mass Assembly (GAMA): The mass-metallicity relationship

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

Principles of supernova-driven winds

The Aquila comparison Project: The Effects of Feedback and Numerical Methods on Simulations of Galaxy Formation

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

The Phoenix Project: the Dark Side of Rich Galaxy Clusters

The UV colours and dust attenuation of Lyman-break galaxies

A Two-Point Correlation Function For Galactic Halo Stars

Galaxy and Mass Assembly (GAMA): The GAMA Galaxy Group Catalogue (G3Cv1)

Galaxy And Mass Assembly (GAMA): the red fraction and radial distribution of satellite galaxies

Galaxy And Mass Assembly (GAMA): The star formation rate dependence of the stellar initial mass function

Galaxy and Mass Assembly (GAMA): ugriz galaxy luminosity functions

Mismatch and Misalignment: Dark Haloes and Satellites of Disc Galaxies

Substructure in the stellar halos of the Aquarius simulations

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

Where will supersymmetric dark matter first be seen?

Assembly History and Structure of Galactic Cold Dark Matter Halos

Galactic stellar haloes in the CDM model

Galaxy And Mass Assembly (GAMA): Galaxies at the faint end of the Halpha luminosity function

Galaxy and Mass Assembly (GAMA): survey diagnostics and core data release

Galaxy And Mass Assembly (GAMA): The input catalogue and star-galaxy separation

Grand unification of AGN activity in the LambdaCDM cosmology

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

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

Substructure Lensing: Effects of Galaxies, Globular Clusters & Satellite Streams

The earliest stars and their relics in the Milky Way

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

The evolution of Lyman-break galaxies in CDM

The Parameter Space of Galaxy Formation

The role of submillimetre galaxies in hierarchical galaxy formation

The spatial distribution of cold gas in hierarchical galaxy formation models

The statistics of the subhalo abundance of dark matter haloes

Effects of Dark Matter Substructures on Gravitational Lensing: Results from the Aquarius Simulations

Galaxy and Mass Assembly (GAMA): Optimal Tiling of Dense Surveys with a Multi-Object Spectrograph

Galaxy morphology in the LambdaCDM cosmology

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

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

The role of black holes in galaxy formation and evolution

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

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

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

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

The 2dF Galaxy Redshift Survey: Hierarchical galaxy clustering

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

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

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

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

Where are the stars?

Stable clustering, the halo model and nonlinear cosmological power spectra

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

Radio sources in the 2dF Galaxy Redshift Survey - II. Local radio luminosity functions for AGN and star-forming galaxies at 1.4 GHz

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

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

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

The 2dF Galaxy Redshift Survey: the luminosity function of cluster galaxies

The 2dF Galaxy Redshift Survey: The population of nearby radio galaxies at the 1 mJy level

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

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

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

Evolution of structure in cold dark matter universes

The seeds of rich galaxy clusters in the Universe

Projection effects in cluster catalogues