The past decade has seen significant progress in understanding galaxy formation and evolution using large-scale cosmological simulations. While these simulations produce galaxies in overall good agreement with observations, they employ different sub-grid models for galaxies and supermassive black holes (BHs). We investigate the impact of the sub-grid models on the BH mass properties of the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA simulations, focusing on the M_BH-M_star relation and the BH mass function. All simulations predict tight M_BH-M_star relations, and struggle to produce the lowest (M_BH< 10^7.5 Msun) in galaxies of M_star~10^10.5-10^11.5 Msun. While the time evolution of the mean M_BH-M_star relation is mild (<1 dex in BH mass for 0<z<5) for all the simulations, its linearity (shape) and normalization varies from simulation to simulation. The strength of SN feedback has a large impact on the linearity and time evolution for M_star<10^10.5 Msun. We find that the low-mass end is a good discriminant of the simulation models, and highlights the need for new observational constraints. At the high-mass end, strong AGN feedback can suppress the time evolution of the relation normalization. Compared with the observations of the local universe, we find an excess of BHs with M_BH>10^9 Msun in most of the simulations. The BH mass function is dominated by efficiently accreting BHs (log10 f_Edd >-2$) at high redshifts, and transitions progressively from the high-mass to the low-mass end to be governed by inactive BHs. The transition time and the contribution of active BHs are different among the simulations, and can be used to evaluate models against observations.

We demonstrate that the UV brightest quasars at z=1-2 live in overdense environments. This is based on an analysis of deep Hubble Space Telescope WFC3 G141 grism spectroscopy of the galaxies along the lines-of-sight to UV luminous quasars in the redshift range z=1-2. This constitutes some of the deepest grism spectroscopy performed by WFC3, with 4 roll angles spread over a year of observations to mitigate the effect of overlapping spectra. Of the 12 quasar fields studied, 8 display evidence for a galaxy overdensity at the redshift of the quasar. One of the overdensities, PG0117+213 at z=1.50, has potentially 36 spectroscopically confirmed members, consisting of 19 with secure redshifts and 17 with single-line redshifts, within a cylinder of radius ~700 kpc. Its halo mass is estimated to be log (M/Msol)=14.7. This demonstrates that spectroscopic and narrow-band observations around distant UV bright quasars may be an excellent route for discovering protoclusters. Our findings agree with previous hints from statistical observations of the quasar population and theoretical works, as feedback regulated black hole growth predicts a correlation between quasar luminosity and halo mass. We also present the high signal-to-noise rest-frame optical spectral and photometric properties of the quasars themselves.

The abundance and distribution of metals in galaxy clusters contains valuable information about their chemical history and evolution. By looking at how metallicity evolves with redshift, it is possible to constrain the different metal production channels. We use the C-EAGLE clusters, a sample of 30 high resolution ($m_{gas} \simeq 1.8\times 10^{6}$ M$_{\odot}$) cluster zoom simulations, to investigate the redshift evolution of metallicity, with particular focus on the cluster outskirts. The early enrichment model, in which the majority of metals are produced in the core of cluster progenitors at high redshift, suggests that metals in cluster outskirts have not significantly evolved since $z=2$. With the C-EAGLE sample, we find reasonable agreement with the early enrichment model as there is very little scatter in the metallicity abundance at large radius across the whole sample, out to at least $z=2$. The exception is Fe for which the radial dependence of metallicity was found to evolve at low redshift as a result of being mainly produced by Type Ia supernovae, which are more likely to be formed at later times than core-collapse supernovae. We also found considerable redshift evolution of metal abundances in the cores of the C-EAGLE clusters which has not been seen in other simulations or observation based metallicity studies. Since we find this evolution to be driven by accretion of low metallicity gas, it suggests that the interaction between outflowing, AGN heated material and the surrounding gas is important for determining the core abundances in clusters.

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.

We address the issue of numerical convergence in cosmological smoothed particle hydrodynamics simulations using a suite of runs drawn from the EAGLE project. Our simulations adopt subgrid models that produce realistic galaxy populations at a fiducial mass and force resolution, but systematically vary the latter in order to study their impact on galaxy properties. We provide several analytic criteria that help guide the selection of gravitational softening for hydrodynamical simulations, and present results from runs that both adhere to and deviate from them. Unlike dark matter-only simulations, hydrodynamical simulations exhibit a strong sensitivity to gravitational softening, and care must be taken when selecting numerical parameters. Our results--which focus mainly on star formation histories, galaxy stellar mass functions and sizes--illuminate three main considerations. First, softening imposes a minimum resolved escape speed, $v_\epsilon$, due to the binding energy between gas particles. Runs that adopt such small softening lengths that $v_\epsilon \gt 10\,{\rm km s^{-1}}$ (the sound speed in ionised $\sim 10^4\,{\rm K}$ gas) suffer from reduced effects of photo-heating. Second, feedback from stars or active galactic nuclei may suffer from numerical over-cooling if the gravitational softening length is chosen below a critical value, $\epsilon_{\rm eFB}$. Third, we note that small softening lengths exacerbate the segregation of stars and dark matter particles in halo centres, often leading to the counter-intuitive result that galaxy sizes {\em increase} as softening is reduced. The structure of dark matter haloes in hydrodynamical runs respond to softening in a way that reflects the sensitivity of their galaxy populations to numerical parameters.

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.

We present an analysis of the gas dynamics of star-forming galaxies at z~1.5 using data from the KMOS Galaxy Evolution Survey (KGES). We quantify the morphology of the galaxies using $HST$ CANDELS imaging parametrically and non-parametrically. We combine the H$\alpha$ dynamics from KMOS with the high-resolution imaging to derive the relation between stellar mass (M$_{*}$) and stellar specific angular momentum (j$_{*}$). We show that high-redshift star-forming galaxies at z~1.5 follow a power-law trend in specific stellar angular momentum with stellar mass similar to that of local late-type galaxies of the form j$_*$$\propto$M$_*^{0.53 \pm 0.10}$. The highest specific angular momentum galaxies are mostly disc-like, although generally, both peculiar morphologies and disc-like systems are found across the sequence of specific angular momentum at a fixed stellar mass. We explore the scatter within the j$_{*}$-M$_{*}$ plane and its correlation with both the integrated dynamical properties of a galaxy (e.g. velocity dispersion, Toomre Q$_{\rm g}$, H$\alpha$ star formation rate surface density $\Sigma_{\rm SFR}$) and its parameterised rest-frame UV/optical morphology (e.g. S\'ersic index, bulge to total ratio, Clumpiness, Asymmetry and Concentration). We establish that the position in the j$_{*}$-M$_{*}$ plane is strongly correlated with the star-formation surface density and the Clumpiness of the stellar light distribution. Galaxies with peculiar rest-frame UV/optical morphologies have comparable specific angular momentum to disc-dominated galaxies of the same stellar mass, but are clumpier and have higher star-formation rate surface densities. We propose that the peculiar morphologies in high--redshift systems are driven by higher star formation rate surface densities and higher gas fractions leading to a more clumpy inter-stellar medium.

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.

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.

We study the impact of numerical parameters on the properties of cold dark matter haloes formed in collisionless cosmological simulations. We quantify convergence in the median spherically-averaged circular velocity profiles for haloes of widely varying particle number, as well as in the statistics of their structural scaling relations and mass functions. In agreement with prior work focused on single haloes, our results suggest that cosmological simulations yield robust halo properties for a wide range of gravitational softening parameters, $\epsilon$, provided: 1) $\epsilon$ is not larger than a "convergence radius", $r_{\rm conv}$, which is dictated by 2-body relaxation and determined by particle number, and 2) a sufficient number of timesteps are taken to accurately resolve particle orbits with short dynamical times. Provided these conditions are met, median circular velocity profiles converge to within $\approx 10$ per cent for radii beyond which the local 2-body relaxation timescale exceeds the Hubble time by a factor $\kappa\equiv t_{\rm relax}/t_{\rm H}\gt 0.177$, with better convergence attained for higher $\kappa$. We provide analytic estimates of $r_{\rm conv}$ that build on previous attempts in two ways: first, by highlighting its explicit (but weak) softening-dependence and, second, by providing a simpler criterion in which $r_{\rm conv}$ is determined entirely by the mean inter-particle spacing, $l$; for example, better than $10$ per cent convergence in circular velocity for $r\gt 0.05\,l$. We show how these analytic criteria can be used to assess convergence in structural scaling relations for dark matter haloes as a function of their mass or maximum circular speed.

We present an analysis of the gas outflow energetics of 529 main-sequence star-forming galaxies at z~1 using KMOS observations of the broad, underlying H-alpha and forbidden lines of [N II] and [S II]. Based on the stacked spectra for a sample with median star-formation rates and stellar masses of SFR ~ 7 Mo/yr and M* = (1.0+/-0.1)x10^10 Mo respectively, we derive a typical mass outflow rate of dM/dt = 1-4 Mo/yr and a mass loading of dM/dt/SFR = 0.2--0.4. The mass loading of the wind does not show a strong trend with star-formation rate over the range SFR ~ 2--20 Mo/yr, although we identify a trend with stellar mass such that dM/dt/SFR ~ M*^(0.26+/-0.07). Finally, we find that the line width of the broad H-alpha increases with disk circular velocity with a sub-linear scaling relation FWHM_broad ~ v^(0.21+/-0.05). As a result of this behavior, in the lowest mass galaxies (M* < 10^10 Mo), a significant fraction of the outflowing gas should have sufficient velocity to escape the gravitational potential of the halo whilst in the highest mass galaxies (M* > 10^10 Mo) most of the gas will be retained, flowing back on to the galaxy disk at later times.

We present recalibrations of the GALFORM semi-analytical model of galaxy formation in a new N-body simulation with the Planck cosmology. The Planck Millennium simulation uses more than 128 billion particles to resolve the matter distribution in a cube of $800$ Mpc on a side, which contains more than 77 million dark matter haloes with mass greater than $2.12 \times 10^{9} h^{-1} {\rm M_{\odot}}$ at the present day. Only minor changes to a very small number of model parameters are required in the recalibration. We present predictions for the atomic hydrogen content (HI) of dark matter halos, which is a key input into the calculation of the HI intensity mapping signal expected from the large-scale structure of the Universe. We find that the HI mass $-$ halo mass relation displays a clear break at the halo mass above which AGN heating suppresses gas cooling, $\approx 3 \times 10^{11} h^{-1} M_{\rm \odot}$. Below this halo mass, the HI content of haloes is dominated by the central galaxy; above this mass it is the combined HI content of satellites that prevails. We find that the HI mass - halo mass relation changes little with redshift up to $z=3$. The bias of HI sources shows a scale dependence that gets more pronounced with increasing redshift.

We present adaptive optics assisted integral field spectroscopy of 34 star-forming galaxies at $z$ = 0.8-3.3 selected from the HiZELS narrow-band survey. We measure the kinematics of the ionised interstellar medium on $\sim$1 kpc scales, and show that the galaxies are turbulent, with a median ratio of rotational to dispersion support of $v$/$\sigma$=0.82$\pm$0.13. We combine the dynamics with high-resolution rest-frame optical imaging and extract emission line rotation curves. We show that high-redshift star-forming galaxies follow a similar power-law trend in specific angular momentum with stellar mass as that of local late type galaxies. We exploit the high resolution of our data and examine the radial distribution of angular momentum within each galaxy by constructing total angular momentum profiles. Although the stellar mass of a typical star-forming galaxy is expected to grow by a factor $\sim$8 in the $\sim$5 Gyrs between $z$$\sim$3.3 and $z$$\sim$0.8, we show that the internal distribution of angular momentum becomes less centrally concentrated in this period i.e the angular momentum grows outwards. To interpret our observations, we exploit the EAGLE simulation and trace the angular momentum evolution of star forming galaxies from $z$$\sim$3 to $z$$\sim$0, identifying a similar trend of decreasing angular momentum concentration. This change is attributed to a combination of gas accretion in the outer disk, and feedback that preferentially arises from the central regions of the galaxy. We discuss how the combination of the growing bulge and angular momentum stabilises the disk and gives rise to the Hubble sequence.

We analyse 850um continuum observations of eight massive X-ray detected galaxy clusters at z~0.8-1.6 taken with SCUBA-2 on the James Clerk Maxwell Telescope. We find an average overdensity of 850um-selected sources of a factor of 4+/-2 per cluster within the central 1Mpc compared to the field. We investigate the multiwavelength properties of these sources and identify 34 infrared counterparts to 26 SCUBA-2 sources. Their colours suggest that the majority of these counterparts are probable cluster members. We use the multi-wavelength far-infrared photometry to measure the total luminosities and total cluster star-formation rates demonstrating that they are roughly three orders of magnitude higher than local clusters. We predict the H-band luminosities of the descendants of our cluster submillimetre galaxies and find that their stellar luminosity distribution is consistent with that of passive elliptical galaxies in z~0 clusters. Together, the faded descendants of the passive cluster population already in place at z~1 and the cluster submillimetre galaxies are able to account for the total luminosity function of early-type cluster galaxies at z~0. This suggests that the majority of the luminous passive population in z~0 clusters are likely to have formed at z>>1 through an extreme, dust-obscured starburst event.

About 35 years ago a class of galaxies with unusually strong Balmer absorption lines and weak emission lines was discovered in distant galaxy clusters. These objects, alternatively referred to as post-starburst, E+A or k+a galaxies, are now known to occur in all environments and at all redshifts, with many exhibiting compact morphologies and low-surface brightness features indicative of past galaxy mergers. They are commonly thought to represent galaxies that are transitioning from blue to red sequence, making them critical to our understanding of the origins of galaxy bimodality. However, recent observational studies have questioned this simple interpretation. From observations alone, it is challenging to disentangle the different mechanisms that lead to the quenching of star formation in galaxies. Here we present examples of three different evolutionary pathways that lead to galaxies with strong Balmer absorption lines in the EAGLE simulation: classical blue-to-red quenching, blue-to-blue cycle and red-to-red rejuvenation. The first two are found in both post-starburst galaxies and galaxies with truncated star formation. Each pathway is consistent with scenarios hypothesised for observational samples. The fact that 'post-starburst' signatures can be attained via various evolutionary channels explains the diversity of observed properties, and lends support to the idea that slower quenching channels are important at low redshift.

We present stellar specific angular momentum $j_*$ measurements of two $z\sim 1.5$ galaxies in the KGES sample and 12 DYNAMO $z\sim 0.1$ analogues of high-redshift galaxies. We combine natural seeing integral field spectroscopic data to trace line emission out to high multiples of effective radius $r_e$, with adaptive optics assisted Keck/OSIRIS observations to trace the rapid rise in rotation curve in the inner regions. Our spaxel-wise integration method gives results that are on average within measurement uncertainty of the traditional rotation curve model method. At $z\sim 0$, combining GMOS and OSIRIS datasets improves the measurement uncertainty in $j_*$ from 13\% (GMOS only) or 16\% (OSIRIS only) to 10\%. At $z\sim 1.5$, systematics allow for at best 20\% uncertainty on $j_*$. DYNAMO analogues of high-$z$ galaxies have low $j_*$ for their stellar mass $M_*$, and low bulge-to-total light ratio $\beta$ for their $j_*/M_*$. The high-$z$ galaxy COSMOS 127977 has $j_*/M_*$ consistent with normal local disk galaxies, while UDS 78317 is consistent with local analogues. However, our high-resolution OSIRIS data reveal that UDS 78317 may be a merging system. We report a relationship between distance to the $\beta-j_*/M_*$ plane and the ratio of velocity dispersion to rotational velocity $\sigma/v_{max}$, where galaxies that deviate more from the plane are more dispersion-dominated due to turbulence. Much of the scatter in $M_*-j_*$ that is not explained by variations in the bulge-to-total ratio or evolution with redshift may be driven by increased turbulence due to star formation, or by treating mergers as rotating disks.

We investigate the connection between the morphology and internal kinematics of the stellar component of central galaxies with mass $M_\star > {10}^{9.5} {\rm M}_\odot$ in the EAGLE simulations. We compare several kinematic diagnostics commonly used to describe simulated galaxies, and find good consistency between them. We model the structure of galaxies as ellipsoids and quantify their morphology via the ratios of their principal axes. We show that the differentiation of blue star-forming and red quiescent galaxies using morphological diagnostics can be achieved with similar efficacy to the use of kinematical diagnostics, but only if one is able to measure both the flattening and the triaxiality of the galaxy. Flattened oblate galaxies exhibit greater rotational support than their spheroidal counterparts, but there is significant scatter in the relationship between morphological and kinematical diagnostics, such that kinematically-similar galaxies can exhibit a broad range of morphologies. The scatter in the relationship between the flattening and the ratio of the rotation and dispersion velocities ($v/\sigma$) correlates strongly with the anisotropy of the stellar velocity dispersion: at fixed $v/\sigma$, flatter galaxies exhibit greater dispersion in the plane defined by the intermediate and major axes than along the minor axis, indicating that the morphology of simulated galaxies is influenced significantly by the structure of their velocity dispersion. The simulations reveal that this anisotropy correlates with the intrinsic morphology of the galaxy's inner dark matter halo, i.e. the halo's morphology that emerges in the absence of dissipative baryonic physics. This implies the existence of a causal relationship between the morphologies of galaxies and that of their host dark matter haloes.

We analyse maps of the spatially-resolved nebular emission of $\approx$1500 star-forming galaxies at $z\approx0.6$-$2.2$ from deep KMOS and MUSE observations to measure the average shape of their rotation curves. We use these to test claims for declining rotation curves at large radii in galaxies at $z\approx1$-$2$ that have been interpreted as evidence for an absence of dark matter. We show that the shape of the average rotation curves, and the extent to which they decline beyond their peak velocities, depends upon the normalisation prescription used to construct the average curve. Normalising in size by the galaxy stellar disk-scale length after accounting for seeing effects ($R_{\rm{d}}^{\prime}$), we construct stacked position-velocity diagrams that trace the average galaxy rotation curve out to $6R_{\rm{d}}^{\prime}$ ($\approx$13 kpc, on average). Combining these curves with average HI rotation curves for local systems, we investigate how the shapes of galaxy rotation curves evolve over $\approx$10 Gyr. The average rotation curve for galaxies binned in stellar mass, stellar surface mass density and/or redshift is approximately flat, or continues to rise, out to at least $6R_{\rm{d}}^{\prime}$. We find a trend between the outer slopes of galaxies' rotation curves and their stellar mass surface densities, with the higher surface density systems exhibiting flatter rotation curves. Drawing comparisons with hydrodynamical simulations, we show that the average shapes of the rotation curves for our sample of massive, star-forming galaxies at $z\approx0$-$2.2$ are consistent with those expected from $\Lambda$CDM theory and imply dark matter fractions within $6R_{\rm{d}}$ of at least $\approx60$ percent.

We investigate the disruption of group and cluster satellite galaxies with total mass (dark matter plus baryons) above 10^10 M_sun in the Hydrangea simulations, a suite of 24 high-resolution cosmological hydrodynamical zoom-in simulations based on the EAGLE model. The simulations predict that ~50 per cent of satellites survive to redshift z = 0, with higher survival fractions in massive clusters than in groups and only small differences between baryonic and pure N-body simulations. For clusters, up to 90 per cent of galaxy disruption occurs in lower-mass sub-groups (i.e., during pre-processing); 96 per cent of satellites in massive clusters that were accreted at z < 2 and have not been pre-processed survive. Of those satellites that are disrupted, only a few per cent merge with other satellites, even in low-mass groups. The survival fraction changes rapidly from less than 10 per cent of those accreted at high z to more than 90 per cent at low z. This shift, which reflects faster disruption of satellites accreted at higher z, happens at lower z for more massive galaxies and those accreted onto less massive haloes. The disruption of satellite galaxies is found to correlate only weakly with their pre-accretion baryon content, star formation rate, and size, so that surviving galaxies are nearly unbiased in these properties. These results suggest that satellite disruption in massive haloes is uncommon, and that it is predominantly the result of gravitational rather than baryonic processes.

We exploit EAGLE, a cosmological hydrodynamical simulation, to reproduce the selection of the observed sub-millimeter (submm) galaxy population by selecting the model galaxies at $z \geq 1$ with mock submm fluxes $S_{850} \geq 1$ mJy. There is a reasonable agreement between the galaxies within this sample and the properties of the observed submm population, such as their star formation rates (SFRs) at $z<3$, redshift distribution and many integrated galaxy properties. We find that the bulk of the $S_{850} \geq 1$ mJy model population is at $z = 2.5$, and that they are massive galaxies ($M_* \sim 10^{11}$ Msol) with high dust masses ($M_{\mathrm{dust}} \sim 10^{8}$ Msol), gas fractions ($f_{\mathrm{gas}} \approx 50$%) and SFRs ($\dot M_* \approx 100$ Msol/yr). They have major and minor merger fractions similar to the general population, suggesting that mergers are not the primary driver of the model submm galaxies. Instead, the $S_{850} \geq 1$ mJy model galaxies yield high SFRs primarily because they maintain a significant gas reservoir as a result of hosting an undermassive black hole. In addition, we find that not all highly star-forming EAGLE galaxies have submm fluxes $S_{850} > 1$ mJy. Thus, we investigate the nature of $z \geq 1$ highly star-forming Submm-Faint galaxies (i.e., $\dot M_* \geq 80$ Msol/yr but $S_{850}< 1$ mJy). We find they are similar to the model submm galaxies; being gas rich and hosting undermassive black holes, however they are typically lower mass ($M_* \sim 10^{10}$ Msol) and are at higher redshifts ($z>4$). These typically higher-$z$ galaxies show stronger evidence for having been triggered by major mergers, and critically, they are likely missed by current submm surveys due to their higher dust temperatures. This suggests a potentially even larger contribution to the SFR density at $z > 3$ from dust-obscured systems than implied by current observations.

The EAGLE simulation suite has previously been used to investigate the relationship between the stellar mass of galaxies, $M_{*}$, and the properties of dark matter halos, using the hydrodynamical reference simulation combined with a dark matter only (DMO) simulation having identical initial conditions. The stellar masses of central galaxies in halos with $M_{\mathrm{200c}} > 10^{11} \mathrm{M_{\odot}}$ were shown to correlate with the DMO halo maximum circular velocity, with $\approx 0.2$ dex of scatter that is uncorrelated with other DMO halo properties. Here we revisit the origin of the scatter in the $M_{*}-V_{\mathrm{max, DMO}}$ relation in EAGLE at $z = 0.1$. We find that the scatter in $M_{*}$ correlates with the mean age of the galaxy stellar population such that more massive galaxies at fixed $V_{\mathrm{max, DMO}}$ are younger. The scatter in the stellar mass and mean stellar population age results from variation in the baryonic mass, $M_{\mathrm{bary}} = M_{\mathrm{gas}} + M_{*}$, of the galaxies' progenitors at fixed halo mass and concentration. At the redshift of peak correlation ($z \approx 1$), the progenitor baryonic mass accounts for $75\%$ of the variance in the $z=0.1$ $M_{*}-V_{\mathrm{max, DMO}}$ relation. The scatter in the baryonic mass, in turn, is primarily set by differences in feedback strength and gas accretion over the course of the evolution of each halo.

We look into the abundance of Dual AGN in the largest hydrodynamical simulation from the EAGLE project. We define a Dual AGN as two active black holes (BHs) with a separation below 30 kpc. We find that only 1 per cent of AGN with $\LhX\geq 10^{42}\ergs$ are part of a Dual AGN system at $z=0.8-1$. During the evolution of a typical binary BH system, the rapid variability of the hard X-ray luminosity on Myr time-scales severely limits the detectability of Dual AGN. To quantify this effect, we calculate a probability of detection, $t_{\rm on}/t_{\rm 30}$, where $t_{\rm 30}$ is the time in which the two black holes are separated at distances below 30 kpc and $t_{\rm on}$, the time that both AGN are visible (e.g. when both AGN have $\LhX\geq 10^{42}\ergs$) in this period. We find that the average fraction of visible Dual systems is 3 per cent. The visible Dual AGN distribution as a function of BH separation presents a pronounced peak at $\sim 20$ kpc that can be understood as a result of the rapid orbital decay of the host galaxies after their first encounter. We also find that $75$ per cent of the host galaxies have recently undergone or are undergoing a merger with stellar mass ratio $\geq 0.1$. Finally, we find that the fraction of visible Dual AGN increases with redshift as found in observations.

We construct Tully-Fisher relations (TFRs), from large samples of galaxies with spatially-resolved H$\alpha$ emission maps from the K-band Multi-Object Spectrograph (KMOS) Redshift One Spectroscopic Survey (KROSS) at $z\approx1$. We compare these to data from the Sydney-Australian-Astronomical-Observatory Multi-object Integral-Field Spectrograph (SAMI) Galaxy Survey at $z\approx0$. We stringently match the data quality of the latter to the former, and apply identical analysis methods and sub-sample selection criteria to both to conduct a direct comparison of the absolute $K$-band magnitude and stellar mass TFRs at $z\approx1$ and $z\approx0$. We find that matching the quality of the SAMI data to that of KROSS results in TFRs that differ significantly in slope, zero-point and (sometimes) scatter in comparison to the corresponding original SAMI relations. These differences are in every case as large or larger than the differences between the KROSS $z\approx1$ and matched SAMI $z\approx0$ relations. Accounting for these differences, we compare the TFRs at $z\approx1$ and $z\approx0$. For disk-like, star-forming galaxies we find no significant difference in the TFR zero-points between the two epochs. This suggests the growth of stellar mass and dark matter in these types of galaxies is intimately linked over this $\approx8$ Gyr period.

Dark matter particles may decay, emitting photons. Drawing on the EAGLE family of hydrodynamic simulations of galaxy formation -- including the APOSTLE and C-EAGLE simulations -- we assess the systematic uncertainties and scatter on the decay flux from different galaxy classes, from Milky Way satellites to galaxy clusters, and compare our results to studies of the 3.55~keV line. We demonstrate that previous detections and non-detections of this line are consistent with a dark matter interpretation. For example, in our simulations the width of the the dark matter decay line for Perseus-analogue galaxy clusters lies in the range 1300-1700~\kms. Therefore, the non-detection of the 3.55~keV line in the centre of the Perseus cluster by the {\it Hitomi} collaboration is consistent with detections by other instruments. We also consider trends with stellar and halo mass and evaluate the scatter in the expected fluxes arising from the anisotropic halo mass distribution and from object-to-object variations. We provide specific predictions for observations with {\it XMM-Newton} and with the planned X-ray telescopes {\it XRISM} and {\it ATHENA}. If future detections of unexplained X-ray lines match our predictions, including line widths, we will have strong evidence that we have discovered the dark matter.

We present the first results from a study of OVI absorption around galaxies at $z<1.44$ using data from a near-infrared grism spectroscopic Hubble Space Telescope Large Program, the Quasar Sightline and Galaxy Evolution (QSAGE) survey. QSAGE is the first grism galaxy survey to focus on the circumgalactic medium at $z\sim1$, providing a blind survey of the galaxy population. Using the first of 12 fields, we provide details of the reduction methods, in particular the handling of the deep grism data which uses multiple position angles to minimise the effects of contamination from overlapping traces. The resulting galaxy sample is H$\alpha$ flux limited ($f({\rm H\alpha}) > 2\times10^{-17}$ erg s$^{-1}$ cm$^{-2}$) at 0.68<z<1.44, corresponding to $\gtrsim0.2-0.8$ M$_\odot$ yr$^{-1}$. We combine the galaxy data with high-resolution STIS and COS spectroscopy of the background quasar to study OVI in the circumgalactic medium. At z>0.68, we find 5 OVI absorption systems along the line of sight with identified galaxies lying at impact parameters of $b\approx100-350$ kpc (proper), whilst we find a further 13 galaxies with no significant associated OVI absorption (i.e. $N({\rm OVI})<10^{13.5-14}$ cm$^{-2}$) in the same impact parameter and redshift range. We find a large scatter in the stellar mass and star-formation rates of the closest galaxies with associated OVI. Whilst one of the OVI absorber systems is found to be associated with a low mass galaxy group at $z\approx1.08$, we infer that the detected OVI absorbers typically lie in the proximity of dark matter halos of masses $10^{11.5} {\rm M_\odot}\lesssim M_{\rm halo}\lesssim10^{12} {\rm M_\odot}$.

We investigate the rapid growth phase of supermassive black holes (BHs) within the hydrodynamical cosmological \eagle simulation. This non-linear phase of BH growth occurs within $\sim$$L_{*}$ galaxies, embedded between two regulatory states of the galaxy host: in sub $L_{*}$ galaxies efficient stellar feedback regulates the gas inflow onto the galaxy and significantly reduces the growth of the central BH, while in galaxies more massive than $L_{*}$ efficient AGN feedback regulates the gas inflow onto the galaxy and curbs further non-linear BH growth. We find evolving critical galaxy and halo mass scales at which rapid BH growth begins. Galaxies in the low-redshift Universe transition into the rapid BH growth phase in haloes that are approximately an order of magnitude more massive than their high-redshift counterparts (\M{200} $\approx 10^{12.4}$~\Msol at $z \approx 0$ decreasing to \M{200} $\approx 10^{11.2}$~\Msol at $z \approx 6$). Instead, BHs enter the rapid growth phase at a fixed critical halo virial temperature ($T_{\mathrm{vir}} \approx 10^{5.6}$~K). We additionally show that major galaxy--galaxy interactions ($\mu \geq \frac{1}{4}$, where $\mu$ is the stellar mass ratio) play a substantial role in triggering the rapid growth phase of BHs in the low-redshift Universe, whilst potentially having a lower influence at high redshift. Approximately 40\% of BHs that initiate the rapid BH growth phase at $z \approx 0$ do so within $\pm 0.5$ dynamical times of a major galaxy--galaxy merger, a fourfold increase above what is expected from the background merger rate. We find that minor mergers ($\frac{1}{10} \leq \mu < \frac{1}{4}$) have a substantially lower influence in triggering the rapid growth phase at all epochs.

Cosmological simulations require the use of a multiple time-stepping scheme. Without such a scheme, cosmological simulations would be impossible due to their high level of dynamic range; over eleven orders of magnitude in density. Such a large dynamic range leads to a range of over four orders of magnitude in time-step, which presents a significant load-balancing challenge. In this work, the extreme adaptivity that cosmological simulations present is tackled in three main ways through the use of the code SWIFT. First, an adaptive mesh is used to ensure that only the relevant particles are interacted in a given time-step. Second, task-based parallelism is used to ensure efficient load-balancing within a single node, using pthreads and SIMD vectorisation. Finally, a domain decomposition strategy is presented, using the graph domain decomposition library METIS, that bisects the work that must be performed by the simulation between nodes using MPI. These three strategies are shown to give SWIFT near-perfect weak-scaling characteristics, only losing 25% performance when scaling from 1 to 4096 cores on a representative problem, whilst being more than 30x faster than the de-facto standard Gadget-2 code.

We use the EAGLE simulations to study the oxygen abundance gradients of gas discs in galaxies within the stellar mass range [10^9.5, 10^10.8]Mo at z=0. The estimated median oxygen gradient is -0.011 (0.002) dex kpc^-1, which is shallower than observed. No clear trend between simulated disc oxygen gradient and galaxy stellar mass is found when all galaxies are considered. However, the oxygen gradient shows a clear correlation with gas disc size so that shallower abundance slopes are found for increasing gas disc sizes. Positive oxygen gradients are detected for ~40 per cent of the analysed gas discs, with a slight higher frequency in low mass galaxies. Galaxies that have quiet merger histories show a positive correlation between oxygen gradient and stellar mass, so that more massive galaxies tend to have shallower metallicity gradients. At high stellar mass, there is a larger fraction of rotational-dominated galaxies in low density regions. At low stellar mass, non-merger galaxies show a large variety of oxygen gradients and morphologies. The normalization of the disc oxygen gradients in non-merger galaxies by the effective radius removes the trend with stellar mass. Conversely, galaxies that experienced mergers show a weak relation between oxygen gradient and stellar mass. Additionally, the analysed EAGLE discs show no clear dependence of the oxygen gradients on local environment, in agreement with current observational findings.

We study the correlation between stellar mass and stellar metallicity in the Evolution and Assembly of GaLaxies and their Environments (EAGLE) suite of cosmological hydrodynamical simulations. At a given stellar mass, simulated galaxies with lower stellar metallicities show, on average, higher gas fractions, higher specific star formation rates and younger stellar populations. Active galactic nuclei feedback seems to play an important role on the determination of the stellar metallicity at high stellar masses. In general, simulated systems follow a well-defined anticorrelation between stellar metallicity and gas fraction, which does not evolve significantly with redshift. All these trends are consistent with previous findings regarding the metallicity of the star-forming gas in EAGLE.

We investigate the effect of the accelerated expansion of the Universe due to a cosmological constant, $\Lambda$, on the cosmic star formation rate. We utilise hydrodynamical simulations from the EAGLE suite, comparing a $\Lambda$CDM Universe to an Einstein-de Sitter model with $\Lambda=0$. Despite the differences in the rate of growth of structure, we find that dark energy, at its observed value, has negligible impact on star formation in the Universe. We study these effects beyond the present day by allowing the simulations to run forward into the future ($t>13.8$ Gyr). We show that the impact of $\Lambda$ becomes significant only when the Universe has already produced most of its stellar mass, only decreasing the total co-moving density of stars ever formed by ${\approx}15\%$. We develop a simple analytic model for the cosmic star formation rate that captures the suppression due to a cosmological constant. The main reason for the similarity between the models is that feedback from accreting black holes dramatically reduces the cosmic star formation at late times. Interestingly, simulations without feedback from accreting black holes predict an upturn in the cosmic star formation rate for $t>15$ Gyr due to the rejuvenation of massive ($ > 10^{11} \mathrm{M}_{\odot}$) galaxies. We briefly discuss the implication of the weak dependence of the cosmic star formation on $\Lambda$ in the context of the anthropic principle.

We investigate the origin of the Hubble sequence by analysing the evolution of the kinematic morphologies of central galaxies in the EAGLE cosmological simulation. By separating each galaxy into disk and spheroidal stellar components and tracing their evolution along the merger tree, we find that the morphology of galaxies follows a common evolutionary trend. We distinguish three phases of galaxy formation. These phases are determined primarily by mass, rather than redshift. For M_star < 10^{9.5} M_sun galaxies grow in a disorganised way, resulting in a morphology that is dominated by random stellar motions. This phase is dominated by in-situ star formation, partly triggered by mergers. In the mass range 10^{9.5} M_sun < M_star < 10^{10.5} M_sun galaxies evolve towards a disk-dominated morphology, driven by in-situ star formation. The central spheroid (i.e. the bulge) at z = 0 consists mostly of stars that formed in-situ, yet the formation of the bulge is to a large degree associated with mergers. Finally, at M_star > 10^{10.5} M_sun growth through in-situ star formation slows down considerably and galaxies transform towards a more spheroidal morphology. This transformation is driven more by the buildup of spheroids than by the destruction of disks. Spheroid formation in these galaxies happens mostly by accretion at large radii of stars formed ex-situ (i.e. the halo rather than the bulge).

Models of the very early universe, including inflationary models, are argued to produce varying universe domains with different values of fundamental constants and cosmic parameters. Using the cosmological hydrodynamical simulation code from the eagle collaboration, we investigate the effect of the cosmological constant on the formation of galaxies and stars. We simulate universes with values of the cosmological constant ranging from Lambda = 0 to Lambda_0 = 300, where Lambda_0 is the value of the cosmological constant in our Universe. Because the global star formation rate in our Universe peaks at t = 3.5 Gyr, before the onset of accelerating expansion, increases in Lambda of even an order of magnitude have only a small effect on the star formation history and efficiency of the universe. We use our simulations to predict the observed value of the cosmological constant, given a measure of the multiverse. Whether the cosmological constant is successfully predicted depends crucially on the measure. The impact of the cosmological constant on the formation of structure in the universe does not seem to be a sharp enough function of Lambda to explain its observed value alone.

We present the first simulated galaxy clusters (M_200 > 10^14 Msun) with both self-interacting dark matter (SIDM) and baryonic physics. They exhibit a greater diversity in both dark matter and stellar density profiles than their counterparts in simulations with collisionless dark matter (CDM), which is generated by the complex interplay between dark matter self-interactions and baryonic physics. Despite variations in formation history, we demonstrate that analytical Jeans modelling predicts the SIDM density profiles remarkably well, and the diverse properties of the haloes can be understood in terms of their different final baryon distributions.

We present a comparison of nine galaxy formation models, eight semi-analytical and one halo occupation distribution model, run on the same underlying cold dark matter simulation (cosmological box of co-moving width 125$h^{-1}$ Mpc, with a dark-matter particle mass of $1.24\times 10^9 h^{-1}$ Msun) and the same merger trees. While their free parameters have been calibrated to the same observational data sets using two approaches, they nevertheless retain some 'memory' of any previous calibration that served as the starting point (especially for the manually-tuned models). For the first calibration, models reproduce the observed z = 0 galaxy stellar mass function (SMF) within 3-{\sigma}. The second calibration extended the observational data to include the z = 2 SMF alongside the z~0 star formation rate function, cold gas mass and the black hole-bulge mass relation. Encapsulating the observed evolution of the SMF from z = 2 to z = 0 is found to be very hard within the context of the physics currently included in the models. We finally use our calibrated models to study the evolution of the stellar-to-halo mass (SHM) ratio. For all models we find that the peak value of the SHM relation decreases with redshift. However, the trends seen for the evolution of the peak position as well as the mean scatter in the SHM relation are rather weak and strongly model dependent. Both the calibration data sets and model results are publicly available.

We present sensitive 870$\mu$m continuum measurements from our ALMA programmes of 114 X-ray selected AGN in the CDF-S and COSMOS fields. We use these observations in combination with data from Spitzer and Herschel to construct a sample of 86 X-ray selected AGN, 63 with ALMA constraints at $z=1.5-3.2$ with stellar mass $>2\times10^{10}M_{\odot}$. We constructed broad-band spectral energy distributions in the infrared band (8 - 1000$\mu$m) and constrain star-formation rates (SFRs) uncontaminated by the AGN. Using a hierarchical Bayesian method that takes into account the information from upper limits, we fit SFR and specific SFR (sSFR) distributions. We explore these distributions as a function of both X-ray luminosity and stellar mass. We compare our measurements to two versions of the EAGLE hydrodynamical simulations: the reference model with AGN feedback and the model without AGN. We find good agreement between the observations and that predicted by the EAGLE reference model for the modes and widths of the sSFR distributions as a function of both X-ray luminosity and stellar mass; however, we found that the EAGLE model without AGN feedback predicts a significantly narrower width when compared to the data. Overall, from the combination of the observations with the model predictions, we conclude that (1) even with AGN feedback, we expect no strong relationship between the sSFR distribution parameters and instantaneous AGN luminosity and (2) a signature of AGN feedback is a broad distribution of sSFRs for all galaxies (not just those hosting an AGN) with stellar masses above $\approx 10^{10}$M$_{\odot}$.

We analyse the velocity dispersion properties of 472 z~0.9 star-forming galaxies observed as part of the KMOS Redshift One Spectroscopic Survey (KROSS). The majority of this sample is rotationally dominated (83 +/- 5% with v_C/sigma_0 > 1) but also dynamically hot and highly turbulent. After correcting for beam smearing effects, the median intrinsic velocity dispersion for the final sample is sigma_0 = 43.2 +/- 0.8 km/s with a rotational velocity to dispersion ratio of v_C/sigma_0 = 2.6 +/- 0.1. To explore the relationship between velocity dispersion, stellar mass, star formation rate and redshift we combine KROSS with data from the SAMI survey (z~0.05) and an intermediate redshift MUSE sample (z~0.5). While there is, at most, a weak trend between velocity dispersion and stellar mass, at fixed mass there is a strong increase with redshift. At all redshifts, galaxies appear to follow the same weak trend of increasing velocity dispersion with star formation rate. Our results are consistent with an evolution of galaxy dynamics driven by disks that are more gas rich, and increasingly gravitationally unstable, as a function of increasing redshift. Finally, we test two analytic models that predict turbulence is driven by either gravitational instabilities or stellar feedback. Both provide an adequate description of the data, and further observations are required to rule out either model.

We use EAGLE to quantify the effect galaxy mergers have on the stellar specific angular momentum of galaxies, $j_{\rm stars}$. We split mergers into: dry (gas-poor)/wet (gas-rich), major/minor, and different spin alignments and orbital parameters. Wet (dry) mergers have an average neutral gas-to-stellar mass ratio of $1.1$ ($0.02$), while major (minor) mergers are those with stellar mass ratios $\ge 0.3$ ($0.1-0.3$). We correlate the positions of galaxies in the $j_{\rm stars}$-stellar mass plane at $z=0$ with their merger history, and find that galaxies of low spins suffered dry mergers, while galaxies of normal/high spins suffered predominantly wet mergers, if any. The radial $j_{\rm stars}$ profiles of galaxies that went through dry mergers are deficient by $\approx 0.3$~dex at $r\lesssim 10\,r_{50}$ (with $r_{50}$ being the half-stellar mass radius), compared to galaxies that went through wet mergers. Studying the merger remnants reveals that dry mergers reduce $j_{\rm stars}$ by $\approx 30$\%, while wet mergers increase it by $\approx 10$\%, on average. The latter is connected to the build-up of the bulge by newly formed stars of high rotational speed. Moving from minor to major mergers accentuates these effects. When the spin vectors of the galaxies prior to the dry merger are misaligned, $j_{\rm stars}$ decreases to a greater magnitude, while in wet mergers co-rotation and high orbital angular momentum efficiently spun-up galaxies. We predict what would be the observational signatures in the $j_{\rm stars}$ profiles driven by dry mergers: (i) shallow radial profiles and (ii) profiles that rise beyond $\approx 10\,r_{50}$, both of which are significantly different from spiral galaxies.

It is now possible for hydrodynamical simulations to reproduce a representative galaxy population. Accordingly, it is timely to assess critically some of the assumptions of traditional semi-analytic galaxy formation models. We use the Eagle simulations to assess assumptions built into the Galform semi-analytic model, focussing on those relating to baryon cycling, angular momentum and feedback. We show that the assumption in Galform that newly formed stars have the same specific angular momentum as the total disc leads to a significant overestimate of the total stellar specific angular momentum of disc galaxies. In Eagle, stars form preferentially out of low specific angular momentum gas in the interstellar medium (ISM) due to the assumed gas density threshold for stars to form, leading to more realistic galaxy sizes. We find that stellar mass assembly is similar between Galform and Eagle but that the evolution of gas properties is different, with various indications that the rate of baryon cycling in Eagle is slower than is assumed in Galform. Finally, by matching individual galaxies between Eagle and Galform, we find that an artificial dependence of AGN feedback and gas infall rates on halo mass doubling events in Galform drives most of the scatter in stellar mass between individual objects. Put together our results suggest that the Galform semi-analytic model can be significantly improved in light of recent advances.

This work presents a study of galactic outflows driven by stellar feedback. We extract main sequence disc galaxies with stellar mass $10^9\le$ M$_{\star}/$M$_{\odot} \le 5.7\times10^{10}$ at redshift $z=0$ from the highest resolution cosmological simulation of the Evolution and Assembly of GaLaxies and their Environments (EAGLE) set. Synthetic gas rotation velocity and velocity dispersion ($\sigma$) maps are created and compared to observations of disc galaxies obtained with the Sydney-AAO Multi-object Integral field spectrograph (SAMI), where $\sigma$-values greater than $150$ km s$^{-1}$ are most naturally explained by bipolar outflows powered by starburst activity. We find that the extension of the simulated edge-on (pixelated) velocity dispersion probability distribution depends on stellar mass and star formation rate surface density ($\Sigma_{\rm SFR}$), with low-M$_{\star}/$low-$\Sigma_{\rm SFR}$ galaxies showing a narrow peak at low $\sigma$ ($\sim30$ km s$^{-1}$) and more active, high-M$_{\star}/$high-$\Sigma_{\rm SFR}$ galaxies reaching $\sigma>150$ km s$^{-1}$. Although supernova-driven galactic winds in the EAGLE simulations may not entrain enough gas with T $<10^5$ K compared to observed galaxies, we find that gas temperature is a good proxy for the presence of outflows. There is a direct correlation between the thermal state of the gas and its state of motion as described by the $\sigma$-distribution. The following equivalence relations hold in EAGLE: $i)$ low-$\sigma$ peak $\,\Leftrightarrow\,$ disc of the galaxy $\,\Leftrightarrow\,$ gas with T $<10^5$ K; $ii)$ high-$\sigma$ tail $\,\Leftrightarrow\,$ galactic winds $\,\Leftrightarrow\,$ gas with T $\ge 10^5$ K.

We use a suite of hydrodynamical cosmological simulations from the Evolution and Assembly of GaLaxies and their Environments (EAGLE) project to investigate the formation of hot hydrostatic haloes and their dependence on feedback mechanisms. We find that the appearance of a strong bimodality in the probability density function (PDF) of the ratio of the radiative cooling and dynamical times for halo gas provides a clear signature of the formation of a hot corona. Haloes of total mass $10^{11.5}-10^{12}{\rm{M}}_{\odot}$ develop a hot corona independent of redshift, at least in the interval $z=0-4$ where the simulation has sufficiently good statistics. We analyse the build up of the hot gas mass in the halo, $M_{\rm{hot}}$, as a function of halo mass and redshift and find that while more energetic galactic winds powered by SNe increases $M_{\rm{hot}}$, AGN feedback reduces it by ejecting gas from the halo. We also study the thermal properties of gas accreting onto haloes and measure the fraction of shock-heated gas as a function of redshift and halo mass. We develop analytic and semianalytic approaches to estimate a `critical halo mass', $M_{\rm{crit}}$, for hot halo formation. We find that the mass for which the heating rate produced by accretion shocks equals the radiative cooling rate, reproduces the mass above which haloes develop a significant hot atmosphere. This yields a mass estimate of $M_{\rm{crit}} \approx 10^{11.7}{\rm{M}}_{\odot}$ at $z=0$, which agrees with the simulation results. The value of $M_{\rm{crit}}$ depends more strongly on the cooling rate than on any of the feedback parameters.

We quantify the correlations between gas-phase and stellar metallicities and global properties of galaxies, such as stellar mass, halo mass, age and gas fraction, in the Evolution and Assembly of GaLaxies and their Environments (EAGLE) suite of cosmological hydrodynamical simulations. The slope of the correlation between stellar mass and metallicity of star-forming (SF) gas ($M_*-Z_{\rm SF,gas}$ relation) depends somewhat on resolution, with the higher-resolution run reproducing a steeper slope. This simulation predicts a non-zero metallicity evolution, increasing by $\approx 0.5$ dex at $\sim 10^9 {\rm M}_{\odot}$ since $z = 3$. The simulated relation between stellar mass, metallicity and star formation rate at $z \lesssim 5$ agrees remarkably well with the observed fundamental metallicity relation. At $M_* \lesssim 10^{10.3} {\rm M}_{\odot}$ and fixed stellar mass, higher metallicities are associated with lower specific star formation rates, lower gas fractions and older stellar populations. On the other hand, at higher $M_*$, there is a hint of an inversion of the dependence of metallicity on these parameters. The fundamental parameter that best correlates with the metal content, in the simulations, is the gas fraction. The simulated gas fraction-metallicity relation exhibits small scatter and does not evolve significantly since $z = 3$. In order to better understand the origin of these correlations, we analyse a set of lower resolution simulations in which feedback parameters are varied. We find that the slope of the simulated $M_*-Z_{\rm SF,gas}$ relation is mostly determined by stellar feedback at low stellar masses ($M_* \lesssim 10^{10} {\rm M}_{\odot}$), and at high masses ($M_* \gtrsim 10^{10} {\rm M}_{\odot}$) by the feedback from active galactic nuclei.

We make use of a catalog of 1600 Pan-STARRS1 groups produced by the probability friends-of-friends algorithm to explore how the galaxy properties, i.e. the specific star formation rate (SSFR) and quiescent fraction, depend on stellar mass and group-centric radius. The work is the extension of Lin et al. (2014). In this work, powered by a stacking technique plus a background subtraction for contamination removal, a finer correction and more precise results are obtained than in our previous work. We find that while the quiescent fraction increases with decreasing group-centric radius the median SSFRs of star-forming galaxies in groups at fixed stellar mass drop slightly from the field toward the group center. This suggests that the major quenching process in groups is likely a fast mechanism. On the other hand, a reduction in SSFRs by ~0.2 dex is seen inside clusters as opposed to the field galaxies. If the reduction is attributed to the slow quenching effect, the slow quenching process acts dominantly in clusters. In addition, we also examine the density-color relation, where the density is defined by using a sixth-nearest neighbor approach. Comparing the quiescent fractions contributed from the density and radial effect, we find that the density effect dominates over the massive group or cluster galaxies, and the radial effect becomes more effective in less massive galaxies. The results support mergers and/or starvation as the main quenching mechanisms in the group environment, while harassment and/or starvation dominate in clusters.

We investigate the evolution of the galaxy Star Formation Rate Function (SFRF) and Cosmic Star Formation Rate Density (CSFRD) of $z\sim 0-8 $ galaxies in the Evolution and Assembly of GaLaxies and their Environments (EAGLE) simulations. In addition, we present a compilation of UV, IR and H$\alpha$ SFRFs and compare these with the predictions from the EAGLE suite of cosmological hydrodynamic simulations. We find that the constraints implied by different indicators are inconsistent with each other for the highest star-forming objects at z < 2, a problem that is possibly related to selection biases and the uncertainties of dust attenuation effects. EAGLE's feedback parameters were calibrated to reproduce realistic galaxy sizes and stellar masses at z = 0.1. In this work we test if and why those choices yield realistic Star Formation Rates (SFRs) for $z \sim 0-8$ as well. We demonstrate that SNe feedback plays a major role at setting the abundance of galaxies at all star-forming regimes, especially at high redshifts. On the contrary, Active Galactic Nuclei (AGN) feedback becomes more prominent at lower redshifts and is a major mechanism that affects only the highest star-forming systems. Furthermore, we find that galaxies with SFR $\sim 1-10 \, {\rm M_{\odot} \, yr^{-1}}$ dominate the CSFRD at redshifts z < 5, while rare high star-forming galaxies (SFR $\sim 10-100 \,{\rm M_{\odot} \, yr^{-1}}$) contribute significantly only briefly around the peak era ($z \sim 2$) and then are quenched by AGN feedback. In the absence of this prescription objects with SFR $\sim 10-100 \,{\rm M_{\odot} \, yr^{-1}}$ would dominate the CSFRD, while the cosmic budget of star formation would be extremely high. Finally, we demonstrate that the majority of the cosmic star formation occurs in relatively rare high mass halos ($ {\rm M_{Halo}} \sim 10^{11-13} \, {\rm M_{\odot}}$) even at the earliest epochs.

We introduce the Cluster-EAGLE (C-EAGLE) simulation project, a set of cosmological hydrodynamical zoom simulations of the formation of $30$ galaxy clusters in the mass range $10^{14}<M_{200}/\mathrm{M}_{\odot}<10^{15.4}$ that incorporates the Hydrangea sample of Bah\'e et al. (2017). The simulations adopt the state-of-the-art EAGLE galaxy formation model, with a gas particle mass of $1.8\times10^{6}\,\mathrm{M}_{\odot}$ and physical softening length of $0.7\,\mathrm{kpc}$. In this paper, we introduce the sample and present the low-redshift global properties of the clusters. We calculate the X-ray properties in a manner consistent with observational techniques, demonstrating the bias and scatter introduced by using estimated masses. We find the total stellar content and black hole masses of the clusters to be in good agreement with the observed relations. However, the clusters are too gas rich, suggesting that the AGN feedback model is not efficient enough at expelling gas from the high-redshift progenitors of the clusters. The X-ray properties, such as the spectroscopic temperature and the soft-band luminosity, and the Sunyaev-Zel'dovich properties are in reasonable agreement with the observed relations. However, the clusters have too high central temperatures and larger-than-observed entropy cores, which is likely driven by the AGN feedback after the cluster core has formed. The total metal content and its distribution throughout the ICM are a good match to the observations.

We describe a new Large Program in progress on the Gemini North and South telescopes: Gemini Observations of Galaxies in Rich Early Environments (GOGREEN). This is an imaging and deep spectroscopic survey of 21 galaxy systems at $1<z<1.5$, selected to span a factor $>10$ in halo mass. The scientific objectives include measuring the role of environment in the evolution of low-mass galaxies, and measuring the dynamics and stellar contents of their host haloes. The targets are selected from the SpARCS, SPT, COSMOS and SXDS surveys, to be the evolutionary counterparts of today's clusters and groups. The new red-sensitive Hamamatsu detectors on GMOS, coupled with the nod-and-shuffle sky subtraction, allow simultaneous wavelength coverage over $\lambda\sim 0.6$--$1.05\mu$m, and this enables a homogeneous and statistically complete redshift survey of galaxies of all types. The spectroscopic sample targets galaxies with AB magnitudes $z^{\prime}<24.25$ and [3.6]$\mu$m$<22.5$, and is therefore statistically complete for stellar masses $M_\ast\gtrsim10^{10.3}M_\odot$, for all galaxy types and over the entire redshift range. Deep, multiwavelength imaging has been acquired over larger fields for most systems, spanning $u$ through $K$, in addition to deep IRAC imaging at 3.6$\mu$m. The spectroscopy is $\sim 50$ per cent complete as of semester 17A, and we anticipate a final sample of $\sim 500$ new cluster members. Combined with existing spectroscopy on the brighter galaxies from GCLASS, SPT and other sources, GOGREEN will be a large legacy cluster and field galaxy sample at this redshift that spectroscopically covers a wide range in stellar mass, halo mass, and clustercentric radius.

Observations with the Chandra X-ray Observatory are used to examine the hot gas properties within a sample of 10 galaxy groups selected from the Galaxy And Mass Assembly survey's optical Friends-of-Friends group catalogue. Our groups have been screened to eliminate spurious and unrelaxed systems, and the effectiveness of this procedure is demonstrated by the detection of intergalactic hot gas in 80 per cent of our sample. However, we find that 9 of the 10 are X-ray underluminous by a mean factor of $\sim$4 compared to typical X-ray-selected samples. Consistent with this, the majority of our groups have gas fractions that are lower and gas entropies somewhat higher than those seen in typical X-ray-selected samples. Two groups, which have high 2{\sigma} lower limits on their gas entropy, are candidates for the population of high-entropy groups predicted by some active galactic nucleus feedback models.

We study present-day galaxy clustering in the EAGLE cosmological hydrodynamical simulation. EAGLE's galaxy formation parameters were calibrated to reproduce the redshift $z=0.1$ galaxy stellar mass function, and the simulation also reproduces galaxy colours well. The simulation volume is too small to correctly sample large-scale fluctuations and we therefore concentrate on scales smaller than a few mega parsecs. We find very good agreement with observed clustering measurements from the Galaxy And Mass Assembly (GAMA) survey, when galaxies are binned by stellar mass, color, or luminosity. However, low-mass red-galaxies are clustered too strongly, which is at least partly due to limited numerical resolution. Apart from this limitation, we conclude that EAGLE galaxies inhabit similar dark matter haloes as observed GAMA galaxies, and that the radial distribution of satellite galaxies as function of stellar mass and colour is similar to that observed as well.

We present mock optical images, broad-band and H$\alpha$ fluxes, and D4000 spectral indices for 30,145 galaxies from the EAGLE hydrodynamical simulation at redshift $z=0.1$, modelling dust with the SKIRT Monte Carlo radiative transfer code. The modelling includes a subgrid prescription for dusty star-forming regions, with both the subgrid obscuration of these regions and the fraction of metals in diffuse interstellar dust calibrated against far-infrared fluxes of local galaxies. The predicted optical colours as a function of stellar mass agree well with observation, with the SKIRT model showing marked improvement over a simple dust screen model. The orientation dependence of attenuation is weaker than observed because EAGLE galaxies are generally puffier than real galaxies, due to the pressure floor imposed on the interstellar medium. The mock H$\alpha$ luminosity function agrees reasonably well with the data, and we quantify the extent to which dust obscuration affects observed H$\alpha$ fluxes. The distribution of D4000 break values is bimodal, as observed. In the simulation, 20$\%$ of galaxies deemed `passive' for the SKIRT model, i.e. exhibiting D4000 $> 1.8$, are classified `active' when ISM dust attenuation is not included. The fraction of galaxies with stellar mass greater than $10^{10}$ M$_\odot$ that are deemed passive is slightly smaller than observed, which is due to low levels of residual star formation in these simulated galaxies. Colour images, fluxes and spectra of EAGLE galaxies are to be made available through the public EAGLE database.

We investigate the relation between kinematic morphology, intrinsic colour and stellar mass of galaxies in the EAGLE cosmological hydrodynamical simulation. We calculate the intrinsic u-r colours and measure the fraction of kinetic energy invested in ordered corotation of 3562 galaxies at z=0 with stellar masses larger than $10^{10}M_{\odot}$. We perform a visual inspection of gri-composite images and find that our kinematic morphology correlates strongly with visual morphology. EAGLE produces a galaxy population for which morphology is tightly correlated with the location in the colour- mass diagram, with the red sequence mostly populated by elliptical galaxies and the blue cloud by disc galaxies. Satellite galaxies are more likely to be on the red sequence than centrals, and for satellites the red sequence is morphologically more diverse. These results show that the connection between mass, intrinsic colour and morphology arises from galaxy formation models that reproduce the observed galaxy mass function and sizes.

We present a clustering comparison of 12 galaxy formation models (including Semi-Analytic Models (SAMs) and Halo Occupation Distribution (HOD) models) all run on halo catalogues and merger trees extracted from a single {\Lambda}CDM N-body simulation. We compare the results of the measurements of the mean halo occupation numbers, the radial distribution of galaxies in haloes and the 2-Point Correlation Functions (2PCF). We also study the implications of the different treatments of orphan (galaxies not assigned to any dark matter subhalo) and non-orphan galaxies in these measurements. Our main result is that the galaxy formation models generally agree in their clustering predictions but they disagree significantly between HOD and SAMs for the orphan satellites. Although there is a very good agreement between the models on the 2PCF of central galaxies, the scatter between the models when orphan satellites are included can be larger than a factor of 2 for scales smaller than 1 Mpc/h. We also show that galaxy formation models that do not include orphan satellite galaxies have a significantly lower 2PCF on small scales, consistent with previous studies. Finally, we show that the 2PCF of orphan satellites is remarkably different between SAMs and HOD models. Orphan satellites in SAMs present a higher clustering than in HOD models because they tend to occupy more massive haloes. We conclude that orphan satellites have an important role on galaxy clustering and they are the main cause of the differences in the clustering between HOD models and SAMs.

We analyze the total and baryonic acceleration profiles of a set of well-resolved galaxies identified in the EAGLE suite of hydrodynamic simulations. Our runs start from the same initial conditions but adopt different prescriptions for unresolved stellar and AGN feedback, resulting in diverse populations of galaxies by the present day. Some of them reproduce observed galaxy scaling relations, while others do not. However, regardless of the feedback implementation, all of our galaxies follow closely a simple relationship between the total and baryonic acceleration profiles, consistent with recent observations of rotationally supported galaxies. The relation has small scatter: different feedback implementations -- which produce different galaxy populations -- mainly shift galaxies along the relation, rather than perpendicular to it. Furthermore, galaxies exhibit a characteristic acceleration, $g_{\dagger}$, above which baryons dominate the mass budget, as observed. These observations, consistent with simple modified Newtonian dynamics, can be accommodated within the standard cold dark matter paradigm.

We introduce the Hydrangea simulations, a suite of 24 cosmological hydrodynamic zoom-in simulations of massive galaxy clusters (M_200c = 10^14-10^15 M_Sun) with baryon particle masses of ~10^6 M_Sun. Designed to study the impact of the cluster environment on galaxy formation, they are a key part of the `Cluster-EAGLE' project (Barnes et al. 2017). They use a galaxy formation model developed for the EAGLE project, which has been shown to yield both realistic field galaxies and hot gas fractions of galaxy groups consistent with observations. The total stellar mass content of the simulated clusters agrees with observations, but central cluster galaxies are too massive, by up to 0.6 dex. Passive satellite fractions are higher than in the field, and at stellar masses Mstar > 10^10 M_Sun this environmental effect is quantitatively consistent with observations. The predicted satellite stellar mass function matches data from local cluster surveys. Normalized to total mass, there are fewer low-mass (Mstar < 10^10 M_Sun) galaxies within the virial radius of clusters than in the field, primarily due to star formation quenching. Conversely, the simulations predict an overabundance of massive galaxies in clusters compared to the field that persists to their far outskirts (> 5r_200c). This is caused by a significantly increased stellar mass fraction of (sub-)haloes in the cluster environment, by up to ~0.3 dex even well beyond r_200c. Haloes near clusters are also more concentrated than equally massive field haloes, but these two effects are largely uncorrelated.

We investigate the connection between the star formation rate (SFR) of galaxies and their central black hole accretion rate (BHAR) using the EAGLE cosmological hydrodynamical simulation. We find, in striking concurrence with recent observational studies, that the <SFR>--BHAR relation for an AGN selected sample produces a relatively flat trend, whilst the <BHAR>--SFR relation for a SFR selected sample yields an approximately linear trend. These trends remain consistent with their instantaneous equivalents even when both SFR and BHAR are time-averaged over a period of 100~Myr. There is no universal relationship between the two growth rates. Instead, SFR and BHAR evolve through distinct paths that depend strongly on the mass of the host dark matter halo. The galaxies hosted by haloes of mass M200 $\lesssim 10^{11.5}$Msol grow steadily, yet black holes (BHs) in these systems hardly grow, yielding a lack of correlation between SFR and BHAR. As haloes grow through the mass range $10^{11.5} \lesssim$ M200 $\lesssim 10^{12.5 }$Msol BHs undergo a rapid phase of non-linear growth. These systems yield a highly non-linear correlation between the SFR and BHAR, which are non-causally connected via the mass of the host halo. In massive haloes (M200 $\gtrsim 10^{12.5}$Msol) both SFR and BHAR decline on average with a roughly constant scaling of SFR/BHAR $\sim 10^{3}$. Given the complexity of the full SFR--BHAR plane built from multiple behaviours, and from the large dynamic range of BHARs, we find the primary driver of the different observed trends in the <SFR>--BHAR and <BHAR>--SFR relationships are due to sampling considerably different regions of this plane.

We investigate the properties of the galaxies that reionized the Universe and the history of cosmic reionization using the "Evolution and Assembly of GaLaxies and their environments" (EAGLE) cosmological hydrodynamical simulations. We obtain the evolution of the escape fraction of ionizing photons in galaxies assuming that galactic winds create channels through which 20~percent of photons escape when the local surface density of star formation is greater than $0.1$ M$_\odot$ yr$^{-1}$ kpc$^{-2}$. Such threshold behaviour for the generation of winds is observed, and the rare local objects which have such high star formation surface densities exhibit high escape fractions of $\sim 10$ percent. In our model the luminosity-weighted mean escape fraction increases with redshift as $\bar f_{\rm esc}=0.045~((1+z)/4)^{1.1}$ at $z>3$, and the galaxy number weighted mean as $\langle f_{\rm esc} \rangle=2.2\times10^{-3}~((1+z)/4)^4$, and becomes constant $\approx0.2$ at redshift $z>10$. The escape fraction evolves as an increasingly large fraction of stars forms above the critical surface density of star formation at earlier times. This evolution of the escape fraction, combined with that of the star formation rate density from EAGLE, reproduces the inferred evolution of the filling factor of ionized regions during the reionization epoch ($6<z<8$), the evolution of the post-reionization ($0\leq z<6$) hydrogen photoionisation rate, and the optical depth due to Thomson scattering of the cosmic microwave background photons measured by the Planck satellite.

We present a MUSE and KMOS dynamical study 405 star-forming galaxies at redshift z=0.28-1.65 (median redshift z=0.84). Our sample are representative of star-forming, main-sequence galaxies, with star-formation rates of SFR=0.1-30Mo/yr and stellar masses M=10^8-10^11Mo. For 49+/-4% of our sample, the dynamics suggest rotational support, 24+/-3% are unresolved systems and 5+/-2% appear to be early-stage major mergers with components on 8-30kpc scales. The remaining 22+/-5% appear to be dynamically complex, irregular (or face-on systems). For galaxies whose dynamics suggest rotational support, we derive inclination corrected rotational velocities and show these systems lie on a similar scaling between stellar mass and specific angular momentum as local spirals with j*=J/M*\propto M^(2/3) but with a redshift evolution that scales as j*\propto M^{2/3}(1+z)^(-1). We identify a correlation between specific angular momentum and disk stability such that galaxies with the highest specific angular momentum, log(j*/M^(2/3))>2.5, are the most stable, with Toomre Q=1.10+/-0.18, compared to Q=0.53+/-0.22 for galaxies with log(j*/M^(2/3))<2.5. At a fixed mass, the HST morphologies of galaxies with the highest specific angular momentum resemble spiral galaxies, whilst those with low specific angular momentum are morphologically complex and dominated by several bright star-forming regions. This suggests that angular momentum plays a major role in defining the stability of gas disks: at z~1, massive galaxies that have disks with low specific angular momentum, appear to be globally unstable, clumpy and turbulent systems. In contrast, galaxies with high specific angular have evolved in to stable disks with spiral structures.

We present dynamical measurements for 586 H-alpha detected star-forming galaxies from the KMOS (K-band Multi-Object Spectrograph) Redshift One Spectroscopic Survey (KROSS). The sample represents typical star-forming galaxies at this redshift (z=0.6-1.0), with a median star formation rate of ~7 Msol/yr and a stellar mass range of log[M/Msol]~9-11. We find that the rotation velocity-stellar mass relationship (the inverse of the Tully-Fisher relationship) for our rotationally-dominated sources (v/sigma>1) has a consistent slope and normalisation as that observed for z=0 disks. In contrast, the specific angular momentum (j; angular momentum divided by stellar mass), is ~0.2-0.3 dex lower on average compared to z=0 disks. The specific angular momentum scales as M^[0.6+/-0.2], consistent with that expected for dark matter (i.e., proportional to M^[2/3]). We find that z~0.9 star-forming galaxies have decreasing specific angular momentum with increasing Sersic index. Visually, the sources with the highest specific angular momentum, for a given mass, have the most disk-dominated morphologies. This implies that an angular momentum-mass-morphology relationship, similar to that observed in local massive galaxies, is already in place by z~1.

We estimate the rate of dark matter scattering in collapsed structures throughout the history of the Universe. If the scattering cross-section is velocity-independent, then the canonical picture is correct that scatterings occur mainly at late times. The scattering rate peaks slightly at redshift z~6, and remains significant today. Half the scatterings occur after z~1, in structures more massive than 10^12 M_sun. Within a factor of two, these numbers are robust to changes in the assumed astrophysics, and the scatterings would be captured in cosmological simulations. However, for particle physics models with a velocity-dependent cross-section (as for Yukawa potential interactions via a massive mediator), the scattering rate peaks before z~20, in objects with mass less than 10^4 M_sun. These precise values are sensitive to the redshift-dependent mass-concentration relation and the small-scale cutoff in the matter power spectrum. In extreme cases, the qualitative effect of early interactions may be reminiscent of warm dark matter and strongly affect the subsequent growth of structure. However, these scatterings are being missed in existing cosmological simulations with limited mass resolution.

Galaxies that have fallen into massive haloes may no longer be able to accrete gas from their surroundings, a process referred to as 'starvation' or 'strangulation' of satellites. We study the environmental dependence of gas accretion onto galaxies using the cosmological, hydrodynamical EAGLE simulation. We quantify the dependence of gas accretion on stellar mass, redshift, and environment, using halo mass and galaxy overdensity as environmental indicators. We find a strong suppression, by many orders of magnitude, of the gas accretion rate in dense environments, primarily for satellite galaxies. This suppression becomes stronger at lower redshift. However, the scatter in accretion rates is very large for satellites. This is (at least in part) due to the variation in halocentric radius, since gas accretion is more suppressed at smaller radii. Central galaxies are influenced less strongly by their environment and exhibit less scatter in their gas accretion rates. The star formation rates of both centrals and satellites show similar behaviour to their gas accretion rates. The relatively small differences between gas accretion and star formation rates demonstrate that galaxies generally exhaust their gas reservoir somewhat faster at higher stellar mass, lower redshift, and in denser environments. We conclude that the environmental suppression of gas accretion could directly result in the quenching of star formation.

We present measurements of star formation rates (SFRs) for dwarf galaxies (M*~10^8.5 Msun ) at z~1 using near-infrared slitless spectroscopy from the Hubble Space Telescope (HST) by targetting and measuring the luminosity of the H{\alpha} emission line. Our sample is derived from the Redshift One LDSS3 Emission Line Survey (ROLES), which used [O II]{\lambda} 3727 as a tracer of star formation to target very low stellar masses down to very low SFRs (~0.1 Msun yr^-1) at this epoch. Dust corrections are estimated using SED-fits and we find, by comparison with other studies using Balmer decrement dust corrections, that we require a smaller ratio between the gas phase and stellar extinction than the nominal Calzetti relation, in agreement with recent findings by other studies. By stacking the WFC3 spectra at the redshifts obtained from ground-based [O II] detections, we are able to push the WFC3 spectra to much lower SFRs and obtain the most complete spectroscopic measurement of the low mass end of the SSFR--mass relation to date. We measure a flatter low mass power law slope (-0.47 +/- 0.04) than found by other (shallower) H{\alpha}-selected samples (\approx -1), although still somewhat steeper than that predicted by the EAGLE simulation (-0.14 +/- 0.05), hinting at possible missing physics not modelled by EAGLE or remaining incompleteness for our H{\alpha} data.

We explore the parameter space of the semi-analytic galaxy formation model GALFORM, studying the constraints imposed by measurements of the galaxy stellar mass function (GSMF) and its evolution. We use the Bayesian Emulator method to quickly eliminate vast implausible volumes of the parameter space and zoom in on the most interesting regions, allowing us to identify a set of models that match the observational data within model uncertainties. We find that the GSMF strongly constrains parameters related to quiescent star formation in discs, stellar and AGN feedback and threshold for disc instabilities, but weakly restricts other parameters. Constraining the model using local data alone does not usually select models that match the evolution of the GSMF well. Nevertheless, we show that a small subset of models provides acceptable match to GSMF data out to redshift 1.5. We explore the physical significance of the parameters of these models, in particular exploring whether the model provides a better description if the mass loading of the galactic winds generated by starbursts ($\beta_{0,\text{burst}}$) and quiescent disks ($\beta_{0,\text{disc}}$) is different. Performing a principal component analysis of the plausible volume of the parameter space, we write a set of relations between parameters obeyed by plausible models with respect to GSMF evolution. We find that while $\beta_{0,\text{disc}}$ is strongly constrained by GSMF evolution data, constraints on $\beta_{0,\text{burst}}$ are weak. Although it is possible to find plausible models for which $\beta_{0,\text{burst}} = \beta_{0,\text{disc}}$, most plausible models have $\beta_{0,\text{burst}}>\beta_{0,\text{disc}}$, implying - for these - larger SN feedback efficiency at higher redshifts.

We use the hydrodynamical EAGLE simulation to study the magnitude and origin of the scatter in the stellar mass - halo mass relation for central galaxies. We separate cause and effect by correlating stellar masses in the baryonic simulation with halo properties in a matched dark matter only (DMO) simulation. The scatter in stellar mass increases with redshift and decreases with halo mass. At $z = 0.1$ it declines from 0.25 dex at $M_{200, \rm DMO} \approx 10^{11}$ M$_{\odot}$ to 0.12 dex at $M_{200, \rm DMO} \approx 10^{13}$ M$_{\odot}$, but the trend is weak above $10^{12}$ M$_{\odot}$. For $M_{200, \rm DMO} < 10^{12.5}$ M$_{\odot}$ up to 0.04 dex of the scatter is due to scatter in the halo concentration. At fixed halo mass, a larger stellar mass corresponds to a more concentrated halo. This is likely because higher concentrations imply earlier formation times and hence more time for accretion and star formation, and/or because feedback is less efficient in haloes with higher binding energies. The maximum circular velocity, $V_{\rm max, DMO}$, and binding energy are therefore more fundamental properties than halo mass, meaning that they are more accurate predictors of stellar mass, and we provide fitting formulae for their relations with stellar mass. However, concentration alone cannot explain the total scatter in the $M_{\rm star} - M_{200, \rm DMO}$ relation, and it does not explain the scatter in $M_{\rm star} -V_{\rm max, DMO}$. Halo spin, sphericity, triaxiality, substructure and environment are also not responsible for the remaining scatter, which thus could be due to more complex halo properties or non-linear/stochastic baryonic effects.

Galaxies fall into two clearly distinct types: `blue-sequence' galaxies that are rapidly forming young stars, and `red-sequence' galaxies in which star formation has almost completely ceased. Most galaxies more massive than $3\times10^{10} M_\odot$ follow the red-sequence while less massive central galaxies lie on the blue sequence. We show that these sequences are created by a competition between star formation-driven outflows and gas accretion on to the supermassive black hole at the galaxy's center. We develop a simple analytic model for this interaction. In galaxies less massive than $3\times10^{10} M_\odot$, young stars and supernovae drive a high entropy outflow that is more buoyant that any diffuse corona. The outflow balances the rate of gas inflow, preventing high gas densities building up in the central regions. More massive galaxies, however, are surrounded by a hot corona. We argue that above a halo mass of $\sim 10^{12} M_\odot$, the supernova-driven outflow is no longer buoyant and star formation is unable to prevent the build up of gas in the central regions. This triggers a strongly non-linear response from the black hole. Its accretion rate rises rapidly, heating the galaxy's corona, disrupting the incoming supply of cool gas and starving the galaxy of the fuel for star formation. The host galaxy makes a transition to the red sequence, and further growth predominantly occurs through galaxy mergers. We show that the analytic model provides a good description of galaxy evolution in the EAGLE hydrodynamic simulations, and demonstrate that, so long as star formation-driven outflows are present, the transition mass scale is almost independent of subgrid parameter choice. The transition mass disappears entirely, however, if star formation driven outflows are absent.

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

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.

We examine the properties of barred disc galaxies in a LCDM cosmological hydrodynamical simulation from the EAGLE project. Our study follows the formation of 269 discs identified at z = 0 in the stellar mass range 10.6 < log Mstr /M < 11. These discs show a wide range of bar strengths, from unbarred discs to weak bars to strongly barred systems (= 20%). Bars in these systems develop after redshift = 1.3, on timescales that depend sen- sitively on the strength of the pattern. Strong bars develop relatively quickly (in a few Gyr, = 10 disc rotation periods) in systems that are disc dominated, gas poor, and have declining rotation curves. Weak bars develop more slowly in systems where the disc is less gravitation- ally important, and are still growing at z = 0. Unbarred galaxies are comparatively gas-rich discs whose rotation speeds do not exceed the maximum circular velocity of the halos they inhabit. Bar lengths compare favourably with observations, ranging from 0.2 to 0.8 times the radius containing 90% of the stars. Bars slow down remarkably quickly as they grow, causing the inner regions of the surrounding dark halo to expand. At z = 0 strong bars have corota- tion radii roughly ten times the bar length. Such slow bars are inconsistent with the few cases where pattern speeds have been measured or inferred observationally, a discrepancy that, if confirmed, might prove a challenge for disc galaxy formation in LCDM.

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.

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$\mu$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 $\mu$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$\%$.

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)

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]

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 $\Lambda$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.

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 $\Lambda$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$. $\Lambda$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.

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.

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.

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.

We study the origin of the stellar $\alpha$-element-to-iron abundance ratio, $[\alpha/\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}$, $[\alpha/\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 $\alpha$-enriched by massive stars. In the absence of feedback from active galactic nuclei (AGN), however, $[\alpha/\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 $\alpha$-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.

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.

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.

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.

To investigate what drives the reversal of the morphology-density relation at intermediate/high redshift, we present a multi-wavelength analysis of 27 dusty starburst galaxies in the massive cluster Cl 0024+17 at z = 0.4. We combine H-alpha dynamical maps from the VLT/FLAMES multi-IFU system with far-infrared imaging using Herschel SPIRE and millimetre spectroscopy from IRAM/NOEMA, in order to measure the dynamics, star formation rates and gas masses of this sample. Most galaxies appear to be rotationally supported, with a median ratio of rotational support to line-of-sight velocity dispersion v/sigma ~ 5 +/- 2, and specific angular momentum lambda_R = 0.83 +/- 0.06 - comparable to field spirals of a similar mass at this redshift. The star formation rates of 3 - 26 M_solar/yr and average 12 CO derived gas mass of 1 x 10^10 M_solar suggest gas depletion timescales of ~ 1Gyr (~ 0.25 of the cluster crossing time). We derive characteristic dust temperatures (mean T_dust = 26 +/- 1 K) consistent with local galaxies of similar far-infrared luminosity, suggesting that the low density gas is yet to be stripped. Taken together, these results suggest that these starbursts have only recently accreted from the field, with star formation rates likely enhanced due to the effects of ram pressure. In order to make the transition to cluster S0s these galaxies must lose ~ 40% of their specific angular momentum. We suggest this must occur > 1 Gyr later, after the molecular gas has been depleted and/or stripped, via multiple tidal interactions with other cluster members.

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 $\mu>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.

We present the stellar mass ($M_{*}$), and K-corrected $K$-band absolute magnitude ($M_{K}$) Tully-Fisher relations (TFRs) for sub-samples of the 584 galaxies spatially resolved in H$\alpha$ emission by the KMOS Redshift One Spectroscopic Survey (KROSS). We model the velocity field of each of the KROSS galaxies and extract a rotation velocity, $V_{80}$ at a radius equal to the major axis of an ellipse containing 80% of the total integrated H$\alpha$ flux. The large sample size of KROSS allowed us to select 210 galaxies with well measured rotation speeds. We extract from this sample a further 56 galaxies that are rotationally supported, using the stringent criterion $V_{80}/\sigma > 3$, where $\sigma$ is the flux weighted average velocity dispersion. We find the $M_{K}$ and $M_{*}$ TFRs for this sub-sample to be $M_{K} / \rm{mag}= (-7.3 \pm 0.9) \times [(\log(V_{80}/\rm{km\ s^{-1}})-2.25]- 23.4 \pm 0.2$ , and $\log(M_{*} / M_{\odot})= (4.7 \pm 0.4) \times [(\log(V_{80}/\rm{km\ s^{-1}}) - 2.25] + 10.0 \pm 0.3$, respectively. We find an evolution of the $M_{*}$ TFR zero-point of $-0.41 \pm 0.08$ dex over the last $\sim $8 billion years. However, we measure no evolution in the $M_{K}$ TFR zero-point over the same period. We conclude that rotationally supported galaxies of a given dynamical mass had less stellar mass at $z \sim 1$ than the present day, yet emitted the same amounts of $K$-band light. The ability of KROSS to differentiate, using integral field spectroscopy with KMOS, between those galaxies that are rotationally supported and those that are not explains why our findings are at odds with previous studies without the same capabilities.

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.

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.

The scaling of disk galaxy rotation velocity with baryonic mass (the "Baryonic Tully-Fisher" relation, BTF) has long confounded galaxy formation models. It is steeper than the M ~ V^3 scaling relating halo virial masses and circular velocities and its zero point implies that galaxies comprise a very small fraction of available baryons. Such low galaxy formation efficiencies may in principle be explained by winds driven by evolving stars, but the tightness of the BTF relation argues against the substantial scatter expected from such vigorous feedback mechanism. We use the APOSTLE/EAGLE simulations to show that the BTF relation is well reproduced in LCDM simulations that match the size and number of galaxies as a function of stellar mass. In such models, galaxy rotation velocities are proportional to halo virial velocity and the steep velocity-mass dependence results from the decline in galaxy formation efficiency with decreasing halo mass needed to reconcile the CDM halo mass function with the galaxy luminosity function. Despite the strong feedback, the scatter in the simulated BTF is smaller than observed, even when considering all simulated galaxies and not just rotationally-supported ones. The simulations predict that the BTF should become increasingly steep at the faint end, although the velocity scatter at fixed mass should remain small. Observed galaxies with rotation speeds below ~40 km/s seem to deviate from this prediction. We discuss observational biases and modeling uncertainties that may help to explain this disagreement in the context of LCDM models of dwarf galaxy formation.

The KMOS Redshift One Spectroscopic Survey (KROSS) is an ESO guaranteed time survey of 795 typical star-forming galaxies in the redshift range z=0.8-1.0 with the KMOS instrument on the VLT. In this paper we present resolved kinematics and star formation rates for 584 z~1 galaxies. This constitutes the largest near-infrared Integral Field Unit survey of galaxies at z~1 to date. We demonstrate the success of our selection criteria with 90% of our targets found to be Halpha emitters, of which 81% are spatially resolved. The fraction of the resolved KROSS sample with dynamics dominated by ordered rotation is found to be 83$\pm$5%. However, when compared with local samples these are turbulent discs with high gas to baryonic mass fractions, ~35%, and the majority are consistent with being marginally unstable (Toomre Q~1). There is no strong correlation between galaxy averaged velocity dispersion and the total star formation rate, suggesting that feedback from star formation is not the origin of the elevated turbulence. We postulate that it is the ubiquity of high (likely molecular) gas fractions and the associated gravitational instabilities that drive the elevated star-formation rates in these typical z~1 galaxies, leading to the ten-fold enhanced star-formation rate density. Finally, by comparing the gas masses obtained from inverting the star-formation law with the dynamical and stellar masses, we infer an average dark matter to total mass fraction within 2.2$r_e$ (9.5kpc) of 65$\pm$12%, in agreement with the results from hydrodynamic simulations of galaxy formation.

We present first results from the KMOS Redshift One Spectroscopic Survey (KROSS), an ongoing large kinematical survey of a thousand, z~1 star forming galaxies, with VLT KMOS. Out of the targeted galaxies (~500 so far), we detect and spatially resolve Ha emission in ~90% and 77% of the sample respectively. Based on the integrated Ha flux measurements and the spatially resolved maps we derive a median star formation rate (SFR) of ~7.0 Msun/yr and a median physical size of <r$_{\rm e}$> = 5.1kpc. We combine the inferred SFRs and effective radii measurements to derive the star formation surface densities ({\Sigma}SFR) and present a "resolved" version of the star formation main sequence (MS) that appears to hold at sub-galactic scales, with similar slope and scatter as the one inferred from galaxy integrated properties. Our data also yield a trend between {\Sigma}SFR and {\Delta}(sSFR) (distance from the MS) suggesting that galaxies with higher sSFR are characterised by denser star formation activity. Similarly, we find evidence for an anti-correlation between the gas phase metallicity (Z) and the {\Delta}(sSFR), suggesting a 0.2dex variation in the metal content of galaxies within the MS and significantly lower metallicities for galaxies above it. The origin of the observed trends between {\Sigma}SFR - {\Sigma}(sSFR) and Z - {\Delta}(sSFR) could be driven by an interplay between variations of the gas fraction or the star formation efficiency of the galaxies along and off the MS. To address this, follow-up observations of the our sample that will allow gas mass estimates are necessary.

We explore the co-evolution of the specific angular momentum of dark matter haloes and the cold baryons that comprise the galaxies within. We study over two thousand central galaxies within the reference cosmological hydrodynamical simulation of the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) project. We employ a methodology within which the evolutionary history of a system is specified by the time-evolving properties of the Lagrangian particles that define it at z=0. We find a strong correlation between the evolution of the specific angular momentum of today's stars (cold gas) and that of the inner (whole) dark matter halo they are associated with. This link is particularly strong for the stars formed before the epoch of maximum expansion and subsequent collapse of the central dark matter halo (turnaround). Spheroids are typically assembled primarily from stars formed prior to turnaround, and are therefore destined to suffer a net loss of angular momentum associated with the strong merging activity during the assembly of the inner dark matter halo. Stellar discs retain their specific angular momentum since they are comprised of stars formed mainly after turnaround, from gas that mostly preserves the high specific angular momentum it acquired by tidal torques during the linear growth of the halo. Since the specific angular momentum loss of the stars is tied to the galaxy's morphology today, it may be possible to use our results to predict, statistically, the assembly history of a halo given the morphology of the galaxy it hosts.

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

We present the first results from the KMOS AGN Survey at High redshift (KASHz), a VLT/KMOS integral-field spectroscopic survey of z>0.6 AGN. We present galaxy-integrated spectra of 89 X-ray AGN (Lx=10^42-10^45 erg/s), for which we observed [O III] (z=1.1-1.7) or Halpha emission (z=0.6-1.1). The targets have X-ray luminosities representative of the parent AGN population and we explore the emission-line luminosities as a function of X-ray luminosity. For the [O III] targets, ~50 per cent have ionised gas velocities indicative of gas that is dominated by outflows and/or highly turbulent material (i.e., overall line-widths >~600 km/s). The most luminous half (i.e., Lx>6x10^43 erg/s) have a >~2 times higher incidence of such velocities. On the basis of our results, we find no evidence that X-ray obscured AGN are more likely to host extreme kinematics than unobscured AGN. Our KASHz sample has a distribution of gas velocities that is consistent with a luminosity-matched sample of z<0.4 AGN. This implies little evolution in the prevalence of ionised outflows, for a fixed AGN luminosity, despite an order-of-magnitude decrease in average star-formation rates over this redshift range. Furthermore, we compare our Halpha targets to a redshift-matched sample of star-forming galaxies and despite a similar distribution of Halpha luminosities and likely star-formation rates, we find extreme ionised gas velocities are up to ~10x more prevalent in the AGN-host galaxies. Our results reveal a high prevalence of extreme ionised gas velocities in high-luminosity X-ray AGN and imply that the most powerful ionised outflows in high-redshift galaxies are driven by AGN activity.

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 \Sigma_HI > 1 M_Sun pc^-2 (by a factor of <~ 2 compared to data from the Bluedisk survey). In the outer region (\Sigma_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)

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" ({\Lambda}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 {\Lambda}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.

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.

We present the evolution of galaxy sizes, from redshift 2 to 0, for actively star forming and passive galaxies in the cosmological hydrodynamical 1003 cMpc3 simulation of the EAGLE project. We find that the sizes increase with stellar mass , but that the relation weakens with increasing redshift. Separating galaxies by their star formation activity, we find that passive galaxies are typically smaller than active galaxies at fixed stellar mass. These trends are consistent with those found in observations and the level of agreement between the predicted and observed size - mass relation is of order 0.1 dex for z < 1 and 0.2-0.3 dex from redshift 1 to 2. We use the simulation to compare the evolution of individual galaxies to that of the population as a whole. While the evolution of the size-stellar mass relation for active galaxies provides a good proxy for the evolution of individual galaxies, the evolution of individual passive galaxies is not well represented by the observed size - mass relation due to the evolving number density of passive galaxies. Observations of z \approx 2 galaxies have revealed an abundance of massive red compact galaxies, that depletes below z \approx 1. We find that a similar population forms naturally in the simulation. Comparing these galaxies to their z = 0 descendants, we find that all compact galaxies grow in size due to the high-redshift stars migrating outwards. Approximately 60% of the compact galaxies increase in size further due to renewed star formation and/or mergers.

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, $\epsilon_{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.

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.

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

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.

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 $\epsilon = 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 $\sigma \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.

We examine the circular velocity profiles of galaxies in {\Lambda}CDM cosmological hydrodynamical simulations from the EAGLE and LOCAL GROUPS projects and compare them with a compilation of observed rotation curves of galaxies spanning a wide range in mass. The shape of the circular velocity profiles of simulated galaxies varies systematically as a function of galaxy mass, but shows remarkably little variation at fixed maximum circular velocity. This is especially true for low-mass dark matter-dominated systems, reflecting the expected similarity of the underlying cold dark matter haloes. This is at odds with observed dwarf galaxies, which show a large diversity of rotation curve shapes, even at fixed maximum rotation speed. Some dwarfs have rotation curves that agree well with simulations, others do not. The latter are systems where the inferred mass enclosed in the inner regions is much lower than expected for cold dark matter haloes and include many galaxies where previous work claims the presence of a constant density "core". The "cusp vs core" issue is thus better characterized as an "inner mass deficit" problem than as a density slope mismatch. For several galaxies the magnitude of this inner mass deficit is well in excess of that reported in recent simulations where cores result from baryon-induced fluctuations in the gravitational potential. We conclude that one or more of the following statements must be true: (i) the dark matter is more complex than envisaged by any current model; (ii) current simulations fail to reproduce the effects of baryons on the inner regions of dwarf galaxies; and/or (iii) the mass profiles of "inner mass deficit" galaxies inferred from kinematic data are incorrect.

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, $\Omega_{\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.

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.

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, $\rho_*(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.

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 $\Lambda$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.

We present a comparison of 14 galaxy formation models: 12 different semi-analytical models and 2 halo-occupation distribution models for galaxy formation based upon the same cosmological simulation and merger tree information derived from it. The participating codes have proven to be very successful in their own right but they have all been calibrated independently using various observational data sets, stellar models, and merger trees. In this paper we apply them without recalibration and this leads to a wide variety of predictions for the stellar mass function, specific star formation rates, stellar-to- halo mass ratios, and the abundance of orphan galaxies. The scatter is much larger than seen in previous comparison studies primarily because the codes have been used outside of their native environment within which they are well tested and calibrated. The purpose of the `nIFTy comparison of galaxy formation models' is to bring together as many different galaxy formation modellers as possible and to investigate a common approach to model calibration. This paper provides a unified description for all participating models and presents the initial, uncalibrated comparison as a baseline for our future studies where we will develop a common calibration framework and address the extent to which that reduces the scatter in the model predictions seen here.

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

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.

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.

We present adaptive optics-assisted integral field spectroscopy around the Ha or Hb lines of 12 gravitationally lensed galaxies obtained with VLT/SINFONI, Keck/OSIRIS and Gemini/NIFS. We combine these data with previous observations and investigate the dynamics and star formation properties of 17 lensed galaxies at z = 1-4. Thanks to gravitational magnification of 1.4 - 90x by foreground clusters, effective spatial resolutions of 40 - 700 pc are achieved. The magnification also allows us to probe lower star formation rates and stellar masses than unlensed samples; our target galaxies feature dust-corrected SFRs derived from Ha or Hb emission of 0.8 - 40Msol/yr, and stellar masses M* ~ 4e8 - 6e10 Msol. All of the galaxies have velocity gradients, with 59% consistent with being rotating discs and a likely merger fraction of 29%, with the remaining 12% classed as 'undetermined.' We extract 50 star-forming clumps with sizes in the range 60pc - 1kpc from the Ha (or Hb) maps, and find that their surface brightnesses and their characteristic luminosities evolve to higher luminosities with redshift. We show that this evolution can be described by fragmentation on larger scales in gas-rich discs, and is likely to be driven by evolving gas fractions.

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.

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.

We present MUSE observations of ESO137-001, a spiral galaxy infalling towards the center of the massive Norma cluster at z~0.0162. During the high-velocity encounter of ESO137-001 with the intracluster medium, a dramatic ram-pressure stripping event gives rise to an extended gaseous tail, traced by our MUSE observations to >30 kpc from the galaxy center. By studying the H-alpha surface brightness and kinematics in tandem with the stellar velocity field, we conclude that ram pressure has completely removed the interstellar medium from the outer disk, while the primary tail is still fed by gas from the inner regions. Gravitational interactions do not appear to be a primary mechanism for gas removal. The stripped gas retains the imprint of the disk rotational velocity to ~20 kpc downstream, without a significant gradient along the tail, which suggests that ESO137-001 is fast moving along a radial orbit in the plane of the sky. Conversely, beyond ~20 kpc, a greater degree of turbulence is seen, with velocity dispersion up to >100 km/s. For a model-dependent infall velocity of ~3000 km/s, we conclude that the transition from laminar to turbulent flow in the tail occurs on timescales >6.5 Myr. Our work demonstrates the terrific potential of MUSE for detailed studies of how ram-pressure stripping operates on small scales, providing a deep understanding of how galaxies interact with the dense plasma of the cluster environment.

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.

We have observed a sample of typical z=1 star forming galaxies, selected from the HiZELS survey, with the new KMOS near-infrared, multi-IFU instrument on the VLT, in order to obtain their dynamics and metallicity gradients. The majority of our galaxies have a metallicity gradient consistent with being flat or negative (i.e. higher metallicity cores than outskirts). Intriguingly, we find a trend between metallicity gradient and specific star formation rate (sSFR), such that galaxies with a high sSFR tend to have relatively metal-poor centres, a result which is strengthened when combined with datasets from the literature. This result appears to explain the discrepancies reported between different high redshift studies and varying claims for evolution. From a galaxy evolution perspective, the trend we see would mean that a galaxy's sSFR is governed by the amount of metal poor gas that can be funnelled into its core, triggered either by merging or through efficient accretion. In fact merging may play a significant role as it is the starburst galaxies at all epochs, which have the more positive metallicity gradients. Our results may help to explain the origin of the fundamental metallicity relation, in which galaxies at a fixed mass are observed to have lower metallicities at higher star formation rates, especially if the metallicity is measured in an aperture encompassing only the central regions of the galaxy. Finally, we note that this study demonstrates the power of KMOS as an efficient instrument for large scale resolved galaxy surveys.

We report the detection at 850um of the central source in SSA22-LAB1, the archetypal Lyman-alpha Blob (LAB), a 100kpc-scale radio-quiet emission-line nebula at z=3.1. The flux density of the source, $S_{850}=4.6\pm1.1$mJy implies the presence of a galaxy, or group of galaxies, with a total luminosity of $L_{\rm IR}\approx10^{12}L_\odot$. The position of an active source at the center of a ~50kpc-radius ring of linearly polarized Ly-alpha emission detected by Hayes et al. (2011) suggests that the central source is leaking Ly-alpha photons preferentially in the plane of the sky, which undergo scattering in HI clouds at large galactocentric radius. The Ly-alpha morphology around the submillimeter detection is reminiscent of biconical outflow, and the average Ly-alpha line profiles of the two `lobes' are dominated by a red peak, expected for a resonant line emerging from a medium with a bulk velocity gradient that is outflowing relative to the line center. Taken together, these observations provide compelling evidence that the central active galaxy (or galaxies) is responsible for a large fraction of the extended Ly-alpha emission and morphology. Less clear is the history of the cold gas in the circumgalactic medium being traced by Ly-alpha: is it mainly pristine material accreting into the halo that has not yet been processed through an interstellar medium (ISM), now being blown back as it encounters an outflow, or does it mainly comprise gas that has been swept-up within the ISM and expelled from the galaxy?

Extending over three Hubble Space Telescope (HST) cycles, the Hubble Frontier Fields (HFF) initiative constitutes the largest commitment ever of HST time to the exploration of the distant Universe via gravitational lensing by massive galaxy clusters. We here present models of the mass distribution in the six HFF cluster lenses, derived from a joint strong- and weak-lensing analysis anchored by a total of 88 multiple-image systems identified in existing HST data. The resulting maps of the projected mass distribution and of the gravitational magnification effectively calibrate the HFF clusters as gravitational telescopes. Allowing the computation of search areas in the source plane, these maps are provided to the community to facilitate the exploitation of forthcoming HFF data for quantitative studies of the gravitationally lensed population of background galaxies. Our models of the gravitational magnification afforded by the HFF clusters allow us to quantify the lensing-induced boost in sensitivity over blank-field observations and predict that galaxies at $z>10$ and as faint as m(AB)=32 will be detectable, up to 2 magnitudes fainter than the limit of the Hubble Ultra Deep Field.

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.

We present the data release of the Gemini-South GMOS spectroscopy in the fields of 11 galaxy groups at $0.8<z<1$, within the COSMOS field. This forms the basis of the Galaxy Environment Evolution Collaboration 2 (GEEC2) project to study galaxy evolution in haloes with $M\sim 10^{13}M_\odot$ across cosmic time. The final sample includes $162$ spectroscopically--confirmed members with $R<24.75$, and is $>50$ per cent complete for galaxies within the virial radius, and with stellar mass $M_{\rm star}>10^{10.3}M_\odot$. Including galaxies with photometric redshifts we have an effective sample size of $\sim 400$ galaxies within the virial radii of these groups. We present group velocity dispersions, dynamical and stellar masses. Combining with the GCLASS sample of more massive clusters at the same redshift we find the total stellar mass is strongly correlated with the dynamical mass, with $\log{M_{200}}=1.20\left(\log{M_{\rm star}}-12\right)+14.07$. This stellar fraction of $~\sim 1$ per cent is lower than predicted by some halo occupation distribution models, though the weak dependence on halo mass is in good agreement. Most groups have an easily identifiable most massive galaxy (MMG) near the centre of the galaxy distribution, and we present the spectroscopic properties and surface brightness fits to these galaxies. The total stellar mass distribution in the groups, excluding the MMG, compares well with an NFW profile with concentration $4$, for galaxies beyond $\sim 0.2R_{200}$. This is more concentrated than the number density distribution, demonstrating that there is some mass segregation.

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.

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.

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.

Using a large sample of field and group galaxies drawn from the Pan-STARRS1 Medium-Deep Survey, we present the specific star formation rate (SSFR) - stellar mass (M*) relation, as well as the quiescent fraction versus M* relation in different environments. We confirm that the fraction of quiescent galaxies is strongly dependent on environment at a fixed stellar mass, but that the amplitude and the slope of the star-forming sequence is similar between the field and groups: in other words, the SSFR-density relation at a fixed stellar mass is primarily driven by the change in the star-forming and quiescent fractions between different environments rather than a global suppression in the star formation rate for the star-forming population. However, when we restrict our sample to the cluster-scale environments ($M>10^{14}M_{solar}$), we find a global reduction in the SSFR of the star forming sequence of $17\%$ at 4$\sigma$ confidence as opposed to its field counterpart. After removing the stellar mass dependence of the quiescent fraction seen in field galaxies, the excess in the quiescent fraction due to the environment quenching in groups and clusters is found to increase with stellar mass. We argue that these results are in favor of galaxy mergers to be the primary environment quenching mechanism operating in galaxy groups whereas strangulation is able to reproduce the observed trend in the environment quenching efficiency and stellar mass relation seen in clusters. Our results also suggest that the relative importance between mass quenching and environment quenching depends on stellar mass -- the mass quenching plays a dominant role in producing quiescent galaxies for more massive galaxies, while less massive galaxies are quenched mostly through the environmental effect, with the transition mass around $1-2\times10^{10}M_{solar}$ in the group/cluster environment. (abridged)

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{\delta} 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).

We present the spatially resolved H-alpha (Ha) dynamics of sixteen star-forming galaxies at z~0.81 using the new KMOS multi-object integral field spectrograph on the ESO VLT. These galaxies were selected using 1.18 um narrow-band imaging from the 10 deg^2 CFHT-HiZELS survey of the SA22hr field, are found in a ~4Mpc over-density of Ha emitters and likely reside in a group/intermediate environment, but not a cluster. We confirm and identify a rich group of star-forming galaxies at z=0.813+-0.003, with thirteen galaxies within 1000 km/s of each other, and 7 within a diameter of 3Mpc. All our galaxies are "typical" star-forming galaxies at their redshift, 0.8+-0.4 SFR*(z=0.8), spanning a range of specific star formation rate of sSFR=0.2-1.1 Gyr^-1 and have a median metallicity very close to solar of 12+log(O/H)=8.62+-0.06. We measure the spatially resolved Ha dynamics of the galaxies in our sample and show that thirteen out of sixteen galaxies can be described by rotating disks and use the data to derive inclination corrected rotation speeds of 50-275 km/s. The fraction of disks within our sample is 75+-8, consistent with previous results based on HST morphologies of Ha selected galaxies at z~1 and confirming that disks dominate the star formation rate density at z~1. Our Ha galaxies are well fitted by the z~1-2 Tully-Fisher relation, confirming the evolution seen in the zero-point. Apart from having, on average, higher stellar masses and lower sSFRs, our group galaxies at z=0.813 present the same mass-metallicity and TF relation as z~1 field galaxies, and are all disk galaxies.

We obtained Subaru FMOS observations of Halpha emitting galaxies selected from the HiZELS narrow-band survey, to investigate the relationship between stellar mass, metallicity and star-formation rate at z = 0.84 - 1.47, for comparison with the Fundamental Metallicity Relation seen at low redshift. Our findings demonstrate, for the first time with a homogeneously selected sample, that a relationship exists for typical star-forming galaxies at z = 1 - 1.5 and that it is surprisingly similar to that seen locally. Therefore, star-forming galaxies at z = 1 - 1.5 are no less metal abundant than galaxies of similar mass and star formation rate (SFR) at z = 0.1, contrary to claims from some earlier studies. We conclude that the bulk of the metal enrichment for this star-forming galaxy population takes place in the 4 Gyr before z = 1.5. We fit a new mass-metallicity-SFR plane to our data which is consistent with other high redshift studies. However, there is some evidence that the mass-metallicity component of this high redshift plane is flattened, at all SFR, compared with z = 0.1, suggesting that processes such as star-formation driven winds, thought to remove enriched gas from low mass halos, are yet to have as large an impact at this early epoch. The negative slope of the SFR-metallicity relation from this new plane is consistent with the picture that the elevation in the SFR of typical galaxies at z > 1 is fuelled by the inflow of metal-poor gas and not major merging.

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.

Chronos is our response to ESA's call for white papers to define the science for the future L2, L3 missions. Chronos targets the formation and evolution of galaxies, by collecting the deepest NIR spectroscopic data, from the formation of the first galaxies at z~10 to the peak of formation activity at z~1-3. The strong emission from the atmospheric background makes this type of survey impossible from a ground-based observatory. The spectra of galaxies represent the equivalent of a DNA fingerprint, containing information about the past history of star formation and chemical enrichment. The proposed survey will allow us to dissect the formation process of galaxies including the timescales of quenching triggered by star formation or AGN activity, the effect of environment, the role of infall/outflow processes, or the connection between the galaxies and their underlying dark matter haloes. To provide these data, the mission requires a 2.5m space telescope optimised for a campaign of very deep NIR spectroscopy. A combination of a high multiplex and very long integration times will result in the deepest, largest, high-quality spectroscopic dataset of galaxies from z=1 to 12, spanning the history of the Universe, from 400 million to 6 billion years after the big bang, i.e. covering the most active half of cosmic history.

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

Mechanical feedback via Active Galactic Nuclei (AGN) jets in the centres of galaxy groups and clusters is a crucial ingredient in current models of galaxy formation and cluster evolution. Jet feedback is believed to regulate gas cooling and thus star formation in the most massive galaxies, but a robust physical understanding of this feedback mode is currently lacking. The large collecting area, excellent spectral resolution and high spatial resolution of Athena+ will provide the breakthrough diagnostic ability necessary to develop this understanding, via: (1) the first kinematic measurements on relevant spatial scales of the hot gas in galaxy, group and cluster haloes as it absorbs the impact of AGN jets, and (2) vastly improved ability to map thermodynamic conditions on scales well-matched to the jets, lobes and gas disturbances produced by them. Athena+ will therefore determine for the first time how jet energy is dissipated and distributed in group and cluster gas, and how a feedback loop operates in group/cluster cores to regulate gas cooling and AGN fuelling. Athena+ will also establish firmly the cumulative impact of powerful radio galaxies on the evolution of baryons from the epoch of group/cluster formation to the present day.

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 (\epsilon_\rho) 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 (\tau) 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 \tau_q<1 Gyr.

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.

We use the HiZELS narrow-band H-alpha survey in combination with CANDELS, UKIDSS and WIRDS near-infrared imaging, to investigate the morphologies, merger rates and sizes of a sample of H-alpha emitting galaxies in the redshift range z=0.40 - 2.23, an epoch encompassing the rise to the peak of the star formation rate density. Merger rates are estimated from space- and ground-based imaging using the M20 coefficient. To account for the increase in the specific star-formation rate (sSFR) of the star forming `main-sequence' with redshift, we normalise the star-formation rate of galaxies at each epoch to the typical value derived from the H-alpha luminosity function. Once this trend in sSFR is removed we see no evidence for an increase in the number density of star-forming galaxies or the merger rate with redshift. We thus conclude that neither is the main driver of the enhanced star-formation rate density at z=1-2, with secular processes such as instabilities within efficiently fuelled, gas-rich discs or multiple minor mergers the most likely alternatives. However, we find that 40-50% of starburst galaxies, those with enhanced specific star formation at their epoch, are major mergers and this fraction is redshift independent. Finally, we find the surprising result that the typical size of a star-forming galaxy of a given mass does not evolve across the redshift range considered, suggesting a universal size-mass relation. Taken in combination, these results indicate a star-forming galaxy population that is statistically similar in physical size, merger rate and mass over the 6 Gyr covered in this study, despite the increase in typical sSFR.

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]

We present the results of a search for molecular gas emission from a star-forming galaxy at z = 4.9. The galaxy benefits from magnification of 22 +/- 5x due to strong gravitational lensing by the foreground cluster MS1358+62. We target the CO(5-4) emission at a known position and redshift from existing Hubble Space Telescope/ACS imaging and Gemini/NIFS [OII]3727 imaging spectroscopy, and obtain a tentative detection at the 4.3sigma level with a flux of 0.104 +/- 0.024Jkm/s. From the CO line luminosity and assuming a CO-to-H2 conversion factor alpha=2, we derive a gas mass M_gas ~ 1^{+1}_{-0.6} x 10^9 M_sun. Combined with the existing data, we derive a gas fraction Mgas/(Mgas + M*) = 0.59^{+0.11}_{-0.06}. The faint line flux of this galaxy highlights the difficulty of observing molecular gas in representative galaxies at this epoch, and suggests that routine detections of similar galaxies in the absence of gravitational lensing will remain challenging even with ALMA in full science operations.

We present HST/WFC3 narrowband imaging of the H-alpha emission in a sample of eight gravitationally-lensed galaxies at z = 1 - 1.5. The magnification caused by the foreground clusters enables us to obtain a median source plane spatial resolution of 360pc, as well as providing magnifications in flux ranging from ~10x to ~50x. This enables us to identify resolved star-forming HII regions at this epoch and therefore study their H-alpha luminosity distributions for comparisons with equivalent samples at z ~ 2 and in the local Universe. We find evolution in the both luminosity and surface brightness of HII regions with redshift. The distribution of clump properties can be quantified with an HII region luminosity function, which can be fit by a power law with an exponential break at some cut-off, and we find that the cut-off evolves with redshift. We therefore conclude that `clumpy' galaxies are seen at high redshift because of the evolution of the cut-off mass; the galaxies themselves follow similar scaling relations to those at z = 0, but their HII regions are larger and brighter and thus appear as clumps which dominate the morphology of the galaxy. A simple theoretical argument based on gas collapsing on scales of the Jeans mass in a marginally unstable disk shows that the clumpy morphologies of high-z galaxies are driven by the competing effects of higher gas fractions causing perturbations on larger scales, partially compensated by higher epicyclic frequencies which stabilise the disk.

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)

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.

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

The observed stellar mass function (SMF) is very different to the halo mass function predicted by Lambda-CDM, and it is widely accepted that this is due to energy feedback from supernovae and black holes. However, the strength and form of this feedback is not understood. In this paper, we use the phenomenological model GALFORM to explore how galaxy formation depends on the strength and halo mass dependence of feedback. We focus on 'expulsion' models in which the wind mass loading, beta, is proportional to 1/\vdisk^n, with n=0,1,2 and contrast these models with the successful Bower et al.\ 2008 model (B8W7). A crucial development is that our code explicitly accounts for the recapture of expelled gas as the system's halo mass (and thus gravitational potential) increases. We find that a model with modest wind speed but high mass loading matches the flat portion of the SMF. When combined with AGN feedback, the model provides a good description of the observed SMF above 10^9 h^-1 Msol. However, in the expulsion models, the brightest galaxies are assembled more recently than in B8W7, and the specific star formation rates of galaxies decrease strongly with decreasing stellar mass. The expulsion models also tend to have a cosmic star formation density that is dominated by lower mass galaxies at z=1-3, and dominated high mass galaxies at low redshift. These trends are in conflict with observational data, but the comparison highlights some deficiencies of the B8W7 model also. The experiments in this paper give us important physical insight to the impact of the feedback process on the formation histories of galaxies, but the strong mass dependence of feedback adopted in B8W7 still appears to provide the most promising description of the observed universe.

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.

We present three dimensional spectroscopy of eleven E+A galaxies, selected for their strong H-delta absorption but weak (or non-existent) [OII]3727 and H-alpha emission. This selection suggests that a recent burst of star-formation was triggered but subsequently abruptly ended. We probe the spatial and spectral properties of both the young (~1Gyr) and old (few Gyr) stellar populations. Using the H-delta equivalent widths we estimate that the burst masses must have been at least 10% by mass (Mburst~10^10Mo), which is also consistent with the star-formation history inferred from the broad-band SEDs. On average the A-stars cover ~33% of the galaxy image, extending over 2-15kpc^2, indicating that the characteristic E+A signature is a property of the galaxy as a whole and not due to a heterogeneous mixture of populations. In approximately half of the sample, we find that the A-stars, nebular emission, and continuum emission are not co-located, suggesting that the newest stars are forming in a different place than those that formed ~1Gyr ago, and that recent star-formation has occurred in regions distinct from the oldest stellar populations. At least ten of the galaxies (91%) have dynamics that class them as "fast rotators" with magnitudes and dynamics comparable to local ellipticals and S0's. We also find a correlation between the spatial extent of the A-stars and dynamics such that the fastest rotators tend to have the most compact A-star populations, providing new constraints on models that aim to explain the transformation of later type galaxies into early types. Finally, we show that there are no obvious differences between the line extents and kinematics of E+A galaxies detected in the radio (AGN) compared to non-radio sources, suggesting that AGN feedback does not play a dramatic role in defining their properties, or that its effects are short.

We present a new cluster detection algorithm designed for the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) survey but with generic application to any multiband data. The method makes no prior assumptions about the properties of clusters other than (a) the similarity in colour of cluster galaxies (the "red sequence") and (b) an enhanced projected surface density. The detector has three main steps: (i) it identifies cluster members by photometrically filtering the input catalogue to isolate galaxies in colour-magnitude space, (ii) a Voronoi diagram identifies regions of high surface density, (iii) galaxies are grouped into clusters with a Friends-of-Friends technique. Where multiple colours are available, we require systems to exhibit sequences in two colours. In this paper we present the algorithm and demonstrate it on two datasets. The first is a 7 square degree sample of the deep Sloan Digital Sky Survey equatorial stripe (Stripe 82), from which we detect 97 clusters with z<=0.6. Benefiting from deeper data, we are 100% complete in the maxBCG optically-selected cluster catalogue (based on shallower single epoch SDSS data) and find an additional 78 previously unidentified clusters. The second dataset is a mock Medium Deep Survey (MDS) Pan-STARRS catalogue, based on the Lambda-CDM model and a semi-analytic galaxy formation recipe. Knowledge of galaxy-halo memberships in the mock allows a quantification of algorithm performance. We detect 305 mock clusters in haloes with mass >10^13 solar masses at z<=0.6 and determine a spurious detection rate of <1%, consistent with tests on the Stripe 82 catalogue. The detector performs well in the recovery of model Lambda-CDM clusters. (abridged)

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.

Galaxy star formation rates (SFRs) are sensitive to the local environment; for example, the high-density regions at the cores of dense clusters are known to suppress star formation. It has been suggested that galaxy transformation occurs largely in groups, which are the intermediate step in density between field and cluster environments. In this paper, we use deep MIPS 24 micron observations of intermediate-redshift (0.3 < z < 0.55) group and field galaxies from the Group Environment and Evolution Collaboration (GEEC) subset of the Second Canadian Network for Observational Cosmology (CNOC2) survey to probe the moderate-density environment of groups, wherein the majority of galaxies are found. The completeness limit of our study is log(L_TIR (L_sun)) > 10.5, corresponding to SFR > 2.7 M_sun/yr. We find that the group and field galaxies have different distributions of morphologies and mass. However, individual group galaxies have star-forming properties comparable to those of field galaxies of similar mass and morphology; that is, the group environment does not appear to modify the properties of these galaxies directly. There is a relatively large number of massive early-type group spirals, along with E/S0 galaxies, that are forming stars above our detection limit. These galaxies account for the nearly comparable level of star-forming activity in groups as compared with the field, despite the differences in mass and morphology distributions between the two environments. The distribution of specific SFRs (SFR/M_*) is shifted to lower values in the groups, reflecting the fact that groups contain a higher proportion of massive and less active galaxies. Considering the distributions of morphology, mass, and SFR, the group members appear to lie between field and cluster galaxies in overall properties.

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.

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.

We present the results of a narrow-band (H2S1, 2.121/0.021um) imaging search with WFCAM/UKIRT for H-alpha emitters around several potential signposts of rare (10^-7-10^-8 Mpc^-3) over-dense regions at z=2.23: an over-density of QSOs (2QZ cluster), a powerful, high-redshift radio galaxy (HzRG), and a concentration of submillimetre galaxies (SMGs) and optically faint radio galaxies (OFRGs). In total, we detect 137 narrow-band emitter candidates down to emission-line fluxes of 0.5-1 x 10^-16 erg s^-1 cm^-2, across a total area of 0.56 sq. degrees (2.1 x 10^5 comoving Mpc at z=2.23) in these fields. The BzK colours of the emitters suggest that at least 80% of our sample are likely to be H-alpha emitters (HAEs) at z=2.23. This is one of the largest HAE samples known at z>2. We find modest (~3-sigma) local over-densities of emitters associated with all the three targets. In the 2QZ cluster field, the emitters show a striking filamentary structure connecting four of the z=2.23 QSOs extending over 30Mpc (comoving). In the HzRG and SMG/OFRG fields, the structures appear to be smaller and seen only in the vicinities of the targets. The K-band magnitudes and the H-alpha equivalent widths of the emitters are weakly correlated with the over-density of the emitters: emitters in over-dense region are more evolved systems compared to those in under-dense regions at z=2.23. We find several examples of extended HAEs in our target fields, including a striking example with a spatial extent of 7.5 arcsec (60 kpc at z=2.23) in the 2QZ field, suggesting that these are relatively common in high-density regions. We conclude that narrow-band H-alpha surveys are efficient routes to map over-dense regions at high-z and thus to understand the relation between the growth of galaxies and their surrounding large-scale structures.

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.

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.

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]

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

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.

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.

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

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

Using the sample from the \it Redshift One LDSS3 Emission line Survey \rm (ROLES), we probe the dependence of star formation rate (SFR) and specific star formation rate (sSFR) as a function of stellar mass $M_*$ and environment as defined by local galaxy density, in the CDFS field. Our spectroscopic sample consists of 312 galaxies with $K_{AB}<24$, corresponding to stellar mass $\log(M_*/M_{\sun})>8.5$, and with [OII] derived star-formation rates SFR$>0.3M_{\sun}/$yr, at $0.889\leq z \leq 1.149$. The results have been compared directly with the Sloan Digital Sky Survey Stripe 82 sample at $0.032\leq z \leq 0.05$. For star-forming galaxies, we confirm that there is little correlation between SFR and density at $z\sim 0$. However, for the lowest mass galaxies in our $z\sim 1$ sample, those with $\log(M_*/M_{\sun})<10$, we find that both the median SFR and specific SFR {\it increase} significantly with increasing local density. The "downsizing" trend for low mass galaxies to be quenched progressively later in time appears to be more pronounced in moderately overdense environments. Overall we find that the evolution of star-formation in galaxies is most strongly driven by their stellar mass, with local galaxy density playing a role that becomes increasingly important for lower mass galaxies.

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

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

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

A generic expectation for gas accreted by high mass haloes is that it is shock heated to the virial temperature of the halo. In low mass haloes, or at high redshift, however, the gas cooling rate is sufficiently rapid that an accretion shock is unlikely to form. Instead, gas can accrete directly into the centre of the halo in a `cold mode' of accretion. Although semi-analytic models have always made a clear distinction between hydrostatic and rapid cooling they have not made a distinction between whether or not an accretion shock forms. Starting from the well-established Galform code, we investigate the effect of explicitly accounting for cold mode accretion using the shock stability model of Birnboim & Dekel. When we modify the code so that there is no effective feedback from galaxy formation, we find that cold mode accretion is the dominant channel for feeding gas into the galaxies at high redshifts. However, this does not translate into a significant difference in the star formation history of the universe compared to the previous code. When effective feedback is included in the model, we find that the the cold mode is much less apparent because of the presence of gas ejected from the galaxy. Thus the inclusion of the additional cold mode physics makes little difference to basic results from earlier semi-analytic models which used a simpler treatment of gas accretion. For more sophisticated predictions of its consequences, we require a better understanding of how the cold mode delivers angular momentum to galaxies and how it interacts with outflows.

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

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.

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.

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]

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.

We report observations of Lyman Alpha Blob 1 (LAB1) in the SSA 22 protocluster region (z=3.09) with the integral-field spectrograph SAURON. We increased the signal-to-noise in the spectra by more than a factor three compared to our previous observations. This allows us to probe the structure of the LAB system in detail, examining its structure in the spatial and wavelength dimensions. We find that the emission from the system comes largely from five distinct blobs. Two of the emission regions are associated with Lyman Break Galaxies, while a third appears to be associated with a heavily obscured submillimeter galaxy. The fourth and fifth components do not appear to be associated with any galaxy despite the deep imaging that is available in this field. If we interpret wavelength shifts in the line centroid as velocity structure in the underlying gas, many of these emission systems show evidence of velocity shear. It remains difficult to distinguish between an underlying rotation of the gas and an outflow driven by the central object. We have examined all of the line profiles for evidence of strong absorption features. While several systems are better fitted by the inclusion of a weak absorption component, we do not see evidence for a large-scale coherent absorption feature such as that seen in LAB2.

(Abridged) We examine the radial entropy distribution and its scaling using 31 nearby galaxy clusters from the Representative XMM-Newton Cluster Structure Survey (REXCESS). The entropy profiles are robustly measured at least out to R_1000 in all systems and out to R_500 in 13 systems. Compared to theoretical expectations, the observed distributions show a radial and mass-dependent excess entropy that is greater and extends to larger radii in lower mass systems. At R_500, the mass dependence and entropy excess are both negligible within the uncertainties. Mirroring this behaviour, the scaling of gas entropy is shallower than self-similar in the inner regions, but steepens with radius, becoming consistent with self-similar at R_500. The dispersion in scaled entropy in the inner regions is linked to the presence of cool cores and dynamical activity; at larger radii the dispersion decreases by a factor of two and the dichotomy between subsamples disappears. Parameterising the profiles with a power law plus constant model, there are two peaks in central entropy K_0; however, we cannot distinguish between a bimodal or a left-skewed distribution. The outer slopes are correlated with system temperature; their distribution is unimodal with a median value of 0.98. Renormalising the dimensionless entropy profiles by the gas mass fraction profile f_gas(< R), leads to a remarkable reduction in the scatter, implying that gas mass fraction variations with radius and mass are the cause of the observed entropy properties. We discuss a tentative scenario to explain the behaviour of the entropy and gas mass fraction in the REXCESS sample, in which extra heating and merger mixing maintains an elevated central entropy level in the majority of the population, and a smaller fraction of systems develops a cool core.

Recent wide-field imaging observations of the X-ray luminous cluster RDCSJ1252.9-2927 at z=1.24 uncovered several galaxy groups that appear to be embedded in filamentary structure extending from the cluster core. We make a spectroscopic study of the galaxies in these groups using GMOS on Gemini-South and FORS2 on VLT with the aim of determining if these galaxies are physically associated to the cluster. We find that three groups contain galaxies at the cluster redshift and that they are probably bound to the cluster. This is the first confirmation of filamentary structure as traced by galaxy groups at z>1. We then use several spectral features in the FORS2 spectra to determine the star formation histories of group galaxies. We find a population of relatively red star-forming galaxies in the groups that are absent from the cluster core. While similarly red star forming galaxies can also be found in the field, the average strength of the hd line is systematically weaker in group galaxies. Interestingly, these groups at z=1.2 are in an environment in which the on-going build-up of red sequence is happening. The unusual line strengths can be explained by star formation that is heavily obscured by dust. We hypothesize that galaxy-galaxy interactions, which is more efficient in the group environment, is the mechanism that drives these dust obscured star formation. The hypothesis can be tested by obtaining spectral observations in the near-IR, high resolution imaging observations and observations in the mid-IR.

X-ray properties of galaxy groups can unlock some of the most challenging research topics in modern extragalactic astronomy: the growth of structure and its influence on galaxy formation. Only with the advent of the Chandra and XMM facilities have X-ray observations reached the depths required to address these questions in a satisfactory manner. Here we present an X-ray imaging study of two patches from the CNOC2 spectroscopic galaxy survey using combined Chandra and XMM data. A state of the art extended source finding algorithm has been applied, and the resultant source catalog, including redshifts from a spectroscopic follow-up program, is presented. The total number of spectroscopically identified groups is 25 spanning a redshift range 0.04-0.79. Approximately 50% of CNOC2 spectroscopically selected groups in the deeper X-ray (RA14h) field are likely X-ray detections, compared to 20% in the shallower (RA21h) field. Statistical modeling shows that this is consistent with expectations, assuming an expected evolution of the Lx-M relation. A significant detection of a stacked shear signal for both spectroscopic and X-ray groups indicates that both samples contain real groups of about the expected mass. We conclude that the current area and depth of X-ray and spectroscopic facilities provide a unique window of opportunity at z~0.4 to test the X-ray appearance of galaxy groups selected in various ways. There is at present no evidence that the correlation between X-ray luminosity and velocity dispersion evolves significantly with redshift, which implies that catalogs based on either method can be fairly compared and modeled.

My aim in this talk is to make clear that there are two sides to galaxy formation: the properties of the galaxies themselves, and the properties of the material that is left over from the galaxy formation process. To date, galaxy formation studies have focused on correctly predicting the properties of galaxies, and I will review the tremendous level of success in this area. However, these models usually ignore the ``flip side'' of galaxy formation: the intergalactic medium and the intra-group/intra-cluster medium (ICM). Yet, Chandra and XMM have given us a good view of the ICM and their results present an equally important challenge for theoretical models. I will show that this challenge is far from easy to meet, but describe the Bower et al 2008 model of galaxy formation which successfully combines both sides of the observational constraints.

We present a detailed study of the spatially resolved kinematics, star-formation and stellar mass in a highly amplified galaxy at z=4.92 behind the lensing cluster MS1358+62. We use the observed optical, near- and mid-infrared imaging from HST ACS & NICMOS and Spitzer IRAC to derive the stellar mass and the Gemini/NIFS IFU to investigate the velocity structure of the galaxy from the nebular [OII] emission. Using a detailed gravitational lens model, we account for lensing amplification factor 12.+/-2.0 and find that this intrinsically L* galaxy has a stellar mass of M*=7+/-2x10^8Mo, a dynamical mass of Mdyn=3+/-1x10^9csc^2(i)Mo (within of 2kpc) and a star-formation rate of 42+/-8Mo/yr. The source-plane UV/optical morphology of this galaxy is dominated by five discrete star-forming regions. Exploiting the dynamical information we derive masses for individual star-forming regions of Mcl~10^(8-9)Mo with sizes of ~200pc. We find that, at a fixed size, the star-formation rate density within these HII regions is approximately two orders of magnitude greater than those observed in local spiral/starburst galaxies, but consistent with the most massive HII regions in the local Universe such as 30Dor. Finally, we compare the spatially resolved nebular emission line velocity with the Ly-alpha and UV-ISM lines and find that this galaxy is surrounded by a galactic scale outflow in which the Ly-alpha appears redshifted by ~150km/s and the UV-ISM lines blue-shifted by -200km/s from the (systemic) nebular emission. The velocity structure of the outflow mirrors that of the nebular emission suggesting the outflow is young (~15Myr), and has yet to burst out of the system. Taken together, these results suggest that this young galaxy is undergoing its first major epoch of mass assembly.

We analyse a 100-ks Chandra observation of the powerful radio galaxy, 4C60.07 at z=3.79. We identify extended X-ray emission with Lx~10^45 erg/s across a ~90-kpc region around the radio galaxy. The energetics of this X-ray halo and its morphological similarity to the radio emission from the galaxy suggest that it arises from inverse Compton (IC) scattering, by relativistic electrons in the radio jets, of Cosmic Microwave Background photons and potentially far-infrared photons from the dusty starbursts around this galaxy. The X-ray emission has a similar extent and morphology to the Ly-alpha halo around the galaxy, suggesting that it may be ionising this halo. Indeed we find that the GHz-radio and X-ray and Ly-alpha luminosities of the halo around 4C60.07 are identical to those of 4C41.17 (also at z=3.8) implying that these three components are linked by a single physical process. This is only the second example of highly-extended IC emission known at z>3, but it underlines the potential importance of IC emission in the formation of the most massive galaxies at high redshifts. In addition, we detect two X-ray luminous active galactic nuclei (AGN) within ~30kpc of the radio galaxy. These two companion AGN imply that the radio and starburst activity in the radio galaxy is triggered through multiple mergers of massive progenitors on a short timescale, ~100Myrs. These discoveries demonstrate the wealth of information which sensitive X-ray observations can yield into the formation of massive galaxies at high redshifts.

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.

We present the results of a 400ks Chandra survey of 29 extended Ly-alpha emitting nebulae (Ly-alpha Blobs, LABs) in the z=3.09 proto-cluster in the SSA22 field. We detect luminous X-ray counterparts in five LABs, implying a large fraction of active galactic nuclei (AGN) in LABs, f_AGN = 17% down to L_2-32keV ~ 10^44 erg/s. All of the AGN appear to be heavily obscured, with spectral indices implying obscuring column densities of N_H > 10^23 cm^-2. The AGN fraction should be considered a lower limit, since several more LABs not detected with Chandra show AGN signatures in their mid-infrared emission. We show that the UV luminosities of the AGN are easily capable of powering the extended Ly-alpha emission via photo-ionization alone. When combined with the UV flux from a starburst component, and energy deposited by mechanical feedback, we demonstrate that 'heating' by a central source, rather than gravitational cooling is the most likely power source of LABs. We argue that all LABs could be powered in this manner, but that the luminous host galaxies are often just below the sensitivity limits of current instrumentation, or are heavily obscured. No individual LABs show evidence for extended X-ray emission, and a stack equivalent to a >9Ms exposure of an average LAB also yields no statistical detection of a diffuse X-ray component. The resulting diffuse X-ray/Ly-alpha luminosity limit implies there is no hot (T>10^7 K) gas component in these halos, and also rules out inverse Compton scattering of cosmic microwave background photons, or local far-infrared photons, as a viable power source for LABs.

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

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.

We use a K-selected (22.5 < K_AB < 24.0) sample of dwarf galaxies (8.4 < log(M*/Msun) < 10) at 0.89<z<1.15 in the Chandra Deep Field South (CDFS) to measure their contribution to the global star-formation rate density (SFRD), as inferred from their [OII] flux. By comparing with [OII]-based studies of higher stellar mass galaxies, we robustly measure a turnover in the [OII] luminosity density at a stellar mass of M~10^10 Msun. By comparison with the [OII]-based SFRD measured from the Sloan Digital Sky Survey we confirm that, while the SFRD of the lowest-mass galaxies changes very little with time, the SFRD of more massive galaxies evolves strongly, such that they dominate the SFRD at z = 1.

(Abridged) The thermal state of the intracluster medium results from a competition between gas cooling and heating. The heating comes from two distinct sources: gravitational heating from the collapse of the dark matter halo and thermal input from galaxy/black hole formation. However, a long standing problem has been that cosmological simulations based on smoothed particle hydrodynamics (SPH) and Eulerian mesh codes predict different results even when cooling and galaxy/black hole heating are switched off. Clusters formed in SPH simulations show near powerlaw entropy profiles, while those formed in mesh simulations develop a core and do not allow gas to reach such low entropies. Since the cooling rate is closely connected to the minimum entropy of the gas, the differences are of potentially key importance. In this paper, we investigate the origin of this discrepancy. By comparing simulations run using the GADGET-2 SPH code and the FLASH adaptive Eulerian mesh code, we show that the discrepancy arises during the idealised merger of two clusters. The difference is not sensitive to the resolution of our simulations, nor is it is due differences in the gravity solvers, Galilean non-invariance of the mesh code, or an effect of unsuitable artificial viscosity in the SPH code. Instead, we find that the difference is inherent to the treatment of eddies and fluid instabilities. These are suppressed in the SPH simulations, while the cluster mergers generate strong vortices in the mesh simulations that efficiently mix the fluid and erase the low entropy gas. Consequently, particles in the SPH simulations retain a close connection to their initial entropy, while this connection is much weaker in the mesh simulations. We discuss the potentially profound implications of these results.

We explore the clustering properties of high redshift dark matter halos, focusing on halos massive enough to host early generations of stars or galaxies at redshift 10 and greater. Halos are extracted from an array of dark matter simulations able to resolve down to the "mini-halo" mass scale at redshifts as high as 30, thus encompassing the expected full mass range of halos capable of hosting luminous objects and sources of reionization. Halo clustering on large-scales agrees with the Sheth, Mo & Tormen halo bias relation within all our simulations, greatly extending the regime where large-scale clustering is confirmed to be "universal" at the 10-20% level (which means, for example, that 3sigma halos of cluster mass at z=0 have the same large-scale bias with respect to the mass distribution as 3sigma halos of galaxy mass at z=10). However, on small-scales, the clustering of our massive halos (> ~10^9 Msun/h) at these high redshifts is stronger than expected from comparisons with small-scale halo clustering extrapolated from lower redshifts. This implies "non-universality" in the scale-dependence of halo clustering, at least for the commonly used parameterizations of the scale-dependence of bias that we consider. We provide a fit for the scale-dependence of bias in our results. This study provides a basis for using extraordinarily high redshift galaxies (redshift ~10) as a probe of cosmology and galaxy formation at its earliest stages. We show also that mass and halo kinematics are strongly affected by finite simulation volumes. This suggests the potential for adverse affects on gas dynamics in hydrodynamic simulations of limited volumes, such as is typical in simulations of the formation of the "first stars", though further study is warranted.

The low redshift Universe (z<~0.5) is not a dull place. Processes leading to the suppression of star formation and morphological transformation are prevalent: this is particularly evident in the dramatic upturn in the fraction of S0-type galaxies in clusters. However, until now, the process and environment of formation has remained unidentified. We present a HST-based morphological analysis of galaxies in the redshift-space selected group and field environments at z~0.4. Groups contain a much higher fraction of S0s at fixed luminosity than the lower density field, with >99.999% confidence. Indeed the S0 fraction in groups is at least as high as in z~0.4 clusters and X-ray selected groups, which have more luminous Intra Group Medium (IGM). An 97% confident excess of S0s at >=0.3Mpc from the group centre at fixed luminosity, tells us that formation is not restricted to, and possibly even avoids, the group cores. Interactions with a bright X-ray emitting IGM cannot be important for the formation of the majority of S0s in the Universe. In contrast to S0s, the fraction of elliptical galaxies in groups at fixed luminosity is similar to the field, whilst the brightest ellipticals are strongly enhanced towards the group centres (>99.999% confidence within 0.3Mpc). We conclude that the group and sub-group environments must be dominant for the formation of S0 galaxies, and that minor mergers, galaxy harassment and tidal interactions are the most likely responsible mechanisms. This has implications not only for the inferred pre-processing of cluster galaxies, but also for the global morphological and star formation budget of galaxies: as hierarchical clustering progresses, more galaxies will be subject to these transformations as they enter the group environment.

We have used the VIMOS IFU to map the properties of the Seyfert 1.9 galaxy LEDA 135736. These maps reveal a number of interesting features including: an Extended Narrow Line Region detectable out to 9 kpc, an area of intense star formation located at a projected distance of 12 kpc from the centre, an elliptical companion galaxy, and kinematic features, aligned along the long-axis of the ENLR, that are consistent with radio jet-driven mass outflow. We propose that the ENLR results from extra-planar gas ionized by the AGN, and that the AGN in turn might be triggered by interaction with the companion galaxy, which can also explain the burst of star formation and morphological features. Only about two percent of the ENLR's kinetic energy is in the mass outflow. We infer from this that the bulk of mechanical energy imparted by the jet is used to heat this gas.

Only ~10% of baryons in the universe are in the form of stars, yet most models of luminous structure formation have concentrated on the properties of the luminous stellar matter. In this paper we focus on the "flip side" of galaxy formation and investigate the properties of the material that is not presently locked up in galaxies. This "by-product" of galaxy formation can be observed as an X-ray emitting plasma (the intracluster medium, hereafter ICM) in groups and clusters, and we present a version of the Durham semi-analytic galaxy formation model GALFORM that allows us to investigate the properties of the ICM. As we would expect on the basis of gravitational scaling arguments, the previous model (presented in Bower et al. 2006) fails to reproduce even the most basic observed properties of the ICM; however, we present a simple modification to the model to allow for heat input into the ICM from the AGN "radio mode" feedback. This heating acts to expel gas from the X-ray luminous central regions of the host halo. With this modification, the model reproduces the observed gas mass fractions and luminosity-temperature relation of groups and clusters. Introducing the heating process into the model requires changes to a number of model parameters in order to retain a good match to the observed galaxy properties. With the revised parameters, the best fitting luminosity function is comparable to that presented in Bower et al. (2006). The new model makes a fundamental step forward, providing a unified model of galaxy and cluster ICM formation. However, the detailed comparison with the data is not completely satisfactory, and we highlight key areas for improvement.

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.

We present a study of the structural and scaling properties of the gas distributions in the intracluster medium (ICM) of 31 nearby (z < 0.2) clusters observed with XMM-Newton, which together comprise the Representative XMM-Newton Cluster Structure Survey (REXCESS). In contrast to previous studies, this sample is unbiased with respect to cluster dynamical state, and it fully samples the cluster X-ray luminosity function. The clusters cover a temperature range of 2.0 -- 8.5 keV and possess a variety of morphologies. The sampling strategy allows us to compare clusters with a wide range of central cooling times on an equal footing. We present non-parametric gas-density profiles out to distances ranging between 0.8 R_500 and 1.5 R_500. The central gas densities differ greatly from system to system, with no clear correlation with system temperature. At intermediate radii the scaled density profiles show much less scatter, with a clear dependence on system temperature, consistent with the presence of an entropy excess as suggested in previous literature. However, at large scaled radii this dependence becomes weaker: clusters with kT > 3 keV scale self-similarly, with no temperature dependence of gas-density normalisation. We find some evidence of a correlation between dynamical state and outer gas density slope, and between dynamical state and both central gas normalisation and cooling time. We find no evidence of a significant bimodality in the distributions of central density, density gradient, or cooling time. Finally, we present the gas mass-temperature relation for the REXCESS sample, which is consistent with the expectation of self-similar scaling modified by the presence of an entropy excess in the inner regions of the cluster, and has a logarithmic intrinsic scatter of ~10%.

We present quantitative morphology measurements of a sample of optically selected group galaxies at 0.3 < z < 0.55 using the Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) and the GIM2D surface brightness--fitting software package. The group sample is derived from the Canadian Network for Observational Cosmology Field Redshift survey (CNOC2) and follow-up Magellan spectroscopy. We compare these measurements to a similarly selected group sample from the Millennium Galaxy Catalogue (MGC) at 0.05 < z < 0.12. We find that, at both epochs, the group and field fractional bulge luminosity (B/T) distributions differ significantly, with the dominant difference being a deficit of disk--dominated (B/T < 0.2) galaxies in the group samples. At fixed luminosity, z=0.4 groups have ~ 5.5 +/- 2 % fewer disk--dominated galaxies than the field, while by z=0.1 this difference has increased to ~ 19 +/- 6 %. Despite the morphological evolution we see no evidence that the group environment is actively perturbing or otherwise affecting the entire existing disk population. At both redshifts, the disks of group galaxies have similar scaling relations and show similar median asymmetries as the disks of field galaxies. We do find evidence that the fraction of highly asymmetric, bulge--dominated galaxies is 6 +/- 3 % higher in groups than in the field, suggesting there may be enhanced merging in group environments. We replicate our group samples at z=0.4 and z=0 using the semi-analytic galaxy catalogues of Bower et al (2006). This model accurately reproduces the B/T distributions of the group and field at z=0.1. However, the model does not reproduce our finding that the deficit of disks in groups has increased significantly since z=0.4.

Luminous red galaxies (LRGs) are much rarer and more massive than L* galaxies. Coupled with their extreme colours, LRGs therefore provide a demanding testing ground for the physics of massive galaxy formation. We present the first self-consistent predictions for the abundance and properties of LRGs in hierarchical structure formation models. We test two published models which use quite different mechanisms to suppress the formation of massive galaxies: the Bower et al. (2006) model, which invokes ``AGN-feedback'' to prevent gas from cooling in massive haloes, and the Baugh et al. (2005) model which relies upon a ``superwind'' to eject gas before it is turned into stars. Without adjusting any parameters, the Bower et al. model gives an excellent match to the observed luminosity function of LRGs in the SDSS (with a median redshift of z=0.24) and to their clustering; the Baugh et al. model is less successful in these respects. Both models fail to match the observed abundance of LRGs at z=0.5 to better than a factor of ~2. In the models, LRGs are typically bulge dominated systems with M* of ~2x10^11 h^{-1} M_sun and velocity dispersions of ~250 km s^{-1}. Around half of the stellar mass in the model LRGs is already formed by z~2.2 and is assembled into one main progenitor by z~1.5; on average, only 25% of the mass of the main progenitor is added after z~1. LRGs are predicted to be found in a wide range of halo masses, a conclusion which relies on properly taking into account the scatter in the formation histories of haloes. Remarkably, we find that the correlation function of LRGs is predicted to be a power law down to small pair separations, in excellent agreement with observational estimates. Neither the Bower et al. nor the Baugh et al. model is able to reproduce the observed radii of LRGs.

(Abridged) X-ray clusters are conventionally divided into two classes: "cool core" (CC) clusters and "non-cool core" (NCC) clusters. Yet relatively little attention has been given to the origins of this dichotomy and, in particular, to the energetics and thermal histories of the two classes. We develop a model for the entropy profiles of clusters starting from the configuration established by gravitational shock heating and radiative cooling. At large radii, gravitational heating accounts for the observed profiles and their scalings well. However, at small and intermediate radii, radiative cooling and gravitational heating cannot be combined to explain the observed profiles of either type of cluster. The inferred entropy profiles of NCC clusters require that material is preheated prior to cluster collapse in order to explain the absence of low entropy (cool) material in these systems. We show that a similar modification is also required in CC clusters in order to match their properties at intermediate radii. In CC clusters, this modification is unstable, and an additional process is required to prevent cooling below a temperature of a few keV. We show that this can be achieved by adding a self-consistent AGN feedback loop in which the lowest-entropy, most rapidly cooling material is heated so that it rises buoyantly to mix with material at larger radii. The resulting model does not require fine tuning and is in excellent agreement with a wide variety of observational data. Some of the other implications of this model are briefly discussed.

The most massive galaxies in the Universe are also the oldest. To overturn this apparent contradiction with hierarchical growth models, we focus on the group scale haloes which host most of these galaxies. A stellar mass selected M_* >~ 2x10^10M_sol sample at z~0.4 is constructed within the CNOC2 redshift survey. A sensitive Mid InfraRed (MIR) IRAC colour is used to isolate passive galaxies. It produces a bimodal distribution, in which passive galaxies (highlighted by morphological early-types) define a tight MIR colour sequence (Infrared Passive Sequence, IPS). This is due to stellar atmospheric emission from old stellar populations. Significantly offset from the IPS are galaxies where reemission by dust boosts emission at 8microns (InfraRed-Excess or IRE galaxies). They include all known morphological late-types. Comparison with EW[OII] shows that MIR colour is highly sensitive to low levels of activity, and allows us to separate dusty-active from passive galaxies. The fraction of IRE galaxies, f(IRE) drops with M_*, such that f(IRE)=0.5 at a ``crossover mass'' of ~1.3x10^11M_sol. Within our optically-defined group sample there is a strong and consistent deficit in f(IRE) at all masses, and most clearly at M_* >~10^11M_sol. Using a mock galaxy catalogue derived from the Millenium Simulation we show that the observed trend of f(IRE) with M_* can be explained if suppression of star formation occurs primarily in the group environment, and particularly for M_*>~10^11M_sol galaxies. In this way, downsizing can be driven solely by structure growth in the Universe.

(Abridged) We demonstrate that the tenet of hierarchical structure growth leads directly to a robust, falsifiable prediction for the correlation between stellar fraction (fstar) and total system mass (M500) of galaxy groups and clusters. This prediction is relatively insensitive to the details of baryonic physics or cosmological parameters. In particular, if the fstar-M500 relation is fixed and does not evolve with redshift, CDM models predict the logarithmic slope of this relation to be b>-0.3. This constraint can be weakened if the fstar-M500 relation evolves strongly, but this implies more stars must be formed in situ in groups at low redshift. Conservatively requiring that at least half the stars in groups were formed by z=1, the constraint from evolution models is b>-0.35. Since the most massive clusters (M500=1E15 Msun) are observed to have fstar=0.01, this means that groups with M500=5E13 Msun must have fstar<0.03. Recent observations by Gonzalez et al. (2007) indicate a much steeper relation, with fstar>0.04 in groups, leading to b=-0.64. If confirmed, this would rule out hierarchical structure formation models: today's clusters could not have been built from today's groups, or even from the higher-redshift progenitors of those groups. We perform a careful analysis of these and other data to identify the most important systematic uncertainties in their measurements. Although correlated uncertainties on stellar and total masses might explain the steep observed relation, the data are only consistent with theory if the observed group masses are systematically underestimated.

In this paper, we present a new implementation of feedback due to active galactic nuclei (AGN) in cosmological simulations of galaxy formation. We assume that a fraction of jet energy, which is generated by an AGN, is transferred to the surrounding gas as thermal energy. Combining a theoretical model of mass accretion onto black holes with a multiphase description of star-forming gas, we self-consistently follow evolution of both galaxies and their central black holes. The novelty in our model is that we consider two distinct accretion modes: standard radiatively efficient thin accretion disks and radiatively inefficient accretion flows which we will generically refer to as RIAFs; motivated by theoretical models for jet production in accretion disks, we assume that only the RIAF is responsible for the AGN feedback. We find that, after an initial episode of bursting star formation, the accretion rate onto the central black hole drops so that the accretion disk switches to a RIAF structure. At this point, the feedback from the AGN becomes efficient and slightly suppresses star formation in the galactic disk and almost completely halts star formation in the bulge. As a result, the nucleus becomes a stochastically fuelled low-luminosity AGN (Seyfert galaxy) with recurrent short-lived episodes of activity after the star bursts. Our model predicts several properties of the low-luminosity AGN including the bolometric luminosity, jet powers, the effect on kpc-scale of the radio jet and the AGN lifetime, which are in broad agreement with observations of Seyfert galaxies and their radio activity. We also find that the mass ratios between the central black hole and the the host spheroid at z = 0 are ~10^{-3} regardless of the strength of either supernova feedback or AGN feedback. (abridged)

A common feature of hierarchical galaxy formation models is the process of "inverse" morphological transformation: a bulge dominated galaxy accretes a gas disk, dramatically reducing the system's bulge-to-disk mass ratio. During their formation, present day galaxies may execute many such cycles across the Hubble diagram. A good candidate for such a "hermaphrodite" galaxy is NGC 3108: a dust-lane early-type galaxy which has a large amount of HI gas distributed in a large scale disk. We present narrow band H_alpha and R-band imaging, and compare the results with the HI distribution. The emission is in two components: a nuclear bar and an extended disk component which coincides with the HI distribution. This suggests that a stellar disk is currently being formed out of the HI gas. The spatial distributions of the H_alpha and HI emission and the HII regions are consistent with a barred spiral structure, extending some 20 kpc in radius. We measure an extinction- corrected SFR of 0.42 Msun/yr. The luminosity function of the HII regions is similar to other spiral galaxies, with a power law index of -2.1, suggesting that the star formation mechanism is similar to other spiral galaxies. We measured the current disk mass and find that it is too massive to have been formed by the current SFR over the last few Gyr. It is likely that the SFR in NGC 3108 was higher in the past. With the current SFR, the disk in NGC 3108 will grow to be ~6.2x10^9 Msun in stellar mass within the next 5.5 Gyr. While this is substantial, the disk will be insignificant compared with the large bulge mass: the final stellar mass disk-to-bulge ratio will be ~0.02. NGC 3108 will fail to transform into anything resembling a spiral without a boost in the SFR and additional supply of gas.

We use a large suite of carefully controlled full hydrodynamic simulations to study the ram pressure stripping of the hot gaseous halos of galaxies as they fall into massive groups and clusters. The sensitivity of the results to the orbit, total galaxy mass, and galaxy structural properties is explored. For typical structural and orbital parameters, we find that ~30% of the initial hot galactic halo gas can remain in place after 10 Gyr. We propose a physically simple analytic model that describes the stripping seen in the simulations remarkably well. The model is analogous to the original formulation of Gunn & Gott (1972), except that it is appropriate for the case of a spherical (hot) gas distribution (as opposed to a face-on cold disk) and takes into account that stripping is not instantaneous but occurs on a characteristic timescale. The model reproduces the results of the simulations to within approximately 10% at almost all times for all the orbits, mass ratios, and galaxy structural properties we have explored. The one exception involves unlikely systems where the orbit of the galaxy is highly non-radial and its mass exceeds about 10% of the group or cluster into which it is falling (in which case the model under-predicts the stripping following pericentric passage). The proposed model has several interesting applications, including modelling the ram pressure stripping of both observed and cosmologically-simulated galaxies and as a way to improve current semi-analytic models of galaxy formation. One immediate consequence is that the colours and morphologies of satellite galaxies in groups and clusters will differ significantly from those predicted with the standard assumption of complete stripping of the hot coronae.

We present a detailed study of the structural properties of four gravitationally lensed disk galaxies at z=1. Modelling the rotation curves on sub-kpc scales we derive the values for the disk mass, the reference dark matter density and core radius, and the angular momentum per unit mass. The derived models suggest that the rotation curve profile and amplitude are best fit with a dark matter component similar to those of local spiral galaxies. The stellar component also has a similar length scale, but with substantially smaller masses than similarly luminous disk galaxies in the local universe. Comparing the average dark matter density inside the optical radius we find that the disk galaxies at z=1 have larger densities (by up to a factor of 7) than similar disk galaxies in the local Universe. Furthermore, the angular momentum per unit mass versus reference velocity is well matched to the local relation, suggesting that the angular momentum of the disk remains constant between high redshifts and the present day. Though statistically limited, these observations point towards a spirals' formation scenario in which stellar disks are slowly grown by the accretion of angular momentum conserving material.

The power of jets from black holes are expected to depend on both the spin of the black hole and the structure of the accretion disk in the region of the last stable orbit. We investigate these dependencies using two different physical models for the jet power: the classical Blandford-Znajek (BZ) model and a hybrid model developed by Meier. In the BZ case, the jets are powered by magnetic fields directly threading the spinning black hole while in the hybrid model, the jet energy is extracted from both the accretion disk as well as the black hole via magnetic fields anchored to the accretion flow inside and outside the hole's ergosphere. The hybrid model takes advantage of the strengths of both the Blandford-Payne and BZ mechanisms, while avoiding the more controversial features of the latter. We develop these models more fully to account for general relativistic effects and to focus on advection-dominated accretion flows (ADAF) for which the jet power is expected to be a significant fraction of the accreted rest mass energy. We apply the models to elliptical galaxies, in order to see if these models can explain the observed correlation between the Bondi accretion rates and the total jet powers. For typical values of the disk viscosity parameter alpha~0.04-0.3 and mass accretion rates consistent with ADAF model expectations, we find that the observed correlation requires j>0.9; i.e., it implies that the black holes are rapidly spinning. Our results suggest that the central black holes in the cores of clusters of galaxies must be rapidly rotating in order to drive jets powerful enough to heat the intracluster medium and quench cooling flows.

The largest uncertainty for cosmological studies using clusters of galaxies is introduced by our limited knowledge of the statistics of galaxy cluster structure, and of the scaling relations between observables and cluster mass. To improve on this situation we have started an XMM-Newton Large Programme for the in-depth study of a representative sample of 33 galaxy clusters, selected in the redshift range z=0.055 to 0.183 from the REFLEX Cluster Survey, having X-ray luminosities above 0.4 X 10^44 h_70^-2 erg s^-1 in the 0.1 - 2.4 keV band. This paper introduces the sample, compiles properties of the clusters, and provides detailed information on the sample selection function. We describe the selection of a nearby galaxy cluster sample that makes optimal use of the XMM-Newton field-of-view, and provides nearly homogeneous X-ray luminosity coverage for the full range from poor clusters to the most massive objects in the Universe. For the clusters in the sample, X-ray fluxes are derived and compared to the previously obtained fluxes from the ROSAT All-Sky Survey. We find that the fluxes and the flux errors have been reliably determined in the ROSAT All-Sky Survey analysis used for the REFLEX Survey. We use the sample selection function documented in detail in this paper to determine the X-ray luminosity function, and compare it with the luminosity function of the entire REFLEX sample. We also discuss morphological peculiarities of some of the sample members. The sample and some of the background data given in this introductory paper will be important for the application of these data in the detailed studies of cluster structure, to appear in forthcoming publications.

The universal baryonic mass fraction (Omega_b/Omega_m) can be sensitively constrained using X-ray observations of galaxy clusters. In this paper, we compare the baryonic mass fraction inferred from measurements of the cosmic microwave background with the gas mass fractions (f_gas) of a large sample of clusters taken from the recent literature. In systems cooler than 4 keV, f_gas declines as the system temperature decreases. However, in higher temperature systems, f_gas(r500) converges to approx. (0.12 +/- 0.02)(h/0.72)^{-1.5}, where the uncertainty reflects the systematic variations between clusters at r500. This is significantly lower than the maximum-likelihood value of the baryon fraction from the recently released WMAP 3-year results. We investigate possible reasons for this discrepancy, including the effects of radiative cooling and non-gravitational heating, and conclude that the most likely solution is that Omega_m is higher than the best-fit WMAP value (we find Omega_m = 0.36^{+0.11}_{-0.08}), but consistent at the 2-sigma level. Degeneracies within the WMAP data require that sigma_8 must also be greater than the maximum likelihood value for consistency between the data sets.

(Abridged) We report a discovery of possible large-scale structures around the RDCS J1252.9-2927 cluster at z=1.24 based on photometric redshifts. We carried out multi-band wide-field imaging with Suprime-Cam on the Subaru Telescope and WFCAM on the United Kingdom Infra-Red Telescope (UKIRT). The distribution of photo-z selected galaxies reveals clumpy structures surrounding the central cluster. We compare the observed structure with an X-ray map and find that two of the four plausible clumps show significant X-ray emissions and one with a marginal detection, which strongly suggest that they are dynamically bound systems. Following the discovery of the possible large-scale structure, we carried out deeper SOFI K_s-band imaging with New Technology Telescope on the four plausible clumps. We construct the optical-to-near-infrared colour-magnitude diagrams of the galaxies in the clumps, and find that the colour-magnitude relation (CMR) of the red galaxies in the clumps is sharply truncated below K_s=22. Interestingly, the main cluster shows a clear relation down to K_s=23 (Lidman et al. 2004). We suggest that galaxies follow the 'environment-dependent down-sizing' evolution. Massive galaxies in high density environments first stop forming stars and become red. Less massive galaxies in less dense environments become red at later times. Based on a few assumptions, we predict that the brightest tip of the CMR appears at z~2.5.

In McCarthy et al. (2007, astro-ph/0701335) we present an analysis of entropy generation during cluster mergers, using SPH simulations. This conference proceeding contribution summarizes some of the most important findings of that paper, without all the technical details. Our main result is that during cluster mergers, entropy is generated in two major shock episodes, that heat the gas from the inside-out, in a way that distributes most of the energy within the more massive clump. This is different from the "onion-skin" model adopted in many simple analytic models of cluster formation.

(Abridged) The thermal history of the intracluster medium (ICM) is complex. Heat input from cluster mergers, AGN, and galaxy winds offsets and may even halt the cooling of the ICM. Consequently, the processes that set the properties of the ICM play a key role in determining how galaxies form. In this paper we focus on the shock heating of the ICM during cluster mergers, with the eventual aim of incorporating this mechanism into semi-analytic models of galaxy formation. We use a suite of hydrodynamic simulations to track the evolution of the ICM in idealised two-body mergers. We find the heating of the ICM can be understood relatively simply by considering the evolution of the gas entropy during the mergers. We examine the processes that generate the entropy in order to understand why previous analytic shock heating models failed. We find that: (1) The energy that is thermalised in the collision greatly exceeds the kinetic energy available when the systems first touch. The smaller system penetrates deep into the potential well before it is disrupted. (2) For unequal mass mergers, most of the energy is thermalised in the more massive component. The heating of the smaller system is minor and its gas sinks to the centre of the final system. (3) The bulk of the entropy generation occurs in two distinct episodes. The first episode occurs following core collision, when a shock wave is generated that propagates outwards from the centre. This causes the combined system to expand rapidly and overshoot hydrostatic equilibrium. The second episode occurs as this material is shock heated as it re-collapses. This revised model for entropy generation significantly improves our physical understanding of cosmological gas simulations.

By combining HST imaging with optical (VIMOS) and near-infrared (SINFONI) integral field spectroscopy we exploit the gravitational potential of a massive, rich cluster at z=0.9 to study the internal properties of a gravitationally lensed galaxy at z=4.88. Using a detailed gravitational lens model of the cluster RCS0224-002 we reconstruct the source-frame morphology of the lensed galaxy on 200pc scales and find an ~L* Lyman-break galaxy with an intrinsic size of only 2.0x0.8kpc, a velocity gradient of <60km/s and an implied dynamical mass of 1.0x10^10Mo within 2kpc. We infer an integrated star-formation rate of just 12+/-2Mo/yr from the intrinsic [OII] emission line flux. The Ly-alpha emission appears redshifted by +200+/-40km/s with respect to the [OII] emission. The Ly-alpha is also significantly more extended than the nebular emission, extending over 11.9x2.4kpc. Over this area, the Ly-alpha centroid varies by less than 10km/s. By examining the spatially resolved structure of the [OII] and asymmetric Ly-alpha emission lines we investigate the nature of this system. The model for local starburst galaxies suggested by Mass-Hesse et al. (2003) provides a good description of our data, and suggests that the galaxy is surrounded by a galactic-scale bi-polar outflow which has recently burst out of the system. The outflow, which appears to be currently located >30kpc from the galaxy, is escaping at a speed of upto ~500km/s. Although the mass of the outflow is uncertain, the geometry and velocity of the outflow suggests that the ejected material is travelling far faster than escape velocity and will travel more than 1Mpc (comoving) before eventually stalling.

Recently, Allen et al. measured a tight correlation between the Bondi accretion rates and jet powers of the nuclei of nearby X-ray luminous elliptical galaxies. We employ two models of jet powering to understand the above correlation and derive constraints on the spin and accretion rate of the central black holes. The first is the Blandford-Znajek model, in which the spin energy of the hole is extracted as jet power; the second model is an hybrid version of the Blandford-Payne and Blandford-Znajek processes, in which the outflow is generated in the inner parts of the accretion disk. We assume advection-dominated accretion flows (ADAF) and account for general relativistic effects. Our modelling implies that for typical values of the disk viscosity parameter alpha~0.01-1 the tight correlation implies the narrow range of spins j~0.7-1 and accretion rates Mdot(R_ms)~(0.01-1) Mdot_Bondi. Our results provide support for the "spin paradigm" scenario and suggest that the central black holes in the cores of clusters of galaxies must be rapidly rotating in order to drive radio jets powerful enough to quench the cooling flows.

We use an array of high-resolution N-body simulations to determine the mass function of dark matter haloes at redshifts 10-30. We develop a new method for compensating for the effects of finite simulation volume that allows us to find an approximation to the true ``global'' mass function. By simulating a wide range of volumes at different mass resolution, we calculate the abundance of haloes of mass 10^{5-12} Msun/h. This enables us to predict accurately the abundance of the haloes that host the sources that reionize the universe. In particular, we focus on the small mass haloes (>~10^{5.5-6} Msun/h) likely to harbour population III stars where gas cools by molecular hydrogen emission, early galaxies in which baryons cool by atomic hydrogen emission at a virial temperature of ~10^4K (10^{7.5-8} Msun/h), and massive galaxies that may be observable at redshift ~10. When we combine our data with simulations that include high mass halos at low redshift, we find that the best fit to the halo mass function depends not only on linear overdensity, as is commonly assumed in analytic models, but also upon the slope of the linear power spectrum at the scale of the halo mass. The Press-Schechter model gives a poor fit to the halo mass function in the simulations at all epochs; the Sheth-Tormen model gives a better match, but still overpredicts the abundance of rare objects at all times by up to 50%. Finally, we consider the consequences of the recently released WMAP 3-year cosmological parameters. These lead to much less structure at high redshift, reducing the number of z=10 ``mini-haloes'' by more than a factor of two and the number of z=30 galaxy hosts by more than four orders of magnitude. Code to generate our best-fit halo mass function may be downloaded from http://icc.dur.ac.uk/Research/PublicDownloads/genmf_readme.html

Models of disk galaxy formation commonly predict the existence of an extended reservoir of hot gas surrounding massive spirals at low redshift. As a test of these models, we have obtained X-ray and optical data of the two massive edge-on spirals NGC 5746 and NGC 5170, in order to investigate the amount and origin of hot gas in their disks and halos. Chandra observations of NGC 5746 reveal evidence for diffuse X-ray emission with a total luminosity of ~7 x 10^39 erg/s surrounding this galaxy out to at least ~20 kpc from the disk, whereas an identical study of the less massive NGC 5170 fails to detect any extraplanar X-ray emission. Unlike the case for other disk galaxies with detected X-ray halos, the halo emission around NGC 5746 is not accompanied by extraplanar H-alpha or radio emission, and there is no evidence for significant nuclear or starburst activity in the disk. In contrast to these other cases, the emission around NGC 5746 therefore appears to arise from the cooling of externally accreted material rather than from disk outflows. To verify this idea, we present results of cosmological simulations of galaxy formation and evolution, showing our observations to be in good agreement with expectations for cosmological accretion, while also confirming that the X-ray halos of other spirals do not fit well into an accretion scenario. We find that the estimated cooling rate of hot halo gas around NGC 5746 would provide sufficient material for star formation in the disk to proceed at its present rate. This lends support to the idea that a supply of hot ambient gas is potentially available as fuel for star formation in massive, nearby spirals, and suggests that accretion of hot gas could be important for maintaining the stellar disks of such galaxies.

We have obtained near-infrared imaging of 58 galaxy groups, in the redshift range 0.1<z<0.6, from the William Herschel Telescope and from the Spitzer IRAC data archive. The groups are selected from the CNOC2 redshift survey, with additional spectroscopy from the Baade telescope (Magellan). Our group samples are statistically complete to K(Vega}=17.7 (INGRID) and [3.6](AB)=19.9 (IRAC). From these data we construct near-infrared luminosity functions, for groups in bins of velocity dispersion, up to 800 km/s, and redshift. The total amount of near-infrared luminosity per group is compared with the dynamical mass, estimated from the velocity dispersion, to compute the mass-to-light ratio, M/Lk. We find that the M/Lk values in these groups are in good agreement with those of their statistical descendants at z=0, with no evidence for evolution beyond that expected for a passively evolving population. There is a trend of M/Lk with group mass, which increases from M/Lk~10 for groups with sigma<250 km/s to M/Lk~100 for sigma=425-800 km/s. This trend is weaker, but still present, if we estimate the total mass from weak lensing measurements. In terms of stellar mass, stars make up ~2 per cent of the mass in the smallest groups, and ~1 per cent in the most massive groups. We also use the near-infrared data to consider the correlations between stellar populations and stellar masses, for group and field galaxies at 0.1<z<0.6. We find that fewer group galaxies show strong [OII] emission compared with field galaxies of the same stellar mass and at the same redshift. We conclude that most of the stellar mass in these groups was already in place by z~0.4, with little environment-driven evolution to the present day.

We analyse a z<0.1 galaxy sample from the Sloan Digital Sky Survey focusing on the variation of the galaxy colour bimodality with stellar mass and projected neighbour density Sigma, and on measurements of the galaxy stellar mass functions. The characteristic mass increases with environmental density from about 10^10.6 Msun to 10^10.9 Msun (Kroupa IMF, H_0=70) for Sigma in the range 0.1--10 per Mpc^2. The galaxy population naturally divides into a red and blue sequence with the locus of the sequences in colour-mass and colour-concentration index not varying strongly with environment. The fraction of galaxies on the red sequence is determined in bins of 0.2 in log Sigma and log mass (12 x 13 bins). The red fraction f_r generally increases continuously in both Sigma and mass such that there is a unified relation: f_r = F(Sigma,mass). Two simple functions are proposed which provide good fits to the data. These data are compared with analogous quantities in semi-analytical models based on the Millennium N-body simulation: the Bower et al. (2006) and Croton et al. (2006) models that incorporate AGN feedback. Both models predict a strong dependence of the red fraction on stellar mass and environment that is qualitatively similar to the observations. However, a quantitative comparison shows that the Bower et al. model is a significantly better match; this appears to be due to the different treatment of feedback in central galaxies.

We examine the prospects for the next generation of surveys aimed at elucidating the nature of dark energy. We review the methods that can be used to determine the redshift evolution of the dark energy equation of state parameter w, highlighting their respective strengths and potential weaknesses. All of the attractive methods require surveys covering more than 5-10,000 square degrees of the sky. We examine the accuracy that each method is likely to deliver within a decade, and discuss the difficulties arising from systematic uncertainties associated with the techniques. We conclude that the proposed photometric and redshift surveys have the potential of delivering measurements of w with percent accuracy at several redshifts out to z ~ 3 . Of particular interest will be the combination of weak lensing and baryonic acoustic oscillations measurements. This exquisite precision is likely to have a fundamental impact on our understanding of the nature of dark energy, providing the necessary guidance for its theoretical explanation.

We report a discovery of proto-cluster candidates around high redshift radio galaxies at z~2.5 on the basis of clear statistical excess of colour-selected galaxies around them seen in the deep near-infrared imaging data obtained with CISCO on Subaru Telescope. We have observed six targets, all at similar redshifts at z~2.5, and our data reach to J=23.5, H=22.6 and K=21.8 (5 sigma) and cover a 1.6' x 1.6' field centered on each radio galaxy. We apply colour cuts in JHKs in order to exclusively search for galaxies located at high redshifts, z>2. Over the magnitude range of 19.5<K<21.5 we see a significant excess of red galaxies with J-K>2.3 by a factor of two around the combined radio galaxies fields compared to those found in the general field of GOODS South. The excess of galaxies around the radio galaxies fields becomes more than factor of three around 19.5<K<20.5 when the two-colours cuts are applied with JHKs. Such overdensity of the colour-selected galaxies suggest that those fields tend to host high density regions at high redshifts, although there seems to be the variety of the density of the colour-selected galaxies in each field. In particular, two radio galaxies fields out of the six observed fields show very strong density excess and these are likely to be proto-clusters associated to the radio galaxies which would evolve into rich clusters of galaxies dominated by old passively evolving galaxies.

Recent observations of the distant Universe suggest that much of the stellar mass of bright galaxies was already in place at $z>1$. This presents a challenge for models of galaxy formation because massive halos are assembled late in hierarchical cosmologies such as cold dark matter (CDM). In this paper, we discuss a new implementation of the Durham semi-analytic model in which feedback due to active galactic nuclei (AGN) is assumed to quench cooling flows in massive halos. This mechanism naturally creates a break in the local galaxy luminosity function at bright magnitudes. The model is implemented within the Millennium N-body simulation; the accurate dark matter merger trees and large number of realizations of the galaxy formation process that the simulation provides results in highly accurate statistics. After adjusting the values of the physical parameters in the model by reference to the properties of local galaxies, we use it to investigate the evolution of the K-band luminosity and galaxy stellar mass functions. We also calculate the volume averaged star formation rate density of the Universe as a function of redshift and the way in which this is apportioned amongst galaxies of different mass. The model robustly predicts a substantial population of massive galaxies out to redshift $z\sim 5$ and a star formation rate density which rises with increasing redshift in objects of all masses. Although observational data on these properties have been cited as evidence for ``anti-hierarchical'' galaxy formation, we find that when AGN feedback is taken into account, the fundamentally hierachical CDM model provides a very good match to these observations.

We exploit the gravitational potential of massive cluster lenses to probe the emission line properties of six z=1 galaxies which appear as highly magnified luminous arcs. Using the GMOS integral field spectrograph together with detailed cluster lens models we reconstruct the intrinsic morphologies and two-dimensional velocity fields in these galaxies on scales corresponds to \~0.5kpc (unlensed) at z=1. Four of the galaxies have stable disk-like kinematics, whilst the other two resemble interacting or starburst galaxies. These galaxies lie close to the mean rest-frame I-band Tully-Fisher relation for nearby spirals suggesting a clear preference for hierarchical growth of structure. In the rest-frame B-band, the observations suggest 0.5+/-0.3mag of brightening, consistent with increased star-formation activity at z=1. However, the galaxies with stable disk kinematics have more slowly rising rotation curves than expected from galaxies with similar surface brightness in the local Universe. We suggest that this may arise because the distant galaxies have lower bulge masses than their local counter-parts. Whilst this study is based on only six galaxies, the gain in flux and in spatial resolution achieved via gravitational magnification provides a much more detailed view of the high redshift Universe than possible with conventional surveys.

Hot gaseous haloes surrounding galaxies and extending well beyond the distribution of stars are a ubiquitous prediction of galaxy formation scenarios. The haloes are believed to consist of gravitationally trapped gas with a temperature of millions of Kelvin. The existence of such hot haloes around massive elliptical galaxies has been established through their X-ray emission. While gas out-flowing from starburst spiral galaxies has been detected, searches for hot haloes around normal, quiescent spiral galaxies have so far failed, casting doubts on the fundamental physics in galaxy formation models. Here we present the first detection of a hot, large-scale gaseous halo surrounding a normal, quiescent spiral galaxy, NGC 5746, alleviating a long-standing problem for galaxy formation models. In contrast to starburst galaxies, where the X-ray halo can be powered by the supernova energy, there is no such power source in NGC 5746. The only compelling explanation is that we are here witnessing a galaxy forming from gradually in-flowing hot and dilute halo gas.

We study how the proportion of star-forming galaxies evolves between z=0.8 and z=0 as a function of galaxy environment, using the [OII] line in emission as a signature of ongoing star formation. Our high-z dataset comprises 16 clusters, 10 groups and another 250 galaxies in poorer groups and the field at z=0.4-0.8 from the ESO Distant Cluster Survey, plus another 9 massive clusters at similar redshifts. As a local comparison, we use samples of galaxy systems selected from the Sloan Digital Sky Survey at 0.04< z < 0.08. At high-z most systems follow a broad anticorrelation between the fraction of star-forming galaxies and the system velocity dispersion. At face value, this suggests that at z=0.4-0.8 the mass of the system largely determines the proportion of galaxies with ongoing star formation. At these redshifts the strength of star formation (as measured by the [OII] equivalent width) in star-forming galaxies is also found to vary systematically with environment. Sloan clusters have much lower fractions of star-forming galaxies than clusters at z=0.4-0.8 and, in contrast with the distant clusters, show a plateau for velocity dispersions $ \ge 550 km s^-1$, where the fraction of galaxies with [OII] emission does not vary systematically with velocity dispersion. We quantify the evolution of the proportion of star-forming galaxies as a function of the system velocity dispersion and find it is strongest in intermediate-mass systems (sigma ~ 500-600 km s^-1 at z=0). To understand the origin of the observed trends, we use the Press-Schechter formalism and the Millennium Simulation and show that galaxy star formation histories may be closely related to the growth history of clusters and groups. We propose a scheme that is able to account for the observed relations between the star-forming fraction and \sigma [abridged].

We model gas cooling in high-resolution N-body simulations in order to investigate the formation of the first generation of stars. We follow a region of a LCDM universe especially selected to contain a rich cluster by the present day. The properties of the dark haloes that form in these sub-solar mass resolution simulations are presented in a companion paper by Gao et al. The first gas clouds able to cool by molecular hydrogen line emission collapse at extremely high redshift, z ~ 47, when the mass of the dark halo is 2.4 x 10^5 Msun/h. By z ~ 30, a substantial population of haloes are capable of undergoing molecular hydrogen cooling although their ability to form stars is dependent on the efficiency of feedback processes such as dissociating Lyman-Werner radiation. The mass of the main halo grows extremely rapidly and, by z ~ 36, its virial temperature has reached 10^4K, at which point gas cooling becomes dominated by more effective atomic line processes. By z ~ 30, a small ``group'' of such potential galaxies will have formed unless prevented from doing so by feedback processes. By this redshift, massive (>~ 100 Msun) population III stars are able to ionise gas well beyond their own host halo and neighbouring HII regions can percolate to form an ionized superbubble. Such patches would be too widely separated to contribute significantly to reionisation at this time. The large number density of early cooling haloes in the pre-reionised universe raises the exciting prospect that this ultra-early generation of stars may be observable as gamma-ray bursts or supernovae.

The evolution of galaxies in groups may have important implications for the global evolution of star formation rate in the Universe, since many processes which operate in groups may suppress star formation, and the fraction of galaxies bound in groups at the present day is as high as ~60%. We present an analysis of our sample of 0.3<=z<=0.55 groups, selected from the CNOC2 redshift survey and supplemented with deep spectroscopy and HST ACS imaging. We find that these groups contain significantly more passive galaxies than the field, with excesses of S0, elliptical and passive spiral galaxy types. The morphological composition is closely matched to that of more massive irregular clusters at a similar epoch. Contrasting with galaxy samples in a variety of environments and epochs, we find that the fraction of passive galaxies (EW[OII]<5A), is strongly evolving in the group environment, with parallel evolution in the (global) field population, whilst little evolution is observed in cluster cores since z~1.

High-velocity galactic outflows, driven by intense bursts of star formation and black hole accretion, are invoked by current theories of galaxy formation to terminate star formation in the most massive galaxies and to deposit heavy elements in the intergalactic medium. From existing observational evidence on high-redshift galaxies, it is unclear whether such outflows are localized to regions of intense star formation just a few kiloparsecs in extent, or whether they instead have a significant impact on the entire galaxy and its surroundings. Here we present two-dimensional spectroscopy of a star-forming galaxy at redshift z=3.09 (seen 11.5 Gyr ago, when the Universe was 20 per cent of its current age): its spatially extended Ly-alpha emission appears to be absorbed by HI in a foreground screen covering the entire galaxy, with a lateral extent of at least 100 kpc and remarkable velocity coherence. It was plausibly ejected from the galaxy during a starburst several 1E8 yr earlier and has subsequently swept up gas from the surrounding intergalactic medium and cooled. This demonstrates the galaxy-wide impact of high-redshift superwinds.

Acoustic oscillations in the baryon-photon fluid leave a signature in the matter power spectrum. The overall shape of the spectrum and the wavelength of the oscillations depend upon the sound horizon scale at recombination. Using the $\Lambda$ cold dark matter Hubble Volume simulation, we show that the imprint of baryons is visible in the power spectrum of cluster-mass dark matter haloes, in spite of significant differences between the halo power spectrum and the prediction of linear perturbation theory. A measurement of the sound horizon scale can constrain the dark energy equation of state. We show that a survey of clusters at intermediate redshift ($ z\sim1 $), like the Sunyaev-Zeldovich survey proposed by the South Pole Telescope or a red sequence photometric survey with VISTA, could potentially constrain the sound horizon scale to an accuracy of $\sim 2%$, in turn fixing the ratio of the pressure of the dark energy to its density ($w$) to better than $\sim 10%$. Our approach does not require knowledge of the cluster mass, unlike those that depend upon the abundance of clusters.

We present optical and near-infrared integral field spectroscopy of the SCUBA galaxy SMM J163650.43+405734.5 (ELAIS N2 850.4) at z=2.385. We combine Ly-alpha and H-alpha emission line maps and velocity structure with high resolution HST ACS and NICMOS imaging to probe the complex dynamics of this vigorous star-burst galaxy. The imaging data shows a complex morphology, consisting of at least three components separated by ~1'' (8kpc) in projection. When combined with the H-alpha velocity field from UKIRT UIST IFU observations we identify two components whose redshifts are coincident with the systemic redshift, measured from previous CO observations, one of which shows signs of AGN activity. A third component is offset by 220+/-50km/s from the systemic velocity. The total star formation rate of the whole system (estimated from the narrow-line H-alpha and uncorrected for reddening) is 340+/-50Mo/yr. The Ly-alpha emission mapped by the GMOS IFU covers the complete galaxy and is offset by +270+/-40km/s from the systemic velocity. This velocity offset is comparable to that seen in rest-frame UV-selected galaxies at similar redshifts and usually interpreted as a star-burst driven wind. The extended structure of the Ly-alpha emission suggests that this wind is not a nuclear phenomenon, but is instead a galactic scale outflow. Our observations suggest that the vigorous activity in N2 850.4 is arising as a result of an interaction between at least two dynamically-distinct components, resulting in a strong starburst, a starburst-driven wind and actively-fuelled AGN activity. [abridged]

We investigate a spectroscopic sample of 48 early-type (E+S0) galaxies in the rich cluster Abell 2390 (z=0.23) and 48 E+S0 galaxies from a previously published survey of Abell 2218 (z=0.18). The A2390 spectra are based on Multi Object Spectroscopy with MOSCA at the Calar Alto 3.5-m, complemented by ground-based and HST F555W and F814W observations. Our investigation spans a broad range in luminosity (-20.5>M_r>-23.0) and a wide FOV of 1.5x1.5 Mpc^2. Since the A2218 and A2390 samples are very similar, we can combine them and analyse a total of 96 E+S0 galaxies at z~0.2. We construct the Faber-Jackson relation (FJR) for all E+S0s and detect a modest luminosity evolution. The average offset from the local FJR in Gunn r-band is (DeltaM_r)=0.32+-0.22^m. Less-massive galaxies show a trend for a larger evolution than more-massive galaxies. HST/WFPC2 surface brightness profile fits were used to derive structural parameters. For 34 E+S0 galaxies which enter the Fundamental Plane (FP) we deduce a mild evolution in Gunn r with a zero-point offset of 0.10+-0.06, corresponding to 0.31+-0.18^m. Elliptical and lenticular galaxies are uniformly distributed along the FP with a similar scatter of 0.1 dex. We find little evidence for differences between the populations of ellipticals and S0s. Lenticulars induce on average a larger evolution of 0.44+-0.18^m than ellipticals with 0.02+-0.21^m. The M/L ratios at z=0.2 are offset by (Deltalog(M/L_r))=-0.12+-0.06 dex compared to those of Coma. Our results can be reconciled with a passive evolution of the stellar populations and a high formation redshift for the bulk of the stars in E+S0 galaxies. However, our findings are also consistent with the hierarchical formation picture for rich clusters. [abridged]

We compare deep Magellan spectroscopy of 26 groups at 0.3 <= z <= 0.55, selected from the Canadian Network for Observational Cosmology 2 field survey (CNOC2), with a large sample of nearby groups from the 2PIGG catalogue (Eke et al., 2004). We find that the fraction of group galaxies with significant [OII] emission (>=5\AA) increases strongly with redshift, from ~29% in 2dFGRS to ~58% in CNOC2, for all galaxies brighter than ~ M*+1.75. This trend is parallel to the evolution of field galaxies, where the equivalent fraction of emission line galaxies increases from ~ 53% to ~ 75%. The fraction of emission-line galaxies in groups is lower than in the field, across the full redshift range, indicating that the history of star formation in groups is influenced by their environment. We show that the evolution required to explain the data is inconsistent with a quiescent model of galaxy evolution; instead, discrete events in which galaxies cease forming stars (truncation events) are required. We constrain the probability of truncation (P_trunc) and find that a high value is required in a simple evolutionary scenario neglecting galaxy mergers (P_trunc >~ 0.3 Gyr^{-1}). However, without assuming significant density evolution, P_trunc is not required to be larger in groups than in the field, suggesting that the environmental dependence of star formation was embedded at redshifts z >~ 0.45.

The evolution of galaxies in groups may have important implications for the evolution of the star formation history of the universe, since many processes which operate in groups may suppress star formation and the fraction of galaxies in bound groups grows rapidly between z=1 and the present day. In this paper, we present an investigation of the properties of galaxies in galaxy groups at intermediate redshift (z ~ 0.4). The groups were selected from the CNOC2 redshift survey as described in Carlberg et al., 2001, with further spectroscopic follow-up undertaken at the Magellan telescope in order to improve the completeness and depth of the sample. We present the data for the individual groups, and find no clear trend in the fraction of passive galaxies with group velocity dispersion and group concentration. We stack the galaxy groups in order to compare the properties of group galaxies with those of field galaxies at the same redshift. The groups contain a larger fraction of passive galaxies than the field, this trend being particularly clear for galaxies brighter than M_{B_J} < -20 in the higher velocity dispersion groups. In addition, we see evidence for an excess of bright passive galaxies in the groups relative to the field. In contrast, the luminosity functions of the star forming galaxies in the groups and the field are consistent. These trends are qualitatively consitent with the differences between group and field galaxies seen in the local universe.

In this paper we present the three dimensional dynamics of the galaxy SDSS J101345.39+011613.66, selected for its unusually strong Balmer absorption lines (Wo(H-delta)=7.5A). Using the GMOS-South IFU in Nod & Shuffle mode we have mapped the continuum and optical absorption lines of this z=0.1055 field galaxy. This galaxy has a disturbed morphology, with a halo of diffuse material distributed asymmetrically toward the north. Using the [OII] emission line (Wo([OII])=4.1A) we find that the gas and hot OB stars are offset from the older stars in the system. The gas also has a spatially extended and elongated morphology with a velocity gradient of 100+/-20km/s across 6kpc in projection. Using the strong H-gamma and H-delta absorption lines we find that the A- stars are widely distributed across the system and are not centrally concentrated arguing that the A-star population has formed in molecular clouds outside the nucleus. By cross correlating the spectra from the datacube with an A-star template we find evidence that the A-star population has a 40km/s shear in the same direction as the gas. The disturbed morphology, strong colour gradients and strong H-delta and H-gamma absorption lines in SDSS J101345.39 argue that this is a recent tidal interaction/merger between a passive elliptical and star-forming galaxy. Although based on a single object, these results show that we can spatially resolve and constrain the dynamics of this short lived (yet important) phase of galaxy formation in which the evolutionary process take galaxies from star-forming to their quiescent end products.

We investigate the evolution of the star formation rate in cluster galaxies. We complement data from the CNOC1 cluster survey (0.15<z<0.6) with measurements from galaxy clusters in the 2dF galaxy redshift survey (0.05<z<0.1) and measurements from recently published work on higher redshift clusters, up to almost z=1. We focus our attention on galaxies in the cluster core, ie. galaxies with r<0.7h^{-1}_{70}Mpc. Averaging over clusters in redshift bins, we find that the fraction of galaxies with strong [OII] emission is < 20% in cluster cores, and the fraction evolves little with redshift. In contrast, field galaxies from the survey show a very strong increase over the same redshift range. It thus appears that the environment in the cores of rich clusters is hostile to star formation at all the redshifts studied. We compare this result with the evolution of the colours of galaxies in cluster cores, first reported by Butcher & Oemler (1984). Using the same galaxies for our analysis of the [OII] emission, we confirm that the fraction of blue galaxies, which are defined as galaxies 0.2 mag bluer in the rest frame B-V than the red sequence of each cluster, increases strongly with redshift. Since the colours of galaxies retain a memory of their recent star formation history, while emission from the [OII] line does not, we suggest that these two results can best be reconciled if the rate at which the clusters are being assembled is higher in the past, and the galaxies from which it is being assembled are typically bluer.

We use a sample of 69726 galaxies from the SDSS to study the variation of the bimodal color-magnitude (CM) distribution with environment. Dividing the galaxy population by environment (Sigma_5) and luminosity (-23<M_r<-17), the u-r color functions are modeled using double-Gaussian functions. This enables a deconvolution of the CM distributions into two populations: red and blue sequences. The changes with increasing environmental density can be separated into two effects: a large increase in the fraction of galaxies in the red distribution, and a small color shift in the CM relations of each distribution. The average color shifts are 0.05+-0.01 and 0.11+-0.02 for the red and blue distributions, respectively, over a factor of 100 in projected neighbor density. The red fraction varies between about 0% and 70% for low-luminosity galaxies and between about 50% and 90% for high-luminosity galaxies. This difference is also shown by the variation of the luminosity functions with environment. We demonstrate that the effects of environment and luminosity can be unified. A combined quantity, Sigma_mod = Sigma_5/Mpc^{-2} + L_r/L_{-20.2}, predicts the fraction of red galaxies, which may be related to the probability of transformation events. Our results are consistent with major interactions (mergers and/or harassment) causing galaxies to transform from the blue to the red distribution. We discuss this and other implications for galaxy evolution from earlier results and model the effect of slow transformations on the color functions.

The observed cooling rate of hot gas in clusters is much lower than that inferred from the gas density profiles. This suggests that the gas is being heated by some source. We use an adaptive-mesh refinement code (FLASH) to simulate the effect of multiple, randomly positioned, injections of thermal energy within 50 kpc of the centre of an initially isothermal cluster with mass M_200=3x10^(14) Msol and kT=3.1 keV. We have performed eight simulations with spherical bubbles of energy generated every 10^8 years, over a total of 1.5 Gyr. Each bubble is created by injecting thermal energy steadily for 10^7 years; the total energy of each bubble ranges from 0.1--3x10^(60) erg, depending on the simulation. We find that 2x10^(60) erg per bubble (corresponding to a average power of 6.3x10^(44) erg/s) effectively balances energy loss in the cluster and prevents the accumulation of gas below kT=1 keV from exceeding the observational limits of 30 Msol/yr. This injection rate is comparable to the radiated luminosity of the cluster, and the required energy and periodic timescale of events are consistent with observations of bubbles produced by central AGN in clusters. The effectiveness of this process depends primarily on the total amount of injected energy and the initial location of the bubbles, but is relatively insensitive to the exact duty cycle of events.

In this letter we present the first 3D spectroscopic study of H-alpha emission in a z ~ 1 field galaxy with an integral field unit. Using the CIRPASS spectrograph on Gemini-South we map the spatial and velocity distribution of H-alpha emission in the z=0.819 galaxy CFRS 22.1313. We detect two H-alpha emitting regions with a velocity separation of 220+/-10 km/s. Combining the 2D map of H-alpha emission with HST F814W imaging, we determine a lower limit of 180+/-20 km/s for the rotation velocity of this M_B(rest) ~ -21 galaxy. We note that our value is significantly higher than the rotation velocity of 120+/-10km/s reported by Barden et al. (2003) for their long-slit spectroscopic study of this galaxy. Our lower limit on the rotation velocity is entirely consistent with no evolution of the rest B-band Tully-Fisher relation. The position of this galaxy relative to the mean rest B-band Tully-Fisher relation of Tully & Pierce (2000) is consistent with brightening of no more than ~ 1mag at z=0.8. A larger integral field unit sample, without the uncertainties inherent to long-slit samples, is needed to accurately determine the evolution of the Tully-Fisher relation out to z ~ 1.

We employ panoramic, multicolour (BRz') and narrow-band Halpha imaging of the cluster Cl0024.0+1652 (z=0.39) from Subaru covering a ~30 arcmin field, to determine cluster membership and star formation rates for a large sample of galaxies across a wide field in the cluster, ~10 Mpc. We detect over 500 galaxies in narrow-band emission, with broad-band colours consistent with them lying at z~0.39. Using this sample we determine the Halpha luminosity function within the cluster and find that its form is approximately independent of local density, and is consistent with that seen in the intermediate redshift field population. This suggests that any density-dependent physical mechanisms which alter the star formation rate must leave the Halpha luminosity function unchanged; this is possible if the time-scale for star formation to cease completely is short compared with a Hubble time. Such short time-scale transformations are also supported by the presence of a population with late-type morphologies but no detectable Halpha emission. The fraction of blue galaxies, and the fraction of galaxies detected in Halpha, decreases strongly with increasing galaxy density in a manner which is qualitatively similar to that seen at lower redshifts. This trend is significantly steeper than the trend with galaxy morphology observed from a panoramic Hubble Space Telescope image of this cluster; this suggests that the physical mechanisms responsible for transformations in morphology and star formation rates may be partially independent. Finally, we compare our data with similar data on clusters spanning a range of redshifts from 0.2<~z<~0.8 and find little evidence for a trend in the total amount of star formation in clusters with redshift.

We analyse the u-r color distribution of 24346 galaxies with Mr<=-18 and z<0.08, drawn from the Sloan Digital Sky Survey first data release, as a function of luminosity and environment. The color distribution is well fit with two Gaussian distributions, which we use to divide the sample into a blue and red population. At fixed luminosity, the mean color of the blue (red) distribution is nearly independent of environment, with a weakly significant (~3sigma) detection of a trend for colors to become redder by 0.1-0.14 (0.03-0.06) mag with a factor ~100 increase in local density, as characterised by the surface density of galaxies within a +/-1000 km/s redshift slice. In contrast, at fixed luminosity the fraction of galaxies in the red distribution is a strong function of local density, increasing from ~10-30 per cent of the population in the lowest density environments, to ~70 per cent at the highest densities. The strength of this trend is similar for both the brightest (-23<Mr<-22) and faintest (-19<Mr<-18) galaxies in our sample. The fraction of red galaxies within the virialised regions of clusters shows no significant dependence on velocity dispersion. Even at the lowest densities explored, a substantial population of red galaxies exists, which might be fossil groups. We propose that most star-forming galaxies today evolve at a rate that is determined primarily by their intrinsic properties, and independent of their environment. Any environmentally triggered transformations from blue to red colors must either occur on a short timescale, or preferentially at high redshift, to preserve the simple Gaussian nature of the color distribution. The mechanism must be effective for both bright and faint galaxies.

We have used the SAURON panoramic integral field spectrograph to study the structure of the Ly-alpha emission-line halo, LAB1, surrounding the sub-millimeter galaxy SMM J221726+0013. This emission-line halo was discovered during a narrow-band imaging survey of the z=3.1 large-scale structure in the SSA22 region. Our observations trace the emission halo out to almost 100 kpc from the sub-millimeter source and identify two distinct Ly-alpha `mini-haloes' around the nearby Lyman-break galaxies. The main emission region has a broad line profile, with variations in the line profile seeming chaotic and lacking evidence for a coherent velocity structure. The data also suggests that Ly-alpha emission is suppressed around the sub-mm source. Interpretation of the line structure needs care because Ly-alpha may be resonantly scattered, leading to complex radiative transfer effects, and we suggest that the suppression in this region arises because of such effects. We compare the structure of the central emission-line halo with local counter-parts, and find that the emission line halo around NGC 1275 in the Perseus cluster may be a good local analogue, although the high redshift halo is factor of ~100 more luminous and appears to have higher velocity broadening. Around the Lyman-break galaxy C15, the emission line is narrower, and a clear shear in the emission wavelength is seen. A plausible explanation for the line profile is that the emission gas is expelled from C15 in a bipolar outflow, similar to that seen in M82.

In this review, we take the reader on a journey. We start by looking at the properties of galaxies in the cores of rich clusters. We have focused on the overall picture: star formation in clusters is strongly suppressed relative to field galaxies at the same redshift. We will argue that the increasing activity and blue populations of clusters with redshift results from a greater level of activity in field galaxies rather than a change in the transformation imposed by the cluster environment. With this in mind, we travel out from the cluster focussing firstly on the properties of galaxies in the outskirts of clusters and then on galaxies in isolated groups. At low redshift, we are able to efficiently probe these environments using the Sloan and 2DF redshift surveys. These allow an accurate comparison of galaxy star formation rates in different regions. The current results show a strong suppression of star formation above a critical threshold in local density. The threshold seems similar regardless of the overall mass of the system. At low redshift at least, only galaxies in close, isolated pairs have their star formation rate boosted above the global average. At higher redshift, work on constructing homogeneous catalogues of galaxies in groups and in the infall regions of clusters is still at an early stage. In the final section, we draw these strands together, summarising what we can deduce about the mechanisms that transform star forming field galaxies into their quiescent cluster counter-parts. We discuss what we can learn about the impact of environment on the global star formation history of the universe.

We analyse the observed correlation between galaxy environment and H-alpha emission line strength, using volume-limited samples and group catalogues of 24968 galaxies drawn from the 2dF Galaxy Redshift Survey (Mb<-19.5) and the Sloan Digital Sky Survey (Mr<-20.6). We characterise the environment by 1) Sigma_5, the surface number density of galaxies determined by the projected distance to the 5th nearest neighbour; and 2) rho1.1 and rho5.5, three-dimensional density estimates obtained by convolving the galaxy distribution with Gaussian kernels of dispersion 1.1 Mpc and 5.5 Mpc, respectively. We find that star-forming and quiescent galaxies form two distinct populations, as characterised by their H-alpha equivalent width, EW(Ha). The relative numbers of star-forming and quiescent galaxies varies strongly and continuously with local density. However, the distribution of EW(Ha) amongst the star-forming population is independent of environment. The fraction of star-forming galaxies shows strong sensitivity to the density on large scales, rho5.5, which is likely independent of the trend with local density, rho1.1. We use two differently-selected group catalogues to demonstrate that the correlation with galaxy density is approximately independent of group velocity dispersion, for sigma=200-1000 km/s. Even in the lowest density environments, no more than ~70 per cent of galaxies show significant H-alpha emission. Based on these results, we conclude that the present-day correlation between star formation rate and environment is a result of short-timescale mechanisms that take place preferentially at high redshift, such as starbursts induced by galaxy-galaxy interactions.

We have been undertaking a programme on the Gemini 8-m telescopes to demonstrate the power of integral field spectroscopy, using the optical GMOS spectrograph, and the new CIRPASS instrument in the near-infrared. Here we present some preliminary results from 3D spectroscopy of extra-galactic objects, mapping the emission lines in a 3CR radio galaxy and in a gravitationally lensed arc, exploring dark matter sub-structure through observations of an Einstein Cross gravitational lens, and the star formation time-scales of young massive clusters in the starburst galaxy NGC 1140.

We investigate the physical mechanisms that shape the luminosity function. Beginning with the mass function of dark matter halos, we show how gas cooling, photoionization, feedback, merging and thermal conduction affect the shape of the luminosity function. We consider three processes whereby supernovae can affect the forming galaxy: (1) reheating of disk gas to the halo temperature; (2) expansion of the diffuse halo gas; (3) expulsion of cold disk gas from the halo. While feedback of form (1) is able to flatten the faint end of the luminosity function, this alone does not produce the sharp cut-off observed at large luminosities. Feedback of form (2) is also unable to solve this problem. The relative paucity of very bright galaxies can only be explained if cooling in massive halos is strongly suppressed. Conduction is a promising mechanism, but an uncomfortably high efficiency is required to suppress cooling to the desired level. If, instead, superwinds are responsible for the lack of bright galaxies, then the total energy budget required to obtain a good match to the galaxy luminosity function greatly exceeds the energy available from supernova explosions. The mechanism is only viable if the formation of central supermassive black holes play a crucial role in limiting the amount of stars that form in the host galaxy. (abridged)

There are three major axes to describe the evolution of galaxies, namely, time (redshift), space (environment) and mass (stellar mass). In this article, one topic each will be presented along these axes. (1) Based on the Subaru wide-field (30') optical imaging of two distant clusters (z=0.55 and 0.83), we show the large scale structures in galaxy distributions on a scale greater than 10Mpc, which serve as an evidence for hierarchical growth of rich clusters of galaxies through assembly of surrounding groups. (2) Based on the deep NIR imaging of high-z clusters at z~1 and 1.5, we show massive galaxies in clusters have assembled most of their mass by z=1.5, which is earlier than the hierarchical model predictions. (3) Based on the Subaru deep and wide optical imaging of Subaru/XMM-Newton Deep Survey Field, we show a deficit of red and faint galaxies and a lack of blue massive galaxies in the high density regions at z~1, which suggest down-sizing in galaxy formation as apparently opposed to the CDM-based bottom-up scenario.

Data from a new, wide field, coincident optical and X-ray survey, the X-ray Dark Cluster Survey (XDCS) are presented. This survey comprises simultaneous and independent searches for clusters of galaxies in the optical and X-ray passbands. Optical cluster detection algorithms implemented on the data are detailed. Two distinct optically selected catalogues are constructed, one based on I-band overdensity, the other on overdensities of colour-selected galaxies. The superior accuracy of the colour-selection technique over that of the single passband method is demonstrated, via internal consistency checks and comparison with external spectroscopic redshift information. This is compared with an X-ray selected cluster catalogue. In terms of gross numbers, the survey yields 185 I-band selected, 290 colour selected and 15 X-ray selected systems, residing in ~11deg^2 of optical + X-ray imaging. The relationship between optical richness/ luminosity and X-ray luminosity is examined, by measuring X-ray luminosities at the positions of our 290 colour-selected systems. Power law correlations between the optical richness/ luminosity versus X-ray luminosity are fitted, both exhibiting approximately 0.2 dex of intrinsic scatter. Interesting outliers in these correlations are discussed in greater detail. Spectroscopic follow up of a subsample of X-ray underluminous systems confirms their reality.

We present the near-infrared (Ks-band) luminosity function of galaxies in two z~1 cluster candidates, 3C336 and Q1335+28. A third cluster, 3C289, was observed but found to be contaminated by a foreground system. Our wide field imaging data reach to Ks=20.5 (5sigma), corresponding to ~M*+2.7 with respect to the passive evolution. The near-infrared luminosity traces the stellar mass of a galaxy due to its small sensitivity to the recent star formation history. Thus the luminosity function can be transformed to the stellar mass function of galaxies using the $J-K$ colours with only a small correction (factor<2) for the effects of on-going star formation. The derived stellar mass function spans a wide range in mass from ~3 x 10^{11}Msun down to ~6 x 10^{9}Msun (set by the magnitude limit). The form of the mass function is very similar to lower redshift counterparts such as that from 2MASS/LCRS clusters (Balogh et al. 2001) and the z=0.31 clusters (Barger et al. 1998). This indicates little evolution of galaxy masses from z=1 to the present-day. Combined with colour data that suggest star formation is completed early (z>>1) in the cluster core, it seems that the galaxy formation processes (both star formation and mass assembly) are strongly accerelated in dense environments and has been largely completed by z=1. We investigate whether the epoch of mass assembly of massive cluster galaxies is earlier than that predicted by the hierarchical galaxy formation models. These models predict the increase of characteristic mass by more than factor ~3 between z=1 and the present day. This seems incompatible with our data.

The most massive galaxies in the present day universe are the giant ellipticals found in the centers of rich clusters. These have old, coeval stellar populations, suggesting they formed at high redshift, and are expected to host supermassive black holes (SMBHs). The recent detection of several high-redshift radio galaxies (HzRGs) at submm wavelengths confirms that indeed some massive galaxies may have formed the bulk of their stellar populations in spectacular dust-enshrouded starbursts at high-z. In this letter we compare sensitive Chandra X-ray images - which identify actively-fueled SMBHs - and submm observations - capable of detecting obscured activity in luminous galaxies at high redshift - of the environments of three HzRGs. These observations exhibit overdensities of X-ray sources in all 3 fields and a close correspondence between the Chandra and SCUBA populations. This suggests that both substantial star formation and nuclear activity may be occuring in these regions. We identify possible pairs of Chandra sources with each of 2 SCUBA sources, suggesting that their ultraluminous activity may be triggered by the interaction of two massive galaxies, each of which hosts an accreting SMBH. The presence of two SMBHs in the resulting remanent is predicted to produce a flattened stellar core in the galaxy, a morphological signature frequently seen in luminous cluster ellipticals. Hence the confirmation of pairs of Chandra sources within individual, luminous SCUBA galaxies would provide additional evidence that these galaxies at z~2-4 are the progenitors of the giant elliptical galaxies found in clusters at the present-day.

In this letter, we exploit the gravitational potential of the rich cluster A2218 as a magnifying glass. We demonstrate that the magnification due to the cluster allows us to observe distant background galaxies at a comparable level of detail to galaxies at z~0.1. Using the GMOS Integral Field Unit on Gemini North we observed the spatially-resolved [OII](3727) emission line spectrum for a lensed disk-galaxy at z=1.034. Using a detailed model for the cluster mass distribution, we are able to correct for the lensing by the cluster and reconstruct the source morphology. We find that the overall magnification is a factor of 4.92 +/- 0.15, and the rest-frame absolute I-band magnitude is M_{I}^{rest} = -22.4 +/- 0.2, where the error bars include conservative estimates of the uncertainty in the source-plane reconstruction. The inclination-corrected circular velocity is 206 +/- 18km/s. The galaxy lies very close to the mean Tully-Fisher relation of present-day spirals. Although our results are based on a single object, they demonstrate that gravitational lensing can be viably used to make detailed studies of the evolution of the structure of distant field galaxies.

We present sensitive, high-resolution, X-ray imaging from Chandra of the high-redshift radio galaxy 4C 41.17 (z=3.8). Our 150-ks Chandra exposure detects strong X-ray emission from a point source coincident with the nucleus of the radio galaxy. In addition we identify extended X-ray emission with a luminosity ~1e45 erg/s covering a 100kpc (15'') diameter region around the radio galaxy. The extended X-ray emission follows the general distribution of radio emission in the radio lobes of this source, and the distribution of a giant Lyman-alpha emission line halo, while the spectrum of the X-ray emission is non-thermal and has a power law index consistent with that of the radio synchrotron. We conclude that the X-ray emission is most likely Inverse-Compton scattering of far-infrared photons from a relativistic electron population probably associated with past and current activity from the central object. Assuming an equipartition magnetic field the CMB energy density at z=3.8 can only account for at most 40% of the Inverse-Compton emission. Published submillimeter maps of 4C 41.17 have detected an apparently extended and extremely luminous far-infrared emission around the radio galaxy. We demonstrate that this photon component and its spatial distribution, in combination with the CMB can reproduce the observed X-ray luminosity. We propose that photo-ionization by these Inverse-Compton X-ray photons plays a significant role in this system, and provides a new physical feedback mechanism to preferentially affect the gas within the most massive halos at high redshift. This is the highest redshift example of extended X-ray emission around a radio galaxy currently known. (Abridged)

The entropy distribution of the intracluster medium and the shape of its confining potential well completely determine the X-ray properties of a relaxed cluster of galaxies, motivating us to explore the origin of intracluster entropy and to describe how it develops in terms of some simple models. We present an analytical model for smooth accretion, including both preheating and radiative cooling, that links a cluster's entropy distribution to its mass accretion history and shows that smooth accretion overproduces the entropy observed in massive clusters by a factor ~2-3, depending on the mass accretion rate. Any inhomogeneity in the accreting gas reduces entropy production at the accretion shock; thus, smoothing of the gas accreting onto a cluster raises its entropy level. Because smooth accretion produces more entropy than hierarchical accretion, we suggest that some of the observed differences between clusters and groups may arise because preheating smooths the smaller-scale lumps of gas accreting onto groups more effectively than it smooths the larger-scale lumps accreting onto clusters. This effect may explain why entropy levels at the outskirts of groups are ~2-3 times larger than expected from self-similar scaling arguments. The details of how the density distribution of accreting gas affects the entropy distribution of a cluster are complex, and we suggest how to explore the relevant physics with numerical simulations.

To examine the evolution of the early-type galaxy population in the rich cluster Abell 2390 at z=0.23 we have gained spectroscopic data of 51 elliptical and lenticular galaxies with MOSCA at the 3.5 m telescope on Calar Alto Observatory. This investigation spans both a broad range in luminosity (-19.3>M_B>-22.3) and uses a wide field of view of 10'x10', therefore the environmental dependence of different formation scenarios can be analysed in detail as a function of radius from the cluster centre. Here we present results on the surface brightness modelling of galaxies where morphological and structural information is available in the F814W filter aboard the Hubble Space Telescope (HST) and investigate for this subsample the evolution of the Fundamental Plane.

We explore the evolution of the early-type galaxy population in the rich cluster Abell 2390 at z=0.23. For this purpose, we have obtained spectroscopic data of 51 elliptical and lenticular galaxies with MOSCA at the 3.5 m telescope on Calar Alto Observatory. As our investigation spans both a broad range in luminosity (-22.3<M_B<-19.3) and a wide field of view (10'x10'), the environmental dependence of different formation scenarios can be analysed in detail as a function of radius from the cluster center. In this paper, we present first results on the Faber-Jackson relation and, for a subsample of 14 galaxies with morphological and structural parameters from HST, we also investigate the evolution of the Kormendy relation and the Fundamental Plane. We find a mild luminosity evolution of the early-type galaxies in Abell 2390: our objects are on average brighter by m_B~0.4 mag.

The globally-averaged star formation rate in the Universe has been steadily declining since at least z~1. This may be due either to very local processes operating within the average galaxy, or to external, environmental effects. Specifically, the build-up of structure may be responsible for terminating star formation in some galaxies and thus decreasing the global average. We summarize our previous and ongoing work to distinguish between these possibilities, by determining the average star formation rate as a function of redshift and environment, out to z=0.5.

We present the first spectroscopic survey of intrinsically low Lx clusters at z>>0, with HST WFPC2 imaging and ground-based spectroscopy. We study 172 confirmed cluster members in a sample of ten clusters at z=0.23-0.3, with Lx<4.E43 h^{-2} ergs/s. The core of each cluster is imaged with WFPC2 in the F702W filter, and the spectroscopic sample is statistically complete to Mr-19.0+5log(h), within an 11' field. The clusters are dynamically well-separated from the surrounding field and the velocity dispersions range from ~350-850 km/s. Emission line galaxies in these clusters are relatively rare, comprising only 22 +/- 4% of the sample. There is no evidence that these emission-line galaxies are dynamically distinct from the majority of the cluster population, though our sample is too small to rule out the ~30% difference that has been observed in more massive clusters. We find eleven galaxies, comprising 6% of the cluster members, which are disk-dominated but show no sign of emission in their spectrum. Most of these are relatively isolated, spiral galaxies with smooth disks. We find no cluster members with a starburst or post-starburst spectrum. The striking similarity between the spectral and morphological properties of galaxies in these clusters and those of galaxies in more massive systems at similar redshifts implies that the physical processes responsible for truncating star formation in galaxies are not restricted to the rare, rich cluster environment, but are viable in much more common environments. In particular, we conclude that ram pressure stripping or cluster-induced starbursts cannot be solely responsible for the low star formation rates in these systems.

We report polarization measurements and photometry for the optical afterglow of the gamma-ray burst GRB 020405. We measured a highly significant 9.9% polarization (in V band) 1.3 days after the burst and argue that it is intrinsic to the GRB. The light curve decay is well fitted by a $t^{-1.72}$ power-law; we do not see any evidence for a break between 1.24 and 4.3 days after the burst. We discuss these measurements in the light of several models of GRB afterglows.

We present a set of cluster models that link the present-day properties of clusters to the processes that govern galaxy formation. These models treat the entropy distribution of the intracluster medium as its most fundamental property. Because convection strives to establish an entropy gradient that rises with radius, the observable properties of a relaxed cluster depend entirely on its dark-matter potential and the entropy distribution of its uncondensed gas. Guided by simulations, we compute the intracluster entropy distribution that arises in the absence of radiative cooling and supernova heating by assuming that the gas-density distribution would be identical to that of the dark matter. The lowest-entropy gas would then fall below a critical entropy threshold at which the cooling time equals a Hubble time. Radiative cooling and whatever feedback is associated with it must modify the entropy of that low-entropy gas, changing the overall entropy distribution function and thereby altering the observable properties of the cluster. Using some phenomenological prescriptions for entropy modification based on the existence of this cooling threshold, we construct a remarkably realistic set of cluster models. The surface-brightness profiles, mass-temperature relation, and luminosity-temperature relation of observed clusters all naturally emerge from these models. By introducing a single adjustable parameter related to the amount of intracluster gas that can cool within a Hubble time, we can also reproduce the observed temperature gradients of clusters and the deviations of cooling-flow clusters from the standard luminosity-temperature relation.

We construct a family of simple analytical models of galaxy clusters at the present epoch and compare its predictions with observational data. We explore two processes that break the self-similarity of galaxy clusters: systematic variation in the dark matter halo concentration and energy injection into the intracluster gas, through their effects on the observed cluster properties. Three observed relations between cluster properties and temperature are employed to constrain the model; mass, slope of gas density profile (beta) and luminosity. The slope of the mass-temperature relation is found to be reproduced by our model when the observed variation in concentration is included, raising the slope from the self-similar prediction of 1.5, to that of the observed relation, ~ 2. Heating of the gas is observed to have little effect on the mass-temperature relation. The mean trend in the beta-temperature relation is reproduced by energy injection of 0.5-0.75 keV per particle, while concentration variation has only a small effect. Excess energies calculated for individual systems from the beta-temperature relation suggest that the lowest mass systems may have excess energies that are biased to lower values by selection effects. The observed properties of the luminosity-temperature relation are reproduced by the combined effects of excess energy and a trend in the dark matter concentration. At high masses the observed variation in dark matter concentration results a slope of ~ 2.7 compared to recent observations in the range 2.6-2.9, whilst the observed steepening in galaxy groups is predicted when heating of 0.5-0.75 keV per particle is included. Hence a combination of energy injection and dark matter concentration variation appears able to account for the mean trends in the observed relations.

We have measured the equivalent width of the H-alpha emission line for 11006 galaxies brighter than M_b=-19 (LCDM) at 0.05<z<0.1 in the 2dF Galaxy Redshift Survey (2dF), in the fields of seventeen known galaxy clusters. The limited redshift range ensures that our results are insensitive to aperture bias, and to residuals from night sky emission lines. We use these measurements to trace mustar, the star formation rate normalized to Lstar, as a function of distance from the cluster centre, and local projected galaxy density. We find that the distribution of mustar steadily skews toward larger values with increasing distance from the cluster centre, converging to the field distribution at distances greater than ~3 times the virial radius. A correlation between star formation rate and local projected density is also found, which is independent of cluster velocity dispersion and disappears at projected densities below ~1 galaxy (brighter than M_b=-19) per Mpc^2. This characteristic scale corresponds approximately to the mean density at the cluster virial radius. The same correlation holds for galaxies more than two virial radii from the cluster centre. We conclude that environmental influences on galaxy properties are not restricted to cluster cores, but are effective in all groups where the density exceeds this critical value. The present day abundance of such systems, and the strong evolution of this abundance, makes it likely that hierarchical growth of structure plays a significant role in decreasing the global average star formation rate. Finally, the low star formation rates well beyond the virialised cluster rule out severe physical processes, such as ram pressure stripping of disk gas, as being completely responsible for the variations in galaxy properties with environment.

(Abridged) We present results of an H-alpha survey in the rich cluster A1689 at z=0.18, using the LDSS++ spectrograph on the AAT. We obtained spectra covering redshifted H-alpha for 522 galaxies brighter than I=22.5, covering a field of 8.7'X8.7'. We detect H-alpha emission in 46 of these galaxies; accounting for selection effects due to sampling and cluster membership, we determine that 24% of cluster members brighter than M_R=-16.5+5log(h) are detected with H-alpha flux greater than 4h^{-2}X 10^{38} ergs/s. This corresponds to a limiting star formation rate of 0.008 h^{-2} Msun/yr, assuming 1 magnitude of dust extinction. From a Hubble Space Telescope mosaic covering 7.5' X 10.0', we determine morphologies for 199 galaxies brighter than I=21, and find that 20% of the cluster members are of type Sa or later. More than 90% of cluster spirals show H-alpha emission, compared with less than 10% of E and S0 galaxies. The cluster H-alpha luminosity function has a low normalisation relative to the z~0.2 field, by ~50%, after accounting for the different fraction of spiral galaxies in the two environments. When compared with local field galaxies, this suggests that star formation activity is suppressed in early-type cluster galaxies, relative to their field counterparts. Our sample includes 29 galaxies previously observed with ISOCAM at 6 and 15 microns. We detect all 15-micron sources at H-alpha, so there is no evidence for any star formation completely hidden at H-alpha. Comparing the 15-micron and H-alpha fluxes, we find evidence that some mid-infrared-detected galaxies could be obscured by as much as 3 magnitudes of extinction at H-alpha, although this depends on the largely unknown contribution from any AGN-heated dust to the mid-infrared flux.

(Abridged) We present a morphological analysis of 17 X-ray selected clusters at z~0.25, imaged uniformly with HST WFPC2. Eight of these clusters have low X-ray luminosities (<10^{44} erg/s), while the remaining nine clusters have L_x>10^{45} ergs/s. The clusters cover a relatively small range in redshift and the data are homogeneous in terms of depth, resolution and rest wavelength observed. We use GIM2D to fit the two dimensional surface brightness profiles of galaxies down to M(702)<=-18.2 (Ho=50; roughly 0.01 L*) with parametric models, and quantify their morphologies using the fractional bulge luminosity (B/T). Within a single WFPC2 image we find that the Low-Lx clusters are dominated by galaxies with low B/T (~0), while the High-Lx clusters are dominated by galaxies with intermediate B/T (~0.4). We test whether this difference could arise from a universal morphology-density relation due to differences in the typical galaxy densities in the two samples. We find that small differences in the B/T distributions of the two samples persist with marginal statistical significance even when we restrict the comparison to galaxies in environments with similar projected local galaxy densities. From the correlations with the bulge and disk luminosity functions, we argue that the global environment affects the population of bulges, over and above trends seen with local density. We conclude that the destruction of disks through ram pressure stripping or harassment is not solely responsible for the morphology-density relation, and that bulge formation is less efficient in low mass clusters.

We present deep, panoramic multi-color imaging of the distant rich cluster A851 (Cl0939+4713, z=0.41) using Suprime-Cam on Subaru. These images cover a 27' field of view, ~11 Mpc at z=0.41, and by exploiting photometric redshifts estimated from our BVRI imaging we can isolate galaxies in a narrow redshift slice at the cluster redshift. Using a sample of ~2700 probable cluster members brighter than 0.02 Lv*, we trace the network of filaments and subclumps around the cluster core. The depth of our observations, combined with the identification of filamentary structure, gives us an unprecedented opportunity to test the influence of the environment on the properties of low luminosity galaxies. We find an abrupt change in the colors of < 0.1 Lv* galaxies at a local density of 100 gal. per sq. Mpc, with the population in lower density regions being predominantly blue, while those in higher density regions are red. The transition in the color-local density behavior occurs at densities corresponding to subclumps within the filaments surrounding the cluster. Identifying the sites where the transition occurs brings us much closer to understanding the mechanisms which are responsible for establishing the present-day relationship between environment and galaxy characteristics.

We present deep, panoramic multi-color imaging of the distant rich cluster A851 (z=0.41) using Suprime-Cam on Subaru. These images cover a 27' field of view (~11 h_{50}^{-1}Mpc), and by exploiting photometric redshifts, we can isolate galaxies in a narrow redshift slice at the cluster redshift. Using a sample of ~2700 probable cluster members (<M*+4), we trace the network of filaments and subclumps around the cluster core. The depth of our observations, combined with the identification of filamentary structure, gives us an unprecedented opportunity to test the influence of the environment on the properties of low luminosity galaxies. We find an abrupt change in the colors of faint galaxies (>M*+1) at a local density of 100 gal. Mpc^{-2}. The transition in the color-local density behavior occurs at densities corresponding to subclumps within the filaments surrounding the cluster. Identifying the sites where the transition occurs brings us much closer to understanding the mechanisms which are responsible for establishing the present-day relationship between environment and galaxy characteristics.

We use a three dimensional hydrodynamical code to simulate the effect of energy injection on cooling flows in the intracluster medium. Specifically, we compare a simulation of a 10$^{15}$ $M_\odot$ cluster with radiative cooling only, with a second simulation in which thermal energy is injected 31 $\kpc$ off-centre, over 64 kpc$^{3}$ at a rate of $4.9\times 10^{44} \ergs$ for 50 Myr. The heat injection forms a hot, low density bubble which quickly rises, dragging behind it material from the cluster core. The rising bubble pushes with it a shell of gas which expands and cools. We find the appearance of the bubble in X-ray temperature and luminosity to be in good qualitative agreement with recent {\it Chandra} observations of cluster cores. Toward the end of the simulation, at 600 Myr, the displaced gas begins to fall back toward the core, and the subsequent turbulence is very efficient at mixing the low and high entropy gas. The result is that the cooling flow is disrupted for up to $\sim 50$ Myr after the injection of energy ceases. Thus, this mechanism provides a very efficient method for regulating cooling flows, if the injection events occur with a 1:1 duty cycle.

We describe an instrument concept that will provide simultaneous spectra for a million spatial samples on the sky. With the proposed field of view and spectral resolution, it will be able to measure redshifts and line strengths for around 2-4000 z~3-7 galaxies in a 16 night campaign. The main science driver is to obtain a complete census of the star formation properties of galaxies with 2.5<z<6.7 as a function of luminosity, environment, morphology and redshift. This survey will also allow us to study the evolution of the 3 dimensional power spectrum and test the connection between these early galaxies and the large scale structure of the Lyman Alpha forest. Additional science drivers include [OII] emitters with 0.13<z<1.52, abundance gradients in nearby galaxies, searches for planetary nebulae in nearby galaxy groups and clusters, lensed emission line objects, and spectroscopic searches for objects with continuum breaks. We expect that this instrument will make an ideal partner to the HST ACS, where ACS morphology and broadband colour information can complement the MEIFU spectroscopy and help remove any degeneracies in line identification.

We present a new photometric catalogue of the Coma galaxy cluster in the Johnson U- and V- bands. We cover an area of 3360arcmin2 of sky, to a depth of V=20 mag in a 13 arcsec diameter aperture, and produce magnitudes for ~1400 extended objects in metric apertures from 8.8 to 26arcsec diameters. The mean internal RMS scatter in the photometry is 0.014 mag in V, and 0.026 mag in U, for V13 < 17 mag. We place new limits on the levels of scatter in the colour--magnitude relation (CMR) in the Coma cluster, and investigate how the slope and scatter of the CMR depends on galaxy morphology, luminosity and position within the cluster. As expected, the lowest levels of scatter are found in the elliptical galaxies, while the late type galaxies have the highest numbers of galaxies bluewards of the CMR. We investigate whether the slope of the CMR is an artifact of colour gradients within galaxies and, show that it persists when the colours are measured within a diameter that scales with galaxy size. Looking at the environmental dependence of the CMR, we find a trend of systematically bluer galaxy colours with increasing projected cluster-centric radius which we associate with a gradient in the mean galactic ages.

Recent measurements of the K-band luminosity function now provide us with strong, reliable constraints on the fraction of baryons which have cooled. Globally, this fraction is only about 5%, and there is no strong evidence that it is significantly higher in clusters. Without an effective sub-grid feedback prescription, the cooled gas fraction in any numerical simulation exceeds these observational constraints, and increases with increasing resolution. This compromises any discussion of galaxy and cluster properties based on results of simulations which include cooling but do not implement an effective feedback mechanism.

(Abridged) We present morphological, photometric, kinematic and line-strength data of 48 early-type members of the rich cluster Abell 2218 at z=0.18. This sample is both larger and spans a wider galaxy luminosity range, down to M*_B+1, than previous studies. We compare A2218 with similar samples drawn from local clusters to identify evolutionary changes between the samples which have occured over the last 3Gyrs. The overall picture is one of little or no evolution in nearly all galaxy parameters. Zeropoint offsets in the Faber--Jackson, Mg_b-sigma and Fundamental Plane relations are all consistent with passively evolving stellar populations. The slopes of these relations have not changed significantly. We do however find a significant spread in the estimated luminosity-weighted ages of the stellar populations in the galaxies, based on line diagnostic diagrams. This age spread is seen in both the disky early-type galaxies (S0) and also the ellipticals. We observe both ellipticals with a strong contribution from a young stellar population and lenticulars dominated by old stellar populations. On average, we find no evidence for systematic differences between the populations of ellipticals and lenticulars. This points to a common formation epoch for the bulk of the stars in most of the early-type galaxies in A2218. This result can be reconciled with the claims of rapid morphological evolution in distant clusters if the suggested transformation from spirals to lenticulars does not involve significant new star formation.

We present a new model for the X-ray properties of the intracluster medium that explicitly includes heating of the gas by the energy released during the evolution of cluster galaxies. We calculate the evolution of clusters by combining the semi-analytic model of galaxy formation of Cole et al. with a simple model for the radial profile of the intracluster gas. We focus on the cluster X-ray luminosity function and on the relation between X-ray temperature and luminosity. We show that cooling alone is not enough to account for the flatness of the observed T-L relation or for the lack of strong redshift evolution in the observed X-ray luminosity function. Gas heating, on the other hand, can solve these two problems: in the \Lambda-CDM cosmology, our model reproduces fairly well the T-L relation and the X-ray luminosity function. Furthermore, it predicts only weak evolution in these two properties out to z=0.5, in agreement with recent observational data. A successful model requires an energy input of 1--2 x 10^49 ergs per solar mass of stars formed. This is comparable to the total energy released by the supernovae associated with the formation of the cluster galaxies. Thus, unless the transfer of supernovae energy to the intracluster gas is very (perhaps unrealistically) efficient, additional sources of energy, such as mechanical energy from AGN winds are required. However, the amplification of an initial energy input by the response of the intracluster medium to protocluster mergers might ease the energy requirements. Our model makes definite predictions for the X-ray properties of groups and clusters at high redshift which may soon be tested with data from the Chandra and Newton satellites.

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.

The observed morphological evolution since z ~ 0.5 within galaxy clusters provides evidence for a hierarchical universe. This evolution is driven by dynamical effects that operate within the cluster environment -- suppression of star-formation by ram-pressure and viscous stripping of the intra-galactic medium and tidal heating of disks by gravitational encounters.

We address the current crucial issues on the formation and evolution of cluster galaxies: ie., connection between the Butcher-Oemler effect, assembly of cluster galaxies, truncation of star formation, and the origin of S0 galaxies. We construct the field corrected colour-magnitude (CM) diagrams for 7 CNOC clusters (0.23<z<0.43) and Coma, and illustrate the evolution of the complete cluster population down to the present-day based on the model in which star formation is truncated when the galaxies infall from the surrounding field. We show that the blue galaxies are incorporated into the present-day tight CM relation as they fade and become redder after the truncation, which is possibly responsible for producing faint S0's (>M*+1). Truncation of star formation is, however, found to be relatively milder (with a time scale of 1 Gyr) than suggested by the viorent processes such as ram-pressure stripping and/or mergers/harassment. The BO effect is after all found to be a combination of three effects, namely, increasing field star formation activity, increasing galaxy infall rate, and the truncatin of star formation after the accretion. Our approach naturally leads to the history of galaxy assembly and `global' star foramtion for `cluster' galaxies.

We present the first 3-dimensional high resolution hydro-dynamical simulations of the interaction between the hot ionised intra-cluster medium and the cold interstellar medium of spiral galaxies. Ram pressure and turbulent/viscous stripping removes 100% of the atomic hydrogen content of luminous galaxies like the Milky Way within 100 million years. These mechanisms naturally account for the morphology of S0 galaxies, the rapid truncation of star formation implied by spectroscopic observations, as well as a host of observational data on the HI morphology of galaxies in clusters.

We have constructed a family of simple models for spiral galaxy evolution to allow us to investigate observational trends in star formation history with galaxy parameters. The models are used to generate broad band colours from which ages and metallicities are derived in the same way as the data. We generate a grid of model galaxies and select only those which lie in regions of parameter space covered by the sample. The data are consistent with the proposition that the star formation history of a region within a galaxy depends primarily on the local surface density of the gas but that one or two additional ingredients are required to fully explain the observational data. The observed age gradients appear steeper than those produced by the density dependent star formation law, indicating that the star formation law or infall history must vary with galactocentric radius. Furthermore, the metallicity-magnitude and age-magnitude correlations are not reproduced by a local density dependence alone. These correlations require one or both of the following: (i) a combination of mass dependent infall and metal enriched outflow, or (ii) a mass dependent galaxy formation epoch. Distinguishing these possibilities on the basis of current data is extremely difficult.

Current theories of galaxy formation predict that spiral galaxies are embedded in a reservoir of hot gas. This gas is able to cool onto the galaxy replenishing cold gas that is consumed by star formation. Estimates of the X-ray luminosity emitted in the cooling region suggest a bolometric luminosity of order 10 x 10^41 ergs/s in massive systems. We have used ROSAT PSPC data to search for extended X-ray emission from the haloes of three nearby, massive, late-type galaxies: NGC 2841, NGC 4594 and NGC 5529. We infer 95 per cent upper limits on the bolometric X-ray luminosities of the haloes of NGC 2841, NGC 4594 and NGC 5529 of 0.4, 1.2 and 3.8 x 10^41 ergs/s respectively. Thus the true luminosity lies well below the straightforward theoretical prediction. We discuss this discrepancy and suggest a number of ways in which the theoretical model can be brought into agreement with the observational results. A possible solution is that the gravitational potentials of the dark matter haloes of these galaxies are weaker than assumed in the current model. Alternatively, the present day accretion may be substantially less than is required on average to build the disk over the Hubble time. Our results are, however, based on only three galaxies, none of which is ideal for this kind of study. A larger data set is required to explore this important problem further.

We have performed deep imaging of a diverse sample of 26 low surface brightness galaxies (LSBGs) in the optical and the near-infrared. Using stellar population synthesis models, we find that it is possible to place constraints on the ratio of young to old stars (which we parameterise in terms of the average age of the galaxy), as well as the metallicity of the galaxy, using optical and near-infrared colours. LSBGs have a wide range of morphologies and stellar populations, ranging from older, high metallicity earlier types to much younger and lower metallicity late type galaxies. Despite this wide range of star formation histories, we find that colour gradients are common in LSBGs. These are most naturally interpreted as gradients in mean stellar age, with the outer regions of LSBGs having younger ages than their inner regions. In an attempt to understand what drives the differences in LSBG stellar populations, we compare LSBG average ages and metallicities with their physical parameters. Strong correlations are seen between a LSBG's star formation history and its K band surface brightness, K band absolute magnitude and gas fraction. These correlations are consistent with a scenario in which the star formation history of a LSBG primarily correlates with its surface density and its metallicity correlates both with its mass and surface density.

We use the $C_{2}4668$, $Fe4383$, $H\gamma_A$ and $H\delta_A$ spectral absorption line indices, together with U- and V-band photometry of 101 galaxies in the Coma cluster, to investigate how mean age and metal abundance correlate with galaxy luminosity. In particular, we use the line index measurements to address the origin of the colour-magnitude relation (CMR). We find that the CMR in Coma is driven primarily by a luminosity-metallicity correlation. We additionally show evidence for a relation between age and luminosity, in the direction predicted by the semi-analytic hierarchical clustering models of Kauffmann & Charlot, but this is only present in the $C_24668$ index models, and could be an effect of the lack of non solar abundance ratios in the Worthey models used. By comparing deviations from the CMR, with deviations in absorption index from analogous `index-magnitude' relations, we find that colour deviations bluewards of the mean relation are strongly correlated with the hydrogen Balmer line series absorption. We show that the properties of these blue galaxies are consistent with the presence of a young stellar population in the galaxies, rather than with a reduced metallicity.

We present results of an extensive study of the X-ray spectral properties of sources detected in the RIXOS survey, that is nearly complete down to a flux limit of 3e-14 cgs (0.5-2 keV). We show that for X-ray surveys containing sources with low count rate spectral slopes estimated using simple hardness ratios in the ROSAT band can be biased. Instead we analyse three-colour X-ray data using statistical techniques appropriate to the Poisson regime which removes the effects of this bias. We have then applied this technique to the RIXOS survey to study the spectral properties of the sample. For the AGN we find an average energy index of 1.05+-0.05 with no evidence for spectral evolution with redshift. Individual AGN are shown to have a range of properties including soft X-ray excesses and intrinsic absorption. Narrow Emission Line Galaxies also seem to fit to a power-law spectrum, which may indicate a non-thermal origin for their X-ray emission. We infer that most of the clusters in the sample have a bremsstrahlung temperature >3 keV, although some show evidence for a cooling flow. The stars deviate strongly from a power-law model but fit to a thermal model. Finally, we have analysed the whole RIXOS sample containing 1762 sources. We find that the mean spectral slope of the sources hardens at lower fluxes in agreement with results from other samples. However, a study of the individual sources demonstrates that the hardening of the mean is caused by the appearance of a population of very hard sources at the lowest fluxes. This has implications for the nature of the soft X-ray background.

We describe similarity solutions that characterize the collapse of collisional gas onto scalefree perturbations in an Einstein-de Sitter universe. We consider the effects of radiative cooling and derive self-similar solutions under the assumption that the cooling function is a power law of density and temperature, $\Lambda(T,\rho) \propto \rho^{3/2} T$. We use these results to test the ability of Smooth Particle Hydrodynamics (SPH) techniques to follow the collapse and accretion of shocked, rapidly cooling gas in a cosmological context. Our SPH code reproduces the analytical results very well in cases that include or exclude radiative cooling. No substantial deviations from the predicted central mass accretion rates or from the temperature, density, and velocity profiles are observed in well resolved regions inside the shock radius. This test problem lends support to the reliability of SPH techniques to model the complex process of galaxy formation.

In this paper, we have extended the entropy-driven model of cluster evolution developed by Bower (1997) in order to be able to predict the evolution of galaxy clusters for a range of cosmological scenarios. We have applied this model to recent measurements of the evolution of the L_x-T normalisation and X-ray luminosity function in order to place constraints on cosmological parameters. We find that these measurements alone do not select a particular cosmological frame-work. An additional constraint is required on the effective slope of the power spectrum to break the degeneracy that exists between this and the background cosmology. We have therefore included a theoretical calculation of the \Omega_0 dependence on the power spectrum, based on the cold dark matter paradigm, which infers \Omega_0<0.55 (0.1<\Omega_0<0.7 for \Omega_0+\Lambda_0=1), at the 95% confidence level. Alternatively, an independent measurement of the slope of the power spectrum from galaxy clustering requires \Omega_0<0.6 (\Omega_0<0.65 for \Omega_0+\Lambda_0=1), again to 95% confidence. The rate of entropy evolution is insensitive to the values of \Omega_0 considered, although is sensitive to changes in the distribution of the intracluster medium.

We use 3-dimensional SPH/N-BODY simulations to study ram pressure stripping of gas from spiral galaxies orbiting in clusters. We find that the analytic expectation of Gunn & Gott (1972) relating the gravitational restoring force provided by the disk to the ram pressure force, provides a good approximation to the radius that gas will be stripped from a galaxy. However, at small radii it is also important to consider the potential provided by the bulge component. A spiral galaxy passing through the core of a rich cluster such as Coma, will have its gaseous disk truncated to $\sim 4$ kpc, thus losing $\sim 80%$ of its diffuse gas mass. The timescale for this to occur is a fraction of a crossing time $\sim 10^7$ years. Galaxies orbiting within poorer clusters, or inclined to the direction of motion through the intra-cluster medium will lose significantly less gas. We conclude that ram-pressure alone is insufficient to account for the rapid and widespread truncation of star-formation observed in cluster galaxies, or the morphological transformation of Sab's to S0's that is necessary to explain the Butcher-Oemler effect.

Near-infrared (NIR) K' images of a sample of five low surface brightness disc galaxies (LSBGs) were combined with optical data, with the aim of constraining their star formation histories. Both red and blue LSBGs were imaged to enable comparison of their stellar populations. For both types of galaxy strong colour gradients were found, consistent with mean stellar age gradients. Very low stellar metallicities were ruled out on the basis of metallicity-sensitive optical-NIR colours. These five galaxies suggest that red and blue LSBGs have very different star formation histories and represent two independent routes to low B band surface brightness. Blue LSBGs are well described by models with low, roughly constant star formation rates, whereas red LSBGs are better described by a `faded disc' scenario.

Photometric redshift classifiers provide a means of estimating galaxy redshifts from observations using a small number of broad-band filters. However, the accuracy with which redshifts can be determined is sensitive to the star formation history of the galaxy, for example the effects of age, metallicity and on-going star formation. We present a photometric classifier that explicitly takes into account the degeneracies implied by these variations, based on the flexible stellar population synthesis code of Kodama & Arimoto. The situation is encouraging since many of the variations in stellar populations introduce colour changes that are degenerate. We use a Bayesian inversion scheme to estimate the likely range of redshifts compatible with the observed colours. When applied to existing multi-band photometry for Abell 370, most of the cluster members are correctly recovered with little field contamination. The inverter is focussed on the recovery of a wide variety of galaxy populations in distant (z~1) clusters from broad band colours covering the 4000 angstrom break. It is found that this can be achieved with impressive accuracy ($|\Delta z| < 0.1$), allowing detailed investigation into the evolution of cluster galaxies with little selection bias.

We present the rest-frame colour-magnitude diagram for 35 early-type galaxies in the Hubble Deep Field with median redshift 0.9. Although with considerable scatter, a red sequence well-described by the passive evolution of an intrinsically old stellar population is observed. Comparison with the passively evolved colour-magnitude relation of the rich Coma cluster ($z=0.023$) indicates that around half of the early-type galaxies in the field at this redshift are as old as those in rich clusters.

This talk investigates the formation of early-type galaxies from a deliberately observational view point. I begin by reviewing the conclusions that can be reached by comparing the detailed properties of galaxies in present-day clusters, focusing on the colour-magnitude relation in particular. The overriding picture is one of homogeneity, implying a remarkable uniformity in the formation of these galaxies. This picture contrasts with the increasing activity seen in clusters as a function of redshift, creating an apparent paradox been the obvious diversity of star formation histories in distant cluster galaxies and their uniformity in local systems. A resolution is feasible so long as star formation occurs over an extended epoch. In addition to placing limits on variations in star formation history the existence of a tight `fundamental relations', such as the colour-magnitude relation, can be used to investigate galaxy mergers and to set limits on the degree to which present-day clusters galaxies are built by combining systems of stars formed in smaller units. The final part of this talk turns to early-type galaxies in the field, and tries to apply the same techniques that have been successful in clusters. This is an emerging field in which appropriate data-sets are only just becoming available; however, comparison of the formation histories of galaxies in a wide variety of environments is key to distinguishing between the Classical and Hierarchical models for galaxy formation.

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.

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-\sigma$ 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.

Observations of the evolution of the galaxy cluster X-ray luminosity function suggest that the entropy of the intra-cluster medium plays a significant role in determining the development of cluster X-ray properties. I present a theoretical framework in which the evolution of the entropy of the central intra-cluster gas is explicitly taken into account. The aim of this work is to develop a theoretical context within which steadily improving measurements of the X-ray luminosities and temperatures of distant galaxy clusters can be interpreted. I discuss the possible range of entropy evolution parameters and relate these to the physical processes heating and cooling the intra-cluster medium. The practical application of this work is demonstrated by combining currently available evolutionary constraints on the X-ray luminosity function and the luminosity--temperature correlation to determine the best-fitting model parameters.

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.

These notes collect together current work on the effect of the environment on galaxy formation and evolution. They are broken into four distinct parts. The first deals with the observational debate surrounding the question of whether galaxies are initially created the same and have only recently been modified by their environment, or whether galaxy formation itself is strongly perturbed by the local neighbourhood. The second section deals with observations of distant galaxy clusters. These provide direct evidence that some galaxies, at least, have changed substantially over the last 5~Gyr. The third and fourth parts look at theoretical work on this theme. They study models for the origin of the density--morphology relation, and the impact of environmentally modulated galaxy formation on the large-scale structure of the universe.

Clusters of galaxies are the largest gravitationally bound systems and therefore provide an important way of studying the formation and evolution of the large scale structure of the Universe. Cluster evolution can be inferred from observations of the X-ray emission of the gas in distant clusters, but interpreting these data is not straightforward. In a simplified view, clusters grow from perturbations in the matter distribution: their intracluster gas is compressed and shock-heated by the gravitational collapse$^{1}$. The resulting X-ray emission is determined by the hydrostatic equilibrium of the gas in the changing gravitational potential. However, if processes such as radiative cooling or pre-collapse heating of the gas are important, then the X-ray evolution will be strongly influenced by the thermal history of the gas. Here we present the first results from a faint flux-limited sample of X-ray selected clusters compiled as part of the ROSAT International X-ray and Optical Survey (RIXOS). Very few distant clusters have been identified. Most importantly, their redshift distribution appears to be inconsistent with simple models based on the evolution of the gravitational potential. Our results suggest that radiative cooling or non-gravitational heating of the intracluster gas must play an important role in the evolution of clusters.