Search

We present the Empirical Dust Attenuation (EDA) framework -- a flexible prescription for assigning realistic dust attenuation to simulated galaxies based on their physical properties. We use the EDA to forward model synthetic observations for three state-of-the-art large-scale cosmological hydrodynamical simulations: SIMBA, IllustrisTNG, and EAGLE. We then compare the optical and UV color-magnitude relations, $(g-r) - M_r$ and $(FUV-NUV)-M_r$, of the simulations to a $M_r < -20$ and UV complete SDSS galaxy sample using likelihood-free inference. Without dust, none of the simulations match observations, as expected. With the EDA, however, we can reproduce the observed color-magnitude with all three simulations. Furthermore, the attenuation curves predicted by our dust prescription are in good agreement with the observed attenuation-slope relations and attenuation curves of star-forming galaxies. However, the EDA does not predict star-forming galaxies with low $A_V$ since simulated star-forming galaxies are intrinsically much brighter than observations. Additionally, the EDA provides, for the first time, predictions on the attenuation curves of quiescent galaxies, which are challenging to measure observationally. Simulated quiescent galaxies require shallower attenuation curves with lower amplitude than star-forming galaxies. The EDA, combined with forward modeling, provides an effective approach for shedding light on dust in galaxies and probing hydrodynamical simulations. This work also illustrates a major limitation in comparing galaxy formation models: by adjusting dust attenuation, simulations that predict significantly different galaxy populations can reproduce the same UV and optical observations., 26 pages, 15 figures

We compare three major large-scale hydrodynamical galaxy simulations (EAGLE, Illustris-TNG, and SIMBA) by forward modeling simulated galaxies into observational space and computing the fraction of isolated and quiescent low mass galaxies as a function of stellar mass. Using SDSS as our observational template, we create mock surveys and synthetic spectroscopic and photometric observations of each simulation, adding realistic noise and observational limits. All three simulations show a decrease in the number of quiescent, isolated galaxies in the mass range $\mathrm{M}_* = 10^{9-10} \ \mathrm{M}_\odot$, in broad agreement with observations. However, even after accounting for observational and selection biases, none of the simulations reproduce the observed absence of quiescent field galaxies below $\mathrm{M}_*=10^{9} \ \mathrm{M}_\odot$. We find that the low mass quiescent populations selected via synthetic observations have consistent quenching timescales, despite apparent variation in the late time star formation histories. The effect of increased numerical resolution is not uniform across simulations and cannot fully mitigate the differences between the simulations and the observations. The framework presented here demonstrates a path towards more robust and accurate comparisons between theoretical simulations and galaxy survey observations, while the quenching threshold serves as a sensitive probe of feedback implementations., Comment: 19 pages, 8 figures. Figure 4 presents the main result. Code used in this work may be accessed at github.com/IQcollaboratory/orchard. Submitted to ApJ

We investigate the lowest-mass quiescent galaxies known to exist in isolated environments ($\mathrm{M^* = 10^{9.0-9.5} \ M_\odot}$; 1.5 Mpc from a more massive galaxy). This population may represent the lowest stellar mass galaxies in which internal feedback quenches galaxy-wide star formation. We present Keck/ESI long-slit spectroscopy for 27 isolated galaxies in this regime: 20 quiescent galaxies and 7 star-forming galaxies. We measure emission line strengths as a function of radius and place galaxies on the Baldwin Phillips Terlevich (BPT) diagram. Remarkably, 16 of 20 quiescent galaxies in our sample host central AGN-like line ratios. Only 5 of these quiescent galaxies were identified as AGN-like in SDSS due to lower spatial resolution and signal-to-noise. We find that many of the quiescent galaxies in our sample have spatially-extended emission across the non-SF regions of BPT-space. When considering only the central 1$^{\prime\prime}$, we identify a tight relationship between distance from the BPT star-forming sequence and host galaxy stellar age as traced by $\mathrm{D_n4000}$, such that older stellar ages are associated with larger distances from the star-forming locus. Our results suggest that the presence of hard ionizing radiation (AGN-like line ratios) is intrinsically tied to the quenching of what may be the lowest-mass self-quenched galaxies., Comment: 16 pages, 9 figures, submitted to ApJ

A tightly correlated star formation rate-stellar mass relation of star forming galaxies, or star-forming sequence (SFS), is a key feature in galaxy property-space that is predicted by modern galaxy formation models. We present a flexible data-driven approach for identifying this SFS over a wide range of star formation rates and stellar masses using Gaussian mixture modeling (GMM). Using this method, we present a consistent comparison of the $z=0$ SFSs of central galaxies in the Illustris, EAGLE, and Mufasa hydrodynamic simulations and the Santa Cruz semi-analytic model (SC-SAM), alongside data from the Sloan Digital Sky Survey. We find, surprisingly, that the amplitude of the SFS varies by up to ${\sim} 0.7\,\mathrm{dex}$ (factor of ${\sim} 5$) among the simulations with power-law slopes range from $0.7$ to $1.2$. In addition to the SFS, our GMM method also identifies sub-components in the star formation rate-stellar mass relation corresponding to star-burst, transitioning, and quiescent sub-populations. The hydrodynamic simulations are similarly dominated by SFS and quiescent sub-populations unlike the SC-SAM, which predicts substantial fractions of transitioning and star-burst galaxies at stellar masses above and below $10^{10} M_\odot$, respectively. All of the simulations also produce an abundance of low-mass quiescent central galaxies in apparent tension with observations. These results illustrate that, even among models that well reproduce many observables of the galaxy population, the $z=0$ SFS and other sub-populations still show marked differences that can provide strong constraints on galaxy formation models., Comment: 28 pages, 15 figures

We investigate the global neutral hydrogen (HI) content of isolated galaxies selected from the SDSS spectroscopic survey with optical evidence of Active Galactic Nuclei (AGN). Our sample includes galaxies with unresolved HI observations from the ALFALFA 70% data release combined with deeper HI observations of low-mass galaxies with 7.0 < log(M*) < 9.5. We examine the HI masses of this sample using the distance from the star-forming sequence on the OIII\Hb and NII\Ha Baldwin Phillips Terlevich (BPT) diagram as a measurement of AGN activity. In agreement with previous studies, we find that, for galaxies with log(M*) > 9.5, AGN activity does not correlate with the global HI content. However, for galaxies with 9.2 < log(M*) < 9.5, we identify a set of objects at large distances from the BPT star-forming sequence and lower than expected HI masses. This gas-depleted sample is red in both g-r and NUV-r colors and compact without distinguishable signs of star formation. This is surprising because the vast majority of isolated galaxies in this stellar mass regime are both star-forming and gas-rich. These galaxies are greater than 1.5 Mpc from any massive galaxy, ruling out environmental processes as a source of the gas-depletion. We suggest that either black hole feedback or shocks from extremely bursty star formation cause the emission lines and have destroyed or otherwise consumed the cold gas., Accepted to ApJ, 14 pages, 7 figures

We present Hα maps at 1 kpc spatial resolution for star-forming galaxies at z ∼ 1, made possible by the Wide Field Camera 3 grism on Hubble Space Telescope (HST). Employing this capability over all five 3D-HST/CANDELS fields provides a sample of 3200 galaxies enabling a division into subsamples based on stellar mass and star formation rate (SFR). By creating deep stacked Hα images, we reach surface brightness limits of 1 × 10-18 erg s-1 cm-2 arcsec-2, allowing us to map the distribution of ionized gas to ∼10 kpc for typical L∗ galaxies at this epoch. We find that the spatial extent of the Hα distribution increases with stellar mass as rHα = 1.5 (M∗/1010 M⊙)0.23 kpc. The Hα emission is more extended than the stellar continuum emission, consistent with inside-out assembly of galactic disks. This effect grows stronger with mass as rHα/r∗ = 1.1(M∗/1010 M⊙)0.054. We map the Hα distribution as a function of SFR(IR+UV) and find evidence for "coherent star formation" across the SFR-M∗ plane: above the main sequence (MS), Hα is enhanced at all radii; below the MS, Hα is depressed at all radii. This suggests that at all masses the physical processes driving the enhancement or suppression of star formation act throughout the disks of galaxies. At high masses (1010.5 ∗/M⊙ 11), above the MS, Hα is particularly enhanced in the center, potentially building bulges and/or supermassive black holes. Below the MS, a strong central dip in the EW(Hα), as well as the inferred specific SFR, appears. Importantly, though, across the entirety of the SFR-M∗ plane, the absolute SFR as traced by Hα is always centrally peaked, even in galaxies below the MS.

We present Ha maps at 1kpc spatial resolution for star-forming galaxies at z~1, made possible by the WFC3 grism on HST. Employing this capability over all five 3D-HST/CANDELS fields provides a sample of 2676 galaxies. By creating deep stacked Halpha (Ha) images, we reach surface brightness limits of 1x10^-18\erg\s\cm^2\arcsec^2, allowing us to map the distribution of ionized gas out to >10kpc for typical L* galaxies at this epoch. We find that the spatial extent of the Ha distribution increases with stellar mass as r(Ha)[kpc]=1.5(Mstars/10^10Msun)^0.23. Furthermore, the Ha emission is more extended than the stellar continuum emission, consistent with inside-out assembly of galactic disks. This effect, however, is mass dependent with r(Ha)/r(stars)=1.1(M/10^10Msun)^0.054, such that at low masses r(Ha)~r(stars). We map the Ha distribution as a function of SFR(IR+UV) and find evidence for `coherent star formation' across the SFR-M plane: above the main sequence, Ha is enhanced at all radii; below the main sequence, Ha is depressed at all radii. This suggests that at all masses the physical processes driving the enhancement or suppression of star formation act throughout the disks of galaxies. It also confirms that the scatter in the star forming main sequence is real and caused by variations in the star formation rate at fixed mass. At high masses (10^10.5, Comment: Submitted to ApJ. Main results: Fig. 6 shows the scale lengths of ionized gas vs. stellar disks. Fig. 9 shows where Halpha in galaxies is enhanced above and suppressed below the star-forming main sequence. Data points for Halpha profiles are available upon request

We present reduced data and data products from the 3D-HST survey, a 248-orbit HST Treasury program. The survey obtained WFC3 G141 grism spectroscopy in four of the five CANDELS fields: AEGIS, COSMOS, GOODS-S, and UDS, along with WFC3 $H_{140}$ imaging, parallel ACS G800L spectroscopy, and parallel $I_{814}$ imaging. In a previous paper (Skelton et al. 2014) we presented photometric catalogs in these four fields and in GOODS-N, the fifth CANDELS field. Here we describe and present the WFC3 G141 spectroscopic data, again augmented with data from GO-1600 in GOODS-N. The data analysis is complicated by the fact that no slits are used: all objects in the WFC3 field are dispersed, and many spectra overlap. We developed software to automatically and optimally extract interlaced 2D and 1D spectra for all objects in the Skelton et al. (2014) photometric catalogs. The 2D spectra and the multi-band photometry were fit simultaneously to determine redshifts and emission line strengths, taking the morphology of the galaxies explicitly into account. The resulting catalog has 98,663 measured redshifts and line strengths down to $JH_{IR}\leq 26$ and 22,548 with $JH_{IR}\leq 24$, where we comfortably detect continuum emission. Of this sample 5,459 galaxies are at $z>1.5$ and 9,621 are at $0.7, Comment: v2 updates Figure 22. Submitted to ApJS. All 2D and 1D spectra, as well as redshifts, line fluxes, and other derived parameters, are publicly available at available at http://3dhst.research.yale.edu