Showing 1–50 of 302 results for author: Faber, S M

Non-Monotonic Relations of Galaxy Star Formation, Radius, and Structure at Fixed Stellar Mass

Authors: Jimena Stephenson, Aldo Rodriguez-Puebla, S. M. Faber, Joel R. Primack, Vladimir Avila-Reese, A. R. Calette, Carlo Cannarozzo, James Kakos, Mariana Cano-Díaz, David C. Koo, Francesco Shankar, D. F. Morell

Abstract: We investigate the relation between galaxy structure and star formation rate (SFR) in a sample of $\sim2.9\times10^{4}$ central galaxies with $z<0.0674$ and axial ratios $b/a>0.5$. The star-forming main sequence (SFMS) shows a bend around the stellar mass of $M_\ast\leq{}M_c=2\times10^{10}{}M_{\odot}$. At $M_\ast\leq{}M_c$ the SFMS follows a power-law $\text{SFR}\propto{}M_\ast^{0.85}$, while at h… ▽ More We investigate the relation between galaxy structure and star formation rate (SFR) in a sample of $\sim2.9\times10^{4}$ central galaxies with $z<0.0674$ and axial ratios $b/a>0.5$. The star-forming main sequence (SFMS) shows a bend around the stellar mass of $M_\ast\leq{}M_c=2\times10^{10}{}M_{\odot}$. At $M_\ast\leq{}M_c$ the SFMS follows a power-law $\text{SFR}\propto{}M_\ast^{0.85}$, while at higher masses it flattens. $M_c$ corresponds to a dark matter halo mass of $M_\text{vir}\sim{}10^{11.8}M_{\odot}$ where virial shocks occurs. Some galaxy structure (e.g., half-light radius, $R_e$) exhibits a non-monotonic dependence across the SFMS at a fixed $M_\ast$. We find $\text{SFR}\propto{R_e^{-0.28}}$ at fixed $M_\ast$, consistent with the global Kennicutt-Schmidt (KS) law. This finding suggests that galaxy sizes contribute to the scatter of the SFMS. However, at $M_\ast>M_c$ the relationship between SFR and $R_e$ diminishes. Low-mass galaxies above the mean of the SFMS have smaller radii, exhibit compact and centrally concentrated profiles resembling green valley (GV) and quiescent galaxies at the same mass, and have higher $M_{\text{H}_2}/M_\text{HI}$. Conversely, those below the SFMS exhibit larger radii, lower densities, have no GV or quiescent counterparts at their mass and have lower $M_{\text{H}_2}/M_\text{HI}$. The above data suggest two pathways for quenching low-mass galaxies, $M_\ast\leq{}M_c$: a fast one that changes the morphology on the SFMS and a slow one that does not. Above $M_c$, galaxies below the SFMS resemble GV and quiescent galaxies structurally, implying that they undergo a structural transformation already within the SFMS. For these massive galaxies, CG are strongly bimodal, with SFMS galaxies exhibiting negative color gradients, suggesting most star formation occurs in their outskirts, maintaining them within the SFMS. △ Less

Submitted 15 April, 2024; originally announced April 2024.

Comments: 30 pages, 28 figures, submitted to MNRAS

Star-forming and Quiescent Central Galaxies Cluster Similarly: Implications for the Galaxy-Halo Connection

Authors: James Kakos, Aldo Rodriguez-Puebla, Joel R. Primack, Sandra M. Faber, David C. Koo, Peter Behroozi, Vladimir Avila-Reese

Abstract: We measure the clustering of low-redshift SDSS galaxies as a function of stellar mass ($10.0<\log(M_*/M_\odot)<11.5$) and specific star formation rate (sSFR) and compare the results to models of the galaxy--halo connection. We find that the auto-correlation functions of central galaxies exhibit little dependence on sSFR, with the well-known stronger clustering of quiescent galaxies mainly attribut… ▽ More We measure the clustering of low-redshift SDSS galaxies as a function of stellar mass ($10.0<\log(M_*/M_\odot)<11.5$) and specific star formation rate (sSFR) and compare the results to models of the galaxy--halo connection. We find that the auto-correlation functions of central galaxies exhibit little dependence on sSFR, with the well-known stronger clustering of quiescent galaxies mainly attributable to satellites. Because halo assembly history is known to affect distinct halo clustering, this result implies that there is little net correlation between halo assembly history and central galaxy sSFR. However, cross-correlations with satellites are stronger for quiescent centrals than star-forming centrals, consistent with quiescent centrals having more satellites in their haloes at fixed $M_*$, as found in SDSS group catalogues. We model the galaxy--halo connection in an $N$-body simulation by assigning sSFRs to central galaxies in three different ways. Two of the models depend on halo assembly history (being based on halo accretion rate or concentration), while the third is independent of halo assembly history (being based on peak halo circular velocity, $V_\text{peak}$, a proxy for halo mass). All three models replicate the observed auto-correlations of central galaxies, while only the $V_\text{peak}$ model reproduces the observed cross-correlations with satellites. This further suggests that the effects of halo assembly history may not be easily seen in auto-correlations of centrals and implies that a more complete understanding of central galaxy clustering may require more than auto-correlations of centrals alone. Additionally, the good agreement with the $V_\text{peak}$ model supports the idea that quiescent galaxies reside in more massive haloes than star-forming galaxies at fixed $M_*$. △ Less

Submitted 2 March, 2024; originally announced March 2024.

Comments: MNRAS submitted, 23 pages

Stellar Populations With Optical Spectra: Deep Learning vs. Popular Spectrum Fitting Codes

Authors: Joanna Woo, Dan Walters, Finn Archinuk, S. M. Faber, Sara L. Ellison, Hossen Teimoorinia, Kartheik Iyer

Abstract: We compare the performance of several popular spectrum fitting codes (Firefly, starlight, pyPipe3D and pPXF), and a deep-learning convolutional neural network (StarNet), in recovering known stellar population properties (mean stellar age, stellar metallicity, stellar mass-to-light ratio M*/L_r and the internal E(B-V)) of simulated galaxy spectra in optical wavelengths. Our mock spectra are constru… ▽ More We compare the performance of several popular spectrum fitting codes (Firefly, starlight, pyPipe3D and pPXF), and a deep-learning convolutional neural network (StarNet), in recovering known stellar population properties (mean stellar age, stellar metallicity, stellar mass-to-light ratio M*/L_r and the internal E(B-V)) of simulated galaxy spectra in optical wavelengths. Our mock spectra are constructed from star-formation histories from the IllustrisTNG100-1 simulation. These spectra mimic the Sloan Digital Sky Survey (SDSS) through a novel method of including the noise, sky residuals and emission lines taken directly from SDSS. We find that StarNet vastly outperforms all conventional codes in both speed and recovery of stellar population properties (error scatter < 0.08 dex, average biases < 0.02 dex for all tested quantities), but it requires an appropriate training set. Of the non-machine-learning codes, pPXF was a factor of 3-4 times faster than the other codes, and was the best in recovering stellar population properties (error scatter of < 0.11 dex, average biases < 0.08 dex). However, the errors and biases are strongly dependent on both true and predicted values of stellar age and metallicity, and signal-to-noise ratio. The biases of all codes can approach 0.15 dex in stellar ages, metallicities and log M*/L_r , but remain < 0.05 for E(B-V). Using unrealistic Gaussian noise in the construction of mock spectra will underestimate the errors in the metallicities by a factor of two or more, and mocks without emission lines will underestimate the errors in stellar age and M*/L_r by a factor of two. △ Less

Submitted 25 April, 2024; v1 submitted 22 January, 2024; originally announced January 2024.

Comments: Accepted to MNRAS

Filaments of The Slime Mold Cosmic Web And How They Affect Galaxy Evolution

Authors: Farhanul Hasan, Joseph N. Burchett, Douglas Hellinger, Oskar Elek, Daisuke Nagai, S. M. Faber, Joel R. Primack, David C. Koo, Nir Mandelker, Joanna Woo

Abstract: We present a novel method for identifying cosmic web filaments using the IllustrisTNG (TNG100) cosmological simulations and investigate the impact of filaments on galaxies. We compare the use of cosmic density field estimates from the Delaunay Tessellation Field Estimator (DTFE) and the Monte Carlo Physarum Machine (MCPM), which is inspired by the slime mold organism, in the DisPerSE structure ide… ▽ More We present a novel method for identifying cosmic web filaments using the IllustrisTNG (TNG100) cosmological simulations and investigate the impact of filaments on galaxies. We compare the use of cosmic density field estimates from the Delaunay Tessellation Field Estimator (DTFE) and the Monte Carlo Physarum Machine (MCPM), which is inspired by the slime mold organism, in the DisPerSE structure identification framework. The MCPM-based reconstruction identifies filaments with higher fidelity, finding more low-prominence/diffuse filaments and better tracing the true underlying matter distribution than the DTFE-based reconstruction. Using our new filament catalogs, we find that most galaxies are located within 1.5-2.5 Mpc of a filamentary spine, with little change in the median specific star formation rate and the median galactic gas fraction with distance to the nearest filament. Instead, we introduce the filament line density, Sigma_fil(MCPM), as the total MCPM overdensity per unit length of a local filament segment, and find that this parameter is a superior predictor of galactic gas supply and quenching. Our results indicate that most galaxies are quenched and gas-poor near high-line density filaments at z<=1. At z=0, quenching in log(M*/Msun)>10.5 galaxies is mainly driven by mass, while lower-mass galaxies are significantly affected by the filament line density. In high-line density filaments, satellites are strongly quenched, whereas centrals have reduced star formation, but not gas fraction, at z<=0.5. We discuss the prospect of applying our new filament identification method to galaxy surveys with SDSS, DESI, Subaru PFS, etc. to elucidate the effect of large-scale structure on galaxy formation. △ Less

Submitted 13 May, 2024; v1 submitted 2 November, 2023; originally announced November 2023.

Comments: Accepted for publication in ApJ. Data available at https://doi.org/10.5281/zenodo.10929036

Galaxies Going Bananas: Inferring the 3D Geometry of High-Redshift Galaxies with JWST-CEERS

Authors: Viraj Pandya, Haowen Zhang, Marc Huertas-Company, Kartheik G. Iyer, Elizabeth McGrath, Guillermo Barro, Steven L. Finkelstein, Martin Kuemmel, William G. Hartley, Henry C. Ferguson, Jeyhan S. Kartaltepe, Joel Primack, Avishai Dekel, Sandra M. Faber, David C. Koo, Greg L. Bryan, Rachel S. Somerville, Ricardo O. Amorin, Pablo Arrabal Haro, Micaela B. Bagley, Eric F. Bell, Emmanuel Bertin, Luca Costantin, Romeel Dave, Mark Dickinson , et al. (31 additional authors not shown)

Abstract: The 3D geometry of high-redshift galaxies remains poorly understood. We build a differentiable Bayesian model and use Hamiltonian Monte Carlo to efficiently and robustly infer the 3D shapes of star-forming galaxies in JWST-CEERS observations with $\log M_*/M_{\odot}=9.0-10.5$ at $z=0.5-8.0$. We reproduce previous results from HST-CANDELS in a fraction of the computing time and constrain the mean e… ▽ More The 3D geometry of high-redshift galaxies remains poorly understood. We build a differentiable Bayesian model and use Hamiltonian Monte Carlo to efficiently and robustly infer the 3D shapes of star-forming galaxies in JWST-CEERS observations with $\log M_*/M_{\odot}=9.0-10.5$ at $z=0.5-8.0$. We reproduce previous results from HST-CANDELS in a fraction of the computing time and constrain the mean ellipticity, triaxiality, size and covariances with samples as small as $\sim50$ galaxies. We find high 3D ellipticities for all mass-redshift bins suggesting oblate (disky) or prolate (elongated) geometries. We break that degeneracy by constraining the mean triaxiality to be $\sim1$ for $\log M_*/M_{\odot}=9.0-9.5$ dwarfs at $z>1$ (favoring the prolate scenario), with significantly lower triaxialities for higher masses and lower redshifts indicating the emergence of disks. The prolate population traces out a ``banana'' in the projected $b/a-\log a$ diagram with an excess of low $b/a$, large $\log a$ galaxies. The dwarf prolate fraction rises from $\sim25\%$ at $z=0.5-1.0$ to $\sim50-80\%$ at $z=3-8$. If these are disks, they cannot be axisymmetric but instead must be unusually oval (triaxial) unlike local circular disks. We simultaneously constrain the 3D size-mass relation and its dependence on 3D geometry. High-probability prolate and oblate candidates show remarkably similar Sérsic indices ($n\sim1$), non-parametric morphological properties and specific star formation rates. Both tend to be visually classified as disks or irregular but edge-on oblate candidates show more dust attenuation. We discuss selection effects, follow-up prospects and theoretical implications. △ Less

Submitted 15 January, 2024; v1 submitted 23 October, 2023; originally announced October 2023.

Comments: Accepted version to appear in ApJ, main body is 36 pages of which ~half are full-page figures

$Σ_{\mathrm{SFR}}$-M* Diagram: A Valuable Galaxy Evolution Diagnostic to Complement (s)SFR-M* Diagrams

Authors: Samir Salim, Sandro Tacchella, Chandler Osborne, S. M. Faber, Janice C. Lee, Sara L. Ellison

Abstract: The specific star formation rate (sSFR) is commonly used to describe the level of galaxy star formation (SF) and to select quenched galaxies. However, being a relative measure of the young-to-old population, an ambiguity in its interpretation may arise because a small sSFR can be either because of a substantial previous mass build up, or because SF is low. We show, using large samples spanning 0 <… ▽ More The specific star formation rate (sSFR) is commonly used to describe the level of galaxy star formation (SF) and to select quenched galaxies. However, being a relative measure of the young-to-old population, an ambiguity in its interpretation may arise because a small sSFR can be either because of a substantial previous mass build up, or because SF is low. We show, using large samples spanning 0 < z < 2, that the normalization of SFR by the physical extent over which SF is taking place (i.e., SFR surface density, $Σ_{\mathrm{SFR}}$) overcomes this ambiguity. $Σ_{\mathrm{SFR}}$ has a strong physical basis, being tied to the molecular gas density and the effectiveness of stellar feedback, so we propose $Σ_{\mathrm{SFR}}$-M* as an important galaxy evolution diagram to complement (s)SFR-M* diagrams. Using the $Σ_{\mathrm{SFR}}$-M* diagram we confirm the Schiminovich et al. (2007) result that the level of SF along the main sequence today is only weakly mass dependent - high-mass galaxies, despite their redder colors, are as active as blue, low-mass ones. At higher redshift, the slope of the "$Σ_{\mathrm{SFR}}$ main sequence" steepens, signaling the epoch of bulge build-up in massive galaxies. We also find that $Σ_{\mathrm{SFR}}$ based on the optical isophotal radius more cleanly selects both the starbursting and the spheroid-dominated (early-type) galaxies than sSFR. One implication of our analysis is that the assessment of the inside-out vs. outside-in quenching scenarios should consider both sSFR and $Σ_{\mathrm{SFR}}$ radial profiles, because ample SF may be present in bulges with low sSFR (red color). △ Less

Submitted 17 October, 2023; originally announced October 2023.

Comments: 16 pages. Accepted to ApJ. Comments on content or relevant missing references welcome

The Evolving Effect Of Cosmic Web Environment On Galaxy Quenching

Authors: Farhanul Hasan, Joseph N. Burchett, Alyssa Abeyta, Douglas Hellinger, Nir Mandelker, Joel R. Primack, S. M. Faber, David C. Koo, Oskar Elek, Daisuke Nagai

Abstract: We investigate how cosmic web structures affect galaxy quenching in the IllustrisTNG (TNG100) cosmological simulations by reconstructing the cosmic web within each snapshot using the DisPerSE framework. We measure the comoving distance from each galaxy with stellar mass $\log(M_{\ast}/\mathrm{M}_{\odot}) \geq 8$ to the nearest node ($d_{\mathrm{node}}$) and the nearest filament spine (… ▽ More We investigate how cosmic web structures affect galaxy quenching in the IllustrisTNG (TNG100) cosmological simulations by reconstructing the cosmic web within each snapshot using the DisPerSE framework. We measure the comoving distance from each galaxy with stellar mass $\log(M_{\ast}/\mathrm{M}_{\odot}) \geq 8$ to the nearest node ($d_{\mathrm{node}}$) and the nearest filament spine ($d_{\mathrm{fil}}$) to study the dependence of both median specific star formation rate (<sSFR>) and median gas fraction (<$f_{\mathrm{gas}}$>) on these distances. We find that the <sSFR> of galaxies is only dependent on cosmic web environment at $z<2$, with the dependence increasing with time. At $z\leq0.5$, $8 \leq \log(M_{\ast}/\mathrm{M}_{\odot}) < 9$ galaxies are quenched at $d_{\mathrm{node}}\lesssim1$~Mpc, and have significantly-suppressed star formation at $d_{\mathrm{fil}}\lesssim1$~Mpc, trends driven mostly by satellite galaxies. At $z\leq1$, in contrast to the monotonic drop in <sSFR> of $\log(M_{\ast}/\mathrm{M}_{\odot}) <10$ galaxies with decreasing $d_{\mathrm{node}}$ and $d_{\mathrm{fil}}$, $\log(M_{\ast}/\mathrm{M}_{\odot}) \geq 10$ galaxies - both centrals and satellites - experience an upturn in <sSFR> at $d_{\mathrm{node}}\lesssim0.2$~Mpc. Much of this cosmic web dependence of star formation activity can be explained by an evolution in $<f_{\mathrm{gas}}>$. Our results suggest that in the past $\sim$10 Gyr, low-mass satellites are quenched by rapid gas stripping in dense environments near nodes and gradual gas starvation in intermediate-density environments near filaments, while at earlier times cosmic web structures efficiently channeled cold gas into most galaxies. State-of-the-art ongoing spectroscopic surveys such as SDSS and DESI, as well as those planned with the Subaru Prime Focus Spectrograph, JWST and Roman, are required to test our predictions against observations. △ Less

Submitted 24 April, 2023; v1 submitted 14 March, 2023; originally announced March 2023.

Comments: Accepted for publication in ApJ. 5 Figures

The Art of Measuring Physical Parameters in Galaxies: A Critical Assessment of Spectral Energy Distribution Fitting Techniques

Authors: Camilla Pacifici, Kartheik G. Iyer, Bahram Mobasher, Elisabete da Cunha, Viviana Acquaviva, Denis Burgarella, Gabriela Calistro Rivera, Adam C. Carnall, Yu-Yen Chang, Nima Chartab, Kevin C. Cooke, Ciaran Fairhurst, Jeyhan Kartaltepe, Joel Leja, Katarzyna Malek, Brett Salmon, Marianna Torelli, Alba Vidal-Garcia, Mederic Boquien, Gabriel G. Brammer, Michael J. I. Brown, Peter L. Capak, Jacopo Chevallard, Chiara Circosta, Darren Croton , et al. (30 additional authors not shown)

Abstract: The study of galaxy evolution hinges on our ability to interpret multi-wavelength galaxy observations in terms of their physical properties. To do this, we rely on spectral energy distribution (SED) models which allow us to infer physical parameters from spectrophotometric data. In recent years, thanks to the wide and deep multi-waveband galaxy surveys, the volume of high quality data have signifi… ▽ More The study of galaxy evolution hinges on our ability to interpret multi-wavelength galaxy observations in terms of their physical properties. To do this, we rely on spectral energy distribution (SED) models which allow us to infer physical parameters from spectrophotometric data. In recent years, thanks to the wide and deep multi-waveband galaxy surveys, the volume of high quality data have significantly increased. Alongside the increased data, algorithms performing SED fitting have improved, including better modeling prescriptions, newer templates, and more extensive sampling in wavelength space. We present a comprehensive analysis of different SED fitting codes including their methods and output with the aim of measuring the uncertainties caused by the modeling assumptions. We apply fourteen of the most commonly used SED fitting codes on samples from the CANDELS photometric catalogs at z~1 and z~3. We find agreement on the stellar mass, while we observe some discrepancies in the star formation rate (SFR) and dust attenuation results. To explore the differences and biases among the codes, we explore the impact of the various modeling assumptions as they are set in the codes (e.g., star formation histories, nebular, dust, and AGN models) on the derived stellar masses, SFRs, and A_V values. We then assess the difference among the codes on the SFR-stellar mass relation and we measure the contribution to the uncertainties by the modeling choices (i.e., the modeling uncertainties) in stellar mass (~0.1dex), SFR (~0.3dex), and dust attenuation (~0.3mag). Finally, we present some resources summarizing best practices in SED fitting. △ Less

Submitted 4 December, 2022; originally announced December 2022.

Comments: 25 pages, 11 figures. Accepted for publication in ApJ

A New Physical Picture for AGNs Lacking Optical Emission Lines

Authors: Christopher J. Agostino, Samir Salim, Sara L. Ellison, Robert W. Bickley, S. M. Faber

Abstract: In this work, we use ~500 low-redshift (z ~ 0.1) X-ray AGNs observed by XMM-Newton and SDSS to investigate the prevalence and nature of AGNs that apparently lack optical emission lines (``optically dull AGNs''). Although 1/4 of spectra appear absorption-line dominated in visual assessment, line extraction with robust continuum subtraction from the MPA/JHU catalog reveals usable [OIII] measurements… ▽ More In this work, we use ~500 low-redshift (z ~ 0.1) X-ray AGNs observed by XMM-Newton and SDSS to investigate the prevalence and nature of AGNs that apparently lack optical emission lines (``optically dull AGNs''). Although 1/4 of spectra appear absorption-line dominated in visual assessment, line extraction with robust continuum subtraction from the MPA/JHU catalog reveals usable [OIII] measurements in 98% of the sample, allowing us to study [OIII]-underluminous AGNs together with more typical AGNs in the context of the L$_{\mathrm{[OIII]}}$--L$_{X}$ relation. We find that ``optically dull AGNs'' do not constitute a distinct population of AGNs. Instead, they are the [OIII]-underluminous tail of a single, unimodal L$_{\mathrm{[OIII]}}$--L$_{X}$ relation that has substantial scatter (0.6 dex). We find the degree to which an AGN is underluminous in [OIII] correlates with the specific SFR or D$_{4000}$ index of the host, which are both linked to the molecular gas fraction. Thus the emerging physical picture for the large scatter seems to involve the gas content of the narrow-line region. We find no significant role for previously proposed scenarios for the presence of optically dull AGNs, such as host dilution or dust obscuration. Despite occasionally weak lines in SDSS spectra, >80% of X-ray AGNs are identified as such with the BPT diagram. >90% are classified as AGNs based only on [NII]/H$α$, providing more complete AGN samples when [OIII] or H$β$ are weak. X-ray AGNs with LINER spectra obey essentially the same \lxo\ relation as Seyfert 2s, suggesting their line emission is produced by AGN activity. △ Less

Submitted 14 November, 2022; originally announced November 2022.

Comments: 21 pages, 14 figures, submitted to ApJ. Comments welcome

The SCUBA-2 Cosmology Legacy Survey: The EGS deep field -- III. The evolution of faint submillimeter galaxies at $z<4$

Authors: L. Cardona-Torres, I. Aretxaga, A. Montaña, J. A. Zavala, S. M. Faber

Abstract: We present a demographic analysis of the physical and morphological properties of $450/850~μ\rm m$-selected galaxies from the deep observations of the SCUBA-2 Cosmology Legacy Survey in the Extended Groth Strip that are detected below the classical submillimeter-galaxy regime ($S_{850 μ\rm m}\lesssim 6~\rm mJy$/beam) and compare them with a sample of optically-selected star-forming galaxies detect… ▽ More We present a demographic analysis of the physical and morphological properties of $450/850~μ\rm m$-selected galaxies from the deep observations of the SCUBA-2 Cosmology Legacy Survey in the Extended Groth Strip that are detected below the classical submillimeter-galaxy regime ($S_{850 μ\rm m}\lesssim 6~\rm mJy$/beam) and compare them with a sample of optically-selected star-forming galaxies detected in the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey in the same field. We derive the evolution of the main sequence of star-forming galaxies, finding a steeper specific star formation rate versus stellar mass at $z>2.5$ than previous studies. Most faint submillimeter-galaxies fall within $3σ$ of the main sequence, but 40~per cent are classified as starbursts. Faint submillimeter galaxies have 50~per cent larger sizes at $2<z<3$ than optically-selected star-forming galaxies of the same mass range. This is also the redshift bin where we find the largest fraction of starbursts, and hence we could be witnessing merging processes, as confirmed by the preference for visual-morphology classifications of these systems as irregular disk galaxies and mergers. Both populations show an increment towards lower redshifts ($z<2$) of their concentration in $H$-band morphology, but faint submillimeter galaxies on average show larger concentration values at later times. These findings support the claim that faint submillimeter galaxies are mostly a population of massive dust-obscured disk-like galaxies that develop larger bulge components at later epochs. While the similarities are great, the median sizes, starburst numbers and $H$-band concentration of faint submillimeter galaxies differ from those of optically-selected star-forming galaxies of the same stellar mass. △ Less

Submitted 10 October, 2022; originally announced October 2022.

Comments: Accepted for publication in MNRAS

A Long Time Ago in a Galaxy Far, Far Away: A Candidate z ~ 12 Galaxy in Early JWST CEERS Imaging

Authors: Steven L. Finkelstein, Micaela B. Bagley, Pablo Arrabal Haro, Mark Dickinson, Henry C. Ferguson, Jeyhan S. Kartaltepe, Casey Papovich, Denis Burgarella, Dale D. Kocevski, Marc Huertas-Company, Kartheik G. Iyer, Rebecca L. Larson, Pablo G. Pérez-González, Caitlin Rose, Sandro Tacchella, Stephen M. Wilkins, Katherine Chworowsky, Aubrey Medrano, Alexa M. Morales, Rachel S. Somerville, L. Y. Aaron Yung, Adriano Fontana, Mauro Giavalisco, Andrea Grazian, Norman A. Grogin , et al. (95 additional authors not shown)

Abstract: We report the discovery of a candidate galaxy with a photo-z of z~12 in the first epoch of the JWST Cosmic Evolution Early Release Science (CEERS) Survey. Following conservative selection criteria we identify a source with a robust z_phot = 11.8^+0.3_-0.2 (1-sigma uncertainty) with m_F200W=27.3, and >7-sigma detections in five filters. The source is not detected at lambda < 1.4um in deep imaging f… ▽ More We report the discovery of a candidate galaxy with a photo-z of z~12 in the first epoch of the JWST Cosmic Evolution Early Release Science (CEERS) Survey. Following conservative selection criteria we identify a source with a robust z_phot = 11.8^+0.3_-0.2 (1-sigma uncertainty) with m_F200W=27.3, and >7-sigma detections in five filters. The source is not detected at lambda < 1.4um in deep imaging from both HST and JWST, and has faint ~3-sigma detections in JWST F150W and HST F160W, which signal a Ly-alpha break near the red edge of both filters, implying z~12. This object (Maisie's Galaxy) exhibits F115W-F200W > 1.9 mag (2-sigma lower limit) with a blue continuum slope, resulting in 99.6% of the photo-z PDF favoring z > 11. All data quality images show no artifacts at the candidate's position, and independent analyses consistently find a strong preference for z > 11. Its colors are inconsistent with Galactic stars, and it is resolved (r_h = 340 +/- 14 pc). Maisie's Galaxy has log M*/Msol ~ 8.5 and is highly star-forming (log sSFR ~ -8.2 yr^-1), with a blue rest-UV color (beta ~ -2.5) indicating little dust though not extremely low metallicity. While the presence of this source is in tension with most predictions, it agrees with empirical extrapolations assuming UV luminosity functions which smoothly decline with increasing redshift. Should followup spectroscopy validate this redshift, our Universe was already aglow with galaxies less than 400 Myr after the Big Bang. △ Less

Submitted 7 September, 2022; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: 17 pages, 6 figures, 3 tables, ApJL in press. Summary of changes from original submission: Improvements in astrometry generated a weak detection in F150W that reduces the photo-z to 11.8 but does not increase the likelihood of lower-z solutions. A full discussion of changes from the original version is available at: https://web.corral.tacc.utexas.edu/ceersdata/papers/Maisie_update.pdf

The Dwarf Galaxy Population at $z\sim 0.7$: A Catalog of Emission Lines and Redshifts from Deep Keck Observations

Authors: John Pharo, Yicheng Guo, Guillermo Barro Calvo, Timothy Carleton, S. M. Faber, Puragra Guhathakurta, Susan A. Kassin, David C. Koo, Jack Lonergan, Teja Teppala, Weichen Wang, Hassen M. Yesuf, Fuyan Bian, Romeel Dave, John C. Forbes, Dusan Keres, Pablo Perez-Gonzalez, Alec Martin, A. J. Puleo, Lauryn Williams, Benjamin Winningham

Abstract: We present a catalog of spectroscopically measured redshifts over $0 < z < 2$ and emission line fluxes for 1440 galaxies. The majority ($\sim$65\%) of the galaxies come from the HALO7D survey, with the remainder from the DEEPwinds program. This catalog includes redshifts for 646 dwarf galaxies with $\log(M_{\star}/M_{\odot}) < 9.5$. 810 catalog galaxies did not have previously published spectrosco… ▽ More We present a catalog of spectroscopically measured redshifts over $0 < z < 2$ and emission line fluxes for 1440 galaxies. The majority ($\sim$65\%) of the galaxies come from the HALO7D survey, with the remainder from the DEEPwinds program. This catalog includes redshifts for 646 dwarf galaxies with $\log(M_{\star}/M_{\odot}) < 9.5$. 810 catalog galaxies did not have previously published spectroscopic redshifts, including 454 dwarf galaxies. HALO7D used the DEIMOS spectrograph on the Keck II telescope to take very deep (up to 32 hours exposure, with a median of $\sim$7 hours) optical spectroscopy in the COSMOS, EGS, GOODS-North, and GOODS-South CANDELS fields, and in some areas outside CANDELS. We compare our redshift results to existing spectroscopic and photometric redshifts in these fields, finding only a 1\% rate of discrepancy with other spectroscopic redshifts. We measure a small increase in median photometric redshift error (from 1.0\% to 1.3\%) and catastrophic outlier rate (from 3.5\% to 8\%) with decreasing stellar mass. We obtained successful redshift fits for 75\% of massive galaxies, and demonstrate a similar 70-75\% successful redshift measurement rate in $8.5 < \log(M_{\star}/M_{\odot}) < 9.5$ galaxies, suggesting similar survey sensitivity in this low-mass range. We describe the redshift, mass, and color-magnitude distributions of the catalog galaxies, finding HALO7D galaxies representative of CANDELS galaxies up to \textit{i}-band magnitudes of 25. The catalogs presented will enable studies of star formation (SF), the mass-metallicity relation, SF-morphology relations, and other properties of the $z\sim0.7$ dwarf galaxy population. △ Less

Submitted 25 July, 2022; v1 submitted 17 March, 2022; originally announced March 2022.

Comments: 23 pages, 19 Figures, updated to version accepted by ApJS

The building up of observed stellar scaling relations of massive galaxies and the connection to black hole growth in the TNG50 simulation

Authors: S. Varma, M. Huertas-Company, A. Pillepich, D. Nelson, V. Rodriguez-Gomez, A. Dekel, S. M. Faber, P. Iglesias-Navarro, D. C. Koo, J. Primack

Abstract: [abridged] We study how mock-observed stellar morphological and structural properties of massive galaxies are built up between $z=0.5$ and $z=3$ in the TNG50 cosmological simulation. We generate mock images with the properties of the CANDELS survey and derive Sersic parameters and optical rest-frame morphologies as usually done in the observations. Overall, the simulation reproduces the observed e… ▽ More [abridged] We study how mock-observed stellar morphological and structural properties of massive galaxies are built up between $z=0.5$ and $z=3$ in the TNG50 cosmological simulation. We generate mock images with the properties of the CANDELS survey and derive Sersic parameters and optical rest-frame morphologies as usually done in the observations. Overall, the simulation reproduces the observed evolution of the abundances of different galaxy morphological types of star-forming and quiescent galaxies. The $\log{M_*}-\log R_e$ and $\log{M_*}-\logΣ_1$ relations of the simulated star-forming and quenched galaxies also match the observed slopes and zeropoints to within 1-$σ$. In the simulation, galaxies increase their observed central stellar mass density ($Σ_1$) and transform in morphology from irregular/clumpy systems to normal Hubble-type systems in the Star Formation Main Sequence at a characteristic stellar mass of $\sim 10^{10.5}~M_\odot$. This morphological transformation is connected to the activity of the central Super Massive Black Holes (SMBHs). At low stellar masses ($10^9$ < $M_*/M_\odot$ < $10^{10}$) SMBHs grow rapidly, while at higher mass SMBHs switch into the kinetic feedback mode and grow more slowly. During this low-accretion phase, SMBH feedback leads to the quenching of star-formation, along with a simultaneous growth in $Σ_1$. More compact massive galaxies grow their SMBHs faster than extended ones of the same mass and end up quenching earlier. In the TNG50 simulation, SMBHs predominantly grow via gas accretion before galaxies quench, and $Σ_1$ increases substantially after SMBH growth slows down. The simulation predicts therefore that quiescent galaxies have higher $Σ_1$ values than star-forming galaxies for the same SMBH mass, which disagrees with alternative models, and may potentially be in tension with some observations. △ Less

Submitted 22 October, 2021; originally announced October 2021.

Comments: accepted for publication in MNRAS

The Baltimore Oriole's Nest: Cool Winds from the Inner and Outer Parts of a Star-Forming Galaxy at z=1.3

Authors: Weichen Wang, Susan A. Kassin, S. M. Faber, David C. Koo, Emily C. Cunningham, Hassen M. Yesuf, Guillermo Barro, Puragra Guhathakurta, Benjamin Weiner, Alexander de la Vega, Yicheng Guo, Timothy M Heckman, Camilla Pacifici, Bingjie Wang, Charlotte Welker

Abstract: Strong galactic winds are ubiquitous at $z\gtrsim 1$. However, it is not well known where inside galaxies these winds are launched from. We study the cool winds ($\sim 10^4$\,K) in two spatial regions of a massive galaxy at $z=1.3$, which we nickname the "Baltimore Oriole's Nest." The galaxy has a stellar mass of $10^{10.3\pm 0.3} M_\odot$, is located on the star-forming main sequence, and has a m… ▽ More Strong galactic winds are ubiquitous at $z\gtrsim 1$. However, it is not well known where inside galaxies these winds are launched from. We study the cool winds ($\sim 10^4$\,K) in two spatial regions of a massive galaxy at $z=1.3$, which we nickname the "Baltimore Oriole's Nest." The galaxy has a stellar mass of $10^{10.3\pm 0.3} M_\odot$, is located on the star-forming main sequence, and has a morphology indicative of a recent merger. Gas kinematics indicate a dynamically complex system with velocity gradients ranging from 0 to 60 $\mathrm{km}\cdot\mathrm{s}^{-1}$. The two regions studied are: a dust-reddened center (Central region), and a blue arc at 7 kpc from the center (Arc region). We measure the \ion{Fe}{2} and \ion{Mg}{2} absorption line profiles from deep Keck/DEIMOS spectra. Blueshifted wings up to 450 km$\cdot$s$^{-1}$ are found for both regions. The \ion{Fe}{2} column densities of winds are $10^{14.7\pm 0.2}\,\mathrm{cm}^{-2}$ and $10^{14.6\pm 0.2}\,\mathrm{cm}^{-2}$ toward the Central and Arc regions, respectively. Our measurements suggest that the winds are most likely launched from both regions. The winds may be driven by the spatially extended star formation, the surface density of which is around 0.2 $M_\odot\,\mathrm{yr}^{-1}\cdot \mathrm{kpc}^{-2}$ in both regions. The mass outflow rates are estimated to be $4\,M_\odot\,\mathrm{yr}^{-1}$ and $3\,M_\odot\,\mathrm{yr}^{-1}$ for the Central and Arc regions, with uncertainties of one order-of-magnitude or more. Findings of this work and a few previous studies suggest that the cool galactic winds at $z\gtrsim 1$ might be commonly launched from the entire spatial extents of their host galaxies due to extended galaxy star formation. △ Less

Submitted 6 April, 2022; v1 submitted 24 September, 2021; originally announced September 2021.

Comments: Accepted to ApJ

What is Important? Morphological Asymmetries are Useful Predictors of Star Formation Rates of Star-forming Galaxies in SDSS Stripe 82

Authors: Hassen M. Yesuf, Luis C. Ho, S. M. Faber

Abstract: Morphology and structure of galaxies reflect their star formation and assembly histories. We use the framework of mutual information ($\mathrm{MI}$) to quantify interdependence among several structural variables and to rank them according to their relevance for predicting specific star formation rate (SSFR) by comparing the $\mathrm{MI}$ of the predictor variables with SSFR and penalizing variable… ▽ More Morphology and structure of galaxies reflect their star formation and assembly histories. We use the framework of mutual information ($\mathrm{MI}$) to quantify interdependence among several structural variables and to rank them according to their relevance for predicting specific star formation rate (SSFR) by comparing the $\mathrm{MI}$ of the predictor variables with SSFR and penalizing variables that are redundant. We apply this framework to study $\sim 3,700$ face-on star-forming galaxies (SFGs) with varying degrees of bulge dominance and central concentration and with stellar mass $M_\star \approx 10^9\,M_\odot - 5\times 10^{11}\,M_\odot$ at redshift $z = 0.02-0.12$. We use the Sloan Digital Sky Survey (SDSS) Stripe 82 deep $i$-band imaging data, which improve measurements of asymmetry and bulge dominance indicators. We find that star-forming galaxies are a multi-parameter family. In addition to $M_\star$, asymmetry emerges as the most powerful predictor of SSFR residuals of SFGs, followed by bulge prominence/concentration. Star-forming galaxies with higher asymmetry and stronger bulges have higher SSFR at a given $M_\star$. The asymmetry reflects both irregular spiral arms and lopsidedness in seemingly isolated SFGs and structural perturbations by galaxy interactions or mergers. △ Less

Submitted 18 September, 2021; originally announced September 2021.

Comments: 28 pages, 14 figures, accepted by ApJ

Physical Drivers of Emission Line Diversity of SDSS Seyfert 2s and LINERs After Removal of Contributions by Star Formation

Authors: Christopher J. Agostino, Samir Salim, S. M. Faber, Stéphanie Juneau, David C. Koo, Yimeng Tang, Yifei Luo, Sofia Quiros, Pin-Song Zhao

Abstract: Ionization sources other than HII regions give rise to the right-hand branch in the standard ([NII]) BPT diagram, populated by Seyfert 2s and LINERs. However, because the majority of Seyfert/LINER hosts are star forming (SF), HII regions contaminate the observed lines to some extent, making it unclear if the position along the branch is merely due to various degrees of mixing between pure Seyfert/… ▽ More Ionization sources other than HII regions give rise to the right-hand branch in the standard ([NII]) BPT diagram, populated by Seyfert 2s and LINERs. However, because the majority of Seyfert/LINER hosts are star forming (SF), HII regions contaminate the observed lines to some extent, making it unclear if the position along the branch is merely due to various degrees of mixing between pure Seyfert/LINER and SF, or whether it reflects the intrinsic diversity of Seyfert/LINER ionizing sources. In this study, we empirically remove SF contributions in ~100,000 Seyfert/LINERs from SDSS using the doppelganger method. We find that mixing is not the principal cause of the extended morphology of the observed branch. Rather, Seyferts/LINERs intrinsically have a wide range of line ratios. Variations in ionization parameter and metallicity can account for much of the diversity of Seyfert/LINER line ratios, but the hardness of ionization field also varies significantly. Furthermore, our k-means classification on seven decontaminated emission lines reveals that LINERs are made up of two populations, which we call soft and hard LINERs. The Seyfert 2s differ from both types of LINERs primarily by higher ionization parameter, whereas the two LINER types mainly differ from each other (and from star-forming regions) in the hardness of the radiation field. We confirm that the [NII] BPT diagram more efficiently identifies LINERs than [SII] and [OI] diagnostics, because in the latter many LINERs, especially soft ones, occupy the same location as pure star-formers, even after the SF has been removed from LINER emission. △ Less

Submitted 17 August, 2021; originally announced August 2021.

Comments: Accepted to ApJ. Main results in Figures 6 and 7. Seyfert/LINER classification (Table 1) available at https://salims.pages.iu.edu/agn . Comments welcome

Implications of Increased Central Mass Surface Densities for the Quenching of Low-mass Galaxies

Authors: Yicheng Guo, Timothy Carleton, Eric F. Bell, Zhu Chen, Avishai Dekel, S. M. Faber, Mauro Giavalisco, Dale D. Kocevski, Anton M. Koekemoer, David C. Koo, Peter Kurczynski, Seong-Kook Lee, F. S. Liu, Casey Papovich, Pablo G. Pérez-González

Abstract: We use the Cosmic Assembly Deep Near-infrared Extragalactic Legacy Survey (CANDELS) data to study the relationship between quenching and the stellar mass surface density within the central radius of 1 kpc ($Σ_1$) of low-mass galaxies (stellar mass $M_* \lesssim 10^{9.5} M_\odot$) at $0.5 \leq z < 1.5$. Our sample is mass complete down to $\sim 10^9 M_\odot$ at $0.5 \leq z < 1.0$. We compare the me… ▽ More We use the Cosmic Assembly Deep Near-infrared Extragalactic Legacy Survey (CANDELS) data to study the relationship between quenching and the stellar mass surface density within the central radius of 1 kpc ($Σ_1$) of low-mass galaxies (stellar mass $M_* \lesssim 10^{9.5} M_\odot$) at $0.5 \leq z < 1.5$. Our sample is mass complete down to $\sim 10^9 M_\odot$ at $0.5 \leq z < 1.0$. We compare the mean $Σ_1$ of star-forming galaxies (SFGs) and quenched galaxies (QGs) at the same redshift and $M_*$. We find that low-mass QGs have higher $Σ_1$ than low-mass SFGs, similar to galaxies above $10^{10} M_\odot$. The difference of $Σ_1$ between QGs and SFGs increases slightly with $M_*$ at $M_* \lesssim 10^{10} M_\odot$ and decreases with $M_*$ at $M_* \gtrsim 10^{10} M_\odot$. The turnover mass is consistent with the mass where quenching mechanisms transition from internal to environmental quenching. At $0.5 \leq z < 1.0$, we find that the $Σ_1$ of galaxies increases by about 0.25 dex in the green valley (i.e., the transitioning region from star forming to fully quenched), regardless of their $M_*$. Using the observed specific star formation rate (sSFR) gradient in the literature as a constraint, we estimate that the quenching timescale (i.e., time spent in the transition) of low-mass galaxies is a few ($\sim4$) Gyrs at $0.5 \leq z < 1.0$. The mechanisms responsible for quenching need to gradually quench star formation in an outside-in way, i.e., preferentially ceasing star formation in outskirts of galaxies while maintaining their central star formation to increase $Σ_1$. An interesting and intriguing result is the similarity of the growth of $Σ_1$ in the green valley between low-mass and massive galaxies, which suggests that the role of internal processes in quenching low-mass galaxies is a question worthy of further investigation. △ Less

Submitted 25 May, 2021; originally announced May 2021.

Comments: 12 pages, 11 figures. Accepted by ApJ

X-ray bubbles in the circumgalactic medium of TNG50 Milky Way- and M31-like galaxies: signposts of supermassive black hole activity

Authors: Annalisa Pillepich, Dylan Nelson, Nhut Truong, Rainer Weinberger, Ignacio Martin-Navarro, Volker Springel, Sandy M. Faber, Lars Hernquist

Abstract: The TNG50 cosmological simulation produces X-ray emitting bubbles, shells, and cavities in the circumgalactic gas above and below the stellar disks of Milky Way- and Andromeda-like galaxies with morphological features reminiscent of the eROSITA and Fermi bubbles in the Galaxy. Two-thirds of the 198 MW/M31 analogues inspected in TNG50 at z=0 show one or more large-scale, coherent features of over-p… ▽ More The TNG50 cosmological simulation produces X-ray emitting bubbles, shells, and cavities in the circumgalactic gas above and below the stellar disks of Milky Way- and Andromeda-like galaxies with morphological features reminiscent of the eROSITA and Fermi bubbles in the Galaxy. Two-thirds of the 198 MW/M31 analogues inspected in TNG50 at z=0 show one or more large-scale, coherent features of over-pressurized gas that impinge into the gaseous halo. Some of the galaxies include a succession of bubbles or shells of increasing size, ranging from a few to many tens of kpc. These are prominent in gas pressure, X-ray emission and gas temperature, and often exhibit sharp boundaries with typical shock Mach numbers of 2-4. The gas in the bubbles outflows with maximum (95th pctl) radial velocities of 100-1500 km/s. TNG50 bubbles expand with speeds as high as 1000-2000 km/s (about 1-2 kpc/Myr), but with a great diversity and with larger bubbles expanding at slower speeds. The bubble gas is at 10^6.4-7.2 K temperatures and is enriched to metallicities of 0.5-2 solar. In TNG50, the bubbles are a manifestation of episodic, kinetic, wind-like energy injections from the supermassive black holes at the galaxy centers that accrete at low Eddington ratios. According to TNG50, X-ray, and possibly gamma-ray, bubbles similar to those observed in the Milky Way should be a frequent feature of disk-like galaxies prior to, or on the verge of, being quenched. They should be within the grasp of eROSITA in the local Universe. △ Less

Submitted 26 October, 2021; v1 submitted 17 May, 2021; originally announced May 2021.

Comments: Accepted on MNRAS. Added details on the SMBH feedback model and a new panel in Fig. 11 with the X-ray luminosities of TNG50 SMBHs. Updated discussion and observability opportunities

A Complete 16 micron-Selected Galaxy Sample at $z\sim1$: Mid-infrared Spectral Energy Distributions

Authors: J. -S. Huang, Y. -S. Dai, S. P. Willner, S. M. Faber, C. Cheng, H. Xu, S. Wu, X. Shao, C. Hao, X. Xia, D. Rigopoulou, M. Pereira Santaella, G. Magdis, I. Cortzen, H. Yan, G. Fazio, P. Assmann, N. Araneda, L. Fan, M. Musin, Z. Wang, K. C. Xu, C. He, A. Esamdin

Abstract: We describe a complete, flux-density-limited sample of galaxies at redshift $0.8 < z < 1.3$ selected at 16 micron. At the selection wavelength near 8 micron rest, the observed emission comes both from dust heated by intense star formation and from active galactic nuclei (AGNs). Fitting the spectral energy distributions (SEDs) of the sample galaxies to local-galaxy templates reveals that more than… ▽ More We describe a complete, flux-density-limited sample of galaxies at redshift $0.8 < z < 1.3$ selected at 16 micron. At the selection wavelength near 8 micron rest, the observed emission comes both from dust heated by intense star formation and from active galactic nuclei (AGNs). Fitting the spectral energy distributions (SEDs) of the sample galaxies to local-galaxy templates reveals that more than half the galaxies have SEDs dominated by star formation. About one sixth of the galaxy SEDs are dominated by an AGN, and nearly all the rest of the SEDs are composite. Comparison with X-ray and far-infrared observations shows that combinations of luminosities at rest-frame 4.5 and 8 micron give good measures of both AGN luminosity and star-formation rate. The sample galaxies mostly follow the established star-forming main sequence for $z=1$ galaxies, but of the galaxies more than 0.5 dex above that main sequence, more than half have AGN-type SEDs. Similarly, the most luminous AGNs tend to have higher star-formation rates than the main sequence value. Galaxies with stellar masses $>$10$^{11}$\,\Msun\ are unlikely to host an AGN. About 1% of the sample galaxies show an SED with dust emission typical of neither star formation nor an AGN. △ Less

Submitted 28 April, 2021; v1 submitted 8 March, 2021; originally announced March 2021.

Comments: 14 pages, 5 tables, 29 figures, ApJ accepted v2 corrects author name formatting

Fast, Slow, Early, Late: Quenching Massive Galaxies at z~0.8

Authors: Sandro Tacchella, Charlie Conroy, S. M. Faber, Benjamin D. Johnson, Joel Leja, Guillermo Barro, Emily C. Cunningham, Alis J. Deason, Puragra Guhathakurta, Yicheng Guo, Lars Hernquist, David C. Koo, Kevin McKinnon, Constance M. Rockosi, Joshua S. Speagle, Pieter van Dokkum, Hassen M. Yesuf

Abstract: We investigate the stellar populations for a sample of 161 massive, mainly quiescent galaxies at $\langle z_{\rm obs} \rangle=0.8$ with deep Keck/DEIMOS rest-frame optical spectroscopy (HALO7D survey). With the fully Bayesian framework Prospector, we simultaneously fit the spectroscopic and photometric data with an advanced physical model (including non-parametric star-formation histories, emissio… ▽ More We investigate the stellar populations for a sample of 161 massive, mainly quiescent galaxies at $\langle z_{\rm obs} \rangle=0.8$ with deep Keck/DEIMOS rest-frame optical spectroscopy (HALO7D survey). With the fully Bayesian framework Prospector, we simultaneously fit the spectroscopic and photometric data with an advanced physical model (including non-parametric star-formation histories, emission lines, variable dust attenuation law, and dust and AGN emission) together with an uncertainty and outlier model. We show that both spectroscopy and photometry are needed to break the dust-age-metallicity degeneracy. We find a large diversity of star-formation histories: although the most massive ($M_{\star}>2\times10^{11}~M_{\odot}$) galaxies formed the earliest (formation redshift of $z_{\rm f}\approx5-10$ with a short star-formation timescale of $τ_{\rm SF}\lesssim1~\mathrm{Gyr}$), lower-mass galaxies have a wide range of formation redshifts, leading to only a weak trend of $z_{\rm f}$ with $M_{\star}$. Interestingly, several low-mass galaxies with have formation redshifts of $z_{\rm f}\approx5-8$. Star-forming galaxies evolve about the star-forming main sequence, crossing the ridgeline several times in their past. Quiescent galaxies show a wide range and continuous distribution of quenching timescales ($τ_{\rm quench}\approx0-5~\mathrm{Gyr}$) with a median of $\langleτ_{\rm quench}\rangle=1.0_{-0.9}^{+0.8}~\mathrm{Gyr}$ and of quenching epochs of $z_{\rm quench}\approx0.8-5.0$ ($\langle z_{\rm quench}\rangle=1.3_{-0.4}^{+0.7}$). This large diversity of quenching timescales and epochs points toward a combination of internal and external quenching mechanisms. In our sample, rejuvenation and "late bloomers" are uncommon. In summary, our analysis supports the "grow & quench" framework and is consistent with a wide and continuously-populated diversity of quenching timescales. △ Less

Submitted 5 January, 2022; v1 submitted 24 February, 2021; originally announced February 2021.

Comments: 30 pages, 12 figures + appendix, accepted by ApJ

Mock Lightcones and Theory Friendly Catalogs for the CANDELS Survey

Authors: Rachel S. Somerville, Charlotte Olsen, L. Y. Aaron Yung, Camilla Pacifici, Henry C. Ferguson, Peter Behroozi, Shannon Osborne, Risa H. Wechsler, Viraj Pandya, Sandra M. Faber, Joel R. Primack, Avishai Dekel

Abstract: We present mock catalogs created to support the interpretation of the CANDELS survey. We extract halos along past lightcones from the Bolshoi Planck dissipationless N-body simulations and populate these halos with galaxies using two different independently developed semi-analytic models of galaxy formation and the empirical model UniverseMachine. Our mock catalogs have geometries that encompass th… ▽ More We present mock catalogs created to support the interpretation of the CANDELS survey. We extract halos along past lightcones from the Bolshoi Planck dissipationless N-body simulations and populate these halos with galaxies using two different independently developed semi-analytic models of galaxy formation and the empirical model UniverseMachine. Our mock catalogs have geometries that encompass the footprints of observations associated with the five CANDELS fields. In order to allow field-to-field variance to be explored, we have created eight realizations of each field. In this paper, we present comparisons with observable global galaxy properties, including counts in observed frame bands, luminosity functions, color-magnitude distributions and color-color distributions. We additionally present comparisons with physical galaxy parameters derived from SED fitting for the CANDELS observations, such as stellar masses and star formation rates. We find relatively good agreement between the model predictions and CANDELS observations for luminosity and stellar mass functions. We find poorer agreement for colors and star formation rate distributions. All of the mock lightcones as well as curated "theory friendly" versions of the observational CANDELS catalogs are made available through a web-based data hub. △ Less

Submitted 29 January, 2021; originally announced February 2021.

Comments: 22 pages, 16 figures. To appear in MNRAS

Radiogenic Heating and its Influence on Rocky Planet Dynamos and Habitability

Authors: Francis Nimmo, Joel Primack, Sandra M. Faber, Enrico Ramirez-Ruiz, Mohammadtaher Safarzadeh

Abstract: The thermal evolution of rocky planets on geological timescales (Gyr) depends on the heat input from the long-lived radiogenic elements potassium, thorium, and uranium. Concentrations of the latter two in rocky planet mantles are likely to vary by up to an order of magnitude between different planetary systems because Th and U, like other heavy r-process elements, are produced by rare stellar proc… ▽ More The thermal evolution of rocky planets on geological timescales (Gyr) depends on the heat input from the long-lived radiogenic elements potassium, thorium, and uranium. Concentrations of the latter two in rocky planet mantles are likely to vary by up to an order of magnitude between different planetary systems because Th and U, like other heavy r-process elements, are produced by rare stellar processes. Here we discuss the effects of these variations on the thermal evolution of an Earth-size planet, using a 1D parameterized convection model. Assuming Th and U abundances consistent with geochemical models of the Bulk Silicate Earth based on chondritic meteorites, we find that Earth had just enough radiogenic heating to maintain a persistent dynamo. According to this model, Earth-like planets of stars with higher abundances of heavy r-process elements, indicated by the relative abundance of europium in their spectra, are likely to have lacked a dynamo for a significant fraction of their lifetimes, with potentially negative consequences for hosting a biosphere. Because the qualitative outcomes of our 1D model are strongly dependent on the treatment of viscosity, further investigations using fully 3D convection models are desirable. △ Less

Submitted 9 November, 2020; originally announced November 2020.

Journal ref: ApJL, 903, L37, 2020

CANDELS Meets GSWLC: Evolution of the Relationship Between Morphology and Star Formation Since z = 2

Authors: Chandler Osborne, Samir Salim, Ivana Damjanov, S. M. Faber, Marc Huertas-Company, David C. Koo, Kameswara Bharadwaj Mantha, Daniel H. McIntosh, Joel R. Primack, Sandro Tacchella

Abstract: Galaxy morphology and its evolution over the cosmic epoch hold important clues for understanding the regulation of star formation (SF). However, studying the relationship between morphology and SF has been hindered by the availability of consistent data at different redshifts. Our sample, combining CANDELS (0.8 < z < 2.5) and the GALEX-SDSS-WISE Legacy Catalog (GSWLC; z ~ 0), has physical paramete… ▽ More Galaxy morphology and its evolution over the cosmic epoch hold important clues for understanding the regulation of star formation (SF). However, studying the relationship between morphology and SF has been hindered by the availability of consistent data at different redshifts. Our sample, combining CANDELS (0.8 < z < 2.5) and the GALEX-SDSS-WISE Legacy Catalog (GSWLC; z ~ 0), has physical parameters derived using consistent SED fitting with flexible dust attenuation laws. We adopt visual classifications from Kartaltepe et al. 2015 and expand them to z ~ 0 using SDSS images matching the physical resolution of CANDELS rest-frame optical images and deep FUV GALEX images matching the physical resolution of the CANDELS rest-frame FUV images. Our main finding is that disks with SF clumps at z ~ 0 make a similar fraction (~15%) of star-forming galaxies as at z ~ 2. The clumpy disk contribution to the SF budget peaks at z ~ 1, rather than z ~ 2, suggesting that the principal epoch of disk assembly continues to lower redshifts. Star-forming spheroids ("blue nuggets"), though less centrally concentrated than quenched spheroids, contribute significantly (~15%) to the SF budget at z ~ 1-2, suggesting that compaction precedes quenching. Among green valley and quiescent galaxies, the pure spheroid fraction drops since z ~ 1, whereas spheroids with disks (S0-like) become dominant. Mergers at or nearing coalescence are enhanced in SFR relative to the main sequence at all redshifts by a factor of ~2, but contribute $\lesssim$5% to the SF budget, with their contribution remaining small above the main sequence. △ Less

Submitted 2 September, 2020; originally announced September 2020.

Comments: 26 pages, 13 figures, 3 tables, accepted for publication in ApJ

Journal ref: The Astrophysical Journal, Volume 902, 77, October 14 2020

The star-forming main sequence and the contribution of dust-obscured star formation since $z\sim4$ from the FUV+IR luminosity functions

Authors: Aldo Rodriguez-Puebla, Vladimir Avila-Reese, Mariana Cano-Diaz, S. M. Faber, Joel R. Primack, Jose Franco, I. Aretxaga, Eder Santiago-Mayoral

Abstract: An analytical approach is proposed to study the evolution of the star-forming galaxy (SFG) main sequence (MS) and the fraction of dust-obscured SF up to $z\sim4$. Far-ultraviolet (FUV) and infrared (IR) star formation rates, SFRs, are described as conditional probability functions of $M_{\ast}$. We convolve them with the galaxy stellar mass function (GSMF) of SFGs to derive the FUV and IR LFs. The… ▽ More An analytical approach is proposed to study the evolution of the star-forming galaxy (SFG) main sequence (MS) and the fraction of dust-obscured SF up to $z\sim4$. Far-ultraviolet (FUV) and infrared (IR) star formation rates, SFRs, are described as conditional probability functions of $M_{\ast}$. We convolve them with the galaxy stellar mass function (GSMF) of SFGs to derive the FUV and IR LFs. The 2 SF modes formalism is used to describe starburst galaxies. By fitting observed FUV and IR LFs, the parametrization of SFR$_{\rm FUV}-M_{\ast}$ and SFR$_{\rm IR}-M_{\ast}$ are constrained. Our derived SFR$_{\rm FUV+IR}-M_{\ast}$ reproduces the evolution of the MS as compared to other observational inferences. At any redshift, we find that the sSFR$_{\rm FUV+IR}-M_{\ast}$ relation for MS SFGs approaches to a power law at the high-mass end. At lower masses, it bends and eventually the slope sign changes from negative to positive at very low masses. At $z\sim0$, this change of sign is at $M_{\ast}\sim5\times10^{8}{\rm M}_{\odot}$ close to dust-obscured SF regime, $M_{\ast}\sim6\times10^{8}{\rm M}_{\odot}$. The slope sign change is related to the knee of the FUV LF. Our derived dust-obscured fractions agree with previous determinations at $0\leq z\leq2.5$. Dust-obscured fractions depend strongly on mass with almost no dependence with redshift at $z\gtrsim1.2$. At $z\lesssim0.75$ high-mass galaxies become more "transparent" compared to their high redshift counterparts. On the opposite, low- and intermediate-mass galaxies have become more obscured by dust. The joint evolution of the GSMF and the FUV and IR LFs is a promising approach to study mass growth and dust formation/destruction mechanisms. △ Less

Submitted 21 October, 2020; v1 submitted 14 August, 2020; originally announced August 2020.

Comments: Submitted to ApJ, 20 pages, 12 figures. Comments are welcome. This is an improved version after the first Referee report. Results remained unchanged

Stellar Masses of Giant Clumps in CANDELS and Simulated Galaxies Using Machine Learning

Authors: M. Huertas-Company, Y. Guo, O. Ginzburg, C. T. Lee, N. Mandelker, M. Metter, J. R. Primack, A. Dekel, D. Ceverino, S. M. Faber, D. C. Koo, A. Koekemoer, G. Snyder, M. Giavalisco, H. Zhang

Abstract: A significant fraction of high redshift star-forming disc galaxies are known to host giant clumps, whose nature and role in galaxy evolution are yet to be understood. In this work we first present a new method based on neural networks to detect clumps in galaxy images. We use this method to detect clumps in the rest-frame optical and UV images of a complete sample of $\sim1500$ star forming galaxi… ▽ More A significant fraction of high redshift star-forming disc galaxies are known to host giant clumps, whose nature and role in galaxy evolution are yet to be understood. In this work we first present a new method based on neural networks to detect clumps in galaxy images. We use this method to detect clumps in the rest-frame optical and UV images of a complete sample of $\sim1500$ star forming galaxies at $1<z<3$ in the CANDELS survey as well as in images from the VELA zoom-in cosmological simulations. We show that observational effects have a dramatic impact on the derived clump properties leading to an overestimation of the clump mass up to a factor of 10, which highlights the importance of fair comparisons between observations and simulations and the limitations of current HST data to study the resolved structure of distant galaxies. After correcting for these effects with a mixture density network, we estimate that the clump stellar mass function follows a power-law down to the completeness limit ($10^{7}$ solar masses) with the majority of the clumps being less massive than $10^9$ solar masses. This is in better agreement with recent gravitational lensing based measurements. The simulations explored in this work overall reproduce the shape of the observed clump stellar mass function and clumpy fractions when confronted under the same conditions, although they tend to lie in the lower limit of the confidence intervals of the observations. This agreement suggests that most of the observed clumps are formed in-situ. △ Less

Submitted 8 September, 2020; v1 submitted 25 June, 2020; originally announced June 2020.

Comments: Accepted for publication in MNRAS - This is the final version

Large Scale Structures in the CANDELS Fields: The Role of the Environment in Star Formation Activity

Authors: Nima Chartab, Bahram Mobasher, Behnam Darvish, Steven L. Finkelstein, Yicheng Guo, Dritan Kodra, Kyoung-Soo Lee, Jeffrey A. Newman, Camilla Pacifici, Casey Papovich, Zahra Sattari, Abtin Shahidi, Mark E. Dickinson, Sandra M. Faber, Henry C. Ferguson, Mauro Giavalisco, Marziye jafariyazani

Abstract: We present a robust method, weighted von Mises kernel density estimation, along with boundary correction to reconstruct the underlying number density field of galaxies. We apply this method to galaxies brighter than $\rm HST/F160w\le 26$ AB mag at the redshift range of $0.4\leq z \leq 5$ in the five CANDELS fields (GOODS-N, GOODS-S, EGS, UDS, and COSMOS). We then use these measurements to explore… ▽ More We present a robust method, weighted von Mises kernel density estimation, along with boundary correction to reconstruct the underlying number density field of galaxies. We apply this method to galaxies brighter than $\rm HST/F160w\le 26$ AB mag at the redshift range of $0.4\leq z \leq 5$ in the five CANDELS fields (GOODS-N, GOODS-S, EGS, UDS, and COSMOS). We then use these measurements to explore the environmental dependence of the star formation activity of galaxies. We find strong evidence of environmental quenching for massive galaxies ($\rm M \gtrsim 10^{11} \rm {M}_\odot$) out to $z\sim 3.5$ such that an over-dense environment hosts $\gtrsim 20\%$ more massive quiescent galaxies compared to an under-dense region. We also find that environmental quenching efficiency grows with stellar mass and reaches $\sim 60\%$ for massive galaxies at $z\sim 0.5$. The environmental quenching is also more efficient in comparison to the stellar mass quenching for low mass galaxies ($\rm M \lesssim 10^{10} \rm {M}_\odot$) at low and intermediate redshifts ($z\lesssim 1.2$). Our findings concur thoroughly with the "over-consumption" quenching model where the termination of cool gas accretion (cosmological starvation) happens in an over-dense environment and the galaxy starts to consume its remaining gas reservoir in depletion time. The depletion time depends on the stellar mass and could explain the evolution of environmental quenching efficiency with the stellar mass. △ Less

Submitted 10 December, 2019; originally announced December 2019.

Comments: 17 pages, 10 figures, 3 tables, accepted for publication in the ApJ. Catalogs of the environment measurements and the density maps will be available with the final published paper in ApJ

The Activation of Galactic Nuclei and Their Accretion Rates are Linked to the Star Formation Rates and Bulge-types of Their Host Galaxies

Authors: Hassen M. Yesuf, S. M. Faber, David C. Koo, Joanna Woo, Joel R. Primack, Yifei Luo

Abstract: We use bulge-type classifications of 809 representative SDSS galaxies by Gadotti (2009) to classify a large sample of galaxies into real bulges (classical or elliptical) and pseudobulges using Random Forest. We use structural and stellar population predictors that can easily be measured without image decomposition. Multiple parameters such as the central mass density with 1 kpc, concentration inde… ▽ More We use bulge-type classifications of 809 representative SDSS galaxies by Gadotti (2009) to classify a large sample of galaxies into real bulges (classical or elliptical) and pseudobulges using Random Forest. We use structural and stellar population predictors that can easily be measured without image decomposition. Multiple parameters such as the central mass density with 1 kpc, concentration index, Sérsic index and velocity dispersion do result in accurate bulge classifications when combined together. We classify $\sim 44,500$ face-on galaxies above stellar mass of 10$^{10}$ M$_\odot$ and redshift $ 0.02 < z < 0.07$ into real bulges or pseudobulges with $93 \pm 2$\% accuracy. We show that $\sim 75 - 90\%$ of AGNs identified by the optical line ratio diagnostic are hosted by real bulges. The pseudobulge fraction significantly decreases with AGN signature as the line ratios change from indicating pure star formation ($\sim 54 \pm 4$ \%), to composite of star formation and AGN ($\sim 18 \pm 3$\%), and to AGN-dominated galaxies ($\sim 5 \pm 3$\%). Using the dust-corrected [\ion{O}{3}] luminosity as an AGN accretion indicator, and the stellar mass and radius as proxies for a black hole mass, we find that AGNs in real bulges have lower Eddington ratios than AGNs in pseudobulges. Real bulges have a wide range of AGN and star formation activities, although most of them are weak AGNs. For both bulge-types, their Eddington ratios are correlated with specific star formation rates (SSFR). Real bulges have lower specific accretion rate but higher AGN fraction than pseudobulges do at similar SSFRs. △ Less

Submitted 8 December, 2019; originally announced December 2019.

Comments: 22 pages, 14 figures, accepted by ApJ

The SFR-radius connection: data and implications for wind strength and halo concentration

Authors: Lin Lin, S. M. Faber, David C. Koo, Samir Salim, Aaron A. Dutton, Jerome J. Fang, Fangzhou Jiang, Cristoph T. Lee, Aldo Rodríguez-Puebla, A. van der Wel, Yicheng Guo, Guillermo Barro, Joel R. Primack, Avishai Dekel, Zhu Chen, Yifei Luo, Viraj Pandya, Rachel S. Somerville, Henry C. Ferguson, Susan Kassin, Anton M. Koekemoer, Norman A. Grogin, Audrey Galametz, P. Santini, Hooshang Nayyeri , et al. (4 additional authors not shown)

Abstract: This paper is one in a series that explores the importance of radius as a second parameter in galaxy evolution. The topic investigated here is the relationship between star formation rate (SFR) and galaxy radius ($R_{\rm e}$) for main-sequence star-forming galaxies. The key observational result is that, over a wide range of stellar mass and redshift in both CANDELS and SDSS, there is little trend… ▽ More This paper is one in a series that explores the importance of radius as a second parameter in galaxy evolution. The topic investigated here is the relationship between star formation rate (SFR) and galaxy radius ($R_{\rm e}$) for main-sequence star-forming galaxies. The key observational result is that, over a wide range of stellar mass and redshift in both CANDELS and SDSS, there is little trend between SFR and $R_{\rm e}$ at fixed stellar mass. The Kennicutt-Schmidt law, or any similar density-related star formation law, then implies that smaller galaxies must have lower gas fractions than larger galaxies (at fixed $M_{\ast}$), and this is supported by observations of local star-forming galaxies. We investigate the implication by adopting the equilibrium "bathtub" model: the ISM gas mass is assumed to be constant over time and the net star formation rate is the difference between the accretion rate of gas onto the galaxy from the halo and the outflow rate due to winds. To match the observed null correlation between SFR and radius, the bathtub model requires that smaller galaxies at fixed mass have weaker galactic winds. Our hypothesis is that galaxies are a 2-dimensional family whose properties are set mainly by halo mass and concentration. Galaxy radius and accretion rate plausibly both depend on halo concentration, which predicts how wind strength should vary with $R_{\rm e}$ and SFR. △ Less

Submitted 20 July, 2020; v1 submitted 24 October, 2019; originally announced October 2019.

Comments: 21 pages, 9 figures, ApJ accepted

Quenching as a Contest between Galaxy Halos and their Central Black Holes

Authors: Zhu Chen, S. M. Faber, David. C. Koo, Rachel S. Somerville, Joel R. Primack, Avishai Dekel, Aldo Rodríguez-Puebla, Yicheng Guo, Guillermo Barro, Dale D. Kocevski, A. van der Wel, Joanna Woo, Eric F. Bell, Jerome J. Fang, Henry C. Ferguson, Mauro Giavalisco, Marc Huertas-Company, Fangzhou Jiang, Susan Kassin, Lin Lin, Fengshan Liu, Yifei Luo, Zhijian Luo, Camilla Pacifici, Viraj Pandya , et al. (5 additional authors not shown)

Abstract: Existing models of galaxy formation have not yet explained striking correlations between structure and star-formation activity in galaxies, notably the sloped and moving boundaries that divide star-forming from quenched galaxies in key structural diagrams. This paper uses these and other relations to ``reverse-engineer'' the quenching process for central galaxies. The basic idea is that star-formi… ▽ More Existing models of galaxy formation have not yet explained striking correlations between structure and star-formation activity in galaxies, notably the sloped and moving boundaries that divide star-forming from quenched galaxies in key structural diagrams. This paper uses these and other relations to ``reverse-engineer'' the quenching process for central galaxies. The basic idea is that star-forming galaxies with larger radii (at a given stellar mass) have lower black-hole masses due to lower central densities. Galaxies cross into the green valley when the cumulative effective energy radiated by their black hole equals $\sim4\times$ their halo-gas binding energy. Since larger-radii galaxies have smaller black holes, one finds they must evolve to higher stellar masses in order to meet this halo-energy criterion, which explains the sloping boundaries. A possible cause of radii differences among star-forming galaxies is halo concentration. The evolutionary tracks of star-forming galaxies are nearly parallel to the green-valley boundaries, and it is mainly the sideways motions of these boundaries with cosmic time that cause galaxies to quench. BH-scaling laws for star-forming, quenched, and green-valley galaxies are different, and most BH mass growth takes place in the green valley. Implications include: the radii of star-forming galaxies are an important second parameter in shaping their black holes; black holes are connected to their halos but in different ways for star-forming, quenched, and green-valley galaxies; and the same BH-halo quenching mechanism has been in place since $z \sim 3$. We conclude with a discussion of black hole-galaxy co-evolution, the origin and interpretation of BH scaling laws. △ Less

Submitted 26 May, 2020; v1 submitted 24 September, 2019; originally announced September 2019.

Comments: 45 pages, 20 figures, 2 tables, accepted for publication in ApJ

Structural and Stellar Population Properties vs. Bulge Types in Sloan Digital Sky Survey Central Galaxies

Authors: Yifei Luo, S. M. Faber, Aldo Rodriguez-Puebla, Joanna Woo, Yicheng Guo, David C. Koo, Joel R. Primack, Zhu Chen, Hassen M. Yesuf, Lin Lin, Guillermo Barro, Jerome J. Fang, Viraj Pandya, M. Huertas-Company, Shude Mao

Abstract: This paper studies pseudo-bulges (P-bulges) and classical bulges (C-bulges) in Sloan Digital Sky Survey central galaxies using the new bulge indicator $ΔΣ_1$, which measures relative central stellar-mass surface density within 1 kpc. We compare $ΔΣ_1$ to the established bulge-type indicator $Δ\langleμ_e\rangle$ from Gadotti (2009) and show that classifying by $ΔΣ_1$ agrees well with… ▽ More This paper studies pseudo-bulges (P-bulges) and classical bulges (C-bulges) in Sloan Digital Sky Survey central galaxies using the new bulge indicator $ΔΣ_1$, which measures relative central stellar-mass surface density within 1 kpc. We compare $ΔΣ_1$ to the established bulge-type indicator $Δ\langleμ_e\rangle$ from Gadotti (2009) and show that classifying by $ΔΣ_1$ agrees well with $Δ\langleμ_e\rangle$. $ΔΣ_1$ requires no bulge-disk decomposition and can be measured on SDSS images out to $z = 0.07$. Bulge types using it are mapped onto twenty different structural and stellar-population properties for 12,000 SDSS central galaxies with masses 10.0 < log $M_*$/$M_{\odot}$ < 10.4. New trends emerge from this large sample. Structural parameters show fairly linear log-log relations vs. $ΔΣ_1$ and $Δ\langleμ_e\rangle$ with only moderate scatter, while stellar-population parameters show a highly non-linear "elbow" in which specific star-formation rate remains roughly flat with increasing central density and then falls rapidly at the elbow, where galaxies begin to quench. P-bulges occupy the low-density end of the horizontal arm of the elbow and are universally star-forming, while C-bulges occupy the elbow and the vertical branch and exhibit a wide range of star-formation rates at fixed density. The non-linear relation between central density and star-formation rate has been seen before, but this mapping onto bulge class is new. The wide range of star-formation rates in C-bulges helps to explain why bulge classifications using different parameters have sometimes disagreed in the past. The elbow-shaped relation between density and stellar indices suggests that central structure and stellar-populations evolve at different rates as galaxies begin to quench. △ Less

Submitted 9 March, 2020; v1 submitted 21 August, 2019; originally announced August 2019.

Comments: 23 pages, 17 figures, published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 493, Issue 2, 2020, p.1686-1707

The CANDELS/SHARDS multi-wavelength catalog in GOODS-N: Photometry, Photometric Redshifts, Stellar Masses, Emission line fluxes and Star Formation Rates

Authors: Guillermo Barro, Pablo G. Perez-Gonzalez, Antonio Cava, Gabriel Brammer, Viraj Pandya, Carmen Eliche Moral, Pilar Esquej, Helena Dominguez-Sanchez, Belen Alcalde Pampliega, Yicheng Guo, Anton M. Koekemoer, Jonathan R. Trump, Matthew L. N. Ashby, Nicolas Cardiel, Marco Castellano, Christopher J. Conselice, Mark E. Dickinson, Timothy Dolch, Jennifer L. Donley, Nestor Espino Briones, Sandra M. Faber, Giovanni G. Fazio, Henry Ferguson, Steve Finkelstein, Adriano Fontana , et al. (30 additional authors not shown)

Abstract: We present a WFC3 F160W ($H$-band) selected catalog in the CANDELS/GOODS-N field containing photometry from the ultraviolet (UV) to the far-infrared (IR), photometric redshifts and stellar parameters derived from the analysis of the multi-wavelength data. The catalog contains 35,445 sources over the 171 arcmin$^{2}$ of the CANDELS F160W mosaic. The 5$σ$ detection limits (within an aperture of radi… ▽ More We present a WFC3 F160W ($H$-band) selected catalog in the CANDELS/GOODS-N field containing photometry from the ultraviolet (UV) to the far-infrared (IR), photometric redshifts and stellar parameters derived from the analysis of the multi-wavelength data. The catalog contains 35,445 sources over the 171 arcmin$^{2}$ of the CANDELS F160W mosaic. The 5$σ$ detection limits (within an aperture of radius 0\farcs17) of the mosaic range between $H=27.8$, 28.2 and 28.7 in the wide, intermediate and deep regions, that span approximately 50\%, 15\% and 35\% of the total area. The multi-wavelength photometry includes broad-band data from UV (U band from KPNO and LBC), optical (HST/ACS F435W, F606W, F775W, F814W, and F850LP), near-to-mid IR (HST/WFC3 F105W, F125W, F140W and F160W, Subaru/MOIRCS Ks, CFHT/Megacam K, and \spitzer/IRAC 3.6, 4.5, 5.8, 8.0 $μ$m) and far IR (\spitzer/MIPS 24$μ$m, HERSCHEL/PACS 100 and 160$μ$m, SPIRE 250, 350 and 500$μ$m) observations. In addition, the catalog also includes optical medium-band data (R$\sim50$) in 25 consecutive bands, $λ=500$ to 950~nm, from the SHARDS survey and WFC3 IR spectroscopic observations with the G102 and G141 grisms (R$\sim210$ and 130). The use of higher spectral resolution data to estimate photometric redshifts provides very high, and nearly uniform, precision from $z=0-2.5$. The comparison to 1,485 good quality spectroscopic redshifts up to $z\sim3$ yields $Δz$/(1+$z_{\rm spec}$)$=$0.0032 and an outlier fraction of $η=$4.3\%. In addition to the multi-band photometry, we release added-value catalogs with emission line fluxes, stellar masses, dust attenuations, UV- and IR- based star formation rates and rest-frame colors. △ Less

Submitted 1 August, 2019; originally announced August 2019.

Comments: The catalogs, images and spectra are available in Rainbow-slicer (https://bit.ly/2OxKKx1 ), navigator (https://bit.ly/2GDS180 ) and MAST (https://bit.ly/2YtoBQ4 ). In addition to the photometry and other added-value catalogs we release UV+IR star formation rates based on Spitzer (MIPS) and Herschel (PACS and SPIRE) in the 5 CANDELS fields

Journal ref: The Astrophysical Journal Supplement Series, Volume 243, Number 2, 2019

Deep ugrizY Imaging and DEEP2/3 Spectroscopy: A Photometric Redshift Testbed for LSST and Public Release of Data from the DEEP3 Galaxy Redshift Survey

Authors: Rongpu Zhou, Michael C. Cooper, Jeffrey A. Newman, Matthew L. N. Ashby, James Aird, Christopher J. Conselice, Marc Davis, Aaron A. Dutton, S. M. Faber, Jerome J. Fang, G. G. Fazio, Puragra Guhathakurta, Dale Kocevski, David C. Koo, Kirpal Nandra, Andrew C. Phillips, David J. Rosario, Edward F. Schlafly, Jonathan R. Trump, Benjamin Weiner, Christopher N. A. Willmer, Renbin Yan

Abstract: We present catalogs of calibrated photometry and spectroscopic redshifts in the Extended Groth Strip, intended for studies of photometric redshifts (photo-z's). The data includes ugriz photometry from CFHTLS and Y-band photometry from the Subaru Suprime camera, as well as spectroscopic redshifts from the DEEP2, DEEP3 and 3D-HST surveys. These catalogs incorporate corrections to produce effectively… ▽ More We present catalogs of calibrated photometry and spectroscopic redshifts in the Extended Groth Strip, intended for studies of photometric redshifts (photo-z's). The data includes ugriz photometry from CFHTLS and Y-band photometry from the Subaru Suprime camera, as well as spectroscopic redshifts from the DEEP2, DEEP3 and 3D-HST surveys. These catalogs incorporate corrections to produce effectively matched-aperture photometry across all bands, based upon object size information available in the catalog and Moffat profile point spread function fits. We test this catalog with a simple machine learning-based photometric redshift algorithm based upon Random Forest regression, and find that the corrected aperture photometry leads to significant improvement in photo-z accuracy compared to the original SExtractor catalogs from CFHTLS and Subaru. The deep ugrizY photometry and spectroscopic redshifts are well-suited for empirical tests of photometric redshift algorithms for LSST. The resulting catalogs are publicly available. We include a basic summary of the strategy of the DEEP3 Galaxy Redshift Survey to accompany the recent public release of DEEP3 data. △ Less

Submitted 12 August, 2019; v1 submitted 19 March, 2019; originally announced March 2019.

Comments: 20 pages, 9 figures; published in MNRAS; for associated catalogs, see http://d-scholarship.pitt.edu/36064/

Journal ref: MNRAS 488, 4565-4584 (2019)

Evolution of the Gas Mass Fraction of Progenitors to Today's Massive Galaxies: ALMA Observations in the CANDELS GOODS-S Field

Authors: Tommy Wiklind, Henry C. Ferguson, Yicheng Guo, David C. Koo, Dale Kocevski, Bahram Mobasher, Gabriel B. Brammer, Susan Kassin, Anton M. Koekemoer, Mauro Giavalisco, Casey Papovich, Swara Ravindranath, Sandra M. Faber, Jonathan Freundlich, Duilia F. de Mello

Abstract: We present an ALMA survey of dust continuum emission in a sample of 70 galaxies in the redshift range z=2-5 selected from the CANDELS GOODS-S field. Multi-Epoch Abundance Matching (MEAM) is used to define potential progenitors of a z = 0 galaxy of stellar mass 1.5 10^11 M_sun. Gas masses are derived from the 850um luminosity. Ancillary data from the CANDELS GOODS-S survey are used to derive the ga… ▽ More We present an ALMA survey of dust continuum emission in a sample of 70 galaxies in the redshift range z=2-5 selected from the CANDELS GOODS-S field. Multi-Epoch Abundance Matching (MEAM) is used to define potential progenitors of a z = 0 galaxy of stellar mass 1.5 10^11 M_sun. Gas masses are derived from the 850um luminosity. Ancillary data from the CANDELS GOODS-S survey are used to derive the gas mass fractions. The results at z<=3 are mostly in accord with expectations: The detection rates are 75% for the z=2 redshift bin, 50% for the z=3 bin and 0% for z>=4. The average gas mass fraction for the detected z=2 galaxies is f_gas = 0.55+/-0.12 and f_gas = 0.62+/-0.15 for the z=3 sample. This agrees with expectations for galaxies on the star-forming main sequence, and shows that gas fractions have decreased at a roughly constant rate from z=3 to z=0. Stacked images of the galaxies not detected with ALMA give upper limits to f_gas of <0.08 and <0.15, for the z=2 and z=3 redshift bins. None of our galaxies in the z=4 and z=5 sample are detected and the upper limit from stacked images, corrected for low metallicity, is f_gas<0.66. We do not think that lower gas-phase metallicities can entirely explain the lower dust luminosities. We briefly consider the possibility of accretion of very low-metallicity gas to explain the absence of detectable dust emission in our galaxies at z>4. △ Less

Submitted 16 March, 2019; originally announced March 2019.

Comments: Accepted for publication in the Astrophysical Journal. 33 pages; 11 figures

Can intrinsic alignments of elongated low-mass galaxies be used to map the cosmic web at high redshift?

Authors: Viraj Pandya, Joel Primack, Peter Behroozi, Avishai Dekel, Haowen Zhang, Elliot Eckholm, Sandra M. Faber, Henry C. Ferguson, Mauro Giavalisco, Yicheng Guo, Nimish Hathi, Dritan Kodra, Anton M. Koekemoer, David C. Koo, Jeffrey Newman, Arjen van der Wel

Abstract: Hubble Space Telescope observations show that low-mass ($M_*=10^9-10^{10}M_{\odot}$) galaxies at high redshift ($z=1.0-2.5$) tend to be elongated (prolate) rather than disky (oblate) or spheroidal. This is explained in zoom-in cosmological hydrodynamical simulations by the fact that these galaxies are forming in cosmic web filaments where accretion happens preferentially along the direction of elo… ▽ More Hubble Space Telescope observations show that low-mass ($M_*=10^9-10^{10}M_{\odot}$) galaxies at high redshift ($z=1.0-2.5$) tend to be elongated (prolate) rather than disky (oblate) or spheroidal. This is explained in zoom-in cosmological hydrodynamical simulations by the fact that these galaxies are forming in cosmic web filaments where accretion happens preferentially along the direction of elongation. We ask whether the elongated morphology of these galaxies allows them to be used as effective tracers of cosmic web filaments at high redshift via their intrinsic alignments. Using mock lightcones and spectroscopically-confirmed galaxy pairs from the CANDELS survey, we test two types of alignments: (1) between the galaxy major axis and the direction to nearby galaxies of any mass, and (2) between the major axes of nearby pairs of low-mass, likely prolate, galaxies. The mock lightcones predict strong signals in 3D real space, 3D redshift space, and 2D projected redshift space for both types of alignments (assuming prolate galaxy orientations are the same as those of their host prolate halos), but we do not detect significant alignment signals in CANDELS observations. However, we show that spectroscopic redshifts have been obtained for only a small fraction of highly elongated galaxies, and accounting for spectroscopic incompleteness and redshift errors significantly degrades the 2D mock signal. This may partly explain the alignment discrepancy and highlights one of several avenues for future work. △ Less

Submitted 12 July, 2019; v1 submitted 25 February, 2019; originally announced February 2019.

Comments: Re-submitted to MNRAS after minor revisions

The structural properties of classical bulges and discs from z~2

Authors: Paola Dimauro, Marc Huertas-Company, Emanuele Daddi, Pablo G. Pérez-González, Mariangela Bernardi, Fernando Caro, Andrea Cattaneo, Boris Häußler, Ulrike Kuchner, Francesco Shankar, Guillermo Barro, Fernando Buitrago, Sandra M. Faber, Dale D. Kocevski, Anton M. Koekemoer, David C. Koo, Simona Mei, Reynier Peletier, Joel Primack, Aldo Rodriguez-Puebla, Mara Salvato, Diego Tuccillo

Abstract: We study the rest-frame optical mass-size relation of bulges and discs from z~2 to z~0 for a complete sample of massive galaxies in the CANDELS fields using 2 component Sérsic models (Dimauro et a. 2018). Discs and star forming galaxies follow similar mass-size relations. The mass-size relation of bulges is less steep than the one of quiescent galaxies (best fit slope of ~ 0.7 for quiescent galaxi… ▽ More We study the rest-frame optical mass-size relation of bulges and discs from z~2 to z~0 for a complete sample of massive galaxies in the CANDELS fields using 2 component Sérsic models (Dimauro et a. 2018). Discs and star forming galaxies follow similar mass-size relations. The mass-size relation of bulges is less steep than the one of quiescent galaxies (best fit slope of ~ 0.7 for quiescent galaxies against ~ 0.4 for bulges). We find little dependence of the structural properties of massive bulges and discs with the global morphology of galaxies (disc vs. bulge dominated) and the star formation activity (star-forming vs. quiescent). This result suggests similar bulge formation mechanisms for most massive galaxies and also that the formation of the bulge component does not significantly affect the disc structure. Our findings pose a challenge to models envisioning multiple channels for massive bulge growth, such as disc instabilities and mergers. △ Less

Submitted 13 February, 2019; v1 submitted 11 February, 2019; originally announced February 2019.

Comments: 20 pages, 15 figures, Submitted to MNRAS

Non-parametric Star Formation History Reconstruction with Gaussian Processes I: Counting Major Episodes of Star Formation

Authors: Kartheik G. Iyer, Eric Gawiser, Sandra M. Faber, Henry C. Ferguson, Anton M. Koekemoer, Camilla Pacifici, Rachel Somerville

Abstract: The star formation histories (SFHs) of galaxies contain imprints of the physical processes responsible for regulating star formation during galaxy growth and quenching. We improve the Dense Basis SFH reconstruction method of Iyer & Gawiser (2017), introducing a nonparametric description of the SFH based on the lookback times at which a galaxy assembles certain quantiles of its stellar mass. The me… ▽ More The star formation histories (SFHs) of galaxies contain imprints of the physical processes responsible for regulating star formation during galaxy growth and quenching. We improve the Dense Basis SFH reconstruction method of Iyer & Gawiser (2017), introducing a nonparametric description of the SFH based on the lookback times at which a galaxy assembles certain quantiles of its stellar mass. The method uses Gaussian Processes to create smooth SFHs that are independent of any functional form, with a flexible number of parameters that is adjusted to extract the maximum possible amount of SFH information from the SEDs being fit. We apply the method to reconstruct the SFHs of 48,791 galaxies with $H<25$ at $0.5 < z < 3.0$ across the five CANDELS fields. Using these SFHs, we study the evolution of galaxies as they grow more massive over cosmic time. We quantify the fraction of galaxies that show multiple major episodes of star formation, finding that the median time between two peaks of star formation is $\sim 0.42_{-0.10}^{+0.15}t_{univ}$ Gyr, where $t_{univ}$ is the age of the universe at a given redshift and remains roughly constant with stellar mass. Correlating SFHs with morphology, we find that studying the median SFHs of galaxies at $0.5<z<1.0$ at the same mass ($10^{10}< M_* < 10^{10.5}M_\odot$) allows us to compare the timescales on which the SFHs decline for different morphological classifications, ranging from $0.60^{-0.54}_{+1.54}$ Gyr for galaxies with spiral arms to $2.50^{-1.50}_{+2.25}$ Gyr for spheroids. The Gaussian Process-based SFH description provides a general approach to reconstruct smooth, nonparametric SFH posteriors for galaxies with a flexible number of parameters that can be incorporated into Bayesian SED fitting codes to minimize the bias in estimating physical parameters due to SFH parametrization. △ Less

Submitted 9 January, 2019; originally announced January 2019.

Comments: 26 pages, 15 figures. Submitted to ApJ

Galaxy inclination and the IRX-beta relation: Effects on UV star-formation rate measurements at intermediate to high redshifts

Authors: Weichen Wang, Susan A. Kassin, Camilla Pacifici, Guillermo Barro, Alexander de la Vega, Raymond C. Simons, S. M. Faber, Brett Salmon, Henry C. Ferguson, Pablo G. Perez-Gonzalez, Gregory F. Snyder, Karl D. Gordon, Zhu Chen, Dritan Kodra

Abstract: At intermediate and high redshifts, measurements of galaxy star-formation rates are usually based on rest-frame ultraviolet (UV) data. A correction for dust attenuation, A_UV, is needed for these measurements. This correction is typically inferred from UV spectral slopes (beta) using an equation known as "Meurer's Relation." In this paper, we study this relation at a redshift of 1.5 using images a… ▽ More At intermediate and high redshifts, measurements of galaxy star-formation rates are usually based on rest-frame ultraviolet (UV) data. A correction for dust attenuation, A_UV, is needed for these measurements. This correction is typically inferred from UV spectral slopes (beta) using an equation known as "Meurer's Relation." In this paper, we study this relation at a redshift of 1.5 using images and photometric measurements in the rest-frame UV (HST) through mid-infrared (Spitzer). It is shown that massive star-forming galaxies (above 10^10 Msun) have dust corrections that are dependent on their inclination to the line-of-sight. Edge-on galaxies have higher A_UV and infrared excess (IRX=L(IR)/L(UV)) than face-on galaxies at a given beta. Interestingly, dust corrections for low-mass star-forming galaxies do not depend on inclination. This is likely because more massive galaxies have more disk-like shapes/kinematics, while low-mass galaxies are more prolate and have more disturbed kinematics. To account for an inclination-dependent dust correction, a modified Meurer's Relation is derived: A_UV=4.43+1.99 beta - 1.73 (b/a-0.67), where b/a is the galaxy axis ratio. This inclination-dependence of A_UV can be explained by a two-component model of the dust distribution inside galaxies. In such a model, the dust attenuation of edge-on galaxies has a higher contribution from a "mixture" component (dust uniformly mixed with stars in the diffuse interstellar medium), and a lower contribution from a "birth cloud" component (near-spherical dust shells surrounding young stars in H II regions) than that of face-on galaxies. The difference is caused by the larger path-lengths through disks at higher inclinations. △ Less

Submitted 8 November, 2018; originally announced November 2018.

Comments: 25 pages, 14 figures. Accepted to ApJ. Figure 4 is the key plot

Evidence of Environmental Quenching at Redshift z ~ 2

Authors: Zhiyuan Ji, Mauro Giavalisco, Christina C. Williams, Sandra M. Faber, Henry C. Ferguson, Yicheng Guo, Teng Liu, Bomee Lee

Abstract: We report evidence of environmental quenching among galaxies at redshift ~ 2, namely the probability that a galaxy quenches its star formation activity is enhanced in the regions of space in proximity of other quenched, more massive galaxies. The effect is observed as strong clustering of quiescent galaxies around quiescent galaxies on angular scales θ< 20 arcsec, corresponding to a proper(comovin… ▽ More We report evidence of environmental quenching among galaxies at redshift ~ 2, namely the probability that a galaxy quenches its star formation activity is enhanced in the regions of space in proximity of other quenched, more massive galaxies. The effect is observed as strong clustering of quiescent galaxies around quiescent galaxies on angular scales θ< 20 arcsec, corresponding to a proper(comoving) scale of 168 (502) kpc at z = 2. The effect is observed only for quiescent galaxies around other quiescent galaxies; the probability to find star-forming galaxies around quiescent or around star-forming ones is consistent with the clustering strength of galaxies of the same mass and at the same redshift, as observed in dedicated studies of galaxy clustering. The effect is mass dependent in the sense that the quenching probability is stronger for galaxies of smaller mass ($\rm{M_*<10^{10} Msun}$) than for more massive ones, i.e. it follows the opposite trend with mass relative to gravitational galaxy clustering. The spatial scale where the effect is observed suggests these environments are massive halos, in which case the observed effect would likely be satellite quenching. The effect is also redshift dependent in that the clustering strength of quiescent galaxies around other quiescent galaxies at z = 1.6 is ~ 1.7 times larger than that of the galaxies with the same stellar mass at z = 2.6. This redshift dependence allows for a crude estimate of the time scale of environmental quenching of low-mass galaxies, which is in the range 1.5 - 4 Gyr, in broad agreement with other estimates and with our ideas on satellite quenching. △ Less

Submitted 11 June, 2018; originally announced June 2018.

Comments: 12 pages, 9 figures, Accepted for publication in ApJ

Journal ref: The Astrophysical Journal 862 (2018) 2-135

The evolution of galaxy shapes in CANDELS: from prolate to oblate

Authors: Haowen Zhang, Joel R. Primack, S. M. Faber, David C. Koo, Avishai Dekel, Zhu Chen, Daniel Ceverino, Yu-Yen Chang, Jerome J. Fang, Yicheng Guo, Lin Lin, Arjen van der Wel

Abstract: We model the projected b/a-log a distributions of CANDELS main sequence star-forming galaxies, where a (b) is the semi-major (semi-minor) axis of the galaxy images. We find that smaller-a galaxies are rounder at all stellar masses M and redshifts, so we include a when analyzing b/a distributions. Approximating intrinsic shapes of the galaxies as triaxial ellipsoids and assuming a multivariate norm… ▽ More We model the projected b/a-log a distributions of CANDELS main sequence star-forming galaxies, where a (b) is the semi-major (semi-minor) axis of the galaxy images. We find that smaller-a galaxies are rounder at all stellar masses M and redshifts, so we include a when analyzing b/a distributions. Approximating intrinsic shapes of the galaxies as triaxial ellipsoids and assuming a multivariate normal distribution of galaxy size and two shape parameters, we construct their intrinsic shape and size distributions to obtain the fractions of prolate, oblate and spheroidal galaxies in each redshift and mass bin. We find that galaxies tend to be prolate at low m and high redshifts, and oblate at high M and low redshifts, qualitatively consistent with van der Wel et al. (2014), implying that galaxies tend to evolve from prolate to oblate. These results are consistent with the predictions from simulations (Ceverino et al. 2015, Tomassetti et al. 2016) that the transition from prolate to oblate is caused by a compaction event at a characteristic mass range, making the galaxy center baryon dominated. We give probabilities of a galaxy's being prolate, oblate or spheroidal as a function of its M, redshift, projected b/a and a, which can facilitate target selections of galaxies with specific shapes at hight redshifts. We also give predicted optical depths of galaxies, which are qualitatively consistent with the expected correlation that AV should be higher for edge-on disk galaxies in each log a slice at low redshift and high mass bins. △ Less

Submitted 31 May, 2018; originally announced May 2018.

Comments: 24 pages, 25 figures, submitted to MNRAS

On the Transition of the Galaxy Quenching Mode at 0.5<z<1 in CANDELS

Authors: F. S. Liu, Meng Jia, Hassen M. Yesuf, S. M. Faber, David C. Koo, Yicheng Guo, Eric F. Bell, Dongfei Jiang, Weichen Wang, Anton M. Koekemoer, Xianzhong Zheng, Jerome J. Fang, Guillermo Barro, Pablo G. Pérez-González, Avishai Dekel, Dale Kocevski, Nimish P. Hathi, Darren Croton, M. Huertas-Company, Xianmin Meng, Wei Tong, Lu Liu

Abstract: We investigate the galaxy quenching process at intermediate redshift using a sample of $\sim4400$ galaxies with $M_{\ast} > 10^{9}M_{\odot}$ between redshift 0.5 and 1.0 in all five CANDELS fields. We divide this sample, using the integrated specific star formation rate (sSFR), into four sub-groups: star-forming galaxies (SFGs) above and below the ridge of the star-forming main sequence (SFMS), tr… ▽ More We investigate the galaxy quenching process at intermediate redshift using a sample of $\sim4400$ galaxies with $M_{\ast} > 10^{9}M_{\odot}$ between redshift 0.5 and 1.0 in all five CANDELS fields. We divide this sample, using the integrated specific star formation rate (sSFR), into four sub-groups: star-forming galaxies (SFGs) above and below the ridge of the star-forming main sequence (SFMS), transition galaxies and quiescent galaxies. We study their $UVI$ ($U-V$ versus $V-I$) color gradients to infer their sSFR gradients out to twice effective radii. We show that on average both star-forming and transition galaxies at all masses are not fully quenched at any radii, whereas quiescent galaxies are fully quenched at all radii. We find that at low masses ($M_{\ast} = 10^{9}-10^{10}M_{\odot}$) SFGs both above and below the SFMS ridge generally have flat sSFR profiles, whereas the transition galaxies at the same masses generally have sSFRs that are more suppressed in their outskirts. In contrast, at high masses ($M_{\ast} > 10^{10.5}M_{\odot}$), SFGs above and below the SFMS ridge and transition galaxies generally have varying degrees of more centrally-suppressed sSFRs relative to their outskirts. These findings indicate that at $z\sim~0.5-1.0$ the main galaxy quenching mode depends on its already formed stellar mass, exhibiting a transition from "the outside-in" at $M_{\ast} \leq 10^{10}M_{\odot}$ to "the inside-out" at $M_{\ast} > 10^{10.5}M_{\odot}$. In other words, our findings support that internal processes dominate the quenching of massive galaxies, whereas external processes dominate the quenching of low-mass galaxies. △ Less

Submitted 1 May, 2018; originally announced May 2018.

Comments: 19 pages, 12 figures, 2 tables, accepted for publication in ApJ

Is the dark-matter halo spin a predictor of galaxy spin and size?

Authors: Fangzhou Jiang, Avishai Dekel, Omer Kneller, Sharon Lapiner, Daniel Ceverino, Joel R. Primack, Sandra M. Faber, Andrea V. Macciò, Aaron Dutton, Shy Genel, Rachel S. Somerville

Abstract: The similarity between the distributions of spins for galaxies ($λ_{\rm g}$) and for dark-matter haloes ($λ_{\rm h}$), indicated both by simulations and observations, is naively interpreted as a one-to-one correlation between the spins of a galaxy and its host halo. This is used to predict galaxy sizes in semi-analytic models via $R_{\rm e}\simeqλ_{\rm h} R_{\rm v}$, with $R_{\rm e}$ the half-mass… ▽ More The similarity between the distributions of spins for galaxies ($λ_{\rm g}$) and for dark-matter haloes ($λ_{\rm h}$), indicated both by simulations and observations, is naively interpreted as a one-to-one correlation between the spins of a galaxy and its host halo. This is used to predict galaxy sizes in semi-analytic models via $R_{\rm e}\simeqλ_{\rm h} R_{\rm v}$, with $R_{\rm e}$ the half-mass radius of the galaxy and $R_{\rm v}$ the halo radius. Utilizing two different suites of zoom-in cosmological simulations, we find that $λ_{\rm g}$ and $λ_{\rm h}$ are in fact only barely correlated, especially at $z\geq 1$. A general smearing of this correlation is expected based on the different spin histories, where the more recently accreted baryons through streams gain and then lose significant angular momentum compared to the gradually accumulated dark matter. Expecting the spins of baryons and dark matter to be correlated at accretion into $R_{\rm v}$, the null correlation at the end reflects an anti-correlation between $λ_{\rm g}/λ_{\rm h}$ and $λ_{\rm h}$, which can partly arise from mergers and a compact star-forming phase that many galaxies undergo. On the other hand, the halo and galaxy spin vectors tend to be aligned, with a median $\cosθ=0.6$-0.7 between galaxy and halo, consistent with instreaming within a preferred plane. The galaxy spin is better correlated with the spin of the inner halo, but this largely reflects the effect of the baryons on the halo. Following the null spin correlation, $λ_{\rm h}$ is not a useful proxy for $R_{\rm e}$. While our simulations reproduce a general relation of the sort $R_{\rm e}=AR_{\rm vir}$, in agreement with observational estimates, the relation becomes tighter with $A=0.02(c/10)^{-0.7}$, where $c$ is the halo concentration, which in turn introduces a dependence on mass and redshift. △ Less

Submitted 19 April, 2018; originally announced April 2018.

Comments: 20 pages, 16 figures, submitted to MNRAS

A catalog of polychromatic bulge-disk decompositions of ~ 17.600 galaxies in CANDELS

Authors: Paola Dimauro, Marc Huertas-Company, Emanuele Daddi, Pablo G. Pérez-González, Mariangela Bernardi, Guillermo Barro, Fernando Buitrago, Fernando Caro, Andrea Cattaneo, Helena Dominguez-Sánchez, Sandra M. Faber, Boris Häußler, Dale D. Kocevski, Anton M. Koekemoer, David C. Koo, Christoph T. Lee, Simona Mei, Berta Margalef-Bentabol, Joel Primack, Aldo Rodriguez-Puebla, Mara Salvato, Francesco Shankar, Diego Tuccillo

Abstract: Understanding how bulges grow in galaxies is critical step towards unveiling the link between galaxy morphology and star-formation. To do so, it is necessary to decompose large sample of galaxies at different epochs into their main components (bulges and disks). This is particularly challenging, especially at high redshifts, where galaxies are poorly resolved. This work presents a catalog of bulge… ▽ More Understanding how bulges grow in galaxies is critical step towards unveiling the link between galaxy morphology and star-formation. To do so, it is necessary to decompose large sample of galaxies at different epochs into their main components (bulges and disks). This is particularly challenging, especially at high redshifts, where galaxies are poorly resolved. This work presents a catalog of bulge-disk decompositions of the surface brightness profiles of ~17.600 H-band selected galaxies in the CANDELS fields (F160W<23, 0<z<2) in 4 to 7 filters covering a spectral range of 430-1600nm. This is the largest available catalog of this kind up to z = 2. By using a novel approach based on deep-learning to select the best model to fit, we manage to control systematics arising from wrong model selection and obtain less contaminated samples than previous works. We show that the derived structural properties are within $\sim10-20\%$ of random uncertainties. We then fit stellar population models to the decomposed SEDs (Spectral Energy Distribution) of bulges and disks and derive stellar masses (and stellar mass bulge-to-total ratios) as well as rest-frame colors (U,V,J) for bulges and disks separately. All data products are publicly released with this paper and through the web page https://lerma.obspm.fr/huertas/form_CANDELS and will be used for scientific analysis in forthcoming works. △ Less

Submitted 30 March, 2018; v1 submitted 27 March, 2018; originally announced March 2018.

Comments: 18 pages, 18 figures, 3 table, submitted to MNRAS after the first referee's report

The VANDELS ESO public spectroscopic survey: observations and first data release

Authors: L. Pentericci, R. J. McLure B. Garilli, O. Cucciati, P. Franzetti, A. Iovino, R. Amorin, M. Bolzonella, A. Bongiorno, A. C. Carnall, M. Castellano, A. Cimatti, M. Cirasuolo, F. Cullen, S. DeBarros, J. S. Dunlop, D. Elbaz, S. Finkelstein, A. Fontana, F. Fontanot, M. Fumana, A. Gargiulo, L. Guaita, W. Hartley, M. Jarvis, S. Juneau , et al. (71 additional authors not shown)

Abstract: This paper describes the observations and the first data release (DR1) of the ESO public spectroscopic survey "VANDELS, a deep VIMOS survey of the CANDELS CDFS and UDS fields". VANDELS' main targets are star-forming galaxies at 2.4<z<5.5 and massive passive galaxies at 1<z<2.5. By adopting a strategy of ultra-long exposure times, from 20 to 80 hours per source, VANDELS is designed to be the deepes… ▽ More This paper describes the observations and the first data release (DR1) of the ESO public spectroscopic survey "VANDELS, a deep VIMOS survey of the CANDELS CDFS and UDS fields". VANDELS' main targets are star-forming galaxies at 2.4<z<5.5 and massive passive galaxies at 1<z<2.5. By adopting a strategy of ultra-long exposure times, from 20 to 80 hours per source, VANDELS is designed to be the deepest ever spectroscopic survey of the high-redshift Universe. Exploiting the red sensitivity of the VIMOS spectrograph, the survey has obtained ultra-deep spectra covering the wavelength 4800-10000 A with sufficient signal-to-noise to investigate the astrophysics of high-redshift galaxy evolution via detailed absorption line studies. The VANDELS-DR1 is the release of all spectra obtained during the first season of observations and includes data for galaxies for which the total (or half of the total) scheduled integration time was completed. The release contains 879 individual objects with a measured redshift and includes fully wavelength and flux-calibrated 1D spectra, the associated error spectra, sky spectra and wavelength-calibrated 2D spectra. We also provide a catalog with the essential galaxy parameters, including spectroscopic redshifts and redshift quality flags. In this paper we present the survey layout and observations, the data reduction and redshift measurement procedure and the general properties of the VANDELS-DR1 sample. We also discuss the spectroscopic redshift distribution, the accuracy of the photometric redshifts and we provide some examples of data products. All VANDELS-DR1 data are publicly available and can be retrieved from the ESO archive. Two further data releases are foreseen in the next 2 years with a final release scheduled for June 2020 which will include improved re-reduction of the entire spectroscopic data set. (abridged) △ Less

Submitted 20 March, 2018; originally announced March 2018.

Comments: Submitted to A&A

The Isophotal Structure of Star-forming Galaxies at $0.5< z <1.8$ in CANDELS: Implications for the Evolution of Galaxy Structure

Authors: Dongfei Jiang, F. S. Liu, Xianzhong Zheng, Hassen M. Yesuf, David C. Koo, S. M. Faber, Yicheng Guo, Anton M. Koekemoer, Weichen Wang, Jerome J. Fang, Guillermo Barro, Meng Jia, Wei Tong, Lu Liu, Xianmin Meng, Dale Kocevski, Elizabeth J. McGrath, Nimish P. Hathi

Abstract: We have measured the radial profiles of isophotal ellipticity ($\varepsilon$) and disky/boxy parameter A$_4$ out to radii of about three times the semi-major axes for $\sim4,600$ star-forming galaxies (SFGs) at intermediate redshifts $0.5<z<1.8$ in the CANDELS/GOODS-S and UDS fields. Based on the average size versus stellar-mass relation in each redshift bin, we divide our galaxies into Small SFGs… ▽ More We have measured the radial profiles of isophotal ellipticity ($\varepsilon$) and disky/boxy parameter A$_4$ out to radii of about three times the semi-major axes for $\sim4,600$ star-forming galaxies (SFGs) at intermediate redshifts $0.5<z<1.8$ in the CANDELS/GOODS-S and UDS fields. Based on the average size versus stellar-mass relation in each redshift bin, we divide our galaxies into Small SFGs (SSFGs), i.e., smaller than average for its mass, and Large SFGs (LSFGs), i.e., larger than average. We find that, at low masses ($M_{\ast} < 10^{10}M_{\odot}$), the SSFGs generally have nearly flat $\varepsilon$ and A$_4$ profiles for both edge-on and face-on views, especially at redshifts $z>1$. Moreover, the median A$_4$ values at all radii are almost zero. In contrast, the highly-inclined, low-mass LSFGs in the same mass-redshift bins generally have monotonically increasing $\varepsilon$ with radius and are dominated by disky values at intermediate radii. These findings at intermediate redshifts imply that low-mass SSFGs are not disk-like, while low-mass LSFGs appear to harbour disk-like components flattened by significant rotation. At high masses ($M_{\ast} > 10^{10}M_{\odot}$), highly-inclined SSFGs and LSFGs both exhibit a general, distinct trend for both $\varepsilon$ and A$_4$ profiles: increasing values with radius at lower radii, reaching maxima at intermediate radii, and then decreasing values at larger radii. Such a trend is more prevalent for more massive ($M_{\ast} > 10^{10.5}M_{\odot}$) galaxies or those at lower redshifts ($z<1.4$). The distinct trend in $\varepsilon$ and A$_4$ can be simply explained if galaxies possess all three components: central bulges, disks in the intermediate regions, and halo-like stellar components in the outskirts. △ Less

Submitted 22 February, 2018; originally announced February 2018.

Journal ref: The Astrophysical Journal, Volume 854, Issue 1, article id. 70, 16 pp. (2018)

Evidence for Merger-Driven Growth in Luminous, High-z, Obscured AGN in the CANDELS/COSMOS Field

Authors: J. L. Donley, J. Kartaltepe, D. Kocevski, M. Salvato, P. Santini, H. Suh, F. Civano, A. M. Koekemoer, J. Trump, M. Brusa, C. Cardamone, A. Castro, M. Cisternas, C. Conselice, D. Croton, N. Hathi, C. Liu, R. A. Lucas, P. Nair, D. Rosario, D. Sanders, B. Simmons, C. Villforth, D. M. Alexander, E. F. Bell , et al. (5 additional authors not shown)

Abstract: While major mergers have long been proposed as a driver of both AGN activity and the M-sigma relation, studies of moderate to high redshift Seyfert-luminosity AGN hosts have found little evidence for enhanced rates of interactions. However, both theory and observation suggest that while these AGN may be fueled by stochastic accretion and secular processes, high-luminosity, high-redshift, and heavi… ▽ More While major mergers have long been proposed as a driver of both AGN activity and the M-sigma relation, studies of moderate to high redshift Seyfert-luminosity AGN hosts have found little evidence for enhanced rates of interactions. However, both theory and observation suggest that while these AGN may be fueled by stochastic accretion and secular processes, high-luminosity, high-redshift, and heavily obscured AGN are the AGN most likely to be merger-driven. To better sample this population of AGN, we turn to infrared selection in the CANDELS/COSMOS field. Compared to their lower-luminosity and less obscured X-ray-only counterparts, IR-only AGN (luminous, heavily obscured AGN) are more likely to be classified as either irregular (50$^{+12}_{-12}$% vs. 9$^{+5}_{-2}$%) or asymmetric (69$^{+9}_{-13}$% vs. 17$^{+6}_{-4}$%) and are less likely to have a spheroidal component (31$^{+13}_{-9}$% vs. 77$^{+4}_{-6}$%). Furthermore, IR-only AGN are also significantly more likely than X-ray-only AGN (75$^{+8}_{-13}$% vs. 31$^{+6}_{-6}$%) to be classified either as interacting or merging in a way that significantly disturbs the host galaxy or disturbed though not clearly interacting or merging, which potentially represents the late stages of a major merger. This suggests that while major mergers may not contribute significantly to the fueling of Seyfert luminosity AGN, interactions appear to play a more dominant role in the triggering and fueling of high-luminosity heavily obscured AGN. △ Less

Submitted 6 December, 2017; originally announced December 2017.

Comments: 13 pages, 7 figures, 1 table, Accepted for publication in the Astrophysical Journal

Clumpy Galaxies in CANDELS. II. Physical Properties of UV-bright Clumps at $0.5\leq z<3$

Authors: Yicheng Guo, Marc Rafelski, Eric F. Bell, Christopher J. Conselice, Avishai Dekel, S. M. Faber, Mauro Giavalisco, Anton M. Koekemoer, David C. Koo, Yu Lu, Nir Mandelker, Joel R. Primack, Daniel Ceverino, Duilia F. de Mello, Henry C. Ferguson, Nimish Hathi, Dale Kocevski, Ray A. Lucas, Pablo G. Pérez-González, Swara Ravindranath, Emmaris Soto, Amber Straughn, Weichen Wang

Abstract: Studying giant star-forming clumps in distant galaxies is important to understand galaxy formation and evolution. At present, however, observers and theorists have not reached a consensus on whether the observed "clumps" in distant galaxies are the same phenomenon that is seen in simulations. In this paper, as a step to establish a benchmark of direct comparisons between observations and theories,… ▽ More Studying giant star-forming clumps in distant galaxies is important to understand galaxy formation and evolution. At present, however, observers and theorists have not reached a consensus on whether the observed "clumps" in distant galaxies are the same phenomenon that is seen in simulations. In this paper, as a step to establish a benchmark of direct comparisons between observations and theories, we publish a sample of clumps constructed to represent the commonly observed "clumps" in the literature. This sample contains 3193 clumps detected from 1270 galaxies at $0.5 \leq z < 3.0$. The clumps are detected from rest-frame UV images, as described in our previous paper. Their physical properties, e.g., rest-frame color, stellar mass (M*), star formation rate (SFR), age, and dust extinction, are measured by fitting the spectral energy distribution (SED) to synthetic stellar population models. We carefully test the procedures of measuring clump properties, especially the method of subtracting background fluxes from the diffuse component of galaxies. With our fiducial background subtraction, we find a radial clump U-V color variation, where clumps close to galactic centers are redder than those in outskirts. The slope of the color gradient (clump color as a function of their galactocentric distance scaled by the semi-major axis of galaxies) changes with redshift and M* of the host galaxies: at a fixed M*, the slope becomes steeper toward low redshift, and at a fixed redshift, it becomes slightly steeper with M*. Based on our SED-fitting, this observed color gradient can be explained by a combination of a negative age gradient, a negative E(B-V) gradient, and a positive specific star formation rate gradient of the clumps. We also find that the color gradients of clumps are steeper than those of intra-clump regions. [Abridged] △ Less

Submitted 5 December, 2017; originally announced December 2017.

Comments: 23 pages, 14 figures. Accepted by ApJ. The clump catalog will be published online in ApJ. It can also be requested through emails to the first author or this link http://faculty.missouri.edu/guoyic/data/clump/

Spatially resolved kinematics in the central 1 kpc of a compact star-forming galaxy at z=2.3 from ALMA CO observations

Authors: G. Barro, M. Kriek, P. G. Perez-Gonzalez, T. Diaz-Santos, S. H. Price, W. Rujopakarn, V. Pandya, D. C. Koo, S. M. Faber, A. Dekel, J. R. Primack, D. D. Kocevski

Abstract: We present high spatial resolution (FWHM$\sim$0.14'') observations of the CO($8-7$) line in GDS-14876, a compact star-forming galaxy at $z=2.3$ with total stellar mass of $\log(M_{\star}/M_{\odot})=10.9$. The spatially resolved velocity map of the inner $r\lesssim1$~kpc reveals a continous velocity gradient consistent with the kinematics of a rotating disk with $v_{\rm rot}(r=1\rm kpc)=163\pm5$ km… ▽ More We present high spatial resolution (FWHM$\sim$0.14'') observations of the CO($8-7$) line in GDS-14876, a compact star-forming galaxy at $z=2.3$ with total stellar mass of $\log(M_{\star}/M_{\odot})=10.9$. The spatially resolved velocity map of the inner $r\lesssim1$~kpc reveals a continous velocity gradient consistent with the kinematics of a rotating disk with $v_{\rm rot}(r=1\rm kpc)=163\pm5$ km s$^{-1}$ and $v_{\rm rot}/σ\sim2.5$. The gas-to-stellar ratios estimated from CO($8-7$) and the dust continuum emission span a broad range, $f^{\rm CO}_{\rm gas}=M_{\rm gas}/M_{\star}=13-45\%$ and $f^{\rm cont}_{\rm gas}=50-67\%$, but are nonetheless consistent given the uncertainties in the conversion factors. The dynamical modeling yields a dynamical mass of$\log(M_{\rm dyn}/M_{\odot})=10.58^{+0.5}_{-0.2}$ which is lower, but still consistent with the baryonic mass, $\log$(M$_{\rm bar}$= M$_{\star}$ + M$^{\rm CO}_{\rm gas}$/M$_{\odot}$)$=11.0$, if the smallest CO-based gas fraction is assumed. Despite a low, overall gas fraction, the small physical extent of the dense, star-forming gas probed by CO($8-7$), $\sim3\times$ smaller than the stellar size, implies a strong concentration that increases the gas fraction up to $f^{\rm CO, 1\rm kpc}_{\rm gas}\sim 85\%$ in the central 1 kpc. Such a gas-rich center, coupled with a high star-formation rate, SFR$\sim$ 500 M$_{\odot}$ yr$^{-1}$, suggests that GDS-14876 is quickly assembling a dense stellar component (bulge) in a strong nuclear starburst. Assuming its gas reservoir is depleted without replenishment, GDS-14876 will quickly ($t_{\rm depl}\sim27$ Myr) become a compact quiescent galaxy that could retain some fraction of the observed rotational support. △ Less

Submitted 4 December, 2017; originally announced December 2017.

Comments: Accepted for Publication in ApJL. Kinematic maps are shown in Figures 2 and 4

CANDELS: Elevated Black Hole Growth in the Progenitors of Compact Quiescent Galaxies at z~2

Authors: Dale D. Kocevski, Guillermo Barro, S. M. Faber, Avishai Dekel, Rachel S. Somerville, Joshua A. Young, Christina C. Williams, Daniel H. McIntosh, Antonis Georgakakis, Guenther Hasinger, Kirpal Nandra, Francesca Civano, David M. Alexander, Omar Almaini, Christopher J. Conselice, Jennifer L. Donley, Harry C. Ferguson, Mauro Giavalisco, Norman A. Grogin, Nimish Hathi, Matthew Hawkins, Anton M. Koekemoer, David C. Koo, Elizabeth J. McGrath, Bahram Mobasher , et al. (9 additional authors not shown)

Abstract: We examine the fraction of massive ($M_{*}>10^{10} M_{\odot}$), compact star-forming galaxies (cSFGs) that host an active galactic nucleus (AGN) at $z\sim2$. These cSFGs are likely the direct progenitors of the compact quiescent galaxies observed at this epoch, which are the first population of passive galaxies to appear in large numbers in the early Universe. We identify cSFGs that host an AGN us… ▽ More We examine the fraction of massive ($M_{*}>10^{10} M_{\odot}$), compact star-forming galaxies (cSFGs) that host an active galactic nucleus (AGN) at $z\sim2$. These cSFGs are likely the direct progenitors of the compact quiescent galaxies observed at this epoch, which are the first population of passive galaxies to appear in large numbers in the early Universe. We identify cSFGs that host an AGN using a combination of Hubble WFC3 imaging and Chandra X-ray observations in four fields: the Chandra Deep Fields, the Extended Groth Strip, and the UKIDSS Ultra Deep Survey field. We find that $39.2^{+3.9}_{-3.6}$\% (65/166) of cSFGs at $1.4<z<3.0$ host an X-ray detected AGN. This fraction is 3.2 times higher than the incidence of AGN in extended star-forming galaxies with similar masses at these redshifts. This difference is significant at the $6.2σ$ level. Our results are consistent with models in which cSFGs are formed through a dissipative contraction that triggers a compact starburst and concurrent growth of the central black hole. We also discuss our findings in the context of cosmological galaxy evolution simulations that require feedback energy to rapidly quench cSFGs. We show that the AGN fraction peaks precisely where energy injection is needed to reproduce the decline in the number density of cSFGs with redshift. Our results suggest that the first abundant population of massive, quenched galaxies emerged directly following a phase of elevated supermassive black hole growth and further hints at a possible connection between AGN and the rapid quenching of star formation in these galaxies. △ Less

Submitted 16 October, 2017; originally announced October 2017.

Comments: Published in the Astrophysical Journal

Journal ref: 2017ApJ...846..112K

Demographics of Star-forming Galaxies since $z\sim2.5$. I. The $UVJ$ Diagram in CANDELS

Authors: Jerome J. Fang, S. M. Faber, David C. Koo, Aldo Rodriguez-Puebla, Yicheng Guo, Guillermo Barro, Peter Behroozi, Gabriel Brammer, Zhu Chen, Avishai Dekel, Henry C. Ferguson, Eric Gawiser, Mauro Giavalisco, Jeyhan Kartaltepe, Dale D. Kocevski, Anton M. Koekemoer, Elizabeth J. McGrath, Daniel McIntosh, Jeffrey A. Newman, Camilla Pacifici, Viraj Pandya, Pablo G. Perez-Gonzalez, Joel R. Primack, Brett Salmon, Jonathan R. Trump , et al. (19 additional authors not shown)

Abstract: This is the first in a series of papers examining the demographics of star-forming galaxies at $0.2<z<2.5$ in CANDELS. We study 9,100 galaxies from GOODS-S and UDS having published values of redshifts, masses, star-formation rates (SFRs), and dust attenuation ($A_V$) derived from UV-optical SED fitting. In agreement with previous works, we find that the $UVJ$ colors of a galaxy are closely correla… ▽ More This is the first in a series of papers examining the demographics of star-forming galaxies at $0.2<z<2.5$ in CANDELS. We study 9,100 galaxies from GOODS-S and UDS having published values of redshifts, masses, star-formation rates (SFRs), and dust attenuation ($A_V$) derived from UV-optical SED fitting. In agreement with previous works, we find that the $UVJ$ colors of a galaxy are closely correlated with its specific star-formation rate (SSFR) and $A_V$. We define rotated $UVJ$ coordinate axes, termed $S_\mathrm{SED}$ and $C_\mathrm{SED}$, that are parallel and perpendicular to the star-forming sequence and derive a quantitative calibration that predicts SSFR from $C_\mathrm{SED}$ with an accuracy of ~0.2 dex. SFRs from UV-optical fitting and from UV+IR values based on Spitzer/MIPS 24 $μ\mathrm{m}$ agree well overall, but systematic differences of order 0.2 dex exist at high and low redshifts. A novel plotting scheme conveys the evolution of multiple galaxy properties simultaneously, and dust growth, as well as star-formation decline and quenching, exhibit "mass-accelerated evolution" ("downsizing"). A population of transition galaxies below the star-forming main sequence is identified. These objects are located between star-forming and quiescent galaxies in $UVJ$ space and have lower $A_V$ and smaller radii than galaxies on the main sequence. Their properties are consistent with their being in transit between the two regions. The relative numbers of quenched, transition, and star-forming galaxies are given as a function of mass and redshift. △ Less

Submitted 6 April, 2018; v1 submitted 15 October, 2017; originally announced October 2017.

Comments: 36 pages, 26 figures, ApJ accepted

Type Ia Supernova Distances at z > 1.5 from the Hubble Space Telescope Multi-Cycle Treasury Programs: The Early Expansion Rate

Authors: Adam G. Riess, Steven A. Rodney, Daniel M. Scolnic, Daniel L. Shafer, Louis-Gregory Strolger, Henry C. Ferguson, Marc Postman, Or Graur, Dan Maoz, Saurabh W. Jha, Bahram Mobasher, Stefano Casertano, Brian Hayden, Alberto Molino, Jens Hjorth, Peter M. Garnavich, David O. Jones, Robert P. Kirshner, Anton M. Koekemoer, Norman A. Grogin, Gabriel Brammer, Shoubaneh Hemmati, Mark Dickinson, Peter M. Challis, Schuyler Wolff , et al. (9 additional authors not shown)

Abstract: We present an analysis of 15 Type Ia supernovae (SNe Ia) at redshift z > 1 (9 at 1.5 < z < 2.3) recently discovered in the CANDELS and CLASH Multi-Cycle Treasury programs using WFC3 on the Hubble Space Telescope. We combine these SNe Ia with a new compilation of 1050 SNe Ia, jointly calibrated and corrected for simulated survey biases to produce accurate distance measurements. We present unbiased… ▽ More We present an analysis of 15 Type Ia supernovae (SNe Ia) at redshift z > 1 (9 at 1.5 < z < 2.3) recently discovered in the CANDELS and CLASH Multi-Cycle Treasury programs using WFC3 on the Hubble Space Telescope. We combine these SNe Ia with a new compilation of 1050 SNe Ia, jointly calibrated and corrected for simulated survey biases to produce accurate distance measurements. We present unbiased constraints on the expansion rate at six redshifts in the range 0.07 < z < 1.5 based only on this combined SN Ia sample. The added leverage of our new sample at z > 1.5 leads to a factor of ~3 improvement in the determination of the expansion rate at z = 1.5, reducing its uncertainty to ~20%, a measurement of H(z=1.5)/H0=2.67 (+0.83,-0.52). We then demonstrate that these six measurements alone provide a nearly identical characterization of dark energy as the full SN sample, making them an efficient compression of the SN Ia data. The new sample of SNe Ia at z > 1 usefully distinguishes between alternative cosmological models and unmodeled evolution of the SN Ia distance indicators, placing empirical limits on the latter. Finally, employing a realistic simulation of a potential WFIRST SN survey observing strategy, we forecast optimistic future constraints on the expansion rate from SNe Ia. △ Less

Submitted 2 October, 2017; originally announced October 2017.

Comments: 14 pages, 5 figures, 7 tables; submitted to ApJ