Search

In theoretical models of galaxy evolution, AGN and star formation (SF) activity are closely linked and AGN feedback is routinely invoked to regulate galaxy growth. In order to constrain such models, we compare the hydrodynamical simulations IllustrisTNG and SIMBA, and the semi-analytical model SAG to the empirical results on AGN and SF at cosmic noon ($0.75 < z < 2.25$) reported in Florez et al. (2020). The empirical results are based on a large mass-complete sample drawn from 93,307 galaxies with and without high X-ray luminosity AGN ($L_X \gtrsim 10^{44}$ erg s$^{-1}$), selected from a 11.8 deg$^2$ area ($\sim 0.18$ Gpc$^3$ comoving volume at $z=0.75-2.25$). The main results of our comparisons are: (i) SAG and IllustrisTNG both qualitatively reproduce the empirical result that galaxies with high X-ray luminosity AGN have higher mean SFR, at a given stellar mass, than galaxies without such AGN. SAG, however, strongly over-produces the number density of high X-ray luminosity AGN by a factor of 10 to 100, while IllustrisTNG shows a lack of high X-ray luminosity AGN at high stellar mass ($M* > 10^{11} \ M_{\odot}$) at $z \sim 2$. (ii) In SIMBA, the mean SFR of galaxies with high X-ray luminosity AGN is lower than the SFR of galaxies without such AGN. Contrary to the data, many high X-ray luminosity AGN in SIMBA have quenched SF, suggesting that AGN feedback, or other feedback modes in galaxies with such AGN, might be too efficient in SIMBA., Comment: 20 pages, 7 figures. Accepted for publication in MNRAS

We present the main sequence for all galaxies and star-forming galaxies for a sample of 28,469 massive ($M_\star \ge 10^{11}$M$_\odot$) galaxies at cosmic noon ($1.5 < z < 3.0$), uniformly selected from a 17.5 deg$^2$ area (0.33 Gpc$^3$ comoving volume at these redshifts). Our large sample allows for a novel approach to investigating the galaxy main sequence that has not been accessible to previous studies. We measure the main sequence in small mass bins in the SFR-M$_{\star}$ plane without assuming a functional form for the main sequence. With a large sample of galaxies in each mass bin, we isolate star-forming galaxies by locating the transition between the star-forming and green valley populations in the SFR-M$_{\star}$ plane. This approach eliminates the need for arbitrarily defined fixed cutoffs when isolating the star-forming galaxy population, which often biases measurements of the scatter around the star-forming galaxy main sequence. We find that the main sequence for all galaxies becomes increasingly flat towards present day at the high-mass end, while the star-forming galaxy main sequence does not. We attribute this difference to the increasing fraction of the collective green valley and quiescent galaxy population from $z=3.0$ to $z=1.5$. Additionally, we measure the total scatter around the star-forming galaxy main sequence and find that it is $\sim0.5-1.0$ dex with little evolution as a function of mass or redshift. We discuss the implications that these results have for pinpointing the physical processes driving massive galaxy evolution., Comment: 16 pages, 10 figures, Accepted for publication in MNRAS

We present the results of a deep K_s-band (2.1 um) imaging survey of the Spitzer/HETDEX Exploratory Large Area (SHELA) field using the NEWFIRM near-infrared (NIR) camera on the KPNO Mayall 4-m telescope. This NEWFIRM HETDEX Survey (NHS) reaches a 5-sigma depth of 22.4 AB mag (2"-diameter apertures corrected to total), is ~50% and 90% complete at K~22.65 and K~22.15, respectively, and covers 22 deg^2 of the 24 deg^2 SHELA Spitzer/ IRAC footprint (within ``Stripe 82' of the Sloan Digital Sky Survey). We present a K_s-band-selected catalog which includes deep ugriz imaging from the Dark Energy Camera and 3.6 and 4.5 um imaging from Spitzer/IRAC, with forced-photometry of 1.7 million sources across 17.5 deg^2. The large area and moderate depth of this catalog enables the study of the most massive galaxies at high redshift, and minimizes uncertainties associated with counting statistics and cosmic variance. As a demonstration, we derive stellar masses (M*) and star-formation rates (SFRs) for candidate galaxies at 3 < z < 5, and select a conservative sample of nine candidate massive (M* > 10^11 Msol) quiescent galaxies, which have measured SFRs significantly below the main-sequence at this redshift. Five are ultra-massive with M* > 10^12, though uncertainties in IRAC blending, gravitational lensing, or AGN emission could result in true masses which are lower. Simulations predict that these galaxies should be extremely rare, thus we discuss what physical processes in models could be altered to allow the formation of such massive quiescent galaxies at such early times., Comment: 38 pages, accepted to the Astrophysical Journal. The full catalog is included as an ancillary file, and is be published with the journal article

We measure the environmental dependence, where environment is defined by the distance to the third nearest neighbor, of multiple galaxy properties inside the Environmental COntext (ECO) catalog. We focus primarily on void galaxies, which we define as the $10 \%$ of galaxies having the lowest local density. We compare the properties of void and non-void galaxies: baryonic mass, color, fractional stellar mass growth rate (FSMGR), morphology, and gas-to-stellar-mass ratio (estimated from a combination of HI data and photometric gas fractions calibrated with the RESOLVE survey). Our void galaxies typically have lower baryonic masses than galaxies in denser environments, and they display the properties expected of a lower mass population: they have more late-types, are bluer, have higher FSMGR, and are more gas rich. We control for baryonic mass and investigate the extent to which void galaxies are different at fixed mass. Void galaxies are bluer, more gas-rich, and more star forming at fixed mass than non-void galaxies, which is a possible signature of galaxy assembly bias. Furthermore, we show that these trends persist even at fixed mass and morphology, and we find that voids host a distinct population of early-types that are bluer and more star-forming than the typical red and quenched early-types. In addition to these empirical observational results, we also present theoretical results from mock catalogs with built-in galaxy assembly bias. We show that a simple matching of galaxy properties to (sub)halo properties, such as mass and age, can recover the observed environmental trends in ECO galaxies., Comment: 18 pages, 11 figures. Accepted for publication in ApJ

We explore the buildup of quiescent galaxies using a sample of 28,469 massive ($M_\star \ge 10^{11}$M$_\odot$) galaxies at redshifts $1.5

We investigate the relation between AGN and star formation (SF) activity at $0.5 < z < 3$ by analyzing 898 galaxies with X-ray luminous AGN ($L_X > 10^{44}$ erg s$^{-1}$) and a large comparison sample of $\sim 320,000$ galaxies without X-ray luminous AGN. Our samples are selected from a large (11.8 deg$^2$) area in Stripe 82 that has multi-wavelength (X-ray to far-IR) data. The enormous comoving volume ($\sim 0.3$ Gpc$^3$) at $0.5 < z < 3$ minimizes the effects of cosmic variance and captures a large number of massive galaxies ($\sim 30,000$ galaxies with $M_* > 10^{11} \ M_{\odot}$) and X-ray luminous AGN. While many galaxy studies discard AGN hosts, we fit the SED of galaxies with and without X-ray luminous AGN with Code Investigating GALaxy Emission (CIGALE) and include AGN emission templates. We find that without this inclusion, stellar masses and star formation rates (SFRs) in AGN host galaxies can be overestimated, on average, by factors of up to $\sim 5$ and $\sim 10$, respectively. The average SFR of galaxies with X-ray luminous AGN is higher by a factor of $\sim 3$ to $10$ compared to galaxies without X-ray luminous AGN at fixed stellar mass and redshift, suggesting that high SFRs and high AGN X-ray luminosities may be fueled by common mechanisms. The vast majority ($> 95 \%$) of galaxies with X-ray luminous AGN at $z=0.5-3$ do not show quenched SF: this suggests that if AGN feedback quenches SF, the associated quenching process takes a significant time to act and the quenched phase sets in after the highly luminous phases of AGN activity., Comment: 27 pages, 15 figures. Accepted for publication in MNRAS

We study the evolution in the number density of the highest mass galaxies over $0.410.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$., Comment: 39 pages, 23 figure, accepted for publication in ApJ

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., Comment: 19 pages, 10 figures, Accepted for publication in MNRAS

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

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., Comment: 25 pages, 13 figures, accepted for publication in ApJ

We present the HI mass inventory for the RESOLVE survey, a volume-limited, multi-wavelength census of >1500 z=0 galaxies spanning diverse environments and complete in baryonic mass down to dwarfs of 10^9 Msun. This first 21cm data release provides robust detections or strong upper limits (1.4M_HI < 5 to 10% of stellar mass M_stars) for 94% of RESOLVE. We examine global atomic gas-to-stellar mass ratios (G/S) in relation to galaxy environment using several metrics: group dark matter halo mass M_h , central/satellite designation, relative mass density of the cosmic web, and distance to nearest massive group. We find that at fixed M_stars, satellites have decreasing G/S with increasing M_h starting clearly at M_h = 10^12 Msun, suggesting the presence of starvation and/or stripping mechanisms associated with halo gas heating in intermediate-mass groups. The analogous relationship for centrals is uncertain because halo abundance matching builds in relationships between central G/S, stellar mass, and halo mass, which depend on the integrated group property used as a proxy for halo mass (stellar or baryonic mass). On larger scales G/S trends are less sensitive to the abundance matching method. At fixed M_h < 10^12 Msun, the fraction of gas-poor centrals increases with large-scale structure density. In overdense regions, we identify a rare population of gas-poor centrals in low-mass (M_h < 10^11.4 Msun) halos primarily located within 1.5 times the virial radius of more massive (M_h > 10^12 Msun) halos, suggesting that gas stripping and/or starvation may be induced by interactions with larger halos or the surrounding cosmic web. We find that the detailed relationship between G/S and environment varies when we examine different subvolumes of RESOLVE independently, which we suggest may be a signature of assembly bias., Comment: 24 pages, 13 figures, accepted for publication in ApJ, data tables available at http://resolve.astro.unc.edu/pages/database.php