We constrain the rest-frame FUV (1546\AA), NUV (2345\AA) and U-band (3690\AA) luminosity functions (LFs) and luminosity densities (LDs) with unprecedented precision from $z\sim0.2$ to $z\sim3$ (FUV, NUV) and $z\sim2$ (U-band). Our sample of over 4.3 million galaxies, selected from the CFHT Large Area $U$-band Deep Survey (CLAUDS) and HyperSuprime-Cam Subaru Strategic Program (HSC-SSP) data lets us probe the very faint regime (down to $M_\mathrm{FUV},M_\mathrm{NUV},M_\mathrm{U} \simeq -15$ at low redshift) while simultaneously detecting very rare galaxies at the bright end down to comoving densities $<10^{-5}$ Mpc$^{-3}$. Our FUV and NUV LFs are well fitted by single Schechter functions, with faint-end slopes that are very stable up to $z\sim2$. We confirm, but self-consistently and with much better precision than previous studies, that the LDs at all three wavelengths increase rapidly with lookback time to $z\sim1$, and then much more slowly at $1<z<2$--$3$. Evolution of the FUV and NUV LFs and LDs at $z<1$ is driven almost entirely by the fading of the characteristic magnitude, $M^\star_{UV}$, while at $z>1$ it is due to the evolution of both $M^\star_{UV}$ and the characteristic number density $\phi^\star_{UV}$. In contrast, the U-band LF has an excess of faint galaxies and is fitted with a double-Schechter form; $M^\star_\mathrm{U}$, both $\phi^\star_\mathrm{U}$ components, and the bright-end slope evolve throughout $0.2<z<2$, while the faint-end slope is constant over at least the measurable $0.05<z<0.6$. We present tables of our Schechter parameters and LD measurements that can be used for testing theoretical galaxy evolution models and forecasting future observations.

Identifying non-contaminated sample of high-redshift galaxies with escaping Lyman continuum (LyC) flux is important for understanding the sources and evolution of cosmic reionization. We present CLAUDS $u$-band photometry of the COSMOS field to probe LyC radiation from spectroscopically confirmed galaxies at $z\geq3.5$ and outside the standard Lyman-break galaxy colour selection expectations. Complementary to the CLAUDS data, we use Subaru multi-filter photometry, \textit {Hubble Space Telescope} ($HST$) multi-filter imaging, and the spectroscopic surveys D10K, VUDS and 3D-HST. We present a sample of Lyman continuum galaxy (LCG) candidates in the redshift range $3.5\lesssim z\lesssim5.1$. Here, we introduce 5 LCG candidates, where two are flagged quality 1 and three quality 2. The estimated $f_{\rm esc}^{\rm abs}$ for quality 1 candidates are in the range $\sim5\% - 73\%$ and $\sim30\% - 93\%$. These estimates are based on our derived parameters from individual galaxies as inputs to a range of BPASS models as well as mean intergalactic medium (IGM) and maximal intergalactic and circumgalactic media (IGM+CGM) transmission. We conclude that our search for LCGs is most likely biased to lines of sight with low HI densities or free from Lyman limit systems. Our two best LCG candidates have EW (Ly$\alpha)\leq50$\AA\ and we find no correlation or anti-correlation between EW (Ly$\alpha$), $f_{\rm esc}^{\rm abs}$, and $R_{\rm obs}$, the ratio of ionizing to non-ionizing observed flux in the measured passbands. Stacking candidates without solid LyC detections ($S/N<3$) results in an estimated $f_{\rm esc}^{\rm abs}$ from galaxies not greater than $1\%$.

We study the environments of a sample of 61 extremely rare z~1.6 Ultra-Massive Passively Evolving Galaxies (UMPEGs: stellar masses M_stars >10^11.5 M_sun) which -- based on clustering analysis presented in Cheema et al. (2020) -- appear to be associated with very massive (M_halo ~ 10^14.1 h^-1 M_sun) dark matter halos that are likely to be the progenitors of z~0 massive (Coma- and Virgo-like) galaxy clusters. We find that UMPEGs on average have fewer than one satellite galaxy with mass ratio M_sat : M_UMPEG >~ 1:5 (i.e., M_sat >~ 10^10.8 M_sun) within 0.5 Mpc; the large mass gap that we observe between the typical UMPEG and its most massive satellite implies that the z~1.6 UMPEGs assembled through major mergers. Using observed satellite counts with merger timescales from the literature, we estimate the growth rate due to mergers with mass ratio of >~ 1:4 to be ~13% Gyr^-1 (with a ~2x systematic uncertainty). This relatively low growth rate is unlikely to significantly affect the shape of the massive end of the stellar mass function, whose evolution must instead be driven by the quenching of new cohorts of ultra-massive star-forming galaxies. However, this growth rate is high enough that, if sustained to z~0, the typical z~1.6 M_UMPEG=10^11.6 M_sun UMPEG can grow into a M_stars~10^12 M_sun brightest cluster galaxy (BCG) of a present-day massive galaxy cluster. Our observations favour a scenario in which our UMPEGs are main-branch progenitors of some of the present-day BCGs that have first assembled through major mergers at high redshifts and grown further through (likely minor) merging at later times.

We use a 27.6 deg$^2$ survey to measure the clustering of $gzK_s$-selected quiescent galaxies at $z\sim1.6$, focusing on ultra-massive quiescent galaxies. We find that $z\sim1.6$ Ultra-Massive Passively Evolving Galaxies (UMPEGs), which have $K_s(AB)<19.75$ (stellar masses of $M_{stars}$ $>\sim 10^{11.4}M_{\odot}$ and mean $<$$M_{stars}$$>$ = $10^{11.5}M_{\odot}$), cluster more strongly than any other known galaxy population at high redshift. Comparing their correlation length, $r_0 = 29.77 \pm 2.75$ $ h^{-1}$Mpc, with the clustering of dark matter halos in the Millennium XXL N-body simulation suggests that these $z\sim1.6$ UMPEGs reside in dark matter halos of mass $M_{h}\sim10^{14.1}h^{-1}M_{\odot}$. Such very massive $z\sim1.6$ halos are associated with the ancestors of $z\sim0$ massive galaxy clusters such as the Virgo and Coma clusters. Given their extreme stellar masses and lack of companions with comparable mass, we surmise that these UMPEGs could be the already-quenched central massive galaxies of their (proto)clusters. We conclude that with only a modest amount of further growth in their stellar mass, $z\sim1.6$ UMPEGs could be the progenitors of some of the massive central galaxies of present-day massive galaxy clusters observed to be already very massive and quiescent near the peak epoch of the cosmic star formation.

The properties of low-mass galaxies hosting central black holes provide clues about the formation and evolution of the progenitors of supermassive black holes. In this letter, we present HSC-XD 52, a spectroscopically confirmed low-mass active galactic nucleus (AGN) at an intermediate redshift of $z\sim0.56$. We detect this object as a very luminous X-ray source coincident with a galaxy observed by the Hyper Suprime-Cam (HSC) as part of a broader search for low-mass AGN. We constrain its stellar mass through spectral energy distribution modeling to be LMC-like at $M_\star \approx 3 \times 10^9 M_\odot$, placing it in the dwarf regime. We estimate a central black hole mass of $M_\mathrm{BH} \sim 10^{6} M_\odot$. With an average X-ray luminosity of $L_X \approx 3.5 \times 10^{43}~\mathrm{erg}~\mathrm{s}^{-1}$, HSC-XD 52 is among the most luminous X-ray selected AGN in dwarf galaxies. The spectroscopic and photometric properties of HSC-XD 52 indicate that it is an intermediate redshift counterpart to local low-mass AGN.

The CFHT Large Area U-band Deep Survey (CLAUDS) uses data taken with the MegaCam mosaic imager on CFHT to produce images of 18.60 deg2 with median seeing of FWHM=0.92 arcsec and to a median depth of U = 27.1 AB (5 sigma in 2 arcsec apertures), with selected areas that total 1.36 deg2 reaching a median depth of U=27.7 AB. These are the deepest U-band images assembled to date over this large an area. These data are located in four fields also imaged to comparably faint levels in grizy and several narrowband filters as part of the Hyper Suprime-Cam (HSC) Subaru Strategic Program (HSC-SSP). These CFHT and Subaru datasets will remain unmatched in their combination of area and depth until the advent of the Large Synoptic Survey Telescope (LSST). This paper provides an overview of the scientific motivation for CLAUDS and gives details of the observing strategy, observations, data reduction, and data merging with the HSC-SSP. Three early applications of these deep data are used to illustrate the potential of the dataset: deep U-band galaxy number counts, z~3 Lyman break galaxy (LBG) selection, and photometric redshifts improved by adding CLAUDS U to the Subaru HSC grizy photometry.

We present James Clerk Maxwell Telescope Submillimetre Common-User Bolometer Array 2 (SCUBA-2) 850 & 450 $\mu$m observations ($\sigma_{850}\sim0.5$ mJy, $\sigma_{450}\sim5$ mJy) of the HS1549+19 and HS1700+64 survey fields containing two of the largest known galaxy over-densities at $z=2.85$ and $2.30$, respectively. We detect 56 sub-millimetre galaxies (SMGs) with SNR $> 4$ over $\sim 50$ arcmin$^2$ at 850 $\mu$m with flux densities of 3 - 17 mJy. The number counts indicate over-densities in the $3$-arcmin diameter core region ($\sim 1.5$ Mpc at $z=2.5$) of $6^{+4}_{-2}\times$ (HS1549) and $4^{+6}_{-2}\times$ (HS1700) compared to blank field surveys. Within these core regions, we spectroscopically confirm that approximately one third of the SMGs lie at the protocluster redshifts for both HS1549 and HS1700. We use statistical identifications of other SMGs in the wider fields to constrain an additional four candidate protocluster members in each system. We combine multi wavelength estimates of the star-formation rates (SFRs) from Lyman-break dropout- and narrowband-selected galaxies, and the SCUBA-2 SMGs, to estimate total SFRs of 12,$500\pm2800$ M$_\odot$ yr$^{-1}$ ($4900\pm1200$ M$_\odot$ yr$^{-1}$) in HS1549 (HS1700), and SFR densities (SFRDs) within the central 1.5-Mpc diameter of each protocluster to be $3000\pm900$ M$_\odot$ yr$^{-1}$ Mpc$^{-3}$ ($1300\pm400$ M$_\odot$ yr$^{-1}$ Mpc$^{-3}$) in the HS1549 (HS1700) protocluster, $\sim10^4\times$ larger than the global SFRDs found at their respective epochs, due to the concentration of star-forming galaxies in the small volume of the dense cluster cores. Our results suggest centrally concentrated starbursts within protoclusters may be a relatively common scenario for the build up of mass in rich clusters assembling at $z\gtrsim2$.

We introduce the largest to date survey of massive quiescent galaxies at redshift z~1.6. With these data, which cover 27.6 deg^2, we can find significant numbers of very rare objects such as ultra-massive quiescent galaxies that populate the extreme massive end of the galaxy mass function, or dense environments that are likely to become present-day massive galaxy clusters. In this paper, the first in a series, we apply our gzK adaptation of the BzK technique to select our z~1.6 galaxy catalog and then study the quiescent galaxy stellar mass function with good statistics over M_stars ~ 10^10.2 - 10^11.7 M_sun --- a factor of 30 in mass --- including 60 ultra-massive z~1.6 quiescent galaxies with M_stars > 10^11.5 M_sun. We find that the stellar mass function of quiescent galaxies at z~1.6 is well represented by the Schechter function over this large mass range. This suggests that the mass quenching mechanism observed at lower redshifts must have already been well established by this epoch, and that it is likely due to a single physical mechanism over a wide range of mass. This close adherence to the Schechter shape also suggests that neither merging nor gravitational lensing significantly affect the observed quenched population. Finally, comparing measurements of M* parameters for quiescent and star-forming populations (ours and from the literature), we find hints of an offset (M*_SF > M*_PE), that could suggest that the efficiency of the quenching process evolves with time.

We investigate the connection between environment and the different quenching channels that galaxies are prone to follow in the rest-frame NUVrK colour diagram, as identified by Moutard et al. (2016b). Namely, the fast quenching channel followed by $young$ low-mass galaxies and the slow quenching channel followed by $old$ high-mass ones. We make use of the >22 deg$^2$ covered the VIPERS Multi-Lambda Survey (VIPERS-MLS) to select a galaxy sample complete down to stellar masses of $M_* > 10^{9.4} M_\odot$ at $z < 0.65$ ($M_* > 10^{8.8} M_\odot$ at $z < 0.5$) and including 33,500 (43,000) quiescent galaxies properly selected at $0.2 < z < 0.65$, while being characterized by reliable photometric redshifts ($\sigma_{\delta z/(1+z)} \leq 0.04$) that we use to measure galaxy local densities. We find that (1) the quiescence of low-mass [$M_* \leq 10^{9.7} M_\odot$] galaxies requires a strong increase of the local density, which confirms the lead role played by environment in their fast quenching and, therefore, confirms that the low-mass upturn observed in the stellar mass function of quiescent galaxies is due to $environmental$ $quenching$. We also observe that (2) the reservoir of low-mass star-forming galaxies located in very dense regions (prone to environmental quenching) has grown between $z \sim 0.6$ and $ z \sim 0.4$ whilst the share of low-mass quiescent galaxies (expected to being environmentally quenched) may have simultaneously increased, which would plead for a rising importance of $environmental$ $quenching$ with cosmic time, compared to $mass$ $quenching$. We finally discuss the composite picture of such environmental quenching of low-mass galaxies and, in particular, how this picture may be consistent with a $delayed$-$then$-$rapid$ quenching scenario.

We present the first results of our pilot study of 8 photometrically selected Lyman continuum (LyC) emitting galaxy candidates from the COSMOS field and focus on their optical emission line ratios. Observations were performed in the H and K bands using the Multi-Object Spectrometer for Infra-Red Exploration (MOSFIRE) instrument at the Keck Observatory, targeting the [OII], H$\beta$, and [OIII] emission lines. We find that photometrically selected LyC emitting galaxy candidates have high ionization parameters, based on their high [OIII]/[OII] ratios (O32), with an average ratio for our sample of 2.5$\pm$0.2. Preliminary results of our companion Low Resolution Imaging Spectrometer (LRIS) observations, targeting LyC and Ly$\alpha$, show that those galaxies with the largest O32 are typically found to also be Ly$\alpha$ emitters. High O32 galaxies are also found to have tentative non-zero LyC escape fractions ($f_{esc}(LyC)$) based on $u$ band photometric detections. These results are consistent with samples of highly ionized galaxies, including confirmed LyC emitting galaxies from the literature. We also perform a detailed comparison between the observed emission line ratios and simulated line ratios from density bounded H$_{\textrm{II}}$ regions modeled using the photoionization code MAPPINGS V. Estimates of $f_{esc}(LyC)$ for our sample fall in the range from 0.0-0.23 and suggest possible tension with published correlations between O32 and $f_{esc}(LyC)$, adding weight to dichotomy of arguments in the literature. We highlight the possible effects of clumpy geometry and mergers that may account for such tension.

Interferometric observations have demonstrated that a significant fraction of single-dish submillimetre (submm) sources are blends of multiple submm galaxies (SMGs), but the nature of this multiplicity, i.e. whether the galaxies are physically associated or chance projections, has not been determined. We performed spectroscopy of 11 SMGs in six multi-component submm sources, obtaining spectroscopic redshifts for nine of them. For an additional two component SMGs, we detected continuum emission but no obvious features. We supplement our observed sources with four sources from the literature. This sample allows us to statistically constrain the physical nature of single-dish submm source multiplicity for the first time. In three [3/7, or 43 (-33/+39) per cent at 95% confidence] of the single-dish sources for which the nature of the blending is unambiguous, the components for which spectroscopic redshifts are available are physically associated, whereas 4/7 [57 (-39/+33) per cent] have at least one unassociated component. When components whose spectra exhibit continuum but no features and for which the photometric redshift is significantly different from the spectroscopic redshift of the other component are also considered, 6/9 [67 (-37/+26) per cent] of the single-dish sources are comprised of at least one unassociated component. The nature of the multiplicity of one single-dish source is ambiguous. We conclude that physically associated systems and chance projections both contribute to the multi-component single-dish submm source population. This result contradicts the conventional wisdom that bright submm sources are solely a result of merger-induced starbursts, as blending of unassociated galaxies is also important.