BZPEERPapers, people, evidence and research paths
Menu

Discover

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

Research tools

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

Network

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

Opportunities

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

Account

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

Source author record

Cedric G. Lacey

Cedric G. Lacey appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.COastro-ph.GAastro-ph.IMastro-phgr-qchep-ph

Catalog footprint

What is connected

36works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

36 published item(s)

preprint2022arXiv

Spin-driven jet feedback in idealised simulations of galaxy groups and clusters

We implement a black hole spin evolution and jet feedback model into SWIFT, a smoothed particle hydrodynamics code. The jet power is determined self-consistently assuming Bondi accretion, using a realistic, spin-dependant efficiency. The jets are launched along the spin axis of the black hole, resulting in natural reorientation and precession. We apply the model to idealised simulations of galaxy groups and clusters, finding that jet feedback successfully quenches gas cooling and star formation in all systems. Our group-size halo ($M_\mathrm{200}=10^{13}$ $\mathrm{M}_\odot$) is quenched by a strong jet episode triggered by a cooling flow, and it is kept quenched by a low-power jet fed from hot halo accretion. In more massive systems ($M_\mathrm{200}\geq 10^{14}$ $\mathrm{M}_\odot$), hot halo accretion is insufficient to quench the galaxies, or to keep them quenched after the first cooling episode. These galaxies experience multiple episodes of gas cooling, star formation and jet feedback. In the most massive galaxy cluster that we simulate ($M_\mathrm{200}=10^{15}$ $\mathrm{M}_\odot$), we find peak cold gas masses of $10^{10}$ $\mathrm{M}_\odot$ and peak star formation rates of a few times $100$ $\mathrm{M}_\odot\mathrm{yr}^{-1}$. These values are achieved during strong cooling flows, which also trigger the strongest jets with peak powers of $10^{47}$ $\mathrm{erg}\hspace{0.3mm}\mathrm{s}^{-1}$. These jets subsequently shut off the cooling flows and any associated star formation. Jet-inflated bubbles draw out low-entropy gas that subsequently forms dense cooling filaments in their wakes, as seen in observations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

AGNs at the cosmic dawn: predictions for future surveys from a $Λ$CDM cosmological model

Telescopes to be launched over the next decade-and-a-half, such as JWST, EUCLID, ATHENA and Lynx, promise to revolutionise the study of the high redshift Universe and greatly advance our understanding of the early stages of galaxy formation. We use a model that follows the evolution of the masses and spins of supermassive black holes (SMBHs) within a semi-analytic model of galaxy formation to make predictions for the Active Galactic Nucleus (AGN) luminosity function at $z\geq7$ in the broadband filters of JWST and EUCLID at near-infrared wavelengths, and ATHENA and Lynx at X-ray energies. The predictions of our model are relatively insensitive to the choice of seed black hole mass, except at the lowest luminosities ($L_{\mathrm{bol}}<10^{43}\mathrm{ergs^{-1}}$) and the highest redshifts ($z>10$). We predict that surveys with these different telescopes will select somewhat different samples of SMBHs, with EUCLID unveiling the most massive, highest accretion rate SMBHs, Lynx the least massive, lowest accretion rate SMBHs, and JWST and ATHENA covering objects inbetween. At $z=7$, we predict that typical detectable SMBHs will have masses, $M_{\mathrm{BH}}\sim10^{5-8}M_{\odot}$, and Eddington normalised mass accretion rates, $\dot{M}/\dot{M}_{\mathrm{Edd}}\sim0.6-2$. The SMBHs will be hosted by galaxies of stellar mass $M_{\star}\sim10^{8-10}M_{\odot}$, and dark matter haloes of mass $M_{\mathrm{halo}}\sim10^{11-12}M_{\odot}$. We predict that the detectable SMBHs at $z=10$ will have slightly smaller black holes, accreting at slightly higher Eddington normalised mass accretion rates, in slightly lower mass host galaxies compared to those at $z=7$, and reside in haloes of mass $M_{\mathrm{halo}}\sim10^{10-11}M_{\odot}$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Lyman-alpha emitters in a cosmological volume II: the impact of the intergalactic medium

In the near future galaxy surveys will target Lyman alpha emitting galaxies (LAEs) to unveil the nature of the dark energy. It has been suggested that the observability of LAEs is coupled to the large scale properties of the intergalactic medium. Such coupling could introduce distortions into the observed clustering of LAEs, adding a new potential difficulty to the interpretation of upcoming surveys. We present a model of LAEs that incorporates Lyman-alpha radiative transfer processes in the interstellar and intergalactic medium. The model is implemented in the GALFORM semi-analytic model of galaxy of formation and evolution. We find that the radiative transfer inside galaxies produces selection effects over galaxy properties. In particular, observed LAEs tend to have low metallicities and intermediate star formation rates. At low redshift we find no evidence of a correlation between the spatial distribution of LAEs and the intergalactic medium properties. However, at high redshift the LAEs are linked to the line of sight velocity and density gradient of the intergalactic medium. The strength of the coupling depends on the outflow properties of the galaxies and redshift. This effect modifies the clustering of LAEs on large scales, adding non linear features. In particular, our model predicts modifications to the shape and position of the baryon acoustic oscillation peak. This work highlights the importance of including radiative transfer physics in the cosmological analysis of LAEs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The evolution of radio jets across cosmic time

We present predictions for the evolution of radio emission from Active Galactic Nuclei (AGNs). We use a model that follows the evolution of Supermassive Black Hole (SMBH) masses and spins, within the latest version of the GALFORM semi-analytic model of galaxy formation. We use a Blandford-Znajek type model to calculate the power of the relativistic jets produced by black hole accretion discs, and a scaling model to calculate radio luminosities. First, we present the predicted evolution of the jet power distribution, finding that this is dominated by objects fuelled by hot halo accretion and an ADAF accretion state for jet powers above $10^{32}\mathrm{W}$ at $z=0$, with the contribution from objects fuelled by starbursts and in a thin disc accretion state being more important for lower jet powers at $z=0$ and at all jet powers at high redshifts ($z\geq3$). We then present the evolution of the jet power density from the model. The model is consistent with current observational estimates of jet powers from radio luminosities, once we allow for the significant uncertainties in these observational estimates. Next, we calibrate the model for radio emission to a range of observational estimates of the $z=0$ radio luminosity function. We compare the evolution of the model radio luminosity function to observational estimates for $0<z<6$, finding that the predicted evolution is similar to that observed. Finally, we explore recalibrating the model to reproduce luminosity functions of core radio emission, finding that the model is in approximate agreement with the observations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A unified multi-wavelength model of galaxy formation

We present a new version of the GALFORM semi-analytical model of galaxy formation. This brings together several previous developments of GALFORM into a single unified model, including a different initial mass function (IMF) in quiescent star formation and in starbursts, feedback from active galactic nuclei supressing gas cooling in massive halos, and a new empirical star formation law in galaxy disks based on their molecular gas content. In addition, we have updated the cosmology, introduced a more accurate treatment of dynamical friction acting on satellite galaxies, and updated the stellar population model. The new model is able to simultaneously explain both the observed evolution of the K-band luminosity function and stellar mass function, and the number counts and redshift distribution of sub-mm galaxies selected at 850 mu. This was not previously achieved by a single physical model within the LambdaCDM framework, but requires having an IMF in starbursts that is somewhat top-heavy. The new model is tested against a wide variety of observational data covering wavelengths from the far-UV to sub-mm, and redshifts from z=0 to z=6, and is found to be generally successful. These observations include the optical and near-IR luminosity functions, HI mass function, fraction of early type galaxies, Tully-Fisher, metallicity-luminosity and size-luminosity relations at z=0, as well as far-IR number counts, and far-UV luminosity functions at z ~ 3-6. [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Constraining SN feedback: a tug of war between reionization and the Milky Way satellites

Theoretical models of galaxy formation based on the cold dark matter cosmogony typically require strong feedback from supernova (SN) explosions in order to reproduce the Milky Way satellite galaxy luminosity function and the faint end of the field galaxy luminosity function. However, too strong a SN feedback also leads to the universe reionizing too late, and the metallicities of Milky Way satellites being too low. The combination of these four observations therefore places tight constraints on SN feedback. We investigate these constraints using the semi-analytical galaxy formation model galform. We find that these observations favour a SN feedback model in which the feedback strength evolves with redshift. We find that, for our best fit model, half of the ionizing photons are emitted by galaxies with rest-frame far-UV absolute magnitudes $M_{\rm AB}(1500{\rm Å})<-17.5$, which implies that already observed galaxy populations contribute about half of the photons responsible for reionization. The $z=0$ descendants of these galaxies are mainly galaxies with stellar mass $M_*>10^{10}\,{\rm M}_{\odot}$ and preferentially inhabit halos with mass $M_{\rm halo}>10^{13}\,{\rm M}_{\odot}$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Galaxies in the EAGLE hydrodynamical simulation and in the Durham and Munich semi-analytical models

We compare global predictions from the EAGLE hydrodynamical simulation, and two semi-analytic (SA) models of galaxy formation, L-GALAXIES and GALFORM. All three models include the key physical processes for the formation and evolution of galaxies and their parameters are calibrated against a small number of observables at $z\approx 0$. The two SA models have been applied to merger trees constructed from the EAGLE dark matter only simulation. We find that at $z\leq 2$, both the galaxy stellar mass functions for stellar masses $M_*<10^{10.5} M_{\odot}$ and the median specific star formation rates (sSFRs) in the three models agree to better than $0.4$~dex. The evolution of the sSFR predicted by the three models closely follows the mass assembly history of dark matter haloes. In both EAGLE and L-GALAXIES there are more central passive galaxies with $M_*<10^{9.5} M_{\odot}$ than in L-GALAXIES. This difference is related to galaxies that have entered and then left a larger halo and which are treated as satellites in GALFORM. In the range $0<z<1$, the slope of the evolution of the star formation rate density in EAGLE is a factor of $\approx 1.5$ steeper than for the two SA models. The median sizes for galaxies with $M_*>10^{9.5} M_{\odot}$ differ in some instances by an order of magnitude, while the stellar mass-size relation in EAGLE is a factor of $\approx 2$ tighter than for the two SA models. Our results suggest the need for a revision of how SA models treat the effect of baryonic self-gravity on the underlying dark matter. The treatment of gas flows in the models needs to be revised based on detailed comparison with observations to understand in particular the evolution of the stellar mass-metallicity relation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Reionisation in sterile neutrino cosmologies

We investigate the process of reionisation in a model in which the dark matter is a warm elementary particle such as a sterile neutrino. We focus on models that are consistent with the dark matter decay interpretation of the recently detected line at 3.5 keV in the X-ray spectra of galaxies and clusters. In warm dark matter models the primordial spectrum of density perturbations has a cut-off on the scale of dwarf galaxies. Structure formation therefore begins later than in the standard cold dark matter (CDM) model and very few objects form below the cut-off mass scale. To calculate the number of ionising photons, we use the Durham semi-analytic model of galaxy formation, GALFORM. We find that even the most extreme 7 keV sterile neutrino we consider is able to reionise the Universe early enough to be compatible with the bounds on the epoch of reionisation from Planck. This, perhaps surprising, result arises from the rapid build-up of high redshift galaxies in the sterile neutrino models which is also reflected in a faster evolution of their far-UV luminosity function between $10>z>7$ than in CDM. The dominant sources of ionising photons are systematically more massive in the sterile neutrino models than in CDM. As a consistency check on the models, we calculate the present-day luminosity function of satellites of Milky Way-like galaxies. When the satellites recently discovered in the DES survey are taken into account, strong constraints are placed on viable sterile neutrino models.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Clustering Evolution of Dusty Star-Forming Galaxies

We present predictions for the clustering of galaxies selected by their emission at far infra-red (FIR) and sub-millimetre wavelengths. This includes the first predictions for the effect of clustering biases induced by the coarse angular resolution of single-dish telescopes at these wavelengths. We combine a new version of the GALFORM model of galaxy formation with a self-consistent model for calculating the absorption and re-emission of radiation by interstellar dust. Model galaxies selected at $850$ $μ$m reside in dark matter halos of mass $M_{\rm halo}\sim10^{11.5}-10^{12}$ $h^{-1}$ M$_{\odot}$, independent of redshift (for $0.2\lesssim z\lesssim4$) or flux (for $0.25\lesssim S_{850μ\rm m}\lesssim4$ mJy). At $z\sim2.5$, the brightest galaxies ($S_{850μ\rm m}>4$ mJy) exhibit a correlation length of $r_{0}=5.5_{-0.5}^{+0.3}$ $h^{-1}$ Mpc, consistent with observations. We show that these galaxies have descendants with stellar masses $M_{\star}\sim10^{11}$ $h^{-1}$ M$_{\odot}$ occupying halos spanning a broad range in mass $M_{\rm halo}\sim10^{12}-10^{14}$ $h^{-1}$ M$_{\odot}$. The FIR emissivity at shorter wavelengths ($250$, $350$ and $500$ $μ$m) is also dominated by galaxies in the halo mass range $M_{\rm halo}\sim10^{11.5}-10^{12}$ $h^{-1}$ M$_{\odot}$, again independent of redshift (for $0.5\lesssim z\lesssim5$). We compare our predictions for the angular power spectrum of cosmic infra-red background anisotropies at these wavelengths with observations, finding agreement to within a factor of $\sim2$ over all scales and wavelengths, an improvement over earlier versions of the model. Simulating images at $850$ $μ$m, we show that confusion effects boost the measured angular correlation function on all scales by a factor of $\sim4$. This has important consequences, potentially leading to inferred halo masses being overestimated by an order of magnitude.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The clustering of dark matter halos: scale-dependent bias on quasi-linear scales

We investigate the spatial clustering of dark matter halos, collapsing from $1-4 σ$ fluctuations, in the redshift range $0 - 5$ using N-body simulations. The halo bias of high redshift halos ($z \geq 2$) is found to be strongly non-linear and scale-dependent on quasi-linear scales that are larger than their virial radii ($0.5-10$ Mpc/h). However, at lower redshifts, the scale-dependence of non-linear bias is weaker and and is of the order of a few percent on quasi-linear scales at $z \sim 0$. We find that the redshift evolution of the scale dependent bias of dark matter halos can be expressed as a function of four physical parameters: the peak height of halos, the non-linear matter correlation function at the scale of interest, an effective power law index of the {\it rms} linear density fluctuations and the matter density of the universe at the given redshift. This suggests that the scale-dependence of halo bias is not a universal function of the dark matter power spectrum, which is commonly assumed. We provide a fitting function for the scale dependent halo bias as a function of these four parameters. Our fit reproduces the simulation results to an accuracy of better than 4 % over the redshift range $0\leq z \leq 5$. We also extend our model by expressing the non-linear bias as a function of the linear matter correlation function. It is important to incorporate our results into the clustering models of dark matter halos at any redshift, including those hosting early generations of stars and galaxies before reionization.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The DESI Experiment Part I: Science,Targeting, and Survey Design

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$α$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The DESI Experiment Part II: Instrument Design

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= λ/Δλ$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The metal enrichment of passive galaxies in cosmological simulations of galaxy formation

Massive early-type galaxies have higher metallicities and higher ratios of $α$ elements to iron than their less massive counterparts. Reproducing these correlations has long been a problem for hierarchical galaxy formation theory, both in semi-analytic models and cosmological hydrodynamic simulations. We show that a simulation in which gas cooling in massive dark haloes is quenched by radio-mode active galactic nuclei (AGNs) feedback naturally reproduces the observed trend between $α$/Fe and the velocity dispersion of galaxies, $σ$. The quenching occurs earlier for more massive galaxies. Consequently, these galaxies complete their star formation before $α$/Fe is diluted by the contribution from type Ia supernovae. For galaxies more massive than $\sim 10^{11}~M_\odot$ whose $α$/Fe correlates positively with stellar mass, we find an inversely correlated mass-metallicity relation. This is a common problem in simulations in which star formation in massive galaxies is quenched either by quasar- or radio-mode AGN feedback. The early suppression of gas cooling in progenitors of massive galaxies prevents them from recapturing enriched gas ejected as winds. Simultaneously reproducing the [$α$/Fe]-$σ$ relation and the mass-metallicity relation is, thus, difficult in the current framework of galaxy formation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A new methodology to test galaxy formation models using the dependence of clustering on stellar mass

We present predictions for the two-point correlation function of galaxy clustering as a function of stellar mass, computed using two new versions of the GALFORM semi-analytic galaxy formation model. These models make use of a high resolution, large volume N-body simulation, set in the WMAP7 cosmology. One model uses a universal stellar initial mass function (IMF), while the other assumes different IMFs for quiescent star formation and bursts. Particular consideration is given to how the assumptions required to estimate the stellar masses of observed galaxies (such as the choice of IMF, stellar population synthesis model and dust extinction) influence the perceived dependence of galaxy clustering on stellar mass. Broad-band spectral energy distribution fitting is carried out to estimate stellar masses for the model galaxies in the same manner as in observational studies. We show clear differences between the clustering signals computed using the true and estimated model stellar masses. As such, we highlight the importance of applying our methodology to compare theoretical models to observations. We introduce an alternative scheme for the calculation of the merger timescales for satellite galaxies in GALFORM, which takes into account the dark matter subhalo information from the simulation. This reduces the amplitude of small-scale clustering. The new merger scheme offers improved or similar agreement with observational clustering measurements, over the redshift range 0 < z < 0.7. We find reasonable agreement with clustering measurements from GAMA, but find larger discrepancies for some stellar mass ranges and separation scales with respect to measurements from SDSS and VIPERS, depending on the GALFORM model used.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Galaxy and Mass Assembly (GAMA): Projected Galaxy Clustering

We measure the projected 2-point correlation function of galaxies in the 180 deg$^2$ equatorial regions of the GAMA II survey, for four different redshift slices between z = 0.0 and z=0.5. To do this we further develop the Cole (2011) method of producing suitable random catalogues for the calculation of correlation functions. We find that more r-band luminous, more massive and redder galaxies are more clustered. We also find that red galaxies have stronger clustering on scales less than ~3 $h^{-1}$ Mpc. We compare to two different versions of the GALFORM galaxy formation model, Lacey et al (in prep.) and Gonzalez-Perez et al. (2014), and find that the models reproduce the trend of stronger clustering for more massive galaxies. However, the models under predict the clustering of blue galaxies, can incorrectly predict the correlation function on small scales and under predict the clustering in our sample of galaxies with ~3$L_r$ . We suggest possible avenues to explore to improve these cluster- ing predictions. The measurements presented in this paper can be used to test other galaxy formation models, and we make the measurements available online to facilitate this.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The environments of high redshift radio galaxies and quasars: probes of protoclusters

We use the GALFORM semi-analytical model to study high density regions traced by radio galaxies and quasars at high redshifts. We explore the impact that baryonic physics has upon the properties of galaxies in these environments. Star-forming emission-line galaxies (Lyα and Hα emitters) are used to probe the environments at high redshifts. Radio galaxies are predicted to be hosted by more massive haloes than quasars, and this is imprinted on the amplitude of galaxy overdensities and cross-correlation functions. We find that Lyα radiative transfer and AGN feedback indirectly affect the clustering on small scales and also the stellar masses, star- formation rates and gas metallicities of galaxies in dense environments. We also investigate the relation between protoclusters associated with radio galaxies and quasars, and their present- day cluster descendants. The progenitors of massive clusters associated with radio galaxies and quasars allow us to determine an average protocluster size in a simple way. Overdensities within the protoclusters are found to correlate with the halo descendant masses. We present scaling relations that can be applied to observational data. By computing projection effects due to the wavelength resolution of modern spectrographs and narrow-band filters we show that the former have enough spectral resolution to map the structure of protoclusters, whereas the latter can be used to measure the clustering around radio galaxies and quasars over larger scales to determine the mass of dark matter haloes hosting them.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Simulated observations of sub-millimetre galaxies: the impact of single-dish resolution and field variance

Recent observational evidence suggests that the coarse angular resolution ($\sim20"$ FWHM) of single-dish telescopes at sub-mm wavelengths has biased the observed galaxy number counts by blending together the sub-mm emission from multiple sub-mm galaxies (SMGs). We use lightcones computed from an updated implementation of the \galform semi-analytic model to generate $50$ mock sub-mm surveys of $0.5$ deg$^2$ at $850$ $μ$m, taking into account the effects of the finite single-dish beam in a more accurate way than has been done previously. We find that blending of SMGs does lead to an enhancement of source extracted number counts at bright fluxes ($S_{\mathrm{850}μ\mathrm{m}}\gtrsim1$ mJy). Typically, $\sim3{-}6$ galaxies contribute $90\%$ of the flux of an $S_{850μ\mathrm{m}}=5$ mJy source and these blended galaxies are physically unassociated. We find that field-to-field variations are comparable to Poisson fluctuations for our $S_{850μ\mathrm{m}}>5$ mJy SMG population, which has a median redshift $z_{50}=2.0$, but are greater than Poisson for the $S_{850μ\mathrm{m}}>1$ mJy population ($z_{50}=2.8$). In a detailed comparison to a recent interferometric survey targeted at single-dish detected sources, we reproduce the difference between single-dish and interferometer number counts and find a median redshift ($z_{50}=2.5$) in excellent agreement with the observed value ($z_{50}=2.5\pm 0.2$). We also present predictions for single-dish survey number counts at $450$ and $1100$ $μ$m, which show good agreement with observational data.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The environments of Ly$α$ blobs I: Wide-field Ly$α$ imaging of TN J1338-1942, a powerful radio galaxy at $z\simeq 4.1$ associated with a giant Ly$α$ nebula

We exploit wide-field Ly$alpha$ imaging with Subaru to probe the environment around TN J1338-1942, a powerful radio galaxy with a >100 kpc Ly$α$ halo at z=4.11. We used a sample of Ly$α$ emitters (LAEs) down to $\log(L_{\rm Lyα} [erg\, s^{-1}])\sim 42.8$ to measure the galaxy density around TNJ1338, compared to a control sample from a blank field taken with the same instrument. We found that TNJ1338 resides in a region with a peak overdensity of $δ_{\rm LAE}=2.8\pm 0.5$ on scales of $8\, h^{-1}\rm Mpc$ (on the sky) and $112\, h^{-1}\rm Mpc$ (line of sight) in comoving coordinates. Adjacent to this overdensity, we found a strong underdensity where virtually no LAEs are detected. We used a semi-analytical model of LAEs derived from the Millennium Simulation to compare our results with theoretical predictions. While the theoretical density distribution is consistent with the blank field, overdense regions such as that around TNJ1338 are very rare, with a number density of $6.4\times 10^{-8}\rm Mpc^{-3}$ (comoving), corresponding to the densest < 0.4 percentile at $z\simeq 4.1$. We also found that the Ly$α$ luminosity function in the TNJ1338 field differs from that in the blank field: the number of bright LAEs ($\log(L_{\rm Lyα}[erg\,s^{-1}]) \gtrsim 43.3$) is enhanced, while the number of fainter LAEs is relatively suppressed. These results suggest that some powerful radio galaxies associated with Ly$α$ nebulae reside in extreme overdensities on $\sim 3$--$6\, \rm Mpc$ scales, where star-formation and AGN activity may be enhanced via frequent galaxy mergers or high rates of gas accretion from the surroundings.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The evolution of the star forming sequence in hierarchical galaxy formation models

It has been argued that the specific star formation rates of star forming galaxies inferred from observational data decline more rapidly below z = 2 than is predicted by hierarchical galaxy formation models. We present a detailed analysis of this problem by comparing predictions from the GALFORM semi-analytic model with an extensive compilation of data on the average star formation rates of star-forming galaxies. We also use this data to infer the form of the stellar mass assembly histories of star forming galaxies. Our analysis reveals that the currently available data favour a scenario where the stellar mass assembly histories of star forming galaxies rise at early times and then fall towards the present day. In contrast, our model predicts stellar mass assembly histories that are almost flat below z = 2 for star forming galaxies, such that the predicted star formation rates can be offset with respect to the observational data by factors of up to 2-3. This disagreement can be explained by the level of coevolution between stellar and halo mass assembly that exists in contemporary galaxy formation models. In turn, this arises because the standard implementations of star formation and supernova feedback used in the models result in the efficiencies of these process remaining approximately constant over the lifetime of a given star forming galaxy. We demonstrate how a modification to the timescale for gas ejected by feedback to be reincorporated into galaxy haloes can help to reconcile the model predictions with the data.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A dynamical model of supernova feedback: gas outflows from the interstellar medium

We present a dynamical model of supernova feedback which follows the evolution of pressurised bubbles driven by supernovae in a multi-phase interstellar medium (ISM). The bubbles are followed until the point of break-out into the halo, starting from an initial adiabatic phase to a radiative phase. We show that a key property which sets the fate of bubbles in the ISM is the gas surface density, through the work done by the expansion of bubbles and its role in setting the gas scaleheight. The multi-phase description of the ISM is essential, and neglecting it leads to order of magnitude differences in the predicted outflow rates. We compare our predicted mass loading and outflow velocities to observations of local and high-redshift galaxies and find good agreement over a wide range of stellar masses and velocities. With the aim of analysing the dependence of the mass loading of the outflow, beta (i.e. the ratio between the outflow and star formation rates), on galaxy properties, we embed our model in the galaxy formation simulation, GALFORM, set in the LCDM framework. We find that a dependence of beta solely on the circular velocity, as is widely assumed in the literature, is actually a poor description of the outflow rate, as large variations with redshift and galaxy properties are obtained. Moreover, we find that below a circular velocity of 80km/s the mass loading saturates. A more fundamental relation is that between beta and the gas scaleheight of the disk, hg, and the gas fraction, fgas, as beta hg^(1.1) fgas^(0.4), or the gas surface density, Σ_g, and the gas fraction, as beta Σ_g^(-0.6) fgas^(0.8). We find that using the new mass loading model leads to a shallower faint-end slope in the predicted optical and near-IR galaxy luminosity functions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Clustering of Extremely Red Objects in Elais-N1 from the UKIDSS DXS with optical photometry from Pan-STARRS1 and Subaru

We measure the angular clustering of 33 415 extremely red objects (EROs) in the Elais-N1 field covering 5.33 deg$^{2}$, which cover the redshift range $z=0.8$ to $2$. This sample was made by merging the UKIDSS Deep eXtragalactic Survey (DXS) with the optical Subaru and Pan-STARRS PS1 datasets. We confirm the existence of a clear break in the angular correlation function at $\sim 0.02^{\circ}$ corresponding to $1 h^{-1}$ Mpc at $z\sim1$. We find that redder or brighter EROs are more clustered than bluer or fainter ones. Halo Occupation Distribution (HOD) model fits imply that the average mass of dark matter haloes which host EROs is over $10^{13} h^{-1} M_{\odot}$ and that EROs have a bias ranging from 2.7 to 3.5. Compared to EROs at $z\sim1.1$, at $z\sim1.5$ EROs have a higher bias and fewer are expected to be satellite galaxies. Furthermore, EROs reside in similar dark matter haloes to those that host $10^{11.0} M_{\odot}<M_{*}<10^{11.5} M_{\odot}$ galaxies. We compare our new measurement and HOD fits with the predictions of the GALFORM semi-analytical galaxy formation model. Overall, the clustering predicted by GALFORM gives an encouraging match to our results. However, compared to our deductions from the measurements, GALFORM puts EROs into lower mass haloes and predicts that a larger fraction of EROs are satellite galaxies. This suggests that the treatment of gas cooling may need to be revised in the model. Our analysis illustrates the potential of clustering analyses to provide observational constraints on theoretical models of galaxy formation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Extending the halo mass resolution of $N$-body simulations

We present a scheme to extend the halo mass resolution of N-body simulations of the hierarchical clustering of dark matter. The method uses the density field of the simulation to predict the number of sub-resolution dark matter haloes expected in different regions. The technique requires as input the abundance of haloes of a given mass and their average clustering, as expressed through the linear and higher order bias factors. These quantities can be computed analytically or, more accurately, derived from a higher resolution simulation as done here. Our method can recover the abundance and clustering in real- and redshift-space of haloes with mass below $\sim 7.5 \times 10^{13}h^{-1}M_{\odot}$ at $z=0$ to better than 10%. We demonstrate the technique by applying it to an ensemble of 50 low resolution, large-volume $N$-body simulations to compute the correlation function and covariance matrix of luminous red galaxies (LRGs). The limited resolution of the original simulations results in them resolving just two thirds of the LRG population. We extend the resolution of the simulations by a factor of 30 in halo mass in order to recover all LRGs. With existing simulations it is possible to generate a halo catalogue equivalent to that which would be obtained from a $N$-body simulation using more than 20 trillion particles; a direct simulation of this size is likely to remain unachievable for many years. Using our method it is now feasible to build the large numbers of high-resolution large volume mock galaxy catalogues required to compute the covariance matrices necessary to analyse upcoming galaxy surveys designed to probe dark energy.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Lightcone mock catalogues from semi-analytic models of galaxy formation - I. Construction and application to the BzK colour selection

We introduce a method for constructing end-to-end mock galaxy catalogues using a semi-analytical model of galaxy formation, applied to the halo merger trees extracted from a cosmological N-body simulation. The mocks that we construct are lightcone catalogues, in which a galaxy is placed according to the epoch at which it first enters the past lightcone of the observer, and incorporate the evolution of galaxy properties with cosmic time. We determine the position between the snapshot outputs at which a galaxy enters the observer's lightcone by interpolation. As an application, we consider the effectiveness of the BzK colour selection technique, which was designed to isolate galaxies in the redshift interval 1.4<z<2.5. The mock catalogue is in reasonable agreement with the observed number counts of all BzK galaxies, as well as with the observed counts of the subsample of BzKs that are star-forming galaxies. We predict that over 75 per cent of the model galaxies with K_{AB}<=23, and 1.4<z<2.5, are selected by the BzK technique. Interloper galaxies, outside the intended redshift range, are predicted to dominate bright samples of BzK galaxies (i.e. with K_{AB}<=21). Fainter K-band cuts are necessary to reduce the predicted interloper fraction. We also show that shallow B-band photometry can lead to confusion in classifying BzK galaxies as being star-forming or passively evolving. Overall, we conclude that the BzK colour selection technique is capable of providing a sample of galaxies that is representative of the 1.4<z<2.5 galaxy population.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Single Colour Diagnostics of the Mass-to-light Ratio: Predictions from Galaxy Formation Models

Accurate galaxy stellar masses are crucial to better understand the physical mechanisms driving the galaxy formation process. We use synthetic star formation and metal enrichment histories predicted by the {\sc galform} galaxy formation model to investigate the precision with which various colours $(m_{a}-m_{b})$ can alone be used as diagnostics of the stellar mass-to-light ratio. As an example, we find that, at $z=0$, the {\em intrinsic} (B$_{f435w}-$V$_{f606w}$) colour can be used to determine the intrinsic rest-frame $V$-band stellar mass-to-light ratio ($\log_{10}Γ_{V}=\log_{10}[(M/M_{\odot})/(L_{V}/L_{V\odot})]$) with a precision of $σ_{lgΓ}\simeq 0.06$ when the initial mass function and redshift are known beforehand. While the presence of dust, assuming a universal attenuation curve, can have a systematic effect on the inferred mass-to-light ratio using a single-colour relation, this is typically small as it is often possible to choose a colour for which the dust reddening vector is approximately aligned with the $(m_{a}-m_{b})-\log_{10}Γ_{V}$ relation. The precision with which the stellar mass-to-light ratio can be recovered using a single colour diagnostic rivals implementations of SED fitting using more information but in which simple parameterisations of the star formation and metal enrichment histories are assumed. To facilitate the wide use of these relations, we provide the optimal observer frame colour to estimate the stellar mass-to-light ratio, along with the associated parameters, as a function of redshift ($0<z<1.5$) for two sets of commonly used filters.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Predictions for the CO emission of galaxies from a coupled simulation of galaxy formation and photon dominated regions

We combine the galaxy formation model GALFORM with the Photon Dominated Region code UCL_PDR to study the emission from the rotational transitions of 12CO (CO) in galaxies from z=0 to z=6 in the Lambda CDM framework. GALFORM is used to predict the molecular (H2) and atomic hydrogen (HI) gas contents of galaxies using the pressure-based empirical star formation relation of Blitz & Rosolowsky. From the predicted H2 mass and the conditions in the interstellar medium, we estimate the CO emission in the rotational transitions 1-0 to 10-9 by applying the UCL_PDR model to each galaxy. We find that deviations from the Milky-Way CO-H2 conversion factor come mainly from variations in metallicity, and in the average gas and star formation rate surface densities. In the local universe, the model predicts a CO(1-0) luminosity function (LF), CO-to-total infrared (IR) luminosity ratios for multiple CO lines and a CO spectral line energy distribution (SLED) which are in good agreement with observations of luminous and ultra-luminous IR galaxies. At high redshifts, the predicted CO SLED of the brightest IR galaxies reproduces the shape and normalization of the observed CO SLED. The model predicts little evolution in the CO-to-IR luminosity ratio for different CO transitions, in good agreement with observations up to z~5. We use this new hybrid model to explore the potential of using colour selected samples of high-redshift star-forming galaxies to characterise the evolution of the cold gas mass in galaxies through observations with the Atacama Large Millimeter Array.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Predictions for the intrinsic UV continuum properties of star forming galaxies and the implications for inferring dust extinction

The observed ultraviolet continuum (UVC) slope is potentially a powerful diagnostic of dust obscuration in star forming galaxies. However, the intrinsic slope is also sensitive to the form of the stellar initial mass function (IMF) and to the recent star formation and metal enrichment histories of a galaxy. Using the galform semi-analytical model of galaxy formation, we investigate the intrinsic distribution of UVC slopes. For star-forming galaxies, we find that the intrinsic distribution of UVC slopes at z=0, parameterised by the power law index beta, has a standard deviation of sigma_beta=0.30. This suggests an uncertainty on the inferred UV attenuation of A_fuv=0.7$ (assuming a Calzetti attenuation curve) for an individual object, even with perfect photometry. Furthermore, we find that the intrinsic UVC slope correlates with star formation rate, intrinsic UV luminosity, stellar mass and redshift. These correlations have implications for the interpretation of trends in the observed UVC slope with these quantities irrespective of the sample size or quality of the photometry. Our results suggest that in some cases the attenuation by dust has been incorrectly estimated.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Simulations and modelling of the ISM in galaxies

The latest observations of molecular gas and the atomic hydrogen content of local and high-redshift galaxies, coupled with how these correlate with star formation activity, have revolutionized our ideas about how to model star formation in a galactic context. A successful theory of galaxy formation has to explain some key facts: (i) high-redshift galaxies have higher molecular gas fractions and star formation rates than local galaxies, (ii) scaling relations show that the atomic-to-stellar mass ratio decreases with stellar mass in the local Universe, and (iii) the global abundance of atomic hydrogen evolves very weakly with time. We review how modern cosmological simulations of galaxy formation attempt to put these pieces together and highlight how approaches simultaneously solving dark matter and gas physics, and approaches first solving the dark matter N-body problem and then dealing with gas physics using semi-analytic models, differ and complement each other. We review the observable predictions, what we think we have learned so far and what still needs to be done in the simulations to allow robust testing by the new observations expected from telescopes such as ALMA, PdBI, LMT, JVLA, ASKAP, MeerKAT, SKA.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The accuracy of the UV continuum as an indicator of the star formation rate in galaxies

The rest-frame intrinsic UV luminosity is often used as an indicator of the instantaneous star formation rate (SFR) in a galaxy. While it is in general a robust indicator of the ongoing star formation activity, the precise value of the calibration relating the UV luminosity to the SFR ($B_ν$), is sensitive to various physical properties, such as the recent star formation and metal enrichment histories, along with the choice of stellar initial mass function. The distribution of these properties for the star-forming galaxy population then suggests that the adoption of a single calibration is not appropriate unless properly qualified with the uncertainties on the calibration. We investigate, with the aid of the {\sc galform} semi-analytic model of galaxy formation, the distribution of UV-SFR calibrations obtained using realistic star formation and metal enrichment histories. At $z=0$, we find that when the initial mass function is fixed (to the Kennicutt IMF), the median calibration is $B_{\rm fuv}=0.9$ where ${\rm SFR}/[{\rm M_{\odot}\,yr^{-1}}]=B_ν\times 10^{-28}\times L_ν/[{\rm ergs\,s^{-1}\,Hz^{-1}}]$. However, the width of the distribution $B_{\rm fuv}$ suggests that for a single object there is around a 20% {\em intrinsic} uncertainty (at $z=0$, rising to $\simeq 30%$ at $z=6$) on the star formation rate inferred from the FUV luminosity without additional constraints on the star formation history or metallicity. We also find that the median value of the calibration $B_{\rm fuv}$ is correlated with the star formation rate and redshift (at $z>3$) raising implications for the correct determination of the star formation rate from the UV.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Can galactic outflows explain the properties of Ly-alpha emitters?

We study the properties of Ly-alpha emitters in a cosmological framework by computing the escape of Ly-alpha photons through galactic outflows. We combine the GALFORM semi-analytical model of galaxy formation with a Monte Carlo Ly-alpha radiative transfer code. The properties of Ly-alpha emitters at 0<z<7 are predicted using two outflow geometries: a Shell of neutral gas and a Wind ejecting material, both expanding at constant velocity. We characterise the differences in the Ly-alpha line profiles predicted by the two outflow geometries in terms of their width, asymmetry and shift from the line centre for a set of outflows with different hydrogen column densities, expansion velocities and metallicities. In general, the Ly-alpha line profile of the Shell geometry is broader and more asymmetric, and the Ly-alpha escape fraction is lower than with the Wind geometry for the same set of parameters. In order to implement the outflow geometries in the semi-analytical model GALFORM, a number of free parameters in the outflow model are set by matching the luminosity function of Ly-alpha emitters over the whole observed redshift range. The models are consistent with the observationally inferred Ly-alpha escape fractions, equivalent width distributions and with the shape of the Ly-alpha line from composite spectra. Interestingly, our predicted UV luminosity function of Ly-alpha emitters and the fraction of Ly-alpha emitters in Lyman-break galaxy samples at high redshift are in partial agreement with observations. Attenuation of the Ly-alpha line by the presence of a neutral intergalactic medium at high redshift could be responsible for this disagreement. We predict that Ly-alpha emitters constitute a subset of the galaxy population with lower metallicities, lower instantaneous star formation rates and larger sizes than the overall population at the same UV luminosity.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Cosmic evolution of the atomic and molecular gas content of galaxies

We study the evolution of the cold gas content of galaxies by splitting the interstellar medium into its atomic and molecular hydrogen components, using the galaxy formation model GALFORM in the LCDM framework. We calculate the molecular-to-atomic hydrogen mass ratio, H2/HI, in each galaxy using two different approaches; the pressure-based empirical relation of Blitz & Rosolowsky and the theoretical model of Krumholz, McKeee & Tumlinson, and apply them to consistently calculate the star formation rates of galaxies. We find that the model based on the Blitz & Rosolowsky law predicts an HI mass function, CO(1-0) luminosity function, correlations between the H2/HI ratio and stellar and cold gas mass, and infrared-CO luminosity relation in good agreement with local and high redshift observations. The HI mass function evolves weakly with redshift, with the number density of high mass galaxies decreasing with increasing redshift. In the case of the H2 mass function, the number density of massive galaxies increases strongly from z=0 to z=2, followed by weak evolution up to z=4. We also find that the H2/HI ratio of galaxies is strongly dependent on stellar and cold gas mass, and also on redshift. The slopes of the correlations between H2/HI and stellar and cold gas mass hardly evolve, but the normalisation increases by up to two orders of magnitude from z=0-8. The strong evolution in the H2 mass function and the H2/HI ratio is primarily due to the evolution in the sizes of galaxies and secondarily, in the gas fractions. The predicted cosmic density evolution of HI agrees with the observed evolution inferred from DLAs, and is dominated by low/intermediate mass halos. We find that previous theoretical studies have largely overestimated the redshift evolution of the global H2/HI ratio due to limited resolution. We predict a maximum of rho_H2/rho_HI~1.2 at z~3.5.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Numerical overcooling in shocks

We present a study of cooling in radiative shocks simulated with smoothed particle hydrodynamics (SPH) and adaptive mesh refinement codes. We obtain a similarity solution for a shock-tube problem in the presence of radiative cooling, and test how well the solution is reproduced in Gadget and Flash. Shock broadening governed by the details of the numerical scheme (artificial viscosity or Riemann solvers) leads to potentially significant overcooling in both codes. We interpret our findings in terms of a resolution criterion, and apply it to realistic simulations of cosmological accretion shocks onto galaxy haloes, cold accretion and thermal feedback from supernovae or active galactic nuclei. To avoid numerical overcooling of accretion shocks onto haloes that should develop a hot corona requires a particle or cell mass resolution of 10^6 M_sun, which is within reach of current state-of-the-art simulations. At this mass resolution, thermal feedback in the interstellar medium of a galaxy requires temperatures of supernova or AGN driven bubbles to be in excess of 10^7 K at densities of n_H=1.0 cm^-3, in order to avoid spurious suppression of the feedback by numerical overcooling.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the impact of empirical and theoretical star formation laws on galaxy formation

We investigate the consequences of applying different star formation laws in the galaxy formation model GALFORM. Three broad star formation laws are implemented: the empirical relations of Kennicutt and Schmidt and Blitz & Rosolowsky and the theoretical model of Krumholz, McKee & Tumlinson. These laws have no free parameters once calibrated against observations of the star formation rate (SFR) and gas surface density in nearby galaxies. We start from published models, and investigate which observables are sensitive to a change in the star formation law, without altering any other model parameters. We show that changing the star formation law (i) does not significantly affect either the star formation history of the universe or the galaxy luminosity functions in the optical and near-IR, due to an effective balance between the quiescent and burst star formation modes; (ii) greatly affects the cold gas contents of galaxies; (iii) changes the location of galaxies in the SFR versus stellar mass plane, so that a second sequence of "passive" galaxies arises, in addition to the known "active" sequence. We show that this plane can be used to discriminate between the star formation laws.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Designing a space-based galaxy redshift survey to probe dark energy

A space-based galaxy redshift survey would have enormous power in constraining dark energy and testing general relativity, provided that its parameters are suitably optimized. We study viable space-based galaxy redshift surveys, exploring the dependence of the Dark Energy Task Force (DETF) figure-of-merit (FoM) on redshift accuracy, redshift range, survey area, target selection, and forecast method. Fitting formulae are provided for convenience. We also consider the dependence on the information used: the full galaxy power spectrum P(k, P(k) marginalized over its shape, or just the Baryon Acoustic Oscillations (BAO). We find that the inclusion of growth rate information (extracted using redshift space distortion and galaxy clustering amplitude measurements) leads to a factor of ~ 3 improvement in the FoM, assuming general relativity is not modified. This inclusion partially compensates for the loss of information when only the BAO are used to give geometrical constraints, rather than using the full P(k) as a standard ruler. We find that a space-based galaxy redshift survey covering ~20,000 deg^2 over 0.5< z < 2 with σ_z/(1+z)<= 0.001 exploits a redshift range that is only easily accessible from space, extends to sufficiently low redshifts to allow both a vast 3-D map of the universe using a single tracer population, and overlaps with ground-based surveys to enable robust modeling of systematic effects. We argue that these parameters are close to their optimal values given current instrumental and practical constraints.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

The Colours of Satellite Galaxies in Groups and Clusters

Current models of galaxy formation predict satellite galaxies in groups and clusters that are redder than observed. We investigate the effect on the colours of satellite galaxies produced by the ram pressure stripping of their hot gaseous atmospheres as the satellites orbit within their parent halo. We incorporate a model of the stripping process based on detailed hydrodynamic simulations within the Durham semi-analytic model of galaxy formation. The simulations show that the environment in groups and clusters is less aggressive than previously assumed. The main uncertainty in the model is the treatment of gas expelled by supernovae. With reasonable assumptions for the stripping of this material, we find that satellite galaxies are able to retain a significant fraction of their hot gas for several Gigayears, thereby replenishing their reservoirs of cold, star forming gas and remaining blue for a relatively long period of time. A bimodal distribution of galaxy colours, similar to that observed in SDSS data, is established and the colours of the satellite galaxies are in good agreement with the data. In addition, our model naturally accounts for the observed dependence of satellite colours on environment, from small groups to high mass clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

The Halo Occupation Distribution and the Physics of Galaxy Formation

The halo occupation distribution (HOD) describes the bias between galaxies and dark matter by specifying (a) the probability P(N|M) that a halo of virial mass M contains N galaxies of a particular class and (b) the relative distributions of galaxies and dark matter within halos. We calculate predicted HODs for a Lambda-CDM cosmological model using an SPH hydrodynamic simulation and a semi-analytic (SA) galaxy formation model. Although the two methods predict different galaxy mass functions, their HOD predictions agree remarkably well. For mass-selected samples, the mean occupation <N(M)> exhibits a sharp cutoff at low halo masses, a slowly rising plateau for <N>~1-2, and a more steeply rising high occupancy regime. At low <N>, the mean pair and triple counts are well below Poisson expectations, with important consequences for small scale behavior of 2- and 3-point correlation functions. The HOD depends strongly on galaxy age, with high mass halos populated mainly by old galaxies and low mass halos by young galaxies. The SPH simulation supports several simplifying assumptions about HOD bias: the most massive galaxy in a halo usually lies close to the center and moves near the halo's mean velocity; satellite galaxies have the same radial profile and velocity dispersion as the dark matter; and the mean occupation at fixed halo mass is independent of the halo's larger scale environment. By applying the SPH and SA HODs to a large volume N-body simulation, we show that both methods predict slight, observable departures from a power-law galaxy correlation function. The predicted HODs are closely tied to the underlying galaxy formation physics, they offer useful guidance to theoretical models of galaxy clustering, and they will be tested empirically by ongoing analyses of galaxy redshift surveys. (Shortened)

0 citations0 reviewsNo code linkedOpen access