Source author record

Richard G. Bower

Richard G. Bower 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.GA astro-ph astro-ph.IM astro-ph.HE physics.hist-ph

Catalog footprint

What is connected

76works

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.

preprint2022arXiv

The importance of black hole repositioning for galaxy formation simulations

Active galactic nucleus (AGN) feedback from accreting supermassive black holes (SMBHs) is an essential ingredient of galaxy formation simulations. The orbital evolution of SMBHs is affected by dynamical friction that cannot be predicted self-consistently by contemporary simulations of galaxy formation in representative volumes. Instead, such simulations typically use a simple "repositioning" of SMBHs, but the effects of this approach on SMBH and galaxy properties have not yet been investigated systematically. Based on a suite of smoothed particle hydrodynamics simulations with the SWIFT code and a Bondi-Hoyle-Lyttleton subgrid gas accretion model, we investigate the impact of repositioning on SMBH growth and on other baryonic components through AGN feedback. Across at least a factor ~1000 in mass resolution, SMBH repositioning (or an equivalent approach) is a necessary prerequisite for AGN feedback; without it, black hole growth is negligible. Limiting the effective repositioning speed to $\lesssim$ 10 km/s delays the onset of AGN feedback and severely limits its impact on stellar mass growth in the centre of massive galaxies. Repositioning has three direct physical consequences. It promotes SMBH mergers and thus accelerates their initial growth. In addition, it raises the peak density of the ambient gas and reduces the SMBH velocity relative to it, giving a combined boost to the accretion rate that can reach many orders of magnitude. Our results suggest that a more sophisticated and/or better calibrated treatment of SMBH repositioning is a critical step towards more predictive galaxy formation simulations.

preprint2020arXiv

Galactic inflow and wind recycling rates in the EAGLE simulations

The role of galactic wind recycling represents one of the largest unknowns in galaxy evolution, as any contribution of recycling to galaxy growth is largely degenerate with the inflow rates of first-time infalling material, and the rates with which outflowing gas and metals are driven from galaxies. We present measurements of the efficiency of wind recycling from the EAGLE cosmological simulation project, leveraging the statistical power of large-volume simulations that reproduce a realistic galaxy population. We study wind recycling at the halo scale, i.e. gas that has been ejected beyond the halo virial radius, and at the galaxy scale, i.e. gas that has been ejected from the ISM to at least $\approx 10 \, \%$ of the virial radius (thus excluding smaller-scale galactic fountains). Galaxy-scale wind recycling is generally inefficient, with a characteristic return timescale that is comparable or longer than a Hubble time, and with an efficiency that clearly peaks at the characteristic halo mass of $M_{200} = 10^{12} \, \mathrm{M_\odot}$. Correspondingly, the majority of gas being accreted onto galaxies in EAGLE is infalling for the first time. At the halo scale, the efficiency of recycling onto haloes differs by orders of magnitude from values assumed by semi-analytic galaxy formation models. Differences in the efficiency of wind recycling with other hydrodynamical simulations are currently difficult to assess, but are likely smaller. We are able to show that the fractional contribution of wind recycling to galaxy growth is smaller in EAGLE than in some other simulations. We find that cumulative first-time gas accretion rates at the virial radius are reduced relative to the expectation from dark matter accretion for haloes with mass, $M_{200} < 10^{12} \, \mathrm{M_\odot}$, indicating efficient preventative feedback on halo scales.

preprint2020arXiv

Galactic outflow rates in the EAGLE simulations

We present measurements of galactic outflow rates from the EAGLE suite of cosmological simulations. We find that gas is removed from the interstellar medium (ISM) of central galaxies with a dimensionless mass loading factor that scales approximately with circular velocity as $V_{\mathrm{c}}^{-3/2}$ in the low-mass regime where stellar feedback dominates. Feedback from active galactic nuclei (AGN) causes an upturn in the mass loading for halo masses $> 10^{12} \, \mathrm{M_\odot}$. We find that more gas outflows through the halo virial radius than is removed from the ISM of galaxies, particularly at low redshift, implying substantial mass loading within the circum-galactic medium (CGM). Outflow velocities span a wide range at a given halo mass/redshift, and on average increase positively with redshift and halo mass up to $M_{200} \sim 10^{12} \, \mathrm{M_\odot}$. Outflows exhibit a bimodal flow pattern on circum-galactic scales, aligned with the galactic minor axis. We present a number of like-for-like comparisons to outflow rates from other recent cosmological hydrodynamical simulations, and show that comparing the propagation of galactic winds as a function of radius reveals substantial discrepancies between different models. Relative to some other simulations, EAGLE favours a scenario for stellar feedback where agreement with the galaxy stellar mass function is achieved by removing smaller amounts of gas from the ISM, but with galactic winds that then propagate and entrain ambient gas out to larger radii.

preprint2019arXiv

How feedback shapes galaxies: an analytic model

We introduce a simple analytic model of galaxy formation that links the growth of dark matter haloes in a cosmological background to the build-up of stellar mass within them. The model aims to identify the physical processes that drive the galaxy-halo co-evolution through cosmic time. The model restricts the role of baryonic astrophysics to setting the relation between galaxies and their haloes. Using this approach, galaxy properties can be directly predicted from the growth of their host dark matter haloes. We explore models in which the effective star formation efficiency within haloes is a function of mass (or virial temperature) and independent of time. Despite its simplicity, the model reproduces self-consistently the shape and evolution of the cosmic star formation rate density, the specific star formation rate of galaxies, and the galaxy stellar mass function, both at the present time and at high redshifts. By systematically varying the effective star formation efficiency in the model, we explore the emergence of the characteristic shape of the galaxy stellar mass function. The origin of the observed double Schechter function at low redshifts is naturally explained by two efficiency regimes in the stellar to halo mass relation, namely, a stellar feedback regulated stage, and a supermassive black hole regulated stage. By providing a set of analytic differential equations, the model can be easily extended and inverted, allowing the roles and impact of astrophysics and cosmology to be explored and understood.

preprint2019arXiv

Hydrostatic mass estimates of massive galaxy clusters: a study with varying hydrodynamics flavours and non-thermal pressure support

We use a set of 45 simulated clusters with a wide mass range ($8\times 10^{13} < M_{500}~[$M$_{\odot}]~< 2\times 10^{15}$) to investigate the effect of varying hydrodynamics flavours on cluster mass estimates. The cluster zooms were simulated using the same cosmological models as the BAHAMAS and C-EAGLE projects, leading to differences in both the hydrodynamic solvers and the subgrid physics but still producing clusters which broadly match observations. At the same mass resolution as BAHAMAS, for the most massive clusters ($M_{500} > 10^{15}$ M$_{\odot}$), we find changes in the SPH method produce the greatest differences in the final halo, while the subgrid models dominate at lower mass. By calculating the mass of all of the clusters using different permutations of the pressure, temperature and density profiles, created with either the true simulated data or mock spectroscopic data, we find that the spectroscopic temperature causes a bias in the hydrostatic mass estimates which increases with the mass of the cluster, regardless of the SPH flavour used. For the most massive clusters, the estimated mass of the cluster using spectroscopic density and temperature profiles is found to be as low as 50 per cent of the true mass compared to $\sim$ 90 per cent for low mass clusters. When including a correction for non-thermal pressure, the spectroscopic hydrostatic mass estimates are less biased on average and the mass dependence of the bias is reduced, although the scatter in the measurements does increase.

preprint2016arXiv

A chronicle of galaxy mass assembly in the EAGLE simulation

We analyse the mass assembly of central galaxies in the EAGLE hydrodynamical simulations. We build merger trees to connect galaxies to their progenitors at different redshifts and characterize their assembly histories by focusing on the time when half of the galaxy stellar mass was assembled into the main progenitor. We show that galaxies with stellar mass $M_*<10^{10.5}M_{\odot}$ assemble most of their stellar mass through star formation in the main progenitor (`in-situ' star formation). This can be understood as a consequence of the steep rise in star formation efficiency with halo mass for these galaxies. For more massive galaxies, however, an increasing fraction of their stellar mass is formed outside the main progenitor and subsequently accreted. Consequently, while for low-mass galaxies the assembly time is close to the stellar formation time, the stars in high-mass galaxies typically formed long before half of the present-day stellar mass was assembled into a single object, giving rise to the observed anti-hierarchical downsizing trend. In a typical present-day $M_*\geq10^{11}M_{\odot}$ galaxy, around $20\%$ of the stellar mass has an external origin. This fraction decreases with increasing redshift. Bearing in mind that mergers only make an important contribution to the stellar mass growth of massive galaxies, we find that the dominant contribution comes from mergers with galaxies of mass greater than one tenth of the main progenitor's mass. The galaxy merger fraction derived from our simulations agrees with recent observational estimates.

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]

preprint2016arXiv

Being WISE II: Reducing the Influence of Star Formation History on the Mass-to-Light Ratio of Quiescent Galaxies

Stellar population synthesis models can now reproduce the photometry of old stellar systems (age $>$ 2 Gyr) in the near-infrared (NIR) bands at 3.4 and 4.6$μ$m (WISE W1 $\&$ W2 or IRAC 1 $\&$ 2). In this paper we derive stellar mass-to-light ratios for these and optical bands, and confirm that the NIR M/L shows dramatically reduced sensitivity to both age and metallicity compared to optical bands, and further, that this behavior leads to significantly more robust stellar masses for quiescent galaxies with [Fe/H] > -0.5 regardless of star formation history (SFH). We then use realistic early-type galaxy SFHs and metallicity distributions from the EAGLE simulations of galaxy formation to investigate two methods to determine the appropriate M/L for a galaxy: 1) We show that the uncertainties introduced by an unknown SFH can be largely removed using a spectroscopically inferred luminosity-weighted age and metallicity for the population to select the appropriate single stellar population (SSP) equivalent M/L. Using this method, the maximum systematic error due to SFH on the M/L of an early-type galaxy is $<$ 4$\%$ at 3.4 $μ$m and typical uncertainties due to errors in the age and metallicity create scatter of $\lesssim$13$\%$. The equivalent values for optical bands are more than 2-3 times greater, even before considering uncertainties associated with internal dust extinction. 2) We demonstrate that if the EAGLE SFHs and metallicities accurately reproduce the true properties of early-type galaxies, the use of an iterative approach to select a mass dependent M/L can provide even more accurate stellar masses for early-type galaxies, with typical uncertainties $<$ 9$\%$.

preprint2016arXiv

Bimodality of low-redshift circumgalactic O VI in non-equilibrium EAGLE zoom simulations

We introduce a series of 20 cosmological hydrodynamical simulations of Lstar (M_200 =10^11.7 - 10^12.3 Msol) and group-sized (M_200 = 10^12.7 - 10^13.3 Msol) haloes run with the model used for the EAGLE project, which additionally includes a non-equilibrium ionization and cooling module that follows 136 ions. The simulations reproduce the observed correlation, revealed by COS-Halos at z~0.2, between O VI column density at impact parameters b < 150 kpc and the specific star formation rate (sSFR=SFR/Mstar) of the central galaxy at z~0.2. We find that the column density of circumgalactic O VI is maximal in the haloes associated with Lstar galaxies, because their virial temperatures are close to the temperature at which the ionization fraction of O VI peaks (T~10^5.5 K). The higher virial temperature of group haloes (> 10^6 K) promotes oxygen to higher ionization states, suppressing the O VI column density. The observed NO VI-sSFR correlation therefore does not imply a causal link, but reflects the changing characteristic ionization state of oxygen as halo mass is increased. In spite of the mass-dependence of the oxygen ionization state, the most abundant circumgalactic oxygen ion in both Lstar and group haloes is O VII; O VI accounts for only 0.1% of the oxygen in group haloes and 0.9-1.3% with Lstar haloes. Nonetheless, the metals traced by O VI absorbers represent a fossil record of the feedback history of galaxies over a Hubble time; their characteristic epoch of ejection corresponds to z > 1 and much of the ejected metal mass resides beyond the virial radius of galaxies. For both Lstar and group galaxies, more of the oxygen produced and released by stars resides in the circumgalactic medium (within twice the virial radius) than in the stars and ISM of the galaxy.

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.

preprint2016arXiv

It's not easy being green: The evolution of galaxy colour in the EAGLE simulation

We examine the evolution of intrinsic u-r colours of galaxies in the EAGLE cosmological hydrodynamical simulations, which has been shown to reproduce the observed redshift z=0.1 colour-magnitude distribution well. The median u-r of star-forming ('blue cloud') galaxies reddens by 1 mag from z=2 to 0 at fixed stellar mass, as their specific star formation rates decrease with time. A red sequence starts to build-up around z=1, due to the quenching of low-mass satellite galaxies at the faint end, and due to the quenching of more massive central galaxies by their active galactic nuclei (AGN) at the bright end. This leaves a dearth of intermediate-mass red sequence galaxies at z=1, which is mostly filled in by z=0. We quantify the time-scales of colour transition due to satellite and AGN quenching, finding that most galaxies spend less than 2 Gyr in the 'green valley'. On examining the trajectories of galaxies in a colour-stellar mass diagram, we identify three characteristic tracks that galaxies follow (quiescently star-forming, quenching and rejuvenating galaxies) and quantify the fraction of galaxies that follow each track.

preprint2016arXiv

Music from the heavens - Gravitational waves from supermassive black hole mergers in the EAGLE simulations

We estimate the expected event rate of gravitational wave signals from mergers of supermassive black holes that could be resolved by a space-based interferometer, such as the Evolved Laser Interferometer Space Antenna (eLISA), utilising the reference cosmological hydrodynamical simulation from the EAGLE suite. These simulations assume a $Λ$CDM cosmogony with state-of-the-art subgrid models for radiative cooling, star formation, stellar mass loss, and feedback from stars and accreting black holes. They have been shown to reproduce the observed galaxy population with unprecedented fidelity. We combine the merger rates of supermassive black holes in EAGLE with the latest phenomenological waveform models to calculate the gravitational waves signals from the intrinsic parameters of the merging black holes. The EAGLE models predict $\sim2$ detections per year by a gravitational wave detector such as eLISA. We find that these signals are largely dominated by mergers between seed mass black holes merging at redshifts between $z\sim2$ and $z\sim1$. In order to investigate the dependence on the assumed black hole seed mass, we introduce an additional model with a black hole seed mass an order of magnitude smaller than in our reference model. We also consider a variation of the reference model where a prescription for the expected delays in the black hole merger timescale has been included after their host galaxies merge. We find that the merger rate is similar in all models, but that the initial black hole seed mass could be distinguished through their detected gravitational waveforms. Hence, the characteristic gravitational wave signals detected by eLISA will provide profound insight into the origin of supermassive black holes and the initial mass distribution of black hole seeds.

preprint2016arXiv

Size matters: abundance matching, galaxy sizes, and the Tully-Fisher relation in EAGLE

The Tully-Fisher relation (TFR) links the stellar mass of a disk galaxy, $M_{\rm str}$, to its rotation speed: it is well approximated by a power law, shows little scatter, and evolves weakly with redshift. The relation has been interpreted as reflecting the mass-velocity scaling ($M\propto V^3$) of dark matter halos, but this interpretation has been called into question by abundance-matching (AM) models, which predict the galaxy-halo mass relation to be non-monotonic and rapidy evolving. We study the TFR of luminous spirals and its relation to AM using the EAGLE set of $Λ$CDM cosmological simulations. Matching both relations requires disk sizes to satisfy constraints given by the concentration of halos and their response to galaxy assembly. EAGLE galaxies approximately match these constraints and show a tight mass-velocity scaling that compares favourably with the observed TFR. The TFR is degenerate to changes in galaxy formation efficiency and the mass-size relation; simulations that fail to match the galaxy stellar mass function may fit the observed TFR if galaxies follow a different mass-size relation. The small scatter in the simulated TFR results because, at fixed halo mass, galaxy mass and rotation speed correlate strongly, scattering galaxies along the main relation. EAGLE galaxies evolve with lookback time following approximately the prescriptions of AM models and the observed mass-size relation of bright spirals, leading to a weak TFR evolution consistent with observation out to $z=1$. $Λ$CDM models that match both the abundance and size of galaxies as a function of stellar mass have no difficulty reproducing the observed TFR and its evolution.

preprint2016arXiv

Supermassive black holes in the EAGLE Universe. Revealing the observables of their growth

We investigate the evolution of supermassive black holes in the `Evolution and Assembly of GaLaxies and their Environments' (EAGLE) cosmological hydrodynamic simulations. The largest of the EAGLE volumes covers a $(100 \,\rm cMpc)^3$ and includes state-of-the-art physical models for star formation and black hole growth that depend only on local gas properties. We focus on the black hole mass function, Eddington ratio distribution and the implied duty cycle of nuclear activity. The simulation is broadly consistent with observational constraints on these quantities. In order to make a more direct comparison with observational data, we calculate the soft and hard X-ray luminosity functions of the active galactic nuclei (AGN). Between redshifts $0$ and $1$, the simulation is in agreement with data. At higher redshifts, the simulation tends to underpredict the luminosities of the brightest observed AGN. This may be due to the limited volume of the simulation, or a fundamental deficiency of the underlying model. It seems unlikely that additional unresolved variability can account for this difference. The simulation shows a similar `downsizing' of the AGN population as seen in observational surveys.

preprint2016arXiv

The brighter galaxies reionised the Universe

Hydrogen in the Universe was (re)ionised between redshifts $z \approx 10$ and $z \approx 6$. The nature of the sources of the ionising radiation is hotly debated, with faint galaxies below current detection limits regarded as prime candidates. Here we consider a scenario in which ionising photons escape through channels punctured in the interstellar medium by outflows powered by starbursts. We take account of the observation that strong outflows occur only when the star formation density is sufficiently high, and estimate the galaxy-averaged escape fraction as a function of redshift and luminosity from the resolved star formation surface densities in the EAGLE cosmological hydrodynamical simulation. We find that the fraction of ionising photons that escape from galaxies increases rapidly with redshift, reaching values of 5-20 percent at $z > 6$, with the brighter galaxies having higher escape fractions. Combining the dependence of escape fraction on luminosity and redshift with the observed luminosity function, we demonstrate that galaxies emit enough ionising photons to match the existing constraints on reionisation while also matching the observed UV-background post-reionisation. Our findings suggest that galaxies above the current Hubble Space Telescope detection limit emit half of the ionising radiation required to reionise the Universe.

preprint2016arXiv

The EAGLE simulations of galaxy formation: public release of halo and galaxy catalogues

We present the public data release of halo and galaxy catalogues extracted from the EAGLE suite of cosmological hydrodynamical simulations of galaxy formation. These simulations were performed with an enhanced version of the GADGET code that includes a modified hydrodynamics solver, time-step limiter and subgrid treatments of baryonic physics, such as stellar mass loss, element-by-element radiative cooling, star formation and feedback from star formation and black hole accretion. The simulation suite includes runs performed in volumes ranging from 25 to 100 comoving megaparsecs per side, with numerical resolution chosen to marginally resolve the Jeans mass of the gas at the star formation threshold. The free parameters of the subgrid models for feedback are calibrated to the redshift z=0 galaxy stellar mass function, galaxy sizes and black hole mass - stellar mass relation. The simulations have been shown to match a wide range of observations for present-day and higher-redshift galaxies. The raw particle data have been used to link galaxies across redshifts by creating merger trees. The indexing of the tree produces a simple way to connect a galaxy at one redshift to its progenitors at higher redshift and to identify its descendants at lower redshift. In this paper we present a relational database which we are making available for general use. A large number of properties of haloes and galaxies and their merger trees are stored in the database, including stellar masses, star formation rates, metallicities, photometric measurements and mock gri images. Complex queries can be created to explore the evolution of more than 10^5 galaxies, examples of which are provided in appendix. (abridged)

preprint2016arXiv

The EAGLE simulations: atomic hydrogen associated with galaxies

We examine the properties of atomic hydrogen (HI) associated with galaxies in the EAGLE simulations of galaxy formation. EAGLE's feedback parameters were calibrated to reproduce the stellar mass function and galaxy sizes at $z=0.1$, and we assess whether this calibration also yields realistic HI properties. We estimate the self-shielding density with a fitting function calibrated using radiation transport simulations, and correct for molecular hydrogen with empirical or theoretical relations. The `standard-resolution' simulations systematically underestimate HI column densities, leading to an HI deficiency in low-mass ($M_\star < 10^{10}M_\odot$) galaxies and poor reproduction of the observed HI mass function. These shortcomings are largely absent from EAGLE simulations featuring a factor of 8 (2) better mass (spatial) resolution, within which the HI mass of galaxies evolves more mildly from $z=1$ to $0$ than in the standard-resolution simulations. The largest-volume simulation reproduces the observed clustering of HI systems, and its dependence on HI-richness. At fixed $M_\star$, galaxies acquire more HI in simulations with stronger feedback, as they become associated with more massive haloes and higher infall rates. They acquire less HI in simulations with a greater star formation efficiency, since the star formation and feedback necessary to balance the infall rate is produced by smaller gas reservoirs. The simulations indicate that the HI of present-day galaxies was acquired primarily by the smooth accretion of ionized, intergalactic gas at $z\simeq1$, which later self-shields, and that only a small fraction is contributed by the reincorporation of gas previously heated strongly by feedback. HI reservoirs are highly dynamic: over $40$ percent of HI associated with $z=0.1$ galaxies is converted to stars or ejected by $z=0$.

preprint2016arXiv

The effect of baryons on redshift space distortions and cosmic density and velocity fields in the EAGLE simulation

We use the EAGLE galaxy formation simulation to study the effects of baryons on the power spectrum of the total matter and dark matter distributions and on the velocity fields of dark matter and galaxies. On scales $k{\stackrel{>}{{}_\sim}} 4{h\,{\rm Mpc}^{-1}}$ the effect of baryons on the amplitude of the total matter power spectrum is greater than $1\%$. The back-reaction of baryons affects the density field of the dark matter at the level of $\sim3\%$ on scales of $1\leq k/({h\,{\rm Mpc}^{-1}})\leq 5$. The dark matter velocity divergence power spectrum at $k{\stackrel{<}{{}_\sim}}0.5{h\,{\rm Mpc}^{-1}}$ is changed by less than $1\%$. The 2D redshift-space power spectrum is affected at the level of $\sim6\%$ at $|\vec{k}|{\stackrel{>}{{}_\sim}} 1{h\,{\rm Mpc}^{-1}}$ (for $μ>0.5$), but for $|\vec{k}|\leq 0.4{h\,{\rm Mpc}^{-1}}$ it differs by less than $1\%$. We report vanishingly small baryonic velocity bias for haloes: the peculiar velocities of haloes with $M_{200}>3\times10^{11}{{\rm M}_{\odot}}$ (hosting galaxies with $M_{*}>10^9{{\rm M}_{\odot}}$) are affected at the level of at most $1~$km/s, which is negligible for $1\%$-precision cosmology. We caution that since EAGLE overestimates cluster gas fractions it may also underestimate the impact of baryons, particularly for the total matter power spectrum. Nevertheless, our findings suggest that for theoretical modelling of redshift space distortions and galaxy velocity-based statistics, baryons and their back-reaction can be safely ignored at the current level of observational accuracy. However, we confirm that the modelling of the total matter power spectrum in weak lensing studies needs to include realistic galaxy formation physics in order to achieve the accuracy required in the precision cosmology era.

preprint2016arXiv

The environmental dependence of HI in galaxies in the EAGLE simulations

We use the EAGLE suite of cosmological hydrodynamical simulations to study how the HI content of present-day galaxies depends on their environment. We show that EAGLE reproduces observed HI mass-environment trends very well, while semi-analytic models typically overpredict the average HI masses in dense environments. The environmental processes act primarily as an on/off switch for the HI content of satellites with stellar mass Mstar>10^9 Msun. At a fixed Mstar, the fraction of HI-depleted satellites increases with increasing host halo mass M200 in response to stronger environmental effects, while at a fixed M200 it decreases with increasing satellite Mstar as the gas is confined by deeper gravitational potentials. HI-depleted satellites reside mostly, but not exclusively, within the virial radius r200 of their host halo. We investigate the origin of these trends by focussing on three environmental mechanisms: ram pressure stripping by the intra-group medium, tidal stripping by the host halo, and satellite-satellite encounters. By tracking back in time the evolution of the HI-depleted satellites, we find that the most common cause of HI removal is satellite encounters. The timescale for HI removal is typically less than 0.5 Gyr. Tidal stripping occurs in halos of M200<10^{14} Msun within half r200, while the other processes act also in more massive halos, generally within r200. Conversely, we find that ram pressure stripping is the most common mechanism that disturbs the HI morphology of galaxies at redshift z=0. This implies that HI removal due to satellite-satellite interactions occurs on shorter timescales than the other processes.

preprint2016arXiv

The fundamental plane of star formation in galaxies revealed by the EAGLE hydrodynamical simulations

We investigate correlations between different physical properties of star-forming galaxies in the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) cosmological hydrodynamical simulation suite over the redshift range $0\le z\le 4.5$. A principal component analysis reveals that neutral gas fraction ($f_{\rm gas, neutral}$), stellar mass ($M_{\rm stellar}$) and star formation rate (SFR) account for most of the variance seen in the population, with galaxies tracing a two-dimensional, nearly flat, surface in the three-dimensional space of $f_{\rm gas, neutral}-M_{\rm stellar}-\rm SFR$ with little scatter. The location of this plane varies little with redshift, whereas galaxies themselves move along the plane as their $f_{\rm gas, neutral}$ and SFR drop with redshift. The positions of galaxies along the plane are highly correlated with gas metallicity. The metallicity can therefore be robustly predicted from $f_{\rm gas, neutral}$, or from the $M_{\rm stellar}$ and SFR. We argue that the appearance of this "fundamental plane of star formation" is a consequence of self-regulation, with the plane's curvature set by the dependence of the SFR on gas density and metallicity. We analyse a large compilation of observations spanning the redshift range $0\lesssim \rm z\lesssim 2.5$, and find that such a plane is also present in the data. The properties of the observed fundamental plane of star formation are in good agreement with EAGLE's predictions.

preprint2016arXiv

The origin of compact galaxies with anomalously high black hole masses

Observations of local galaxies harbouring supermassive black holes (BHs) of anomalously high mass, M_BH, relative to their stellar mass, M_star, appear to be at odds with simple models of the co-evolution between galaxies and their central BHs. We study the origin of such outliers in a Lambda cold dark matter context using the EAGLE cosmological, hydrodynamical simulation. We find 15 'M_BH(M_star)-outlier' galaxies, defined as having M_BH more than 1.5 dex above the median M_BH(M_star) relation in the simulation, M_{BH,med}. All M_BH(M_star)-outliers are satellite galaxies, typically with M_star ~ 10^10 M_sun and M_BH ~ 10^8 M_sun. They have all become outliers due to a combination of tidal stripping of their outer stellar component acting over several Gyr and early formation times leading to rapid BH growth at high redshift, with the former mechanism being most important for 67 per cent of these outliers. The same mechanisms also cause the M_BH(M_star)-outlier satellites to be amongst the most compact galaxies in the simulation, making them ideal candidates for ultracompact dwarf galaxy progenitors. The 10 most extreme central galaxies found at z=0 (with log_{10}(M_BH/M_{BH,med}) in [1.2, 1.5]) grow rapidly in M_BH to lie well above the present-day M_BH-M_star relation at early times (z > 2), and either continue to evolve parallel to the z=0 relation or remain unchanged until the present day, making them 'relics' of the high-redshift universe. This high-z formation mechanism may help to explain the origin of observed M_BH(M_star)-outliers with extended dark matter haloes and undisturbed morphologies.

preprint2016arXiv

The origin of the $α$-enhancement of massive galaxies

We study the origin of the stellar $α$-element-to-iron abundance ratio, $[α/\mathrm{Fe}]_{\ast}$, of present-day central galaxies, using cosmological, hydrodynamical simulations from the Evolution and Assembly of GaLaxies and their Environments (EAGLE) project. For galaxies with stellar masses of $M_{\ast} > 10^{10.5}$ M$_{\odot}$, $[α/\mathrm{Fe}]_{\ast}$ increases with increasing galaxy stellar mass and age. These trends are in good agreement with observations of early-type galaxies, and are consistent with a `downsizing' galaxy formation scenario: more massive galaxies have formed the bulk of their stars earlier and more rapidly, hence from an interstellar medium that was mostly $α$-enriched by massive stars. In the absence of feedback from active galactic nuclei (AGN), however, $[α/\mathrm{Fe}]_{\ast}$ in $M_{\ast} > 10^{10.5}$ M$_{\odot}$ galaxies is roughly constant with stellar mass and decreases with mean stellar age, extending the trends found for lower-mass galaxies in both simulations with and without AGN. We conclude that AGN feedback can account for the $α$-enhancement of massive galaxies, as it suppresses their star formation, quenching more massive galaxies at earlier times, thereby preventing the iron from longer-lived intermediate-mass stars (supernova Type Ia) from being incorporated into younger stars.

preprint2016arXiv

The origin of the enhanced metallicity of satellite galaxies

Observations of galaxies in the local Universe have shown that both the ionized gas and the stars of satellites are more metal-rich than of equally massive centrals. To gain insight into the connection between this metallicity enhancement and other differences between centrals and satellites, such as their star formation rates, gas content, and growth history, we study the metallicities of >3600 galaxies with M_star > 10^10 M_sun in the cosmological hydrodynamical EAGLE 100 Mpc `Reference' simulation, including ~1500 in the vicinity of galaxy groups and clusters (M_200 >= 10^13 M_sun). The simulation predicts excess gas and stellar metallicities in satellites consistent with observations, except for stellar metallicities at M_star <~ 10^10.2 M_sun where the predicted excess is smaller than observed. The exact magnitude of the effect depends on galaxy selection, aperture, and on whether the metallicity is weighted by stellar mass or luminosity. The stellar metallicity excess in clusters is also sensitive to the efficiency scaling of star formation feedback. We identify stripping of low-metallicity gas from the galaxy outskirts, as well as suppression of metal-poor inflows towards the galaxy centre, as key drivers of the enhancement of gas metallicity. Stellar metallicities in satellites are higher than in the field as a direct consequence of the more metal-rich star forming gas, whereas stripping of stars and suppressed stellar mass growth, as well as differences in accreted vs. in-situ star formation between satellites and the field, are of secondary importance.

preprint2015arXiv

Baryon effects on the internal structure of LCDM halos in the EAGLE simulations

We investigate the internal structure and density profiles of halos of mass $10^{10}-10^{14}~M_\odot$ in the Evolution and Assembly of Galaxies and their Environment (EAGLE) simulations. These follow the formation of galaxies in a $Λ$CDM Universe and include a treatment of the baryon physics thought to be relevant. The EAGLE simulations reproduce the observed present-day galaxy stellar mass function, as well as many other properties of the galaxy population as a function of time. We find significant differences between the masses of halos in the EAGLE simulations and in simulations that follow only the dark matter component. Nevertheless, halos are well described by the Navarro-Frenk-White (NFW) density profile at radii larger than ~5% of the virial radius but, closer to the centre, the presence of stars can produce cuspier profiles. Central enhancements in the total mass profile are most important in halos of mass $10^{12}-10^{13}M_\odot$, where the stellar fraction peaks. Over the radial range where they are well resolved, the resulting galaxy rotation curves are in very good agreement with observational data for galaxies with stellar mass $M_*<5\times10^{10}M_\odot$. We present an empirical fitting function that describes the total mass profiles and show that its parameters are strongly correlated with halo mass.

preprint2015arXiv

Colours and luminosities of z=0.1 simulated galaxies in the EAGLE simulations

We calculate the colours and luminosities of redshift z = 0.1 galaxies from the EAGLE simulation suite using the GALAXEV population synthesis models. We take into account obscuration by dust in birth clouds and diffuse ISM using a two-component screen model, following the prescription of Charlot and Fall. We compare models in which the dust optical depth is constant to models where it depends on gas metallicity, gas fraction and orientation. The colours of EAGLE galaxies for the more sophisticated models are in broad agreement with those of observed galaxies. In particular, EAGLE produces a red sequence of passive galaxies and a blue cloud of star forming galaxies, with approximately the correct fraction of galaxies in each population and with g-r colours within 0.1 magnitudes of those observed. Luminosity functions from UV to NIR wavelengths differ from observations at a level comparable to systematic shifts resulting from a choice between Petrosian and Kron photometric apertures. Despite the generally good agreement there are clear discrepancies with observations. The blue cloud of EAGLE galaxies extends to somewhat higher luminosities than in the data, consistent with the modest underestimate of the passive fraction in massive EAGLE galaxies. There is also a moderate excess of bright blue galaxies compared to observations. The overall level of agreement with the observed colour distribution suggests that EAGLE galaxies at z = 0.1 have ages, metallicities and levels of obscuration that are comparable to those of observed galaxies.

preprint2015arXiv

Evidence for a change in the dominant satellite galaxy quenching mechanism at z=1

We present an analysis of galaxies in groups and clusters at $0.8<z<1.2$, from the GCLASS and GEEC2 spectroscopic surveys. We compute a "conversion fraction" $f_{\rm convert}$ that represents the fraction of galaxies that were prematurely quenched by their environment. For massive galaxies, $M_{\rm star}>10^{10.3}M_\odot$, we find $f_{\rm convert}\sim 0.4$ in the groups and $\sim 0.6$ in the clusters, similar to comparable measurements at $z=0$. This means the time between first accretion into a more massive halo and final star formation quenching is $t_p\sim 2$ Gyr. This is substantially longer than the estimated time required for a galaxy's star formation rate to become zero once it starts to decline, suggesting there is a long delay time during which little differential evolution occurs. In contrast with local observations we find evidence that this delay timescale may depend on stellar mass, with $t_p$ approaching $t_{\rm Hubble}$ for $M_{\rm star}\sim 10^{9.5}M_\odot$. The result suggests that the delay time must not only be much shorter than it is today, but may also depend on stellar mass in a way that is not consistent with a simple evolution in proportion to the dynamical time. Instead, we find the data are well-matched by a model in which the decline in star formation is due to "overconsumption", the exhaustion of a gas reservoir through star formation and expulsion via modest outflows in the absence of cosmological accretion. Dynamical gas removal processes, which are likely dominant in quenching newly accreted satellites today, may play only a secondary role at $z=1$.

preprint2015arXiv

Intrinsic alignments of galaxies in the EAGLE and cosmo-OWLS simulations

We report results for the alignments of galaxies in the EAGLE and cosmo-OWLS simulations as a function of galaxy separation and halo mass. The combination of these hydro-cosmological simulations enables us to span four orders of magnitude in halo mass ($10.7<log_{10}(M_{200}/[h^{-1}M_\odot])<15$) and a large range of separations ($-1<log_{10}(r/[h^{-1}Mpc])< 2$). We focus on two classes of alignments: the orientations of galaxies with respect to either the directions to, or the orientations of, surrounding galaxies. We find that the strength of the alignment is a strongly decreasing function of the distance between galaxies. The orientation-direction alignment can remain significant up to ~100 Mpc, for galaxies hosted by the most massive haloes in our simulations. Galaxies hosted by more massive subhaloes show stronger alignment. At a fixed halo mass, more aspherical or prolate galaxies exhibit stronger alignments. The spatial distribution of satellites is anisotropic and significantly aligned with the major axis of the main host halo. The major axis of satellite galaxies, when all stars are considered, are preferentially aligned towards the centre of the main host halo. The predicted projected direction-orientation alignment, $ε_{g+}(r_{p})$, is in broad agreement with recent observations when only stars within the typical observable extent of a galaxy are used to define galaxy orientations. We find that the orientation-orientation alignment is weaker than the orientation-direction alignment on all scales. Overall, the strength of galaxy alignments depends strongly on the subset of stars that are used to measure the orientations of galaxies and it is always weaker than the alignment of the dark matter haloes. Thus, alignment models that use halo orientation as a direct proxy for galaxy orientation will overestimate the impact of intrinsic alignments on weak lensing analyses.

preprint2015arXiv

Molecular hydrogen abundances of galaxies in the EAGLE simulations

We investigate the abundance of galactic molecular hydrogen (H$_2$) in the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) cosmological hydrodynamic simulations. We assign H$_2$ masses to gas particles in the simulations in post-processing using two different prescriptions that depend on the local dust-to-gas ratio and the interstellar radiation field. Both result in H$_2$ galaxy mass functions that agree well with observations in the local and high-redshift Universe. The simulations reproduce the observed scaling relations between the mass of H$_2$ and the stellar mass, star formation rate and stellar surface density. Towards high edshifts, galaxies in the simulations display larger H$_2$ mass fractions, and correspondingly lower H$_2$ depletion timescales, also in good agreement with observations. The comoving mass density of H$_2$ in units of the critical density, $Ω_{\rm H_2}$, peaks at $z\approx 1.2-1.5$, later than the predicted peak of the cosmic star formation rate activity, at $z\approx 2$. This difference stems from the decrease in gas metallicity and increase in interstellar radiation field with redshift, both of which hamper H$_2$ formation. We find that the cosmic H$_2$ budget is dominated by galaxies with $M_{\rm H_2}>10^9\,\rm M_{\odot}$, star formation rates $>10\,\rm M_{\odot}\,\rm yr^{-1}$ and stellar masses $M_{\rm stellar}>10^{10}\,\rm M_{\odot}$, which are readily observable in the optical and near-IR. The match between the H$_2$ properties of galaxies that emerge in the simulations and observations is remarkable, particularly since H$_2$ observations were not used to adjust parameters in EAGLE.

preprint2015arXiv

Recycled stellar ejecta as fuel for star formation and implications for the origin of the galaxy mass-metallicity relation

We use cosmological, hydrodynamical simulations from the EAGLE and OWLS projects to assess the significance of recycled stellar ejecta as fuel for star formation. The fractional contributions of stellar mass loss to the cosmic star formation rate (SFR) and stellar mass densities increase with time, reaching $35 \%$ and $19 \%$, respectively, at $z=0$. The importance of recycling increases steeply with galaxy stellar mass for $M_{\ast} < 10^{10.5}$ M$_{\odot}$, and decreases mildly at higher mass. This trend arises from the mass dependence of feedback associated with star formation and AGN, which preferentially suppresses star formation fuelled by recycling. Recycling is more important for satellites than centrals and its contribution decreases with galactocentric radius. The relative contribution of AGB stars increases with time and towards galaxy centers. This is a consequence of the more gradual release of AGB ejecta compared to that of massive stars, and the preferential removal of the latter by star formation-driven outflows and by lock up in stellar remnants. Recycling-fuelled star formation exhibits a tight, positive correlation with galaxy metallicity, with a secondary dependence on the relative abundance of alpha elements (which are predominantly synthesized in massive stars), that is insensitive to the subgrid models for feedback. Hence, our conclusions are directly relevant for the origin of the mass-metallicity relation and metallicity gradients. Applying the relation between recycling and metallicity to the observed mass-metallicity relation yields our best estimate of the mass-dependent contribution of recycling. For centrals with a mass similar to that of the Milky Way, we infer the contributions of recycled stellar ejecta to the SFR and stellar mass to be $35 \%$ and $20 \%$, respectively.

preprint2015arXiv

The alignment and shape of dark matter, stellar, and hot gas distributions in the EAGLE and cosmo-OWLS simulations

We report the alignment and shape of dark matter, stellar, and hot gas distributions in the EAGLE and cosmo-OWLS simulations. The combination of these state-of-the-art hydro-cosmological simulations enables us to span four orders of magnitude in halo mass ($11 < log_{10}(M_{200}/ [h^{-1}M_\odot]) < 15$), a wide radial range ($-2.3 < log_{10}(r/[h^{-1}Mpc ]) < 1.3$) and redshifts $0 < z < 1$. The shape parameters of the dark matter, stellar and hot gas distributions follow qualitatively similar trends: they become more aspherical (and triaxial) with increasing halo mass, radius and redshift. We measure the misalignment of the baryonic components (hot gas and stars) of galaxies with their host halo as a function of halo mass, radius, redshift, and galaxy type (centrals vs satellites and early- vs late-type). Overall, galaxies align well with the local distribution of the total (mostly dark) matter. However, the stellar distributions on galactic scales exhibit a median misalignment of about 45-50 degrees with respect to their host haloes. This misalignment is reduced to 25-30 degrees in the most massive haloes ($13 < log_{10}(M_{200}/ [h^{-1}M_\odot ]) < 15$). Half of the disc galaxies in the EAGLE simulations have a misalignment angle with respect to their host haloes larger than 40 degrees. We present fitting functions and tabulated values for the probability distribution of galaxy-halo misalignment to enable a straightforward inclusion of our results into models of galaxy formations based on purely collisionless N-body simulations.

preprint2015arXiv

The APOSTLE simulations: solutions to the Local Group's cosmic puzzles

The Local Group of galaxies offer some of the most discriminating tests of models of cosmic structure formation. For example, observations of the Milky Way (MW) and Andromeda satellite populations appear to be in disagreement with N-body simulations of the "Lambda Cold Dark Matter" (ΛCDM) model: there are far fewer satellite galaxies than substructures in cold dark matter halos (the "missing satellites" problem); dwarf galaxies seem to avoid the most massive substructures (the "too-big-to-fail" problem); and the brightest satellites appear to orbit their host galaxies on a thin plane (the "planes of satellites" problem). Here we present results from APOSTLE (A Project Of Simulating The Local Environment), a suite of cosmological hydrodynamic simulations of twelve volumes selected to match the kinematics of the Local Group (LG) members. Applying the Eagle code to the LG environment, we find that our simulations match the observed abundance of LG galaxies, including the satellite galaxies of the MW and Andromeda. Due to changes to the structure of halos and the evolution in the LG environment, the simulations reproduce the observed relation between stellar mass and velocity dispersion of individual dwarf spheroidal galaxies without necessitating the formation of cores in their dark matter profiles. Satellite systems form with a range of spatial anisotropies, including one similar to that of the MW, confirming that such a configuration is not unexpected in ΛCDM. Finally, based on the observed velocity dispersion, size, and stellar mass, we provide new estimates of the maximum circular velocity for the halos of nine MW dwarf spheroidals.

preprint2015arXiv

The distribution of atomic hydrogen in EAGLE galaxies: morphologies, profiles, and HI holes

We compare the mass and internal distribution of atomic hydrogen (HI) in 2200 present-day central galaxies with M_star > 10^10 M_Sun from the 100 Mpc EAGLE Reference simulation to observational data. Atomic hydrogen fractions are corrected for self-shielding using a fitting formula from radiative transfer simulations and for the presence of molecular hydrogen using an empirical or a theoretical prescription from the literature. The resulting neutral hydrogen fractions, M_(HI+H2)/M_star, agree with observations to better than 0.1 dex for galaxies with M_star between 10^10 and 10^11 M_Sun. Our fiducial, empirical H2 model based on gas pressure results in galactic HI mass fractions, M_HI/M_star, that agree with observations from the GASS survey to better than 0.3 dex, but the alternative theoretical H2 formula leads to a negative offset in M_HI/M_star of up to 0.5 dex. Visual inspection reveals that most HI disks in simulated HI-rich galaxies are vertically disturbed, plausibly due to recent accretion events. Many galaxies (up to 80 per cent) contain spuriously large HI holes, which are likely formed as a consequence of the feedback implementation in EAGLE. The HI mass-size relation of all simulated galaxies is close to (but 16 per cent steeper than) observed, and when only galaxies without large holes in the HI disc are considered, the agreement becomes excellent (better than 0.1 dex). The presence of large HI holes also makes the radial HI surface density profiles somewhat too low in the centre, at Σ_HI > 1 M_Sun pc^-2 (by a factor of <~ 2 compared to data from the Bluedisk survey). In the outer region (Σ_HI < 1 M_Sun pc^-2), the simulated profiles agree quantitatively with observations. Scaled by HI size, the simulated profiles of HI-rich (M_HI > 10^9.8 M_Sun) and control galaxies (10^9.1 M_Sun > M_HI > 10^9.8 M_Sun) follow each other closely, as observed. (Abridged)

preprint2015arXiv

The distribution of neutral hydrogen around high-redshift galaxies and quasars in the EAGLE simulation

The observed high covering fractions of neutral hydrogen (HI) with column densities above $\sim 10^{17} \rm{cm}^{-2}$ around Lyman-Break Galaxies (LBGs) and bright quasars at redshifts z ~ 2-3 has been identified as a challenge for simulations of galaxy formation. We use the EAGLE cosmological, hydrodynamical simulation, which has been shown to reproduce a wide range of galaxy properties and for which the subgrid feedback was calibrated without considering gas properties, to study the distribution of HI around high-redshift galaxies. We predict the covering fractions of strong HI absorbers ($N_{\rm{HI}} \gtrsim 10^{17} \rm{cm}^{-2}$) inside haloes to increase rapidly with redshift but to depend only weakly on halo mass. For massive ($M_{200} \gtrsim 10^{12} {\rm M_{\odot}}$) halos the covering fraction profiles are nearly scale-invariant and we provide fitting functions that reproduce the simulation results. While efficient feedback is required to increase the HI covering fractions to the high observed values, the distribution of strong absorbers in and around halos of a fixed mass is insensitive to factor of two variations in the strength of the stellar feedback. In contrast, at fixed stellar mass the predicted HI distribution is highly sensitive to the feedback efficiency. The fiducial EAGLE simulation reproduces both the observed global column density distribution function of HI and the observed radial covering fraction profiles of strong HI absorbers around LBGs and bright quasars.

preprint2015arXiv

The EAGLE simulations of galaxy formation: calibration of subgrid physics and model variations

We present results from thirteen cosmological simulations that explore the parameter space of the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) simulation project. Four of the simulations follow the evolution of a periodic cube L = 50 cMpc on a side, and each employs a different subgrid model of the energetic feedback associated with star formation. The relevant parameters were adjusted so that the simulations each reproduce the observed galaxy stellar mass function at z = 0.1. Three of the simulations fail to form disc galaxies as extended as observed, and we show analytically that this is a consequence of numerical radiative losses that reduce the efficiency of stellar feedback in high-density gas. Such losses are greatly reduced in the fourth simulation - the EAGLE reference model - by injecting more energy in higher density gas. This model produces galaxies with the observed size distribution, and also reproduces many galaxy scaling relations. In the remaining nine simulations, a single parameter or process of the reference model was varied at a time. We find that the properties of galaxies with stellar mass <~ M* (the "knee" of the galaxy stellar mass function) are largely governed by feedback associated with star formation, while those of more massive galaxies are also controlled by feedback from accretion onto their central black holes. Both processes must be efficient in order to reproduce the observed galaxy population. In general, simulations that have been calibrated to reproduce the low-redshift galaxy stellar mass function will still not form realistic galaxies, but the additional requirement that galaxy sizes be acceptable leads to agreement with a large range of observables.

preprint2015arXiv

The EAGLE simulations of galaxy formation: the importance of the hydrodynamics scheme

We present results from a subset of simulations from the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) suite in which the formulation of the hydrodynamics scheme is varied. We compare simulations that use the same subgrid models without re-calibration of the parameters but employing the standard GADGET flavour of smoothed particle hydrodynamics (SPH) instead of the more recent state-of-the-art ANARCHY formulation of SPH that was used in the fiducial EAGLE runs. We find that the properties of most galaxies, including their masses and sizes, are not significantly affected by the details of the hydrodynamics solver. However, the star formation rates of the most massive objects are affected by the lack of phase mixing due to spurious surface tension in the simulation using standard SPH. This affects the efficiency with which AGN activity can quench star formation in these galaxies and it also leads to differences in the intragroup medium that affect the X-ray emission from these objects. The differences that can be attributed to the hydrodynamics solver are, however, likely to be less important at lower resolution. We also find that the use of a time step limiter is important for achieving the feedback efficiency required to match observations of the low-mass end of the galaxy stellar mass function.

preprint2015arXiv

The effect of baryons on the inner density profiles of rich clusters

We use the "Evolution and assembly of galaxies and their environments" (EAGLE) cosmological simulation to investigate the effect of baryons on the density profiles of rich galaxy clusters. We focus on EAGLE clusters with $M_{200}>10^{14}~M_\odot$ of which we have six examples. The central brightest cluster galaxies (BCGs) in the simulation have steep stellar density profiles, $ρ_*(r) \propto r^{-3}$. Stars dominate the mass density for $r<10~\rm{kpc}$, and, as a result, the total mass density profiles are steeper than the Navarro-Frenk-White (NFW) profile, in remarkable agreement with observations. The dark matter halo itself closely follows the NFW form at all resolved radii ($r\gtrsim3.0~\rm{kpc}$). The EAGLE BCGs have similar surface brightness and line-of-sight velocity dispersion profiles as the BCGs in the sample of Newman et al., which have the most detailed measurements currently available. After subtracting the contribution of the stars to the central density, Newman et al. infer significantly shallower slopes than the NFW value, in contradiction with the EAGLE results. We discuss possible reasons for this discrepancy, and conclude that an inconsistency between the kinematical model adopted by Newman et al. for their BCGs, which assumes isotropic stellar orbits, and the kinematical structure of the EAGLE BCGs, in which the orbital stellar anisotropy varies with radius and tends to be radially biased, could explain at least part of the discrepancy.

preprint2015arXiv

The offsets between galaxies and their dark matter in Lambda-CDM

We use the "Evolution and Assembly of GaLaxies and their Environments" ( EAGLE ) suite of hydrodynamical cosmological simulations to measure offsets between the centres of stellar and dark matter components of galaxies. We find that the vast majority (>95%) of the simulated galaxies display an offset smaller than the gravitational softening length of the simulations (Plummer-equivalent $ε= 700$ pc), both for field galaxies and satellites in clusters and groups. We also find no systematic trailing or leading of the dark matter along a galaxy's direction of motion. The offsets are consistent with being randomly drawn from a Maxwellian distribution with $σ\leq 196$ pc. Since astrophysical effects produce no feasible analogues for the $1.62^{+0.47}_{-0.49}$ kpc offset recently observed in Abell 3827, the observational result is in tension with the collisionless cold dark matter model assumed in our simulations.

preprint2014arXiv

A Medieval Multiverse: Mathematical Modelling of the 13th Century Universe of Robert Grosseteste

In his treatise on light, written in about 1225, Robert Grosseteste describes a cosmological model in which the Universe is created in a big-bang like explosion and subsequent condensation. He postulates that the fundamental coupling of light and matter gives rises to the material body of the entire cosmos. Expansion is arrested when matter reaches a minimum density and subsequent emission of light from the outer region leads to compression and rarefaction of the inner bodily mass so as to create nine celestial spheres, with an imperfect residual core. In this paper we reformulate the Latin description in terms of a modern mathematical model. The equations which describe the coupling of light and matter are solved numerically, subject to initial conditions and critical criteria consistent with the text. Formation of a universe with a non-infinite number of perfected spheres is extremely sensitive to the initial conditions, the intensity of the light and the transparency of these spheres. In this "medieval multiverse", only a small range of opacity and initial density profiles lead to a stable universe with nine perfected spheres. As in current cosmological thinking, the existence of Grosseteste's universe relies on a very special combination of fundamental parameters.

preprint2014arXiv

Bent by baryons: the low mass galaxy-halo relation

The relation between galaxies and dark matter halos is of vital importance for evaluating theoretical predictions of structure formation and galaxy formation physics. We show that the widely used method of abundance matching based on dark matter only simulations fails at the low mass end because two of its underlying assumptions are broken: only a small fraction of low mass (below 10^9.5 solar masses) halos host a visible galaxy, and halos grow at a lower rate due to the effect of baryons. In this regime, reliance on dark matter only simulations for abundance matching is neither accurate nor self-consistent. We find that the reported discrepancy between observational estimates of the halo masses of dwarf galaxies and the values predicted by abundance matching does not point to a failure of LCDM, but simply to a failure to account for baryonic effects. Our results also imply that the Local Group contains only a few hundred observable galaxies in contrast with the thousands of faint dwarfs that abundance matching would suggest. We show how relations derived from abundance matching can be corrected, so that they can be used self-consistently to calibrate models of galaxy formation.

preprint2014arXiv

Local Group galaxies emerge from the dark

The "Lambda Cold Dark Matter" (LCDM) model of cosmic structure formation is eminently falsifiable: once its parameters are fixed on large scales, it becomes testable in the nearby Universe. Observations within our Local Group of galaxies, including the satellite populations of the Milky Way and Andromeda, appear to contradict LCDM predictions: there are far fewer satellite galaxies than dark matter halos (the "missing satellites" problem), galaxies seem to avoid the largest substructures (the "too big to fail" problem), and the brightest satellites appear to orbit their host galaxies on a thin plane (the "planes of satellites" problem). We present results from the first hydrodynamic simulations of the Local Group that match the observed abundance of galaxies. We find that when baryonic and dark matter are followed simultaneously in the context of a realistic galaxy formation model, all three "problems" are resolved within the LCDM paradigm.

preprint2014arXiv

Overconsumption, outflows and the quenching of satellite galaxies

The baryon cycle of galaxies is a dynamic process involving the intake, consumption and ejection of vast quantities of gas. In contrast, the conventional picture of satellite galaxies has them methodically turning a large gas reservoir into stars until this reservoir is forcibly removed due to external ram pressure. This picture needs revision. Our modern understanding of the baryon cycle suggests that in some regimes the simple interruption of the fresh gas supply may quench satellite galaxies long before stripping events occur, a process we call overconsumption. We compile measurements from the literature of observed satellite quenching times at a range of redshifts to determine if satellites are principally quenched through orbit-based gas stripping events -- either direct stripping of the disk (ram pressure stripping) or the extended gas halo (strangulation) -- or from internally-driven star formation outflows via overconsumption. The observed timescales show significant deviation from the evolution expected for gas stripping mechanisms and suggest that either ram pressure stripping is much more efficient at high redshift, or that secular outflows quench satellites before orbit-based stripping occurs. Given the strong redshift evolution of star formation rates, at high redshift (z > 1.5) even moderate outflow rates will lead to extremely short quenching times with the expectation that such satellites will be quenched almost immediately following the cessation of cosmological inflow, regardless of stripping events. Observations of high redshift satellites give an indirect but sensitive measure of the outflow rate with current measurements suggesting that outflows are no larger than 2.5 times the star formation rate for galaxies with a stellar mass of 10^{10.5} solar masses.

preprint2014arXiv

The chosen few: the low mass halos that host faint galaxies

Since reionization prevents star formation in most halos below 3 x 10^9 solar masses, dwarf galaxies only populate a fraction of existing dark matter halos. We use hydrodynamic cosmological simulations of the Local Group to study the discriminating factors for galaxy formation in the early Universe and connect them to the present-day properties of galaxies and halos. A combination of selection effects related to reionization, and the subsequent evolution of halos in different environments, introduces strong biases between the population of halos that host dwarf galaxies, and the total halo population. Halos that host galaxies formed earlier and are more concentrated. In addition, halos more affected by tidal stripping are more likely to host a galaxy for a given mass or maximum circular velocity, vmax, today. Consequently, satellite halos are populated more frequently than field halos, and satellite halos of 10^8 - 10^9 solar masses or vmax of 12 - 20 km/s, similar to the Local Group dwarf spheroidals, have experienced a greater than average reduction in both mass and vmax after infall. They are on closer, more radial orbits with higher infall velocities and earlier infall times. Together, these effects make dwarf galaxies highly biased tracers of the underlying dark matter distribution.

preprint2014arXiv

The EAGLE project: Simulating the evolution and assembly of galaxies and their environments

We introduce the Virgo Consortium's EAGLE project, a suite of hydrodynamical simulations that follow the formation of galaxies and black holes in representative volumes. We discuss the limitations of such simulations in light of their finite resolution and poorly constrained subgrid physics, and how these affect their predictive power. One major improvement is our treatment of feedback from massive stars and AGN in which thermal energy is injected into the gas without the need to turn off cooling or hydrodynamical forces, allowing winds to develop without predetermined speed or mass loading factors. Because the feedback efficiencies cannot be predicted from first principles, we calibrate them to the z~0 galaxy stellar mass function and the amplitude of the galaxy-central black hole mass relation, also taking galaxy sizes into account. The observed galaxy mass function is reproduced to $\lesssim 0.2$ dex over the full mass range, $10^8 < M_*/M_\odot \lesssim 10^{11}$, a level of agreement close to that attained by semi-analytic models, and unprecedented for hydrodynamical simulations. We compare our results to a representative set of low-redshift observables not considered in the calibration, and find good agreement with the observed galaxy specific star formation rates, passive fractions, Tully-Fisher relation, total stellar luminosities of galaxy clusters, and column density distributions of intergalactic CIV and OVI. While the mass-metallicity relations for gas and stars are consistent with observations for $M_* \gtrsim 10^9 M_\odot$, they are insufficiently steep at lower masses. The gas fractions and temperatures are too high for clusters of galaxies, but for groups these discrepancies can be resolved by adopting a higher heating temperature in the subgrid prescription for AGN feedback. EAGLE constitutes a valuable new resource for studies of galaxy formation.

preprint2014arXiv

The metallicity of galactic winds

The abundance evolution of galaxies depends critically on the balance between the mixing of metals in their interstellar medium, the inflow of new gas and the outflow of enriched gas. We study these processes in gas columns perpendicular to a galactic disk using sub-parsec resolution simulations that track stellar ejecta with the Flash code. We model a simplified interstellar medium stirred and enriched by supernovae and their progenitors. We vary the density distribution of the gas column and integrate our results over an exponential disk to predict wind and ISM enrichment properties for disk galaxies. We find that winds from more massive galaxies are hotter and more highly enriched, in stark contrast to that which is often assumed in galaxy formation models. We use these findings in a simple model of galactic enrichment evolution, in which the metallicity of forming galaxies is the result of accretion of nearly pristine gas and outflow of enriched gas along an equilibrium sequence. We compare these predictions to the observed mass-metallicity relation, and demonstrate how the galaxy's gas fraction is a key controlling parameter. This explains the observed flattening of the mass-metallicity relation at higher stellar masses.

preprint2013arXiv

Do group dynamics play a role in the evolution of member galaxies?

We examine galaxy groups from the present epoch to z = 1 to explore the impact of group dynamics on galaxy evolution. We use group catalagues from the Sloan Digital Sky Survey (SDSS), the Group Environment and Evolution Collaboration (GEEC) and the high redshift GEEC2 sample to study how the observed member properties depend on galaxy stellar mass, group dynamical mass and dynamical state of the host group. We find a strong correlation between the fraction of non-star-forming (quiescent) galaxies and galaxy stellar mass, but do not detect a significant difference in the quiescent fraction with group dynamical mass, within our sample halo mass range of 10^13-10^14.5 M_sun, or with dynamical sate. However, at a redshift of approximately 0.4 we do see some evidence that the quiescent fraction in low mass galaxies (log(M_star/M_sun) < 10.5) is lower in groups with substructure. Additionally, our results show that the fraction of groups with non-Gaussian velocity distributions increases with redshift to roughly z = 0.4, while the amount of detected substructure remains constant to z = 1. Based on these results, we conclude that for massive galaxies (log(M_star/M_sun_ > 10.5), evolution is most strongly correlated to the stellar mass of a galaxy with little or no additional effect related to either the group dynamical mass or dynamical state. For low mass galaxies, we do see some evidence of a correlation between the quiescent fraction and the amount of detected substructure, highlighting the need to probe further down the stellar mass function to elucidate the role of the environment in galaxy evolution.

preprint2013arXiv

Efficient satellite quenching at z~1 from the GEEC2 spectroscopic survey of galaxy groups

We present deep GMOS-S spectroscopy for 11 galaxy groups at 0.8<z<1.0, for galaxies with r_{AB}<24.75. Our sample is highly complete (>66%) for eight of the eleven groups. Using an optical-NIR colour-colour diagram, the galaxies in the sample were separated with a dust insensitive method into three categories: passive (red), star-forming (blue), and intermediate (green). The strongest environmental dependence is observed in the fraction of passive galaxies, which make up only ~20 per cent of the field in the mass range 10^{10.3}<M_{star}/M_\odot<10^{11.0} but are the dominant component of groups. If we assume that the properties of the field are similar to those of the `pre-accreted' population, the environment quenching efficiency (ε_ρ) is defined as the fraction of field galaxies required to be quenched in order to match the observed red fraction inside groups. The efficiency obtained is ~0.4, similar to its value in intermediate-density environments locally. While green (intermediate) galaxies represent ~20 per cent of the star-forming population in both the group and field, at all stellar masses, the average sSFR of the group population is lower by a factor of ~3. The green population does not show strong H-delta absorption that is characteristic of starburst galaxies. Finally, the high fraction of passive galaxies in groups, when combined with satellite accretion models, require that most accreted galaxies have been affected by their environment. Thus, any delay between accretion and the onset of truncation of star formation (τ) must be <2 Gyr, shorter than the 3-7 Gyr required to fit data at z=0. The relatively small fraction of intermediate galaxies requires that the actual quenching process occurs quickly, with an exponential decay timescale of τ_q<1 Gyr.

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.

preprint2013arXiv

Star Formation and Environmental Quenching of GEEC2 Group Galaxies at z~1

We present new analysis from the GEEC2 spectroscopic survey of galaxy groups at $0.8<z<1$. Our previous work revealed an intermediate population between the star-forming and quiescent sequences and a strong environmental dependence in the fraction of quiescent galaxies. Only $\sim5$ per cent of star-forming galaxies in both the group and field sample show a significant enhancement in star formation, which suggests that quenching is the primary process in the transition from the star-forming to the quiescent state. To model the environmental quenching scenario, we have tested the use of different exponential quenching timescales and delays between satellite accretion and the onset of quenching. We find that with no delay, the quenching timescale needs to be long in order to match the observed quiescent fraction, but then this model produces too many intermediate galaxies. Fixing a delay time of 3 Gyr, as suggested from the local universe, produces too few quiescent galaxies. The observed fractions are best matched with a model that includes a delay that is proportional to the dynamical time and a rapid quenching timescale ($\sim0.25$ Gyr), but this model also predicts intermediate galaxies Hδ strength higher than that observed. Using stellar synthesis models, we have tested other scenarios, such as the rejuvenation of star formation in early-type galaxies and a portion of quenched galaxies possessing residual star formation. If environment quenching plays a role in the GEEC2 sample, then our work suggests that only a fraction of intermediate galaxies may be undergoing this transition and that quenching occurs quite rapidly in satellite galaxies ($\lesssim0.25$ Gyr).

preprint2012arXiv

Exploring the Diversity of Groups at 0.1<z<0.8 with X-ray and Optically Selected Samples

We present the global group properties of two samples of galaxy groups containing 39 high quality X-ray selected systems and 38 optically (spectroscopically) selected systems in coincident spatial regions at 0.12<z<0.79. Only nine optical systems are associable with X-ray systems. We discuss the confusion inherent in the matching of both galaxies to extended X-ray emission and of X-ray emission to already identified optical systems. Extensive spectroscopy has been obtained and the resultant redshift catalog and group membership are provided here. X-ray, dynamical, and total stellar masses of the groups are also derived and presented. We explore the effects of applying three different kinds of radial cut to our systems: a constant cut of 1 Mpc and two r200 cuts, one based on the velocity dispersion of the system and the other on the X-ray emission. We find that an X-ray based r200 results in less scatter in scaling relations and less dynamical complexity as evidenced by results of the Anderson-Darling and Dressler-Schectman tests, indicating that this radius tends to isolate the virialized part of the system. The constant and velocity dispersion based cuts can overestimate membership and can work to inflate velocity dispersion and dynamical and stellar mass. We find Lx-sigma and Mstellar-Lx scaling relations for X-ray and optically selected systems are not dissimilar. The mean fraction of mass found in stars for our systems is approximately 0.014 with a logarithmic standard deviation of 0.398 dex. We also define and investigate a sample of groups which are X-ray underluminous given the total group stellar mass. For these systems the fraction of stellar mass contributed by the most massive galaxy is typically lower than that found for the total population of groups implying that there may be less IGM contributed from the most massive member in these systems. (Abridged)

preprint2012arXiv

How supernova explosions power galactic winds

Feedback from supernovae is an essential aspect of galaxy formation. In order to improve subgrid models of feedback we perform a series of numerical experiments to investigate how supernova explosions power galactic winds. We use the Flash hydrodynamic code to model a simplified ISM, including gravity, hydrodynamics, radiative cooling above 10,000 K, and star formation that reproduces the Kennicutt-Schmidt relation. By simulating a small patch of the ISM in a tall box perpendicular to the disk, we obtain sub-parsec resolution allowing us to resolve individual supernova events and we investigate how the wind properties depend on those of the ISM and the galaxy. We find that outflows are more efficient in disks with lower surface densities or gas fractions. A simple model in which the warm cloudy medium is the barrier that limits the expansion of blast waves reproduces the scaling of outflow properties with disk parameters at high star formation rates. The scaling we find sets the investigation of galaxy winds on a new footing, providing a physically motivated sub-grid description of winds that can be implemented in cosmological hydrodynamic simulations and phenomenological models. [Abridged]

preprint2012arXiv

Substructure in the Most Massive GEEC Groups: Field-like Populations in Dynamically Active Groups

The presence of substructure in galaxy groups and clusters is believed to be a sign of recent galaxy accretion and can be used not only to probe the assembly history of these structures, but also the evolution of their member galaxies. Using the Dressler-Shectman (DS) Test, we study substructure in a sample of intermediate redshift (z ~ 0.4) galaxy groups from the Group Environment and Evolution Collaboration (GEEC) group catalog. We find that 4 of the 15 rich GEEC groups, with an average velocity dispersion of ~525 km s-1, are identified as having significant substructure. The identified regions of localized substructure lie on the group outskirts and in some cases appear to be infalling. In a comparison of galaxy properties for the members of groups with and without substructure, we find that the groups with substructure have a significantly higher fraction of blue and star-forming galaxies and a parent colour distribution that resembles that of the field population rather than the overall group population. In addition, we observe correlations between the detection of substructure and other dynamical measures, such as velocity distributions and velocity dispersion profiles. Based on this analysis, we conclude that some galaxy groups contain significant substructure and that these groups have properties and galaxy populations that differ from groups with no detected substructure. These results indicate that the substructure galaxies, which lie preferentially on the group outskirts and could be infalling, do not exhibit signs of environmental effects, since little or no star-formation quenching is observed in these systems.

preprint2012arXiv

The slowly evolving role of environment in a spectroscopic survey of star formation in Mstar > 5E8 Msun galaxies since z=1

We present a deep [OII] emission line survey of faint galaxies (22.5<KAB<24) in the Chandra Deep Field South and the FIRES field. With these data we measure the star formation rate (SFR) in galaxies in the stellar mass range 8.85 < log(M*/Msun) < 9.5 at 0.62<z<0.885, to a limit of SFR = 0.1Msun/yr. The presence of a massive cluster (MS1054-03) in the FIRES field, and of significant large scale structure in the CDFS field, allows us to study the environmental dependence of SFRs amongst this population of low-mass galaxies. Comparing our results with more massive galaxies at this epoch, with our previous survey (ROLES) at the higher redshift z=1, and with SDSS Stripe 82 data, we find no significant evolution of the stellar mass function of star-forming galaxies between z=0 and z=1, and no evidence that its shape depends on environment. The correlation between specific star formation rate (sSFR) and stellar mass at z=0.75 has a power-law slope of beta=-0.2, with evidence for a steeper relation at the lowest masses. The normalization of this correlation lies as expected between that corresponding to z=1 and the present day. The global SFR density is consistent with an evolution of the form (1+z)^2 over 0<z<1, with no evidence for a dependence on stellar mass. The sSFR of these star-forming galaxies at z=0.75 does not depend upon the density of their local environment. Considering just high-density environments, the low-mass end of the sSFR-M* relation in our data is steeper than that in Stripe 82 at z=0, and shallower than that measured by ROLES at z=1. Evolution of low-mass galaxies in dense environments appears to be more rapid than in the general field.

preprint2011arXiv

4098 galaxy clusters to z~0.6 in the Sloan Digital Sky Survey equatorial Stripe 82

We present a catalogue of 4098 photometrically selected galaxy clusters with a median redshift <z> = 0.32 in the 270 square degree 'Stripe 82' region of the Sloan Digital Sky Survey (SDSS), covering the celestial equator in the Southern Galactic Cap (-50 < RA < 59 deg, |Dec| < 1.25 deg). Owing to the multi-epoch SDSS coverage of this region, the ugriz photometry is ~2 magnitudes deeper than single scans within the main SDSS footprint. We exploit this to detect clusters of galaxies using an algorithm that searches for statistically significant overdensities of galaxies in a Voronoi tessellation of the projected sky. 32% of the clusters have at least one member with a spectroscopic redshift from existing public data (SDSS Data Release 7, 2SLAQ & WiggleZ), and the remainder have a robust photometric redshift (accurate to ~5-9% at the median redshift of the sample). The weighted average of the member galaxies' redshifts provides a reasonably accurate estimate of the cluster redshift. The cluster catalogue is publicly available for exploitation by the community to pursue a range of science objectives. In addition to the cluster catalogue, we provide a linked catalogue of 18,295 V<21 mag quasar sight-lines with impact parameters within <3 Mpc of the cluster cores selected from the catalogue of Veron et al. (2010). The background quasars cover 0.25 < z < 2, where MgII absorption-line systems associated with the clusters are detectable in optical spectra.

preprint2011arXiv

A spectroscopic measurement of galaxy formation timescales with ROLES

We present measurements of the specific star-formation rate (SSFR)-stellar mass relation for star-forming galaxies. Our deep spectroscopic samples are based on the Redshift One LDSS3 Emission line Survey, ROLES, and European Southern Observatory, ESO, public spectroscopy at z=1, and on the Sloan Digital Sky Survey (SDSS) at z=0.1. These datasets cover an equally deep mass range of 8.5<~log(M*/Msun)<~11 at both epochs. We find that the SSFR--mass relation evolves in a way which is remarkably independent of stellar mass, as we previously found for the star-formation rate density (SFRD)--mass relation. At higher masses, such as those probed by previous surveys, the evolution in SSFR--mass is almost independent of stellar mass. At higher masses (log(M*/Msun)>10) the shapes of the cumulative cosmic SFRDs are very similar at both z=0.1 and z=1.0, both showing 70% of the total SFRD above a mass of log(M*/Msun)>10. Mass functions are constructed for star-forming galaxies and found to evolve by only <35% between z=1 and z=0.1 over the whole mass range. The evolution is such that the mass function decreases with increasing cosmic time, confirming that galaxies are leaving the star-forming sequence/blue cloud. The observational results are extended to z~2 by adding two recent Lyman break galaxy samples, and data at these three epochs (z=0.1, 1, 2) are compared with the GALFORM semi-analytic model of galaxy formation. GALFORM predicts an overall SFR density (SFRD) as a function of stellar mass in reasonable agreement with the observations. The star formation timescales inferred from 1/SSFR also give reasonable overall agreement, with the agreement becoming worse at the lowest and highest masses. [abridged]

preprint2011arXiv

Direct observational evidence for a large transient galaxy population in groups at 0.85<z<1

(abridged) We introduce our survey of galaxy groups at 0.85<z<1, as an extension of the Group Environment and Evolution Collaboration (GEEC). Here we present the first results, based on Gemini GMOS-S nod-and-shuffle spectroscopy of seven galaxy groups selected from spectroscopically confirmed, extended XMM detections in COSMOS. In total we have over 100 confirmed group members, and four of the groups have >15 members. The dynamical mass estimates are in good agreement with the masses estimated from the X-ray luminosity, with most of the groups having 13<log(Mdyn/Msun)<14. Our spectroscopic sample is statistically complete for all galaxies with Mstar>1E10.1 Msun, and for blue galaxies we sample masses as low as Mstar=1E8.8 Msun. Like lower-redshift groups, these systems are dominated by red galaxies, at all stellar masses Mstar>1E10.1 Msun. Few group galaxies inhabit the "blue cloud" that dominates the surrounding field; instead, we find a large and possibly distinct population of galaxies with intermediate colours. The "green valley" that exists at low redshift is instead well-populated in these groups, containing ~30 per cent of galaxies. These do not appear to be exceptionally dusty galaxies, and about half show prominent Balmer-absorption lines. Furthermore, their HST morphologies appear to be intermediate between those of red-sequence and blue-cloud galaxies of the same stellar mass. We postulate that these are a transient population, migrating from the blue cloud to the red sequence, with a star formation rate that declines with an exponential timescale 0.6 Gyr< tau < 2 Gyr. Their prominence among the group galaxy population, and the marked lack of blue, star-forming galaxies, provides evidence that the group environment either directly reduces star formation in member galaxies, or at least prevents its rejuvenation during the normal cycle of galaxy evolution.

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.

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.

preprint2011arXiv

The local star-formation rate density: assessing calibrations using [OII], Ha and UV luminosities

We explore the use of simple star-formation rate (SFR) indicators (such as may be used in high-redshift galaxy surveys) in the local Universe using [OII], Ha, and u-band luminosities from the deeper 275 deg^2 Stripe 82 subsample of the Sloan Digital Sky Survey (SDSS) coupled with UV data from the Galaxy Evolution EXplorer satellite (GALEX). We examine the consistency of such methods using the star-formation rate density (SFRD) as a function of stellar mass in this local volume, and quantify the accuracy of corrections for dust and metallicity on the various indicators. Rest-frame u-band promises to be a particularly good SFR estimator for high redshift studies since it does not require a particularly large or sensitive extinction correction, yet yields results broadly consistent with more observationally expensive methods. We suggest that the [OII]-derived SFR, commonly used at higher redshifts (z~1), can be used to reliably estimate SFRs for ensembles of galaxies, but for high mass galaxies (log(M*/Msun)>10), a larger correction than is typically used is required to compensate for the effects of metallicity dependence and dust extinction. We provide a new empirical mass-dependent correction for the [OII]-SFR.

preprint2011arXiv

The population of Milky Way satellites in the LambdaCDM cosmology

We present a model for the satellites of the Milky Way in which galaxy formation is followed using semi-analytic techniques applied to the six high-resolution N-body simulations of galactic halos of the Aquarius project. The model, calculated using the Galform code, incorporates improved treatments of the relevant physics in the LambdaCDM cosmogony, particularly a self-consistent calculation of reionization by UV photons emitted by the forming galaxy population, including the progenitors of the central galaxy. Along the merger tree of each halo, the model calculates gas cooling (by Compton scattering off cosmic microwave background photons, molecular hydrogen and atomic processes), gas heating (from hydrogen photoionization and supernova energy), star formation and evolution. The evolution of the intergalactic medium is followed simultaneously with that of the galaxies. Star formation in the more massive progenitor subhalos is suppressed primarily by supernova feedback, while for smaller subhalos it is suppressed primarily by photoionization due to external and internal sources. The model is constrained to match a wide range of properties of the present day galaxy population as a whole, but at high redshift it requires an escape fraction of UV photons near unity in order completely to reionize the universe by redshift z ~ 8. In the most successful model the local sources photoionize the pre-galactic region completely by z ~ 10. In addition to the luminosity function of Milky Way satellites, the model matches their observed luminosity-metallicity relation, their radial distribution and the inferred values of the mass within 300 pc, which in the models increase slowly but significantly with luminosity. There is a large variation in satellite properties from halo to halo, with the luminosity function, for example, varying by a factor of ~ 2 among the six simulations.

preprint2010arXiv

Galaxy Formation Spanning Cosmic History

Over the past several decades, galaxy formation theory has met with significant successes. In order to test current theories thoroughly we require predictions for as yet unprobed regimes. To this end, we describe a new implementation of the Galform semi-analytic model of galaxy formation. Our motivation is the success of the model described by Bower et al. in explaining many aspects of galaxy formation. Despite this success, the Bower et al. model fails to match some observational constraints and certain aspects of its physical implementation are not as realistic as we would like. The model described in this work includes substantially updated physics, taking into account developments in our understanding over the past decade, and removes certain limiting assumptions made by this (and most other) semi-analytic models. This allows it to be exploited reliably in high-redshift and low mass regimes. Furthermore, we have performed an exhaustive search of model parameter space to find a particular set of model parameters which produce results in good agreement with a wide range of observational data (luminosity functions, galaxy sizes and dynamics, clustering, colours, metal content) over a wide range of redshifts. This model represents a solid basis on which to perform calculations of galaxy formation in as yet unprobed regimes.

preprint2010arXiv

Gas expulsion by quasar-driven winds as a solution to the over-cooling problem in galaxy groups and clusters

Galaxy groups are not scaled down versions of massive galaxy clusters - the hot gas in groups (known as the intragroup medium, IGrM for short) is, on average, less dense than the intracluster medium, implying that one or more non-gravitational processes (e.g., radiative cooling, star formation, and/or feedback) has had a relatively larger effect on groups. In the present study, we compare a number of cosmological hydrodynamic simulations that form part of the OverWhelmingly Large Simulations project to isolate and quantify the effects of cooling and feedback from supernovae (SNe) and active galactic nuclei (AGN) on the gas. This is achieved by comparing Lagrangian thermal histories of the gas in the different runs, which were all started from identical initial conditions. While radiative cooling, star formation, and SN feedback are all necessary ingredients, only runs that also include AGN feedback are able to successfully reproduce the optical and X-ray properties of groups and low-mass clusters. We isolate how, when, and exactly what gas is heated by AGN. Interestingly, we find that the gas that constitutes the present-day IGrM is that which was not strongly heated by AGN. Instead, the low median density/high median entropy of the gas in present-day groups is achieved by the ejection of lower entropy gas from low-mass progenitor galaxies at high redshift (primarily 2 < z < 4). This corresponds to the epoch when supermassive black holes accreted most of their mass, typically at a rate that is close to the Eddington limit (i.e., when the black holes are in a `quasar mode').

preprint2010arXiv

Mass transport by buoyant bubbles in galaxy clusters

We investigate the effect of three important processes by which AGN-blown bubbles transport material: drift, wake transport and entrainment. The first of these, drift, occurs because a buoyant bubble pushes aside the adjacent material, giving rise to a net upward displacement of the fluid behind the bubble. For a spherical bubble, the mass of upwardly displaced material is roughly equal to half the mass displaced by the bubble, and should be ~ 10^{7-9} solar masses depending on the local ICM and bubble parameters. We show that in classical cool core clusters, the upward displacement by drift may be a key process in explaining the presence of filaments behind bubbles. A bubble also carries a parcel of material in a region at its rear, known as the wake. The mass of the wake is comparable to the drift mass and increases the average density of the bubble, trapping it closer to the cluster centre and reducing the amount of heating it can do during its ascent. Moreover, material dropping out of the wake will also contribute to the trailing filaments. Mass transport by the bubble wake can effectively prevent the build-up of cool material in the central galaxy, even if AGN heating does not balance ICM cooling. Finally, we consider entrainment, the process by which ambient material is incorporated into the bubble. Abridged

preprint2010arXiv

The case for AGN feedback in galaxy groups

[Abridged] The relatively recent insight that energy input from supermassive black holes (BHs) can have a substantial effect on the star formation rates (SFRs) of galaxies motivates us to examine its effects on the scale of galaxy groups. At present, groups contain most of the galaxies and a significant fraction of the overall baryon content of the universe. To explore the effects of BH feedback on groups, we analyse two high resolution cosmological hydro simulations from the OverWhelmingly Large Simulations project. While both include galactic winds driven by supernovae, only one includes feedback from BHs. We compare the properties of the simulated groups to a wide range of observational data, including hot gas radial profiles and gas mass fractions (fgas), luminosity-mass-temperature (L-M-T) scaling relations, K-band luminosity of the group and its central brightest galaxy (CBG), SFRs and ages of the CBG, and gas/stellar metallicities. Both runs yield entropy profiles similar to the data, while the run without AGN feedback yields highly peaked temperature profiles, in discord with the observations. Energy input from BHs significantly reduces fgas for groups with masses less than ~10^14 Msun, yielding fgas-T and L-T relations that are in agreement with the data. The run without AGN feedback suffers from the well known overcooling problem; the resulting K-band luminosities are much larger than observed. By contrast, the run that includes BH feedback yields K-band luminosities and CBG SFRs and ages in agreement with current estimates. Both runs yield very similar gas-phase metallicities that match X-ray data, but they predict very different stellar metallicities. Based on the above, galaxy groups provide a compelling case that BH feedback is a crucial ingredient in the formation of massive galaxies.

preprint2010arXiv

The Dawn of the Red: Star formation histories of group galaxies over the past 5 billion years

We examine the star formation properties of group and field galaxies in two surveys, the Sloan Digital Sky Survey (SDSS; at z ~ 0.08) and the Group Environment and Evolution Collaboration (GEEC; at z ~ 0.4). Using UV imaging from the GALEX space telescope, along with optical and, for GEEC, near infrared photometry, we compare the observed spectral energy distributions to large suites of stellar population synthesis models. This allows us to accurately determine star formation rates and stellar masses. We find that star forming galaxies of all environments undergo a systematic lowering of their star formation rate between z=0.4 and z=0.08 regardless of mass. Nonetheless, the fraction of passive galaxies is higher in groups than the field at both redshifts. Moreover, the difference between the group and field grows with time and is mass-dependent, in the sense the the difference is larger at low masses. However, the star formation properties of star forming galaxies, as measured by their average specific star formation rates, are consistent within the errors in the group and field environment at fixed redshift. The evolution of passive fraction in groups between z=0.4 and z=0 is consistent with a simple accretion model, in which galaxies are environmentally affected 3 Gyrs after falling into a ~ 10E13 Msun group. This long timescale appears to be inconsistent with the need to transform galaxies quickly enough to ensure that star forming galaxies appear similar in both the group and field, as observed.

preprint2010arXiv

The mass assembly of galaxy groups and the evolution of the magnitude gap

We investigate the assembly of groups and clusters of galaxies using the Millennium dark matter simulation and the associated gas simulations and semi-analytic catalogues of galaxies. In particular, in order to find an observable quantity that could be used to identify early-formed groups, we study the development of the difference in magnitude between their brightest galaxies to assess the use of magnitude gaps as possible indicators. We select galaxy groups and clusters at redshift z=1 with dark matter halo mass M(R200) > 1E13/h Msun, and trace their properties until the present time (z=0). We consider only the systems with X-ray luminosity L_X> 0.25E42/h^2 erg/s at z=0. While it is true that a large magnitude gap between the two brightest galaxies of a particular group often indicates that a large fraction of its mass was assembled at an early epoch, it is not a necessary condition. More than 90% of fossil groups defined on the basis of their magnitude gaps (at any epoch between 0<z<1) cease to be fossils within 4 Gyr, mostly because other massive galaxies are assembled within their cores, even though most of the mass in their haloes might have been assembled at early times. We show that, compared to the conventional definition of fossil galaxy groups based on the magnitude gap Delta m(12)> 2 (in the R-band, within 0.5R200 of the centre of the group), an alternative criterion Delta m(14)>2.5 (within the same radius) finds 50% more early-formed systems, and those that on average retain their fossil phase longer. However, the conventional criterion performs marginally better at finding early-formed groups at the high-mass end of groups. Nevertheless, both criteria fail to identify a majority of the early-formed systems.

preprint2010arXiv

The Redshift One LDSS-3 Emission line Survey (ROLES) II: Survey method and z~1 mass-dependent star-formation rate density

Motivated by suggestions of 'cosmic downsizing', in which the dominant contribution to the cosmic star formation rate density (SFRD) proceeds from higher to lower mass galaxies with increasing cosmic time, we describe the design and implementation of the Redshift One LDSS3 Emission line Survey (ROLES). ROLES is a K-selected (22.5 < K_AB < 24.0) survey for dwarf galaxies [8.5<log(M*/Msun)< 9.5] at 0.89 < z < 1.15 drawn from two extremely deep fields (GOODS-S and MS1054-FIRES). Using the [OII]3727 emission line, we obtain redshifts and star-formation rates (SFRs) for star-forming galaxies down to a limit of ~0.3 Msun/yr. We present the [OII] luminosity function measured in ROLES and find a faint end slope of alpha_faint ~ -1.5, similar to that measured at z~0.1 in the SDSS. By combining ROLES with higher mass surveys, we measure the SFRD as a function of stellar mass using [OII] (with and without various empirical corrections), and using SED-fitting to obtain the SFR from the rest-frame UV luminosity for galaxies with spectroscopic redshifts. Our best estimate of the corrected [OII]-SFRD and UV SFRD both independently show that the SFRD evolves equally for galaxies of all masses between z~1 and z~0.1. The exact evolution in normalisation depends on the indicator used, with the [OII]-based estimate showing a change of a factor of ~2.6 and the UV-based a factor of ~6. We discuss possible reasons for the discrepancy in normalisation between the indicators, but note that the magnitude of this uncertainty is comparable to the discrepancy between indicators seen in other z~1 works. Our result that the shape of the SFRD as a function of stellar mass (and hence the mass range of galaxies dominating the SFRD) does not evolve between z~1 and z~0.1 is robust to the choice of indicator. [abridged]

preprint2010arXiv

The stellar and hot gas content of low-mass galaxy clusters

We analyse the stellar and hot gas content of 18 nearby, low-mass galaxy clusters, detected in redshift space and selected to have a dynamical mass 3E14<M/Msun<6E14, as measured from the 2dF Galaxy Redshift Survey. We combine X-ray measurements from both Chandra and XMM with ground-based near-infrared observations from CTIO, AAT and CFHT to compare the mass in hot gas and stars to the dynamical mass and state of the clusters. Only 13 of the clusters are detected in X-ray emission, and for these systems we find that a range of 7-20 per cent of their baryonic mass, and <3 per cent of their dynamical mass, is detected in starlight, similar to what is observed in more massive clusters. In contrast, the five undetected clusters are underluminous in X-ray emission, by up to a factor 10, given their stellar mass. Although the velocity distribution of cluster members in these systems is indistinguishable from a Gaussian, all show subtle signs of being unrelaxed: either they lack a central, dominant galaxy, or the bright galaxy distribution is less concentrated and/or more elongated than the rest of the sample. Thus we conclude that low-mass clusters and groups selected from the velocity distribution of their galaxies exhibit a dichotomy in their hot gas properties. Either they are detected in X-ray, in which case they generally lie on the usual scaling relations, or they are completely undetected in X-ray emission. The non-detections may be partly related to the apparently young dynamical state of the clusters, but it remains a distinct possibility that some of these systems are exceptionally devoid of hot emitting gas as the result of its expulsion or rarefaction.

preprint2009arXiv

Hot gas halos around disk galaxies: Confronting cosmological simulations with observations

Models of disk galaxy formation commonly predict the existence of an extended reservoir of accreted hot gas surrounding massive spirals at low redshift. As a test of these models, we use X-ray and H-alpha data of the two massive, quiescent edge-on spirals NGC 5746 and NGC 5170 to investigate the amount and origin of any hot gas in their halos. Contrary to our earlier claim, the Chandra analysis of NGC 5746, employing more recent calibration data, does not reveal any significant evidence for diffuse X-ray emission outside the optical disk, with a 3-sigma upper limit to the halo X-ray luminosity of 4e39 erg/s. An identical study of the less massive NGC 5170 also fails to detect any extraplanar X-ray emission. By extracting hot halo properties of disk galaxies formed in cosmological hydrodynamical simulations, we compare these results to expectations for cosmological accretion of hot gas by spirals. For Milky Way-sized galaxies, these high-resolution simulations predict hot halo X-ray luminosities which are lower by a factor of ~2 compared to our earlier results reported by Toft et al. (2002). We find the new simulation predictions to be consistent with our observational constraints for both NGC 5746 and NGC 5170, while also confirming that the hot gas detected so far around more actively star-forming spirals is in general probably associated with stellar activity in the disk. Observational results on quiescent disk galaxies at the high-mass end are nevertheless providing powerful constraints on theoretical predictions, and hence on the assumed input physics in numerical studies of disk galaxy formation and evolution.

preprint2009arXiv

Large-scale Structure and Dynamics of the Most X-ray Luminous Galaxy Cluster Known -- RX J1347-1145

We present photometric, spectroscopic and weak lensing analysis of the large-scale structure and dynamics of the most X-ray luminous galaxy cluster known, RX J1347-1145, at z=0.451. We spectroscopically confirmed 47 new members with LDSS3 on the Magellan telescope. Together with previously known members, we measure a new velocity dispersion of 1163(+/-97) km/s. The mass inferred from our velocity dispersion is M200 = 1.16^{+0.32}_{-0.27}x10^{15} solar mass, with r200=1.85Mpc, under the assumption of a singular isothermal sphere. We also present a weak lensing analysis using deep CFHT data on this cluster, and find a deprojected mass of 1.47^{+0.46}_{-0.43}x10^{15} solar mass within r200, in excellent agreement with our dynamical estimate. Thus, our new dynamical mass estimate is consistent with that from weak lensing and X-ray studies in the literature, resolving a previously claimed discrepancy. We photometrically detect and spectroscopically confirm another massive cluster with sigma=780(+/-100) km/s and M200=3.4^{+1.4}_{-1.1}x10^{14} solar mass ~7Mpc south-west of RX J1347-1145, which we refer to as RXJ1347-SW. Our spectroscopic survey reveals a possible excess of galaxies in velocity space in the region between RX J1347-1145 and RXJ1347-SW; comparing with simulations, this excess appears consistent with that expected from a large filamentary structure traced by galaxies connecting these two clusters.

preprint2009arXiv

The accretion of galaxies into groups and clusters

We use the galaxy stellar mass and halo merger tree information from the semi-analytic model galaxy catalogue of Font et al. (2009) to examine the accretion of galaxies into a large sample of groups and clusters, covering a wide range in halo mass (10E12.9 to 10E15.3 Msun/h), and selected from each of four redshift epochs (z=0, 0.5, 1.0 and 1.5). We find that clusters at all examined redshifts have accreted a significant fraction of their final galaxy populations through galaxy groups. A 10E14.5 Msun/h mass cluster at z=0 has, on average, accreted ~ 40% of its galaxies (Mstellar > 10E9 Msun/h) from halos with masses greater than 10E13 Msun/h. Further, the galaxies which are accreted through groups are more massive, on average, than galaxies accreted through smaller halos or from the field population. We find that at a given epoch, the fraction of galaxies accreted from isolated environments is independent of the final cluster or group mass. In contrast, we find that observing a cluster of the same halo mass at each redshift epoch implies different accretion rates of isolated galaxies, from 5-6 % per Gyr at z=0 to 15% per Gyr at z=1.5. We find that combining the existence of a Butcher Oemler effect at z=0.5 and the observations that galaxies within groups display significant environmental effects with galaxy accretion histories justifies striking conclusions. Namely, that the dominant environmental process must begin to occur in halos of 10E12 -- 10E13 Msun/h, and act over timescales of > 2 Gyrs. This argues in favor of a mechanism like "strangulation", in which the hot halo of a galaxy is stripped upon infalling into a more massive halo . This simple model predicts that by z=1.5 galaxy groups and clusters will display little to no environmental effects.

preprint2009arXiv

The colour of galaxies in distant groups

(Abridged) We present new optical and near-infrared imaging for a sample of 98 spectroscopically-selected galaxy groups at 0.25<z<0.55. We measure accurate colours for group members and the surrounding field population, statistically complete above a stellar mass limit of M=1E10 Msun. The overall colour distribution is bimodal in both the field and group samples; but at fixed luminosity the fraction of group galaxies populating the red peak is larger, by 20+/-7 per cent, than that of the field. In particular, group members with early-type morphologies, as identified in HST imaging, exhibit a tight red sequence, similar to that seen for more massive clusters. We show that approximately 20-30 per cent of galaxies on the red sequence may be dust-reddened galaxies with non-negligible star formation and early-spiral morphologies. This is true of both the field and group sample, and shows little dependence on near infrared luminosity. Thus, the fraction of bright group members with no sign of star formation or AGN activity is 54+/-6 per cent. Our field sample, which includes galaxies in all environments, contains 35+/-3 per cent of such inactive galaxies, consistent with the amount expected if all such galaxies are located in groups and clusters. This reinforces our earlier conclusions, that dense environments at z<0.5 are associated with a premature cessation of star formation in some galaxies; in particular we find no evidence for significantly enhanced star formation in these environments. Simple galaxy formation models predict a quenching of star formation in groups that is too efficient, overpopulating the red sequence. Attempts to fix this by increasing the timescale of this quenching equally for all group members distorts the colour distribution in a way that is inconsistent with observations.

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.

preprint2000arXiv

A low global star formation rate in the rich galaxy cluster AC114 at z=0.32

We present the results of a wide-field survey for H-alpha emitting galaxies in the cluster AC114 at z=0.32. Spectra centred on H-alpha at the cluster redshift have been obtained for 586 galaxies to I~22 out to a radius around 2Mpc. At most, only ~10% of these were found to be H-alpha-emitting cluster members. These objects are predominantly blue and of late-type spiral morphology, consistent with them hosting star formation. However, ~65% of the cluster members classified morphologically as spirals (with HST), have no detectable H-alpha emission; star-formation and morphological evolution in cluster galaxies appear to be largely decoupled. Changes in the H-alpha detection rate and the strength of H-alpha emission with environment are found to be weak within the region studied. Star formation within the cluster members is also found to be strongly and uniformly suppressed, with the rates inferred from the H-alpha emission not exceeding 4Mo/yr, and AC114's H-alpha luminosity function being an order of magnitude below that observed for field galaxies at the same redshift. None of the galaxies detected have the high star formation rates associated with `starburst' galaxies; however, this may still be reconciliable with the known (8+/-3%) fraction of `post-starburst' galaxies within AC114, given the poorly determined but short lifetimes of starbursts and the possibility that much of the associated star formation is obscured by dust.

preprint1998arXiv

The Colour-Magnitude Relation as a Constraint on the Formation of Rich Cluster Galaxies

In this paper, we examine the role that the colour-magnitude relation (CMR) can play in constraining the formation history of rich cluster galaxies. Firstly, we consider the colour evolution of galaxies after star formation ceases. We show that the scatter of the CMR places a strong constraint on the spread in age of the bulk of the stellar population. However, although the bulk of stars must be formed in a short period, continuing formation of stars in a fraction of the galaxies is not so strongly constrained. We examine a model in which star formation occurs over an extended period of time in most galaxies. An extension of this type of star formation history allows us to reconcile the small present-day scatter of the CMR with the observed blue galaxy fractions of intermediate redshift galaxy clusters. Secondly, the CMR can also be used to constrain the degree of merging between pre-existing stellar systems. This test relies on the slope of the CMR. We show that random mergers between galaxies very rapidly remove any well-defined CMR. However, we prefer to examine the merger process using a self-consistent merger tree. In such a model there are two effects: massive galaxies preferentially merge with systems of similar mass; and the rate of mass growth is considerably smaller than for the random merger case. As a result of both of these effects, the CMR persists through a larger number of merger steps. The passive evolution of galaxy colours and their averaging in dissipationless mergers provide opposing constraints on the formation of cluster galaxies in a hierarchical model; but at the level of current constraints, a compromise solution appears possible.

preprint1997arXiv

A Dynamical Study of Optically Selected Distant Clusters

We present a programme of spectroscopic observations of galaxies in a sample of optically-selected clusters taken from the catalogue of Couch et al (1991). Previous ROSAT observations of these clusters have shown them to have lower X-ray luminosities, given their optical richness, than might be expected on the basis of local samples. In the present paper we extend this work by determining velocity dispersions of a subsample of the clusters. We confirm the dynamical reality of all but one of the original sample, and find velocity dispersions comparable with present-day clusters of equivalent comoving space density. Thus, in the context of the $L_X-σ$ relation for present-day clusters, there is evidence for a higher velocity dispersion at fixed X-ray luminosity. A key question is whether the high velocity dispersions are indicative of the gravitational potential. If they are, the X-ray luminosities measured in Bower et al., 1994 (Paper I), would then imply an implausibly low efficiency of X-ray generation. Alternatively, the discrepancy could be explained if the clusters were systems of lower virial temperature, in which the apparent velocity dispersion is inflated by an infalling, unrelaxed halo. This might result either from an increase with redshift in the infall rate for clusters, or from the preferential selection of clusters embedded in filaments oriented along the line of sight. Since clusters with similar properties can be found in local optically selected catalogues, we suggest that the latter explanation is more likely.

preprint1996arXiv

A Weak Lensing Survey in the Fields of z~1 Luminous Radio Sources

In this paper we present weak lensing observations of the fields around 8 z~1 luminous radio sources. These data are searched for the lensing signatures of clusters that are either physically associated with the radio objects, or are foreground systems projected along the line of sight. The radio sources were all imaged with WFPC-2/HST providing high quality shape information on large numbers of faint galaxies around them. Statistical analysis of the coherent shear field visible in the shapes of the faint galaxies indicates that we have detected a weak lensing signal close to one of the targets, 3C336 at z=0.927, with a high level of confidence. A second, independent WFPC-2 observation of this target reinforces this detection. Our results support the earlier suggestion of weak lensing in this field by Fort et al (1996) using ground-based data. We also combined the shear distributions in the remaining 7 field to improve our sensitivity to weak shear signals from any structure typically associated with these sources. We find no detectable signal and estimate an upper limit on the maximum shear allowed by our observations. Using an N(z) estimated from lensing analyses we convert our observed lensing signal and limits into estimates of the masses of the various structures. We suggest that further lensing observations of distant radio sources and their host environments may allow the cluster L_X-mass relationship to be mapped at high-z. This is crucial for interpreting the results of the next generation of deep X-ray surveys, and thus constraining the redshift evolution of the cluster mass function out to z=1.

Richard G. Bower

What is connected

Connect this record

See the researcher in context

Building this map preview

76 published item(s)

The importance of black hole repositioning for galaxy formation simulations

Galactic inflow and wind recycling rates in the EAGLE simulations

Galactic outflow rates in the EAGLE simulations

How feedback shapes galaxies: an analytic model

Hydrostatic mass estimates of massive galaxy clusters: a study with varying hydrodynamics flavours and non-thermal pressure support

A chronicle of galaxy mass assembly in the EAGLE simulation

A unified multi-wavelength model of galaxy formation

Being WISE II: Reducing the Influence of Star Formation History on the Mass-to-Light Ratio of Quiescent Galaxies

Bimodality of low-redshift circumgalactic O VI in non-equilibrium EAGLE zoom simulations

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

It's not easy being green: The evolution of galaxy colour in the EAGLE simulation

Music from the heavens - Gravitational waves from supermassive black hole mergers in the EAGLE simulations

Size matters: abundance matching, galaxy sizes, and the Tully-Fisher relation in EAGLE

Supermassive black holes in the EAGLE Universe. Revealing the observables of their growth

The brighter galaxies reionised the Universe

The EAGLE simulations of galaxy formation: public release of halo and galaxy catalogues

The EAGLE simulations: atomic hydrogen associated with galaxies

The effect of baryons on redshift space distortions and cosmic density and velocity fields in the EAGLE simulation

The environmental dependence of HI in galaxies in the EAGLE simulations

The fundamental plane of star formation in galaxies revealed by the EAGLE hydrodynamical simulations

The origin of compact galaxies with anomalously high black hole masses

The origin of the $α$-enhancement of massive galaxies

The origin of the enhanced metallicity of satellite galaxies

Baryon effects on the internal structure of LCDM halos in the EAGLE simulations

Colours and luminosities of z=0.1 simulated galaxies in the EAGLE simulations

Evidence for a change in the dominant satellite galaxy quenching mechanism at z=1

Intrinsic alignments of galaxies in the EAGLE and cosmo-OWLS simulations

Molecular hydrogen abundances of galaxies in the EAGLE simulations

Recycled stellar ejecta as fuel for star formation and implications for the origin of the galaxy mass-metallicity relation

The alignment and shape of dark matter, stellar, and hot gas distributions in the EAGLE and cosmo-OWLS simulations

The APOSTLE simulations: solutions to the Local Group's cosmic puzzles

The distribution of atomic hydrogen in EAGLE galaxies: morphologies, profiles, and HI holes

The distribution of neutral hydrogen around high-redshift galaxies and quasars in the EAGLE simulation

The EAGLE simulations of galaxy formation: calibration of subgrid physics and model variations

The EAGLE simulations of galaxy formation: the importance of the hydrodynamics scheme

The effect of baryons on the inner density profiles of rich clusters

The offsets between galaxies and their dark matter in Lambda-CDM

A Medieval Multiverse: Mathematical Modelling of the 13th Century Universe of Robert Grosseteste

Bent by baryons: the low mass galaxy-halo relation

Local Group galaxies emerge from the dark

Overconsumption, outflows and the quenching of satellite galaxies

The chosen few: the low mass halos that host faint galaxies

The EAGLE project: Simulating the evolution and assembly of galaxies and their environments

The metallicity of galactic winds

Do group dynamics play a role in the evolution of member galaxies?

Efficient satellite quenching at z~1 from the GEEC2 spectroscopic survey of galaxy groups

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

Star Formation and Environmental Quenching of GEEC2 Group Galaxies at z~1

Exploring the Diversity of Groups at 0.1<z<0.8 with X-ray and Optically Selected Samples

How supernova explosions power galactic winds

Substructure in the Most Massive GEEC Groups: Field-like Populations in Dynamically Active Groups

The slowly evolving role of environment in a spectroscopic survey of star formation in Mstar > 5E8 Msun galaxies since z=1

4098 galaxy clusters to z~0.6 in the Sloan Digital Sky Survey equatorial Stripe 82

A spectroscopic measurement of galaxy formation timescales with ROLES

Direct observational evidence for a large transient galaxy population in groups at 0.85<z<1

Numerical overcooling in shocks

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

The local star-formation rate density: assessing calibrations using [OII], Ha and UV luminosities

The population of Milky Way satellites in the LambdaCDM cosmology

Galaxy Formation Spanning Cosmic History

Gas expulsion by quasar-driven winds as a solution to the over-cooling problem in galaxy groups and clusters

Mass transport by buoyant bubbles in galaxy clusters

The case for AGN feedback in galaxy groups

The Dawn of the Red: Star formation histories of group galaxies over the past 5 billion years

The mass assembly of galaxy groups and the evolution of the magnitude gap

The Redshift One LDSS-3 Emission line Survey (ROLES) II: Survey method and z~1 mass-dependent star-formation rate density

The stellar and hot gas content of low-mass galaxy clusters

Hot gas halos around disk galaxies: Confronting cosmological simulations with observations

Large-scale Structure and Dynamics of the Most X-ray Luminous Galaxy Cluster Known -- RX J1347-1145

The accretion of galaxies into groups and clusters

The colour of galaxies in distant groups

The Colours of Satellite Galaxies in Groups and Clusters

A low global star formation rate in the rich galaxy cluster AC114 at z=0.32

The Colour-Magnitude Relation as a Constraint on the Formation of Rich Cluster Galaxies