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 explore how the estimated star formation rate (SFR) of a sample of isolated, massive dusty star-forming galaxies at early cosmic epochs ($1.5 < z < 3.5$) changes when their ultraviolet (UV) to near-infrared (NIR) spectral energy distribution is extended to longer wavelengths by adding far-infrared/sub-millimeter data to trace the reprocessed radiation from dust heated by young massive stars. We use large-area surveys with multi-wavelength datasets that include DECam UV-to-optical, VICS82 NIR, Spitzer-IRAC NIR, and Herschel-SPIRE far-infrared/sub-millimeter data. We find that the inclusion of far-infrared/sub-millimeter data leads to SFRs that span $\sim$100-3500 $M_{\odot} yr^{-1}$ and are higher than the extinction-corrected UV-based SFR by an average factor of $\sim$3.5, and by a factor of over 10 in many individual galaxies. Our study demonstrates the importance of far-IR/sub-millimeter data for deriving accurate SFRs in massive dusty galaxies at early epochs, and underscores the need for next-generation far-IR/sub-millimeter facilities with high sensitivity, field of view, and angular resolution.

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.

Observations have shown that supermassive black holes in nearby elliptical galaxies correlate tightly with the stellar velocity dispersion (the $\MBH - \sigma$ relation) and the stellar mass (the $\MBH - \Mhost$ relation) of their host spheroids. However, the origin of these correlations remains ambiguous. In a previous paper by Zhu et al., we proposed a model which links the M-$\sigma$ relation to the the dynamical state of the system and the $\MBH - \Mhost$ relation to the self-regulation of galaxy growth. To test this model, we compile a sample of observed galaxies with different properties and examine the dependence of the above correlations on these parameters. We find that galaxies that satisfy the the $\MBH - \sigma$ correlation appear to have reached virial equilibrium, as indicated by the ratio between kinetic energy and gravitational potential, 2K/U $\sim$ 1. Furthermore, the ratio of black hole accretion rate to star formation rate remains nearly constant, BHAR /SFR $\sim$ $10^{-3}$, in active galaxies over a wide range of mass in the redshift range z=0 - 3. These results confirm our theoretical model that the observed correlations have different origins: the $\MBH - \sigma$ relation may result from galaxy relaxation, while the $\MBH$ - $\Mhost$ relation may be due to self-regulated black hole accretion and star formation in galaxies.

It is well established that supermassive black holes in nearby elliptical galaxies correlate tightly with the kinematic property ($\mbhsigma$ correlation) and stellar mass ($\mbhhost$ correlation) of their host spheroids. However, it is not clear what the relations would be at the low-mass end, and how they evolve. Here, we investigate these relations in low-mass systems ($\MBH \sim \rm{10^{6}- 10^{8}}\, \Msun$) using the Aquila Simulation, a high-resolution cosmological hydrodynamic simulation which follows the formation and evolution of stars and black holes in a Milky Way-size galaxy and its substructures. We find a number of interesting results on the origin and evolution of the scaling relations in these systems: (1) there is a strong redshift evolution in the $\mbhsigma$ relation, but a much weaker one in the $\mbhhost$ relation; (2) there is a close link between the $\mbhsigma$ relation and the dynamical state of the system -- the galaxies that fall on the observed correlation appear to have reached virial equilibrium. (3) the star formation and black hole growth are self-regulated in galaxies -- the ratio between black hole accretion rate and star formation rate remains nearly constant in a wide redshift span $z = 0-6$. These findings suggest that the observed correlations have different origins: the $\mbhsigma$ relation may be the result of virial equilibrium, while the $\mbhhost$ relation may the result of self-regulated star formation and black hole growth in galaxies.