Source author record

Cedric Lacey

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

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

14 published item(s)

preprint2022arXiv

The spatial distribution of satellites in galaxy clusters

The planar distributions of satellite galaxies around the Milky Way and Andromeda have been extensively studied as potential challenges to the standard cosmological model. Using the Sloan Digital Sky Survey and the Millennium simulation we extend such studies to the satellite galaxies of massive galaxy clusters. We find that both observations and simulations of galaxy clusters show an excess of anisotropic satellite distributions. On average, satellites in clusters have a higher degree of anisotropy than their counterparts in Milky-Way-mass hosts once we account for the difference in their radial distributions. The normal vector of the plane of satellites is strongly aligned with the host halo's minor axis, while the alignment with the large-scale structure is weak. At fixed cluster mass, the degree of anisotropy is higher at higher redshift. This reflects the highly anisotropic nature of satellites accretion points, a feature that is partly erased by the subsequent orbital evolution of the satellites. We also find that satellite galaxies are mostly accreted singly so group accretion is not the explanation for the high flattening of the planes of satellites.

preprint2020arXiv

Qwind code release: a non-hydrodynamical approach to modelling line-driven winds in active galactic nuclei

Ultraviolet (UV) line driven winds may be an important part of the active galactic nucleus (AGN) feedback process, but understanding their impact is hindered by the complex nature of the radiation hydrodynamics. Instead, we have taken the approach pioneered by Risaliti & Elvis, calculating only ballistic trajectories from radiation forces and gravity, but neglecting gas pressure. We have completely re-written their QWIND code using more robust algorithms, and can now quickly model the acceleration phase of these winds for any AGN spectral energy distribution spanning UV and X-ray wavebands. We demonstrate the code using an AGN with black hole mass $10^8\, M_\odot$ emitting at half the Eddington rate and show that this can effectively eject a wind with velocities $\simeq (0.1-0.2)\, c$. The mass loss rates can be up to $\simeq 0.3 M_\odot$ per year, consistent with more computationally expensive hydrodynamical simulations, though we highlight the importance of future improvements in radiation transfer along the multiple different lines of sight illuminating the wind. The code is fully public, and can be used to quickly explore the conditions under which AGN feedback can be dominated by accretion disc winds.

preprint2015arXiv

The evolution of the stellar mass versus halo mass relationship

We present an analysis of the predictions made by the Galform semi-analytic galaxy formation model for the evolution of the relationship between stellar mass and halo mass. We show that for the standard implementations of supernova feedback and gas reincorporation used in semi-analytic models, this relationship is predicted to evolve weakly over the redshift range 0<z<4. Modest evolution in the median stellar mass versus halo mass (SHM) relationship implicitly requires that, at fixed halo mass, the efficiency of stellar mass assembly must be almost constant with cosmic time. We show that in our model, this behaviour can be understood in simple terms as a result of a constant efficiency of gas reincorporation, and an efficiency of SNe feedback that is, on average, constant at fixed halo mass. We present a simple explanation of how feedback from active galactic nuclei (AGN) acts in our model to introduce a break in the SHM relation whose location is predicted to evolve only modestly. Finally, we show that if modifications are introduced into the model such that, for example, the gas reincorporation efficiency is no longer constant, the median SHM relation is predicted to evolve significantly over 0<z<4. Specifically, we consider modifications that allow the model to better reproduce either the evolution of the stellar mass function or the evolution of average star formation rates inferred from observations.

preprint2015arXiv

The Star Formation and AGN luminosity relation: Predictions from a semi-analytical model

In a Universe where AGN feedback regulates star formation in massive galaxies, a strong correlation between these two quantities is expected. If the gas causing star formation is also responsible for feeding the central black hole, then a positive correlation is expected. If powerful AGNs are responsible for the star formation quenching, then a negative correlation is expected. Observations so far have mainly found a mild correlation or no correlation at all (i.e. a flat relation between star formation rate (SFR) and AGN luminosity), raising questions about the whole paradigm of "AGN feedback". In this paper, we report the predictions of the GALFORM semi-analytical model, which has a very strong coupling between AGN activity and quenching of star formation. The predicted SFR-AGN luminosity correlation appears negative in the low AGN luminosity regime, where AGN feedback acts, but becomes strongly positive in the regime of the brightest AGN. Our predictions reproduce reasonably well recent observations by Rosario et al., yet there is some discrepancy in the normalisation of the correlation at low luminosities and high redshifts. Though this regime could be strongly influenced by observational biases, we argue that the disagreement could be ascribed to the fact that GALFORM neglects AGN variability effects. Interestingly, the galaxies that dominate the regime where the observations imply a weak correlation are massive early-type galaxies that are subject to AGN feedback. Nevertheless, these galaxies retain high enough molecular hydrogen contents to maintain relatively high star formation rates and strong infrared emission.

preprint2014arXiv

Evolution of the dust emission of massive galaxies up to z=4 and constraints on their dominant mode of star formation

We aim to measure the average dust and molecular gas content of massive star-forming galaxies ($\rm > 3 \times 10^{10}\,M_\odot$) up to z=4 in the COSMOS field to determine if the intense star formation observed at high redshift is induced by major mergers or caused by large gas reservoirs. Firstly, we measured the evolution of the average spectral energy distributions as a function of redshift using a stacking analysis of Spitzer, Herschel, LABOCA, and AzTEC data for two samples of galaxies: normal star-forming objects and strong starbursts, as defined by their distance to the main sequence. We found that the mean intensity of the radiation field $< U >$ heating the dust (strongly correlated with dust temperature) increases with increasing redshift up to z$\sim$4 in main-sequence galaxies. We can reproduce this evolution with simple models that account for the decrease of the gas metallicity with redshift. No evolution of $< U >$ with redshift is found in strong starbursts. We then deduced the evolution of the molecular gas fraction (defined here as $\rm M_{\rm mol}/(M_{\rm mol}+M_\star)$) with redshift and found a similar, steeply increasing trend for both samples. At z$\sim$4, this fraction reaches $\sim$60%. The average position of the main-sequence galaxies is on the locus of the local, normal star-forming disks in the integrated Schmidt-Kennicutt diagram (star formation rate versus mass of molecular gas), suggesting that the bulk of the star formation up to z=4 is dominated by secular processes.

preprint2014arXiv

Galaxy And Mass Assembly (GAMA): The dependence of the galaxy luminosity function on environment, redshift and colour

We use 80922 galaxies in the Galaxy And Mass Assembly (GAMA) survey to measure the galaxy luminosity function (LF) in different environments over the redshift range 0.04<z<0.26. The depth and size of GAMA allows us to define samples split by colour and redshift to measure the dependence of the LF on environment, redshift and colour. We find that the LF varies smoothly with overdensity, consistent with previous results, with little environmental dependent evolution over the last 3 Gyrs. The modified GALFORM model predictions agree remarkably well with our LFs split by environment, particularly in the most overdense environments. The LFs predicted by the model for both blue and red galaxies are consistent with GAMA for the environments and luminosities at which such galaxies dominate. Discrepancies between the model and the data seen in the faint end of the LF suggest too many faint red galaxies are predicted, which is likely to be due to the over-quenching of satellite galaxies. The excess of bright blue galaxies predicted in underdense regions could be due to the implementation of AGN feedback not being sufficiently effective in the lower mass halos.

preprint2014arXiv

Herschel-ATLAS/GAMA:How does the far-IR luminosity function depend on galaxy group properties?

We use the Herschel ATLAS (H-ATLAS) Phase I data to study the conditional luminosity function of far-IR (250 um) selected galaxies in optically-selected galaxy groups from the GAMA spectroscopic survey, as well as environmental effects on the far-IR-to-optical colour. We applied two methods, which gave consistent results for the far-IR conditional luminosity functions. The direct matching method matches H-ATLAS sources to GAMA/SDSS galaxies, then links the optical counterparts to GAMA groups. The stacking method counts the number of far-IR sources within the projected radii of GAMA groups, subtracting the local background. We investigated the dependence of the far-IR (250 um) luminosity function on group mass in the range 10d12 < M_h < 10d14Msun/h and on redshift in the range 0 < z < 0.4, using a sample of 3000 groups containing H-ATLAS sources with GAMA redshifts over an area of 126 sq. deg. We find that the characteristic 250 um luminosity, L^*(250), increases with group mass up to M_h ~ 10d13Msun/h, but is roughly constant above this, while it increases with redshift at high group masses, but less so at low masses. We also find that the group far-IR luminosity-to-mass ratio L(250)/M_h increases with redshift and is higher in low-mass groups. We estimate that around 70% of the 250 um luminosity density in the local universe is contributed by groups with M_h > 10d12Msun/h. We also find that the far-IR-to-optical colours of H-ATLAS galaxies are independent of group mass over the range 10d12 < M_h <10d14Msun/h in the local universe. We also compare our observational results with recent semi-analytical models, and find that none of these galaxy formation model can reproduce the conditional far-IR luminosity functions of galaxy groups.

preprint2012arXiv

Effects of cosmological model assumptions on galaxy redshift survey measurements

The clustering of galaxies observed in future redshift surveys will provide a wealth of cosmological information. Matching the signal at different redshifts constrains the dark energy driving the acceleration of the expansion of the Universe. In tandem with these geometrical constraints, redshift-space distortions (RSD) depend on the build up of large-scale structure. As pointed out by many authors measurements of these effects are intrinsically coupled. We investigate this link, and argue that it strongly depends on the cosmological assumptions adopted when analysing data. Using representative assumptions for the parameters of the Euclid survey in order to provide a baseline future experiment, we show how the derived constraints change due to different model assumptions. We argue that even the assumption of a Friedman-Robertson-Walker (FRW) space-time is sufficient to reduce the importance of the coupling to a significant degree. Taking this idea further, we consider how the data would actually be analysed and argue that we should not expect to be able to simultaneously constrain multiple deviations from the standard $Λ$CDM model. We therefore consider different possible ways in which the Universe could deviate from the $Λ$CDM model, and show how the coupling between geometrical constraints and structure growth affects the measurement of such deviations.

preprint2010arXiv

The galaxies that reionized the Universe

The Durham GALFORM semi-analytical galaxy formation model has been shown to reproduce the observed rest-frame 1500Å luminosity function of galaxies well over the whole redshift range z=5-10. We show that in this model, this galaxy population also emits enough ionizing photons to reionize the Universe by redshift z=10, assuming a modest escape fraction of 20 per cent. The bulk of the ionizing photons is produced in faint galaxies during starbursts triggered by galaxy mergers. The bursts introduce a dispersion up to ~ 5 dex in galaxy ionizing luminosity at a given halo mass. Almost 90 per cent of the ionizing photons emitted at z=10 are from galaxies below the current observational detection limit at that redshift. Photo-ionization suppression of star formation in these galaxies is unlikely to affect this conclusion significantly, because the gas that fuels the starbursts has already cooled out of their host halos. The galaxies that dominate the ionizing emissivity at z=10 are faint, with M_{1500, AB} ~ -16, have low star formation rates, \dot{M_{*}} ~ 0.06 h^{-1} M_sun yr^{-1}, and reside in halos of mass M ~ 10^9 h^{-1} M_sun.

preprint2010arXiv

The Herschel-SPIRE Legacy Survey (HSLS): the scientific goals of a shallow and wide submillimeter imaging survey with SPIRE

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

preprint2008arXiv

The Properties of Sub-mm Galaxies in Hierarchical Models

We use the combined GALFORM semi-analytical model of galaxy formation and GRASIL spectrophotometric code to investigate the properties of galaxies selected via their sub-mm emission. Our fiducial model has previously been shown to fit the properties of local ULIRGs, as well as the number counts of faint sub-mm galaxies. Here, we test the model in detail by comparing the SEDs and stellar, dynamical, gas and halo masses of sub-mm galaxies against observational data. We precisely mimic the sub-mm and radio selection function of the observations and show that the predicted far-infrared properties of model galaxies with S_850>5mJy and S_1.4>30uJy are in good agreement with observations. Although the dust emission model does not assume a single dust temperature, the far-infrared SEDs are well described by single component modified black-body spectrum with characteristic temperature 32+/-5K. We also find evidence that the observations may have uncovered evolution in the far-infrared--radio relation in ULIRGs out to z~2. We show that the predicted redshift distribution of sub-mm galaxies provides a reasonable fit to the observational data with a median redshift z=2.0, with the radio-selected subset predicted to make up approximately 75% of the population. However, the predicted K-band and mid-infrared (3--8um) flux densities of the sub-mm galaxies (and LBGs) are up to a factor 10x fainter than observed. This discrepancy may indicate that the stellar masses of the sub-mm galaxies in the model are too low: M~10^10Mo, while observations suggest more massive systems, M~10^11Mo. Finally, we discuss the potential modifications to the models which may improve the fit to the observational data. [Abridged]

preprint2004arXiv

The abundance of Lyman-(α) emitters in hierarchical models

We present predictions for the abundance of Ly-$α$ emitters in hierarchical structure formation models. We use the {\tt GALFORM} semi-analytical model to explore the impact on the predicted counts of varying assumptions about the escape fraction of Ly-$α$ photons, the redshift at which the universe reionised and the cosmological density parameter. A model with a fixed escape fraction gives a remarkably good match to the observed counts over a wide redshift interval. We present predictions for the expected counts in a typical observation with the Multi Unit Spectroscopic Explorer instrument proposed for the Very Large Telescope.

preprint1996arXiv

The Structure of Dark Matter Haloes in Hierarchical Clustering Models

We use a set of large cosmological N-body simulations to study the internal structure of dark matter haloes which form in scale-free models. We find that the radius r_178 corresponding to a mean interior overdensity of 178 accurately delineates the quasi-static halo interior from the surrounding infalling material, in agreement with the simple spherical collapse model. The interior velocity dispersion correlates with mass, again in good agreement with the spherical collapse model. Interior to the virial radius r_178, the spherical averaged density, circular velocity and velocity dispersion profiles are well fit by a simple 2-parameter analytic model proposed by Navarro etal (1995). This model has density going as 1/r at small radii, steepening to 1/r^3 at large radii, and fits our haloes to the resolution limit of the simulations. The two model parameters, scalelength and mass, are tightly correlated. Lower mass haloes are more centrally concentrated, and so have scalelengths which are a smaller fraction of their virial radius than those of their higher mass counterparts. This reflects the earlier formation times of low mass haloes. The haloes are moderately aspherical, with typical axial ratios 1:0.8:0.65 at their virial radii. These shapes are maintained by an anisotropic velocity dispersion tensor. The median value of the spin parameter is lambda=0.04, with a weak trend for lower lambda at higher halo mass. We also investigate how the halo properties depend on the algorithm used to identify them, using both friends-of-friends and spherical overdensity methods. We find that for groups selected at mean overdensities 100-400 by either method the properties are insensitive to how the haloes are selected, if the halo centre is taken as the position of the most bound particle

preprint1994arXiv

Merger rates in hierarchical models of galaxy formation. II: Comparison with N-body simulations

We have made a detailed comparison of the results of large N-body simulations with the analytical description of the merging histories of dark matter halos presented in Lacey & Cole 1993, which is based on an extension of the Press- Schechter method (Bond etal 1991,Bower 1991). We find the analytical predictions for the halo mass function, merger rates and formation times to be remarkably accurate. The analytical model is however expected to apply for arbitrary Omega and general power spectra. Dark matter halos were identified in the simulations using two different methods. For halos selected at mean overdensities 100-200, the analytical mass function was found to provide a good fit to the simulations with a collapse threshold close to that predicted by the spherical collapse model, with a typical error of <30% over a range of 1000 in mass, which is the full dynamic range of our N-body simulations. This was insensitive to the type of filtering used. Over a range of 100 - 1000 in mass, there was also good agreement with the analytical predictions for merger rates, including their dependence on the masses of the two halos involved and the time interval being considered, and for formation times, including the dependence on halo mass and formation epoch. The analytical Press-Schechter mass function and its extension to halo lifetimes and merger rates thus provide a very useful description of the growth of dark matter halos through hierarchical clustering, and should providea valuable tool in studies of the formation and evolution of galaxies and galaxy clusters.