Showing 1–32 of 32 results for author: Kawinwanichakij, L

A massive galaxy that formed its stars at $z \sim 11$

Authors: Karl Glazebrook, Themiya Nanayakkara, Corentin Schreiber, Claudia Lagos, Lalitwadee Kawinwanichakij, Colin Jacobs, Harry Chittenden, Gabriel Brammer, Glenn G. Kacprzak, Ivo Labbe, Danilo Marchesini, Z. Cemile Marsan, Pascal A. Oesch, Casey Papovich, Rhea-Silvia Remus, Kim-Vy H. Tran, James Esdaile, Angel Chandro Gomez

Abstract: The formation of galaxies by gradual hierarchical co-assembly of baryons and cold dark matter halos is a fundamental paradigm underpinning modern astrophysics and predicts a strong decline in the number of massive galaxies at early cosmic times. Extremely massive quiescent galaxies (stellar masses $>10^{11}$ M$_\odot$) have now been observed as early as 1-2 billions years after the Big Bang; these… ▽ More The formation of galaxies by gradual hierarchical co-assembly of baryons and cold dark matter halos is a fundamental paradigm underpinning modern astrophysics and predicts a strong decline in the number of massive galaxies at early cosmic times. Extremely massive quiescent galaxies (stellar masses $>10^{11}$ M$_\odot$) have now been observed as early as 1-2 billions years after the Big Bang; these are extremely constraining on theoretical models as they form 300-500 Myr earlier and only some models can form massive galaxies this early. Here we report on the spectroscopic observations with the James Webb Space Telescope of a massive quiescent galaxy ZF-UDS-7329 at redshift 3.205 $\pm$ 0.005 that eluded deep ground-based spectrscopy, is significantly redder than typical and whose spectrum reveals features typical of much older stellar populations. Detailed modelling shows the stellar population formed around 1.5 billion years earlier in time (z ~ 11) at an epoch when dark matter halos of sufficient hosting mass have not yet assembled in the standard scenario. This observation may point to the presence of undetected populations of early galaxies and the possibility of significant gaps in our understanding of early stellar populations, galaxy formation and/or the nature of dark matter. △ Less

Submitted 3 May, 2024; v1 submitted 10 August, 2023; originally announced August 2023.

Comments: 21 pages, 6 figures, updated to reflect accepted version. v3 with figure caption clarified

Journal ref: Nature 2024

The eROSITA Final Equatorial-Depth Survey (eFEDS): Host-galaxy Demographics of X-ray AGNs with Subaru Hyper Suprime-Cam

Authors: Junyao Li, John D. Silverman, Andrea Merloni, Mara Salvato, Johannes Buchner, Andy Goulding, Teng Liu, Riccardo Arcodia, Johan Comparat, Xuheng Ding, Kohei Ichikawa, Masatoshi Imanishi, Toshihiro Kawaguchi, Lalitwadee Kawinwanichakij, Yoshiki Toba

Abstract: We investigate the physical properties, such as star-forming activity, disk vs. bulge nature, galaxy size, and obscuration of 3796 X-ray selected AGNs at $0.2<z<0.8$ in the eFEDS field. Using Subaru Hyper Suprime-Cam imaging data in the $grizy$ bands for SRG/eROSITA-detected AGNs, we measure the structural parameters for AGN host galaxies by performing a 2D AGN-host image decomposition. We then co… ▽ More We investigate the physical properties, such as star-forming activity, disk vs. bulge nature, galaxy size, and obscuration of 3796 X-ray selected AGNs at $0.2<z<0.8$ in the eFEDS field. Using Subaru Hyper Suprime-Cam imaging data in the $grizy$ bands for SRG/eROSITA-detected AGNs, we measure the structural parameters for AGN host galaxies by performing a 2D AGN-host image decomposition. We then conduct spectral energy distribution fitting to derive stellar mass and rest-frame colors for AGN hosts. We find that (1) AGNs can contribute significantly to the total optical light down to ${\rm log}\,L_{\rm X}\sim 42.5\ \rm erg\,s^{-1}$, thus ignoring the AGN component can significantly bias the structural measurements; (2) AGN hosts are predominately star-forming galaxies at ${\rm log}\,\mathcal{M}_\star \lesssim 11.3\ M_\odot$; (3) the bulk of AGNs (64%) reside in galaxies with significant stellar disks, while their host galaxies become increasingly bulge dominated and quiescent at ${\rm log}\,\mathcal{M}_\star \gtrsim 11.0\ M_\odot$; (4) the size-stellar mass relation of AGN hosts tends to lie between that of inactive star-forming and quiescent galaxies, suggesting that the physical mechanism responsible for building the central stellar density also efficiently fuel the black hole growth; (5) the hosts of X-ray unobscured AGNs are biased towards face-on systems and the average $E(B-V)/N_{\rm H}$ is similar to the galactic dust-to-gas ratio, suggesting that some of the obscuration of the nuclei could come from galaxy-scale gas and dust, which may partly account for (up to 30%) the deficiency of star-forming disks as host galaxies for the most massive AGNs. These results are consistent with a scenario in which the black hole and galaxy grow in mass while transform in structure and star-forming activity, as desired to establish the local scaling relations. △ Less

Submitted 23 February, 2023; originally announced February 2023.

Comments: 15 pages, 13 figures, submitted to MNRAS

The Spitzer-HETDEX Exploratory Large Area Survey. IV. Model-Based Multi-wavelength Photometric Catalog

Authors: Gene C. K. Leung, Steven Finkelstein, John Weaver, Casey Papovich, Rebecca Larson, Katherine Chworowsky, Robin Ciardullo, Eric Gawiser, Caryl Gronwall, Shardha Jogee, Lalitwadee Kawinwanichakij, Rachel Somerville, Isak Wold, L. Y. Aaron Yung

Abstract: We present a 0.3--4.5 $μ$m 16-band photometric catalog for the Spitzer/HETDEX Exploratory Large-Area (SHELA) survey. SHELA covers a $\sim 27$ deg$^2$ field within the footprint of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). Here we present new DECam imaging and a $rizK_s$-band-selected catalog of four million sources extracted using a fully model-based approach. We validate our pho… ▽ More We present a 0.3--4.5 $μ$m 16-band photometric catalog for the Spitzer/HETDEX Exploratory Large-Area (SHELA) survey. SHELA covers a $\sim 27$ deg$^2$ field within the footprint of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). Here we present new DECam imaging and a $rizK_s$-band-selected catalog of four million sources extracted using a fully model-based approach. We validate our photometry by comparing with the model-based DECam Legacy Survey. We analyze the differences between model-based and aperture photometry by comparing with the previous SHELA catalog, finding that our model-based photometry can measure point sources to fainter fluxes and better capture the full emission of resolved sources. The catalog is $80\%$ ($50\%$) complete at $riz \sim 24.7$ ($25.1$) AB mag, and the optical photometry reaches a $5σ$ depth of $\sim 25.5$ AB mag. We measure photometric redshifts and achieve $1σ$ scatter of $Δz/(1+z)$ of 0.04 with available spectroscopic redshifts at $0 \le z \le 1$. This large area, multi-wavelength photometric catalog, combined with spectroscopic information from HETDEX, will enable a wide range of extragalactic science investigations. △ Less

Submitted 2 January, 2023; originally announced January 2023.

Comments: 22 pages, 12 figures, submitted to ApJ

A population of faint, old, and massive quiescent galaxies at 3 < z < 4 revealed by JWST NIRSpec Spectroscopy

Authors: Themiya Nanayakkara, Karl Glazebrook, Colin Jacobs, Lalitwadee Kawinwanichakij, Corentin Schreiber, Gabriel Brammer, James Esdaile, Glenn G. Kacprzak, Ivo Labbe, Claudia Lagos, Danilo Marchesini, Z. Cemile Marsan, Pascal A. Oesch, Casey Papovich, Rhea-Silvia Remus, Kim-Vy H. Tran

Abstract: Here we present a sample of 12 massive quiescent galaxy candidates at z~3-4 observed with the James Webb Space Telescope (JWST) Near Infrared Spectrograph (NIRSpec). These galaxies were pre-selected from the Hubble Space Telescope imaging and 10 of our sources were unable to be spectroscopically confirmed by ground based spectroscopy. By combining spectroscopic data from NIRSpec with multi-wavelen… ▽ More Here we present a sample of 12 massive quiescent galaxy candidates at z~3-4 observed with the James Webb Space Telescope (JWST) Near Infrared Spectrograph (NIRSpec). These galaxies were pre-selected from the Hubble Space Telescope imaging and 10 of our sources were unable to be spectroscopically confirmed by ground based spectroscopy. By combining spectroscopic data from NIRSpec with multi-wavelength imaging data from the JWST Near Infrared Camera (NIRCam), we analyse their stellar populations and their formation histories. We find that all of our galaxies classify as quiescent based on the reconstruction of their star formation histories but show a variety of quenching timescales and ages. All our galaxies are massive ($\sim0.1-1.2 \times 10^{11} M\odot$), with masses comparable to massive galaxies in the local Universe. We find that the oldest galaxy in our sample formed $\sim1.0\times10^{11} M\odot$ of mass within the first few hundred million years of the Universe and has been quenched for more than a billion years by the time of observation at z$\sim$3.2 ($\sim$2 billion years after the Big Bang). Our results point to very early formation of massive galaxies requiring a high conversion rate of baryons to stars in the early Universe. △ Less

Submitted 14 February, 2024; v1 submitted 22 December, 2022; originally announced December 2022.

Comments: Published in Nature Scientific Reports

Journal ref: Scientific Reports volume 14, Article number: 3724 (2024)

Dark matter halos of luminous AGNs from galaxy-galaxy lensing with the HSC Subaru Strategic Program

Authors: Wentao Luo, John D. Silverman, Surhud More, Andy Goulding, Hironao Miyatake, Takahiro Nishimichi, Chiaki Hikage, Lalitwadee Kawinwanichakij, Junyao Li, Xiangchong Li, Elinor Medezinski, Masamune Oguri, Taira Oogi, Cristobal Sifon

Abstract: We assess the dark matter halo masses of luminous AGNs over the redshift range 0.2 to 1.2 using galaxy-galaxy lensing based on imaging data from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). We measure the weak lensing signal of a sample of 8882 AGNs constructed using HSC and WISE photometry. The lensing detection around AGNs has a signal-to-noise ratio of 15. As expected, we find that… ▽ More We assess the dark matter halo masses of luminous AGNs over the redshift range 0.2 to 1.2 using galaxy-galaxy lensing based on imaging data from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). We measure the weak lensing signal of a sample of 8882 AGNs constructed using HSC and WISE photometry. The lensing detection around AGNs has a signal-to-noise ratio of 15. As expected, we find that the lensing mass profile is consistent with that of massive galaxies ($M_{*}\sim 10.8~M_\odot$). Surprisingly, the lensing signal remains unchanged when the AGN sample is split into low and high stellar mass hosts. Specifically, we find that the excess surface density (ESD) of AGNs, residing in galaxies with high stellar masses, significantly differs from that of the control sample. We further fit a halo occupation distribution model to the data to infer the posterior distribution of parameters including the average halo mass. We find that the characteristic halo mass of the full AGN population lies near the knee ($\rm log(M_h/h^{-1}M_{\odot})=12.0$) of the stellar-to-halo mass relation (SHMR). Illustrative of the results given above, the halo masses of AGNs residing in host galaxies with high stellar masses (i.e., above the knee of the SHMR) falls below the calibrated SHMR while the halo mass of the low stellar mass sample is more consistent with the established SHMR. These results indicate that massive halos with higher clustering bias tend to suppress AGN activity, probably due to the lack of available gas. △ Less

Submitted 11 April, 2022; v1 submitted 7 April, 2022; originally announced April 2022.

Comments: 17 pages,10 figures, accepted for publication in the Astrophysical Journal

AGN and Star Formation at Cosmic Noon: Comparison of Data to Theoretical Models

Authors: Jonathan Florez, Shardha Jogee, Yuchen Guo, Sofía A. Cora, Rainer Weinberger, Romeel Davé, Lars Hernquist, Mark Vogelsberger, Robin Ciardullo, Steven L. Finkelstein, Caryl Gronwall, Lalitwadee Kawinwanichakij, Gene C. K. Leung, Stephanie LaMassa, Casey Papovich, Matthew L. Stevans, Isak Wold

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

Submitted 22 September, 2021; originally announced September 2021.

Comments: 20 pages, 7 figures. Accepted for publication in MNRAS

Hyper Suprime-Cam Subaru Strategic Program: A Mass-Dependent Slope of the Galaxy Size-Mass Relation at $z<1$

Authors: Lalitwadee Kawinwanichakij, John D. Silverman, Xuheng Ding, Angelo George, Ivana Damjanov, Marcin Sawicki, Masayuki Tanaka, Dan S. Taranu, Simon Birrer, Song Huang, Junyao Li, Masato Onodera, Takatoshi Shibuya, Naoki Yasuda

Abstract: We present the galaxy size-mass ($R_{e}-M_{\ast}$) distributions using a stellar-mass complete sample of $\sim1.5$ million galaxies, covering $\sim100$ deg$^2$, with $\log(M_{\ast}/M_{\odot})>10.2~(9.2)$ over the redshift range $0.2<z<1.0$ $(z<0.6)$ from the second public data release of the Hyper Suprime-Cam Subaru Strategic Program. We confirm that, at fixed redshift and stellar mass over the ra… ▽ More We present the galaxy size-mass ($R_{e}-M_{\ast}$) distributions using a stellar-mass complete sample of $\sim1.5$ million galaxies, covering $\sim100$ deg$^2$, with $\log(M_{\ast}/M_{\odot})>10.2~(9.2)$ over the redshift range $0.2<z<1.0$ $(z<0.6)$ from the second public data release of the Hyper Suprime-Cam Subaru Strategic Program. We confirm that, at fixed redshift and stellar mass over the range of $\log(M_{\ast}/M_{\odot})<11$, star-forming galaxies are on average larger than quiescent galaxies. The large sample of galaxies with accurate size measurements, thanks to the excellent imaging quality, also enables us to demonstrate that the $R_{e}-M_{\ast}$ relations of both populations have a form of broken power-law, with a clear change of slopes at a pivot stellar mass $M_{p}$. For quiescent galaxies, below an (evolving) pivot mass of $\log(M_{p}/M_{\odot})=10.2-10.6$ the relation follows $R_{e}\propto M_{\ast}^{0.1}$; above $M_{p}$ the relation is steeper and follows $R_{e}\propto M_{\ast}^{0.6-0.7}$. For star-forming galaxies, below $\log(M_{p}/M_{\odot})\sim10.7$ the relation follows $R_{e}\propto M_{\ast}^{0.2}$; above $M_{p}$ the relation evolves with redshift and follows $R_{e}\propto M_{\ast}^{0.3-0.6}$. The shallow power-law slope for quiescent galaxies below $M_{p}$ indicates that large low-mass quiescent galaxies have sizes similar to those of their counterpart star-forming galaxies. We take this as evidence that large low-mass quiescent galaxies have been recently quenched (presumably through environment-specific process) without significant structural transformation. Interestingly, the pivot stellar mass of the $R_{e}-M_{\ast}$ relations coincides with mass at which half of the galaxy population is quiescent, implied that the pivot mass represents the transition of galaxy growth from being dominated by in-situ star formation to being dominated by (dry) mergers. △ Less

Submitted 20 September, 2021; originally announced September 2021.

Comments: Accepted for publication in ApJ. 45 pages, 17 figures

A Census of the Bright z=8.5-11 Universe with the Hubble and Spitzer Space Telescopes in the CANDELS Fields

Authors: Steven L. Finkelstein, Micaela Bagley, Mimi Song, Rebecca Larson, Casey Papovich, Mark Dickinson, Keely Finkelstein, Anton M. Koekemoer, Norbert Pirzkal, Rachel S. Somerville, L. Y. Aaron Yung, Peter Behroozi, Harry Ferguson, Mauro Giavalisco, Norman Grogin, Nimish Hathi, Taylor Hutchison, Intae Jung, Dale Kocevski, Lalitwadee Kawinwanichakij, Sofia Rojas-Ruiz, Russell Ryan Jr., Gregory F. Snyder, Sandro Tacchella

Abstract: We present the results from a new search for candidate galaxies at z ~ 8.5-11 discovered over the 850 arcmin^2 area probed by the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS). We use a photometric redshift selection including both Hubble and Spitzer Space Telescope photometry to robustly identify galaxies in this epoch at F160W < 26.6. We use a detailed vetting procedur… ▽ More We present the results from a new search for candidate galaxies at z ~ 8.5-11 discovered over the 850 arcmin^2 area probed by the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS). We use a photometric redshift selection including both Hubble and Spitzer Space Telescope photometry to robustly identify galaxies in this epoch at F160W < 26.6. We use a detailed vetting procedure, including screening for persistence, stellar contamination, inclusion of ground-based imaging, and followup space-based imaging to build a robust sample of 11 candidate galaxies, three presented here for the first time. The inclusion of Spitzer/IRAC photometry in the selection process reduces contamination, and yields more robust redshift estimates than Hubble alone. We constrain the evolution of the rest-frame ultraviolet luminosity function via a new method of calculating the observed number densities without choosing a prior magnitude bin size. We find that the abundance at our brightest probed luminosities (M_UV=-22.3) is consistent with predictions from simulations which assume that galaxies in this epoch have gas depletion times at least as short as those in nearby starburst galaxies. Due to large Poisson and cosmic variance uncertainties we cannot conclusively rule out either a smooth evolution of the luminosity function continued from z=4-8, or an accelerate decline at z > 8. We calculate that the presence of seven galaxies in a single field (EGS) is an outlier at the 2-sigma significance level, implying the discovery of a significant overdensity. These scenarios will be imminently testable to high confidence within the first year of observations of the James Webb Space Telescope. △ Less

Submitted 19 November, 2021; v1 submitted 25 June, 2021; originally announced June 2021.

Comments: 47 figures, 25 pages, 10 tables, Accepted to the Astrophysical Journal. Photometric catalogs available via email request to lead author

Optical Spectroscopy of Dual Quasar Candidates from the Subaru HSC-SSP program

Authors: Shenli Tang, John D. Silverman, Xuheng Ding, Junyao Li, Khee-Gan Lee, Michael A. Strauss, Andy Goulding, Malte Schramm, Lalitwadee Kawinwanichakij, J. Xavier Prochaska, Joseph F. Hennawi, Masatoshi Imanishi, Kazushi Iwasawa, Yoshiki Toba, Issha Kayo, Masamune Oguri, Yoshiki Matsuoka, Kohei Ichikawa, Tilman Hartwig, Nobunari Kashikawa, Toshihiro Kawaguchi, Kotaro Kohno, Yuichi Matsuda, Tohru Nagao, Yoshiaki Ono , et al. (8 additional authors not shown)

Abstract: We report on a spectroscopic program to search for dual quasars using Subaru Hyper Suprime-Cam (HSC) images of SDSS quasars which represent an important stage during galaxy mergers. Using Subaru/FOCAS and Gemini-N/GMOS, we identify three new physically associated quasar pairs having projected separations less than 20 kpc, out of 26 observed candidates. These include the discovery of the highest re… ▽ More We report on a spectroscopic program to search for dual quasars using Subaru Hyper Suprime-Cam (HSC) images of SDSS quasars which represent an important stage during galaxy mergers. Using Subaru/FOCAS and Gemini-N/GMOS, we identify three new physically associated quasar pairs having projected separations less than 20 kpc, out of 26 observed candidates. These include the discovery of the highest redshift ($z=3.1$) quasar pair with a separation $<$ 10 kpc. Based on the sample acquired to date, the success rate of identifying physically associated dual quasars is $19\%$ when excluding stars based on their HSC colors. Using the full sample of six spectroscopically confirmed dual quasars, we find that the black holes in these systems have black hole masses ($M_{BH} \sim 10^{8-9}M_{\odot}$) similar to single SDSS quasars as well as their bolometric luminosities and Eddington ratios. We measure the stellar mass of their host galaxies based on 2D image decomposition of the five-band ($grizy$) optical emission and assess the mass relation between supermassive black holes (SMBHs) and their hosts. Dual SMBHs appear to have elevated masses relative to their host galaxies. Thus mergers may not necessarily align such systems onto the local mass relation, as suggested by the Horizon-AGN simulation. This study suggests that dual luminous quasars are triggered prior to the final coalescence of the two SMBHs, resulting in early mass growth of the black holes relative to their host galaxies. △ Less

Submitted 21 May, 2021; originally announced May 2021.

Comments: 36 pages, 21 figures, 4 tables

The Sizes of Quasar Host Galaxies with the Hyper Suprime-Cam Subaru Strategic Program

Authors: Junyao Li, John D. Silverman, Xuheng Ding, Michael A. Strauss, Andy Goulding, Simon Birrer, Hassen M. Yesuf, Yongquan Xue, Lalitwadee Kawinwanichakij, Yoshiki Matsuoka, Yoshiki Toba, Tohru Nagao, Malte Schramm, Kohei Inayoshi

Abstract: The relationship between quasars and their host galaxies provides clues on how supermassive black holes (SMBHs) and massive galaxies are jointly assembled. To elucidate this connection, we measure the structural and photometric properties of the host galaxies of ~5000 SDSS quasars at 0.2<z<1 using five-band (grizy) optical imaging from the Hyper Suprime-Cam Subaru Strategic Program. An automated a… ▽ More The relationship between quasars and their host galaxies provides clues on how supermassive black holes (SMBHs) and massive galaxies are jointly assembled. To elucidate this connection, we measure the structural and photometric properties of the host galaxies of ~5000 SDSS quasars at 0.2<z<1 using five-band (grizy) optical imaging from the Hyper Suprime-Cam Subaru Strategic Program. An automated analysis tool is used to forward-model the blended emission of the quasar as characterized by the point spread function and the underlying host galaxy as a two-dimensional Sersic profile. In agreement with previous studies, quasars are preferentially hosted by massive star-forming galaxies with disk-like light profiles. Furthermore, we find that the size distribution of quasar hosts is broad at a given stellar mass and the average values exhibit a size-stellar mass relation as seen with inactive galaxies. In contrast, the sizes of quasar hosts are more compact than inactive star-forming galaxies on average, but not as compact as quiescent galaxies of similar stellar masses. This is true irrespective of quasar properties including bolometric luminosity, Eddington ratio, and black hole mass. These results are consistent with a scenario in which galaxies are concurrently fueling a SMBH and building their stellar bulge from a centrally-concentrated gas reservoir. Alternatively, quasar hosts may be experiencing a compaction process in which stars from the disk and inflowing gas are responsible for growing the bulge. In addition, we confirm that the host galaxies of type-1 quasars have a bias of being closer towards face-on, suggesting that galactic-scale dust can contribute to obscuring the broad-line region. △ Less

Submitted 13 May, 2021; originally announced May 2021.

Comments: 34 pages, 22 figures, 3 tables, resubmitted to ApJ after addressing referee's comments; catalog will be available after the paper is published

The Shape and Scatter of The Galaxy Main Sequence for Massive Galaxies at Cosmic Noon

Authors: Sydney Sherman, Shardha Jogee, Jonathan Florez, Steven L. Finkelstein, Robin Ciardullo, Isak Wold, Matthew L. Stevans, Lalitwadee Kawinwanichakij, Casey Papovich, Caryl Gronwall

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

Submitted 8 May, 2021; originally announced May 2021.

Comments: 16 pages, 10 figures, Accepted for publication in MNRAS

The NEWFIRM HETDEX Survey: Photometric Catalog and a Conservative Sample of Massive Quiescent Galaxies at z=3-5 over 17.5 deg^2 in the SHELA Field

Authors: Matthew L. Stevans, Steven L. Finkelstein, Lalitwadee Kawinwanichakij, Isak Wold, Casey Papovich, Rachel S. Somerville, L. Y. Aaron Yung, Sydney Sherman, Robin Ciardullo, Romeel Dave, Jonathan Florez, Caryl Gronwall, Shardha Jogee

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

Submitted 26 March, 2021; originally announced March 2021.

Comments: 38 pages, accepted to the Astrophysical Journal. The full catalog is included as an ancillary file, and is be published with the journal article

Investigating The Growing Population of Massive Quiescent Galaxies at Cosmic Noon

Authors: Sydney Sherman, Shardha Jogee, Jonathan Florez, Matthew L. Stevans, Lalitwadee Kawinwanichakij, Isak Wold, Steven L. Finkelstein, Casey Papovich, Robin Ciardullo, Caryl Gronwall, Sofía A. Cora, Tomás Hough, Cristian A. Vega-Martínez

Abstract: We explore the buildup of quiescent galaxies using a sample of 28,469 massive ($M_\star \ge 10^{11}$M$_\odot$) galaxies at redshifts $1.5<z<3.0$, drawn from a 17.5 deg$^2$ area (0.33 Gpc$^3$ comoving volume at these redshifts). This allows for a robust study of the quiescent fraction as a function of mass at $1.5<z<3.0$ with a sample $\sim$40 times larger at log($M_{\star}$/$\rm M_{\odot}$)… ▽ More We explore the buildup of quiescent galaxies using a sample of 28,469 massive ($M_\star \ge 10^{11}$M$_\odot$) galaxies at redshifts $1.5<z<3.0$, drawn from a 17.5 deg$^2$ area (0.33 Gpc$^3$ comoving volume at these redshifts). This allows for a robust study of the quiescent fraction as a function of mass at $1.5<z<3.0$ with a sample $\sim$40 times larger at log($M_{\star}$/$\rm M_{\odot}$)$\ge11.5$ than previous studies. We derive the quiescent fraction using three methods: specific star-formation rate, distance from the main sequence, and UVJ color-color selection. All three methods give similar values at $1.5<z<2.0$, however the results differ by up to a factor of two at $2.0<z<3.0$. At redshifts $1.5 < z < 3.0$ the quiescent fraction increases as a function of stellar mass. By $z=2$, only 3.3 Gyr after the Big Bang, the universe has quenched $\sim$25% of $M_\star = 10^{11}$M$_\odot$ galaxies and $\sim$45% of $M_\star = 10^{12}$M$_\odot$ galaxies. We discuss physical mechanisms across a range of epochs and environments that could explain our results. We compare our results with predictions from hydrodynamical simulations SIMBA and IllustrisTNG and semi-analytic models (SAMs) SAG, SAGE, and Galacticus. The quiescent fraction from IllustrisTNG is higher than our empirical result by a factor of $2-5$, while those from SIMBA and the three SAMs are lower by a factor of $1.5-10$ at $1.5<z<3.0$. △ Less

Submitted 9 October, 2020; originally announced October 2020.

Comments: 22 pages, 14 figures, Accepted for publication in MNRAS

Exploring AGN and Star Formation Activity of Massive Galaxies at Cosmic Noon

Authors: Jonathan Florez, Shardha Jogee, Sydney Sherman, Matthew L. Stevans, Steven L. Finkelstein, Casey Papovich, Lalitwadee Kawinwanichakij, Robin Ciardullo, Caryl Gronwall, C. Megan Urry, Allison Kirkpatrick, Stephanie M. LaMassa, Tonima Tasnim Ananna, Isak Wold

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

Submitted 29 July, 2020; v1 submitted 23 July, 2020; originally announced July 2020.

Comments: 27 pages, 15 figures. Accepted for publication in MNRAS

On the (Lack of) Evolution of the Stellar Mass Function of Massive Galaxies from $z$=1.5 to 0.4

Authors: Lalitwadee Kawinwanichakij, Casey Papovich, Robin Ciardullo, Steven L. Finkelstein, Matthew L. Stevans, Isak G. Wold, Shardha Jogee, Sydney Sherman, Jonathan Florez, Caryl Gronwall

Abstract: We study the evolution in the number density of the highest mass galaxies over $0.4<z<1.5$ (covering 9 Gyr). We use the Spitzer/HETDEX Exploratory Large-Area (SHELA) Survey, which covers 17.5 $\mathrm{deg}^2$ with eight photometric bands spanning 0.3-4.5 $μ$m within the SDSS Stripe 82 field. This size produces the lowest counting uncertainties and cosmic variance yet for massive galaxies at… ▽ More We study the evolution in the number density of the highest mass galaxies over $0.4<z<1.5$ (covering 9 Gyr). We use the Spitzer/HETDEX Exploratory Large-Area (SHELA) Survey, which covers 17.5 $\mathrm{deg}^2$ with eight photometric bands spanning 0.3-4.5 $μ$m within the SDSS Stripe 82 field. This size produces the lowest counting uncertainties and cosmic variance yet for massive galaxies at $z\sim1.0$. We study the stellar mass function (SMF) for galaxies with $\log(M_\ast/M_\odot)>10.3$ using a forward-modeling method that fully accounts for statistical and systematic uncertainties on stellar mass. From $z$=0.4 to 1.5 the massive end of the SMF shows minimal evolution in its shape: the characteristic mass ($M^\ast$) evolves by less than 0.1 dex ($\pm$0.05 dex); the number density of galaxies with $\log (M_\ast/M_\odot) >11$ stays roughly constant at $\log (n/\mathrm{Mpc}^{-3})$ $\simeq$ $-$3.4 ($\pm$0.05), then declines to $\log (n/\mathrm{Mpc}^{-3})$=$-$3.7 ($\pm$0.05) at $z$=1.5. We discuss the uncertainties in the SMF, which are dominated by assumptions in the star formation history and details of stellar population synthesis models for stellar mass estimations. For quiescent galaxies, the data are consistent with no (or slight) evolution ($\lesssim0.1$ dex) in the characteristic mass nor number density from $z\sim 1.5$ to the present. This implies that any mass growth (presumably through "dry" mergers) of the quiescent massive galaxy population must balance the rate of mass losses from late-stage stellar evolution and the formation of quenching galaxies from the star-forming population. We provide a limit on this mass growth from $z=1.0$ to 0.4 of $ΔM_\ast/M_\ast\leq$ 45% (i.e., $\simeq0.16$ dex) for quiescent galaxies more massive than $10^{11}$ $M_\odot$. △ Less

Submitted 17 February, 2020; originally announced February 2020.

Comments: 39 pages, 23 figure, accepted for publication in ApJ

The Atacama Cosmology Telescope: SZ-based masses and dust emission from IR-selected cluster candidates in the SHELA survey

Authors: Brittany J. Fuzia, Lalitwadee Kawinwanichakij, Nicola Mehrtens, Simone Aiola, Nicholas Battaglia, Robin Ciardullo, Mark Devlin, Steven L. Finkelstein, Megan Gralla, Matt Hilton, Kevin M. Huffenberger, John P. Hughes, Shardha Jogee, Felipe A. Maldonado, Lyman A. Page, Casey Papovich, Bruce Partridge, Eli Rykoff, Neelima Sehgal, Cristobal Sifon, Suzanne T. Staggs, Edward Wollack

Abstract: We examine the stacked thermal Sunyaev-Zel\text{'}dovich (SZ) signals for a sample of galaxy cluster candidates from the Spitzer-HETDEX Exploratory Large Area (SHELA) Survey, which are identified in combined optical and infrared SHELA data using the redMaPPer algorithm. We separate the clusters into three richness bins, with average photometric redshifts ranging from 0.70 to 0.80. The richest bin… ▽ More We examine the stacked thermal Sunyaev-Zel\text{'}dovich (SZ) signals for a sample of galaxy cluster candidates from the Spitzer-HETDEX Exploratory Large Area (SHELA) Survey, which are identified in combined optical and infrared SHELA data using the redMaPPer algorithm. We separate the clusters into three richness bins, with average photometric redshifts ranging from 0.70 to 0.80. The richest bin shows a clear temperature decrement at 148 GHz in the Atacama Cosmology Telescope data, which we attribute to the SZ effect. All richness bins show an increment at 220 GHz, which we attribute to dust emission from cluster galaxies. We correct for dust emission using stacked profiles from Herschel Stripe 82 data, and allow for synchrotron emission using stacked profiles created by binning source fluxes from NVSS data. We see dust emission in all three richness bins, but can only confidently detect the SZ decrement in the highest richness bin, finding $M_{500}$ = $8.7^{+1.7}_{-1.3} \times 10^{13} M_\odot$. Neglecting the correction for dust depresses the inferred mass by 26 percent, indicating a partial fill-in of the SZ decrement from thermal dust and synchrotron emission by the cluster member galaxies. We compare our corrected SZ masses to two redMaPPer mass--richness scaling relations and find that the SZ mass is lower than predicted by the richness. We discuss possible explanations for this discrepancy, and note that the SHELA richnesses may differ from previous richness measurements due to the inclusion of IR data in redMaPPer. △ Less

Submitted 26 January, 2020; originally announced January 2020.

Comments: 14 pages, 7 figures. Submitted to MNRAS

Exploring the High-Mass End of the Stellar Mass Function of Star Forming Galaxies at Cosmic Noon

Authors: Sydney Sherman, Shardha Jogee, Jonathan Florez, Matthew L. Stevans, Lalitwadee Kawinwanichakij, Isak Wold, Steven L. Finkelstein, Casey Papovich, Viviana Acquaviva, Robin Ciardullo, Carly Gronwall, Zacharias Escalante

Abstract: We present the high-mass end of the galaxy stellar mass function using the largest sample to date (5,352) of star-forming galaxies with $M_{\star} > 10^{11} M_{\odot}$ at cosmic noon, $1.5 < z < 3.5$. This sample is uniformly selected across 17.2 deg$^2$ ($\sim$0.44 Gpc$^3$ comoving volume from $1.5 < z < 3.5$), mitigating the effects of cosmic variance and encompassing a wide range of environment… ▽ More We present the high-mass end of the galaxy stellar mass function using the largest sample to date (5,352) of star-forming galaxies with $M_{\star} > 10^{11} M_{\odot}$ at cosmic noon, $1.5 < z < 3.5$. This sample is uniformly selected across 17.2 deg$^2$ ($\sim$0.44 Gpc$^3$ comoving volume from $1.5 < z < 3.5$), mitigating the effects of cosmic variance and encompassing a wide range of environments. This area, a factor of 10 larger than previous studies, provides robust statistics at the high-mass end. Using multi-wavelength data in the Spitzer/HETDEX Exploratory Large Area (SHELA) footprint we find that the SHELA footprint star-forming galaxy stellar mass function is steeply declining at the high-mass end probing values as high as $\sim$$10^{-4}$ Mpc$^3$/dex and as low as $\sim$5$\times$$10^{-8}$ Mpc$^3$/dex across a stellar mass range of log($M_\star$/$M_\odot$) $\sim$ 11 - 12. We compare our empirical star-forming galaxy stellar mass function at the high mass end to three types of numerical models: hydrodynamical models from IllustrisTNG, abundance matching from the UniverseMachine, and three different semi-analytic models (SAMs; SAG, SAGE, GALACTICUS). At redshifts $1.5 < z < 3.5$ we find that results from IllustrisTNG and abundance matching models agree within a factor of $\sim$2 to 10, however the three SAMs strongly underestimate (up to a factor of 1,000) the number density of massive galaxies. We discuss the implications of these results for our understanding of galaxy evolution. △ Less

Submitted 15 November, 2019; originally announced November 2019.

Comments: 19 pages, 10 figures, Accepted for publication in MNRAS

ALMA 200-parsec Resolution Imaging of Smooth Cold Dusty Disks in Typical $z \sim 3$ Star-Forming Galaxies

Authors: W. Rujopakarn, E. Daddi, G. H. Rieke, A. Puglisi, M. Schramm, P. G. Pérez-González, G. E. Magdis, S. Alberts, F. Bournaud, D. Elbaz, M. Franco, L. Kawinwanichakij, K. Kohno, D. Narayanan, J. D. Silverman, T. Wang, C. C. Williams

Abstract: We present high-fidelity, 30 milliarcsecond (200-pc) resolution ALMA rest-frame 240 $μ$m observations of cold dust emission in three typical main-sequence star-forming galaxies (SFGs) at $z \sim 3$ in the Hubble Ultra-Deep Field (HUDF). The cold dust is distributed within the smooth disk-like central regions of star formation $1 - 3$ kpc in diameter, despite their complex and disturbed rest-frame… ▽ More We present high-fidelity, 30 milliarcsecond (200-pc) resolution ALMA rest-frame 240 $μ$m observations of cold dust emission in three typical main-sequence star-forming galaxies (SFGs) at $z \sim 3$ in the Hubble Ultra-Deep Field (HUDF). The cold dust is distributed within the smooth disk-like central regions of star formation $1 - 3$ kpc in diameter, despite their complex and disturbed rest-frame UV and optical morphologies. No dust substructures or clumps are seen down to $\simeq 1- 3$ $M_\odot$yr$^{-1}$ (1$σ$) per 200-pc beam. No dust emission is observed at the locations of UV-emitting clumps, which lie $\simeq 2-10$ kpc from the bulk of star formation. Clumpy substructures can contribute no more than $1-7$% of the total star formation in these galaxies (3$σ$ upper limits). The lack of star-forming substructures in our HUDF galaxies is to be contrasted with the multiple substructures characteristic of submillimeter-selected galaxies (SMGs) at the same cosmic epoch, particularly the far-IR-bright SMGs with similarly high-fidelity ALMA observations of Hodge et al. (2019). Individual star-forming substructures in these SMGs contain $\sim10-30$% of their total star formation. A substructure in these SMGs is often comparably bright in the far-infrared as (or in some cases brighter than) our typical SFGs, suggesting that these SMGs originate from a class of disruptive event involving multiple objects at the scale of our HUDF galaxies. The scale of the disruptive event found in our main-sequence SFGs, characterized by the lack of star-forming substructures at our resolution and sensitivity, could be less violent, e.g., gas-rich disk instability or minor mergers. △ Less

Submitted 4 August, 2019; v1 submitted 9 April, 2019; originally announced April 2019.

Comments: 22 pages, 7 figures; ApJ, accepted

The Spitzer-HETDEX Exploratory Large Area Survey II: Dark Energy Camera and Spitzer/IRAC Multiwavelength Catalog

Authors: Isak Wold, Lalitwadee Kawinwanichakij, Matthew Stevans, Steven Finkelstein, Casey Papovich, Yaswant Devarakonda, Robin Ciardullo, John Feldmeier, Jonathan Florez, Eric Gawiser, Caryl Gronwall, Shardha Jogee, Jennifer Marshall, Sydney Sherman, Heath Shipley, Rachel Somerville, Francisco Valdes, Gregory Zeimann

Abstract: We present the $ugriz$-band Dark Energy Camera (DECam) plus 3.6 and 4.5 $μ$m IRAC catalogs for the Spitzer/HETDEX Exploratory Large-Area (SHELA) survey. SHELA covers $\sim24$ deg$^{2}$ of the Sloan Digital Sky Survey (SDSS) Stripe 82 region, with seven bandpasses spanning a wavelength range of 0.35 to 4.5 $μ$m. SHELA falls within the footprint of the Hobby-Eberly Telescope Dark Energy Experiment (… ▽ More We present the $ugriz$-band Dark Energy Camera (DECam) plus 3.6 and 4.5 $μ$m IRAC catalogs for the Spitzer/HETDEX Exploratory Large-Area (SHELA) survey. SHELA covers $\sim24$ deg$^{2}$ of the Sloan Digital Sky Survey (SDSS) Stripe 82 region, with seven bandpasses spanning a wavelength range of 0.35 to 4.5 $μ$m. SHELA falls within the footprint of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX), which will provide spectroscopic redshifts for $\sim200{,}000$ Ly$α$ emitters at $1.9<z<3.5$ and also for $\sim200{,}000$ [OII] emitters at $z<0.5$. SHELA's deep, wide-area multiwavelength images combined with HETDEX's spectroscopic information, will facilitate many extragalactic studies, including measuring the evolution of galaxy stellar mass, halo mass, and environment from $1.5<z<3.5$. Here we present $riz$-band selected $ugriz$-band DECam catalogs that reach a $5σ$ depth of $\sim24.5$ AB mag (for point sources with an aperture that encloses $70\%$ of the total flux) and cover $17.5$ deg$^{2}$ of the overall SHELA field. We validate our DECam catalog by comparison to the DECam Legacy Survey (DECaLS) DR5 and the Dark Energy Survey (DES) DR1. We perform IRAC forced photometry with The Tractor image modeling code to measure 3.6 and 4.5 $μ$m fluxes for all objects within our DECam catalog. We demonstrate the utility of our catalog by computing galaxy number counts and estimating photometric redshifts. Our photometric redshifts recover the available $\left\langle z \right\rangle = 0.33 $ SDSS spectroscopic redshifts with a $1σ$ scatter in $Δz/(1 +z)$ of 0.04. △ Less

Submitted 11 December, 2018; originally announced December 2018.

Comments: 20 pages, 14 figures, accepted by ApJS

ZFOURGE: Using Composite Spectral Energy Distributions to Characterize Galaxy Populations at 1<z<4

Authors: Ben Forrest, Kim-Vy H. Tran, Adam Broussard, Jonathan H. Cohn, Robert C. Kennicutt Jr., Casey Papovich, Rebecca Allen, Michael Cowley, Karl Glazebrook, Glenn G. Kacprzak, Lalitwadee Kawinwanichakij, Themiya Nanayakkara, Brett Salmon, Lee R. Spitler, Caroline M. S. Straatman

Abstract: We investigate the properties of galaxies as they shut off star formation over the 4 billion years surrounding peak cosmic star formation. To do this we categorize $\sim7000$ galaxies from $1<z<4$ into $90$ groups based on the shape of their spectral energy distributions (SEDs) and build composite SEDs with $R\sim 50$ resolution. These composite SEDs show a variety of spectral shapes and also show… ▽ More We investigate the properties of galaxies as they shut off star formation over the 4 billion years surrounding peak cosmic star formation. To do this we categorize $\sim7000$ galaxies from $1<z<4$ into $90$ groups based on the shape of their spectral energy distributions (SEDs) and build composite SEDs with $R\sim 50$ resolution. These composite SEDs show a variety of spectral shapes and also show trends in parameters such as color, mass, star formation rate, and emission line equivalent width. Using emission line equivalent widths and strength of the 4000Å break, $D(4000)$, we categorize the composite SEDs into five classes: extreme emission line, star-forming, transitioning, post-starburst, and quiescent galaxies. The transitioning population of galaxies show modest H$α$ emission ($EW_{\rm REST}\sim40$Å) compared to more typical star-forming composite SEDs at $\log_{10}(M/M_\odot)\sim10.5$ ($EW_{\rm REST}\sim80$Å). Together with their smaller sizes (3 kpc vs. 4 kpc) and higher Sérsic indices (2.7 vs. 1.5), this indicates that morphological changes initiate before the cessation of star formation. The transitional group shows a strong increase of over one dex in number density from $z\sim3$ to $z\sim1$, similar to the growth in the quiescent population, while post-starburst galaxies become rarer at $z\lesssim1.5$. We calculate average quenching timescales of 1.6 Gyr at $z\sim1.5$ and 0.9 Gyr at $z\sim2.5$ and conclude that a fast quenching mechanism producing post-starbursts dominated the quenching of galaxies at early times, while a slower process has become more common since $z\sim2$. △ Less

Submitted 10 July, 2018; originally announced July 2018.

Comments: Accepted for publication in The Astrophysical Journal

Bridging Star-Forming Galaxy and AGN Ultraviolet Luminosity Functions at $z=4$ with the SHELA Wide-Field Survey

Authors: Matthew L. Stevans, Steven L. Finkelstein, Isak Wold, Lalitwadee Kawinwanichakij, Casey Papovich, Sydney Sherman, Robin Ciardullo, Jonathan Florez, Caryl Gronwall, Shardha Jogee, Rachel S. Somerville, L. Y. Aaron Yung

Abstract: We present a joint analysis of the rest-frame ultraviolet (UV) luminosity functions of continuum-selected star-forming galaxies and galaxies dominated by active galactic nuclei (AGNs) at $z \sim$ 4. These 3,740 $z \sim$ 4 galaxies are selected from broad-band imaging in nine photometric bands over 18 deg$^2$ in the \textit{Spitzer}/HETDEX Exploratory Large Area Survey (SHELA) field. The large area… ▽ More We present a joint analysis of the rest-frame ultraviolet (UV) luminosity functions of continuum-selected star-forming galaxies and galaxies dominated by active galactic nuclei (AGNs) at $z \sim$ 4. These 3,740 $z \sim$ 4 galaxies are selected from broad-band imaging in nine photometric bands over 18 deg$^2$ in the \textit{Spitzer}/HETDEX Exploratory Large Area Survey (SHELA) field. The large area and moderate depth of our survey provide a unique view of the intersection between the bright end of the galaxy UV luminosity function (M$_{AB}<-$22) and the faint end of the AGN UV luminosity function. We do not separate AGN-dominated galaxies from star-formation-dominated galaxies, but rather fit both luminosity functions simultaneously. These functions are best fit with a double power-law (DPL) for both the galaxy and AGN components, where the galaxy bright-end slope has a power-law index of $-3.80\pm0.10$, and the corresponding AGN faint-end slope is $α_{AGN} = -1.49^{+0.30}_{-0.21}$. We cannot rule out a Schechter-like exponential decline for the galaxy UV luminosity function, and in this scenario the AGN luminosity function has a steeper faint-end slope of $-2.08^{+0.18}_{-0.11}$. Comparison of our galaxy luminosity function results with a representative cosmological model of galaxy formation suggests that the molecular gas depletion time must be shorter, implying that star formation is more efficient in bright galaxies at $z=4$ than at the present day. If the galaxy luminosity function does indeed have a power-law shape at the bright end, the implied ionizing emissivity from AGNs is not inconsistent with previous observations. However, if the underlying galaxy distribution is Schechter, it implies a significantly higher ionizing emissivity from AGNs at this epoch. △ Less

Submitted 13 June, 2018; originally announced June 2018.

Comments: 25 pages, 13 figures, accepted for publication in ApJ

The Effects of Environment on the Evolution of the Galaxy Stellar Mass Function

Authors: Casey Papovich, Lalitwadee Kawinwanichakij, Ryan Quadri, Karl Glazebrook, Ivo Labbe, Kim-Vy Tran, Ben Forrest, Glenn G. Kacprzak, Lee R. Spitler, Caroline S. Straatman, Adam Tomczak

Abstract: We study the effects of galaxy environment on the evolution of the stellar-mass function (SMF) over 0.2 < z < 2.0 using the FourStar Galaxy Evolution (ZFOURGE) survey and NEWFIRM Medium-Band Survey (NMBS) down to the stellar-mass completeness limit, log M / Msun > 9.0 (9.5) at z = 1.0 (2.0). We compare the SMFs for quiescent and star-forming galaxies in the highest and lowest environments using a… ▽ More We study the effects of galaxy environment on the evolution of the stellar-mass function (SMF) over 0.2 < z < 2.0 using the FourStar Galaxy Evolution (ZFOURGE) survey and NEWFIRM Medium-Band Survey (NMBS) down to the stellar-mass completeness limit, log M / Msun > 9.0 (9.5) at z = 1.0 (2.0). We compare the SMFs for quiescent and star-forming galaxies in the highest and lowest environments using a density estimator based on the distance to the galaxies' third-nearest neighbors. For star-forming galaxies, at all redshifts there are only minor differences with environment in the shape of the SMF. For quiescent galaxies, the SMF in the lowest densities shows no evolution with redshift, other than an overall increase in number density (phi*) with time. This suggests that the stellar-mass dependence of quenching in relatively isolated galaxies is both universal and does not evolve strongly. While at z >~ 1.5 the SMF of quiescent galaxies is indistinguishable in the highest and lowest densities, at lower redshifts it shows a rapidly increasing number density of lower-mass galaxies, log M / Msun ~= 9-10. We argue this evolution can account for all the redshift evolution in the shape of the total quiescent-galaxy SMF. This evolution in the quiescent-galaxy SMF at higher redshift (z > 1) requires an environmental-quenching efficiency that decreases with decreasing stellar mass at 0.5 < z < 1.5 or it would overproduce the number of lower-mass quiescent galaxies in denser environments. This requires a dominant environment process such as starvation combined with rapid gas depletion and ejection at z > 0.5 - 1.0 for galaxies in our mass range. The efficiency of this process decreases with redshift allowing other processes (such as galaxy interactions and ram-pressure stripping) to become more important at later times, z < 0.5. △ Less

Submitted 15 January, 2018; originally announced January 2018.

Comments: Accepted for publication in the Astrophysical Journal, 19 pages

Effect of local environment and stellar mass on galaxy quenching and morphology at $0.5<z<2.0$

Authors: Lalitwadee Kawinwanichakij, Casey Papovich, Ryan F. Quadri, Karl Glazebrook, Glenn G. Kacprzak, Rebecca J. Allen, Eric F. Bell, Darren J. Croton, Avishai Dekel, Henry C. Ferguson, Ben Forrest, Norman A. Grogin, Yicheng Guo, Dale D. Kocevski, Anton M. Koekemoer, Ivo Labbé, Ray A. Lucas, Themiya Nanayakkara, Lee R. Spitler, Caroline M. S. Straatman, Kim-Vy H. Tran, Adam Tomczak, Pieter van Dokkum

Abstract: We study galactic star-formation activity as a function of environment and stellar mass over 0.5<z<2.0 using the FourStar Galaxy Evolution (ZFOURGE) survey. We estimate the galaxy environment using a Bayesian-motivated measure of the distance to the third nearest neighbor for galaxies to the stellar mass completeness of our survey, $\log(M/M_\odot)>9 (9.5)$ at z=1.3 (2.0). This method, when applie… ▽ More We study galactic star-formation activity as a function of environment and stellar mass over 0.5<z<2.0 using the FourStar Galaxy Evolution (ZFOURGE) survey. We estimate the galaxy environment using a Bayesian-motivated measure of the distance to the third nearest neighbor for galaxies to the stellar mass completeness of our survey, $\log(M/M_\odot)>9 (9.5)$ at z=1.3 (2.0). This method, when applied to a mock catalog with the photometric-redshift precision ($σ_z / (1+z) \lesssim 0.02$), recovers galaxies in low- and high-density environments accurately. We quantify the environmental quenching efficiency, and show that at z> 0.5 it depends on galaxy stellar mass, demonstrating that the effects of quenching related to (stellar) mass and environment are not separable. In high-density environments, the mass and environmental quenching efficiencies are comparable for massive galaxies ($\log (M/M_\odot)\gtrsim$ 10.5) at all redshifts. For lower mass galaxies ($\log (M/M)_\odot) \lesssim$ 10), the environmental quenching efficiency is very low at $z\gtrsim$ 1.5, but increases rapidly with decreasing redshift. Environmental quenching can account for nearly all quiescent lower mass galaxies ($\log(M/M_\odot) \sim$ 9-10), which appear primarily at $z\lesssim$ 1.0. The morphologies of lower mass quiescent galaxies are inconsistent with those expected of recently quenched star-forming galaxies. Some environmental process must transform the morphologies on similar timescales as the environmental quenching itself. The evolution of the environmental quenching favors models that combine gas starvation (as galaxies become satellites) with gas exhaustion through star-formation and outflows ("overconsumption"), and additional processes such as galaxy interactions, tidal stripping and disk fading to account for the morphological differences between the quiescent and star-forming galaxy populations. △ Less

Submitted 20 September, 2017; v1 submitted 12 June, 2017; originally announced June 2017.

Comments: 29 pages, 15 figure, accepted for publication in ApJ

The FourStar Galaxy Evolution Survey (ZFOURGE): ultraviolet to far-infrared catalogs, medium-bandwidth photometric redshifts with improved accuracy, stellar masses, and confirmation of quiescent galaxies to z~3.5

Authors: Caroline M. S. Straatman, Lee R. Spitler, Ryan F. Quadri, Ivo Labbe, Karl Glazebrook, S. Eric Persson, Casey Papovich, Kim-Vy H. Tran, Gabriel B. Brammer, Michael Cowley, Adam Tomczak, Themiya Nanayakkara, Leo Alcorn, Rebecca Allen, Adam Broussard, Pieter van Dokkum, Ben Forrest, Josha van Houdt, Glenn G. Kacprzak, Lalitwadee Kawinwanichakij, Daniel D. Kelson, Janice Lee, Patrick J. McCarthy, Nicola Mehrtens, Andrew Monson , et al. (4 additional authors not shown)

Abstract: The FourStar galaxy evolution survey (ZFOURGE) is a 45 night legacy program with the FourStar near-infrared camera on Magellan and one of the most sensitive surveys to date. ZFOURGE covers a total of $400\ \mathrm{arcmin}^2$ in cosmic fields CDFS, COSMOS and UDS, overlapping CANDELS. We present photometric catalogs comprising $>70,000$ galaxies, selected from ultradeep $K_s$-band detection images… ▽ More The FourStar galaxy evolution survey (ZFOURGE) is a 45 night legacy program with the FourStar near-infrared camera on Magellan and one of the most sensitive surveys to date. ZFOURGE covers a total of $400\ \mathrm{arcmin}^2$ in cosmic fields CDFS, COSMOS and UDS, overlapping CANDELS. We present photometric catalogs comprising $>70,000$ galaxies, selected from ultradeep $K_s$-band detection images ($25.5-26.5$ AB mag, $5σ$, total), and $>80\%$ complete to $K_s<25.3-25.9$ AB. We use 5 near-IR medium-bandwidth filters ($J_1,J_2,J_3,H_s,H_l$) as well as broad-band $K_s$ at $1.05\ - 2.16\ μm$ to $25-26$ AB at a seeing of $\sim0.5$". Each field has ancillary imaging in $26-40$ filters at $0.3-8\ μm$. We derive photometric redshifts and stellar population properties. Comparing with spectroscopic redshifts indicates a photometric redshift uncertainty $σ_z={0.010,0.009}$, and 0.011 in CDFS, COSMOS, and UDS. As spectroscopic samples are often biased towards bright and blue sources, we also inspect the photometric redshift differences between close pairs of galaxies, finding $σ_{z,pairs}= 0.01-0.02$ at $1<z<2.5$. We quantify how $σ_{z,pairs}$ depends on redshift, magnitude, SED type, and the inclusion of FourStar medium bands. $σ_{z,pairs}$ is smallest for bright, blue star-forming samples, while red star-forming galaxies have the worst $σ_{z,pairs}$. Including FourStar medium bands reduces $σ_{z,pairs}$ by 50\% at $1.5<z<2.5$. We calculate SFRs based on ultraviolet and ultradeep far-IR $Spitzer$/MIPS and Herschel/PACS data. We derive rest-frame $U-V$ and $V-J$ colors, and illustrate how these correlate with specific SFR and dust emission to $z=3.5$. We confirm the existence of quiescent galaxies at $z\sim3$, demonstrating their SFRs are suppressed by $>\times15$. △ Less

Submitted 30 August, 2016; v1 submitted 26 August, 2016; originally announced August 2016.

Comments: 39 pages, 35 figures, accepted for publication in ApJ, full data release available at http://zfourge.tamu.edu/

The Spitzer-HETDEX Exploratory Large-Area Survey

Authors: Casey Papovich, H. V. Shipley, N. Mehrtens, C. Lanham, M. Lacy, R. Ciardullo, S. L. Finkelstein, R. Bassett, P Behroozi, G. A. Blanc, R. S. de Jong, D. L. DePoy, N. Drory, E. Gawiser, K. Gebhardt, C. Gronwall, G. J. Hill, U. Hopp, S. Jogee, L. Kawinwanichakij, J. L. Marshall, E. McLinden, E. Mentuch Cooper, R. S. Somerville, M. Steinmetz , et al. (5 additional authors not shown)

Abstract: We present post-cryogenic Spitzer imaging at 3.6 and 4.5 micron with the Infrared Array Camera (IRAC) of the Spitzer/HETDEX Exploratory Large-Area (SHELA) survey. SHELA covers $\sim$deg$^2$ of the Sloan Digital Sky Survey "Stripe 82" region, and falls within the footprints of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) and the Dark Energy Survey. The HETDEX blind R $\sim$ 800 spectr… ▽ More We present post-cryogenic Spitzer imaging at 3.6 and 4.5 micron with the Infrared Array Camera (IRAC) of the Spitzer/HETDEX Exploratory Large-Area (SHELA) survey. SHELA covers $\sim$deg$^2$ of the Sloan Digital Sky Survey "Stripe 82" region, and falls within the footprints of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) and the Dark Energy Survey. The HETDEX blind R $\sim$ 800 spectroscopy will produce $\sim$ 200,000 redshifts from the Lyman-$α$ emission for galaxies in the range 1.9 < z < 3.5, and an additional $\sim$200,000 redshifts from the [OII] emission for galaxies at z < 0.5. When combined with deep ugriz images from the Dark Energy Camera, K-band images from NEWFIRM, and other ancillary data, the IRAC photometry from Spitzer will enable a broad range of scientific studies of the relationship between structure formation, galaxy stellar mass, halo mass, AGN, and environment over a co-moving volume of $\sim$0.5 Gpc$^3$ at 1.9 < z < 3.5. Here, we discuss the properties of the SHELA IRAC dataset, including the data acquisition, reduction, validation, and source catalogs. Our tests show the images and catalogs are 80% (50%) complete to limiting magnitudes of 22.0 (22.6) AB mag in the detection image, which is constructed from the weighted sum of the IRAC 3.6 and 4.5 micron images. The catalogs reach limiting sensitivities of 1.1 $μ$Jy at both 3.6 and 4.5 micron (1$σ$, for R=2 arcsec circular apertures). As a demonstration of science, we present IRAC number counts, examples of highly temporally variable sources, and galaxy surface density profiles of rich galaxy clusters. In the spirit of Spitzer Exploratory programs we provide all images and catalogs as part of the publication. △ Less

Submitted 17 March, 2016; originally announced March 2016.

Comments: Accepted for publication in ApJS. 30 pages, emulateapj format. Many figures. All data products, including images and catalogs available http://irsa.ipac.caltech.edu/data/SPITZER/SHELA

UV to IR Luminosities and Dust Attenuation Determined from ~4000 K-Selected Galaxies at 1<z<3 in the ZFOURGE Survey

Authors: Ben Forrest, Kim-Vy H. Tran, Adam R. Tomczak, Adam Broussard, Ivo Labbé, Casey Papovich, Mariska Kriek, Rebecca J. Allen, Michael Cowley, Mark Dickinson, Karl Glazebrook, Josha van Houdt, Hanae Inami, Glenn G. Kacprzak, Lalitwadee Kawinwanichakij, Daniel Kelson, Patrick J. McCarthy, Andrew Monson, Glenn Morrison, Themiya Nanayakkara, S. Eric Persson, Ryan F. Quadri, Lee R. Spitler, Caroline Straatman, Vithal Tilvi

Abstract: We build a set of composite galaxy SEDs by de-redshifting and scaling multi-wavelength photometry from galaxies in the ZFOURGE survey, covering the CDFS, COSMOS, and UDS fields. From a sample of ~4000 K_s-band selected galaxies, we define 38 composite galaxy SEDs that yield continuous low-resolution spectra (R~45) over the rest-frame range 0.1-4 um. Additionally, we include far infrared photometry… ▽ More We build a set of composite galaxy SEDs by de-redshifting and scaling multi-wavelength photometry from galaxies in the ZFOURGE survey, covering the CDFS, COSMOS, and UDS fields. From a sample of ~4000 K_s-band selected galaxies, we define 38 composite galaxy SEDs that yield continuous low-resolution spectra (R~45) over the rest-frame range 0.1-4 um. Additionally, we include far infrared photometry from the Spitzer Space Telescope and the Herschel Space Observatory to characterize the infrared properties of our diverse set of composite SEDs. From these composite SEDs we analyze the rest-frame UVJ colors, as well as the ratio of IR to UV light (IRX) and the UV slope ($β$) in the IRX$-β$ dust relation at 1<z<3. Blue star-forming composite SEDs show IRX and $β$ values consistent with local relations; dusty star-forming galaxies have considerable scatter, as found for local IR bright sources, but on average appear bluer than expected for their IR fluxes. We measure a tight linear relation between rest-frame UVJ colors and dust attenuation for star-forming composites, providing a direct method for estimating dust content from either (U-V) or (V-J) rest-frame colors for star-forming galaxies at intermediate redshifts. △ Less

Submitted 2 February, 2016; originally announced February 2016.

Comments: 6 pages, 4 figures. Accepted for publication in ApJ Letters

Satellite Quenching and Galactic Conformity at 0.3 < z < 2.5

Authors: Lalitwadee Kawinwanichakij, Ryan F. Quadri, Casey Papovich, Glenn G. Kacprzak, Ivo Labbé, Lee R. Spitler, Caroline Straatman, Kim-Vy Tran, Rebecca Allen, Peter S. Behroozi, Michael Cowley, Avishai Dekel, Karl Glazebrook, William G. Hartley, Daniel D. Kelson, David C. Koo, Seong-Kook Lee, Yu Lu, Themiya Nanayakkara, Eric Persson, Joel R. Primack, Vithal Tilvi, Adam R. Tomczak, Pieter van Dokkum

Abstract: We measure the evolution of the quiescent fraction and quenching efficiency of satellites around star-forming and quiescent central galaxies with stellar mass $\log(M_{\mathrm{cen}}/M_{\odot})>10.5$ at $0.3<z<2.5$. We combine imaging from three deep near-infrared-selected surveys (ZFOURGE/CANDELS, UDS, and UltraVISTA), which allows us to select a stellar-mass complete sample of satellites with… ▽ More We measure the evolution of the quiescent fraction and quenching efficiency of satellites around star-forming and quiescent central galaxies with stellar mass $\log(M_{\mathrm{cen}}/M_{\odot})>10.5$ at $0.3<z<2.5$. We combine imaging from three deep near-infrared-selected surveys (ZFOURGE/CANDELS, UDS, and UltraVISTA), which allows us to select a stellar-mass complete sample of satellites with $\log(M_{\mathrm{sat}}/M_{\odot})>9.3$. Satellites for both star-forming and quiescent central galaxies have higher quiescent fractions compared to field galaxies matched in stellar mass at all redshifts. We also observe "galactic conformity": satellites around quiescent centrals are more likely to be quenched compared to the satellites around star-forming centrals. In our sample, this conformity signal is significant at $\gtrsim3σ$ for $0.6<z<1.6$, whereas it is only weakly significant at $0.3<z<0.6$ and $1.6<z<2.5$. Therefore, conformity (and therefore satellite quenching) has been present for a significant fraction of the age of the universe. The satellite quenching efficiency increases with increasing stellar mass of the central, but does not appear to depend on the stellar mass of the satellite to the mass limit of our sample. When we compare the satellite quenching efficiency of star-forming centrals with stellar masses 0.2 dex higher than quiescent centrals (which should account for any difference in halo mass), the conformity signal decreases, but remains statistically significant at $0.6<z<0.9$. This is evidence that satellite quenching is connected to the star-formation properties of the central as well as to the mass of the halo. We discuss physical effects that may contribute to galactic conformity, and emphasize that they must allow for continued star-formation in the central galaxy even as the satellites are quenched. △ Less

Submitted 9 November, 2015; originally announced November 2015.

Comments: Submitted to ApJ after first round of referee comments, 22 pages, 12 figures

The sizes of massive quiescent and star forming galaxies at z~4 with ZFOURGE and CANDELS

Authors: Caroline M. S. Straatman, Ivo Labbe, Lee R. Spitler, Karl Glazebrook, Adam Tomczak, Rebecca Allen, Gabriel B. Brammer, Michael Cowley, Pieter van Dokkum, Glenn G. Kacprzak, Lalit Kawinwanichakij, Nicola Mehrtens, Themiya Nanayakkara, Casey Papovich, S. Eric Persson, Ryan F. Quadri, Glen Rees, Vithal Tilvi, Kim-Vy Tran, Katherine E. Whitaker

Abstract: We study the rest-frame ultra-violet sizes of massive (~0.8 x 10^11 M_Sun) galaxies at 3.4<z<4.2, selected from the FourStar Galaxy Evolution Survey (ZFOURGE), by fitting single Sersic profiles to HST/WFC3/F160W images from the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS). Massive quiescent galaxies are very compact, with a median circularized half-light radius r_e = 0.… ▽ More We study the rest-frame ultra-violet sizes of massive (~0.8 x 10^11 M_Sun) galaxies at 3.4<z<4.2, selected from the FourStar Galaxy Evolution Survey (ZFOURGE), by fitting single Sersic profiles to HST/WFC3/F160W images from the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS). Massive quiescent galaxies are very compact, with a median circularized half-light radius r_e = 0.63 +/- 0.18 kpc. Removing 5/16 (31%) sources with signs of AGN activity does not change the result. Star-forming galaxies have r_e = 2.0 +/- 0.60 kpc, 3.2 +/- 1.3 x larger than quiescent galaxies. Quiescent galaxies at z~4 are on average 6.0 +\- 0.17 x smaller than at z~0 and 1.9 +/- 0.7 x smaller than at z~2. Star-forming galaxies of the same stellar mass are 2.4 +/- 0.7 x smaller than at z~0. Overall, the size evolution at 0<z<4 is well described by a powerlaw, with r_e = 5.08 +/- 0.28 (1+z)^(-1.44+/-0.08) kpc for quiescent and r_e = 6.02 +/- 0.28 (1+z)^(-0.72+/-0.05) kpc for star-forming galaxies. Compact star-forming galaxies are rare in our sample: we find only 1/14 (7%) with r_e / (M / 10^11 M_Sun)^0.75 < 1.5, whereas 13/16 (81%) of the quiescent galaxies is compact. The number density of compact quiescent galaxies at z~4 is 1.8 +/- 0.8 x 10^-5 Mpc^-3 and increases rapidly, by >5 x, between 2<z<4. The paucity of compact star-forming galaxies at z~4 and their large rest-frame ultra-violet median sizes suggest that the formation phase of compact cores is very short and/or highly dust obscured. △ Less

Submitted 3 June, 2015; originally announced June 2015.

Comments: 9 pages, 4 figures. Accepted for publication in ApjL

The distribution of satellites around massive galaxies at 1<z<3 in ZFOURGE/CANDELS: dependence on star formation activity

Authors: Lalitwadee Kawinwanichakij, Casey Papovich, Ryan F. Quadri, Kim-Vy H. Tran, Lee R. Spitler, Glenn G. Kacprzak, Ivo Labbe, Caroline M. S. Straatman, Karl Glazebrook, Rebecca Allen, Michael Cowley, Romeel Davé, Avishai Dekel, Henry C. Ferguson, W. G Hartley, Anton M. Koekemoer, David C. Koo, Yu Lu, Nicola Mehrtens, Themiya Nanayakkara, S. Eric Persson, Glen Rees, Brett Salmon, Vithal Tilvi, Adam R. Tomczak , et al. (1 additional authors not shown)

Abstract: We study the statistical distribution of satellites around star-forming and quiescent central galaxies at 1<z<3 using imaging from the FourStar Galaxy Evolution Survey (ZFOURGE) and the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). The deep near-IR data select satellites down to $\log(M/M_\odot)>9$ at z<3. The radial satellite distribution around centrals is consistent with a… ▽ More We study the statistical distribution of satellites around star-forming and quiescent central galaxies at 1<z<3 using imaging from the FourStar Galaxy Evolution Survey (ZFOURGE) and the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). The deep near-IR data select satellites down to $\log(M/M_\odot)>9$ at z<3. The radial satellite distribution around centrals is consistent with a projected NFW profile. Massive quiescent centrals, $\log(M/M_\odot)>10.78$, have $\sim$2 times the number of satellites compared to star-forming centrals with a significance of 2.7$σ$ even after accounting for differences in the centrals' stellar-mass distributions. We find no statistical difference in the satellite distributions of intermediate-mass quiescent and star-forming centrals, $10.48<\log(M/M_\odot)<10.78$. Comparing to the Guo2011 semi-analytic model, the excess number of satellites indicates that quiescent centrals have halo masses 0.3 dex larger than star-forming centrals, even when the stellar-mass distributions are fixed. We use a simple toy model that relates halo mass and quenching, which roughly reproduces the observed quenched fractions and the differences in halo mass between star-forming and quenched galaxies only if galaxies have a quenching probability that increases with halo mass from $\sim$0 for $\log(M_h/M_\odot)\sim$11 to $\sim$1 for $\log(M_h/M_\odot)\sim$13.5. A single halo-mass quenching threshold is unable to reproduce the quiescent fraction and satellite distribution of centrals. Therefore, while halo quenching may be an important mechanism, it is unlikely to be the only factor driving quenching. It remains unclear why a high fraction of centrals remain star-forming even in relatively massive halos. △ Less

Submitted 7 July, 2014; v1 submitted 23 June, 2014; originally announced June 2014.

Comments: 19 pages, 17 figures, accepted by ApJ. Information on ZFOURGE can be found at http://zfourge.tamu.edu

Journal ref: 2014 ApJ 792 103

Exploring the z=3-4 massive galaxy population with ZFOURGE: the prevalence of dusty and quiescent galaxies

Authors: Lee R. Spitler, Caroline M. S. Straatman, Ivo Labbe, Karl Glazebrook, Kim-Vy H. Tran, Glenn G. Kacprzak, Ryan F. Quadri, Casey Papovich, S. Eric Persson, Pieter van Dokkum, Rebecca Allen, Lalitwadee Kawinwanichakij, Daniel D. Kelson, Patrick J. McCarthy, Nicola Mehrtens, Andrew J. Monson, Themiya Nanayakkara, Glen Rees, Vithal Tilvi, Adam R. Tomczak

Abstract: Our understanding of the redshift $z>3$ galaxy population relies largely on samples selected using the popular "dropout" technique, typically consisting of UV-bright galaxies with blue colors and prominent Lyman breaks. As it is currently unknown if these galaxies are representative of the massive galaxy population, we here use the FourStar Galaxy Evolution (ZFOURGE) Survey to create a stellar mas… ▽ More Our understanding of the redshift $z>3$ galaxy population relies largely on samples selected using the popular "dropout" technique, typically consisting of UV-bright galaxies with blue colors and prominent Lyman breaks. As it is currently unknown if these galaxies are representative of the massive galaxy population, we here use the FourStar Galaxy Evolution (ZFOURGE) Survey to create a stellar mass-limited sample at $z=3-4$. Uniquely, ZFOURGE uses deep near-infrared medium-bandwidth filters to derive accurate photometric redshifts and stellar population properties. The mass-complete sample consists of 57 galaxies with log M $>10.6$, reaching below $M^{\star}$ at $z=3-4$. On average, the massive $z=3-4$ galaxies are extremely faint in the observed optical with median $R_{tot}^{AB}=27.48\pm0.41$ (restframe $M_{1700}=-18.05\pm0.37$). They lie far below the UV luminosity-stellar mass relation for Lyman break galaxies and are about $\sim100\times$ fainter at the same mass. The massive galaxies are red ($R-Ks_{AB}=3.9\pm0.2$; restframe UV-slope $β=-0.2\pm0.3$) likely from dust or old stellar ages. We classify the galaxy SEDs by their restframe $U-V$ and $V-J$ colors and find a diverse population: $46^{+6+10}_{-6-17}$% of the massive galaxies are quiescent, $54^{+8+17}_{-8-10}$% are dusty star-forming galaxies, and only $14^{+3+10}_{-3-4}$% resemble luminous blue star forming Lyman break galaxies. This study clearly demonstrates an inherent diversity among massive galaxies at higher redshift than previously known. Furthermore,we uncover a reservoir of dusty star-forming galaxies with $4\times$ lower specific star-formation rates compared to submillimeter-selected starbursts at $z>3$. With $5\times$ higher numbers, the dusty galaxies may represent a more typical mode of star formation compared to submillimeter-bright starbursts. △ Less

Submitted 5 May, 2014; originally announced May 2014.

Comments: Accepted for publication in ApJ Letters. See http://zfourge.tamu.edu for more details about ZFOURGE

Journal ref: 2014 ApJ 787 L36

A substantial population of massive quiescent galaxies at z~4 from ZFOURGE

Authors: Caroline M. S. Straatman, Ivo Labbé, Lee R. Spitler, Rebecca Allen, Bruno Altieri, Gabriel B. Brammer, Mark Dickinson, Pieter van Dokkum, Hanae Inami, Karl Glazebrook, Glenn G. Kacprzak, Lalit Kawinwanichakij, Daniel D. Kelson, Patrick J. McCarthy, Nicola Mehrtens, Andy Monson, David Murphy, Casey Papovich, S. Eric Persson, Ryan Quadri, Glen Rees, Adam Tomczak, Kim-Vy H. Tran, Vithal Tilvi

Abstract: We report the likely identification of a substantial population of massive M~10^11M_Sun galaxies at z~4 with suppressed star formation rates (SFRs), selected on rest-frame optical to near-IR colors from the FourStar Galaxy Evolution Survey. The observed spectral energy distributions show pronounced breaks, sampled by a set of near-IR medium-bandwidth filters, resulting in tightly constrained photo… ▽ More We report the likely identification of a substantial population of massive M~10^11M_Sun galaxies at z~4 with suppressed star formation rates (SFRs), selected on rest-frame optical to near-IR colors from the FourStar Galaxy Evolution Survey. The observed spectral energy distributions show pronounced breaks, sampled by a set of near-IR medium-bandwidth filters, resulting in tightly constrained photometric redshifts. Fitting stellar population models suggests large Balmer/4000AA breaks, relatively old stellar populations, large stellar masses and low SFRs, with a median specific SFR of 2.9+/-1.8 x 10^-11/yr. Ultradeep Herschel/PACS 100micron, 160micron and Spitzer/MIPS 24micron data reveal no dust-obscured SFR activity for 15/19 (79%) galaxies. Two far-IR detected galaxies are obscured QSOs. Stacking the far-IR undetected galaxies yields no detection, consistent with the SED fit, indicating independently that the average specific SFR is at least 10x smaller than of typical star-forming galaxies at z~4. Assuming all far-IR undetected galaxies are indeed quiescent, the volume density is 1.8+/-0.7 x 10^-5Mpc^-3 to a limit of log10M/M_Sun>10.6, which is 10x and 80x lower than at z = 2 and z = 0.1. They comprise a remarkably high fraction (~35%) of z~4 massive galaxies, suggesting that suppression of star formation was efficient even at very high redshift. Given the average stellar age of 0.8Gyr and stellar mass of 0.8x10^11M_Sun, the galaxies likely started forming stars before z =5, with SFRs well in excess of 100M_Sun/yr, far exceeding that of similarly abundant UV-bright galaxies at z>4. This suggests that most of the star-formation in the progenitors of quiescent z~4 galaxies was obscured by dust. △ Less

Submitted 13 February, 2014; v1 submitted 17 December, 2013; originally announced December 2013.

Comments: 8 pages, 4 figures, 1 table, accepted by ApjL

Galaxy Stellar Mass Functions from ZFOURGE/CANDELS: An Excess of Low-Mass Galaxies Since z=2 and the Rapid Buildup of Quiescent Galaxies

Authors: Adam R. Tomczak, Ryan F. Quadri, Kim-Vy H. Tran, Ivo Labbe, Caroline M. S. Straatman, Casey Papovich, Karl Glazebrook, Rebecca Allen, Gabriel B. Brammer, Glenn G. Kacprzak, Lalitwadee Kawinwanichakij, Daniel D. Kelson, Patrick J. McCarthy, Nicola Mehrtens, Andrew J. Monson, S. Eric Persson, Lee R. Spitler, Vithal Tilvi, Pieter van Dokkum

Abstract: Using observations from the FourStar Galaxy Evolution Survey (ZFOURGE), we obtain the deepest measurements to date of the galaxy stellar mass function at 0.5 < z < 2.5. ZFOURGE provides well-constrained photometric redshifts made possible through deep medium-bandwidth imaging at 1-2um . We combine this with HST imaging from the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS), al… ▽ More Using observations from the FourStar Galaxy Evolution Survey (ZFOURGE), we obtain the deepest measurements to date of the galaxy stellar mass function at 0.5 < z < 2.5. ZFOURGE provides well-constrained photometric redshifts made possible through deep medium-bandwidth imaging at 1-2um . We combine this with HST imaging from the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS), allowing for the efficient selection of both blue and red galaxies down to stellar masses ~10^9.5 Msol at z ~ 2.5. The total surveyed area is 316 arcmin^2 distributed over three independent fields. We supplement these data with the wider and shallower NEWFIRM Medium-Band Survey (NMBS) to provide stronger constraints at high masses. Several studies at z<=1 have revealed a steepening of the slope at the low-mass end of the stellar mass function (SMF), leading to an upturn at masses <10^10 Msol that is not well-described by a standard single-Schechter function. We find evidence that this feature extends to at least z ~ 2, and that it can be found in both the star-forming and quiescent populations individually. The characteristic mass (M*) and slope at the lowest masses (alpha) of a double-Schechter function fit to the SMF stay roughly constant at Log(M/Msol) ~ 10.65 and ~-1.5 respectively. The SMF of star-forming galaxies has evolved primarily in normalization, while the change in shape is relatively minor. This is not the case for quiescent galaxies: the depth of our imaging allows us to show for the first time significantly more evolution at Log(M/Msol) < 10.5 than at higher masses. We find that the total mass density (down to 10^9 Msol) in star-forming galaxies has increased by a factor of ~2.2 since z ~ 2.5, whereas in quiescent galaxies it has increased by a factor of ~12 . △ Less

Submitted 3 February, 2014; v1 submitted 23 September, 2013; originally announced September 2013.

Comments: 17 pages, 13 figures. Accepted for publication in the ApJ