Search

We exploit the unprecedented depth of integral field data from the KMOS Ultra-deep Rotational Velocity Survey (KURVS) to analyse the strong (H$\alpha$) and forbidden ([NII], [SII]) emission line ratios in 22 main-sequence galaxies at $z\approx1.5$. Using the [NII]/H$\alpha$ emission-line ratio we confirm the presence of the stellar mass $-$ gas-phase metallicity relation at this epoch, with galaxies exhibiting on average 0.13$\pm$0.04 dex lower gas-phase metallicity (12+log(O/H)$_{\rm M13}$=8.40$\pm$0.03) for a given stellar mass ($\log_{10}$($M_{\rm *}$[$M_{\odot}$]=10.1$\pm$0.1) than local main-sequence galaxies. We determine the galaxy-integrated [SII] doublet ratio, with a median value of [SII]$\lambda$6716/$\lambda$6731=1.26$\pm$0.14 equivalent to an electron density of log$_{10}$($n_{\rm e}$[cm$^{-3}$])=1.95$\pm$0.12. Utilising CANDELS $HST$ multi-band imaging we define the pixel surface-mass and star-formation rate density in each galaxy and spatially resolve the fundamental metallicity relation at $z\approx1.5$, finding an evolution of 0.05$\pm$0.01 dex compared to the local relation. We quantify the intrinsic gas-phase metallicity gradient within the galaxies using the [NII]/H$\alpha$ calibration, finding a median annuli-based gradient of $\Delta$Z/$\Delta$R=$-$0.015$\pm$0.005 dex kpc$^{-1}$. Finally we examine the azimuthal variations in gas-phase metallicity, which show a negative correlation with the galaxy integrated star-formation rate surface density ($r_{\rm s}$=$-$0.40, $p_{\rm s}$=0.07) but no connection to the galaxies kinematic or morphological properties nor radial variations in stellar mass surface density or star formation rate surface density. This suggests both the radial and azimuthal variations in interstellar medium properties are connected to the galaxy integrated density of recent star formation., Comment: Submitted for publication in MNRAS. Comments welcome. 18 pages

In this paper we present data from 72 low redshift, hard X-ray selected AGN taken from the {\it Swift}-BAT 58 month catalogue. We utilise spectral energy distribution fitting to the optical to IR photometry in order to estimate host galaxy properties. We compare this observational sample to a volume and flux matched sample of AGN from the EAGLE hydrodynamical simulations in order to verify how accurately the simulations can reproduce observed AGN host galaxy properties. After correcting for the known +0.2 dex offset in the SFRs between EAGLE and previous observations, we find agreement in the SFR and X-ray luminosity distributions; however we find that the stellar masses in EAGLE are $0.2 - 0.4$ dex greater than the observational sample, which consequently leads to lower sSFRs. We compare these results to our previous study at high redshift, finding agreement in both the observations and simulations, whereby the widths of sSFR distributions are similar ($\sim0.4-0.6$ dex) and the median of the SFR distributions lie below the star forming main sequence by $\sim0.3-0.5$ dex across all samples. We also use EAGLE to select a sample of AGN host galaxies at high and low redshift and follow their characteristic evolution from $z=8$ to $z=0$. We find similar behaviour between these two samples, whereby star formation is quenched when the black hole goes through its phase of most rapid growth. Utilising EAGLE we find that 23\% of AGN selected at $z\sim0$ are also AGN at high redshift, and that their host galaxies are among the most massive objects in the simulation. Overall we find EAGLE reproduces the observations well, with some minor inconsistencies ($\sim$ 0.2 dex in stellar masses and $\sim$ 0.4 dex in sSFRs)., Comment: 14 pages, 6 figures, 3 tables, accepted for publication in MNRAS

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., Comment: Published in MNRAS. This is the author accepted version. A small number of typos were corrected after journal proof stage

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., Comment: Accepted for publication in MNRAS. Comments welcome. 25 pages, 14 figures. The full pdf, which includes an additional 17 pages of supplementary figures is available on this http://astro.dur.ac.uk/KROSS/publications.html

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., Comment: Accepted for publication in MNRAS. 14 pages, 9 figures, 3 tables

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., Comment: 31 pages, 18 figures including appendices. Accepted for publication in MNRAS 2019 March 12. Updated text and figures

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., Comment: 32 pages, 21 figures including appendices. Accepted for publication in MNRAS after addressing final comments from referee. Additional appendices added, plus some updates and additions to plots and text

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., Comment: Accepted for publication in MNRAS

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.680.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}$., Comment: 22 pages, 18 figures, 4 tables, submitted to MNRAS 14th Sept 2018

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}$., Comment: Accepted for publication in MNRAS

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., Comment: 18 pages, 12 figures, accepted for publication in MNRAS

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., Comment: 19 pages, 11 figures, accepted to ApJ

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., Comment: 21 pages, 12 figures, accepted for publication in MNRAS

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., Comment: Accepted for publication in MNRAS. Data available at http://astro.dur.ac.uk/KROSS

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., Comment: Accepted for publication in MNRAS. Table of values and extended figures available at: http://astro.dur.ac.uk/KROSS/

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., Comment: Accepted for publication in MNRAS. Main paper: 12 pages, 7 figures. Plus appendix

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., Comment: Accepted for publication in MNRAS. Main paper: 20 pages, 14 figures. Plus appendix

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., Comment: 14 pages, 10 figures

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., Comment: 26 pages, 11 figures. Accepted for publication in MNRAS

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., Comment: Accepted for publication by MNRAS; 18 pages, 7 figures; comments welcome

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?, Comment: Accepted to ApJ

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), Comment: 19 pages, 10 figures, ApJ accepted

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., Comment: 24 pages, 14 figures. Accepted to MNRAS for publication. Updated version with typos corrected

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., Comment: Supporting paper for the science theme The Hot and Energetic Universe to be implemented by the Athena+ X-ray observatory (http://www.the-athena-x-ray-observatory.eu). 9 pages, 7 figures

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., Comment: 14 pages, 4 figures, published in ApJL

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., Comment: 16 pages, 9 figures, 2 tables, accepted for publication in MNRAS

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., Comment: 27 pages (31 color figures). Accepted for publication in MNRAS. E-mail to obtain Model catalogues

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., Comment: 16 pages, 13 figures. Accepted to MNRAS

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), Comment: 22 pages, 17 figures. Accepted for publication in MNRAS. ORCA cluster catalogues available at http://orca.dur.ac.uk/

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., Comment: Accepted to ApJ. 44 pages, 13 figures

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., Comment: 23 pages, 17 figures, submitted to MNRAS

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

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., Comment: 33 pages, 15 figures. Accepted by MNRAS

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., Comment: 9 pages, 5 figures, accepted for publication in MNRAS

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., Comment: Accepted for publication in MNRAS after correcting for revisions suggested by the referee

(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., Comment: 14 pages, 9 figures, A&A, in press. Final version corrects Eqn. 3 and columns 5-8 in Table 1 and has additional minor text clarifications

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., Comment: 12 pages, 12 figures. Accepted for publication in A&A

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., Comment: To appear in proceedings of the conference "The Monster's Fiery Breath: Feedback in Galaxies, Groups, and Clusters", June 2009, Madison Wisconsin

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., Comment: 15 pages, ApJ in press

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., Comment: ApJL, in press

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., Comment: 9 pages, 4 figures. Accepted for publication in the Astrophysical Journal. This version updated to accepted version

(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., Comment: 19 pages, 10 figures, Accepted for publication in MNRAS

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., Comment: 13 pages, 6 figures. Accepted for publication in ApJ

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., Comment: Accepted for publication in MNRAS Letters

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., Comment: 15 pages, 6 figures. Accepted for publication in ApJ

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%., Comment: 20 pages, 21 figures. Accepted for publication in A&A. v3 uses updated data to correct some minor mistakes, text is unchanged

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., Comment: 19 pages; MNRAS in press

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., Comment: 6 Pages, 3 figures, Submitted to MNRAS

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., Comment: 17 pages, 17 figures; to appear in A&A. A pdf version with full-quality figures can be found at ftp://ftp.xray.mpe.mpg.de/people/gwp/xmmlp/xmmlp.pdf

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., Comment: Accepted for publication in MNRAS