We present new estimates of the luminosity function (LF) and star formation rate density (SFRD) for an H$\alpha$ selected sample at $z\sim0.62$ from the Deep And Wide Narrow-band (DAWN) survey. Our results are based on a new H$\alpha$ sample in the extended COSMOS region (compared to Coughlin et al. 2018) with the inclusion of flanking fields, resulting in a total area coverage of $\sim$1.5 deg$^2$. A total of 241 H$\alpha$ emitters were selected based on robust selection criteria using spectro-photometric redshifts and broadband color-color classification. We explore the effect of different dust correction prescriptions by calculating the LF and SFRD using a constant dust extinction correction, A{$_{\textrm{H}\alpha}=1$} mag, a luminosity-dependent correction, and a stellar-mass dependent correction. The resulting H$\alpha$ LFs are well fitted using Schechter functions with best-fit parameters: L$^*=10^{42.24}$ erg s$^{-1}$, $\phi^*=10^{-2.85}$ Mpc$^{-3}$, $\alpha = -1.62$ for constant dust correction, L$^*=10^{42.31}$ erg s$^{-1}$, $\phi^*=10^{-2.8}$ Mpc$^{-3}$, $\alpha=-1.39$ for luminosity-dependent dust correction, and L$^*=10^{42.36}$ erg s$^{-1}$, $\phi^*=10^{-2.91}$ Mpc$^{-3}$, $\alpha = -1.48$, for stellar mass-dependent dust correction. The deep and wide nature of the DAWN survey effectively samples H$\alpha$ emitters over a wide range of luminosities, thereby providing better constraints on both the faint and bright end of the LF. Also, the SFRD estimates $\rho_{\textrm{SFR}}=10^{-1.39}$ M$_{\odot}$yr$^{-1}$Mpc$^{-3}$ (constant dust correction), $\rho_{\textrm{SFR}}=10^{-1.47}$ M$_{\odot}$yr$^{-1}$Mpc$^{-3}$ (luminosity-dependent dust correction), and $\rho_{\textrm{SFR}}=10^{-1.49}$ M$_{\odot}$yr$^{-1}$Mpc$^{-3}$ (stellar mass-dependent dust correction) are in good agreement with the evolution of SFRD across redshifts ($0 < z < 2$) seen from previous H$\alpha$ surveys.

We present our measurements of the H$\alpha$, [OIII], and [OII] luminosity functions as part of the Lyman Alpha Galaxies at Epoch of Reionization (LAGER) survey using our samples of 1577 $z = 0.47$ H$\alpha$-, 3933 $z = 0.93$ [OIII]-, and 5367 $z = 1.59$ [OII]-selected emission line galaxies in a single 3 deg$^2$ CTIO/Blanco DECam pointing of the COSMOS field. Our observations reach 5$\sigma$ depths of $8.2\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$ and comoving volumes of $(1-7)\times10^{5}$ Mpc$^3$ making our survey one of the deepest narrowband surveys. We measure the observed luminosity functions and find best-fits of $\phi^\star = 10^{-3.16\pm0.09}$ Mpc$^{-3}$ and $L^\star = 10^{41.72\pm0.09}$ erg s$^{-1}$ for H$\alpha$, $\phi^\star = 10^{-2.16^{+0.10}_{-0.12}}$ Mpc$^{-3}$ and $L^\star = 10^{41.38^{+0.07}_{-0.06}}$ erg s$^{-1}$ for [OIII], and $\phi^\star = 10^{-1.97^{+0.07}_{-0.07}}$ Mpc$^{-3}$ and $L^\star = 10^{41.66\pm0.03}$ erg s$^{-1}$ for [OII], with $\alpha$ fixed to $-1.75$, $-1.6$, and $-1.3$, respectively. An excess of bright $> 10^{42}$ erg s$^{-1}$ [OIII] emitters is observed and may be due to AGN contamination. Dust corrections are applied assuming $A_{\rm{H}\alpha} = 1$ mag. We also design our own empirical rest-frame $g - r$ calibration using SDSS DR12 data, test it against our $z = 0.47$ H$\alpha$ emitters with $z$COSMOS $1$D spectra, and calibrate it for $(g - r)$ between $-0.8$ and $1.3$ mag. Dust and AGN-corrected star formation rate densities (SFRDs) are measured as $\log_{10}\rho_{\rm{SFR}}/(\rm{M}_\odot\ \rm{yr}^{-1}\ \rm{Mpc}^{-3}) = -1.63\pm0.04$, $-1.07\pm0.06$, and $-0.90\pm0.10$ for H$\alpha$, [OIII], and [OII], respectively. We find our [OIII] and [OII] samples fully trace cosmic star formation activity at their respective redshifts in comparison to multi-wavelength SFRDs, while the H$\alpha$ sample traces $\sim 70$ percent of the total $z = 0.47$ SFRD.

We present the largest-ever sample of 79 Ly$\alpha$ emitters (LAEs) at $z\sim$ 7.0 selected in the COSMOS and CDFS fields of the LAGER project (the Lyman Alpha Galaxies in the Epoch of Reionization). Our newly amassed ultradeep narrowband exposure and deeper/wider broadband images have more than doubled the number of LAEs in COSMOS, and we have selected 30 LAEs in the second field CDFS. We detect two large-scale LAE-overdense regions in the COSMOS that are likely protoclusters at the highest redshift to date. We perform injection and recovery simulations to derive the sample incompleteness. We show significant incompleteness comes from blending with foreground sources, which however has not been corrected in LAE luminosity functions in {the} literature. The bright end bump in the Ly$\alpha$ luminosity function in COSMOS is confirmed with 6 (2 newly selected) luminous LAEs (L$_{Ly\alpha}$ $>$ 10$^{43.3}$ erg s$^{-1}$). Interestingly, the bump is absent in CDFS, in which only one luminous LAE is detected. Meanwhile, the faint end luminosity functions from the two fields well agree with each other. The 6 luminous LAEs in COSMOS coincide with 2 LAE-overdense regions, while such regions are not seen in CDFS. The bright-end luminosity function bump could be attributed to ionized bubbles in a patchy reionization. It appears associated with cosmic overdensities, thus supports an inside-out reionization topology at $z$ $\sim$ 7.0, i.e., the high density peaks were ionized earlier compared to the voids. An average neutral hydrogen fraction of $x_{HI}$ $\sim$ 0.2 -- 0.4 is derived at $z\sim$ 7.0 based on the cosmic evolution of the Ly$\alpha$ luminosity function.

Studying the internal processes of individual galaxies at kilo-parsec scales is crucial in enhancing our understanding of galaxy formation and evolution processes. In this work, we investigate the distribution of star formation rate (SFR), specific SFR (sSFR), and dust attenuation across individual galaxies for a sample of 32 galaxies selected from the MUSE-Wide Survey at 0.1 $< \textit{z} <$ 0.42 with a dynamic range in stellar masses between $10^{7.7}$ and $10^{10.3}$ M$_{\odot}$. We take advantage of the high spatial resolution of the MUSE integral-field spectrograph and measure reliable spatially resolved H$\alpha$ and H$\beta$ emission line maps for individual galaxies. We also derive resolved stellar mass, SFR and dust maps using pixel-by-pixel SED fitting on high resolution multi-band HST/ACS and HST/WFC3 data from the CANDELS survey. By combining these, we analyze the radial profile of various physical parameters across these galaxies. We observe a radial dependence in both stellar and nebular color excess profiles peaking at the inner regions of galaxies. We also find the color excess profiles to most strongly correlate with the integrated sSFRs of galaxies. The median sSFR$_{\mathrm{H}\alpha}$ radial profiles of galaxies in our sample show a 0.8 dex increase from the central regions outward. This increase compared to the almost flat median radial profile of sSFR$_{\mathrm{SED}}$, which traces longer timescales of star formation, is in favor of the inside-out quenching of star formation. We bring further evidence for this quenching scenario from the locus of different subregions of galaxies on the SFR-M$_{*}$ and sSFR-M$_{*}$ relations.

We present the design for the first narrowband filter NB964 for the Dark Energy Camera (DECam), which is operated on the 4m Blanco Telescope at the Cerro Tololo Inter-American Observatory. The NB964 filter profile is essentially defined by maximizing the power of searching for Lyman alpha emitting galaxies (LAEs) in the epoch of reionization, with the consideration of the night sky background in the near-infrared and the DECam quantum efficiency. The NB964 filter was manufactured by Materion in 2015. It has a central wavelength of 964.2 nm and a full width at half maximum (FWHM) of 9.2 nm. An NB964 survey named LAGER (Lyman Alpha Galaxies in the Epoch of Reionization) has been ongoing since December 2015. Here we report results of lab tests, on-site tests and observations with the NB964 filter. The excellent performances of this filter ensure that the LAGER project is able to detect LAEs at z~7 with a high efficiency.

We investigate the clustering and halo properties of $\sim 5000$ Ly$\alpha$-selected emission line galaxies (LAEs) from the Slicing COSMOS 4K (SC4K) and from archival NB497 imaging of SA22 split in 15 discrete redshift slices between $z \sim 2.5 - 6$. We measure clustering lengths of $r_0 \sim 3 - 6\ h^{-1}$ Mpc and typical halo masses of $\sim 10^{11}$ M$_\odot$ for our narrowband-selected LAEs with typical $L_{\rm{Ly}\alpha} \sim 10^{42 - 43}$ erg s$^{-1}$. The intermediate band-selected LAEs are observed to have $r_0 \sim 3.5 - 15\ h^{-1}$ Mpc with typical halo masses of $\sim 10^{11 - 12}$ M$_\odot$ and typical $L_{\rm{Ly}\alpha} \sim 10^{43 - 43.6}$ erg s$^{-1}$. We find a strong, redshift-independent correlation between halo mass and Ly$\alpha$ luminosity normalized by the characteristic Ly$\alpha$ luminosity, $L^\star(z)$. The faintest LAEs ($L \sim 0.1\ L^\star(z)$) typically identified by deep narrowband surveys are found in $10^{10}$ M$_\odot$ halos and the brightest LAEs ($L \sim 7\ L^\star(z)$) are found in $\sim 5 \times 10^{12}$ M$_\odot$ halos. A dependency on the rest-frame 1500 \AA~UV luminosity, M$_\rm{UV}$, is also observed where the halo masses increase from $10^{11}$ to $10^{13}$ M$_\odot$ for M$_\rm{UV} \sim -19$ to $-23.5$ mag. Halo mass is also observed to increase from $10^{9.8}$ to $10^{12.3}$ M$_\odot$ for dust-corrected UV star formation rates from $\sim 0.6$ to $10$ M$_\odot$ yr$^{-1}$ and continues to increase up to $10^{13.5}$ M$_\odot$ in halo mass, where the majority of those sources are AGN. All the trends we observe are found to be redshift-independent. Our results reveal that LAEs are the likely progenitors of a wide range of galaxies depending on their luminosity, from dwarf-like, to Milky Way-type, to bright cluster galaxies. LAEs therefore provide unique insight into the early formation and evolution of the galaxies we observe in the local Universe.

We present and explore deep narrow- and medium-band data obtained with the Subaru and the Isaac Newton telescopes in the ~2 deg$^2$ COSMOS field. We use these data as an extremely wide, low-resolution (R~20-80) IFU survey to slice through the COSMOS field and obtain a large sample of ~4000 Lyman-$\alpha$ (Lya) emitters from z~2 to z~6 in 16 redshift slices (SC4K). We present new Lya luminosity functions (LFs) covering a co-moving volume of ~10$^8$Mpc$^3$. SC4K extensively complements ultra-deep surveys, jointly covering over 4 dex in Lya luminosity and revealing a global (2.5<z<6) synergy LF with $\alpha=-1.93\pm0.12$, $\log_{10}\Phi^*=-3.45^{+0.22}_{-0.29}$ Mpc$^{-3}$ and $\log_{10}L^*=42.93^{+0.15}_{-0.11}$ erg/s. The Schechter component of the Lya LF reveals a factor ~5 rise in $L^*$ and a ~7x decline in $\Phi^*$ from z~2 to z~6. The data reveal an extra power-law (or Schechter) component above L~10$^{43.3}$ erg/s at z~2.2-3.5 and we show that it is partially driven by X-ray and radio AGN, as their Lya LF resembles the excess. The power-law component vanishes and/or is below our detection limits above z>3.5, likely linked with the evolution of the AGN population. The Lya luminosity density rises by a factor ~2 from z~2 to z~3 but is then found to be roughly constant (~$1.1\times10^{40}$ erg s$^{-1}$ Mpc$^{-3}$) to z~6, despite the ~0.7 dex drop in UV luminosity density. The Lya/UV luminosity density ratio rises from $4\pm1$% to $30\pm6$% from z~2.2 to z~6. Our results imply a rise of a factor of ~2 in the global ionisation efficiency ($\xi_{\rm ion}$) and a factor ~$4\pm1$ in the Lya escape fraction from z~2 to z~6, hinting for evolution in both the typical burstiness/stellar populations and even more so in the typical ISM conditions allowing Ly$\alpha$ photons to escape.

We present new results from near-infrared spectroscopy with Keck/MOSFIRE of [OIII]-selected galaxies at $z\sim3.2$. With our $H$ and $K$-band spectra, we investigate the interstellar medium (ISM) conditions, such as ionization states and gas metallicities. [OIII] emitters at $z\sim3.2$ show a typical gas metallicity of $\mathrm{12+log(O/H) = 8.07\pm0.07}$ at $\mathrm{log(M_*/M_\odot) \sim 9.0-9.2}$ and $\mathrm{12+log(O/H) = 8.31\pm0.04}$ at $\mathrm{log(M_*/M_\odot) \sim 9.7-10.2}$ when using the empirical calibration method. We compare the [OIII] emitters at $z\sim3.2$ with UV-selected galaxies and Ly$\alpha$ emitters at the same epoch and find that the [OIII]-based selection does not appear to show any systematic bias in the selection of star-forming galaxies. Moreover, comparing with star-forming galaxies at $z\sim2$ from literature, our samples show similar ionization parameters and gas metallicities as those obtained by the previous studies using the same calibration method. We find no strong redshift evolution in the ISM conditions between $z\sim3.2$ and $z\sim2$. Considering that the star formation rates at a fixed stellar mass also do not significantly change between the two epochs, our results support the idea that the stellar mass is the primary quantity to describe the evolutionary stages of individual galaxies at $z>2$.

We investigate the clustering properties of $\sim 7000$ H$\beta$+[OIII] and [OII] narrowband-selected emitters at $z \sim 0.8 - 4.7$ from the High-$z$ Emission Line Survey. We find clustering lengths, $r_0$, of $1.5 - 4.0h^{-1}$ Mpc and minimum dark matter halo masses of $10^{10.7 - 12.1}\rm{M}_\odot$ for our $z = 0.8 - 3.2$ H$\beta$+[OIII] emitters and $r_0 \sim 2.0 - 8.3h^{-1}$ Mpc and halo masses of $10^{11.5 - 12.6}\rm{M}_\odot$ for our $z = 1.5 - 4.7$ [OII] emitters. We find $r_0$ to strongly increase both with increasing line luminosity and redshift. By taking into account the evolution of the characteristic line luminosity, $L^\star(z)$, and using our model predictions of halo mass given $r_0$, we find a strong, redshift-independent increasing trend between $L/L^\star(z)$ and minimum halo mass. The faintest H$\beta$+[OIII] emitters are found to reside in $10^{9.5}\rm{M}_\odot$ halos and the brightest emitters in $10^{13.0}\rm{M}_\odot$ halos. For [OII] emitters, the faintest emitters are found in $10^{10.5} \rm{M}_\odot$ halos and the brightest emitters in $10^{12.6}\rm{M}_\odot$ halos. A redshift-independent stellar mass dependency is also observed where the halo mass increases from $10^{11}\rm{M}_\odot$ to $10^{12.5} \rm{M}_\odot$ for stellar masses of $10^{8.5}\rm{M}_\odot$ to $10^{11.5}\rm{M}_\odot$, respectively. We investigate the interdependencies of these trends by repeating our analysis in a $L_\textrm{line} - \rm{M}_\textrm{star}$ grid space for our most populated samples (H$\beta$+[OIII] $z = 0.84$ and [OII] $z = 1.47$) and find that the line luminosity dependency is stronger than the stellar mass dependency on halo mass. For $L > L^\star$ emitters at all epochs, we find a relatively flat trend with halo masses of $10^{12.5 - 13}\rm{M}_\odot$ which may be due to quenching mechanisms in massive halos which is consistent with a transitional halo mass predicted by models.

We present a multi-wavelength photometric catalog in the COSMOS field as part of the observations by the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). The catalog is based on Hubble Space Telescope Wide Field Camera 3 (HST/WFC3) and Advanced Camera for Surveys (ACS) observations of the COSMOS field (centered at RA: $10^h00^m28^s$, Dec:$+02^{\circ}12^{\prime}21^{\prime\prime}$). The final catalog has 38671 sources with photometric data in forty two bands from UV to the infrared ($\rm \sim 0.3-8\,\mu m$). This includes broad-band photometry from the HST, CFHT, Subaru, VISTA and Spitzer Space Telescope in the visible, near infrared and infrared bands along with intermediate and narrow-band photometry from Subaru and medium band data from Mayall NEWFIRM. Source detection was conducted in the WFC3 F160W band (at $\rm 1.6\,\mu m$) and photometry is generated using the Template FITting algorithm. We further present a catalog of the physical properties of sources as identified in the HST F160W band and measured from the multi-band photometry by fitting the observed spectral energy distributions of sources against templates.

We study the production rate of ionizing photons of a sample of 588 H$\alpha$ emitters (HAEs) and 160 Lyman-$\alpha$ emitters (LAEs) at $z=2.2$ in the COSMOS field in order to assess the implied emissivity from galaxies, based on their UV luminosity. By exploring the rest-frame Lyman Continuum (LyC) with GALEX/$NUV$ data, we find f$_{\rm esc} < 2.8\, (6.4)$% through median (mean) stacking. By combining the H$\alpha$ luminosity density with IGM emissivity measurements from absorption studies, we find a globally averaged $\langle$f$_{\rm esc}\rangle$ of $5.9^{+14.5}_{-4.2}$ % at $z=2.2$ if we assume HAEs are the only source of ionizing photons. We find similarly low values of the global $\langle$f$_{\rm esc}\rangle$ at $z\approx3-5$, also ruling out a high $\langle$f$_{\rm esc}\rangle$ at $z<5$. These low escape fractions allow us to measure $\xi_{ion}$, the number of produced ionizing photons per unit UV luminosity, and investigate how this depends on galaxy properties. We find a typical $\xi_{ion} \approx 10^{24.77\pm0.04}$ Hz erg$^{-1}$ for HAEs and $\xi_{ion} \approx 10^{25.14\pm0.09}$ Hz erg$^{-1}$ for LAEs. LAEs and low mass HAEs at $z=2.2$ show similar values of $\xi_{ion}$ as typically assumed in the reionization era, while the typical HAE is three times less ionizing. Due to an increasing $\xi_{ion}$ with increasing EW(H$\alpha$), $\xi_{ion}$ likely increases with redshift. This evolution alone is fully in line with the observed evolution of $\xi_{ion}$ between $z\approx2-5$, indicating a typical value of $\xi_{ion} \approx 10^{25.4}$ Hz erg$^{-1}$ in the reionization era.

We investigate the properties of $z$=2.23 H$\alpha$ and [OIII]$\lambda$5007 emitters using the narrow-band-selected samples obtained from the High-$z$ Emission Line Survey (HiZELS: Sobral et al. 2013). We construct two samples of the H$\alpha$ and [OIII] emitters and compare their integrated physical properties. We find that the distribution of stellar masses, dust extinction, star formation rates (SFRs), and specific SFRs, is not statistically different between the two samples. When we separate the full galaxy sample into three subsamples according to the detections of the H$\alpha$ and/or [OIII] emission lines, most of the sources detected with both H$\alpha$ and [OIII] show ${\rm log(sSFR_{UV})}$$\gtrsim$-9.5. The comparison of the three subsamples suggests that sources with strong [OIII] line emission tend to have the highest star-forming activity out all galaxies that we study. We argue that the [OIII] emission line can be used as a tracer of star-forming galaxies at high redshift, and that it is especially useful to investigate star-forming galaxies at $z$$>$3, for which H$\alpha$ emission is no longer observable from the ground.

We investigate the properties of $\sim7000$ narrow-band selected galaxies with strong H$\beta+$[OIII] and [OII] nebular emission lines from the High-$z$ Emission Line Survey (HiZELS) between $z \sim 0.8 - 5.0$. Our sample covers a wide range in stellar mass ($M\sim10^{7.5 - 12.0}$ M$_\odot$), rest-frame equivalent widths (EW$_\mathrm{rest}\sim 10 - 10^5$ \AA), and line luminosities ($L\sim10^{40.5 - 43.2}$ erg s$^{-1}$). We measure the H$\beta+$[OIII] and [OII]-selected stellar mass functions out to $z \sim 3.5$ and find that both $M_\star$ and $\phi_\star$ increases with cosmic time, which may be due to the [OIII] selection including an increasing fraction of AGN at lower redshifts. The [OII]-selected stellar mass functions show a constant $M_\star\sim10^{11.6}$ M$_\odot$ and a strong, increasing evolution with cosmic time in $\phi_\star$ in line with H$\alpha$ studies. We also investigate the EW$_\mathrm{rest}$ evolution as a function of redshift with a fixed mass range (10$^{9.5 - 10.0}$ M$_\odot$) and find an increasing trend best represented by $(1+z)^{3.81\pm0.14}$ and $(1+z)^{2.72\pm0.19}$ up to $z\sim 2$ and $\sim 3$ for H$\beta$+[OIII] and [OII] emitters, respectively. This is the first time that the EW evolution has been directly measured for H$\beta+$[OIII] and [OII] emitters up to these redshifts. There is evidence for a slower evolution for $z > 2$ in the H$\beta+$[OIII] EW and a decreasing trend for $z > 3$ in the [OII] EW evolution, which would imply low [OII] EW at the highest redshifts and higher [OIII]/[OII] line ratios. This suggests that the ionization parameter at higher redshift may be significantly higher than the local Universe. Our results set the stage for future near-IR space-based spectroscopic surveys to test our extrapolated predictions and also produce $z > 5$ measurements to constrain the high-$z$ end of the EW$_\mathrm{rest}$ and [OIII]/[OII] evolution.

We investigate the evolution of the H$\beta$+[OIII] and [OII] luminosity functions from $z \sim 0.8$ to $\sim5$ in four redshift slices per emission line using data from the High-{\it z} Emission Line Survey (HiZELS). This is the first time that the H$\beta$+[OIII] and [OII] luminosity functions have been studied at these redshifts in a self-consistent analysis. This is also the largest sample of [OII] and H$\beta$+[OIII] emitters (3475 and 3298 emitters, respectively) in this redshift range, with large co-moving volumes $\sim 1 \times 10^6$ Mpc$^{-3}$ in two independent volumes (COSMOS and UDS), greatly reducing the effects of cosmic variance. The emitters were selected by a combination of photometric redshift and color-color selections, as well as spectroscopic follow-up, including recent spectroscopic observations using DEIMOS and MOSFIRE on the Keck Telescopes and FMOS on Subaru. We find a strong increase in $L_\star$ and a decrease in $\phi_\star$ for both H$\beta$+[OIII] and [OII] emitters. We derive the [OII] star-formation history of the Universe since $z\sim5$ and find that the cosmic SFRD rises from $z \sim 5$ to $\sim 3$ and then drops towards $z \sim 0$. We also find that our star-formation history is able to reproduce the evolution of the stellar mass density up to $z\sim 5$ based only on a single tracer of star-formation. When comparing the H$\beta$+[OIII] SFRDs to the [OII] and H$\alpha$ SFRD measurements in the literature, we find that there is a remarkable agreement, suggesting that the H$\beta$+[OIII] sample is dominated by star-forming galaxies at high-$z$ rather than AGNs.

We present results from the largest contiguous narrow-band survey in the near-infrared. We have used WIRCam/CFHT and the lowOH2 filter (1.187$\pm$0.005 $\mu$m) to survey ~10 deg$^2$ of contiguous extragalactic sky in the SA22 field. A total of ~6000 candidate emission-line galaxies are found. We use deep ugrizJK data to obtain robust photometric redshifts. We combine our data with the High-redshift Emission Line Survey (HiZELS), explore spectroscopic surveys (VVDS, VIPERS) and obtain our own spectroscopic follow-up with KMOS, FMOS and MOSFIRE to derive large samples of high-redshift emission-line selected galaxies: 3471 H$\alpha$ emitters at z=0.8, 1343 [OIII]+H$\beta$ emitters at z=1.4 and 572 [OII] emitters at z=2.2. We probe co-moving volumes of >10$^6$ Mpc$^3$ and find significant over-densities, including an 8.5$\sigma$ (spectroscopically confirmed) over-density of H$\alpha$ emitters at z=0.81. We derive H$\alpha$, [OIII]+H$\beta$ and [OII] luminosity functions at z=0.8,1.4,2.2, respectively, and present implications for future surveys such as Euclid. Our uniquely large volumes/areas allow us to sub-divide the samples in thousands of randomised combinations of areas and provide a robust empirical measurement of sample/cosmic variance. We show that surveys for star-forming/emission-line galaxies at a depth similar to ours can only overcome cosmic-variance (errors <10%) if they are based on volumes >5x10$^{5}$ Mpc$^{3}$; errors on $L^*$ and $\phi^*$ due to sample (cosmic) variance on surveys probing ~10$^4$ Mpc$^{3}$ and ~10$^5$ Mpc$^{3}$ are typically very high: ~300% and ~40-60%, respectively.

We use spitzer-IRAC data to identify near-infrared counterparts to submillimeter galaxies detected with Herschel-SPIRE at 250um in the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). Using a likelihood ratio analysis we identify 146 reliable IRAC counterparts to 123 SPIRE sources out of the 159. We find that, compared to the field population, the SPIRE counterparts occupy a distinct region of 3.6 and 4.5um color-magnitude space, and we use this property to identify a further 23 counterparts to 13 SPIRE sources. The IRAC identification rate of 86% is significantly higher than those that have been demonstrated with wide-field ground-based optical and near-IR imaging of Herschel fields. We estimate a false identification rate of 3.6%, corresponding to 4 to 5 sources. Among the 73 counterparts that are undetected in SDSS, 57 have both 3.6 and 4.5um coverage. Of these 43 have [3.6] - [4.5]> 0 indicating that they are likely to be at z > 1.4. Thus, ~ 40% of identified SPIRE galaxies are likely to be high redshift (z > 1.4) sources. We discuss the statistical properties of the IRAC-identified SPIRE galaxy sample including far-IR luminosities, dust temperatures, star-formation rates, and stellar masses. The majority of our detected galaxies have 10^10 to 10^11 L_sun total IR luminosities and are not intense starbursting galaxies as those found at z ~ 2, but they have a factor of 2 to 3 above average specific star-formation rates compared to near-IR selected galaxy samples.

This paper describes the Hubble Space Telescope imaging data products and data reduction procedures for the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). This survey is designed to document the evolution of galaxies and black holes at $z\sim1.5-8$, and to study Type Ia SNe beyond $z>1.5$. Five premier multi-wavelength sky regions are selected, each with extensive multiwavelength observations. The primary CANDELS data consist of imaging obtained in the Wide Field Camera 3 / infrared channel (WFC3/IR) and UVIS channel, along with the Advanced Camera for Surveys (ACS). The CANDELS/Deep survey covers \sim125 square arcminutes within GOODS-N and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of \sim800 square arcminutes across GOODS and three additional fields (EGS, COSMOS, and UDS). We summarize the observational aspects of the survey as motivated by the scientific goals and present a detailed description of the data reduction procedures and products from the survey. Our data reduction methods utilize the most up to date calibration files and image combination procedures. We have paid special attention to correcting a range of instrumental effects, including CTE degradation for ACS, removal of electronic bias-striping present in ACS data after SM4, and persistence effects and other artifacts in WFC3/IR. For each field, we release mosaics for individual epochs and eventual mosaics containing data from all epochs combined, to facilitate photometric variability studies and the deepest possible photometry. A more detailed overview of the science goals and observational design of the survey are presented in a companion paper.

We present physical properties of two submillimeter selected gravitationally lensed sources, identified in the Herschel Astrophysical Terahertz Large Area Survey. These submillimeter galaxies (SMGs) have flux densities > 100 mJy at 500 um, but are not visible in existing optical imaging. We fit light profiles to each component of the lensing systems in Spitzer IRAC 3.6 and 4.5 um data and successfully disentangle the foreground lens from the background source in each case, providing important constraints on the spectral energy distributions (SEDs) of the background SMG at rest-frame optical-near-infrared wavelengths. The SED fits show that these two SMGs have high dust obscuration with Av ~4 to 5 and star formation rates of ~100 M_sun/yr. They have low gas fractions and low dynamical masses compared to 850 um selected galaxies.

The extragalactic background light at far-infrared wavelengths originates from optically-faint, dusty, star-forming galaxies in the universe with star-formation rates at the level of a few hundred solar masses per year. Due to the relatively poor spatial resolution of far-infrared telescopes, the faint sub-millimetre galaxies are challenging to study individually. Instead, their average properties can be studied using statistics such as the angular power spectrum of the background intensity variations. A previous attempt at measuring this power spectrum resulted in the suggestion that the clustering amplitude is below the level computed with a simple ansatz based on a halo model. Here we report a clear detection of the excess clustering over the linear prediction at arcminute angular scales in the power spectrum of brightness fluctuations at 250, 350, and 500 microns. From this excess, we find that sub-millimetre galaxies are located in dark matter halos with a minimum mass of log[M_min/M_sun ]= 11.5^+0.7_-0.2 at 350 microns. This minimum dark matter halo mass corresponds to the most efficient mass scale for star formation in the universe, and is lower than that predicted by semi-analytical models for galaxy formation.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

We measure the angular correlation function, w(theta), from 0.5 to 30 arcminutes of detected sources in two wide fields of the Herschel Multi-tiered Extragalactic Survey (HerMES). Our measurements are consistent with the expected clustering shape from a population of sources that trace the dark matter density field, including non-linear clustering at arcminute angular scales arising from multiple sources that occupy the same dark matter halos. By making use of the halo model to connect the spatial clustering of sources to the dark matter halo distribution, we estimate source bias and halo occupation number for dusty sub-mm galaxies at z ~ 2. We find that sub-mm galaxies with 250 micron flux densities above 30 mJy reside in dark matter halos with mass above (5\pm4) x 10^12 M_sun, while (14\pm8)% of such sources appear as satellites in more massive halos.