In a multi-wavelength survey of $13$ quasars at $5.8\lesssim z\lesssim6.5$, that were pre-selected to be potentially young, we find five objects with extremely small proximity zone sizes that may imply UV-luminous quasar lifetimes of $\lesssim 100,000$ years. Proximity zones are regions of enhanced transmitted flux in the vicinity of the quasars that are sensitive to the quasars' lifetimes because the intergalactic gas has a finite response time to their radiation. We combine sub-mm observations from the Atacama Large Millimetre Array (ALMA) and the NOrthern Extended Millimeter Array (NOEMA), as well as deep optical and near-infrared spectra from medium-resolution spectrograph on the Very Large Telescope (VLT) and on the Keck telescopes, in order to identify and characterize these new young quasars, which provide valuable clues about the accretion behavior of supermassive black holes (SMBHs) in the early universe, and pose challenges on current black hole formation models to explain the rapid formation of billion solar mass black holes. We measure the quasars' systemic redshifts, black hole masses, Eddington ratios, emission line luminosities, and star formation rates of their host galaxies. Combined with previous results we estimate the fraction of young objects within the high-redshift quasar population at large to be $5\%\lesssim f_{\rm young}\lesssim 10\%$. One of the young objects, PSO J158-14, shows a very bright dust continuum flux ($F_{\rm cont}=3.46\pm 0.02\,\rm mJy$), indicating a highly star-bursting host galaxy with a star formation rate of approximately $1420\,M_{\odot}\,\rm yr^{-1}$.

We present deep (9 hours) Gemini-N/GNIRS near-infrared spectroscopic observations of ULAS J1342+0928, a luminous quasar at z=7.54. Various broad emission lines were detected, as well as the underlying continuum and iron forests over the rest-frame wavelength 970-2930A. There is a clear trend that higher-ionization emission lines show larger blueshifts with CIV1549 exhibiting 5510^{+240}_{-110} km s-1 blueshift with respect to the systematic redshift from the far-infrared [CII] 158um emission line. Those high ionization lines have wide profiles with FWHM more than 10000 km s-1. A modest blueshift of 340^{+110}_{-80} km s-1 is also seen in MgII, the lowest ionization line identified in the spectrum. The updated MgII-based black hole mass of M_BH=9.1_{-1.3}^{+1.4} x 10^8 M_sun and the Eddington ratio of L_bol/L_Edd=1.1_{-0.2}^{+0.2} confirm that ULAS J1342+0928 is powered by a massive and actively accreting black hole. There is no significant difference in the emission line ratios such as SiIV/CIV and AlIII/CIV when compared to lower-redshift quasars in a similar luminosity range, which suggests early metal pollution of the broad-line-region clouds. This trend also holds for the FeII/MgII line ratio, known as a cosmic clock that traces the iron enrichment in the early universe. Different iron templates and continuum fitting ranges were used to explore how the FeII/MgII measurement changes as a function of spectral modeling. Quasars at even higher redshift or at fainter luminosity range (L_bol<10^46 erg s-1) are needed to probe the sites of early metal enrichment and a corresponding change in the FeII/MgII ratio.

We report the discovery of a luminous quasar, J1007+2115 at $z=7.515$ ("P\={o}niu\={a}'ena"), from our wide-field reionization-era quasar survey. J1007+2115 is the second quasar now known at $z>7.5$, deep into the reionization epoch. The quasar is powered by a $(1.5\pm0.2)\times10^9$ $M_{\odot}$ supermassive black hole (SMBH), based on its broad MgII emission-line profile from Gemini and Keck near-IR spectroscopy. The SMBH in J1007+2115 is twice as massive as that in quasar J1342+0928 at $z=7.54$, the current quasar redshift record holder. The existence of such a massive SMBH just 700 million years after the Big Bang significantly challenges models of the earliest SMBH growth. Model assumptions of Eddington-limited accretion and a radiative efficiency of 0.1 require a seed black hole of $\gtrsim 10^{4}$ $M_{\odot}$ at $z=30$. This requirement suggests either a massive black hole seed as a result of direct collapse or earlier periods of rapid black hole growth with hyper-Eddington accretion and/or a low radiative efficiency. We measure the damping wing signature imprinted by neutral hydrogen absorption in the intergalactic medium (IGM) on J1007+2115's Ly$\alpha$ line profile, and find that it is weaker than that of J1342+0928 and two other $z\gtrsim7$ quasars. We estimate an IGM volume-averaged neutral fraction $\langle x\rm_{HI}\rangle=0.39^{+0.22}_{-0.13}$. This range of values suggests a patchy reionization history toward different IGM sightlines. We detect the 158 $\mu$m [C II] emission line in J1007+2115 with ALMA; this line centroid yields a systemic redshift of $z=7.5149\pm0.0004$ and indicates a star formation rate of $\sim210$ $M_{\odot}$ yr$^{-1}$ in its host galaxy.

Luminous $z\ge7$ quasars provide direct probes of the evolution of supermassive black holes (SMBHs) and the intergalactic medium (IGM) during the epoch of reionization (EoR). The Ly$\alpha$ damping wing absorption imprinted by neutral hydrogen in the IGM can be detected in a single EoR quasar spectrum, allowing the measurement of the IGM neutral fraction towards that line of sight. However, damping wing features have only been detected in two $z>7$ quasars in previous studies. In this paper, we present new high quality optical and near-infrared spectroscopy of the $z=7.00$ quasar DES J025216.64--050331.8 obtained with Keck/NIRES and Gemini/GMOS. By using the MgII single-epoch virial method, we find that it hosts a $\rm (1.39\pm0.16) \times10^{9} ~M_\odot$ SMBH accreting at an Eddington ratio of $\lambda_{\rm Edd}=0.7\pm0.1$, consistent with the values seen in other luminous $z\sim 7$ quasars. Furthermore, the Ly$\alpha$ region of the spectrum exhibits a strong damping wing absorption feature. The lack of associated metal absorption in the quasar spectrum indicates that this absorption is imprinted by a neutral IGM. Using a state-of-the-art model developed by Davies et al., we measure a volume-averaged neutral hydrogen fraction at $z=7$ of $\langle x_{\rm HI} \rangle = 0.70^{+0.20}_{-0.23} (^{+0.28}_{-0.48})$ within 68% (95%) confidence intervals when marginalizing over quasar lifetimes of $10^3\le t_{\rm Q}\le10^8$ yr. This is the highest IGM neutral fraction yet measured using reionization-era quasar spectra.

We identify a strong Ly-alpha damping wing profile in the spectrum of the quasar P183+05 at z=6.4386. Given the detection of several narrow metal absorption lines at z=6.40392, the most likely explanation for the absorption profile is that it is due to a damped Ly-alpha system. However, in order to match the data a contribution of an intergalactic medium 5-38% neutral or additional weaker absorbers near the quasar is also required. The absorption system presented here is the most distant damped Ly-alpha system currently known. We estimate an HI column density ($10^{20.68\pm0.25}\,$cm$^{-2}$), metallicity ([O/H]$=-2.92\pm 0.32$), and relative chemical abundances of a system consistent with a low-mass galaxy during the first Gyr of the universe. This object is among the most metal-poor damped Ly-alpha systems known and, even though it is observed only ~850 Myr after the big bang, its relative abundances do not show signatures of chemical enrichment by Population III stars.

The discovery of quasars few hundred megayears after the Big Bang represents a major challenge to our understanding of black holes and galaxy formation and evolution. Their luminosity is produced by extreme gas accretion onto black holes, which already reached masses of 10$^9$ M$_\odot$ by z ~ 6. Simultaneously, their host galaxies form hundreds of stars per year, using up gas in the process. To understand which environments are able to sustain the rapid formation of these extreme sources we started a VLT/MUSE effort aimed at characterizing the surroundings of a sample of 5.7 < z < 6.6 quasars dubbed: the Reionization Epoch QUasar InvEstigation with MUSE (REQUIEM) survey. We here present results of our searches for extended Ly-Alpha halos around the first 31 targets observed as part of this program. Reaching 5-sigma surface brightness limits of 0.1-1.1 x 10$^{-17}$ erg/s/cm$^2$/arcsec$^2$ over a 1 arcsec$^2$ aperture, we were able to unveil the presence of 12 Ly-Alpha nebulae, 8 of which are newly discovered. The detected nebulae show a variety of emission properties and morphologies with luminosities ranging from 8 x 10$^{42}$ to 2 x 10$^{44}$ erg/s, FWHMs between 300 and 1700 km/s, sizes < 30 pkpc, and redshifts consistent with those of the quasar host galaxies. As the first statistical and homogeneous investigation of the circum-galactic medium of massive galaxies at the end of the reionization epoch, the REQUIEM survey enables the study of the evolution of the cool gas surrounding quasars in the first 3 Gyr of the Universe. A comparison with the extended Ly-Alpha emission observed around bright (M$_{1450}$ < -25 mag) quasars at intermediate redshift indicates little variations on the properties of the cool gas from z ~ 6 to z ~ 3 followed by a decline in the average surface brightness down to z ~ 2.

We present ALMA observations of the [CII] 158 $\mu$m fine structure line and dust continuum emission from two quasars, SDSS J104433.04-012502.2 and SDSS J012958.51-003539.7, at z=5.78. The ALMA observations at 0.2'' resolution map the dust and gas on kpc scales. The spatially resolved emission shows a similar trend of decreasing [CII]-FIR ratios with increasing FIR surface brightnesses as was found in the infrared luminous galaxies with intense star formation. We confirm the velocity gradients of [CII] emission found previously in SDSS J0129-0035. No clear evidence of order motion is detected in SDSS J1044-0125. The velocity maps and PV-diagrams also suggest turbulent gas clumps in both objects. We tentatively detect a [CII] peak offset 4.9 kpc to the East of SDSS J1044-0125. This may be associated with an infalling companion, or node of gas outflow. All these results suggest significant dynamical evolution of the ISM in the nuclear region of these young quasar-starburst systems. We fit the velocity map of the [CII] emission from SDSS J0129-0035 with a rotating disk model. The result suggests a face-on system with an inclination angle of 16 +/- 20 degree and constrains the lower limit of the host galaxy dynamical mass to be 2.6x10^10 Msun within the [CII] emitting region. It is likely that SDSS J0129-0035, as well as other young quasars with super massive black hole masses on order of 10^7 to 10^8 Msun, falls close to the black hole and host galaxy mass relation defined by local galaxies.

We present ALMA 0.28"x0.20" (1.4kpc x 1.0kpc) resolution observations of the [CII] 158um line and underlying dust continuum in the host galaxy of the most distant quasar currently known, ULAS J1342+0928 at z=7.5413. Both the [CII] and continuum are detected and spatially resolved. The [CII] line is ~1.5 times more extended than the continuum emission, showing an elongated and complex structure of approximately 3.2kpc x 6.4kpc. Two separate peaks are clearly seen (>6sigma each) in three 100 km/s width [CII] channel maps. The velocity field of the [CII] gas does not show evidence of a coherent rotation field but rather chaotic motion, reminiscent of an on-going merger.

We report on ~0.35"(~2 kpc) resolution observations of the [CII] and dust continuum emission from five z>6 quasar host-companion galaxy pairs obtained with the Atacama Large Millimeter/submillimeter Array. The [CII] emission is resolved in all galaxies, with physical extents of 3.2-5.4 kpc. The dust continuum is on-average 40% more compact, which results in larger [CII] deficits in the center of the galaxies. However, the measured [CII] deficits are fully consistent with those found at lower redshifts. Four of the galaxies show [CII] velocity fields that are consistent with ordered rotation, while the remaining six galaxies show no clear velocity gradient. All galaxies have high (~80-200 km/s) velocity dispersions, consistent with the interpretation that the interstellar medium (ISM) of these high redshift galaxies is turbulent. By fitting the galaxies with kinematic models, we estimate the dynamical mass of these systems, which range between (0.3 -> 5.4) x 1E10 Msun. For the three closest separation galaxy pairs, we observe dust and [CII] emission from gas in between and surrounding the galaxies, which is an indication that tidal interactions are disturbing the gas in these systems. Although gas exchange in these tidal interactions could power luminous quasars, the existence of quasars in host galaxies without nearby companions suggests that tidal interactions are not the only viable method for fueling their active centers. These observations corroborate the assertion that accreting supermassive black holes do not substantially contribute to the [CII] and dust continuum emission of the quasar host galaxies, and showcase the diverse ISM properties of galaxies when the universe was less than one billion years old.

Massive, quiescent galaxies are already observed at redshift $z\sim4$, i.e. $\sim$1.5 Gyr after the Big Bang. Current models predict them to be formed via massive, gas-rich mergers at $z>6$. Recent ALMA observations of the cool gas and dust in $z\gtrsim$6 quasars have discovered [CII]- and far infrared-bright galaxies adjacent to several quasars. In this work, we present sensitive imaging and spectroscopic follow-up observations, with HST/WFC3, Spitzer/IRAC, VLT/MUSE, Magellan/FIRE and LBT/LUCI-MODS, of ALMA-detected, dust-rich companion galaxies of four quasars at $z\gtrsim 6$, specifically acquired to probe their stellar content and unobscured star formation rate. Three companion galaxies do not show significant emission in the observed optical/IR wavelength range. The photometric limits suggest that these galaxies are highly dust-enshrouded, with unobscured star formation rates SFR$_{\rm UV}<$few M$_{\odot}$ yr$^{-1}$, and a stellar content of $M_{*}<$10$^{10}$ M$_{\odot}$ yr$^{-1}$. However, the companion to PJ167-13 shows bright rest-frame UV emission (F140W AB = 25.48). Its SED resembles that of a star-forming galaxy with a total SFR$\sim$50 M$_{\odot}$ yr$^{-1}$ and $M_{*}\sim 9 \times 10^{9}$ M$_{\odot}$. All the companion sources are consistent with residing on the galaxy main sequence at $z\sim$6. Additional, deeper data from future facilities, such as JWST, are needed in order to characterize these gas-rich sources in the first Gyr of cosmic history.

We present IRAM/NOEMA, JCMT/SCUBA-2 and VLA observations of the most distant known gravitationally lensed quasar J0439+1634 at z = 6.5. We detect strong dust emission, [CII] 158 $\mu$m, [CI] 369 $\mu$m, [OI] 146 $\mu$m, CO(6-5), CO(7-6), CO(9-8), CO(10-9), H$_{\rm 2}$O $3_{\rm 1,2}-2_{\rm 2,1}$, and H$_{\rm 2}$O $3_{\rm 2,1}-3_{\rm 1,2}$ lines as well as a weak radio continuum. The strong [CII] line yields a systemic redshift of the host galaxy to be z=6.5188$\pm$0.0001. The magnification makes J0439+1634 the far-infrared (FIR) brightest quasar at z > 6 known, with the brightest [CII] line yet detected at this redshift. The FIR luminosity is (3.4$\pm$0.2)$\times$10$^{13}$ $\mu^{-1}$ $L_{\odot}$, where $\mu$ $\sim$ 2.6 - 6.6 is the magnification of the host galaxy, estimated based on the lensing configuration from HST imaging. We estimate the dust mass to be (2.2$\pm$0.1)$\times$10$^{9}$ $\mu^{-1}$ $M_{\odot}$. The CO Spectral Line Energy Distribution using four CO lines is best fit by a two-component model of the molecular gas excitation. The estimates of molecular gas mass derived from CO lines and atomic carbon mass are consistent, in the range of 3.9 - 8.9 $\times$10$^{10} \mu^{-1}$ $M_{\odot}$. The [CII]/[CI], [CII]/CO, and [OI]/[CII] line luminosity ratios suggest a photodissociation region model with more than one component. The ratio of H$_{\rm 2}$O $3_{\rm 2,1}-3_{\rm 1,2}$ line luminosity to $L_{\rm TIR}$ is consistent with values in local and high redshift ultra-/hyper-luminous infrared galaxies. The VLA observations reveal an unresolved radio continuum source, and indicate that J0439+1634 is a radio quiet quasar with R = 0.05 - 0.17.

We use ALMA observations of the host galaxy of the quasar ULAS-J1342+0928 at z=7.54 to study the dust continuum and far infrared lines emitted from its interstellar medium. The Rayleigh-Jeans tail of the dust continuum is well sampled with eight different spectral setups, and from a modified black body fit we obtain an emissivity coefficient of beta=1.85+-0.3. Assuming a standard dust temperature of 47 K we derive a dust mass of Mdust=0.35x10^8 Msol and a star formation rate of 150+-30 Msol/yr. We have >4sigma detections of the [CII]_158, [OIII]_88 and [NII]_122 atomic fine structure lines and limits on the [CI]_369, [OI]_146 and [NII]_205 emission. We also report multiple limits of CO rotational lines with Jup >= 7, as well as a tentative 3.3sigma detection of the stack of four CO lines (Jup=11, 10, 8, 7). We find line deficits that are in agreement with local ultra luminous infrared galaxies. Comparison of the [NII]_122 and [CII]_158 lines indicates that the [CII]_158 emission arises predominantly from the neutral medium, and we estimate that the photo-disassociation regions in J1342+0928 have densities <~5x10^4cm^-3. The data suggest that ~16% of hydrogen is in ionized form and that the HII regions have high electron densities of n_e>180 cm^-3. Our observations favor a low gas-to-dust ratio of <100 and a metallicity of the interstellar medium comparable to the Solar value. All the measurements presented here suggest that the host galaxy of J1342+0928 is highly enriched in metal and dust, despite being observed just 680 Myr after the Big Bang.

We present deep MUSE observations of five quasars within the first Gyr of the Universe ($z\gtrsim6$), four of which display extended Ly$\alpha$ halos. After PSF-subtraction, we reveal halos surrounding two quasars for the first time, as well as confirming the presence of two more halos for which tentative detections exist in long-slit spectroscopic observations and narrow-band imaging. The four Ly$\alpha$ halos presented here are diverse in morphology and size, they each display spatial asymmetry, and none are centred on the position of the quasar. Spectra of the diffuse halos demonstrate that none are dramatically offset in velocity from the systemic redshift of the quasars ($\Delta$ v $< 200$ kms$^{-1}$), however each halo shows a broad Ly$\alpha$ line, with a velocity width of order $\sim1000$ kms$^{-1}$. Total Ly$\alpha$ luminosities range between $\sim$ $2 \times 10^{43}$ erg s$^{-1}$ and $\sim$ $2 \times 10^{44}$ erg s$^{-1}$, reaching maximum radial extents of $13 - 30$ pkpc from the quasar positions. We find larger sizes and higher Ly$\alpha$ luminosities than previous literature results at this redshift, but find no correlation between the quasar properties and the Ly$\alpha$ halo, suggesting that the detected emission is most closely related to the physical properties of the circum-galactic medium

We present kpc-scale ALMA and HST imaging of the quasar PJ308-21 at $z$=$6.2342$, tracing dust, gas (via the [CII] 158 $\mu$m line) and young stars. At a resolution of $\sim0.3''$ ($\approx1.7$ kpc), the system is resolved over $>4''$ ($>$20 kpc). In particular, it features a main component, identified to be the quasar host galaxy, centered on the accreting supermassive black hole; and two other extended components on the West and East side, one redshifted and the other blueshifted relative to the quasar. The [CII] emission of the entire system stretches over $>$1500 km/s along the line of sight. All the components of the system are observed in dust, [CII], and rest-frame UV emission. The inferred [CII] luminosities [(0.9-4.6)$\times 10^9$ L$_\odot$], dust luminosities [(0.15-2.6)$\times10^{12}$ L$_\odot$], and rest-frame UV luminosities [(6.6-15)$\times10^{10}$ L$_\odot$], their ratios, and the implied gas/dust masses and star formation rates [11-290 M$_\odot$ yr$^{-1}$] are typical of high-redshift star-forming galaxies. A toy model of a single satellite galaxy that is tidally stripped by the interaction with the quasar host galaxy can account for the observed velocity and spatial extent of the two extended components. An outflow interpretation of the unique features in PJ308-21 is not supported by the data. PJ308-21 is thus one of the earliest galaxy mergers imaged at cosmic dawn.

We report high spatial resolution (~0.076", 410pc) Atacama Large Millimeter/submillimeter Array imaging of the dust continuum and the ionized carbon line [CII] in a luminous quasar host galaxy at z=6.6, 800 million years after the big bang. Based on previous studies, this galaxy hosts a ~1x10^9 M_sun black hole and has a star-formation rate of ~1500 M_sun/yr. The unprecedented high resolution of the observations reveals a complex morphology of gas within 3kpc of the accreting central black hole. The gas has a high velocity dispersion with little ordered motion along the line of sight, as would be expected from gas accretion that has yet to settle in a disk. In addition, we find the presence of [CII] cavities in the gas distribution (with diameters of ~0.5kpc), offset from the central black hole. This unique distribution and kinematics cannot be explained by a simple model. Plausible scenarios are that the gas is located in a truncated or warped disk, or the holes are created by interactions with nearby galaxies or due to energy injection into the gas. In the latter case, the energy required to form the cavities must originate from the central active galactic nucleus, as the required energy far exceeds the energy output expected from supernovae. This energy input into the gas, however, does not inhibit the high rate of star-formation. Both star-formation and black hole activity could have been triggered by interactions with satellite galaxies; our data reveal three additional companions detected in [CII] emission around the quasar.

We present sub-kpc-scale mapping of the 870 $\mu$m ALMA continuum emission in six luminous ($L_{\rm IR}~\sim~5~\times10^{12}$ L$_{\odot}$) submillimeter galaxies (SMGs) from the ALESS survey of the Extended Chandra Deep Field South. Our high-fidelity 0.07$''$-resolution imaging ($\sim$500 pc) reveals robust evidence for structures with deconvolved sizes of $\lesssim$0.5-1 kpc embedded within (dominant) exponential dust disks. The large-scale morphologies of the structures within some of the galaxies show clear curvature and/or clump-like structures bracketing elongated nuclear emission, suggestive of bars, star-forming rings, and spiral arms. In this interpretation, the ratio of the `ring' and `bar' radii (1.9$\pm$0.3) agrees with that measured for such features in local galaxies. These potential spiral/ring/bar structures would be consistent with the idea of tidal disturbances, with their detailed properties implying flat inner rotation curves and Toomre-unstable disks (Q<1). The inferred one-dimensional velocity dispersions ($\sigma_{\rm r}\lesssim$ 70-160 km s$^{-1}$) are marginally consistent with the limits implied if the sizes of the largest structures are comparable to the Jeans length. We create maps of the star formation rate density ($\Sigma_{\rm SFR}$) on $\sim$500 pc scales and show that the SMGs are able to sustain a given (galaxy-averaged) $\Sigma_{\rm SFR}$ over much larger physical scales than local (ultra-)luminous infrared galaxies. However, on 500 pc scales, they do not exceed the Eddington limit set by radiation pressure on dust. If confirmed by kinematics, the potential presence of non-axisymmetric structures would provide a means for net angular momentum loss and efficient star formation, helping to explain the very high star formation rates measured in SMGs.

The discovery of luminous quasars at redshifts up to 7.5 demonstrates the existence of several billion M_sun supermassive black holes (SMBHs) less than a billion years after the Big Bang. They are accompanied by intense star formation in their host galaxies, pinpointing sites of massive galaxy assembly in the early universe, while their absorption spectra reveal an increasing neutral intergalactic medium (IGM) at the epoch of reionization. Extrapolating from the rapid evolution of the quasar density at z=5-7, we expect that there is only one luminous quasar powered by a billion M_sun SMBH in the entire observable universe at z~9. In the next decade, new wide-field, deep near-infrared (NIR) sky surveys will push the redshift frontier to the first luminous quasars at z~9-10; the combination with new deep X-ray surveys will probe fainter quasar populations that trace earlier phases of SMBH growth. The identification of these record-breaking quasars, and the measurements of their BH masses and accretion properties require sensitive spectroscopic observations with next generation of ground-based and space telescopes at NIR wavelengths. High-resolution integral-field spectroscopy at NIR, and observations at millimeter and radio wavelengths, will together provide a panchromatic view of the quasar host galaxies and their galactic environment at cosmic dawn, connecting SMBH growth with the rise of the earliest massive galaxies. Systematic surveys and multiwavelength follow-up observations of the earliest luminous quasars will strongly constrain the seeding and growth of the first SMBHs in the universe, and provide the best lines of sight to study the history of reionization.

We report initial results from a large Gemini program to observe z>~5.7 quasars with GNIRS near-IR spectroscopy. Our sample includes 50 quasars with simultaneous ~0.85-2.5 micron spectra covering the rest-frame ultraviolet and major broad emission lines from Ly-alpha to MgII. We present spectral measurements for these quasars and compare to their lower-redshift counterparts at z=1.5-2.3. We find that when quasar luminosity is matched, there are no significant differences between the rest-UV spectra of z>~5.7 quasars and the low-z comparison sample. High-z quasars have similar continuum and emission line properties and occupy the same region in the black hole mass and luminosity space as the comparison sample, accreting at an average Eddington ratio of ~0.3. There is no evidence for super-Eddington accretion or hypermassive (>10^10 Msun) black holes within our sample. We find a mild excess of quasars with weak CIV lines relative to the control sample. Our results, corroborating earlier studies but with better statistics, demonstrate that these high-z quasars are already mature systems of accreting supermassive black holes operating with the same physical mechanisms as those at lower redshifts.

Luminous high-redshift quasars can be used to probe of the intergalactic medium (IGM) in the early universe because their UV light is absorbed by the neutral hydrogen along the line of sight. They help us to measure the neutral hydrogen fraction of the high-z universe, shedding light on the end of reionization epoch. In this paper, we present a discovery of a new quasar (PSO J006.1240+39.2219) at redshift $z=6.61\pm0.02$ from Panoramic Survey Telescope & Rapid Response System 1. Including this quasar, there are nine quasars above $z>6.5$ up to date. The estimated continuum brightness is $M_\text{1450}$=$-25.96\pm0.08$. PSO J006.1240+39.2219 has a strong Ly~$\alpha$ emission compared with typical low-redshift quasars, but the measured near-zone region size is $R_\text{NZ}=3.2\pm1.1$ proper megaparsecs, which is consistent with other quasars at z$\sim$6.

We present ALMA band 8 observations of the [OIII] 88um line and the underlying thermal infrared continuum emission in the z=6.08 quasar CFHQS J2100-1715 and its dust-obscured starburst companion galaxy (projected distance: ~60 kpc). Each galaxy hosts dust-obscured star formation at rates > 100 M_sun/yr, but only the quasar shows evidence for an accreting 10^9 M_sun black hole. Therefore we can compare the properties of the interstellar medium in distinct galactic environments in two physically associated objects, ~1 Gyr after the Big Bang. Bright [OIII] 88um emission from ionized gas is detected in both systems; the positions and line-widths are consistent with earlier [CII] measurements, indicating that both lines trace the same gravitational potential on galactic scales. The [OIII] 88um/FIR luminosity ratios in both sources fall in the upper range observed in local luminous infrared galaxies of similar dust temperature, although the ratio of the quasar is smaller than in the companion. This suggests that gas ionization by the quasar (expected to lead to strong optical [OIII] 5008A emission) does not dominantly determine the quasar's FIR [OIII] 88um luminosity. Both the inferred number of photons needed for the creation of O++ and the typical line ratios can be accounted for without invoking extreme (top-heavy) stellar initial mass functions in the starbursts of both sources.

Despite extensive efforts, only two quasars have been found at $z>7$ to date due to a combination of low spatial density and high contamination from more ubiquitous Galactic cool dwarfs in quasar selection. This limits our current knowledge of the super-massive black hole (SMBH) growth mechanism and reionization history. In this letter, we report the discovery of a luminous quasar at $z=7.021$, DELS J003836.10$-$152723.6 (hereafter J0038$-$1527), selected using photometric data from DESI Legacy imaging Survey (DELS), Pan-STARRS1 (PS1) imaging Survey, as well as Wide-field Infrared Survey Explore ($WISE$) mid-infrared all-sky survey. With an absolute magnitude of $M_{1450}$=$-$27.1 and bolometric luminosity of $L_{\rm Bol}$=5.6$\times$10$^{13}$ $L_\odot$, J0038$-$1527 is the most luminous quasar known at $z>7$. Deep optical to near infrared spectroscopic observations suggest that J0038-1527 hosts a 1.3 billion solar mass BH accreting at the Eddington limit, with an Eddington ratio of 1.25$\pm$0.19. The CIV broad emission line of J0038$-$1527 is blue-shifted by more than 3000 km s$^{-1}$ to the systemic redshift. More detailed investigations of the high quality spectra reveal three extremely high velocity CIV broad absorption lines (BALs) with velocity from 0.08 to 0.14 times the speed of light and total balnicity index of more than 5000 km s$^{-1}$, suggesting the presence of relativistic outflows. J0038$-$1527 is the first quasar found at the epoch of reionization (EoR) with such strong outflows and provides a unique laboratory to investigate AGN feedback on the formation and growth of the most massive galaxies in the early universe.

We present Atacama Large Millimeter Array 1mm observations of the rest-frame far-infrared (FIR) dust continuum in 27 quasars at redshifts 6.0 < z < 6.7. We detect FIR emission at >3sigma in all quasar host galaxies with flux densities at ~1900GHz in the rest-frame of 0.12 < S_rest,1900GHz < 5.9mJy, with a median (mean) flux density of 0.88mJy (1.59mJy). The implied FIR luminosities range from L_FIR = (0.27-13)x10^12 L_sun, with 74% of our quasar hosts having L_FIR > 10^12 L_sun. The estimated dust masses are M_dust = 10^7-10^9 M_sun. If the dust is heated only by star formation, then the star formation rates in the quasar host galaxies are between 50 and 2700 M_sun/yr. In the framework of the host galaxy-black hole coevolution model a correlation between ongoing black hole growth and star formation in the quasar host galaxy would be expected. However, combined with results from the literature to create a luminosity-limited quasar sample, we do not find a strong correlation between quasar UV luminosity (a proxy for ongoing black hole growth) and FIR luminosity (star formation in the host galaxy). The absence of such a correlation in our data does not necessarily rule out the coevolution model, and could be due to a variety of effects (including different timescales for black hole accretion and FIR emission).

Strong gravitational lensing provides a powerful probe of the physical properties of quasars and their host galaxies. A high fraction of the most luminous high-redshift quasars was predicted to be lensed due to magnification bias. However, no multiple imaged quasar was found at z>5 in previous surveys. We report the discovery of J043947.08+163415.7, a strongly lensed quasar at z=6.51, the first such object detected at the epoch of reionization, and the brightest quasar yet known at z>5. High-resolution HST imaging reveals a multiple imaged system with a maximum image separation theta ~ 0.2", best explained by a model of three quasar images lensed by a low luminosity galaxy at z~0.7, with a magnification factor of ~50. The existence of this source suggests that a significant population of strongly lensed, high redshift quasars could have been missed by previous surveys, as standard color selection techniques would fail when the quasar color is contaminated by the lensing galaxy.

This is the third paper in a series aims at finding reionzation-era quasars with the combination of DESI Legacy imaging Surveys (DELS) and near-infrared imaging surveys, such as the UKIRT Hemisphere Survey (UHS), as well as the Wide-field Infrared Survey Explore ($WISE$) mid-infrared survey. In this paper, we describe the updated quasar candidate selection procedure, report the discovery of 16 quasars at $6.4\lesssim z \lesssim6.9$ from area of $\sim$13,020 deg$^2$, and present the quasar luminosity function (QLF) at $z\sim6.7$. The measured QLF follows $\Phi(L_{1450})\propto L_{1450}^{-2.35}$ in the magnitude range $27.6<M_{1450}<-25.5$. We determine the quasar comoving spatial density at $\langle z \rangle$=6.7 and $M_{1450}<-26.0$ to be $\rm 0.39\pm0.11 Gpc^{-3}$ and find that the exponential density evolution parameter to be $k=-0.78\pm0.18$ from $z\sim6$ to $z\sim6.7$, corresponding to a rapid decline by a factor of $\sim 6$ per unit redshift towards earlier epoch, a rate significantly faster than that at $z\sim 3- 5$. The cosmic time between $z\sim6$ and $z\sim6.7$ is only 121 Myrs. The quasar density declined by a factor of more than three within such short time requires that SMBHs must grow rapidly or they are less radiatively efficient at higher redshifts. We measured quasar comoving emissivity at $z\sim6.7$ which indicate that high redshift quasars are highly unlikely to make a significant contribution to hydrogen reionization. The broad absorption line (BAL) quasar fraction at $z\gtrsim6.5$ is measured to be $\gtrsim$22%. In addition, we also report the discovery of additional five quasars at $z\sim6$ in the appendix.

Bright high-redshift quasars ($z>6$), hosting supermassive black holes ($\rm M_{BH} > 10^8 M_{\odot}$), are expected to reside in massive host galaxies embedded within some of the earliest and most massive galaxy overdensities. We analyze 1.2\,mm ALMA dust continuum maps of 35 bright quasars at $6<z<7$ and search the primary beam for excess dust continuum emission from sources with L$_{\rm IR} \gtrsim{ 10^{12}}$\,L$_{\odot}$ as evidence for early protoclusters. We compare the detection rates of continuum sources at $\geq5\sigma$ significance in the fields surrounding the quasars (A$_{\rm eff}$ = 4.3\,arcmin$^2$) with millimeter number counts in blank field surveys. We discover 15 millimeter sources in the fields excluding the quasars themselves, corresponding to an overdensity $\delta_{\rm gal} \equiv (N_{\rm gal} - N_{\rm exp})/N_{\rm exp} = -0.07\pm0.56$, consistent with no detected overdensity of dusty galaxies within 140 physical kpc of the quasars. However, the apparent lack of continuum overdensity does not negate the hypothesis that quasars live in overdense environments, as evidenced by strong [CII] overdensities found on the same scales to similarly-selected quasars. The small field of view of ALMA could miss a true overdensity if it exists on scales larger than 1\,cMpc, if the quasar is not centered in the overdensity, or if quasar feedback plays a role close to the quasar, but it is most likely that the large line of sight volume probed by a continuum survey will wash out a true overdensity signal. We discuss the necessary factors in determining the bias with which dusty star-forming galaxies trace true dark matter overdensities in order to improve upon overdensity searches in dust continuum.

We present further results of a survey for absorption line systems in the spectra of four high redshift quasars (5.79 $\le$ z$_{\textrm{em}}$ $\le$ 6.13) obtained with the ESO Very Large Telescope X-Shooter. We identify 36 $\textrm{CIV}$ and 7 $\textrm{SiIV}$ systems with a $\ge$ 5$\sigma$ significance. The highest redshift $\textrm{CIV}$ and $\textrm{SiIV}$ absorbers identified in this work are at z = 5.80738 $\pm$ 0.00017 and z = 5.77495 $\pm$ 0.00038, respectively. We compute the comoving mass density of $\textrm{SiIV}$ ($\Omega_{\textrm{SiIV}}$) and find that it evolves from $\Omega_{\textrm{SiIV}}$ = 4.3$^{+2.1}_{-2.1}$ $\times$10$^{-9}$ at <z> = 5.05 to $\Omega_{\textrm{SiIV}}$ = 1.4$^{+0.6}_{-0.4}$ $\times$10$^{-9}$ at <z> = 5.66. We also measure $\Omega_{\textrm{CIV}}$ = 1.6$^{+0.4}_{-0.1}$ $\times$10$^{-8}$ at <z> = 4.77 and $\Omega_{\textrm{CIV}}$ = 3.4$^{+1.6}_{-1.1}$ $\times$10$^{-9}$ at <z> = 5.66. We classify our $\textrm{CIV}$ absorber population by the presence of associated $\textit{low}$ and/or $\textit{high ionisation}$ systems and compute their velocity width ($\Delta$v$_{90}$). We find that all $\textrm{CIV}$ systems with $\Delta$v$_{90}$ > 200 kms$^{-1}$ have associated $\textit{low ionisation}$ systems. We investigate two such systems, separated by 550 physical kpc along a line of sight, and find it likely that they are both tracing a multi-phase medium where hot and cold gas is mixing at the interface between the CGM and IGM. We further discuss the \textrm{MgII} systems presented in a previous work and we identify 5 $\textrm{SiII}$, 10 $\textrm{AlII}$, 12 $\textrm{FeII}$, 1 $\textrm{CII}$, 7 $\textrm{MgI}$ and 1 $\textrm{CaII}$ associated transitions. We compute the respective comoving mass densities in the redshift range 2 to 6, as allowed by the wavelength coverage.

We present the discovery of the radio-loud quasar PSO J352.4034-15.3373 at z=5.84 pm 0.02. This quasar is the radio brightest source known, by an order of magnitude, at z~6 with a flux density in the range of 8-100 mJy from 3GHz to 230MHz and a radio loudness parameter R>~1000. This source provides an unprecedented opportunity to study powerful jets and radio-mode feedback at the highest redshifts, and presents the first real chance to probe deep into the neutral intergalactic medium by detecting 21 cm absorption at the end of cosmic reionization.

We present high angular resolution imaging ($23.9 \times 11.3$ mas, $138.6 \times 65.5$ pc) of the radio-loud quasar PSO~J352.4034$-$15.3373 at $z=5.84$ with the Very Long Baseline Array (VLBA) at 1.54 GHz. This quasar has the highest radio-to-optical flux density ratio at such a redshift, making it the radio-loudest source known to date at $z \sim 6$. The VLBA observations presented here resolve this quasar into multiple components with an overall linear extent of 1.62 kpc ($0\rlap{.}{''}28$) and with a total flux density of $6.57 \pm 0.38$ mJy, which is about half of the emission measured at a much lower angular resolution. The morphology of the source is comparable with either a radio core with a one-sided jet, or a compact or a medium-size Symmetric Object (CSO/MSO). If the source is a CSO/MSO, and assuming an advance speed of $0.2c$, then the estimated kinematic age is $\sim 10^4$ yr.

Quasars are the most luminous non-transient objects known and as a result they enable studies of the Universe at the earliest cosmic epochs. Despite extensive efforts, however, the quasar ULAS J1120+0641 at z=7.09 has remained the only one known at z>7 for more than half a decade. Here we report observations of the quasar ULAS J134208.10+092838.61 (hereafter J1342+0928) at redshift z=7.54. This quasar has a bolometric luminosity of 4e13 times the luminosity of the Sun and a black hole mass of 8e8 solar masses. The existence of this supermassive black hole when the Universe was only 690 million years old---just five percent of its current age---reinforces models of early black-hole growth that allow black holes with initial masses of more than about 1e4 solar masses or episodic hyper-Eddington accretion. We see strong evidence of absorption of the spectrum of the quasar redwards of the Lyman alpha emission line (the Gunn-Peterson damping wing), as would be expected if a significant amount (more than 10 per cent) of the hydrogen in the intergalactic medium surrounding J1342+0928 is neutral. We derive a significant fraction of neutral hydrogen, although the exact fraction depends on the modelling. However, even in our most conservative analysis we find a fraction of more than 0.33 (0.11) at 68 per cent (95 per cent) probability, indicating that we are probing well within the reionization epoch of the Universe.

We use ALMA observations of four sub-millimetre galaxies (SMGs) at $z\sim2-3$ to investigate the spatially resolved properties of the inter-stellar medium (ISM) at scales of 1--5 kpc (0.1--0.6$''$). The velocity fields of our sources, traced by the $^{12}$CO($J$=3-2) emission, are consistent with disk rotation to first order, implying average dynamical masses of $\sim$3$\times10^{11}$M$_{\odot}$ within two half-light radii. Through a Bayesian approach we investigate the uncertainties inherent to dynamically constraining total gas masses. We explore the covariance between the stellar mass-to-light ratio and CO-to-H$_{2}$ conversion factor, $\alpha_{\rm CO}$, finding values of $\alpha_{\rm CO}=1.1^{+0.8}_{-0.7}$ for dark matter fractions of 15 \%. We show that the resolved spatial distribution of the gas and dust continuum can be uncorrelated to the stellar emission, challenging energy balance assumptions in global SED fitting. Through a stacking analysis of the resolved radial profiles of the CO(3-2), stellar and dust continuum emission in SMG samples, we find that the cool molecular gas emission in these sources (radii $\sim$5--14 kpc) is clearly more extended than the rest-frame $\sim$250 $\mu$m dust continuum by a factor $>2$. We propose that assuming a constant dust-to-gas ratio, this apparent difference in sizes can be explained by temperature and optical-depth gradients alone. Our results suggest that caution must be exercised when extrapolating morphological properties of dust continuum observations to conclusions about the molecular gas phase of the ISM.

We present a 45 ks Chandra observation of the quasar ULAS J1342+0928 at z=7.54. We detect 14.0^{+4.8}_{-3.7} counts from the quasar in the observed-frame energy range 0.5-7.0 keV (6-sigma detection), representing the most distant non-transient astronomical source identified in X-rays to date. The present data are sufficient only to infer rough constraints on the spectral parameters. We find an X-ray hardness ratio of HR = -0.51^{+0.26}_{-0.28} between the 0.5-2.0 keV and 2.0-7.0 keV ranges and derive a power-law photon index of Gamma = 1.95^{+0.55}_{-0.53}. Assuming a typical value for high-redshift quasars of Gamma = 1.9, ULAS J1342+0928 has a 2-10 keV rest-frame X-ray luminosity of L_{2-10} = 11.6^{+4.3}_{-3.5} x 10^{44} erg/s. Its X-ray-to-optical power-law slope is alpha_{OX}=-1.67^{+0.16}_{-0.10}, consistent with the general trend indicating that the X-ray emission in the most bolometrically powerful quasars is weaker relative to their optical emission.

Quasars (QSOs) hosting supermassive black holes are believed to reside in massive halos harboring galaxy overdensities. However, many observations revealed average or low galaxy densities around $z\gtrsim6$ QSOs. This could be partly because they measured galaxy densities in only tens of arcmin$^2$ around QSOs and might have overlooked potential larger scale galaxy overdensities. Some previous studies also observed only Lyman break galaxies (LBGs, massive older galaxies) and missed low mass young galaxies like Ly$\alpha$ emitters (LAEs) around QSOs. Here we present observations of LAE and LBG candidates in $\sim700$ arcmin$^2$ around a $z=6.61$ luminous QSO using Subaru Telescope Suprime-Cam with narrow/broadbands. We compare their sky distributions, number densities and angular correlation functions with those of LAEs/LBGs detected in the same manner and comparable data quality in our control blank field. In the QSO field, LAEs and LBGs are clustering in 4-20 comoving Mpc angular scales, but LAEs show mostly underdensity over the field while LBGs are forming $30\times60$ comoving Mpc$^2$ large scale structure containing 3-$7\sigma$ high density clumps. The highest density clump includes a bright (23.78 mag in the narrowband) extended ($\gtrsim 16$ kpc) Ly$\alpha$ blob candidate, indicative of a dense environment. The QSO could be part of the structure but is not located exactly at any of the high density peaks. Near the QSO, LAEs show underdensity while LBGs average to $4\sigma$ excess densities compared to the control field. If these environments reflect halo mass, the QSO may not be in the most massive halo, but still in a moderately massive one.

During reionization, neutral hydrogen in the intergalactic medium (IGM) imprints a damping wing absorption feature on the spectrum of high-redshift quasars. A detection of this signature provides compelling evidence for a significantly neutral Universe, and enables measurements of the hydrogen neutral fraction $x_{\rm HI}(z)$ at that epoch. Obtaining reliable quantitative constraints from this technique, however, is challenging due to stochasticity induced by the patchy inside-out topology of reionization, degeneracies with quasar lifetime, and the unknown unabsorbed quasar spectrum close to rest-frame Ly$\alpha$. We combine a large-volume semi-numerical simulation of reionization topology with 1D radiative transfer through high-resolution hydrodynamical simulations of the high-redshift Universe to construct models of quasar transmission spectra during reionization. Our state-of-the-art approach captures the distribution of damping wing strengths in biased quasar halos that should have reionized earlier, as well as the erosion of neutral gas in the quasar environment caused by its own ionizing radiation. Combining this detailed model with our new technique for predicting the quasar continuum and its associated uncertainty, we introduce a Bayesian statistical method to jointly constrain the neutral fraction of the Universe and the quasar lifetime from individual quasar spectra. We apply this methodology to the spectra of the two highest redshift quasars known, ULAS J1120+0641 and ULAS J1342+0928, and measured volume-averaged neutral fractions $\langle x_{\rm HI} \rangle(z=7.09)=0.48^{+0.26}_{-0.26}$ and $\langle x_{\rm HI} \rangle(z=7.54)=0.60^{+0.20}_{-0.23}$ (posterior medians and 68% credible intervals) when marginalized over quasar lifetimes of $10^3 \leq t_{\rm q} \leq 10^8$ years.

Measuring the proximity effect and the damping wing of intergalactic neutral hydrogen in quasar spectra during the epoch of reionization requires an estimate of the intrinsic continuum at rest-frame wavelengths $\lambda_{\rm rest}\sim1200$-$1260$ {\AA}. In contrast to previous works which used composite spectra with matched spectral properties or explored correlations between parameters of broad emission lines, we opted for a non-parametric predictive approach based on principal component analysis (PCA) to predict the intrinsic spectrum from the spectral properties at redder (i.e. unabsorbed) wavelengths. We decomposed a sample of $12764$ spectra of $z\sim2$-$2.5$ quasars from SDSS/BOSS into 10 red-side ($1280$ {\AA} $<\lambda_{\rm rest}<2900$ {\AA}) and 6 blue-side ($1180$ {\AA} $<\lambda_{\rm rest}<1280$ {\AA}) PCA basis spectra, and constructed a projection matrix to predict the blue-side coefficients from a fit to the red-side spectrum. We found that our method predicts the blue-side continuum with $\sim6$-$12\%$ precision and $\lesssim1\%$ bias by testing on the full training set sample. We then computed predictions for the blue-side continua of the two quasars currently known at $z>7$: ULAS J1120+0641 ($z=7.09$) and ULAS J1342+0928 ($z=7.54$). Both of these quasars are known to exhibit extreme emission line properties, so we individually calibrated the precision of the continuum predictions from similar quasars in the training set. We find that both $z>7$ quasars, and in particular ULAS J1342+0928, show signs of damping wing-like absorption at wavelengths redward of Ly$\alpha$.

We present a survey of the [CII] 158 $\mu$m line and underlying far-infrared (FIR) dust continuum emission in a sample of 27 z>6 quasars using the Atacama Large Millimeter Array (ALMA) at ~1" resolution. The [CII] line was significantly detected (at >5-sigma) in 23 sources (85%). We find typical line luminosities of $L_{\rm [CII]}=10^{9-10}$ L$_\odot$, and an average line width of ~385 km/s. The [CII]-to-far-infrared luminosity ratio ([CII]/FIR) in our sources span one order of magnitude, highlighting a variety of conditions in the star-forming medium. Four quasar host galaxies are clearly resolved in their [CII] emission on a few kpc scales. Basic estimates of the dynamical masses of the host galaxies give masses between $2\times10^{10}$ and $2\times10^{11}$ M$_\odot$, i.e., more than an order of magnitude below what is expected from local scaling relations, given the available limits on the masses of the central black holes ($>3\times10^8$ M$_\odot$, assuming Eddington-limited accretion). In stacked ALMA [CII] spectra of individual sources in our sample, we find no evidence of a deviation from a single Gaussian profile. The quasar luminosity does not strongly correlate with either the [CII] luminosity or equivalent width. This survey (with typical on-source integration times of 8 min) showcases the unparalleled sensitivity of ALMA at millimeter wavelengths, and offers a unique reference sample for the study of the first massive galaxies in the universe.

We present IRAM/NOEMA and JVLA observations of the quasar J1342+0928 at z=7.54 and report detections of copious amounts of dust and [CII] emission in the interstellar medium (ISM) of its host galaxy. At this redshift, the age of the universe is 690 Myr, about 10% younger than the redshift of the previous quasar record holder. Yet, the ISM of this new quasar host galaxy is significantly enriched by metals, as evidenced by the detection of the [CII] 158micron cooling line and the underlying far-infrared (FIR) dust continuum emission. To the first order, the FIR properties of this quasar host are similar to those found at a slightly lower redshift (z~6), making this source by far the FIR-brightest galaxy known at z>7.5. The [CII] emission is spatially unresolved, with an upper limit on the diameter of 7 kpc. Together with the measured FWHM of the [CII] line, this yields a dynamical mass of the host of <1.5x10^11 M_sun. Using standard assumptions about the dust temperature and emissivity, the NOEMA measurements give a dust mass of (0.6-4.3)x10^8 M_sun. The brightness of the [CII] luminosity, together with the high dust mass, imply active ongoing star formation in the quasar host. Using [CII]-SFR scaling relations, we derive star formation rates of 85-545 M_sun/yr in the host, consistent with the values derived from the dust continuum. Indeed, an episode of such past high star formation is needed to explain the presence of ~10^8 M_sun of dust implied by the observations.

We present final statistics from a survey for intervening MgII absorption towards 100 quasars with emission redshifts between $z=3.55$ and $z=7.08$. Using infrared spectra from Magellan/FIRE, we detect 279 cosmological MgII absorbers, and confirm that the incidence rate of $W_r>0.3 \AA$ MgII absorption per comoving path length does not evolve measurably between $z=0.25$ and $z=7$. This is consistent with our detection of seven new MgII systems at $z>6$, a redshift range that was not covered in prior searches. Restricting to relatively strong MgII systems ($W_r>1$\AA), there is significant evidence for redshift evolution. These systems roughly double in number density between $z=0$ and $z=2$-$3$, but decline by an order of magnitude from this peak by $z\sim 6$. This evolution mirrors that of the global star formation rate density, which could reflect a connection between star formation feedback and strong MgII absorbers. We compared our results to the Illustris cosmological simulation at $z=2$-$4$ by assigning absorption to catalogued dark-matter halos and by direct extraction of spectra from the simulation volume. To reproduce our results using the halo catalogs, we require circumgalactic (CGM) MgII envelopes within halos of progressively smaller mass at earlier times. This occurs naturally if we define the lower integration cutoff using SFR rather than mass. MgII profiles calculated directly from the Illustris volume yield far too few strong absorbers. We argue that this arises from unresolved phase space structure of CGM gas, particularly from turbulent velocities on sub-mesh scales. The presence of CGM MgII at $z>6$-- just $\sim 250$ Myr after the reionization redshift implied by Planck--suggests that enrichment of intra-halo gas may have begun before the presumed host galaxies' stellar populations were mature and dynamically relaxed. [abridged]

Quasars are galaxies hosting accreting supermassive black holes; due to their brightness, they are unique probes of the early universe. To date, only few quasars have been reported at $z > 6.5$ ($<$800 Myr after the Big Bang). In this work, we present six additional $z \gtrsim 6.5$ quasars discovered using the Pan-STARRS1 survey. We use a sample of 15 $z \gtrsim 6.5$ quasars to perform a homogeneous and comprehensive analysis of this highest-redshift quasar population. We report four main results: (1) the majority of $z\gtrsim$6.5 quasars show large blueshifts of the broad CIV 1549\AA$\,$emission line compared to the systemic redshift of the quasars, with a median value $\sim$3$\times$ higher than a quasar sample at $z\sim$1; (2) we estimate the quasars' black hole masses (M$\rm_{BH}\sim$0.3$-$5 $\times$ 10$^{9}$ M$_{\odot}$) via modeling of the MgII 2798\AA$\,$emission line and rest-frame UV continuum; we find that quasars at high redshift accrete their material (with $\langle (L_{\mathrm{bol}}/L_{\mathrm{Edd}}) \rangle = 0.39$) at a rate comparable to a luminosity-matched sample at lower$-$redshift, albeit with significant scatter ($0.4$ dex); (3) we recover no evolution of the FeII/MgII abundance ratio with cosmic time; (4) we derive near zone sizes; together with measurements for $z\sim6$ quasars from recent work, we confirm a shallow evolution of the decreasing quasar near zone sizes with redshift. Finally, we present new millimeter observations of the [CII] 158 $\mu$m emission line and underlying dust continuum from NOEMA for four quasars, and provide new accurate redshifts and [CII]/infrared luminosities estimates. The analysis presented here shows the large range of properties of the most distant quasars.

We utilize the Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT) to search for extended Lyman-Alpha emission around the z~6.6 QSO J0305-3150. After carefully subtracting the point-spread-function, we reach a nominal 5-sigma surface brightness limit of SB = 1.9x10$^{-18}$ erg/s/cm$^2$/arcsec$^2$ over a 1 arcsec$^2$ aperture, collapsing 5 wavelength slices centered at the expected location of the redshifted Lyman-Alpha emission (i.e. at 9256 Ang.). Current data suggest the presence (5-sigma, accounting for systematics) of a Lyman-Alpha nebula that extends for 9 kpc around the QSO. This emission is displaced and redshifted by 155 km/s with respect to the location of the QSO host galaxy traced by the [CII] emission line. The total luminosity is L = 3.0x10$^{42}$ erg/s. Our analysis suggests that this emission is unlikely to rise from optically thick clouds illuminated by the ionizing radiation of the QSO. It is more plausible that the Lyman-Alpha emission is due to fluorescence of the highly ionized optically thin gas. This scenario implies a high hydrogen volume density of n$_H$ ~ 6 cm$^{-3}$. In addition, we detect a Lyman-Alpha emitter (LAE) in the immediate vicinity of the QSO: i.e., with a projected separation of 12.5 kpc and a line-of-sight velocity difference of 560 km/s. The luminosity of the LAE is L = 2.1x10$^{42}$ erg/s and its inferred star-formation-rate is SFR ~ 1.3 M$_\odot$/yr. The probability of finding such a close LAE is one order of magnitude above the expectations based on the QSO-galaxy cross-correlation function. This discovery is in agreement with a scenario where dissipative interactions favour the rapid build-up of super-massive black holes at early Cosmic times.

We present the results of a survey for intervening MgII absorbers in the redshift range z $\simeq$2-6 in the foreground of four high redshift quasar spectra, 5.79$\le z_{em}\le$6.133, obtained with the ESO VLT X-shooter. We identify 24 absorbers at $\ge 5\sigma$ significance in the equivalent width range 0.117$\le W_{2796}\le 3.655\unicode{xC5}$ with the highest redshift absorber at $z=4.89031\pm4\times10^{-5}$. For weak ($W_{2796}<0.3\unicode{xC5}$) systems, we measure an incidence rate $dN/dz$=1.35$\pm$0.58 at <z>=2.34 and find that it almost doubles to $dN/dz$=2.58$\pm$0.67 by <z>=4.81. Weak absorbers exceeds the number expected from an exponential fit to stronger systems ($W_{2796}>0.3\unicode{xC5}$). We find that there must be significant evolution in the absorption halo properties of MgII absorbers with $W_{2796} >0.1\unicode{xC5}$ by <z>=4.77 and/or that they are associated with galaxies with luminosities beyond the limits of the current luminosity function at z $\sim$5. We find that the incidence rate of strong MgII absorbers ($W_{2796}>1.0\unicode{xC5}$) can be explained if they are associated with galaxies with $L\ge0.29L_{\ast}$ and/or their covering fraction increases. If they continue to only be associated with galaxies with $L\ge0.50L_{\ast}$ then their physical cross section ($\sigma_{phys}$) increases from 0.015 Mpc$^2$ at z=2.3 to 0.041 Mpc$^2$ at <z>=4.77. We measure $\Omega_{MgII}$=2.1$^{+6.3}_{-0.6}\times10^{-8}$, 1.9$^{+2.9}_{-0.2}\times10^{-8}$, 3.9$^{+7.1}_{-2.4}\times10^{-7}$ at <z>=2.48, 3.41, 4.77, respectively. At <z>=4.77, $\Omega_{MgII}$ exceeds the value expected from $\Omega_{HI}$ estimated from the global metallicity of DLAs at z $\simeq$4.85 by a factor of $\sim$44 suggesting that either MgII absorbers trace both ionised and neutral gas and/or are more metal rich than the average DLA at this redshift.

We present ALMA band 3 observations of the CO(6-5), CO(7-6), and [CI] 369micron emission lines in three of the highest redshift quasar host galaxies at 6.6<z<6.9. These measurements constitute the highest-redshift CO detections to date. The target quasars have previously been detected in [CII] 158micron emission and the underlying far-infrared (FIR) dust continuum. We detect (spatially unresolved, at a resolution of >2", or >14kpc) CO emission in all three quasar hosts. In two sources, we detect the continuum emission around 400micron (rest-frame), and in one source we detect [CI] at low significance. We derive molecular gas reservoirs of (1-3)x10^10 M_sun in the quasar hosts, i.e. approximately only 10 times the mass of their central supermassive black holes. The extrapolated [CII]-to-CO(1-0) luminosity ratio is 2500-4200, consistent with measurements in galaxies at lower redshift. The detection of the [CI] line in one quasar host galaxy and the limit on the [CI] emission in the other two hosts enables a first characterization of the physical properties of the interstellar medium in z~7 quasar hosts. In the sources, the derived global CO/[CII]/[CI] line ratios are consistent with expectations from photodissociation regions (PDR), but not X-ray dominated regions (XDR). This suggest that quantities derived from the molecular gas and dust emission are related to ongoing star-formation activity in the quasar hosts, providing further evidence that the quasar hosts studied here harbor intense starbursts in addition to their active nucleus.

The existence of massive ($10^{11}$ solar masses) elliptical galaxies by redshift z~4 (when the Universe was 1.5 billion years old) necessitates the presence of galaxies with star-formation rates exceeding 100 solar masses per year at z>6 (corresponding to an age of the Universe of less than 1 billion years). Surveys have discovered hundreds of galaxies at these early cosmic epochs, but their star-formation rates are more than an order of magnitude lower. The only known galaxies with very high star-formation rates at z>6 are, with only one exception, the host galaxies of quasars, but these galaxies also host accreting supermassive (more than $10^9$ solar masses) black holes, which probably affect the properties of the galaxies. Here we report observations of an emission line of singly ionized carbon ([CII] at a wavelength of 158 micrometres) in four galaxies at z>6 that are companions of quasars, with velocity offsets of less than 600 kilometers per second and linear offsets of less than 600 kiloparsecs. The discovery of these four galaxies was serendipitous; they are close to their companion quasars and appear bright in the far-infrared. On the basis of the [CII] measurements, we estimate star-formation rates in the companions of more than 100 solar masses per year. These sources are similar to the host galaxies of the quasars in [CII] brightness, linewidth and implied dynamical masses, but do not show evidence for accreting supermassive black holes. Similar systems have previously been found at lower redshift. We find such close companions in four out of twenty-five z>6 quasars surveyed, a fraction that needs to be accounted for in simulations. If they are representative of the bright end of the [CII] luminosity function, then they can account for the population of massive elliptical galaxies at z~4 in terms of cosmic space density.

We present a search for metal absorption line systems at the highest redshifts to date using a deep (30h) VLT/X-Shooter spectrum of the z = 7.084 quasi-stellar object (QSO) ULAS J1120+0641. We detect seven intervening systems at z > 5.5, with the highest-redshift system being a C IV absorber at z = 6.51. We find tentative evidence that the mass density of C IV remains flat or declines with redshift at z < 6, while the number density of C II systems remains relatively flat over 5 < z < 7. These trends are broadly consistent with models of chemical enrichment by star formation-driven winds that include a softening of the ultraviolet background towards higher redshifts. We find a larger number of weak ( W_rest < 0.3A ) Mg II systems over 5.9 < z < 7.0 than predicted by a power-law fit to the number density of stronger systems. This is consistent with trends in the number density of weak Mg II systems at z = 2.5, and suggests that the mechanisms that create these absorbers are already in place at z = 7. Finally, we investigate the associated narrow Si IV, C IV, and N V absorbers located near the QSO redshift, and find that at least one component shows evidence of partial covering of the continuum source.

Using MUSE on the ESO-VLT, we obtained a 4 hour exposure of the z=3.12 radio galaxy MRC0316-257. We detect features down to ~10^-19 erg/s/cm^2/arcsec^2 with the highest surface brightness regions reaching more than a factor of 100 higher. We find Ly-alpha emission out to ~250 kpc in projection from the active galactic nucleus (AGN). The emission shows arc-like morphologies arising at 150-250 kpc from the nucleus in projection with the connected filamentary structures reaching down into the circum-nuclear region. The most distant arc is offset by 700 km/s relative to circum-nuclear HeII 1640 emission, which we assume to be at the systemic velocity. As we probe emission closer to the nucleus, the filamentary emission narrows in projection on the sky, the relative velocity decreases to ~250 km/s, and line full-width at half maximum range from 300-700 km/s. From UV line ratios, the emission on scales of 10s of kpc from the nucleus along a wide angle in the direction of the radio jets is clearly excited by the radio jets and ionizing radiation of the AGN. Assuming ionization equilibrium, the more extended emission outside of the axis of the jet direction would require 100% or more illumination to explain the observed surface brightness. High speed (>300 km/s) shocks into rare gas would provide sufficiently high surface brightness. We discuss the possibility that the arcs of Ly-alpha emission represent accretion shocks and the filamentary emission represent gas flows into the halo, and compare our results with gas accretion simulations.

We present a 30h integration Very Large Telescope X-shooter spectrum of the Lyman series forest towards the $z = 7.084$ quasar ULAS J1120+0641. The only detected transmission at ${\rm S/N}>5$ is confined to seven narrow spikes in the Ly$\alpha$ forest, over the redshift range $5.858<z<6.122$, just longward of the wavelength of the onset of the Ly$\beta$ forest. There is also a possible detection of one further unresolved spike in the Ly$\beta$ forest at $z=6.854$, with ${\rm S/N}=4.5$. We also present revised Hubble Space Telescope F814W photometry of the source. The summed flux from the transmission spikes is in agreement with the F814W photometry, so all the transmission in the Lyman series forest may have been detected. There is a Gunn-Peterson (GP) trough in the Ly$\alpha$ forest from $z=6.122$ all the way to the quasar near zone at $z=7.04$. The trough, of comoving length $240\,h^{-1}$Mpc, is over twice as long as the next longest known GP trough. We combine the spectroscopic and photometric results to constrain the evolution of the Ly$\alpha$ effective optical depth with redshift, $\tau_{\rm GP}^{\rm eff}$ extending a similar analysis by Simpson et al. We find $\tau_{\rm GP}^{\rm eff} \propto (1+z)^{\xi}$ where $\xi = 11.2^{+0.4}_{-0.6}$, for $z > 5.5$. The data nevertheless provide only a weak limit on the volume-weighted hydrogen intergalactic (IGM) neutral fraction at $z\sim 6.5$, $x_{{\rm HI}} > 10^{-4}$, similar to limits at redshift $z\sim6$ from less distant quasars. The new observations cannot extend measurements of the neutral fraction of the IGM to higher values because absorption in the Ly$\alpha$ forest is already saturated near $z\sim6$. For higher neutral fractions, other methods such as measuring the red damping wing of the IGM will be required.

We present ALMA observations of the [CII] fine structure line and the underlying far-infrared (FIR) dust continuum emission in J1120+0641, the most distant quasar currently known (z=7.1). We also present observations targeting the CO(2-1), CO(7-6) and [CI] 369 micron lines in the same source obtained at the VLA and PdBI. We find a [CII] line flux of F_[CII]=1.11+/-0.10 Jy km/s and a continuum flux density of S_227GHz=0.53+/-0.04 mJy/beam, consistent with previous unresolved measurements. No other source is detected in continuum or [CII] emission in the field covered by ALMA (~25"). At the resolution of our ALMA observations (0.23", or 1.2 kpc, a factor ~70 smaller beam area compared to previous measurements), we find that the majority of the emission is very compact: a high fraction (~80%) of the total line and continuum flux is associated with a region 1-1.5 kpc in diameter. The remaining ~20% of the emission is distributed over a larger area with radius <4 kpc. The [CII] emission does not exhibit ordered motion on kpc-scales: applying the virial theorem yields an upper limit on the dynamical mass of the host galaxy of (4.3+/-0.9)x10^10 M_sun, only ~20x higher than the central black hole. The other targeted lines (CO(2-1), CO(7-6) and [CI]) are not detected, but the limits of the line ratios with respect to the [CII] emission imply that the heating in the quasar host is dominated by star formation, and not by the accreting black hole. The star-formation rate implied by the FIR continuum is 105-340 M_sun/yr, with a resulting star-formation rate surface density of ~100-350 M_sun/yr/kpc^2, well below the value for Eddington-accretion-limited star formation.

Bright quasars, observed when the Universe was less than one billion years old (z>5.5), are known to host massive black holes (~10$^{9}$ M$_{\odot}$), and are thought to reside in the center of massive dark matter overdensities. In this picture, overdensities of galaxies are expected around high redshift quasars. However, observations based on the detection of Lyman Break Galaxies (LBGs) around these quasars do not offer a clear picture: this may be due to the uncertain redshift constraints of LBGs, which are selected through broad-band filters only. To circumvent such uncertainties, we here perform a search for Lyman Alpha Emitting galaxies (LAEs) in the field of the quasar PSO J215.1512-16.0417 at z~5.73, through narrow band, deep imaging with FORS2 at the VLT. We study an area of 37 arcmin$^{2}$, i.e. ~206 comoving Mpc$^{2}$ at the redshift of the quasar. We find no evidence for an overdensity of LAEs in the quasar field with respect to blank field studies. Possible explanations for these findings include that our survey volume is too small, or that the strong ionizing radiation from the quasar hinders galaxy formation in its immediate proximity. Another possibility is that these quasars are not situated in the dense environments predicted by some simulations.

We present the results of a new spectroscopic survey for dusty intervening absorption systems, particularly damped Ly$\alpha$ absorbers (DLAs), towards reddened quasars. The candidate quasars are selected from mid-infrared photometry from the Wide-field Infrared Survey Explorer combined with optical and near-infrared photometry. Out of 1073 candidates, we secure low-resolution spectra for 108 using the Nordic Optical Telescope on La Palma, Spain. Based on the spectra, we are able to classify 100 of the 108 targets as quasars. A large fraction (50 %) is observed to have broad absorption lines (BALs). Moreover, we find 6 quasars with strange breaks in their spectra, which are not consistent with regular dust reddening. Using template fitting we infer the amount of reddening along each line of sight ranging from A(V)$\approx$0.1 mag to 1.2 mag (assuming an SMC extinction curve). In four cases, the reddening is consistent with dust exhibiting the 2175{\AA} feature caused by an intervening absorber, and for two of these, a MgII absorption system is observed at the best-fit absorption redshift. In the rest of the cases, the reddening is most likely intrinsic to the quasar. We observe no evidence for dusty DLAs in this survey. However, the large fraction of BAL quasars hampers the detection of absorption systems. Out of the 50 non-BAL quasars only 28 have sufficiently high redshift to detect Ly$\alpha$ in absorption.

Luminous quasars at z>5.6 can be studied in detail with the current generation of telescopes and provide us with unique information on the first gigayear of the universe. Thus far these studies have been statistically limited by the number of quasars known at these redshifts. Such quasars are rare and therefore wide-field surveys are required to identify them and multiwavelength data are needed to separate them efficiently from their main contaminants, the far more numerous cool dwarfs. In this paper, we update and extend the selection for z~6 quasars presented in Banados et al. (2014) using the Pan-STARRS1 (PS1) survey. We present the PS1 distant quasar sample, which currently consists of 124 quasars in the redshift range 5.6<z<6.7 that satisfy our selection criteria. Seventy-seven of these quasars have been discovered with PS1, and 63 of them are newly identified in this paper. We present composite spectra of the PS1 distant quasar sample. This sample spans a factor of ~20 in luminosity and shows a variety of emission line properties. The number of quasars at z>5.6 presented in this work almost double the quasars previously known at these redshifts, marking a transition phase from studies of individual sources to statistical studies of the high-redshift quasar population, which was impossible with earlier, smaller samples.

We present the serendipitous discovery of a projected pair of quasi-stellar objects (QSOs) with an angular separation of $\Delta\theta =4.50$ arcsec. The redshifts of the two QSOs are widely different: one, our programme target, is a QSO with a spectrum consistent with being a narrow line Seyfert 1 AGN at $z=2.05$. For this target we detect Lyman-$\alpha$, \ion{C}{4}, and \ion{C}{3]}. The other QSO, which by chance was included on the spectroscopic slit, is a Type 1 QSO at a redshift of $z=1.68$, for which we detect \ion{C}{4}, \ion{C}{3]} and \ion{Mg}{2}. We compare this system to previously detected projected QSO pairs and find that only about a dozen previously known pairs have smaller angular separation.

We analyze the near-infrared to UV data of 16 quasars with redshifts ranging from 0.71 $<$ $z$ $<$ 2.13 to investigate dust extinction properties. The sample presented in this work is obtained from the High $A_V$ Quasar (HAQ) survey. The quasar candidates were selected from the Sloan Digital Sky Survey (SDSS) and the UKIRT Infrared Deep Sky Survey (UKIDSS), and follow-up spectroscopy was carried out at the Nordic Optical Telescope (NOT) and the New Technology Telescope (NTT). To study dust extinction curves intrinsic to the quasars, from the HAQ survey we selected 16 cases where the Small Magellanic Cloud (SMC) law could not provide a good solution to the spectral energy distributions (SEDs). We derived the extinction curves using Fitzpatrick & Massa 1986 (FM) law by comparing the observed SEDs to the combined quasar template from Vanden Berk et al. 2001 and Glikman et al. 2006. The derived extinction, $A_V$, ranges from 0.2-1.0 mag. All the individual extinction curves of our quasars are steeper ($R_V=2.2$-2.7) than that of the SMC, with a weighted mean value of $R_V=2.4$. We derive an `average quasar extinction curve' for our sample by fitting SEDs simultaneously by using the weighted mean values of the FM law parameters and a varying $R_V$. The entire sample is well fit with a single best-fit value of $R_V=2.2\pm0.2$. The `average quasar extinction curve' deviates from the steepest Milky Way and SMC extinction curves at a confidence level $\gtrsim95\%$. Such steep extinction curves suggest a significant population of silicates to produce small dust grains. Moreover, another possibility could be that the large dust grains may have been destroyed by the activity of the nearby active galactic nuclei (AGN), resulting in steep extinction curves.

We present ALMA detections of the [CII] 158 micron emission line and the underlying far-infrared continuum of three quasars at 6.6<z<6.9 selected from the VIKING survey. The [CII] line fluxes range between 1.6-3.4 Jy km/s ([CII] luminosities ~(1.9-3.9)x10^9 L_sun). We measure continuum flux densities of 0.56-3.29 mJy around 158 micron (rest-frame), with implied far-infrared luminosities between (0.6-7.5)x10^12 L_sun and dust masses M_d=(0.7-24)x10^8 M_sun. In one quasar we derive a dust temperature of 30^+12_-9 K from the continuum slope, below the canonical value of 47 K. Assuming that the [CII] and continuum emission are powered by star formation, we find star-formation rates from 100-1600 M_sun/yr based on local scaling relations. The L_[CII]/L_FIR ratios in the quasar hosts span a wide range from (0.3-4.6)x10^-3, including one quasar with a ratio that is consistent with local star-forming galaxies. We find that the strength of the L_[CII] and 158 micron continuum emission in z>~6 quasar hosts correlate with the quasar's bolometric luminosity. In one quasar, the [CII] line is significantly redshifted by ~1700 km/s with respect to the MgII broad emission line. Comparing to values in the literature, we find that, on average, the MgII is blueshifted by 480 km/s (with a standard deviation of 630 km/s) with respect to the host galaxy redshift, i.e. one of our quasars is an extreme outlier. Through modeling we can rule out a flat rotation curve for our brightest [CII] emitter. Finally, we find that the ratio of black hole mass to host galaxy (dynamical) mass is higher by a factor 3-4 (with significant scatter) than local relations.

We present the results of our first year of quasar search in the on-going ESO public Kilo Degree Survey (KiDS) and VISTA Kilo-Degree Infrared Galaxy (VIKING) surveys. These surveys are among the deeper wide-field surveys that can be used to uncovered large numbers of z~6 quasars. This allows us to probe a more common population of z~6 quasars that is fainter than the well-studied quasars from the main Sloan Digital Sky Survey. From this first set of combined survey catalogues covering ~250 deg^2 we selected point sources down to Z_AB=22 that had a very red i-Z (i-Z>2.2) colour. After follow-up imaging and spectroscopy, we discovered four new quasars in the redshift range 5.8<z<6.0. The absolute magnitudes at a rest-frame wavelength of 1450 A are between -26.6 < M_1450 < -24.4, confirming that we can find quasars fainter than M^*, which at z=6 has been estimated to be between M^*=-25.1 and M^*=-27.6. The discovery of 4 quasars in 250 deg^2 of survey data is consistent with predictions based on the z~6 quasar luminosity function. We discuss various ways to push the candidate selection to fainter magnitudes and we expect to find about 30 new quasars down to an absolute magnitude of M_1450=-24. Studying this homogeneously selected faint quasar population will be important to gain insight into the onset of the co-evolution of the black holes and their stellar hosts.

The [CII] 158$\mu$m fine-structure line is known to trace regions of active star formation and is the main coolant of the cold, neutral atomic medium. In this \textit{Letter}, we report a strong detection of the [CII] line in the host galaxy of the brightest quasar known at $z>6.5$, the Pan-STARRS1 selected quasar PSO J036.5078+03.0498 (hereafter P036+03), using the IRAM NOEMA millimeter interferometer. Its [CII] and total far-infrared luminosities are $(5.8 \pm 0.7) \times 10^9 \,L_\odot$ and $(7.6\pm1.5) \times 10^{12}\,L_\odot$, respectively. This results in a $L_{[CII]} /L_{TIR}$ ratio of $\sim 0.8\times 10^{-3}$, which is at the high end for those found for active galaxies, though it is lower than the average found in typical main sequence galaxies at $z\sim 0$. We also report a tentative additional line which we identify as a blended emission from the $3_{22} - 3_{13}$ and $5_{23} - 4_{32}$ H$_2$O transitions. If confirmed, this would be the most distant detection of water emission to date. P036+03 rivals the current prototypical luminous J1148+5251 quasar at $z=6.42$, in both rest-frame UV and [CII] luminosities. Given its brightness and because it is visible from both hemispheres (unlike J1148+5251), P036+03 has the potential of becoming an important laboratory for the study of star formation and of the interstellar medium only $\sim 800\,$Myr after the Big Bang.

Radio-loud Active Galactic Nuclei at z~2-4 are typically located in dense environments and their host galaxies are among the most massive systems at those redshifts, providing key insights for galaxy evolution. Finding radio-loud quasars at the highest accessible redshifts (z~6) is important to study their properties and environments at even earlier cosmic time. They would also serve as background sources for radio surveys intended to study the intergalactic medium beyond the epoch of reionization in HI 21 cm absorption. Currently, only five radio-loud ($R=f_{\nu,5{\rm GHz}}/f_{\nu,4400\AA}>10$) quasars are known at z~6. In this paper we search for 5.5 < z < 7.2 quasars by cross-matching the optical Pan-STARRS1 and radio FIRST surveys. The radio information allows identification of quasars missed by typical color-based selections. While we find no good 6.4 < z <7.2 quasar candidates at the sensitivities of these surveys, we discover two new radio-loud quasars at z~6. Furthermore, we identify two additional z~6 radio-loud quasars which were not previously known to be radio-loud, nearly doubling the current z~6 sample. We show the importance of having infrared photometry for z>5.5 quasars to robustly classify them as radio-quiet or radio-loud. Based on this, we reclassify the quasar J0203+0012 (z=5.72), previously considered radio-loud, to be radio-quiet. Using the available data in the literature, we constrain the radio-loud fraction of quasars at z~6, using the Kaplan--Meier estimator, to be $8.1^{+5.0}_{-3.2}\%$. This result is consistent with there being no evolution of the radio-loud fraction with redshift, in contrast to what has been suggested by some studies at lower redshifts.

We present deep MUSE integral-field unit (IFU) spectroscopic observations of the giant (~150 x 80 kpc) Ly-alpha halo around the z=4.1 radio galaxy TNJ J1338-1942. This 9-hr observation maps the two-dimensional kinematics of the Ly-alpha emission across the halo. We identify two HI absorbers which are seen against the Ly-alpha emission, both of which cover the full 150 x 80 kpc extent of the halo and so have covering fractions ~1. The stronger and more blue-shifted absorber (dv~1200 km/s) has dynamics that mirror that of the underlying halo emission and we suggest that this high column material (n(HI) ~ 10^19.4 /cm^2), which is also seen in CIV absorption, represents an out-flowing shell that has been driven by the AGN (or star formation) within the galaxy. The weaker (n(HI)~10^14 /cm^2) and less blue shifted (dv~500 km/s) absorber most likely represents material in the cavity between the out-flowing shell and the Ly-alpha halo. We estimate that the mass in the shell must be of order 10^10 Msol -- a significant fraction of the ISM from a galaxy at z=4. The large scales of these coherent structures illustrate the potentially powerful influence of AGN feedback on the distribution and energetics of material in their surroundings. Indeed, the discovery of high-velocity (~1000 km/s), group-halo-scale (i.e. >150 kpc) and mass-loaded winds in the vicinity of the central radio source are broadly in agreement with the requirements of models that invoke AGN-driven outflows to regulate star formation and black-hole growth in massive galaxies at early times.

Luminous distant quasars are unique probes of the high redshift intergalactic medium (IGM) and of the growth of massive galaxies and black holes in the early universe. Absorption due to neutral Hydrogen in the IGM makes quasars beyond a redshift of z~6.5 very faint in the optical $z$-band, thus locating quasars at higher redshifts require large surveys that are sensitive above 1 micron. We report the discovery of three new z>6.5 quasars, corresponding to an age of the universe of <850 Myr, selected as z-band dropouts in the Pan-STARRS1 survey. This increases the number of known z>6.5 quasars from 4 to 7. The quasars have redshifts of z=6.50, 6.52, and 6.66, and include the brightest z-dropout quasar reported to date, PSO J036.5078+03.0498 with M_1450=-27.4. We obtained near-infrared spectroscopy for the quasars and from the MgII line we estimate that the central black holes have masses between 5x10^8 and 4x10^9 M_sun, and are accreting close to the Eddington limit (L_Bol/L_Edd=0.13-1.2). We investigate the ionized regions around the quasars and find near zone radii of R_NZ=1.5-5.2 proper Mpc, confirming the trend of decreasing near zone sizes with increasing redshift found for quasars at 5.7<z<6.4. By combining R_NZ of the PS1 quasars with those of 5.7<z<7.1 quasars in the literature, we derive a luminosity corrected redshift evolution of R_NZ,corrected=(7.2+/-0.2)-(6.1+/-0.7)x(z-6) Mpc. However, the large spread in R_NZ in the new quasars implies a wide range in quasar ages and/or a large variation in the neutral Hydrogen fraction along different lines of sight.

Quasi-stellar objects (QSOs) whose spectral energy distributions (SEDs) are reddened by dust either in their host galaxies or in intervening absorber galaxies are to a large degree missed by optical color selection criteria like the one used by the Sloan Digital Sky Survey (SDSS). To overcome this bias against red QSOs, we employ a combined optical and near-infrared color selection. In this paper, we present a spectroscopic follow-up campaign of a sample of red candidate QSOs which were selected from the SDSS and the UKIRT Infrared Deep Sky Survey (UKIDSS). The spectroscopic data and SDSS/UKIDSS photometry are supplemented by mid-infrared photometry from the Wide-field Infrared Survey Explorer. In our sample of 159 candidates, 154 (97%) are confirmed to be QSOs. We use a statistical algorithm to identify sightlines with plausible intervening absorption systems and identify nine such cases assuming dust in the absorber similar to Large Magellanic Cloud sightlines. We find absorption systems toward 30 QSOs, 2 of which are consistent with the best-fit absorber redshift from the statistical modeling. Furthermore, we observe a broad range in SED properties of the QSOs as probed by the rest-frame 2 {\mu}m flux. We find QSOs with a strong excess as well as QSOs with a large deficit at rest-frame 2 {\mu}m relative to a QSO template. Potential solutions to these discrepancies are discussed. Overall, our study demonstrates the high efficiency of the optical/near-infrared selection of red QSOs.

We report the discovery of an extremely long ($\sim$110 Mpc/$h$) and dark ($\tau_{\rm eff} \gtrsim 7$) Ly$\alpha$ trough extending down to $z \simeq 5.5$ towards the $z_{\rm em} \simeq 6.0$ quasar ULAS J0148+0600. We use these new data in combination with Ly$\alpha$ forest measurements from 42 quasars at $4.5 \le z_{\rm em} \le 6.4$ to conduct an updated analysis of the line-of-sight variance in the intergalactic Ly$\alpha$ opacity over $4 \le z \le 6$. We find that the scatter in transmission among lines of sight near $z \sim 6$ significantly exceeds theoretical expectations for either a uniform ultraviolet background (UVB) or simple fluctuating UVB models in which the mean free path to ionizing photons is spatially invariant. The data, particularly near $z \simeq 5.6$-5.8, instead require fluctuations in the volume-weighted hydrogen neutral fraction that are a factor of 3 or more beyond those expected from density variations alone. We argue that these fluctuations are most likely driven by large-scale variations in the mean free path, consistent with expectations for the final stages of inhomogeneous hydrogen reionization. Even by $z \simeq 5.6$, however, a large fraction of the data are consistent with a uniform UVB, and by $z \sim 5$ the data are fully consistent with opacity fluctuations arising solely from the density field. This suggests that while reionization may be ongoing at $z \sim 6$, it has fully completed by $z \sim 5$.

We present new observations of the highest-redshift quasar known, ULAS J1120+0641, redshift $z=7.084$, obtained in the optical, at near-, mid-, and far-infrared wavelengths, and in the sub-mm. We combine these results with published X-ray and radio observations to create the multiwavelength spectral energy distribution (SED), with the goals of measuring the bolometric luminosity $L_{\rm bol}$, and quantifying the respective contributions from the AGN and star formation. We find three components are needed to fit the data over the wavelength range $0.12-1000\,\mu$m: the unobscured quasar accretion disk and broad-line region, a dusty clumpy AGN torus, and a cool 47K modified black body to characterise star formation. Despite the low signal-to-noise ratio of the new long-wavelength data, the normalisation of any dusty torus model is constrained within $\pm40\%$. We measure a bolometric luminosity $L_{\rm bol}=2.6\pm0.6\times10^{47}\,$erg$\,$s$^{-1}=6.7 \pm 1.6\times10^{13}L_{\odot}$, to which the three components contribute $31\%,32\%,3\%$, respectively, with the remainder provided by the extreme UV $<0.12\,\mu$m. We tabulate the best-fit model SED. We use local scaling relations to estimate a star formation rate (SFR) in the range $60-270\,{\rm M}_\odot$/yr from the [C$\,{\scriptsize \rm II}$] line luminosity and the $158\,\mu$m continuum luminosity. An analysis of the equivalent widths of the [C$\,{\scriptsize \rm II}$] line in a sample of $z>5.7$ quasars suggests that these indicators are promising tools for estimating the SFR in high-redshift quasars in general. At the time observed the black hole was growing in mass more than 100 times faster than the stellar bulge, relative to the mass ratio measured in the local universe, i.e. compared to ${M_{\rm BH}}/{M_{\rm bulge}} \simeq 1.4\times10^{-3}$, for ULAS J1120+0641 we measure ${\dot{M}_{\rm BH}}/{\dot{M}_{\rm bulge}} \simeq 0.2$.

We present APEX LABOCA 870 micron observations of the field around the high-redshift radio galaxy MRC1138-262 at z=2.16. We detect 16 submillimeter galaxies in this ~140 square arcmin bolometer map with flux densities in the range 3-11 mJy. The raw number counts indicate a density of submillimeter galaxies (SMGs) that is up to four times that of blank field surveys. Based on an exquisite multiwavelength database, including VLA 1.4 GHz radio and infrared observations, we investigate whether these sources are members of the protocluster structure at z=2.2. Using Herschel PACS+SPIRE and Spitzer MIPS photometry, we derive reliable far-infrared photometric redshifts for all sources. Follow-up VLT ISAAC and SINFONI near-infrared spectra confirm that four of these SMGs have redshifts of z=2.2. We also present evidence that another SMG in this field, detected earlier at 850 micron, has a counterpart that exhibits Halpha and CO(1-0) emission at z=2.15. Including the radio galaxy and two SMGs with far-IR photometric redshifts at z=2.2, we conclude that at least eight submm sources are part of the protocluster at z=2.16 associated with the radio galaxy MRC1138-262. We measure a star formation rate density SFRD ~1500 Msun yr^-1 Mpc^-3, four magnitudes higher than the global SFRD of blank fields at this redshift. Strikingly, these eight sources are concentrated within a region of 2 Mpc (the typical size of clusters in the local universe) and are distributed within the filaments traced by the Halpha emitters at z=2.2. This concentration of massive, dusty starbursts is not centered on the submillimeter-bright radio galaxy which could support the infalling of these sources into the cluster center. Approximately half (6/11) of the SMGs that are covered by the Halpha imaging data are associated with Halpha emitters, demonstrating the potential of tracing SMG counterparts with this population (abridged).

We present the analysis of optical and near-infrared spectra of the only four $z>6.5$ quasars known to date, discovered in the UKIDSS-LAS and VISTA-VIKING surveys. Our data-set consists of new VLT/X-Shooter and Magellan/FIRE observations. These are the best optical/NIR spectroscopic data that are likely to be obtained for the $z>6.5$ sample using current $6$ - $10$ m facilities. We estimate the black hole mass, the Eddington ratio, and the SiIV/CIV, CIII]/CIV, and FeII/MgII emission-line flux ratios. We perform spectral modeling using a procedure that allows us to derive a probability distribution for the continuum components and to obtain the quasar properties weighted upon the underlying distribution of continuum models. The $z>6.5$ quasars show the same emission properties as their counterparts at lower redshifts. The $z>6.5$ quasars host black holes with masses of $\sim 10^9$ M$_{\odot}$ that are accreting close to the Eddington luminosity ($\langle{\rm log} (L_{\rm Bol}/L_{\rm Edd})\rangle= -0.4\pm0.2$), in agreement with what has been observed for a sample of $4.0<z<6.5$ quasars. By comparing the SiIV/CIV and CIII]/CIV flux ratios with the results obtained from luminosity-matched samples at $z\sim6$ and $2\leq z\leq4.5$, we find no evidence of evolution of the line ratios with cosmic time. We compare the measured FeII/MgII flux ratios with those obtained for a sample of $4.0<z<6.4$ sources. The two samples are analyzed using a consistent procedure. There is no evidence that the FeII/MgII flux ratio evolves between $z=7$ and $z=4$. Under the assumption that the FeII/MgII traces the Fe/Mg abundance ratio, this implies the presence of major episodes of chemical enrichment in the quasar hosts in the first $\sim0.8$ Gyr after the Big Bang.

We present optical and near-infrared imaging of the field of the z=7.0842 quasar ULAS J112001.48+064124.3 taken with the Hubble Space Telescope. We use these data to search for galaxies that may be physically associated with the quasar, using the Lyman break technique, and find three such objects, although the detection of one in Spitzer Space Telescope imaging strongly suggests it lies at z~2. This is consistent with the field luminosity function and indicates that there is no excess of >L* galaxies within 1 Mpc of the quasar. A detection of the quasar shortward of the Ly-alpha line is consistent with the previously observed evolution of the intergalactic medium at z>5.5.

High-redshift quasars are currently the only probes of the growth of supermassive black holes and potential tracers of structure evolution at early cosmic time. Here we present our candidate selection criteria from the Panoramic Survey Telescope & Rapid Response System 1 and follow-up strategy to discover quasars in the redshift range 5.7<z<6.2. With this strategy we discovered eight new 5.7<z<6.0 quasars, increasing the number of known quasars at z>5.7 by more than 10%. We additionally recovered 18 previously known quasars. The eight quasars presented here span a large range of luminosities (-27.3 < M_{1450} < -25.4; 19.6 < z_ps1 < 21.2) and are remarkably heterogeneous in their spectral features: half of them show bright emission lines whereas the other half show a weak or no Ly$\alpha$ emission line (25% with rest-frame equivalent width of the Ly$\alpha$ + Nv line lower than 15{\AA}). We find a larger fraction of weak-line emission quasars than in lower redshift studies. This may imply that the weak-line quasar population at the highest redshifts could be more abundant than previously thought. However, larger samples of quasars are needed to increase the statistical significance of this finding.

We present X-ray imaging and spectroscopy of the redshift z=7.084 radio-quiet quasar ULAS J112001.48+064124.3 obtained with Chandra and XMM-Newton. The quasar is detected as a point source with both observatories. The Chandra observation provides a precise position, confirming the association of the X-ray source and the quasar, while a sufficient number of photons is detected in the XMM-Newton observation to yield a meaningful X-ray spectrum. In the XMM-Newton observation the quasar has a 2-10 keV luminosity of 4.7+-0.9 times 10^44 erg/s and a spectral slope alpha = 1.6+0.4/-0.3 (where f_nu is proportional to nu^-alpha). The quasar appears to have dimmed in the 15 months between the two observations, with a 2-10 keV luminosity of $1.8+1.0/-0.7 times 10^45 erg/s during the Chandra observation. We derive optical to X-ray spectral slopes alpha_OX of 1.76+-0.07 and 1.54+0.09/-0.08 at the times of the XMM-Newton and Chandra observations respectively, consistent with the range of alpha_OX found in other quasars of comparable ultraviolet luminosity. The very soft X-ray spectrum suggests that the quasar is accreting above the Eddington rate, L/L_Edd = 5+15/-4, compared to L/L_Edd = 1.2+0.6/-0.5 derived from the rest-frame ultraviolet. Super-Eddington accretion would help to reduce the discrepancy between the age of the quasar implied by the small size of the ionized near zone in which it sits (<10^7 years), and the characteristic e-folding time (2.5 times 10^7 years if L/L_Edd=2). Such super-Eddington accretion would also alleviate the challenging constraints on the seed black hole mass provided that the quasar has been rapidly accreting throughout its history. The remnant of an individual population III star is a plausible progenitor if an average L/L_Edd>1.46 has been maintained over the quasar's lifetime.

Studying quasars at the highest redshifts can constrain models of galaxy and black hole formation, and it also probes the intergalactic medium in the early universe. Optical surveys have to date discovered more than 60 quasars up to z~6.4, a limit set by the use of the z-band and CCD detectors. Only one z>6.4 quasar has been discovered, namely the z=7.08 quasar ULAS J1120+0641, using near-infrared imaging. Here we report the discovery of three new z>6.4 quasars in 332 square degrees of the Visible and Infrared Survey Telescope for Astronomy Kilo-degree Infrared Galaxy (VIKING) survey, thus extending the number from 1 to 4. The newly discovered quasars have redshifts of z=6.60, 6.75, and 6.89. The absolute magnitudes are between -26.0 and -25.5, 0.6-1.1 mag fainter than ULAS J1120+0641. Near-infrared spectroscopy revealed the MgII emission line in all three objects. The quasars are powered by black holes with masses of ~(1-2)x10^9 M_sun. In our probed redshift range of 6.44<z<7.44 we can set a lower limit on the space density of supermassive black holes of \rho(M_BH>10^9 M_sun) > 1.1x10^(-9) Mpc^(-3). The discovery of three quasars in our survey area is consistent with the z=6 quasar luminosity function when extrapolated to z~7. We do not find evidence for a steeper decline in the space density of quasars with increasing redshift from z=6 to z=7.

We present 1-2 GHz Very Large Array A-configuration continuum observations on the highest redshift quasar known to date, the $z=7.085$ quasar ULAS J112001.48+064124.3. The results show no radio continuum emission at the optical position of the quasar or its vicinity at a level of $\geq 3\sigma$ or $23.1 \mu$Jy beam$^{-1}$. This $3\sigma$ limit corresponds to a rest frame 1.4 GHz luminosity density limit of $L_{\nu,1.4\,GHz} < 1.76 \times 10^{24}$ W Hz$^{-1}$ for a spectral index of $\alpha=0$, and $L_{\nu,1.4\,GHz} < 1.42 \times 10^{25}$ W Hz$^{-1}$ for a spectral index of $\alpha=-1$. The rest-frame 1.4 GHz luminosity limits are $L_{rad} < 6.43 \times 10^6 L_{\odot}$ and $L_{\rm rad} < 5.20 \times 10^7 L_{\odot}$ for $\alpha=0$ and $\alpha=-1$, respectively. The derived limits for the ratio of the rest frame 1.4 GHz luminosity density to the $B$-band optical luminosity density are $R\rlap{}_{1.4}^{*} < 0.53$ and $< 4.30$ for the above noted spectral indices, respectively. Given our upper limits on the radio continuum emission and the radio-to-optical luminosity ratio, we conclude that this quasar is radio-quiet and located at the low end of the radio quiet distribution of high redshift ($z \gtrsim 6$) quasars.

This paper presents the first results of a far-infrared search for protocluster-associated galaxy overdensities using the SPIRE instrument on-board the {\it Herschel} Space Observatory. Large ($\sim$400 arcmin$^{2}$) fields surrounding 26 powerful high-redshift radio galaxies ($2.0 < z < 4.1$; $L_{\rm 500 MHz} > 10^{28.5}$ WHz$^{-1}$) are mapped at 250, 350 and 500\mic to give a unique wide-field sample. On average the fields have a higher than expected, compared to blank fields, surface density of 500\mic sources within 6 comoving Mpc of the radio galaxy. The analysis is then restricted to potential protocluster members only, which are identified using a far-infrared colour selection; this reveals significant overdensities of galaxies in 2 fields, neither of which are previously known protoclusters. The probability of finding 2 overdensities of this size by chance, given the number of fields observed is $5 \times 10^{-4}$. Overdensities here exist around radio galaxies with $L_{\rm 500 MHz} \gtrsim 10^{29}$ WHz$^{-1}$ and $z < 3$. The radial extent of the average far-infrared overdensity is found to be $\sim$6 comoving Mpc. Comparison with predictions from numerical simulations shows that the overdensities are consistent with having masses $> 10^{14}$Msolar. However, the large uncertainty in the redshift estimation means that it is possible that these far-infrared overdensities consist of several structures across the redshift range searched.

High-redshift quasars are believed to reside in massive halos in the early universe and should therefore be located in fields with overdensities of galaxies, which are thought to evolve into galaxy clusters seen in the local universe. However, despite many efforts, the relationship between galaxy overdensities and z~6 quasars is ambiguous. This can possibly be attributed to the difficulty of finding galaxies with accurate redshifts in the vicinity of z~6 quasars. So far, overdensity searches around z~6 quasars have been based on studies of Lyman break galaxies (LBGs), which probe a redshift range of Delta z ~ 1. This range is large enough to select galaxies that may not be physically related to the quasar. We use deep narrow- and broadband imaging to study the environment of the z=5.72 quasar ULAS J0203+0012. The redshift range probed by our narrow-band selection of Lyman alpha emitters (LAEs) is Delta z ~ 0.1, which is significantly narrower than the LBG searches. This is the first time that LAEs were searched for near a z~6 quasar, in an effort to provide clues about the environments of quasars at the end of the epoch of reionization. We find no enhancement of LAEs in the surroundings of ULAS J0203+0012 in comparison with blank fields. We explore different explanations and interpretations for this non-detection of a galaxy overdensity, including that (1) the strong ionization from the quasar may prevent galaxy formation in its immediate vicinity and (2) high-redshift quasars may not reside in the center of the most massive dark matter halos.

We present the results of a search for red QSOs using a selection based on optical imaging from SDSS and near-infrared imaging from UKIDSS. For a sample of 58 candidates 46 (79%) are confirmed to be QSOs. The QSOs are predominantly dust-reddened except a handul at redshifts z>3.5. The dust is most likely located in the QSO host galaxies. 4 (7%) of the candidates turned out to be late-type stars, and another 4 (7%) are compact galaxies. The remaining 4 objects we could not identify. In terms of their optical spectra the QSOs are similar to the QSOs selected in the FIRST-2MASS red Quasar survey except they are on average fainter, more distant and only two are detected in the FIRST survey. We estimate the amount of extinction using the SDSS QSO template reddened by SMC-like dust. It is possible to get a good match to the observed (restframe ultraviolet) spectra, but for nearly all the reddened QSOs it is not possible to match the near-IR photometry from UKIDSS. The likely reasons are that the SDSS QSO template is too red at optical wavelengths due to contaminating host galaxy light and that the assumed SMC extinction curve is too shallow. Our survey has demonstrated that selection of QSOs based on near-IR photometry is an efficent way to select QSOs, including reddened QSOs, with only small contamination from late-type stars and compact galaxies. This will be useful with ongoing and future wide-field near-IR surveys such as the VISTA and EUCLID surveys. [Abridged]

We present a new sample of purely near infra-red selected KVega<16.5 extremely red ((J-K)>2.5) quasar candidates at z~2 from 900 sq deg of data in the UKIDSS Large Area Survey (LAS). Five of these are spectroscopically confirmed to be heavily reddened Type 1 AGN with broad emission lines bringing our total sample of reddened quasars from the UKIDSS-LAS to 12 at z=1.4-2.7. At these redshifts, Ha (6563A) is in the K-band. However, the mean Ha equivalent width of the reddened quasars is only ten per cent larger than that of the optically selected population and cannot explain the extreme colours. Instead, dust extinction of Av~2-6 mags is required to reproduce the continuum colours of our sources. This is comparable to the dust extinctions seen in submillimetre galaxies at similar redshifts. We argue that the AGN are likely being observed in a relatively short-lived breakout phase when they are expelling gas and dust following a massive starburst, subsequently turning into UV-luminous quasars. Some of our quasars show direct evidence for strong outflows (v~800-1000 km/s) affecting the Ha line consistent with this scenario. We predict that a larger fraction of reddened quasar hosts are likely to be submillimetre bright compared to the UV-luminous quasar population. We use our sample to place new constraints on the fraction of obscured Type 1 AGN likely to be missed in optical surveys. Taken at face-value our findings suggest that the obscured fraction depends on quasar luminosity. The space density of obscured quasars is ~5x that inferred for UV-bright quasars from the SDSS luminosity function at Mi<-30 but seems to drop at lower luminosities even accounting for various sources of incompleteness in our sample (abridged).

We search for extended Ly-alpha emission around two z>6 quasars, SDSS J1030+0524 (z=6.309) and SDSS J1148+5251 (z=6.419) using WFC3 narrow-band filters on board the Hubble Space Telescope. For each quasar, we collected two deep, narrow-band images, one sampling the Ly-alpha line+continuum at the quasar redshifts and one of the continuum emission redwards of the line. After carefully modeling the Point Spread Function, we find no evidence for extended Ly-alpha emission. These observations set 2-sigma limits of L(Ly-alpha, extended) < 3.2 x 10^{44} erg/s for J1030+0524 and L(Ly-alpha, extended) < 2.5 x 10^{44} erg/s for J1148+5251. Given the star formation rates typically inferred from (rest-frame) far-infrared measurements of z~6 quasars, these limits are well below the intrinsic bright Ly-alpha emission expected from the recombination of gas photoionized by the quasars or by the star formation in the host galaxies, and point towards significant Ly-alpha suppression or dust attenuation. However, small extinction values have been observed along the line of sight to the nuclei, thus reddening has to be coupled with other mechanisms for Ly-alpha suppression (e.g., resonance scattering). No Ly-alpha emitting companions are found, down to a 5-sigma sensitivity of ~ 1 x 10^{-17} erg/s/cm^2/arcsec^2 (surface brightness) and ~ 5 x 10^{-17} erg/s/cm^2 (assuming point sources).

Using the IRAM Plateau de Bure Interferometer, we report the detection of the 158 micron [CII] emission line and underlying dust continuum in the host galaxy of the quasar ULAS J112001.48+064124.3 (hereafter J1120+0641) at z=7.0842+/-0.0004. This is the highest redshift detection of the [CII] line to date, and allows us to put first constraints on the physical properties of the host galaxy. The [CII] line luminosity is (1.2+/-0.2)x10^9 Lsun, which is a factor ~4 lower than observed in a luminous quasar at z=6.42 (SDSS J1148+5251). The underlying far-infrared (FIR) continuum has a flux density of 0.61+/-0.16 mJy, similar to the average flux density of z~6 quasars that were not individually detected in the rest-frame FIR. Assuming the FIR luminosity of L_FIR = 5.8x10^11-1.8x10^12 Lsun is mainly powered by star-formation, we derive a star-formation rate in the range 160-440 Msun/yr and a total dust mass in the host galaxy of 6.7x10^7-5.7x10^8 Msun (both numbers have significant uncertainties given the unknown nature of dust at these redshifts). The [CII] line width of sigma_V=100+/-15 km/s is among the smallest observed when compared to the molecular line widths detected in z~6 quasars. Both the [CII] and dust continuum emission are spatially unresolved at the current angular resolution of 2.0x1.7 arcsec^2 (corresponding to 10x9 kpc^2 at the redshift of J1120+0641).

We present spectroscopic follow-up observations of Lyman Break Galaxies (LBGs) selected in the field surrounding the radio galaxy MRC0316-257 at z~3.13 (0316). Robust spectroscopic redshifts are determined for 20 out of 24 objects. Three of the spectroscopically confirmed galaxies have 3.12<z<3.13 indicating that these objects reside in a protocluster structure previously found around the radio galaxy. An additional 5 objects are found 1600 km/s blue-shifted with respect to the main protocluster structure. This is in addition to three [OIII] emitters found at this redshift in a previous study. This is further evidence that a structure exists directly in front of the 0316 protocluster. We estimate that the foreground structure is responsible for half of the surface overdensity of LBGs found in the field as a whole. The foreground structure is associated with a strong surface density peak 1.4 Mpc to the North-West of the radio galaxy and a 2D Kolmogorov-Smirnov test indicates that the spatial distributions of the 0316 and foreground galaxies differ at the 3 sigma level. In addition, we compare the properties of protocluster, foreground structure and field galaxies, but we find no significant differences. In terms of the nature of the two structures, a merger scenario is a possible option. Simple merger dynamics indicates that the observed relative velocity of 1600 km/s can be reproduced if the two structures have masses of ~5x10^14 Msun and have starting separations of around 2.5 to 3 Mpc. It is also possible that the foreground structure is unrelated to the 0316 protocluster in which case the two structures will not interact before z=0.

The European Southern Observatory's (ESO) Visible and Infrared Survey Telescope for Astronomy (VISTA) is a 4-m class survey telescope for wide-field near-infrared imaging. VISTA is currently running a suite of six public surveys, which will shortly deliver their first Europe wide public data releases to ESO. The VISTA Kilo-degree Infrared Galaxy Survey (VIKING) forms a natural intermediate between current wide shallow, and deeper more concentrated surveys, by targeting two patches totalling 1500 sq.deg in the northern and southern hemispheres with measured 5-sigma limiting depths of Z ~ 22.4, Y ~ 21.4, J ~ 20.9, H ~ 19.9 and Ks ~19.3 (Vega). This architecture forms an ideal working parameter space for the discovery of a significant sample of 6.5 <= z <= 7.5 quasars. In the first data release priority has been placed on small areas encompassing a number of fields well sampled at many wavelengths, thereby optimising science gains and synergy whilst ensuring a timely release of the first products. For rare object searches e.g. high-z quasars, this policy is not ideal since photometric selection strategies generally evolve considerably with the acquisition of data. Without a reasonably representative data set sampling many directions on the sky it is not clear how a rare object search can be conducted in a highly complete and efficient manner. In this paper, we alleviate this problem by supplementing initial data with a realistic model of the spatial, luminosity and colour distributions of sources known to heavily contaminate photometric quasar selection spaces, namely dwarf stars of spectral type M, L and T. We use this model along with a subset of available data to investigate contamination of quasar selection space by cool stars and galaxies and lay down a set of benchmark selection constraints that limit contamination to reasonable levels whilst maintaining high completeness...

Aims. Ly-alpha emitters (LAEs) can be detected out to very high redshifts during the epoch of reionization. The evolution of the LAE luminosity function with redshift is a direct probe of the Ly-alpha transmission of the intergalactic medium (IGM), and therefore of the IGM neutral-hydrogen fraction. Measuring the Ly-alpha luminosity function (LF) of LAEs at redshift z = 7.7 therefore allows us to constrain the ionizing state of the Universe at this redshift. Methods. We observed three 7.5'x7.5' fields with the HAWK-I instrument at the VLT with a narrow band filter centred at 1.06 $\mu$m and targeting LAEs at redshift z ~ 7.7. The fields were chosen for the availability of multiwavelength data. One field is a galaxy cluster, the Bullet Cluster, which allowed us to use gravitational amplification to probe luminosities that are fainter than in the field. The two other fields are subareas of the GOODS Chandra Deep Field South and CFHTLS-D4 deep field. We selected z=7.7 LAE candidates from a variety of colour criteria, in particular from the absence of detection in the optical bands. Results. We do not find any LAE candidates at z = 7.7 in ~2.4 x 10^4 Mpc^3 down to a narrow band AB magnitude of ~ 26, which allows us to infer robust constraints on the Ly-alpha LAE luminosity function at this redshift. Conclusions. The predicted mean number of objects at z = 6.5, derived from somewhat different LFs of Hu et al. (2010), Ouchi et al. (2010), and Kashikawa et al. (2011) are 2.5, 13.7, and 11.6, respectively. Depending on which of these LFs we refer to, we exclude a scenario with no evolution from z = 6.5 to z = 7.7 at 85% confidence without requiring a strong change in the IGM Ly-alpha transmission, or at 99% confidence with a significant quenching of the IGM Ly-alpha transmission, possibly from a strong increase in the high neutral-hydrogen fraction between these two redshifts.

The quasar ULAS J1120+0641 at redshift z=7.085 has a highly ionised near zone which is smaller than those around quasars of similar luminosity at z~6. The spectrum also exhibits evidence for a damping wing extending redward of the systemic Lya redshift. We use radiative transfer simulations in a cosmological context to investigate the implications for the ionisation state of the inhomogeneous IGM surrounding this quasar. Our simulations show that the transmission profile is consistent with an IGM in the vicinity of the quasar with a volume averaged HI fraction of f_HI>0.1 and that ULAS J1120+0641 has been bright for 10^6--10^7 yr. The observed spectrum is also consistent with smaller IGM neutral fractions, f_HI ~ 10^-3--10-4, if a damped Lya system in an otherwise highly ionised IGM lies within 5 proper Mpc of the quasar. This is, however, predicted to occur in only ~5 per cent of our simulated sight-lines for a bright phase of 10^6--10^7 yr. Unless ULAS J1120+0641 grows during a previous optically obscured phase, the low age inferred for the quasar adds to the theoretical challenge of forming a 2x10^9 M_sol black hole at this high redshift.

The intergalactic medium was not completely reionized until approximately a billion years after the Big Bang, as revealed by observations of quasars with redshifts of less than 6.5. It has been difficult to probe to higher redshifts, however, because quasars have historically been identified in optical surveys, which are insensitive to sources at redshifts exceeding 6.5. Here we report observations of a quasar (ULAS J112001.48+064124.3) at a redshift of 7.085, which is 0.77 billion years after the Big Bang. ULAS J1120+0461 had a luminosity of 6.3x10^13 L_Sun and hosted a black hole with a mass of 2x10^9 M_Sun (where L_Sun and M_Sun are the luminosity and mass of the Sun). The measured radius of the ionized near zone around ULAS J1120+0641 is 1.9 megaparsecs, a factor of three smaller than typical for quasars at redshifts between 6.0 and 6.4. The near zone transmission profile is consistent with a Ly alpha damping wing, suggesting that the neutral fraction of the intergalactic medium in front of ULAS J1120+0641 exceeded 0.1.

We present the first results obtained using a tunable narrowband filter in the search for high-z protoclusters. Using the recently commissioned red tunable filter on the Gran Telescopio Canarias we have searched for Lya emitters in a 75 arcmin^2 field centered on the z=4.413 radio galaxy 6C0140+326. With three different wavelength tunings we find a total of 27 unique candidate Lya emitters. The availability of three different wavelength tunings allows us to make estimates of the redshifts for each of the objects. It also allows us to separate a possible protocluster from structure in the immediate foreground. This division shows that the foreground region contains significantly fewer Lya emitters. Also, the spatial distribution of the objects in the protocluster field deviates from a random distribution at the 2.5 sigma level. The observed redshift distribution of the emitters is different from the expected distribution of a blank field at the ~3 sigma level, with the Lya emitters concentrated near the radio galaxy at z>4.38. The 6C0140+326 field is denser by a factor of 9+/-5 than a blank field, and the number density of Lya emitters close to the radio galaxy is similar to that of the z~4.1 protocluster around TNJ1338-1942. We thus conclude that there is an overdensity of Lya emitters around the radio galaxy 6C0140+326. This is one of few known overdensities at such a high redshift.

The environment of the high-z radio galaxy PKS 1138-262 at z~2.2 is a prime example of a forming galaxy cluster. We use deep SINFONI data to perform a detailed study of the kinematics of the galaxies within 60 kpc of the radio core and we link this to the kinematics of the protocluster on the megaparsec scale. Identification of optical emission lines shows that 11 galaxies are at the redshift of the protocluster. The density of line emitters is more than an order of magnitude higher in the core of the protocluster than the larger scale environment. This implies a matter overdensity in the core of delta_m~70 which is similar to the outskirts of local galaxy clusters. The velocity distribution of the confirmed satellite galaxies shows a broad, double-peaked velocity structure with sigma=1360+/-206 km/s. A similar broad, double-peaked distribution was found in a previous study targeting the large scale protocluster structure, indicating that a common process is acting on both small and large scales. Including all spectroscopically confirmed protocluster galaxies, a velocity dispersion of 1013+/-87 km/s is found. We show that the protocluster has likely decoupled from the Hubble flow and is a dynamically evolved structure. Comparison to the Millenium simulation indicates that the protocluster velocity distribution is consistent with that of the most massive haloes at z~2, but we rule out that the protocluster is a fully virialized structure based on dynamical arguments and its X-ray luminosity. Comparison to merging haloes in the Millennium simulation shows that the structure as observed in and around the Spiderweb galaxy is best interpreted as being the result of a merger between two massive haloes. We propose that this merger can result in an increase in star formation and AGN activity in the protocluster core and is possibly an important stage in the evolution of massive cD galaxies.

This is a study of H-alpha emitters in two dense galaxy proto-clusters surrounding radio galaxies at z~2. We show that the proto-cluster surrounding MRC 1138-262 contains 14+/-2 times more H-alpha candidates than the average field (9 sigma significance), and the z=2.35 radio galaxy 4C+10.48 is surrounded by 12+/-2 times more emitters than the field (5 sigma), so it is also likely to reside in a dense proto-cluster environment. We compared these H-alpha emitters, situated in dense environments, to a control field sample selected from 3 separate fields forming a total area of 172 arcmin^2. We constructed and compared H-alpha and rest-frame R continuum luminosity functions of the emitters in both environments. The star formation density is on average 13 times greater in the proto-clusters than the field at z~2, so the total star formation rate within the central 1.5Mpc of the proto-clusters exceeds 3000Msun/yr. However, we found no significant difference in the shape of the H-alpha luminosity functions, implying that environment does not substantially affect the strength of the H-alpha line from strongly star forming galaxies. The proto-cluster emitters are typically 0.8mag brighter in rest-frame R continuum than field emitters, implying they are twice as massive as their field counterparts at the same redshift. We also show the proto-cluster galaxies have lower specific star formation rates than field galaxies, meaning the emitters in the dense environments formed more of their stars earlier than the field galaxies. We conclude that galaxy growth in the early Universe was accelerated in dense environments, and that cluster galaxies differed from field galaxies even before the cluster had fully formed.

We present the results of a deep, near-infrared, narrow band imaging survey at a central wavelength of 1.062 microns (FWHM=0.01 microns) in the GOODS-South field using the ESO VLT instrument, HAWK-I. The data are used to carry out the highest redshift search for [OII]3727 emission line galaxies to date. The images reach an emission line flux limit (5 sigma) of 1.5 x 10^-17 erg cm^-2 s^-1, additionally making the survey the deepest of its kind at high redshift. In this paper we identify a sample of [OII]3727 emission line objects at redshift z~1.85 in a co-moving volume of ~4100 Mpc^3. Objects are selected using an observed equivalent width (EW_obs) threshold of EW_obs = 50 angstroms. The sample is used to derive the space density and constrain the luminosity function of [OII] emitters at z=1.85. We find that the space density of objects with observed [OII] luminosities in the range log(L_[OII]) > 41.74 erg s^-1 is log(rho)=-2.45+/-0.14 Mpc^-3, a factor of 2 greater than the observed space density of [OII] emitters reported at z~1.4. After accounting for completeness and assuming an internal extinction correction of A_Halpha=1 mag (equivalent to A_[OII]=1.87), we report a star formation rate density of rho* ~0.38+/-0.06 Msun yr^-1 Mpc^-3. We independently derive the dust extinction of the sample using 24 micron fluxes and find a mean extinction of A_[OII]=0.98+/-0.11 magnitudes (A_Halpha=0.52). This is significantly lower than the A_Halpha=1 (A[OII]=1.86) mag value widely used in the literature. Finally we incorporate this improved extinction correction into the star formation rate density measurement and report rho*~0.24+/-0.06 Msun yr^-1 Mpc^-3.

High redshift quasars (HZQs) with redshifts of z >~ 6 are so rare that any photometrically-selected sample of sources with HZQ-like colours is likely to be dominated by Galactic stars and brown dwarfs scattered from the stellar locus. It is impractical to reobserve all such candidates, so an alternative approach was developed in which Bayesian model comparison techniques are used to calculate the probability that a candidate is a HZQ, P_q, by combining models of the quasar and star populations with the photometric measurements of the object. This method was motivated specifically by the large number of HZQ candidates identified by cross-matching the UKIRT Infrared Deep Sky Survey (UKIDSS) Large Area Survey (LAS) to the Sloan Digital Sky Survey (SDSS): in the ~1900 deg^2 covered by the LAS in the UKIDSS Seventh Data Release (DR7) there are ~10^3 real astronomical point-sources with the measured colours of the target quasars, of which only ~10 are expected to be HZQs. Applying Bayesian model comparison to the sample reveals that most sources with HZQ-like colours have P_q <~ 0.1 and can be confidently rejected without the need for any further observations. In the case of the UKIDSS DR7 LAS, there were just 88 candidates with P_q >= 0.1; these object were prioritized for reobservation by ranking according to P_q (and their likely redshift, which was also inferred from the photometric data). Most candidates were rejected after one or two (moderate depth) photometric measurements by recalculating P_q using the new data. That left seven confirmed HZQs, three of which were previously identified in the SDSS and four of which were new UKIDSS discoveries. The high efficiency of this Bayesian selection method suggests that it could usefully be extended to other HZQ surveys (e.g. searches by Pan-STARRS or VISTA) as well as to other searches for rare objects.

We report the discovery of a significant excess of candidate Halpha emitters (HAEs) in the field of the radio galaxy 4C 23.56 at z=2.483. Using the MOIRCS near-infrared imager on the Subaru Telescope we found 11 candidate emission-line galaxies to a flux limit of ~7.5 10^-17 erg s-1 cm-2, which is about 5 times excess from the expected field counts with ~3-sigma significance. Three of these are spectroscopically confirmed as redshifted Halpha at z=2.49. The distribution of candidate emitters on the sky is tightly confined to a 1.2-Mpc-radius area at z=2.49, locating 4C 23.56 at the western edge of the distribution. Analysis of the deep Spitzer MIPS 24 mu m imaging shows that there is also an excess of faint MIPS sources. All but two of the 11 HAEs are also found in the MIPS data. The inferred star-formation rate (SFR) of the HAEs based on the extinction-corrected Halpha luminosity (median SFR >~100 M_solar yr-1) is similar to those of HAEs in random fields at z~2. On the other hand, the MIPS-based SFR for the HAEs is on average 3.6 times larger, suggesting the existence of the star-formation significanly obscured by dust. The comparison of the Halpha-based star-formation activities of the HAEs in the 4C 23.56 field to those in another proto-cluster around PKS 1138-262 at z=2.16 reveals that the latter tend to have fainter Halpha emission despite similar K-band magnitudes. This suggests that star-formation may be suppressed in the PKS 1138-262 protocluster relative to the 4C 23.56 protocluster. This difference among the HAEs in the two proto-clusters at z > 2 may imply that some massive cluster galaxies are just forming at these epochs with some variation among clusters.

We study the environments of 6 radio galaxies at 2.2 < z < 2.6 using wide-field near-infrared images. We use colour cuts to identify galaxies in this redshift range, and find that three of the radio galaxies are surrounded by significant surface overdensities of such galaxies. The excess galaxies that comprise these overdensities are strongly clustered, suggesting they are physically associated. The colour distribution of the galaxies responsible for the overdensity are consistent with those of galaxies that lie within a narrow redshift range at z ~ 2.4. Thus the excess galaxies are consistent with being companions of the radio galaxies. The overdensities have estimated masses in excess of 10^14 solar masses, and are dense enough to collapse into virizalised structures by the present day: these structures may evolve into groups or clusters of galaxies. A flux-limited sample of protocluster galaxies with K < 20.6 mag is derived by statistically subtracting the fore- and background galaxies. The colour distribution of the protocluster galaxies is bimodal, consisting of a dominant blue sequence, comprising 77 +/- 10% of the galaxies, and a poorly populated red sequence. The blue protocluster galaxies have similar colours to local star-forming irregular galaxies (U -V ~ 0.6), suggesting most protocluster galaxies are still forming stars at the observed epoch. The blue colours and lack of a dominant protocluster red sequence implies that these cluster galaxies form the bulk of their stars at z < 3.

We present a study of protoclusters associated with high redshift radio galaxies. We imaged MRC1017-220 (z=1.77) and MRC0156-252 (z=2.02) using the near-infrared wide-field (7.5'x7.5') imager VLT/HAWK-I in the Y, H and Ks bands. We present the first deep Y-band galaxy number counts within a large area (200 arcmin2). We then develop a purely near-infrared colour selection technique to isolate galaxies at 1.6<z<3 that may be associated with the two targets, dividing them into (i) red passively evolving or dusty star-forming galaxies or (ii) blue/star-formation dominated galaxies with little or no dust. Both targeted fields show an excess of star-forming galaxies with respect to control fields. No clear overdensity of red galaxies is detected in the surroundings of MRC1017-220 although the spatial distribution of the red galaxies resembles a filament-like structure within which the radio galaxy is embedded. In contrast, a significant overdensity of red galaxies is detected in the field of MRC0156-252, ranging from a factor of 2-3 times the field density at large scales (2.5Mpc, angular distance) up to a factor of 3-4 times the field density within a 1Mpc radius of the radio galaxy. Half of these red galaxies have colours consistent with red sequence models at z~2, with a large fraction being bright (Ks<21.5, i.e. massive). In addition, we also find a small group of galaxies within 5" of MRC0156-252 suggesting that the radio galaxy has multiple companions within ~50 kpc. We conclude that the field of MRC0156-252 shows many remarkable similarities with the well-studied protocluster surrounding PKS1138-262 (z=2.16) suggesting that MRC0156-252 is associated with a galaxy protocluster at z~2.

We report on the environmental dependence of galaxy properties at z=2.15. We construct multi-band photometric data sets in the (proto-)cluster PKS1138-26 field and in the GOODS field. We then fit spectral energy distributions of the galaxies with model templates generated with the latest stellar population synthesis code and derive physical properties of galaxies from the fits. To quantify the environmental dependence of galaxy properties, a special care is taken of systematic errors -- we use data sets that have almost the same wavelength samplings, use the same code to fit SEDs with the same set of templates, and compare relative differences between the two samples. We find that the PKS1138 galaxies have similar ages, shorter star formation time scales, lower star formation rates, and weaker dust extinction compared to the GOODS galaxies at z~2. This trend is similar to that observed locally, suggesting that the environmental dependence of galaxy properties is already partly in place as early as z=2.15. We show that the PKS1138 galaxies assemble the bulk of their masses ~1 Gyr earlier than field galaxies, i.e., the galaxy formation depends on environment. Galaxy mergers should frequently occur during the first collapse of clusters and they might play an important role in driving the observed environmental dependence of galaxy properties at z=2.15.

We present a population study of several types of galaxies within the protocluster surrounding the radio galaxy MRC0316-257 at z~3.1. In addition to the known population of Ly_alpha emitters (LAEs) and [OIII] emitters, we use colour selection techniques to identify protocluster candidates that are Lyman break galaxies (LBG) and Balmer break galaxies (BBGs). The radio galaxy field contains an excess of LBG candidates, with a surface density 1.6\pm0.3 times larger than found for comparable blank fields. This surface overdensity corresponds to an LBG volume overdensity of ~8\pm4. The BBG photometric redshift distribution peaks at the protocluster's redshift, but we detect no significant surface overdensity of BBG. This is not surprising because a volume overdensity similar to the LBGs would have resulted in a surface density of ~1.2 that found in the blank field. This could not have been detected in our sample. Masses and star formation rates of the candidate protocluster galaxies are determined using SED fitting. These properties are not significantly different from those of field galaxies. The galaxies with the highest masses and star formation rates are located near the radio galaxy, indicating that the protocluster environment influences galaxy evolution at z~3. We conclude that the protocluster around MRC0316-257 is still in the early stages of formation.

We present a spectroscopic campaign to follow-up red colour-selected candidate massive galaxies in two high redshift proto-clusters surrounding radio galaxies. We observed a total of 57 galaxies in the field of MRC0943-242 (z=2.93) and 33 in the field of PKS1138-262 (z=2.16) with a mix of optical and near-infrared multi-object spectroscopy. We confirm two red galaxies in the field of PKS1138-262 at the redshift of the radio galaxy. Based on an analysis of their spectral energy distributions, and their derived star formation rates from the H-alpha and 24um flux, one object belongs to the class of dust-obscured star-forming red galaxies, while the other is evolved with little ongoing star formation. This result represents the first red and mainly passively evolving galaxy to be confirmed as companion galaxies in a z>2 proto-cluster. Both red galaxies in PKS1138-262 are massive, of the order of 4-6x10^11 M_Sol. They lie along a Colour-Magnitude relation which implies that they formed the bulk of their stellar population around z=4. In the MRC0943-242 field we find no red galaxies at the redshift of the radio galaxy but we do confirm the effectiveness of our JHK_s selection of galaxies at 2.3<z<3.1, finding that 10 out of 18 (56%) of JHK_s-selected galaxies whose redshifts could be measured fall within this redshift range. We also serendipitously identify an interesting foreground structure of 6 galaxies at z=2.6 in the field of MRC0943-242. This may be a proto-cluster itself, but complicates any interpretation of the red sequence build-up in MRC0943-242 until more redshifts can be measured.

Optical and near-infrared (NIR) spectra are presented for ULAS J131911.29+095951.4 (hereafter ULAS J1319+0950), a new redshift z = 6.127 +/- 0.004 quasar discovered in the Third Data Release (DR3) of the UKIRT Infrared Deep Sky Survey (UKIDSS). The source has Y_Vega = 19.10 +/- 0.03, corresponding to M_1450,AB = -27.12, which is comparable to the absolute magnitudes of the z ~= 6 quasars discovered in the SDSS. ULAS J1319+0950 was, in fact, registered by SDSS as a faint source with z_AB = 20.1 +/- 0.1, just below the signal--to--noise ratio limit of the SDSS high-redshift quasar survey. The faint z-band magnitude is a consequence of the weak Ly alpha/N V emission line, which has a rest-frame equivalent width of only ~20 A and provides only a small boost to the z-band flux. Nevertheless, there is no evidence from this UKIDSS-based search for a significant population of high-redshift quasars with weak emission lines. The Ly alpha optical depth to ULAS J1319+0950 is consistent with that measured towards similarly distant SDSS quasars, implying that results from optical- and NIR-selected quasars may be combined. Also presented is a new NIR-spectrum of the previously discovered UKIDSS quasar ULAS J020332.38+001229.2, which reveals the object to be a broad absorption line quasar. The new spectrum shows that the emission line previously identified as Ly alpha is actually N V, leading to a revised redshift of z = 5.72, rather than z = 5.86 as previously estimated.

We present new information on galaxies in the vicinity of luminous radio galaxies and quasars at z=4,5,6. These fields were previously found to contain overdensities of Lyman Break Galaxies (LBGs) or spectroscopic Lyman alpha emitters. We use HST and Spitzer data to infer stellar masses, and contrast our results with large samples of LBGs in more average environments as probed by the Great Observatories Origins Deep Survey (GOODS). The following results were obtained. First, LBGs in both overdense regions and in the field at z=4-5 lie on a very similar sequence in a z'-[3.6] versus [3.6] color-magnitude diagram. This is interpreted as a sequence in stellar mass (log[M*/Msun] = 9-11) in which galaxies become increasingly red due to dust and age as their star formation rate (SFR) increases. Second, the two radio galaxies are among the most massive objects (log[M*/Msun]~11) known to exist at z~4-5, and are extremely rare based on the low number density of such objects as estimated from the ~25x larger area GOODS survey. We suggest that the presence of these massive galaxies and supermassive black holes has been boosted through rapid accretion of gas or merging inside overdense regions. Third, the total stellar mass found in the z=4 ``proto-cluster'' TN1338 accounts for <30% of the stellar mass on the cluster red sequence expected to have formed at z>4, based on a comparison with the massive X-ray cluster Cl1252 at z=1.2. Although future near-infrared observations should determine whether any massive galaxies are currently being missed, one possible explanation for this mass difference is that TN1338 evolves into a smaller cluster than Cl1252. This raises the interesting question of whether the most massive protocluster regions at z>4 remain yet to be discovered.

We present the discovery of two brown dwarfs in the UKIRT Infrared Deep Sky Survey (UKIDSS) Deep Extragalactic Survey (DXS) Data Release 2. Both objects were selected photometrically from six square degrees in DXS for their blue J-K colour and the lack of optical counterparts in the Sloan Digital Sky Survey (SDSS) Stripe 82. Additional optical photometry provided by the Canada-France-Hawaii Telescope Legacy Survey (CFHT-LS) corroborated the possible substellarity of these candidates. Subsequent methane imaging of UDXS J221611.51+003308.1 and UDXS J221903.10+002418.2, has confirmed them as T7$\pm$1 and T6$\pm$1 dwarfs at photometric distances of 81 (52-118 pc) and 60 (44-87 pc; 2 sigma confidence level). A similar search in the second data release of the Ultra Deep Survey over a smaller area (0.77 square degree) and shallower depth didn't return any late-T dwarf candidate. The numbers of late-T dwarfs in our study are broadly in line with a declining mass function when considering the current area and depth of the DXS and UDS. These brown dwarfs are the first discovered in the VIMOS 4 field and among the few T dwarfs found in pencil-beam surveys. They are valuable to investigate the scale height of T dwarfs.

The UKIRT Infrared Deep Sky Survey (UKIDSS) Large Area Survey (LAS) has the necessary combination of filters (Y, J, H and K), depth (Y <~ 20.2) and area coverage (~4000 deg^2) to detect several redshift z >~ 6.4 quasars. The Third Data Release (DR3) included ~1000 deg^2 of LAS observations which have so far yielded two previously known z ~= 6 quasars and two new discoveries: ULAS J0203+0012, at z = 5.72; and ULAS J1319+0950, at z = 6.13.

We report the discovery of three very late T dwarfs in the UKIRT Infrared Deep Sky Survey (UKIDSS) Third Data Release: ULAS J101721.40+011817.9 (ULAS1017), ULAS J123828.51+095351.3 (ULAS1238) and ULAS J133553.45+113005.2 (ULAS1335).We detail optical and near-infrared photometry for all three sources, and mid-infrared photometry for ULAS1335. We use near-infrared spectra of each source to assign spectral types T8p (ULAS1017), T8.5 (ULAS1228) and T9 (ULAS1335) to these objects. We estimate that ULAS1017 has 750 < Teff < 850K, and 5.0 < log g < 5.5, assuming solar metallicity, an age of 1.6-15 Gyr, a mass of 33-70 MJ and lies at a distance of 31-54 pc. We extend the unified scheme of Burgasser et al. (2006) to the type T9 and suggest the inclusion of the WJ index to replace the now saturated J-band indices. ULAS1335 is the same spectral type as ULAS J003402.77-005206.7 and CFBDS J005910.90-011401.3. Comparison of model spectra with that of ULAS1335 suggest a temperature below 600K. We find ULAS1335 to be extremely red in near to mid-infrared colours, with H-[4.49]=4.34+/-0.04. This is the reddest near to mid-infrared colour yet observed for a T dwarf, which supports Teff < 600K, and we estimate Teff ~550-600K for ULAS1335. We estimate that ULAS1335 has an age of 0.6-5.3 Gyr, a mass of 15-31 MJ and lies at a distance of 8-12 pc.

Venemans et al. (2005) found evidence for an overdensity of Ly-alpha emission line galaxies associated with the radio galaxy MRC 0316-257 at z=3.13 indicating the presence of a massive protocluster. Here, we present the results of a search for additional star-forming galaxies and AGN within the protocluster. Narrow-band infrared imaging was used to select candidate [O III] emitters in a 1.1 x 1.1 Mpc^2 region around the radio galaxy. Thirteen candidates have been detected. Four of these are among the previously confirmed sample of Ly-alpha galaxies, and an additional three have been confirmed through follow-up infrared spectroscopy. The three newly confirmed objects lie within a few hundred km/s of each other, but are blueshifted with respect to the radio galaxy and Ly-alpha emitters by ~2100 km/s. Although the sample is currently small, our results indicate that the radio-selected protocluster is forming at the centre of a larger, ~60 co-moving Mpc super-structure. On the basis of an HST/ACS imaging study we calculate dust-corrected star-formation rates and investigate morphologies and sizes of the [O III] candidate emitters. From a comparison of the star formation rate derived from UV-continuum and [O III] emission, we conclude that at least two of the [O III] galaxies harbour an AGN which ionized the O+ gas.

We present the discovery of fifteen new T2.5-T7.5 dwarfs (with estimated distances between ~24-93pc, identified in the first three main data releases of the UKIRT Infrared Deep Sky Survey. This brings the total number of T dwarfs discovered in the Large Area Survey (to date) to 28. These discoveries are confirmed by near infrared spectroscopy, from which we derive spectral types on the unified scheme of Burgasser et al. (2006). Seven of the new T dwarfs have spectral types of T2.5-T4.5, five have spectral types of T5-T5.5, one is a T6.5p, and two are T7-7.5. We assess spectral morphology and colours to identify T dwarfs in our sample that may have non-typical physical properties (by comparison to solar neighbourhood populations). The colours of the full sample of LAS T dwarfs show a possible trend to bluer Y-J with decreasing effective temperature beyond T8. By accounting for the main sources of incompleteness (selection, follow-up and spatial) as well as the effects of unresolved binarity and Malmquist bias, we estimate that there are 17+-4 >=T4 dwarfs in the J<=19 volume of the LAS second data release. Comparing this to theoretical predictions is most consistent with a sub-stellar mass function exponent alpha between -1.0 and 0. This is consistent with the latest 2MASS/SDSS constraint (which is based on lower number statistics), and is significantly lower than the alpha~1.0 suggested by L dwarf field populations, possibly a result of the lower mass range probed by the T dwarf class.

We present new constraints on the evolution of the early-type galaxy color-magnitude relation (CMR) based on deep near-infrared imaging of a galaxy protocluster at z=2.16 obtained using NICMOS on-board the Hubble Space Telescope. This field contains a spectroscopically confirmed space-overdensity of Lyman-alpha and H-alpha emitting galaxies which surrounds the powerful radio galaxy MRC 1138-262. Using these NICMOS data we identify a significant surface-overdensity (= 6.2x) of red J-H galaxies in the color-magnitude diagram (when compared with deep NICMOS imaging from the HDF-N and UDF). The optical-NIR colors of these prospective red-sequence galaxies indicate the presence of on-going dust-obscured star-formation or recently formed (<~ 1.5 Gyr)stellar populations in a majority of the red galaxies. We measure the slope and intrinsic scatter of the CMR for three different red galaxy samples selected by a wide color cut, and using photometric redshifts both with and without restrictions on rest-frame optical morphology. In all three cases both the rest-frame $U-B$ slope and intrinsic color scatter are considerably higher than corresponding values for lower redshift galaxy clusters. These results suggest that while some relatively quiescent galaxies do exist in this protocluster both the majority of the galaxy population and hence the color-magnitude relation are still in the process of forming, as expected.

We present deep HST/ACS observations in g,r,i,z towards the z=4.1 radio galaxy TN J1338-1942 and its overdensity of >30 spectroscopically confirmed Lya emitters (LAEs). We select 66 g-band dropouts to z=27, 6 of which are also a LAE. Although our color-color selection results in a relatively broad redshift range centered on z=4.1, the field of TN J1338-1942 is richer than the average field at the >5 sigma significance, based on a comparison with GOODS. The angular distribution is filamentary with about half of the objects clustered near the radio galaxy, and a small, excess signal (2 sigma) in the projected pair counts at separations of <10" is interpreted as being due to physical pairs. The LAEs are young (a few x 10^7 yr), small (<r_50> = 0.13") galaxies, and we derive a mean stellar mass of ~10^8-9 Msun based on a stacked K-band image. We determine star formation rates, sizes, morphologies, and color-magnitude relations of the g-dropouts and find no evidence for a difference between galaxies near TN J1338-1942 and in the field. We conclude that environmental trends as observed in clusters at much lower redshift are either not yet present, or are washed out by the relatively broad selection in redshift. The large galaxy overdensity, its corresponding mass overdensity and the sub-clustering at the approximate redshift of TN J1338-1942 suggest the assemblage of a >10^14 Msun structure, confirming that it is possible to find and study cluster progenitors in the linear regime at z>4.

(abridged) We report the discovery of a very cool brown dwarf, ULAS J003402.77-005206.7 (ULAS J0034-00), identified in UKIDSS DR1. We provide optical, near-infrared, and mid-infrared photometry of the source, and two near-infrared spectra. Comparing the spectral energy distribution of ULAS J0034-00 to that of the T8 brown dwarf 2MASS J0415-09, the latest-type and coolest well-studied brown dwarf to date, with Teff~750 K, we find evidence that ULAS J0034-00 is significantly cooler. First, the measured values of the near-infrared absorption spectral indices imply a later classification, of T8.5. Second, the H-[4.49] colour provides an empirical estimate of the temperature of 540<Teff<660 K (+/-2sig range). Third, the J- and H-band peaks are somewhat narrower in ULAS J0034-00, and detailed comparison against spectral models calibrated to 2MASS J0415-09 yields an estimated temperature lower by 60-120 K relative to 2MASS J0415-09 i.e. 630<Teff<690 K (+/-2sig), and lower gravity or higher metallicity according to the degenerate combination -0.5<delta(log g-2[m/H])<-0.25 (+/-2sig). Combining these estimates, and considering systematics, it is likely the temperature lies in the range 600<Teff<700 K. Despite the low inferred Teff we find no evidence for strong absorption by NH3 over the wavelength range 1.51-1.56 um. Evolutionary models imply that the mass and age are in the ranges 15-36 M(Jup) and 0.5-8 Gyr, respectively. The measured proper motion, of (0.37+/-0.07)arcsec/yr, combined with the photometrically estimated distance of 14-22 pc, implies a tangential velocity of ~30 km/s. ULAS J0034-00 is significantly bluer than 2MASS J0415-09 in Y-J, so future searches should allow for the possibility that cooler T dwarfs are bluer still.

We explore the evolved galaxy population in the proto-clusters around four high-z radio galaxies at 2<~z<~3 based on wide-field near-infrared imaging. Three of the four fields are known proto-clusters as demonstrated by overdensities of line emitting galaxies at the same redshifts as the radio galaxies. We imaged the fields of three targets (PKS1138-262, USS0943-242 and MRC0316-257) to a depth of Ks~22 (5sigma) over a 4'x7' area centered on the radio galaxies with a new wide-field NIR camera, MOIRCS, on the Subaru Telescope. Another target (USS1558-003) was observed with SOFI on the NTT to a depth of Ks=20.5 over a 5'x5' area. We apply colour cuts in J-Ks and/or JHKs in order to exclusively search for galaxies located at high redshifts: z>2. To the 5sigma limiting magnitudes, we see a significant excess of NIR selected galaxies by a factor of two to three compared to those found in the field of GOODS-South. The spatial distribution of these NIR selected galaxies is not uniform and traces structures similar to those of emission line galaxies, although the samples of NIR selected galaxies and emitters show little overlap. We focus on the NIR colour-magnitude sequence of the evolved population and find that the bright-end (M_{stars}>10^{11}Msun) of the red sequence is well populated by z~2 but much less so in the z~3 proto-clusters. This may imply that the bright-end of the colour-magnitude sequence first appeared between z=3 and 2, an era coinciding with the appearance of submm galaxies and the peak of the cosmic star formation rate. Our observations show that during the same epoch, massive galaxies are forming in high density environments by vigorous star formation and assembly.

We present the results of an optical and near-IR spectroscopic study of giant nebular emission line halos associated with three z > 3 radio galaxies, 4C 41.17, 4C 60.07 and B2 0902+34. Previous deep narrow band Ly-alpha imaging had revealed complex morphologies with sizes up to 100 kpc), possibly connected to outflows and AGN feedback from the central regions. The outer regions of these halos show quiet kinematics with typical velocity dispersions of a few hundred km/s, and velocity shears that can mostly be interpreted as being due to rotation. The inner regions show shocked cocoons of gas closely associated with the radio lobes. These display disturbed kinematics and have expansion velocities and/or velocity dispersions >1000 km/s. The core region is chemically evolved, and we also find spectroscopic evidence for the ejection of enriched material in 4C 41.17 up to a distance of approximately 60 kpc along the radio-axis. The dynamical structures traced in the Ly-alpha line are, in most cases, closely echoed in the Carbon and Oxygen lines. This shows that the Ly-alpha line is produced in a highly clumped medium of small filling factor, and can therefore be used as a tracer of the dynamics of high-z radio galaxies (HzRGs). We conclude that these HzRGs are undergoing a final jet-induced phase of star formation with ejection of most of their interstellar medium before becoming "red and dead" Elliptical galaxies.

We present the initial results from our search for high redshift, z > 6, quasars using near infrared data from the UKIRT Infrared Deep Sky Survey (UKIDSS) Large Area Survey (LAS). Our analysis of 106 deg^2 of sky from Data Release 1 (DR1) has resulted in the discovery of ULAS J020332.38+001229.2, a luminous (J_AB=20.0, J_Vega=19.1, M_1450=-26.2) quasar at z=5.86. Following candidate selection from the combined IR and optical catalogue data and stacking of multiple epoch Sloan Digital Sky Survey (SDSS) data, we have obtained optical spectroscopy for the only two z > 6 quasar candidates. The VLT FORS2 spectrum of ULAS J020332.38+001229.2 shows broad Ly alpha + NV1240 emission at ~8350 A and an abrupt continuum break due to absorption by the Ly alpha forest. The quasar is not present in the SDSS DR5 catalogue and the continuum spectral index of alpha=-1.4 (F_nu nu^alpha) is redder than a composite of SDSS quasars at similar redshifts (alpha=-0.5). The discovery of one z~6 quasar in ~100 deg^2 in a complete sample within our selection criteria down to a median depth of Y_AB=20.4 (7 sigma) is consistent with existing SDSS results. We describe our survey methodology, including the use of optical data from the SDSS and the highly effective procedures developed to isolate the very small surface density of high-probability quasar candidates.

[Abridged] We present the results of a large program conducted with the Very Large Telescope and Keck telescope to search for forming clusters of galaxies near powerful radio galaxies at 2.0 < z < 5.2. We obtained narrow- and broad-band images of nine radio galaxies and their surroundings. The imaging was used to select candidate Lyman alpha emitting galaxies in ~3x3 Mpc^2 areas near the radio galaxies. A total of 337 candidate emitters were found with a rest-frame Lyman alpha equivalent width of EW_0 > 15 A and Sigma = EW_0/Delta EW_0 > 3. Follow-up spectroscopy confirmed 168 Lyman alpha emitters near eight radio galaxies. The success rate of our selection procedure is 91%. At least six of our eight fields are overdense in Lyman alpha emitters by a factor 3-5. Also, the emitters show significant clustering in velocity space. In the overdense fields, the width of the velocity distributions of the emitters is a factor 2-5 smaller than the width of the narrow-band filters. Taken together, we conclude that we have discovered six forming clusters of galaxies (protoclusters). We estimate that roughly 75% of powerful (L_2.7GHz > 10^33 erg/s/Hz/sr) high redshift radio galaxies reside in a protocluster, with a sizes of at least 1.75 Mpc. We estimate that the protoclusters have masses in the range 2-9 x 10^14 Msun and they are likely to be progenitors of present-day (massive) clusters of galaxies. For the first time, we have been able to estimate the velocity dispersion of cluster progenitors from z~5 to ~2. The velocity dispersion of the emitters increases with cosmic time, in agreement with the dark matter velocity dispersion in numerical simulations of forming massive clusters.

We present a deep image of the radio galaxy MRC 1138-262 taken with the Hubble Space Telescope (HST) at a redshift of z = 2.2. The galaxy is known to have properties of a cD galaxy progenitor and be surrounded by a 3 Mpc-sized structure, identified with a protocluster. The morphology shown on the new deep HST/ACS image is reminiscent of a spider's web. More than 10 individual clumpy features are observed, apparently star-forming satellite galaxies in the process of merging with the progenitor of a dominant cluster galaxy 11 Gyr ago. There is an extended emission component, implying that star formation was occurring over a 50 times 40 kpc region at a rate of more than 100 M_sun/yr. A striking feature of the newly named ``Spiderweb galaxy'' is the presence of several faint linear galaxies within the merging structure. The dense environments and fast galaxy motions at the centres of protoclusters may stimulate the formation of these structures, which dominate the faint resolved galaxy populations in the Hubble Ultra Deep Field. The new image provides a unique testbed for simulations of forming dominant cluster galaxies.

We present broad-band imaging with the Subaru Telescope of a 25'x25' field surrounding the radio galaxy TN J1338-1942 at redshift z=4.1. The field contains excesses of Lyman-alpha emitters (LAEs) and Lyman break galaxies (LBGs) identified with a protocluster surrounding the radio galaxy. Our new wide-field images provide information about the boundary of the protocluster and its surroundings. There are 874 candidate LBGs within our field, having redshifts in the range z=3.5-4.5. An examination of the brightest of these (with i'< 25.0) shows that the most prominent concentration coincides with the previously discovered protocluster. The diameter of this galaxy overdensity corresponds to ~2 Mpc at z=4, consistent with the previous estimation using LAEs. Several other concentrations of LBGs are observed in the field, some of which may well be physically connected with the z=4.1 protocluster. The observed structure in the smoothed LBG distribution can be explained as the projection of large-scale structure, within the redshift range z=3.5-4.5, comprising compact overdensities and prominent larger voids. If the 5-8 observed compact overdensities are associated with protoclusters, the observed protocluster volume density is ~5x10^-6 Mpc^-3, similar to the volume density of rich clusters in the local Universe.

We present results from a new method of exploring the distant Universe. We use 3-D spectroscopy to sample a large cosmological volume at a time when the Universe was less than 3 billion years old to investigate the evolution of star-formation activity. Within this study we also discovered a high redshift type-II quasar which would not have been identified with imaging studies alone. This highlights the crucial role that integral-field spectroscopy may play in surveying the distant Universe in the future.

We present HST/ACS observations of the most distant radio galaxy known, TN J0924-2201 at z=5.2. This radio galaxy has 6 spectroscopically confirmed Lya emitting companion galaxies, and appears to lie within an overdense region. The radio galaxy is marginally resolved in i_775 and z_850 showing continuum emission aligned with the radio axis, similar to what is observed for lower redshift radio galaxies. Both the half-light radius and the UV star formation rate are comparable to the typical values found for Lyman break galaxies at z~4-5. The Lya emitters are sub-L* galaxies, with deduced star formation rates of 1-10 Msun/yr. One of the Lya emitters is only detected in Lya. Based on the star formation rate of ~3 Msun/yr calculated from Lya, the lack of continuum emission could be explained if the galaxy is younger than ~2 Myr and is producing its first stars. Observations in V_606, i_775, and z_850 were used to identify additional Lyman break galaxies associated with this structure. In addition to the radio galaxy, there are 22 V-break (z~5) galaxies with z_850<26.5 (5sigma), two of which are also in the spectroscopic sample. We compare the surface density of 2/arcmin^2 to that of similarly selected V-dropouts extracted from GOODS and the UDF Parallel fields. We find evidence for an overdensity to very high confidence (>99%), based on a counts-in-cells analysis applied to the control field. The excess is suggestive of the V-break objects being associated with a forming cluster around the radio galaxy.

Observations of the radio galaxy MRC 0316-257 at z=3.13 and the surrounding field are presented. Using narrow- and broad-band imaging obtained with the VLT, 77 candidate Ly-alpha emitters with a rest-frame equivalent width of > 15 A were selected in a ~7'x7' field around the radio galaxy. Spectroscopy of 40 candidate emitters resulted in the discovery of 33 emission line galaxies of which 31 are Ly-alpha emitters with redshifts similar to that of the radio galaxy, while the remaining two galaxies turned out to be [OII] emitters. The Ly-alpha profiles have widths (FWHM) in the range of 120-800 km/s, with a median of 260 km/s. Where the signal-to-noise was large enough, the Ly-alpha profiles were found to be asymmetric, with apparent absorption troughs blueward of the profile peaks, indicative of absorption along the line of sight of an HI mass of at least 2x10^2 - 5x10^4 M_sun. The properties of the Ly-alpha galaxies (faint, blue and small) are consistent with young star forming galaxies which are still nearly dust free. The volume density of Ly-alpha emitting galaxies in the field around MRC 0316-257 is a factor of 3.3+0.5-0.4 larger compared with the density of field Ly-alpha emitters at that redshift. The velocity distribution of the spectroscopically confirmed emitters has a FWHM of 1510 km/s, which is substantially smaller than the width of the narrow-band filter (FWHM ~ 3500 km/s). The peak of the velocity distribution is located within 200 km/s of the redshift of the radio galaxy. We conclude that the confirmed Ly-alpha emitters are members of a protocluster of galaxies at z~3.13. The size of the protocluster is larger than 3.3x3.3 Mpc^2. The mass of this structure is estimated to be > 3-6x10^14 M_sun and could be the progenitor of a cluster of galaxies similar to e.g. the Virgo cluster. (Abridged)

We present deep optical imaging of the z=4.1 radio galaxy TN J1338--1942 obtained using the ACS on-board HST. The radio galaxy is known to reside within a large galaxy overdensity (both in physical extent and density contrast). There is good evidence that this `protocluster' region is the progenitor of a present-day rich galaxy cluster. TN J1338 is the dominant galaxy in the protocluster, in terms of size and luminosity and therefore seems destined to evolve into the brightest cluster galaxy. The high spatial-resolution ACS images reveal several kpc-scale features within and around the radio galaxy. The continuum light is aligned with the radio axis and is resolved into two clumps in the i-band and z-band bands. These components have luminosities ~10^9 L_sun and sizes of a few kpc. The estimated star-formation rate for the whole radio galaxy is ~200 M_sun/yr. A simple model in which the jet has triggered star-formation in these continuum knots is consistent with the available data. An unusual feature is seen in Lyman-alpha emission. A wedge-shaped extension emanates from the radio galaxy perpendicularly to the radio axis. This `wedge' naturally connects to the surrounding, asymmetric, large-scale (~100 kpc) Lyman-alpha halo. We posit that the wedge is a starburst-driven superwind, associated with the first major epoch of formation of the brightest cluster galaxy. The shock and wedge are examples of feedback processes due to both AGN and star-formation in the earliest stages of massive galaxy formation.

We have used VIsible MultiObject Spectrograph Integral-Field Unit (VIMOS-IFU) observations centred on a radio galaxy at z=2.9 to search for Ly-alpha emitters within a comoving volume of ~ 10^4 Mpc^3. We find 14 Ly-alpha emitters with flux > 1.4 x 10^-20 Wm^-2, yielding a comoving space density of 0.0018 +/- 0.0006 Mpc^-3. We fit a Schechter luminosity function which agrees well with previous studies both at similar redshift (z ~ 3.4) and higher redshift (z ~ 5.7). We therefore find no evidence for evolution in the properties of Ly-alpha emitters between 3 < z < 6, although our sample is small. By summing the star-formation rates of the individual Ly-alpha emitters we find a total cosmic star-formation rate density of rho_SFR = 6.7 +/- 0.5 x 10^-3 M_solar yr^-1 Mpc^-3. Integrating over the luminosity function for the combined Ly-alpha surveys at z ~ 3.4 and accounting for the difference in obscuration between the Ly-alpha line and the UV-continuum yields an estimate of rho_SFR ~ 2.2 x 10^-2 M_solar yr^-1 Mpc^-3, in line with previous multi-colour and narrow-band surveys of high-redshift star-forming galaxies. The detection of high-redshift emission-line galaxies in our volumetric search shows that the unique capabilities of wide-field integral-field spectroscopy are well suited in searching for high-redshift galaxies in a relatively unbiased manner.

We present the results of narrow-band and broad-band imaging with the Very Large Telescope of the field surrounding the radio galaxy TN J0924-2201 at z = 5.2. Fourteen candidate Ly alpha emitters with a rest-frame equivalent width of > 20 A were detected. Spectroscopy of 8 of these objects showed that 6 have redshifts similar to that of the radio galaxy. The density of emitters at the redshift of the radio galaxy is estimated to be a factor 1.5-6.2 higher than in the field, and comparable to the density of Ly alpha emitters in radio galaxy protoclusters at z = 4.1, 3.1 and 2.2. The Ly alpha emitters near TN J0924-2201 could therefore be part of a structure that will evolve into a massive cluster. These observations confirm that substantial clustering of Ly alpha emitters occurs at z > 5 and support the idea that radio galaxies pinpoint high density regions in the early Universe.

The clustering properties of clusters, galaxies and AGN as a function of redshift are briefly discussed. It appears that extremely red objects at z ~ 1, and objects with J-K > 1.7 and photometric redshifts 2 < z_phot < 4 are highly clustered, indicating that a majority of these objects constitutes the progenitors of nearby ellipticals. Similarly clustered seem luminous radio galaxies at z\sim 1, indicating that these objects comprise a short lived phase in the lifetime of these red objects. The high level of clustering furthermore suggests that distant powerful radio galaxies (e.g. z>2) might be residing in the progenitors of nearby clusters -- proto-clusters. A number of observational projects targetting fields with distant radio galaxies, including studies of Lya and Ha emitters, Lyman break galaxies and (sub)mm and X-ray emitters, all confirm that such radio galaxies are located in such proto-clusters. Estimates of the total mass of the proto-clusters are similar to the masses of local clusters. If the total star formation rate which we estimate for the entire proto-clusters is sustained up to z~1, the metals in the hot cluster gas of local clusters can easily be accounted for.

We present a 1.2 mm (250 GHz) map obtained with MAMBO on the IRAM 30m telescope of the central 25 arcmin^2 of the proto-cluster surrounding the z=4.1 radio galaxy TN J1338-1942. The map reaches a 1sigma sensitivity of 0.6 mJy in the central area, increasing to 1.2 mJy at the edges. We detect 10 candidate mm sources, of which 8 are also detected in a deep VLA 1.4 GHz map and/or a VLT R-band image. Three sources have a flux density S_{1.2 mm}>4.0 mJy, representing a 7sigma overdensity compared to random field surveys, which predict only 1 such source in our map area. We obtained SCUBA/JCMT 850 um and 450 um photometry of six radio/optically identified MAMBO sources, confirming 5 of them with S/N>4. Radio-to-mm and mm-to-submm redshift estimators cannot put strong constraints on the redshifts of these MAMBO sources, but 9 of them are consistent within the uncertainties (mean Delta z=+2.6) with z=4.1. One faint MAMBO source is possibly identified with an extremely red object (R-K=6.1) at a likely spectroscopic redshift z=1.18. The four brightest MAMBO sources are all located north of the radio galaxy, while the densest area of companion Ly-alpha excess and Lyman break galaxies is to the southeast. None of the 14 spectroscopically confirmed Ly-alpha emitters in the MAMBO field are detected at 1.2 mm; their average 1.2 mm flux density is <S_{\rm 1.2mm}>=0.25+-0.24 mJy. If the mm sources lie at z=4.1, none of them show excess Ly-alpha emission in our narrow-band images. Both populations thus show no apparent overlap, possibly due to dust quenching the Ly-alpha emission. If the mm sources are part of the proto-cluster, our results suggest that galaxies with star formation rates of a few 1000 M_{Sun} yr^{-1} could be spread throughout the proto-cluster over projected scales of at least 2 Mpc.

Recently, we conducted a Very Large Telescope (VLT) large program to search for forming clusters by looking for overdensities of Ly-alpha emitters around high redshift radio galaxies. In total seven proto-clusters were discovered, including a proto-cluster around the radio galaxy MRC 0316--257 at z ~ 3.13. This structure has an excess of Ly-alpha emitters by a factor of 3 as compared to the field, and the derived mass is 2-5 x 10^14 M_sun. The Ly-alpha emitters in the proto-cluster are on average bluer than Lyman Break Galaxies (LBGs). Also, the galaxies are faint (sub L_*) and small (half light radii < 1.7 kpc, which is smaller than the average size of LBGs). This might indicate that, at least a fraction of, Ly-alpha emitters could be young (~ 10^6 yr), nearly dust-free, forming galaxies.

We have carried out narrow band imaging targeted at the Lya line of emitters associated with radio galaxies at 2 < z < 5. Subsequent spectroscopy has confirmed the identity of > 120 Lya emitters and led to the discovery of six proto-clusters up to z = 4.1 with velocity dispersions between 300 and 1000 km/s. The number of emission line galaxies in the observed fields is four to fifteen times higher than in blank fields. These results strongly support the idea that high redshift radio galaxies are the progenitors of central brightest cluster galaxies located in the progenitors of clusters of galaxies.

High redshift radio galaxies are great cosmological tools for pinpointing the most massive objects in the early Universe: massive forming galaxies, active super--massive black holes and proto--clusters. We report on deep narrow--band imaging and spectroscopic observations of several z > 2 radio galaxy fields to investigate the nature of giant Ly-alpha nebulae centered on the galaxies and to search for over--dense regions around them. We discuss the possible implications for our understanding of the formation and evolution of massive galaxies and galaxy clusters.

Imaging and spectroscopy with the Very Large Telescope have revealed 20 Lyman-alpha emitters within a projected distance of 1.3 Mpc and 600 km/s of the luminous radio galaxy TN J1338-1942 at z=4.1. Compared to the field density of Lyman-alpha emitters, this implies an overdensity on the order of 15. The structure has a projected size of at least 2.7 Mpc x 1.8 Mpc and a velocity dispersion of 325 km/s, which makes it the most distant structure known. Using the galaxy overdensity and assuming a bias parameter b = 3 - 5, the mass is estimated to be ~10^15 M_sun. The radio galaxy itself is surrounded by an uniquely asymmetric Lyman-alpha halo. Taken together with our previous data on PKS 1138-262 at z ~ 2.16, these results suggest that luminous radio sources are excellent tracers of high density regions in the early Universe, which evolve into present-day clusters. The statistics of bright radio sources and of concentrations in the Lyman break galaxy population are consistent with the picture that each of those concentrations harbours an active or passive luminous radio source.

We have recently discovered a forming cluster around the radio galaxy MRC 1138-262 at redshift 2.2. Besides the population of Ly alpha emitting galaxies that have been confirmed spectroscopically, we have detected many candidate H alpha emitters that seem to have a different spatial distribution from the other galaxies: they are more clustered towards the center of the cluster and seem to be distributed along the same direction as the radio source. We present here the characteristics of the Ly alpha and H alpha emitters and study the nature of these populations.

High redshift radio galaxies are amongst the most massive galaxies in the early Universe and have properties expected from central galaxies in forming clusters. We are carrying out an observational programme on the VLT to find and study galaxy proto clusters around radio galaxies at redshifts 2 < z < 4. First, we use narrow band imaging to select candidate galaxies which show excess Lyman alpha emission at redshifts similar to the central radio galaxy. Then, we use multi object spectroscopy to confirm the redshifts of these candidates and measure the velocity dispersion of the cluster members. Our goal is to observe a sample of about 10 targets and investigate galaxy overdensities as a function of redshift. Here, we report on the current progress of the programme and show some preliminary results which include the discovery of a structure of galaxies at redshift 4.1.

Recent observations of the environments of gamma-ray bursts (GRBs) favour massive stars as their progenitors, which are likely to be surrounded by gas and dust. The visibility of the optical and UV emission of a GRB are expected to depend on the characteristics of both the dust and the GRB emission itself. A reasonable distribution of surrounding dust is capable of absorbing all the optical and UV emission of the optical flash and afterglow of a GRB, unless the optical flash has a peak isotropic luminosity L_peak > 10^49 erg/s. This means that dark bursts should exist and these bursts will have to be studied at infrared rather than optical wavelengths. In this paper details will be given about the infrared GRB dust emission. The reprocessed dust emission peaks at a rest-frame wavelength of about 8 micron. Forthcoming space telescopes, in particular the IRAC camera aboard the Space Infrared Telescope Facility could detect this emission out to a redshift of about 2. However, an accurate position of the GRB afterglow must be provided for this emission to be identified, because the light curve of the reprocessed dust emission does not vary on time-scales less than several years.