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.

We present a study of the galaxy environment of 9 strong HI+CIV absorption line systems ($16.2<{\rm log}(N({\rm HI}))<21.2$) spanning a wide range in metallicity at $z\sim4-5$, using MUSE integral field and X-Shooter spectroscopic data collected in a $z\approx 5.26$ quasar field. We identify galaxies within a 250 kpc and $\pm1000$ km s$^{-1}$ window for 6 out of the 9 absorption systems, with 2 of the absorption line systems showing multiple associated galaxies within the MUSE field of view. The space density of Ly$\alpha$ emitting galaxies (LAEs) around the HI and CIV systems is $\approx10-20$ times the average sky density of LAEs given the flux limit of our survey, showing a clear correlation between the absorption and galaxy populations. Further, we find that the strongest CIV systems in our sample are those that are most closely aligned with galaxies in velocity space, i.e. within velocities of $\pm500$ km s$^{-1}$. The two most metal poor systems lie in the most dense galaxy environments, implying we are potentially tracing gas that is infalling for the first time into star-forming groups at high redshift. Finally, we detect an extended Ly$\alpha$ nebula around the $z\approx 5.26$ quasar, which extends up to $\approx50$ kpc at the surface brightness limit of $3.8 \times 10^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$. After scaling for surface brightness dimming, we find that this nebula is centrally brighter, having a steeper radial profile than the average for nebulae studied at $z\sim3$ and is consistent with the mild redshift evolution seen from $z\approx 2$.

We present the design, methods, and first results of the MUSE Analysis of Gas around Galaxies (MAGG) survey, a large programme on the Multi Unit Spectroscopic Explorer (MUSE) instrument at the Very Large Telescope (VLT) which targets 28 z > 3.2 quasars to investigate the connection between optically-thick gas and galaxies at z~3-4. MAGG maps the environment of 52 strong absorption line systems at z > 3, providing the first statistical sample of galaxies associated with gas-rich structures in the early Universe. In this paper, we study the galaxy population around a very metal poor gas cloud at z~3.5 towards the quasar J124957.23-015928.8. We detect three Lyman alpha emitters within <200km/s of the cloud redshift, at projected separations <185 kpc (physical). The presence of star-forming galaxies near a very metal-poor cloud indicates that metal enrichment is still spatially inhomogeneous at this redshift. Based on its very low metallicity and the presence of nearby galaxies, we propose that the most likely scenario for this LLS is that it lies within a filament which may be accreting onto a nearby galaxy. Taken together with the small number of other LLSs studied with MUSE, the observations to date show a range of different environments near strong absorption systems. The full MAGG survey will significantly expand this sample and enable a statistical analysis of the link between gas and galaxies to pin down the origin of these diverse environments at z~3-4.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

We present first results from Multi Unit Spectroscopic Explorer (MUSE) observations at the Very Large Telescope in the MUSE Ultra Deep Field (MUDF), a $\approx 1.2\times 1.4$ arcmin$^2$ region for which we are collecting $\approx$200 hours of integral field spectroscopy. The $\approx 40$-hour observation completed to date reveals the presence of a group of three Ly$\alpha$ nebulae associated with a bright quasar pair at $z\simeq3.23$ with projected separation of $\approx 500\rm~kpc$. Two of the nebulae are physically associated with the quasars which are likely powering the Ly$\alpha$ emission, and extend for $\gtrsim 100~\rm kpc$ at a surface brightness level of $\approx 6\times 10^{-19}~\rm erg~s^{-1}~cm^{-2}~arcsec^{-2}$. A third smaller ($\approx$35 kpc) nebula lies at a velocity offset of $\approx 1550$ km s$^{-1}$. Despite their clustered nature, the two large nebulae have properties similar to those observed in isolated quasars and exhibit no sharp decline in flux at the current depth, suggesting an even more extended distribution of gas around the quasars. We interpret the shape and the alignment of the two brighter nebulae as suggestive of the presence of an extended structure connecting the two quasar host galaxies, as seen for massive galaxies forming within gas-rich filaments in cosmological simulations.

The DESI Legacy Imaging Surveys are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image approximately 14,000 deg^2 of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamically adjusting the exposure times and pointing selection during observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering a catalog, constructed by using a probabilistic inference-based approach to estimate source shapes and brightnesses. The catalog includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4, 4.6, 12 and 22 micorons) observed by the Wide-field Infrared Survey Explorer (WISE) satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project.

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

Cosmological simulations predict that a significant fraction of the low-$z$ baryon budget resides in large-scale filaments in the form of a diffuse plasma at temperatures $T \sim 10^{5} - 10^{7}$ K. However, direct observation of this so-called warm-hot intergalactic medium (WHIM) has been elusive. In the $\Lambda$CDM paradigm, galaxy clusters correspond to the nodes of the cosmic web at the intersection of several large-scale filamentary threads. In previous work, we used HST/COS data to conduct the first survey of broad HI Ly$\alpha$ absorbers (BLAs) potentially produced by WHIM in inter-cluster filaments. We targeted a single QSO, namely Q1410, whose sight-line intersects $7$ independent inter-cluster axes at impact parameters $<3$ Mpc (co-moving), and found a tentative excess of a factor of ${\sim}4$ with respect to the field. Here, we further investigate the origin of these BLAs by performing a blind galaxy survey within the Q1410 field using VLT/MUSE. We identified $77$ sources and obtained the redshifts for $52$ of them. Out of the total sample of $7$ BLAs in inter-cluster axes, we found $3$ without any galaxy counterpart to stringent luminosity limits ($\sim 4 \times 10^{8}$ L$_{\odot}$ ${\sim} 0.01$ L$_{*}$), providing further evidence that these BLAs may represent genuine WHIM detections. We combined this sample with other suitable BLAs from the literature and inferred the corresponding baryon mean density for these filaments in the range $\Omega^{\rm fil}_{\rm bar}= 0.02-0.04$. Our rough estimates are consistent with the predictions from numerical simulations but still subject to large systematic uncertainties, mostly from the adopted geometry, ionization corrections and density profile.

We recently discovered that the luminous radio-quiet QSO LBQS0302-0019 at z=3.286 is likely accompanied by an obscured AGN at 20kpc projected distance, which we dubbed Jil. It represents the tightest candidate obscured/unobscured dual AGN system at z>3. To verify the dual AGN scenario we obtained deep $K_s$ band (rest-frame $V$ band) imaging with the VLT/HAWK-I+GRAAL instrument at 0.4" resolution during science verification in January 2018. Indeed, we detect the individual host galaxies of the QSO and Jil with estimated stellar masses of $\log(M_\star/M_{\odot})=11.4\pm0.5$ and $\log(M_\star/M_{\odot})=10.9\pm0.5$, respectively. Near-IR spectra obtained with VLT-KMOS reveal a clear [OIII] $\lambda$5007 line detection at the location of Jil which does not contribute significantly to the $K_s$ band flux. Both observations therefore corroborate the dual AGN scenario. A comparison to Illustris cosmological simulations suggests a parent halo mass of $\log(M_\mathrm{halo}/M_{\odot})=13.2\pm0.5$ for this interacting galaxy system, corresponding to a very massive dark matter halo at that epoch.

We present observations of escaping Lyman Continuum (LyC) radiation from 34 massive star-forming galaxies and 12 weak AGN with reliably measured spectroscopic redshifts at $z$$\simeq$2.3-4.1. We analyzed Hubble Space Telescope Wide Field Camera 3 (WFC3) mosaics of the Early Release Science field in three UVIS filters to sample the rest-frame LyC over this redshift range. With our best current assessment of the WFC3 systematics, we provide 1$\sigma$ upper limits for the average LyC emission of galaxies at $z$=2.35, 2.75, and 3.60 to $\sim$28.5, 28.1, and 30.7 mag in image stacks of 11-15 galaxies in the WFC3/UVIS F225W, F275W, and F336W, respectively. The LyC flux of weak AGN at $z$=2.62 and 3.32 are detected at 28.3 and 27.4 mag with SNRs of $\sim$2.7 and 2.5 in F275W and F336W for stacks of 7 and 3 AGN, respectively, while AGN at $z$=2.37 are constrained to $\gtrsim$27.9 mag at 1$\sigma$ in a stack of 2 AGN. The stacked AGN LyC light profiles are flatter than their corresponding non-ionizing UV continuum profiles out to radii of r$\lesssim$0."9, which may indicate a radial dependence of porosity in the ISM. With synthetic stellar SEDs fit to UV continuum measurements longwards of Ly$\alpha$ and IGM transmission models, we constrain the absolute LyC escape fractions to $f_{\rm esc}^{\rm abs}$$\simeq$$22^{+44}_{-22}$% at $z$=2.35 and $\lesssim$55% at $z$=2.75 and 3.60, respectively. All available data for galaxies, including published work, suggests a more sudden increase of $f_{\rm esc}$ with redshift at $z$$\simeq$2. Dust accumulating in (massive) galaxies over cosmic time correlates with increased HI column density, which may lead to reducing $f_{\rm esc}$ more suddenly at $z$$\lesssim$2. This may suggest that star-forming galaxies collectively contributed to maintaining cosmic reionization at redshifts $z$$\gtrsim$2-4, while AGN likely dominated reionization at $z$$\lesssim$2.

We present new MUSE observations of a galaxy group probed by a background quasar. The quasar sightline passes between multiple $z=0.28$ galaxies, whilst showing at the same redshift low ionised metal line species, including Ca II, Mg I, Mg II and Fe II. Based on the galaxy redshifts measured from the MUSE data, we estimate the galaxies to be part of a small galaxy group with a halo mass of $\approx6\times10^{12}$ M$_{\odot}$. We use the MUSE data to reveal the two dimensional dynamical properties of the gas and stars in the group galaxies, and relate these to the absorber kinematics. With these data we consider a number of scenarios for the nature of the gas probed by the sightline absorbers: a co-rotating gas halo associated with a single galaxy within the group; outflowing material from a single group member powered by recent star-formation; and cool dense gas associated with an intra-group medium. We find that the dynamics, galaxy impact parameters, star-formation rates, and the absorber strength suggest the cool gas can not be clearly associated with any single galaxy within the group. Instead we find that the observations are consistent with a superposition of cool gas clouds originating with the observed galaxies as they fall into the group potential, and are now likely in the process of forming the intra-group medium.

We present an analysis of the spatial distribution of gas and galaxies using new X-Shooter observations of $z\sim3-4$ quasars. Adding the X-Shooter data to an existing dataset of high resolution quasar spectroscopy, we use a total sample of 29 quasars alongside $\sim1700$ Lyman Break Galaxies in the redshift range $2<z<3.5$. Analysing the Ly$\alpha$ forest auto-correlation function using the full quasar sample, we find $s_0=0.081\pm0.006h^{-1}$Mpc. We then investigate the clustering and dynamics of Ly$\alpha$ forest absorbers around $z\sim3$ LBGs. From the redshift-space cross-correlation, we find $s_0=0.27\pm0.14h^{-1}$Mpc, with power-law slope $\gamma=1.1\pm0.2$. We make a first analysis of the dependence of this clustering length on absorber strength based on cuts in the sightline transmitted flux, finding a clear preference for stronger absorption features to be more strongly clustered around the galaxy population than weaker absorption features. Further, we calculate the projected correlation function, finding $r_0=0.24\pm0.04h^{-1}$Mpc (assuming a fixed slope $\gamma=1.1$). Taking this as the underlying real-space clustering, we fit the 2D cross-correlation function with a dynamical model incorporating the infall parameter and the peculiar velocity, finding $\beta_{\rm F}=1.02\pm0.22$ and $240\pm60$ km s$^{-1}$ respectively. This result shows a significant detection of gas infall relative to the galaxy population, whilst the measured velocity dispersion is consistent with the velocity uncertainties on the galaxy redshifts. We evaluate the Cauchy-Schwarz inequality between the galaxy-galaxy, absorber-absorber, and galaxy-absorber correlation functions, finding a result significantly less than unity: $\xi_{\rm ag}^2/(\xi_{\rm gg}\xi_{\rm aa})=0.25\pm0.14$, implying that galaxies and Ly$\alpha$ absorbers do not linearly trace the underlying dark matter distribution in the same way.

We present new results on the auto- and cross-correlation functions of galaxies and OVI absorbers in a $\sim 18~\textrm{Gpc}^3$ comoving volume at $z < 1$. We use a sample of 51,296 galaxies and 140 OVI absorbers in the column density range $13 \lesssim \log N \lesssim 15$ to measure two-point correlation functions in the two dimensions transverse and orthogonal to the line-of-sight $\xi(r_{\perp}, r_{\parallel})$. We furthermore infer the corresponding 'real-space' correlation functions, $\xi(r)$, by projecting $\xi(r_{\perp}, r_{\parallel})$ along $r_{\parallel}$, and assuming a power-law form, $\xi(r) = (r / r_0)^{-\gamma}$. Comparing the results from the absorber-galaxy cross-correlation function, $\xi_{\textrm{ag}}$, the galaxy auto-correlation function, $\xi_{\textrm{gg}}$, and the absorber auto-correlation function, $\xi_{\textrm{aa}}$, we constrain the statistical connection between galaxies and the metal-enriched intergalactic medium as a function of star-formation activity. We also compare these results to predictions from the EAGLE cosmological hydrodynamical simulation and find a reasonable agreement. We find that: (i) OVI absorbers show very little velocity dispersion with respect to galaxies on $\sim$ Mpc scales, likely $\lesssim$ 100 \kms; (ii) OVI absorbers and galaxies may not linearly trace the same underlying distribution of matter in general. In particular, our results demonstrate that OVI absorbers are less clustered, and potentially more extended around galaxies than galaxies are around themselves; (iii) On $\gtrsim 100$ kpc scales, the likelihood of finding OVI absorbers around star-forming galaxies is similar to the likelihood of finding OVI absorbers around non star-forming galaxies (abridged)

We present a new redshift survey, the 2dF Quasar Dark Energy Survey pilot (2QDESp), which consists of ${\approx}10000$ quasars from ${\approx}150$ deg$^2$ of the southern sky, based on VST-ATLAS imaging and 2dF/AAOmega spectroscopy. Combining our optical photometry with the WISE (W1,W2) bands we can select essentially contamination free quasar samples with $0.8{<}z{<}2.5$ and $g{<}20.5$. At fainter magnitudes, optical UVX selection is still required to reach our $g{\approx}22.5$ limit. Using both these techniques we observed quasar redshifts at sky densities up to $90$ deg$^{-2}$. By comparing 2QDESp with other surveys (SDSS, 2QZ and 2SLAQ) we find that quasar clustering is approximately luminosity independent, with results for all four surveys consistent with a correlation scale of $r_{0}{=}6.1{\pm}0.1 \: h^{-1}$Mpc, despite their decade range in luminosity. We find a significant redshift dependence of clustering, particularly when BOSS data with $r_{0}{=}7.3{\pm}0.1 \: h^{-1}$Mpc are included at $z{\approx}2.4$. All quasars remain consistent with having a single host halo mass of ${\approx}2{\pm}1{\times}10^{12} \: h^{-1}M_\odot$. This result implies that either quasars do not radiate at a fixed fraction of the Eddington luminosity or AGN black hole and dark matter halo masses are weakly correlated. No significant evidence is found to support fainter, X-ray selected quasars at low redshift having larger halo masses as predicted by the `hot halo' mode AGN model of Fanidakis et al. 2013. Finally, although the combined quasar sample reaches an effective volume as large as that of the original SDSS LRG sample, we do not detect the BAO feature in these data.

We present a survey of $z\sim3$ Ly$\alpha$ emitters (LAEs) within the fields of the VLT LBG Redshift Survey. The data encompasses 5 independent survey fields co-spatial with spectroscopic LBG data and covering a larger total area than previously analysed for LAE number counts and clustering. This affords an improved analysis over previous work by minimising the effects of cosmic variance and allowing the cross-clustering analysis of LAEs and LBGs. Our photometric sample consists of $\approx600$ LAE candidates, over an area of 1.07~deg$^2$, with equivalent widths of $\gtrsim65$~\AA\ and a flux limit of $\approx2\times10^{-17}$~erg~cm$^{-2}$~s$^{-1}$. From spectroscopic follow-up, we measured a success rate of $78\pm18\%$. We find the $R$-band continuum luminosity function to be $\sim10\times$ lower than the luminosity function of LBGs at this redshift, consistent with previous studies. Exploiting the large area of the survey, we estimate the LAE auto-correlation function and find a clustering length of $r_0=2.86\pm0.33~h^{-1}$~Mpc, low compared to the $z\sim3$ LBG population, but somewhat higher than previous LAE measurements. This corresponds to a median halo mass of $M_{\rm DM}=10^{11.0\pm0.3}~h^{-1}~$M$_{\odot}$. We present an analysis of clustering length versus continuum magnitude and find that the measurements for LAEs and LBGs are consistent at faint magnitudes. Our combined dataset of LAEs and LBGs allows us to measure, for the first time, the LBG-LAE cross-correlation, finding a clustering length of $r_0=3.29\pm0.57~h^{-1}$~Mpc and a LAE halo mass of $10^{11.1\pm0.4}~h^{-1}$~M$_{\odot}$. Overall, we conclude that LAEs inhabit primarily low mass halos, but form a relatively small proportion of the galaxy population found in such halos.

Using data from four deep fields (COSMOS, AEGIS, ECDFS, and CDFN), we study the correlation between the position of galaxies in the star formation rate (SFR) versus stellar mass plane and local environment at $z<1.1$. To accurately estimate the galaxy SFR, we use the deepest available Spitzer/MIPS 24 and Herschel/PACS datasets. We distinguish group environments ( $M_{halo}\sim$10$^{12.5-14.2}$$M_{\odot}$) based on the available deep X-ray data and lower halo mass environments based on the local galaxy density. We confirm that the Main Sequence (MS) of star forming galaxies is not a linear relation and there is a flattening towards higher stellar masses ( $M_*>10^{10.4-10.6}$ $M_{\odot}$), across all environments. At high redshift ( $0.5<z<1.1$ ), the MS varies little with environment. At low redshift ( $0.15<z<0.5$ ), group galaxies tend to deviate from the mean MS towards the region of quiescence with respect to isolated galaxies and less-dense environments. We find that the flattening of the MS toward low SFR is due to an increased fraction of bulge dominated galaxies at high masses. Instead, the deviation of group galaxies from the MS at low redshift is caused by a large fraction of red disk dominated galaxies which are not present in the lower density environments. Our results suggest that above a mass threshold ( $\sim10^{10.4}-10^{10.6}$$M_{\odot}$ ) stellar mass, morphology and environment act together in driving the evolution of the SF activity towards lower level. The presence of a dominating bulge and the associated quenching processes are already in place beyond $z\sim$1. The environmental effects appear, instead, at lower redshifts and have a long time-scale.

Modern analyses of structure formation predict a universe tangled in a 'cosmic web' of dark matter and diffuse baryons. These theories further predict that at low-z, a significant fraction of the baryons will be shock-heated to $T \sim 10^{5}-10^{7}$K yielding a warm-hot intergalactic medium (WHIM), but whose actual existence has eluded a firm observational confirmation. We present a novel experiment to detect the WHIM, by targeting the putative filaments connecting galaxy clusters. We use HST/COS to observe a remarkable QSO sightline that passes within $\Delta d = 3$ Mpc from the 7 inter-cluster axes connecting 7 independent cluster-pairs at redshifts $0.1 \le z \le 0.5$. We find tentative excesses of total HI, narrow HI (NLA; Doppler parameters $b<50$ km/s), broad HI (BLA; $b \ge 50$ km/s) and OVI absorption lines within rest-frame velocities of $\Delta v \lesssim 1000$ km/s from the cluster-pairs redshifts, corresponding to $\sim 2$, $\sim 1.7$, $\sim 6$ and $\sim 4$ times their field expectations, respectively. Although the excess of OVI likely comes from gas close to individual galaxies, we conclude that most of the excesses of NLAs and BLAs are truly intergalactic. We find that the covering fractions, $f_c$, of BLAs close to cluster-pairs are $\sim 4-7$ times higher than the random expectation (at the $\sim 2 \sigma$ c.l.), whereas the $f_c$ of NLAs and OVI are not significantly enhanced. We argue that a larger relative excess of BLAs compared to those of NLAs close to cluster-pairs may be a signature of the WHIM in inter-cluster filaments. By extending the present analysis to tens of sightlines our experiment offers a promising route to detect the WHIM.

We study the evolution of the total star formation (SF) activity, total stellar mass and halo occupation distribution in massive halos by using one of the largest X-ray selected sample of galaxy groups with secure spectroscopic identification in the major blank field surveys (ECDFS, CDFN, COSMOS, AEGIS). We provide an accurate measurement of SFR for the bulk of the star-forming galaxies using very deep mid-infrared Spitzer MIPS and far-infrared Herschel PACS observations. For undetected IR sources, we provide a well-calibrated SFR from SED fitting. We observe a clear evolution in the level of SF activity in galaxy groups. The total SF activity in the high redshift groups (0.5<z<1.1) is higher with respect to the low redshift (0.15<z<0.5) sample at any mass by 0.8+/-0.12 dex. A milder difference (0.35+/-0.1 dex) is observed between the low redshift bin and the groups at z~0. We show that the level of SF activity is declining more rapidly in the more massive halos than in the more common lower mass halos. We do not observe any evolution in the halo occupation distribution and total stellar mass- halo mass relations in groups. The picture emerging from our findings suggests that the galaxy population in the most massive systems is evolving faster than galaxies in lower mass halos, consistently with a "halo downsizing" scenario.

We use observations and simulation to study the relationship between star-forming galaxies and the intergalactic medium at z~3. The observed galaxy sample is based on spectroscopic redshift data from a combination of the VLT LBG Redshift Survey and Keck observations in fields centred on bright z>3 QSOs, whilst the simulation data is taken from GIMIC. In the simulation, we find that the dominant peculiar velocities are in the form of large-scale coherent motions of gas and galaxies. Gravitational infall of galaxies towards one another is also seen. At smaller scales, the peculiar velocities in the simulation over-predict the difference between the simulated real- and z-space galaxy correlation functions. Peculiar velocity pairs separated by <1Mpc/h have a smaller dispersion and explain the z-space correlation function better. The Ly{\alpha} auto- and cross-correlation functions in the GIMIC simulation show infall smaller than implied by previous work. This reduced infall may be due to the galaxy wide outflows implemented in the simulation. The main challenge in comparing these simulated results with the observed correlation functions comes from the presence of velocity errors for the observed LBGs which dominate at ~1Mpc/h scales. When these are taken into account, the observed LBG correlation function is well matched by a simulated $M_*>10^9M_\odot$ galaxy sample. The simulated cross-correlation shows similar neutral gas densities around galaxies as are seen in the observations. The simulated and observed Ly{\alpha} z-space autocorrelation functions agree well with each other. Our overall conclusion is that gas and galaxy peculiar velocities are towards the low end of expectation. Finally, little direct evidence is seen in either simulation or observations for high transmission near galaxies due to feedback. (Abridged)

We present an analysis of the clustering of galaxies from z ~ 2 to the present day using the WIRCam Deep Survey (WIRDS). WIRDS combines deep near-infrared data with the deep optical data from the CFHTLS Deep fields, providing a photometric data-set over an effective area of 2.4 sq. deg., from which accurate photometric redshifts and stellar masses can be estimated. We use the data to calculate the angular correlation function for galaxy samples split by star-formation activity, stellar mass and redshift. We estimate the real-space clustering for each sample, determining clustering lengths and power-law slopes. For galaxies selected by constant mass, we find that the clustering scale shows no evolution up to z ~ 2. Splitting the galaxy sample by mass, we see that higher mass galaxies have larger clustering scales at all redshifts. We use our results to test the GALFORM semi-analytical galaxy formation model and find the two are consistent. We split the galaxy population into passive and star-forming populations and find that the passive galaxy population shows a significantly larger clustering scale at all redshifts than the star-forming population below masses of ~$10^{11}M_\odot/h$, showing that even at z ~ 2 passive galaxies exist in denser environments than the bulk of the star-forming galaxy population. For star-forming galaxies with stellar masses $>10^{11}M_\odot/h$, we find a clustering strength of ~8Mpc/h across all redshifts, comparable to the measurements for the passive population. Also, for star-forming galaxies we see that clustering strength increases for higher stellar mass systems, however there is little sign of a mass dependence in passive galaxies. Finally, we investigate the connection between galaxy stellar mass and dark matter halo mass, showing a clear correlation between the two in both the WIRDS data and the GALFORM predictions.

We present a catalogue of 2135 galaxy redshifts from the VLT LBG Redshift Survey (VLRS), a spectroscopic survey of z ~ 3 galaxies in wide fields centred on background quasi-stellar objects. We have used deep optical imaging to select galaxies via the Lyman-break technique. Spectroscopy of the Lyman-break galaxies (LBGs) was then made using the Visible Multi-Object Spectrograph (VIMOS), giving a mean redshift of z=2.79. We analyse the clustering properties of the VLRS sample and also of the VLRS sample combined with the smaller area Keck-based survey of Steidel et al. From the semiprojected correlation function, wp({\sigma}) we find that the results are well fit with a single power-law model, with clustering scale lengths of r0=3.46+-0.41 and 3.83+-0.24 Mpc/h, respectively. We note that the corresponding combined {\xi}(r) slope is flatter than for local galaxies at {\gamma} = 1.5-1.6 rather than {\gamma}=1.8. This flat slope is confirmed by the z-space correlation function, {\xi}(s), and in the range 10<s<100 Mpc/h the VLRS shows ~2.5{\sigma} excess over the {\Lambda} cold dark matter. This excess may be consistent with recent evidence for non-Gaussianity in clustering results at z~1. We then analyse the LBG z-space distortions using the 2D correlation function, {\xi}({\sigma}, {\pi}), finding for the combined sample a large-scale infall parameter of $\beta$ = 0.38+-0.19 and a velocity dispersion of 420km/s. Based on our measured {\beta}, we are able to determine the gravitational growth rate, finding a value of f(z = 3)=0.99+-0.50 (or f{\sigma}8 = 0.26+-0.13), which is the highest redshift measurement of the growth rate via galaxy clustering and is consistent with {\Lambda}CDM. Finally, we constrain the mean halo mass for the LBG population, finding that the VLRS and combined sample suggest mean halo masses of log(MDM/Msun) = 11.57+-0.15 and 11.73+-0.07, respectively.

We present the results of a search for extended X-ray sources and their corresponding galaxy groups from 800-ks Chandra coverage of the All-wavelength Extended Groth Strip International Survey (AEGIS). This yields one of the largest X-ray selected galaxy group catalogs from a blind survey to date. The red-sequence technique and spectroscopic redshifts allow us to identify 100$%$ of reliable sources, leading to a catalog of 52 galaxy groups. The groups span the redshift range $z\sim0.066-1.544$ and virial mass range $M_{200}\sim1.34\times 10^{13}-1.33\times 10^{14}M_\odot$. For the 49 extended sources which lie within DEEP2 and DEEP3 Galaxy Redshift Survey coverage, we identify spectroscopic counterparts and determine velocity dispersions. We select member galaxies by applying different cuts along the line of sight or in projected spatial coordinates. A constant cut along the line of sight can cause a large scatter in scaling relations in low-mass or high-mass systems depending on the size of cut. A velocity dispersion based virial radius can more overestimate velocity dispersion in comparison to X-ray based virial radius for low mass systems. There is no significant difference between these two radial cuts for more massive systems. Independent of radial cut, overestimation of velocity dispersion can be created in case of existence of significant substructure and also compactness in X-ray emission which mostly occur in low mass systems. We also present a comparison between X-ray galaxy groups and optical galaxy groups detected using the Voronoi-Delaunay method (VDM) for DEEP2 data in this field.

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 present the results of a search for the most luminous star-forming galaxies at redshifts z~6 based on CFHT Legacy Survey data. We identify a sample of 40 Lyman break galaxies brighter than magnitude z'=25.3 across an area of almost 4 square degrees. Sensitive spectroscopic observations of seven galaxies provide redshifts for four, of which only two have moderate to strong Lyman alpha emission lines. All four have clear continuum breaks in their spectra. Approximately half of the Lyman break galaxies are spatially resolved in 0.7 arcsec seeing images, indicating larger sizes than lower luminosity galaxies discovered with the Hubble Space Telescope, possibly due to on-going mergers. The stacked optical and infrared photometry is consistent with a galaxy model with stellar mass ~ 10^{10} solar masses. There is strong evidence for substantial dust reddening with a best-fit A_V=0.7 and A_V>0.48 at 2 sigma confidence, in contrast to the typical dust-free galaxies of lower luminosity at this epoch. The spatial extent and spectral energy distribution suggest that the most luminous z~6 galaxies are undergoing merger-induced starbursts. The luminosity function of z=5.9 star-forming galaxies is derived. This agrees well with previous work and shows strong evidence for an exponential decline at the bright end, indicating that the feedback processes which govern the shape of the bright end are occurring effectively at this epoch.

[Abridged] Sub-mm observations of the William Herschel Deep Field using LABOCA revealed possible counterparts for 2 X-ray absorbed QSOs. The aim here is to exploit EVLA imaging at 8.4 GHz to establish the QSOs as radio/sub-mm sources. The challenge in reducing the EVLA data was the presence of a strong 4C source in the field. A new calibration algorithm was applied to the data to subtract it. The resulting thermal noise limited radio map covers the 16'x16' Extended WHDF. It contains 41 sources above a 4-sigma limit, 17 of which have primary beam corrected flux. The radio observations show that the absorbed AGN with LABOCA detections are coincident with radio sources, confirming the tendency for X-ray absorbed AGN to be sub-mm bright. These sources show strong ultraviolet excess (UVX) suggesting the nuclear sightline is gas- but not dust-absorbed. Of the 3 remaining LABOCA sources within the ~5' half-power beam width, 1 is identified with a faint nuclear X-ray/radio source in a nearby galaxy, 1 with a faint radio source and 1 is unidentified in any other band. More generally, differential radio source counts are in good agreement with previous observations, showing at S<50 micro-Jy a significant excess over a pure AGN model. In the full area, of 10 sources fainter than this limit, 6 have optical counterparts of which 3 are UVX (i.e. likely QSOs) including the 2 absorbed quasar LABOCA sources. The other faint radio counterparts are not UVX but are only slightly less blue and likely to be star-forming/merging galaxies, predominantly at lower luminosities and redshifts. The 4 faint, optically unidentified radio sources may be either dust obscured QSOs or galaxies. These high-z obscured AGN and lower-z star-forming populations are thus the main candidates to explain the observed excess in faint source counts and hence the excess radio background found previously by the ARCADE2 experiment.

We present a new near-infrared imaging survey in the four CFHTLS deep fields: the WIRCam Deep Survey (WIRDS). WIRDS comprises extremely deep, high quality (FWHM ~0.6") J, H and K imaging covering a total effective area of 2.1 deg^2 and reaching AB 50% completeness limits of ~24.5. We combine our images with the CFHTLS to create a unique eight-band ugrizJHK photometric catalogues in the CFHTLS deep fields; these four separate fields allow us to make a robust estimate of the effect of cosmic variance for all our measurements. We use these catalogues to estimate precise photometric redshifts, galaxy types and stellar masses for a unique sample of ~1.8 million galaxies. Our JHK number counts are consistent with previous studies. We apply the BzK selection to our gzK filter set and find that the star forming BzK selection successfully selects 76% of star-forming galaxies in the redshift range 1.4<z<2.5 in our photometric catalogue. The passive BzK selection returns 52% of the passive 1.4<z<2.5 population identified in the photometric catalogue. We present the galaxy stellar mass function as a function of redshift up to z=2 and present fits using double Schechter functions. A mass-dependent evolution of the mass function is seen with the numbers of galaxies with masses of log(M)<10.75 still evolving at z<1, but galaxies of higher mass reaching their present day numbers by z~0.8-1. This is consistent with the present picture of downsizing in galaxy evolution. We compare our results with the predictions of the GALFORM semi-analytical galaxy formation model and find that the simulations provide a relatively successful fit to the observed mass functions at intermediate masses (i.e. 10<log(M)<11). However, the GALFORM results under-predict the mass function at low masses, whilst the fit as a whole degrades beyond redshifts of z~1.2.

[Abridged] We investigate the global galaxy evolution over 12 Gyr (0.05<z<4.5), from the star formation rate density (SFRD), combining the VVDS Deep (17.5<=I<=24.0) and Ultra-Deep (23.00<=i<=24.75) surveys. We obtain a single homogeneous spectroscopic redshift sample, totalizing about 11000 galaxies. We estimate the rest-frame FUV luminosity function (LF) and luminosity density (LD), extract the dust attenuation of the FUV radiation using SED fitting, and derive the dust-corrected SFRD. We find a constant and flat faint-end slope alpha in the FUV LF at z<1.7. At z>1.7, we set alpha steepening with (1+z). The absolute magnitude M*_FUV brightens in the entire range 0<z<4.5, and at z>2 it is on average brighter than in the literature, while phi* is smaller. Our total LD shows a peak at z=2, present also when considering all sources of uncertainty. The SFRD history peaks as well at z=2. It rises by a factor of 6 during 2 Gyr (from z=4.5 to z=2), and then decreases by a factor of 12 during 10 Gyr down to z=0.05. This peak is mainly produced by a similar peak within the population of galaxies with -21.5<=M_FUV<=-19.5 mag. As times goes by, the total SFRD is dominated by fainter and fainter galaxies. The presence of a clear peak at z=2 and a fast rise at z>2 of the SFRD is compelling for models of galaxy formation. The mean dust attenuation A_FUV of the global galaxy population rises by 1 mag during 2 Gyr from z=4.5 to z=2, reaches its maximum at z=1 (A_FUV=2.2 mag), and then decreases by 1.1 mag during 7 Gyr down to z=0. The dust attenuation maximum is reached 2 Gyr after the SFRD peak, implying a contribution from the intermediate-mass stars to the dust production at z<2.

We investigate the contribution made by active galactic nuclei (AGN) to the high-redshift, luminous, submillimetre (submm) source population using deep (< 2 mJy/beam) Large Apex Bolometer Camera (LABOCA) 870 um observations within the William Herschel Deep Field (WHDF). This submm data complements previously obtained Chandra X-ray data of the field, from which AGN have been identified with the aid of follow-up optical spectra. From the LABOCA data, we detect 11 submm sources (based on a detection threshold of 3.2 sigma) with estimated fluxes of > 3 mJy/beam. Of the 11 identified submm sources, we find that 2 coincide with observed AGN and that, based on their hardness ratios, both of these AGN appear to be heavily obscured. We perform a stacking of the submm data around the AGN, which we group by estimated column density, and find that only the obscured (N_H > 10^22 cm^2) AGN show significant associated submm emission. These observations support the previous findings of Page et al and Hill et al that obscured AGN preferentially show submm emission. Hill et al have argued that, in this case, the contribution to the observed submm emission (and thus the submm background) from AGN heating of the dust in these sources may be higher than previously thought.

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 present the initial imaging and spectroscopic data acquired as part of the VLT VIMOS Lyman-break galaxy Survey. UBR (or UBVI) imaging covers five 36'x36' fields centred on bright z>3 QSOs, allowing ~21,000 2<z<3.5 galaxy candidates to be selected using the Lyman-break technique. We performed spectroscopic follow-up using VIMOS, measuring redshifts for 1020 z>2 LBGs and 10 z>2 QSOs from a total of 19 VIMOS pointings. From the galaxy spectra, we observe a 625+-510 km/s velocity offset between the ISM absorption and Ly-alpha emission line redshifts. Using the photometric and spectroscopic catalogues, we have analysed the galaxy clustering at z~3. In the photometric case, the angular correlation function, w(theta), is well fit by a double power-law with clustering scale-length, r_0 = 3.19+0.32-0.54 Mpc/h for r < 1 Mpc/h and r_0 = 4.59+0.31-0.33 Mpc/h at larger scales. Using the redshift sample we estimate the semi-projected correlation function, w_p(sigma) and find r_0 = 3.67+0.23-0.24 Mpc/h for the VLT sample and r_0 = 3.98+0.14-0.15 Mpc/h for a combined VLT+Keck sample. From the z-space correlation functions and assuming the above xi(r) models, we find that the combined VLT and Keck surveys require a galaxy velocity dispersion, <w_z^2>^1/2 ~ 700 km/s, higher than the ~400 km/s found by previous authors. We also measure a value for the gravitational growth rate parameter of beta(z=3) = 0.48+-0.17, implying a low value for the bias of b = 2.06+1.1-0.5. This value is consistent with the galaxy clustering amplitude which gives b = 2.22+-0.16, assuming the standard cosmology, implying that the evolution of the gravitational growth rate is also consistent with Einstein gravity. We have compared our LBG clustering amplitudes with lower redshift measurements and find that the clustering strength is not inconsistent with that of low-redshift L* spirals for simple 'long-lived' galaxy models.

We use CFHTLS deep optical data, WIRCam Deep Survey (WIRDS) NIR data and XMM data to identify z>1.1 clusters in the CFHTLS D1 and D4 fields. Counterparts to such clusters can not be identified without deep NIR data and as such the total of =1deg2 of J , H & Ks band imaging provided by WIRDS is an indispensable tool in such work. Using public XMM X-ray data, we identify extended X-ray sources in the two fields. The resulting catalogue of extended X-ray sources was analyzed for optical/NIR counterparts, using a red-sequence algorithm. Redshifts of candidate groups and clusters were estimated using the median photometric redshifts of detected counterparts and where available spectroscopic data. Additionally, we surveyed X-ray point sources for potential group systems at the limit of our detection range in the X-ray data. A catalogue of z > 1.1 cluster candidates in the two fields has been compiled and cluster masses, radii and temperatures have been estimated using the scaling relations. The catalogue consists of 15 z > 1.1 candidates. Three of the detections are previously published extended X-ray sources. Of note is JKSC 041 for which we identify possible structures at z = 0.8, z = 0.96, z = 1.13 and z = 1.49. We also make an independent detection of the massive cluster, XMMXCS J2215.9-1738. We use the z > 1.1 catalogue to compare the cluster number counts in these fields with models based on WMAP 7-year cosmology and find that the models slightly over-predict the observations, whilst at z>1.5 we do not detect any clusters. We note that cluster number counts at z > 1.1 are highly sensitive to the cosmological model, however a significant reduction in present statistical (due to available survey area) and systematic (due to cluster scaling relations) uncertainties is required in order to confidently constrain cosmological parameters using cluster number counts at high redshift.

Aims. The aim of this work is to study the contribution of the Ly-a emitters (LAE) to the star formation rate density (SFRD) of the Universe in the interval 2<z<6.6. Methods. We assembled a sample of 217 LAE from the Vimos-VLT Deep Survey (VVDS) with secure spectroscopic redshifts in the redshift range 2 < z < 6.62 and fluxes down to F=1.5x10^18 erg/s/cm^2. 133 LAE are serendipitous identifications in the 22 arcmin^2 total slit area surveyed with the VVDS-Deep and the 3.3 arcmin^2 from the VVDS Ultra-Deep survey, and 84 are targeted identifications in the 0.62 deg^2 surveyed with the VVDS-DEEP and 0.16 deg^2 from the Ultra-Deep survey. Among the serendipitous targets we estimate that 90% of the emission lines are most probably Ly-a, while the remaining 10% could be either [OII]3727 or Ly-a. We computed the LF and derived the SFRD from LAE at these redshifts. Results. The VVDS-LAE sample reaches faint line fluxes F(Lya) = 1.5x1^18 erg/s/cm^2 (corresponding to L(Lya)=10^41 erg/s at z~3) enabling the faint end slope of the luminosity function to be constrained to a=-1.6+-0.12 at redshift z~2.5 and to a=-1.78+0.1-0.12 at z=4, placing on firm statistical grounds trends found in previous LAE studies, and indicating that sub-L* LAE contribute significantly to the SFRD. The projected number density and volume density of faint LAE in 2<z<6.6 with F>1.5x10^18 erg/s/cm^2 are 33 galaxies/arcmin^2 and 4x10^-2 Mpc^-3, respectively. We find that the the observed luminosity function of LAE does not evolve from z=2 to z=6. This implies that, after correction for the redshift-dependent IGM absorption, the intrinsic LF must have evolved significantly over 3 Gyr. The SFRD from LAE contributes to about 20% of the SFRD at z =2-3, while the LAE appear to be the dominant source of star formation producing ionizing photons in the early universe z>5-6, becoming equivalent to that of Lyman Break galaxies.

We have measured redshifts for 243 z ~3 quasars in nine VLT VIMOS LBG redshift survey areas, each of which is centred on a known bright quasar. Using spectra of these quasars, we measure the cross-correlation between neutral hydrogen gas causing the Lya forest and 1020 Lyman-break galaxies at z ~3. We find an increase in neutral hydrogen absorption within 5 h^-1 Mpc of a galaxy in agreement with the results of Adelberger et al. (2003, 2005). The Lya-LBG cross-correlation can be described by a power-law on scales larger than 3 h^-1 Mpc. When galaxy velocity dispersions are taken into account our results at smaller scales (<2 h^-1 Mpc) are also in good agreement with the results of Adelberger et al. (2005). There is little immediate indication of a region with a transmission spike above the mean IGM value which might indicate the presence of star-formation feedback. To measure the galaxy velocity dispersions, which include both intrinsic LBG velocity dispersion and redshift errors, we have used the LBG-LBG redshift space distortion measurements of Bielby et al. (2010). We find that the redshift-space transmission spike implied in the results of Adelberger et al. (2003) is too narrow to be physical in the presence of the likely LBG velocity dispersion and is likely to be a statistical fluke. Nevertheless, neither our nor previous data can rule out the presence of a narrow, real-space transmission spike, given the evidence of the increased Lya absorption surrounding LBGs which can mask the spike's presence when convolved with a realistic LBG velocity dispersion. Finally, we identify 176 CIV systems in the quasar spectra and find an LBG-CIV correlation strength on scales of 10 h^-1 Mpc consistent with the relation measured at ~Mpc scales.

We investigate the use of simple colour cuts applied to the SDSS optical imaging to perform photometric selections of emission line galaxies out to z<1. From colour-cuts using the SDSS g, r and i bands, we obtain mean photometric redshifts of z=0.32+-0.08, z=0.44+-0.12 and z=0.65+-0.21. We further calibrate our high redshift selection using spectroscopic observations with the AAOmega spectrograph on the 4m Anglo-Australian Telescope (AAT), observing ~50-200 galaxy candidates in 4 separate fields. With just 1-hour of integration time and with seeing of ~1.6", we successfully determined redshifts for ~65% of the targeted candidates. We calculate the angular correlation functions of the samples and find correlation lengths of r0=2.64 h-1 Mpc, r0=3.62 h-1 Mpc and r0=5.88 h-1 Mpc for the low, mid and high redshift samples respectively. Comparing these results with predicted dark matter clustering, we estimate the bias parameter for each sample to be b=0.70, b=0.92 and b=1.46. We calculate the 2-point redshift-space correlation function at z~0.6 and find a clustering amplitude of s0=6.4 h-1 Mpc. Finally, we use our photometric sample to search for the Integrated Sachs-Wolfe signal in the WMAP 5yr data. We cross-correlate our three redshift samples with the WMAP W, V, Q and K bands and find an overall trend for a positive signal similar to that expected from models. However, the signal in each is relatively weak. Combining all three galaxy samples we find a signal of wTg(<100')=0.20+-0.12 microK in the WMAP W-band, a significance of 1.7sigma.

We first show that the new WMAP 3 year data confirm the detection by Myers et al 2004 of an extended SZ signal centred on 606 Abell (ACO) clusters with richness class, R>1. Our results also show SZ decrements around APM and 2MASS groups at increased significance than previously detected. We then follow the approach of Lieu et al 2006 and compare the stacked WMAP results for the decrement in 31 clusters with ROSAT X-ray profiles where Lieu et al found on average less SZ decrement in the WMAP 1 year data than predicted. We confirm that in the 3 year data these same clusters again show less SZ decrement than the X-ray data predicts. We then analysed the WMAP results for the 38 X-ray clusters with OVRO/BIMA measured SZ decrements as presented by Bonamente et al 2006. We again find that the average decrement is measured to be significantly less (5.5 sigma) than predicted by the Chandra X-ray data. Thus while we confirm the original detection of an extended SZ effect by Myers et al 2004, these X-ray comparisons may now suggest that the central SZ amplitudes detected by WMAP may actually be lower than expected. One possible explanation is that there is contamination of the WMAP SZ signal by radio sources in the clusters but we argue that this appears implausible. We then consider the possibility that the SZ decrement has been lensed away by foreground galaxy groups. Such a model predicts that the SZ decrement should depend on cluster redshift. A reduction in the SZ decrement with redshift is suggested from the ACO cluster sample and also from comparing the samples of Lieu et al 2006 and Bonamente et al 06. However, the mass power-spectrum would require a far higher amplitude than currently expected if lensing was to explain the SZ deficit in high redshift clusters.