We study the evolution in the number density of the highest mass galaxies over $0.4<z<1.5$ (covering 9 Gyr). We use the Spitzer/HETDEX Exploratory Large-Area (SHELA) Survey, which covers 17.5 $\mathrm{deg}^2$ with eight photometric bands spanning 0.3-4.5 $\mu$m within the SDSS Stripe 82 field. This size produces the lowest counting uncertainties and cosmic variance yet for massive galaxies at $z\sim1.0$. We study the stellar mass function (SMF) for galaxies with $\log(M_\ast/M_\odot)>10.3$ using a forward-modeling method that fully accounts for statistical and systematic uncertainties on stellar mass. From $z$=0.4 to 1.5 the massive end of the SMF shows minimal evolution in its shape: the characteristic mass ($M^\ast$) evolves by less than 0.1 dex ($\pm$0.05 dex); the number density of galaxies with $\log (M_\ast/M_\odot) >11$ stays roughly constant at $\log (n/\mathrm{Mpc}^{-3})$ $\simeq$ $-$3.4 ($\pm$0.05), then declines to $\log (n/\mathrm{Mpc}^{-3})$=$-$3.7 ($\pm$0.05) at $z$=1.5. We discuss the uncertainties in the SMF, which are dominated by assumptions in the star formation history and details of stellar population synthesis models for stellar mass estimations. For quiescent galaxies, the data are consistent with no (or slight) evolution ($\lesssim0.1$ dex) in the characteristic mass nor number density from $z\sim 1.5$ to the present. This implies that any mass growth (presumably through "dry" mergers) of the quiescent massive galaxy population must balance the rate of mass losses from late-stage stellar evolution and the formation of quenching galaxies from the star-forming population. We provide a limit on this mass growth from $z=1.0$ to 0.4 of $\Delta M_\ast/M_\ast\leq$ 45% (i.e., $\simeq0.16$ dex) for quiescent galaxies more massive than $10^{11}$ $M_\odot$.

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.

We present the high-mass end of the galaxy stellar mass function using the largest sample to date (5,352) of star-forming galaxies with $M_{\star} > 10^{11} M_{\odot}$ at cosmic noon, $1.5 < z < 3.5$. This sample is uniformly selected across 17.2 deg$^2$ ($\sim$0.44 Gpc$^3$ comoving volume from $1.5 < z < 3.5$), mitigating the effects of cosmic variance and encompassing a wide range of environments. This area, a factor of 10 larger than previous studies, provides robust statistics at the high-mass end. Using multi-wavelength data in the Spitzer/HETDEX Exploratory Large Area (SHELA) footprint we find that the SHELA footprint star-forming galaxy stellar mass function is steeply declining at the high-mass end probing values as high as $\sim$$10^{-4}$ Mpc$^3$/dex and as low as $\sim$5$\times$$10^{-8}$ Mpc$^3$/dex across a stellar mass range of log($M_\star$/$M_\odot$) $\sim$ 11 - 12. We compare our empirical star-forming galaxy stellar mass function at the high mass end to three types of numerical models: hydrodynamical models from IllustrisTNG, abundance matching from the UniverseMachine, and three different semi-analytic models (SAMs; SAG, SAGE, GALACTICUS). At redshifts $1.5 < z < 3.5$ we find that results from IllustrisTNG and abundance matching models agree within a factor of $\sim$2 to 10, however the three SAMs strongly underestimate (up to a factor of 1,000) the number density of massive galaxies. We discuss the implications of these results for our understanding of galaxy evolution.

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.

The epoch of reionization (6 < z < 10) marks the period in our universe when the first large galaxies grew to fruition, and began to affect the universe around them. Massive stars, and potentially accreting supermassive black holes, filled the universe with ionizing radiation, burning off the haze of neutral gas that had filled the intergalactic medium (IGM) since recombination (z~1000). The evolution of this process constrains key properties of these earliest luminous sources, thus observationally constraining reionization is a key science goal for the next decade. The measurement of Lyman-alpha emission from photometrically-identified galaxies is a highly constraining probe of reionization, as a neutral IGM will resonantly scatter these photons, reducing detectability. While significant work has been done with 8-10m telescopes, these observations require extremely large telescopes (ELTs); the flux limits available from today's 10m class telescopes are sufficient for only the brightest known galaxies (m < 26). Ultra-deep surveys with the Giant Magellan Telescope (GMT) and Thirty Meter Telescope (TMT) will be capable of detecting Lyman-alpha emission from galaxies 2-3 magnitudes fainter than today's deepest surveys. Wide-field fiber spectroscopy on the GMT combined with narrow-field AO-assisted slit spectroscopy on the TMT will be able to probe the expected size of ionized bubbles throughout the epoch of reionization, following up degree scale deep imaging surveys with the Wide Field Infrared Space Telescope. These data will provide the first resolved Lyman-alpha-based maps of the ionized intergalactic medium throughout the epoch of reionization, constraining models of both the temporal and spatial evolution of this phase change.

Realizing the utility of Ly$\alpha$ emission to trace the evolution of the intergalactic medium (IGM) during the epoch of reionization requires deep spectroscopy across the boundary of optical and near-infrared (NIR) spectrographs at $z\sim7.2$ when Ly$\alpha$ emission is at $\sim$1$\mu$m. Our Texas Spectroscopic Search for Ly$\alpha$ Emission at the End of Reionization includes 18 nights of deep spectroscopic observations using the Keck DEIMOS (optical) and MOSFIRE (NIR) spectrographs. Within this dataset we observe Ly$\alpha$ emission from 183 photometric-redshift selected galaxies at $z =$ 5.5 - 8.3 from the Cosmic Assembly Near infrared Deep Extragalactic Legacy Survey (CANDELS). Our overlapping MOSFIRE observations, over 84 galaxies total, provide the deepest NIR spectroscopic data yet obtained for Ly$\alpha$ from galaxies $z > 7$, with $>16$ hr integration time for four observed galaxies. Here we analyze these four targets, and we report the discovery of a new $z = 7.60$ Ly$\alpha$ detection as well as provide an updated observation of the previously confirmed $z=7.51$ Ly$\alpha$ emission from Finkelstein et al. (2013) with a $\sim$3$\times$ longer exposure time. Our analysis of these Ly$\alpha$ emission line profiles reveal a significant asymmetric shape. The two detected Ly$\alpha$ emission lines from bright sources ($M_{\text{UV}}<-20.25$) could imply that these bright galaxies inhabit ionized bubbles in a partially neutral IGM, although deeper exposures may yet reveal Ly$\alpha$ emission in the fainter sources.

We present the $ugriz$-band Dark Energy Camera (DECam) plus 3.6 and 4.5 $\mu$m IRAC catalogs for the Spitzer/HETDEX Exploratory Large-Area (SHELA) survey. SHELA covers $\sim24$ deg$^{2}$ of the Sloan Digital Sky Survey (SDSS) Stripe 82 region, with seven bandpasses spanning a wavelength range of 0.35 to 4.5 $\mu$m. SHELA falls within the footprint of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX), which will provide spectroscopic redshifts for $\sim200{,}000$ Ly$\alpha$ emitters at $1.9<z<3.5$ and also for $\sim200{,}000$ [OII] emitters at $z<0.5$. SHELA's deep, wide-area multiwavelength images combined with HETDEX's spectroscopic information, will facilitate many extragalactic studies, including measuring the evolution of galaxy stellar mass, halo mass, and environment from $1.5<z<3.5$. Here we present $riz$-band selected $ugriz$-band DECam catalogs that reach a $5\sigma$ depth of $\sim24.5$ AB mag (for point sources with an aperture that encloses $70\%$ of the total flux) and cover $17.5$ deg$^{2}$ of the overall SHELA field. We validate our DECam catalog by comparison to the DECam Legacy Survey (DECaLS) DR5 and the Dark Energy Survey (DES) DR1. We perform IRAC forced photometry with The Tractor image modeling code to measure 3.6 and 4.5 $\mu$m fluxes for all objects within our DECam catalog. We demonstrate the utility of our catalog by computing galaxy number counts and estimating photometric redshifts. Our photometric redshifts recover the available $\left\langle z \right\rangle = 0.33 $ SDSS spectroscopic redshifts with a $1\sigma$ scatter in $\Delta z/(1 +z)$ of 0.04.

We present a joint analysis of the rest-frame ultraviolet (UV) luminosity functions of continuum-selected star-forming galaxies and galaxies dominated by active galactic nuclei (AGNs) at $z \sim$ 4. These 3,740 $z \sim$ 4 galaxies are selected from broad-band imaging in nine photometric bands over 18 deg$^2$ in the \textit{Spitzer}/HETDEX Exploratory Large Area Survey (SHELA) field. The large area and moderate depth of our survey provide a unique view of the intersection between the bright end of the galaxy UV luminosity function (M$_{AB}<-$22) and the faint end of the AGN UV luminosity function. We do not separate AGN-dominated galaxies from star-formation-dominated galaxies, but rather fit both luminosity functions simultaneously. These functions are best fit with a double power-law (DPL) for both the galaxy and AGN components, where the galaxy bright-end slope has a power-law index of $-3.80\pm0.10$, and the corresponding AGN faint-end slope is $\alpha_{AGN} = -1.49^{+0.30}_{-0.21}$. We cannot rule out a Schechter-like exponential decline for the galaxy UV luminosity function, and in this scenario the AGN luminosity function has a steeper faint-end slope of $-2.08^{+0.18}_{-0.11}$. Comparison of our galaxy luminosity function results with a representative cosmological model of galaxy formation suggests that the molecular gas depletion time must be shorter, implying that star formation is more efficient in bright galaxies at $z=4$ than at the present day. If the galaxy luminosity function does indeed have a power-law shape at the bright end, the implied ionizing emissivity from AGNs is not inconsistent with previous observations. However, if the underlying galaxy distribution is Schechter, it implies a significantly higher ionizing emissivity from AGNs at this epoch.

The distribution of Ly$\alpha$ emission is an presently accessible method for studying the state of the intergalactic medium (IGM) into the reionization era. We carried out deep spectroscopic observations in order to search for Ly$\alpha$ emission from galaxies with photometric redshifts $z$ = 5.5 - 8.3 selected from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). Utilizing data from the Keck/DEIMOS spectrograph, we explore a wavelength coverage of Ly$\alpha$ emission at $z$ ~ 5 - 7 with four nights of spectroscopic observations for 118 galaxies, detecting five emission lines with ~ 5$\sigma$ significance: three in the GOODS-N and two in the GOODS-S field. We constrain the equivalent width (EW) distribution of Ly$\alpha$ emission by comparing the number of detected objects with the expected number constructed from detailed simulations of mock emission lines that account for the observational conditions (e.g., exposure time, wavelength coverage, and sky emission) and galaxy photometric redshift probability distribution functions. The Ly$\alpha$ EW distribution is well described by an exponential form, $\text{dN/dEW}\propto \text{exp(-EW/}W_0)$, characterized by the $e$-folding scale ($W_0$) of ~ 60 - 100$\AA$ at 0.3 < $z$ < 6. By contrast, our measure of the Ly$\alpha$ EW distribution at 6.0 < $z$ < 7.0 rejects a Ly$\alpha$ EW distribution with $W_0$ > 36.4$\AA$ (125.3$\AA$) at 1$\sigma$ (2$\sigma$) significance. This provides additional evidence that the EW distribution of Ly$\alpha$ declines at $z$ > 6, suggesting an increasing fraction of neutral hydrogen in the IGM at that epoch.

We present a flux-limited sample of $z\sim0.3$ Ly$\alpha$ emitters (LAEs) from Galaxy Evolution Explorer (GALEX) grism spectroscopic data. The published GALEX $z\sim0.3$ LAE sample is pre-selected from continuum-bright objects and thus is biased against high equivalent width (EW) LAEs. We remove this continuum pre-selection and compute the EW distribution and the luminosity function of the Ly$\alpha$ emission line directly from our sample. We examine the evolution of these quantities from $z\sim0.3$ to $2.2$ and find that the EW distribution shows little evidence for evolution over this redshift range. As shown by previous studies, the Ly$\alpha$ luminosity density from star-forming galaxies declines rapidly with declining redshift. However, we find that the decline in Ly$\alpha$ luminosity density from $z=2.2$ to $z=0.3$ may simply mirror the decline seen in the H$\alpha$ luminosity density from $z=2.2$ to $z=0.4$, implying little change in the volumetric Ly$\alpha$ escape fraction. Finally, we show that the observed Ly$\alpha$ luminosity density from AGNs is comparable to the observed Ly$\alpha$ luminosity density from star-forming galaxies at $z=0.3$. We suggest that this significant contribution from AGNs to the total observed Ly$\alpha$ luminosity density persists out to $z\sim2.2$.

We report the detection of the most luminous high-redshift Lyman Alpha Emitting galaxy (LAE) yet seen, with log L(Ly alpha) = 43.9 ergs/s. The galaxy -- COSMOS Lyman alpha 1, or COLA1 -- was detected in a search for ultra-luminous LAEs with Hyper Suprime-Cam on the Subaru telescope. It was confirmed to lie at z = 6.593 based on a Lyman alpha line detection obtained from followup spectroscopy with the DEIMOS spectrograph on Keck2. COLA1 is the first very high-redshift LAE to show a multi-component Lyman alpha line profile with a blue wing, which suggests that it could lie in a highly ionized region of the intergalactic medium and could have significant infall. If this interpretation is correct, then ultra-luminous LAEs like COLA1 offer a unique opportunity to determine the properties of the HII regions around these galaxies which will help in understanding the ionization of the z ~ 7 intergalactic medium.

We construct a flux-limited sample of 135 candidate z~1 Lya emitters (LAEs) from Galaxy Evolution Explorer (GALEX) grism data using a new data cube search method. These LAEs have luminosities comparable to those at high redshifts and lie within a 7 Gyr gap present in existing LAE samples. We use archival and newly obtained optical spectra to verify the UV redshifts of these LAEs. We use the combination of the GALEX UV spectra, optical spectra, and X-ray imaging data to estimate the active galactic nucleus (AGN) fraction and its dependence on Lya luminosity. We remove the AGNs and compute the luminosity function (LF) from 60 z~1 LAE galaxies. We find that the best fit LF implies a luminosity density increase by a factor of ~1.5 from z~0.3 to z~1 and ~20 from z~1 to z~2. We find a z~1 volumetric Lya escape fraction of 0.7+/-0.4%.

We present 1.4 GHz catalogs for the cluster fields Abell 370 and Abell 2390 observed with the Very Large Array. These are two of the deepest radio images of cluster fields ever taken. The Abell 370 image covers an area of 40'x40' with a synthesized beam of ~1.7" and a noise level of ~5.7 uJy near field center. The Abell 2390 image covers an area of 34'x34' with a synthesized beam of ~1.4" and a noise level of ~5.6 uJy near field center. We catalog 200 redshifts for the Abell 370 field. We construct differential number counts for the central regions (radius < 16') of both clusters. We find that the faint (S_1.4GHz < 3 mJy) counts of Abell 370 are roughly consistent with the highest blank field number counts, while the faint number counts of Abell 2390 are roughly consistent with the lowest blank field number counts. Our analyses indicate that the number counts are primarily from field radio galaxies. We suggest that the disagreement of our counts can be largely attributed to cosmic variance.

Recent cosmological modeling efforts have shown that a local underdensity on scales of a few hundred Mpc (out to z ~ 0.1), could produce the apparent acceleration of the expansion of the universe observed via type Ia supernovae. Several studies of galaxy counts in the near-infrared (NIR) have found that the local universe appears under-dense by ~25-50% compared with regions a few hundred Mpc distant. Galaxy counts at low redshifts sample primarily L ~ L* galaxies. Thus, if the local universe is under-dense, then the normalization of the NIR galaxy luminosity function (LF) at z>0.1 should be higher than that measured for z<0.1. Here we present a highly complete (> 90%) spectroscopic sample of 1436 galaxies selected in the H-band to study the normalization of the NIR LF at 0.1<z<0.3 and address the question of whether or not we reside in a large local underdensity. We find that for the combination of our six fields, the product phi* L* at 0.1 < z < 0.3 is ~ 30% higher than that measured at lower redshifts. While our statistical errors in this measurement are on the ~10% level, we find the systematics due to cosmic variance may be larger still. We investigate the effects of cosmic variance on our measurement using the COSMOS cone mock catalogs from the Millennium simulation and recent empirical estimates. We find that our survey is subject to systematic uncertainties due to cosmic variance at the 15% level ($1 sigma), representing an improvement by a factor of ~ 2 over previous studies in this redshift range. We conclude that observations cannot yet rule out the possibility that the local universe is under-dense at z<0.1.

We describe a method for obtaining a flux-limited sample of Ly-alpha emitters from GALEX grism data. We show that the multiple GALEX grism images can be converted into a three-dimensional (two spatial axes and one wavelength axis) data cube. The wavelength slices may then be treated as narrowband images and searched for emission-line galaxies. For the GALEX NUV grism data, the method provides a Ly-alpha flux-limited sample over the redshift range z=0.67-1.16. We test the method on the Chandra Deep Field South field, where we find 28 Ly-alpha emitters with faint continuum magnitudes (NUV>22) that are not present in the GALEX pipeline sample. We measure the completeness by adding artificial emitters and measuring the fraction recovered. We find that we have an 80% completeness above a Ly-alpha flux of 10^-15 erg/cm^2/s. We use the UV spectra and the available X-ray data and optical spectra to estimate the fraction of active galactic nuclei in the selection. We report the first detection of a giant Ly-alpha blob at z<1, though we find that these objects are much less common at z=1 than at z=3. Finally, we compute limits on the z~1 Ly-alpha luminosity function and confirm that there is a dramatic evolution in the luminosity function over the redshift range z=0-1.

There is increasing evidence of a connection between AGN activity and galaxy evolution. To obtain further insight into this potentially important evolutionary phase, we analyse the properties of quasar host galaxies. In this paper, we present a population synthesis modeling technique for off-axis spectra, the results of which constrain host colour and the stellar ages of luminous quasars (M_V(nuc)<-23). Our technique is similar to well established quiescent-galaxy models, modified to accommodate scattered nuclear light (a combination of atmospheric, instrumental and host galaxy scattered light) observed off axis. In our model, subtraction of residual scattered quasar light is performed, while simultaneously modeling the constituent stellar populations of the host galaxy. The reliability of this technique is tested via a Monte-Carlo routine in which the correspondence between synthetic spectra with known parameters and the model output is determined. Application of this model to a preliminary sample of 10 objects is presented and compared to previous studies. Spectroscopic data was obtained via long-slit and integral-field unit observations on the Keck and WIYN telescopes. We confirm that elliptical quasar hosts are distinguishable (bluer) from inactive ellipticals in rest frame B-V colour. Additionally, we note a trend for radio luminous (L_5GHz > 10^40 erg s^-1) quasars to be located in redder host galaxies in comparison to their less luminous radio counterparts. While the host colour and age of our radio luminous sample is in close proximity to the green valley, our radio faint sample is consistent with quiescent star-forming galaxies. However, further observations are needed to confirm these results. Finally, we discuss future applications for our technique on a larger sample of objects being obtained via SALT and WIYN telescope observing campaigns.

We present new multi-wavelength observations of the first submm-selected galaxy, SMM J02399-0136 at z=2.8. These observations include mapping of the CO J=1-0 emission using elements of the Expanded VLA, as well as high-resolution 1.4-GHz imaging and optical/IR data from the VLA, HST, Spitzer and Keck. Together these new data provide fundamental insights into the mass and distribution of stars, gas and dust within this archetypal SMG. The CO J=1-0 emission, with its minimal excitation and density requirements, traces the bulk of the metal-rich molecular gas, and reveals a mass of ~10^11 M(sun), extending over ~5" (~25 kpc in the source plane), although there is tentative evidence that it may be significantly larger. Our data suggest that three or more distinct structures are encompassed by this molecular gas reservoir, including the BAL quasar from which the redshift of the SMG was initially determined. In particular, the new rest-frame near-IR observations identify a massive, obscured, starburst coincident with a previously known Ly-alpha cloud. This starburst dominates the far-IR emission from the system and requires a re-assessment of previous claims that the gas reservoir resides in a massive, extended disk around the BAL QSO. Instead it appears that SMM J02399-0136 comprises a merger between a far-IR-luminous, but highly obscured starburst, the BAL QSO host and a faint third component. Our findings suggest that this SMG and its immediate environment mark a vast and complex galactic nursery and that detailed studies of other SMGs are likely to uncover a similarly rich diversity of properties.