Looking for bright galaxies born in the early universe is fundamental to investigating the Epoch of Reionization, the era when the first stars and galaxies ionized the intergalactic medium. We utilize Hubble Space Telescope pure parallel imaging to select galaxy candidates at a time 500 to 650 million years after the Big Bang, which corresponds to redshifts z ~ 8-10. These data come from the Brightest of Reionizing Galaxies Survey (BoRG) Cycle 22 dataset, which consists of pure-parallel imaging in ~ 90 different lines of sight that sum up to an area of ~ 420 arcmin^2. This survey uses five filters and has the advantage (compared to the Cycle 21 BoRG program) of including imaging in the JH140 band, covering continuous wavelengths from the visible to near-infrared (lambda = 0.35um - 1.7um). This allows us to perform reliable selection of galaxies at z>8 using the photometric redshift technique. We use these galaxy candidates to constrain the bright end of the rest-frame ultraviolet luminosity function in this epoch. These candidates are excellent targets for follow-up observations, particularly with the James Webb Space Telescope.

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

We present an ultraviolet spectroscopic survey of strong H I absorbers in the intergalactic medium, probing their evolution over the last 6-7 Gyr at redshifts $0.24 \leq z \leq 0.84$. We measure column densities $N_{\rm HI} \,( {\rm cm}^{-2})$ from the pattern of Lyman-series absorption lines and flux decrement at the Lyman limit (LL) when available. We analyzed 220 H I absorbers in ultraviolet spectra of 102 active galactic nuclei (AGN) taken by the Cosmic Origins Spectrograph aboard the Hubble Space Telescope with G130M/G160M gratings (1134-1795 \AA). For 158 absorbers with $\log N_{\rm HI} \geq 15$, the mean frequency is $d {\cal N}/dz = 4.95 \pm 0.39$ over pathlength $\Delta z = 31.94$ ($0.24 \leq z \leq 0.84)$. We identify 8 Lyman Limit Systems (LLS, $\log N_{\rm HI} \geq 17.2$) and 54 partial systems (pLLS) with $16.0 \leq \log N_{\rm HI} < 17.2$. Toward 159 AGN between $0.01 < z_{\rm abs} < 0.84$ with $\Delta z \approx 48$, we find four damped Ly$\alpha$ absorbers (DLA) with $(d {\cal N}/dz)_{\rm DLA} = 0.083^{+0.066}_{-0.040}$ at $\langle z \rangle = 0.18$. The mean LLS frequency between $z = 0.24-0.48$ is $(d {\cal N}/dz)_{\rm LLS} = 0.36^{+0.20}_{-0.13}$ fitted to $N(z) = (0.25^{+0.13}_{-0.09})(1+z)^{1.14}$. For 54 pLLS we find $(d {\cal N}/dz)_{\rm pLLS} = 1.69\pm0.23$ at $\langle z \rangle = 0.39$, a frequency consistent with gaseous halo sizes $R \approx 100 h^{-1}~{\rm kpc}$ for ($0.3-3L^*$) galaxies. A maximum-likelihood analysis yields a distribution $f(N,z) = C_0 N^{-\beta} (1+z)^{\gamma}$ with $\beta = 1.48 \pm 0.05$ and $\gamma = 1.14^{+0.88}_{-0.89}$ for $15 \leq \log N_{\rm HI} \leq 17.5$. The far-UV opacity gradient is $d \tau_{\rm eff} / dz \approx (0.444)(1+z)^{1.14}$ over the range $15 \leq \log N_{\rm HI} \leq 17$, implying mean LyC optical depth $\tau_{\rm eff} \approx 0.3-0.5$ toward sources at $z = 1-2$.

We use high-quality, medium-resolution {\it Hubble Space Telescope}/Cosmic Origins Spectrograph (\HST/COS) observations of 82 UV-bright AGN at redshifts $z_{AGN}<0.85$ to construct the largest survey of the low-redshift intergalactic medium (IGM) to date: 5343 individual extragalactic absorption lines in HI and 25 different metal-ion species grouped into 2610 distinct redshift systems at $z_{abs}<0.75$ covering total redshift pathlengths $\Delta z_{HI}=21.7$ and $\Delta z_{OVI}=14.5$. Our semi-automated line-finding and measurement technique renders the catalog as objectively-defined as possible. The cumulative column-density distribution of HI systems can be parametrized $dN(>N)/dz=C_{14}(N/10^{14} cm^{-2})^{-(\beta-1)}$, with $C_{14}=25\pm1$ and $\beta=1.65\pm0.02$. This distribution is seen to evolve both in amplitude, $C_{14}\sim(1+z)^{2.0\pm0.1}$, and slope $\beta(z)=1.73-0.26 z$ for $z<0.47$. We observe metal lines in 427 systems, and find that the fraction of IGM absorbers detected in metals is strongly dependent on N_{HI}. The distribution of OVI absorbers appear to evolve in the same sense as the Lya forest. We calculate contributions to $\Omega_b$ from different components of the low-$z$ IGM and determine the Lya decrement as a function of redshift. IGM absorbers are analyzed via a two-point correlation function (TPCF) in velocity space. We find substantial clustering of \HI\ absorbers on scales of $\Delta v=50-300$ km/s with no significant clustering at $\Delta v>1000$ km/s. Splitting the sample into strong and weak absorbers, we see that most of the clustering occurs in strong, $N_{HI}>10^{13.5} cm^{-2}$, metal-bearing IGM systems. The full catalog of absorption lines and fully-reduced spectra is available via MAST as a high-level science product at http://archive.stsci.edu/prepds/igm/.

The rest-frame ultraviolet (UV) spectra of active galactic nuclei (AGNs) are important diagnostics of both accretion disk physics and their contribution to the metagalactic ionizing UV background. Though the mean AGN spectrum is well characterized with composite spectra at wavelengths greater than 912 Angstroms, the shorter-wavelength extreme-UV (EUV) remains poorly studied. In this third paper in a series on the spectra of AGNs, we combine 11 new spectra taken with the Cosmic Origins Spectrograph on the Hubble Space Telescope with archival spectra to characterize the typical EUV spectral slope of AGNs from $\lambda_{\rm rest}\sim 850~{\rm Angstroms}$ down to $\lambda_{\rm rest}\sim 425~{\rm Angstroms}$. Parameterizing this slope as a power law, we obtain $F_\nu\propto \nu^{ -0.72\pm 0.26}$, but we also discuss the limitations and systematic uncertainties of this model. We identify broad emission features in this spectral region, including emission due to ions of O, Ne, Mg, and other species, and we limit the intrinsic HeI 504 Angstrom photoelectric absorption edge opacity to $\tau_{\rm HeI}<0.047$.

The ionizing fluxes from quasars and other active galactic nuclei (AGN) are critical for interpreting their emission-line spectra and for photoionizing and heating the intergalactic medium (IGM). Using far-ultraviolet spectra from the Cosmic Origins Spectrograph (COS) on the Hubble Space Telescope (HST), we directly measure the rest-frame ionizing continua and emission lines for 159 AGN at redshifts 0.001 < z_AGN < 1.476 and construct a composite spectrum from 475-1875A. We identify the underlying AGN continuum and strong EUV emission lines from ions of oxygen, neon, and nitrogen after masking out absorption lines from the HI Lya forest, 7 Lyman-limit systems (N_HI > 10^17.2 cm^-2) and 214 partial Lyman-limit systems (15.0 < log N_HI < 17.2). The 159 AGN exhibit a wide range of FUV/EUV spectral shapes, F_nu \propto nu^(alpha_nu), typically with -2 < alpha_nu < 0 and no discernible continuum edges at 912A (H I) or 504A (He I). The composite rest-frame continuum shows a gradual break at 1000 A, with mean spectral index alpha_nu = -0.83 +/- 0.09 in the FUV (1200-2000A) steepening to alpha_nu = -1.41 +/- 0.15 in the EUV (500-1000A). We discuss the implications of the UV flux turnovers and lack of continuum edges for the structure of accretion disks, AGN mass inflow rates, and luminosities relative to Eddington values.

The ionizing fluxes from quasars and other active galactic nuclei (AGN) are critical for interpreting the emission-line spectra of AGN and for photoionization and heating of the intergalactic medium. Using ultraviolet spectra from the Cosmic Origins Spectrograph (COS) on the Hubble Space Telescope (HST), we have directly measured the rest-frame ionizing continua and emission lines for 22 AGN. Over the redshift range 0.026 < z < 1.44, COS samples the Lyman continuum and many far-UV emission lines (Lya 1216, C IV 1549, Si IV/OIV] 1400, N V 1240, O VI 1035). Strong EUV emission lines with 14-22 eV excitation energies (Ne VIII 770,780, Ne V 569, O II 834, O III 833, 702, O IV 788,608,554, O V 630, N III 685) suggest the presence of hot gas in the broad emission-line region. The rest-frame continuum, F_nu ~ nu^{alpha_nu}, shows a break at wavelengths below 1000 A, with spectral index alpha_nu = -0.68 +/- 0.14 in the FUV (1200-2000 A) steepening to alpha_nu = -1.41 +/- 0.21 in the EUV (500-1000 A). The COS EUV index is similar to that of radio-quiet AGN in the 2002 HST/FOS survey (alpha_nu = -1.57 +/- 0.17). We see no Lyman edge (tau_HI < 0.03) or He I 584 emission in the AGN composite. Our 22 AGN exhibit a substantial range of FUV/EUV spectral indices and a correlation with AGN luminosity and redshift, likely due to observing below the 1000 A break.

We report ultraviolet spectra of Galactic high-velocity clouds (HVCs) in Complex C, taken by the Cosmic Origins Spectrograph (COS) on the Hubble Space Telescope (HST), together with new 21-cm spectra from the Green Bank Telescope. The wide spectral coverage and higher S/N, compared to previous HST spectra, provide better velocity definition of the HVC absorption, additional ionization species, and improved abundances in this halo gas. Complex C has a metallicity of 0.1-0.3 solar and a wide range of ions, suggesting dynamical and thermal interactions with hot gas in the Galactic halo. Spectra in the COS medium-resolution G130M (1133-1468 A) and G160M (1383-1796 A) gratings detect ultraviolet absorption lines from 8 elements in low ionization stages (O I, N I, C II, S II, Si II, Al II, Fe II, P II) and 3 elements in intermediate and high-ionization states (Si III, Si IV, C IV, N V). Our four AGN sight lines toward Mrk 817, Mrk 290, Mrk 876, and PG1259+593 have high-velocity H I and O VI column densities, log N_HI = 19.39-20.05 and log N_OVI = 13.58-14.10, with substantial amounts of kinematically associated photoionized gas. The high-ion abundance ratios are consistent with cooling interfaces between photoionized gas and collisionally ionized gas: N(C IV)/N(O VI) = 0.3-0.5, N(Si IV)/N(O VI) = 0.05-0.11, N(N V)/N(O VI) = 0.07-0.13, and N(Si IV)/N(Si III) = 0.2.