Showing 1–50 of 73 results for author: Straatman, C

Spatially resolved mock observations of stellar kinematics: full radiative transfer treatment of simulated galaxies

Authors: Daniela Barrientos Acevedo, Arjen van der Wel, Maarten Baes, Robert J. J Grand, Anand Utsav Kapoor, Peter Camps, Anna de Graaff, Caroline M. S. Straatman, Rachel Bezanson

Abstract: We present a framework to build realistic mock spectroscopic observations for state-of-the-art hydrodynamical simulations, using high spectral resolution stellar population models and full radiative transfer treatment with SKIRT. As a first application we generate stellar continuum mock observations for the Auriga cosmological zoom simulations emulating integral-field observations from the SAMI Ga… ▽ More We present a framework to build realistic mock spectroscopic observations for state-of-the-art hydrodynamical simulations, using high spectral resolution stellar population models and full radiative transfer treatment with SKIRT. As a first application we generate stellar continuum mock observations for the Auriga cosmological zoom simulations emulating integral-field observations from the SAMI Galaxy survey. We perform spectral fitting on our synthetic cubes and compute the resulting rotation velocity ($V_{\rm{rot}}$) and velocity dispersion within 1$R_{\text{e}}$ ($σ_{\text{e}}$) for a sub-set of the Auriga sample. We find that the kinematics produced by Auriga are in good agreement with the observations from the SAMI Galaxy survey after taking into account the effects of dust and the systematics produced by the observation limitations. We also explore the effects of seeing convolution, inclination, and attenuation on the line-of-sight velocity distribution. For highly inclined galaxies, these effects can lead to an artificial decrease in the measured $V/σ$ by nearly a factor two (after inclination correction). We also demonstrate the utility of our method for high-redshift galaxies by emulating spatially resolved continuum spectra from the LEGA-C survey and, looking forward, E-ELT HARMONI. Our framework represents a crucial link between the ground truth for stellar populations and kinematics in simulations and the observed stellar continuum observations at low and high redshift. △ Less

Submitted 13 June, 2023; originally announced June 2023.

Comments: 17 pages, 11 figures. Accepted for publication in MNRAS

Galaxy kinematics and mass estimates at z~1 from ionised gas and stars

Authors: Hannah Übler, Natascha M. Förster Schreiber, Arjen van der Wel, Rachel Bezanson, Sedona H. Price, Francesco D'Eugenio, Emily Wisnioski, Reinhard Genzel, Linda J. Tacconi, Stijn Wuyts, Thorsten Naab, Dieter Lutz, Caroline M. S. Straatman, T. Taro Shimizu, Ric Davies, Daizhong Liu, J. Trevor Mendel

Abstract: We compare ionised gas and stellar kinematics of 16 star-forming galaxies ($\log(M_\star/M_\odot)=9.7-11.2$, SFR=6-86 $M_\odot/yr$) at $z\sim1$ using near-infrared integral field spectroscopy (IFS) of H$α$ emission from the KMOS$^{\rm 3D}$ survey and optical slit spectroscopy of stellar absorption and gas emission from the LEGA-C survey. H$α$ is dynamically colder than stars, with higher disc rota… ▽ More We compare ionised gas and stellar kinematics of 16 star-forming galaxies ($\log(M_\star/M_\odot)=9.7-11.2$, SFR=6-86 $M_\odot/yr$) at $z\sim1$ using near-infrared integral field spectroscopy (IFS) of H$α$ emission from the KMOS$^{\rm 3D}$ survey and optical slit spectroscopy of stellar absorption and gas emission from the LEGA-C survey. H$α$ is dynamically colder than stars, with higher disc rotation velocities (by ~45 per cent) and lower disc velocity dispersions (by a factor ~2). This is similar to trends observed in the local Universe. We find higher rotational support for H$α$ relative to [OII], potentially explaining systematic offsets in kinematic scaling relations found in the literature. Regarding dynamical mass measurements, for six galaxies with cumulative mass profiles from Jeans Anisotropic Multi-Gaussian Expansion (JAM) models the H$α$ dynamical mass models agree remarkably well out to ~10 kpc for all but one galaxy (average $ΔM_{\rm dyn}(R_{e,\rm F814W})<0.1$ dex). Simpler dynamical mass estimates based on integrated stellar velocity dispersion are less accurate (standard deviation 0.24 dex). Differences in dynamical mass estimates are larger, for example, for galaxies with stronger misalignments of the H$α$ kinematic major axis and the photometric position angle, highlighting the added value of IFS observations for dynamics studies. The good agreement between the JAM models and the dynamical models based on H$α$ kinematics at $z\sim1$ corroborates the validity of dynamical mass measurements from H$α$ IFS observations also for higher redshift rotating disc galaxies. △ Less

Submitted 12 December, 2023; v1 submitted 6 October, 2022; originally announced October 2022.

Comments: 29 pages, 16 figures, 4 tables; accepted for publication in MNRAS

The Velocity Dispersion Function for Massive Quiescent and Star-Forming Galaxies at 0.6 $<$ z $\leq$ 1.0

Authors: Lance Taylor, Rachel Bezanson, Arjen van der Wel, Alan Pearl, Eric F. Bell, Francesco D'Eugenio, Marijn Franx, Michael V. Maseda, Adam Muzzin, David Sobral, Caroline Straatman, Katherine E. Whitaker, Po-Feng Wu

Abstract: We present the first direct spectroscopic measurement of the stellar velocity dispersion function (VDF) for massive quiescent and star-forming galaxies at $0.6 < z \leq 1.0$. For this analysis we use individual measurements of stellar velocity dispersion from high-S/N spectra from the public Large Early Galaxy Astrophysics Census (LEGA-C) survey. We report a remarkable stability of the VDF for bot… ▽ More We present the first direct spectroscopic measurement of the stellar velocity dispersion function (VDF) for massive quiescent and star-forming galaxies at $0.6 < z \leq 1.0$. For this analysis we use individual measurements of stellar velocity dispersion from high-S/N spectra from the public Large Early Galaxy Astrophysics Census (LEGA-C) survey. We report a remarkable stability of the VDF for both quiescent and star-forming galaxies within this redshift range, though we note the presence of weak evolution in the number densities of star-forming galaxies. We compare both VDFs with previous direct and inferred measurements at local and intermediate redshifts, with the caveat that previous measurements of the VDF for star-forming galaxies are poorly constrained at all epochs. We emphasize that this work is the first to directly push to low-stellar velocity dispersion ($σ_\star > 100$ km s$^{-1}$) and extend to star-forming galaxies. We are largely consistent with the high-sigma tail measured from BOSS, and we find that the VDF remains constant from the median redshift of LEGA-C, $z\sim0.8$, to the present day. △ Less

Submitted 30 September, 2022; originally announced October 2022.

Comments: 11 pages, 3 figures, accepted for publication in ApJ

The Mass Scale of High-Redshift Galaxies: Virial Mass Estimates Calibrated with Stellar Dynamical Models from LEGA-C

Authors: Arjen van der Wel, Josha van Houdt, Rachel Bezanson, Marijn Franx, Francesco D'Eugenio, Caroline Straatman, Eric F. Bell, Adam Muzzin, David Sobral, Michael V. Maseda, Anna de Graaff, Bradford P. Holden

Abstract: Dynamical models for $673$ galaxies at $z=0.6-1.0$ with spatially resolved (long-slit) stellar kinematic data from LEGA-C are used to calibrate virial mass estimates defined as $M_{\rm{vir}}=K σ'^2_{\star,\rm{int}} R$, with $K$ a scaling factor, $σ'_{\star,\rm{int}}$ the spatially-integrated stellar velocity second moment from the LEGA-C survey and $R$ the effective radius measured from a Sérsic p… ▽ More Dynamical models for $673$ galaxies at $z=0.6-1.0$ with spatially resolved (long-slit) stellar kinematic data from LEGA-C are used to calibrate virial mass estimates defined as $M_{\rm{vir}}=K σ'^2_{\star,\rm{int}} R$, with $K$ a scaling factor, $σ'_{\star,\rm{int}}$ the spatially-integrated stellar velocity second moment from the LEGA-C survey and $R$ the effective radius measured from a Sérsic profile fit to HST imaging. The sample is representative for $M_{\star}>3\times10^{10}~M_{\odot}$ and includes all types of galaxies, irrespective of morphology and color. We demonstrate that using $R=R_{\rm{sma}}$~(the semi-major axis length of the ellipse that encloses 50\% of the light) in combination with an inclination correction on $σ'_{\star,\rm{int}}$~produces an unbiased $M_{\rm{vir}}$. We confirm the importance of projection effects on $σ'_{\star,\rm{int}}$ by showing the existence of a similar residual trend between virial mass estimates and inclination for the nearby early-type galaxies in the ATLAS$^{\rm{3D}}$~survey. Also, as previously shown, when using a Sérsic profile-based $R$ estimate, then a Sérsic index-dependent correction to account for non-homology in the radial profiles is required. With respect to analogous dynamical models for low-redshift galaxies from the ATLAS$^{\rm{3D}}$~survey we find a systematic offset of 0.1 dex in the calibrated virial constant for LEGA-C, which may be due to physical differences between the galaxy samples or an unknown systematic error. Either way, with our work we establish a common mass scale for galaxies across 8 Gyr of cosmic time with a systematic uncertainty of at most 0.1 dex. △ Less

Submitted 26 August, 2022; originally announced August 2022.

Comments: Published in ApJ

Journal ref: The Astrophysical Journal, 2022, 936, 1

Beyond UVJ: Color Selection of Galaxies in the JWST Era

Authors: Jacqueline Antwi-Danso, Casey Papovich, Joel Leja, Danilo Marchesini, Z. Cemile Marsan, Nicholas S. Martis, Ivo Labbé, Adam Muzzin, Karl Glazebrook, Caroline M. S. Straatman, Kim-Vy H. Tran

Abstract: We present a new rest-frame color-color selection method using "synthetic $u_s-g_s$ and $g_s-i_s$'', $(ugi)_s$ colors to identify star-forming and quiescent galaxies. Our method is similar to the widely-used $U-V$ versus $V-J$ ($UVJ$) diagram. However, $UVJ$ suffers known systematics. Spectroscopic campaigns have shown that $UVJ$-selected quiescent samples at $z \gtrsim 3$ include $\sim 10-30\%$ c… ▽ More We present a new rest-frame color-color selection method using "synthetic $u_s-g_s$ and $g_s-i_s$'', $(ugi)_s$ colors to identify star-forming and quiescent galaxies. Our method is similar to the widely-used $U-V$ versus $V-J$ ($UVJ$) diagram. However, $UVJ$ suffers known systematics. Spectroscopic campaigns have shown that $UVJ$-selected quiescent samples at $z \gtrsim 3$ include $\sim 10-30\%$ contamination from galaxies with dust-obscured star formation and strong emission lines. Moreover, at $z>3$, $UVJ$ colors are extrapolated because the rest-frame J-band shifts beyond the coverage of the deepest bandpasses at $< 5~μm$ (typically $Spitzer$/IRAC 4.5 $μm$ or future $JWST$/NIRCam observations). We demonstrate that $(ugi)_s$ offers improvements to $UVJ$ at $z>3$, and can be applied to galaxies in the $JWST$ era. We apply $(ugi)_s$ selection to galaxies at $0.5<z<6$ from the (observed) 3D-HST and UltraVISTA catalogs, and to the (simulated) JAGUAR catalogs. We show that extrapolation can affect $(V-J)_0$ color by up to 1 magnitude, but changes $(u_s-g_s)_0$ color by $\leq$ 0.2 mag, even at $z\simeq 6$. While $(ugi)_s$-selected quiescent samples are comparable to $UVJ$ in completeness (both achieve $\sim$85-90% at $z=3-3.5$), $(ugi)_s$ reduces contamination in quiescent samples by nearly a factor of two, from $\simeq$35% to $\simeq$17% at $z=3$, and from $\simeq $60% to $\simeq $33% at $z=6$. This leads to improvements in the true-to-false-positive ratio (TP/FP), where we find TP/FP $\gtrsim$ 2.2 for $(ugi)_s$ at $z \simeq 3.5 - 6$, compared to TP/FP $<$ 1 for $UVJ$-selected samples. This indicates that contaminants will outnumber true quiescent galaxies in $UVJ$ at these redshifts, while $(ugi)_s$ will provide higher-fidelity samples. △ Less

Submitted 16 November, 2022; v1 submitted 14 July, 2022; originally announced July 2022.

Comments: 33 pages, 20 figures, accepted for publication in ApJ

Low star-formation activity and low gas content of quiescent galaxies at $z=$ 3.5-4.0 constrained with ALMA

Authors: Tomoko L. Suzuki, Karl Glazebrook, Corentin Schreiber, Tadayuki Kodama, Glenn G. Kacprzak, Roger Leiton, Themiya Nanayakkara, Pascal A. Oesch, Casey Papovich, Lee Spitler, Caroline M. S. Straatman, Kim-Vy Tran, Tao Wang

Abstract: The discovery in deep near-infrared surveys of a population of massive quiescent galaxies at $z>3$ has given rise to the question of how they came to be quenched so early in the history of the Universe. Measuring their molecular gas properties can distinguish between physical processes where they stop forming stars due to a lack of fuel versus those where star-formation efficiency is reduced and t… ▽ More The discovery in deep near-infrared surveys of a population of massive quiescent galaxies at $z>3$ has given rise to the question of how they came to be quenched so early in the history of the Universe. Measuring their molecular gas properties can distinguish between physical processes where they stop forming stars due to a lack of fuel versus those where star-formation efficiency is reduced and the gas is retained. We conducted Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of four quiescent galaxies at $z=$ 3.5-4.0 found by the Fourstar Galaxy Evolution Survey (ZFOURGE) and a serendipitous optically dark galaxy at $z=3.71$. We aim to investigate the presence of dust-obscured star-formation and their gas content by observing the dust continuum emission at Band-7 and the atomic carbon [C I]($^3P_1$-$^3P_0$) line at 492.16 GHz. Among the four quiescent galaxies, only one source is detected in the dust continuum at $λ_{\rm obs} = 870 {\rm μm}$. The sub-mm observations confirm their passive nature, and all of them are located more than four times below the main sequence of star-forming galaxies at $z=3.7$. None of the targets are detected in [C I], constraining their gas mass fractions to be $<$ 20%. These gas mass fractions are more than three times lower than the scaling relation for star-forming galaxies at $z=3.7$. These results support scenarios where massive galaxies at $z=$ 3.5-4.0 quench by consuming/expelling all the gas rather than by reducing the efficiency of the conversion of their gas into stars. △ Less

Submitted 28 June, 2022; originally announced June 2022.

Comments: 13 pages, 5 figures, 1 table, Accepted for publication in ApJ

LEGA-C: analysis of dynamical masses from ionized gas and stellar kinematics at z~0.8

Authors: Caroline M. S. Straatman, Arjen van der Wel, Josha van Houdt, Rachel Bezanson, Eric F. Bell, Pieter van Dokkum, Francesco D'Eugenio, Marijn Franx, Anna Gallazzi, Anna de Graaff, Michael Maseda, Sharon E. Meidt, Adam Muzzin, David Sobral, Po-Feng Wu

Abstract: We compare dynamical mass estimates based on spatially extended stellar and ionized gas kinematics ($\mathrm{M_{dyn,*}}$ and $\mathrm{M_{dyn,eml}}$, respectively) of 157 star forming galaxies at $0.6\leq z<1$. Compared to $z\sim0$, these galaxies have enhanced star formation rates, with stellar feedback likely affecting the dynamics of the gas. We use LEGA-C DR3, the highest redshift dataset provi… ▽ More We compare dynamical mass estimates based on spatially extended stellar and ionized gas kinematics ($\mathrm{M_{dyn,*}}$ and $\mathrm{M_{dyn,eml}}$, respectively) of 157 star forming galaxies at $0.6\leq z<1$. Compared to $z\sim0$, these galaxies have enhanced star formation rates, with stellar feedback likely affecting the dynamics of the gas. We use LEGA-C DR3, the highest redshift dataset providing sufficiently deep measurements of a $K_s-$band limited sample. For $\mathrm{M_{dyn,*}}$ we use Jeans Anisotropic Multi-Gaussian Expansion models. For $\mathrm{M_{dyn,eml}}$ we first fit a custom model of a rotating exponential disk with uniform dispersion, whose light is projected through a slit and corrected for beam smearing. We then apply an asymmetric drift correction based on assumptions common in the literature to the fitted kinematic components to obtain the circular velocity, assuming hydrostatic equilibrium. Within the half-light radius, $\mathrm{M_{dyn,eml}}$ is on average lower than $\mathrm{M_{dyn,*}}$, with a mean offset of $-0.15\pm0.016$ dex and galaxy-to-galaxy scatter of $0.19$ dex, reflecting the combined random uncertainty. While data of higher spatial resolution are needed to understand this small offset, it supports the assumption that the galaxy-wide ionized gas kinematics do not predominantly originate from disruptive events such as star formation driven outflows. However, a similar agreement can be obtained without modeling from the integrated emission line dispersions for axis ratios $q<0.8$. This suggests that our current understanding of gas kinematics is not sufficient to efficiently apply asymmetric drift corrections to improve dynamical mass estimates compared to observations lacking the $S/N$ required for spatially extended dynamics. △ Less

Submitted 11 March, 2022; originally announced March 2022.

Comments: 19 pages, 7 figures, accepted for publication by ApJ, in press

The LEGA-C of nature and nurture in stellar populations of galaxies at z~0.6-1.0: D4000 and H-delta reveal different assembly histories for quiescent galaxies in different environments

Authors: David Sobral, Arjen van der Wel, Rachel Bezanson, Eric Bell, Adam Muzzin, Francesco D'Eugenio, Behnam Darvish, Anna Gallazzi, Po-Feng Wu, Michael Maseda, Jorryt Matthee, Ana Paulino-Afonso, Caroline Straatman, Pieter van Dokkum

Abstract: Galaxy evolution is driven by a variety of physical processes which are predicted to proceed at different rates for different dark matter haloes and environments across cosmic times. A record of this evolution is preserved in galaxy stellar populations, which we can access using absorption-line spectroscopy. Here we explore the large LEGA-C survey (DR3) to investigate the role of the environment a… ▽ More Galaxy evolution is driven by a variety of physical processes which are predicted to proceed at different rates for different dark matter haloes and environments across cosmic times. A record of this evolution is preserved in galaxy stellar populations, which we can access using absorption-line spectroscopy. Here we explore the large LEGA-C survey (DR3) to investigate the role of the environment and stellar mass on stellar populations at z~0.6-1.0 in the COSMOS field. Leveraging the statistical power and depth of LEGA-C, we reveal significant gradients in D4000 and H-delta equivalent widths (EWs) distributions over the stellar mass vs environment 2D spaces for the massive galaxy population (M>10^10 M$_{\odot}$) at z~0.6-1.0. D4000 and H-delta EWs primarily depend on stellar mass, but they also depend on environment at fixed stellar mass. By splitting the sample into centrals and satellites, and in terms of star-forming galaxies and quiescent galaxies, we reveal that the significant environmental trends of D4000 and H-delta EW when controlling for stellar mass are driven by quiescent galaxies. Regardless of being centrals or satellites, star-forming galaxies reveal D4000 and H-delta EWs which depend strongly on their stellar mass and are completely independent of the environment at 0.6<z<1.0. The environmental trends seen for satellite galaxies are fully driven by the trends that hold only for quiescent galaxies, combined with the strong environmental dependency of the quiescent fraction at fixed stellar mass. Our results are consistent with recent predictions from simulations that point towards massive galaxies forming first in over-densities or the most compact dark matter haloes. △ Less

Submitted 15 December, 2021; originally announced December 2021.

Comments: Accepted for publication in ApJ. 18 pages, 9 figures

Ubiquitous [OII] Emission in Quiescent Galaxies at z ~ 0.85

Authors: Michael V. Maseda, Arjen van der Wel, Marijn Franx, Eric F. Bell, Rachel Bezanson, Adam Muzzin, David Sobral, Francesco D'Eugenio, Anna Gallazzi, Anna de Graaff, Joel Leja, Caroline Straatman, Katherine E. Whitaker, Christina C. Williams, Po-Feng Wu

Abstract: Using deep rest-frame optical spectroscopy from the Large Early Galaxy Astrophysical Census (LEGA-C) survey, conducted using VIMOS on the ESO Very Large Telescope, we systematically search for low-ionization [OII] 3726,3729 emission in the spectra of a mass-complete sample of z~0.85 galaxies. Intriguingly, we find that 59 percent of UVJ-quiescent (i.e. non star-forming) galaxies in the sample have… ▽ More Using deep rest-frame optical spectroscopy from the Large Early Galaxy Astrophysical Census (LEGA-C) survey, conducted using VIMOS on the ESO Very Large Telescope, we systematically search for low-ionization [OII] 3726,3729 emission in the spectra of a mass-complete sample of z~0.85 galaxies. Intriguingly, we find that 59 percent of UVJ-quiescent (i.e. non star-forming) galaxies in the sample have ionized gas, as traced by [OII] emission, detected above our completeness limit of 1.5 Angstroms. The median stacked spectrum of the lowest equivalent width quiescent galaxies also shows [OII] emission. The overall fraction of sources with [OII] above our equivalent width limit is comparable to what we find in the low-redshift Universe from GAMA and MASSIVE, except perhaps at the highest stellar masses (log Mstar/Msol > 11.5). However, stacked spectra for the individual low-equivalent width systems uniquely indicates ubiquitous [OII] emission in the higher-z LEGA-C sample, with typical [OII] luminosities per unit stellar mass that are a factor of 3 larger than the lower-z GAMA sample. Star formation in these otherwise quiescent galaxies could play a role in producing the [OII] emission at higher-z, although it is unlikely to provide the bulk of the ionizing photons. More work is required to fully quantify the contributions of evolved stellar populations or active galactic nuclei to the observed spectra. △ Less

Submitted 30 September, 2021; originally announced October 2021.

Comments: ApJ accepted: 16 pages, 11 figures

Introducing the FLAMINGOS-2 Split-K Medium Band Filters: The Impact on Photometric Selection of High-z Galaxies in the FENIKS-pilot survey

Authors: James Esdaile, Ivo Labbe, Karl Glazebrook, Jacqueline Antwi-Danso, Casey Papovich, Edward Taylor, Z. Cemile Marsan, Adam Muzzin, Caroline M. S. Straatman, Danilo Marchesini, Ruben Diaz, Lee Spitler, Kim-Vy H. Tran, Stephen Goodsell

Abstract: Deep near-infrared photometric surveys are efficient in identifying high-redshift galaxies, however they can be prone to systematic errors in photometric redshift. This is particularly salient when there is limited sampling of key spectral features of a galaxy's spectral energy distribution (SED), such as for quiescent galaxies where the expected age-sensitive Balmer/4000 A break enter the $K$-ban… ▽ More Deep near-infrared photometric surveys are efficient in identifying high-redshift galaxies, however they can be prone to systematic errors in photometric redshift. This is particularly salient when there is limited sampling of key spectral features of a galaxy's spectral energy distribution (SED), such as for quiescent galaxies where the expected age-sensitive Balmer/4000 A break enter the $K$-band at $z>4$. With single filter sampling of this spectral feature, degeneracies between SED models and redshift emerge. A potential solution to this comes from splitting the $K$-band into multiple filters. We use simulations to show an optimal solution is to add two medium-band filters, $K_\mathrm{blue}$ ($λ_\mathrm{cen}$=2.06 $μ$m, $Δλ$=0.25 $μ$m) and $K_\mathrm{red}$ ($λ_\mathrm{cen}$=2.31 $μ$m, $Δλ$=0.27 $μ$m), that are complementary to the existing $K_\mathrm{s}$ filter. We test the impact of the $K$-band filters with simulated catalogues comprised of galaxies with varying ages and signal-to-noise. The results suggest that the $K$-band filters do improve photometric redshift constraints on $z>4$ quiescent galaxies, increasing precision and reducing outliers by up to 90$\%$. We find that the impact from the $K$-band filters depends on the signal-to-noise, the redshift and the SED of the galaxy. The filters we designed were built and used to conduct a pilot of the FLAMINGOS-2 Extra-galactic Near-Infrared $K$-band Split (FENIKS) survey. While no new $z>4$ quiescent galaxies are identified in the limited area pilot, the $K_\mathrm{blue}$ and $K_\mathrm{red}$ filters indicate strong Balmer/4000 A breaks in existing candidates. Additionally we identify galaxies with strong nebular emission lines, for which the $K$-band filters increase photometric redshift precision and in some cases indicate extreme star-formation. △ Less

Submitted 14 September, 2021; originally announced September 2021.

Comments: 27 pages, 18 figures. Accepted for publication in the Astronomical Journal

Towards precise galaxy evolution: a comparison between spectral indices of $z\sim1$ galaxies in the IllustrisTNG simulation and the LEGA-C survey

Authors: Po-Feng Wu, Dylan Nelson, Arjen van der Wel, Annalisa Pillepich, Stefano Zibetti, Rachel Bezanson, Francesco D'Eugenio, Anna Gallazzi, Camilla Pacifici, Caroline M. S. Straatman, Ivana Barišić, Eric F. Bell, Michael V. Maseda, Adam Muzzin, David Sobral, Katherine E. Whitaker

Abstract: We present the first comparison of observed stellar continuum spectra of high-redshift galaxies and mock galaxy spectra generated from hydrodynamical simulations. The mock spectra are produced from the IllustrisTNG TNG100 simulation combined with stellar population models and take into account dust attenuation and realistic observational effects (aperture effects and noise). We compare the simulat… ▽ More We present the first comparison of observed stellar continuum spectra of high-redshift galaxies and mock galaxy spectra generated from hydrodynamical simulations. The mock spectra are produced from the IllustrisTNG TNG100 simulation combined with stellar population models and take into account dust attenuation and realistic observational effects (aperture effects and noise). We compare the simulated $D_n4000$ and EW(H$δ$) of galaxies with $10.5 \leq \log(M_\ast/M_\odot) \leq 11.5$ at $0.6 \leq z \leq 1.0$ to the observed distributions from the LEGA-C survey. TNG100 globally reproduces the observed distributions of spectral indices, implying that the age distribution of galaxies in TNG100 is generally realistic. Yet there are small but significant differences. For old galaxies, TNG100 shows small $D_n4000$ when compared to LEGA-C, while LEGA-C galaxies have larger EW(H$δ$) at fixed $D_n4000$. There are several possible explanations: 1) LEGA-C galaxies have overall older ages combined with small contributions (a few percent in mass) from younger ($<1$~Gyr) stars, while TNG100 galaxies may not have such young sub-populations; 2) the spectral mismatch could be due to systematic uncertainties in the stellar population models used to convert stellar ages and metallicities to observables. In conclusion, the latest cosmological galaxy formation simulations broadly reproduce the global age distribution of galaxies at $z\sim1$ and, at the same time, the high quality of the latest observed and simulated datasets help constrain stellar population synthesis models as well as the physical models underlying the simulations. △ Less

Submitted 26 September, 2021; v1 submitted 23 August, 2021; originally announced August 2021.

Comments: 17 pages, 12 figures, accepted by AJ

Stellar Dynamical Models for 797 $z\sim 0.8$ Galaxies from LEGA-C

Authors: Josha van Houdt, Arjen van der Wel, Rachel Bezanson, Marijn Franx, Francesco D'Eugenio, Ivana Barisic, Eric F. Bell, Anna Gallazzi, Anna de Graaff, Michael V. Maseda, Camilla Pacifici, Jesse van de Sande, David Sobral, Caroline Straatman, Po-Feng Wu

Abstract: We present spatially resolved stellar kinematics for 797 $z=0.6-1$ galaxies selected from the LEGA-C survey and construct axisymmetric Jeans models to quantify their dynamical mass and degree of rotational support. The survey is $K_s$-band selected, irrespective of color or morphological type, and allows for a first assessment of the stellar dynamical structure of the general $L^*$ galaxy populati… ▽ More We present spatially resolved stellar kinematics for 797 $z=0.6-1$ galaxies selected from the LEGA-C survey and construct axisymmetric Jeans models to quantify their dynamical mass and degree of rotational support. The survey is $K_s$-band selected, irrespective of color or morphological type, and allows for a first assessment of the stellar dynamical structure of the general $L^*$ galaxy population at large lookback time. Using light profiles from Hubble Space Telescope imaging as a tracer, our approach corrects for observational effects (seeing convolution and slit geometry), and uses well-informed priors on inclination, anisotropy and a non-luminous mass component. Tabulated data include total mass estimates in a series of spherical apertures (1, 5, and 10 kpc; 1$\times$ and 2$\times$\re), as well as rotational velocities, velocity dispersions and anisotropy. We show that almost all star-forming galaxies and $\sim$50\% of quiescent galaxies are rotation-dominated, with deprojected $V/σ\sim1-2$. Revealing the complexity in galaxy evolution, we find that the most massive star-forming galaxies are among the most rotation-dominated, and the most massive quiescent galaxies among the least rotation-dominated galaxies. These measurements set a new benchmark for studying galaxy evolution, using stellar dynamical structure for galaxies at large lookback time. Together with the additional information on stellar population properties from the LEGA-C spectra, the dynamical mass and $V/σ$ measurements presented here create new avenues for studying galaxy evolution at large lookback time. △ Less

Submitted 18 August, 2021; originally announced August 2021.

Comments: Accepted for publication in ApJ. Data table will be published with journal article and is now available upon request

The Large Early Galaxy Astrophysics Census (LEGA-C) Data Release 3: 3000 High-Quality Spectra of $K_s$-selected galaxies at $z>0.6$

Authors: Arjen van der Wel, Rachel Bezanson, Francesco D'Eugenio, Caroline Straatman, Marijn Franx, Josha van Houdt, Michael V. Maseda, Anna Gallazzi, Po-Feng Wu, Camilla Pacifici, Ivana Barisic, Gabriel B. Brammer, Juan Carlos Munoz-Mateos, Sarah Vervalcke, Stefano Zibetti, David Sobral, Anna de Graaff, Joao Calhau, Yasha Kaushal, Adam Muzzin, Eric F. Bell, Pieter G. van Dokkum

Abstract: We present the third and final data release of the Large Early Galaxy Astrophysics Census (LEGA-C), an ESO/VLT public spectroscopic survey targeting $0.6 < z < 1.0$, Ks-selected galaxies. The data release contains 3528 spectra with measured stellar velocity dispersions and stellar population properties, a 25-fold increase in sample size compared to previous work. This $K_s$-selected sample probes… ▽ More We present the third and final data release of the Large Early Galaxy Astrophysics Census (LEGA-C), an ESO/VLT public spectroscopic survey targeting $0.6 < z < 1.0$, Ks-selected galaxies. The data release contains 3528 spectra with measured stellar velocity dispersions and stellar population properties, a 25-fold increase in sample size compared to previous work. This $K_s$-selected sample probes the galaxy population down to $\sim0.3 L^*$, for all colors and morphological types. Along with the spectra we publish a value-added catalog with stellar and ionized gas velocity dispersions, stellar absorption line indices, emission line fluxes and equivalent widths, complemented with structural parameters measured from HST/ACS imaging. With its combination of high precision and large sample size, LEGA-C provides a new benchmark for galaxy evolution studies. △ Less

Submitted 2 August, 2021; originally announced August 2021.

Comments: Accepted for publication in ApJS. Data available at http://archive.eso.org/cms/eso-archive-news/Third-and-final-release-of-the-Large-Early-Galaxy-Census-LEGA-C-Spectroscopic-Public-Survey-published.html or https://users.ugent.be/\string~avdrwel/research.html\#legac

The Fundamental Plane in the LEGA-C Survey: unraveling the $M/L$ variations of massive star-forming and quiescent galaxies at $z\sim0.8$

Authors: Anna de Graaff, Rachel Bezanson, Marijn Franx, Arjen van der Wel, Bradford Holden, Jesse van de Sande, Eric F. Bell, Francesco D'Eugenio, Michael V. Maseda, Adam Muzzin, David Sobral, Caroline M. S. Straatman, Po-Feng Wu

Abstract: We explore the connection between the kinematics, structures and stellar populations of massive galaxies at $0.6<z<1.0$ using the Fundamental Plane (FP). Combining stellar kinematic data from the Large Early Galaxy Astrophysics Census (LEGA-C) survey with structural parameters measured from deep Hubble Space Telescope imaging, we obtain a sample of 1419 massive ($\log(M_*/M_\odot) >10.5$) galaxies… ▽ More We explore the connection between the kinematics, structures and stellar populations of massive galaxies at $0.6<z<1.0$ using the Fundamental Plane (FP). Combining stellar kinematic data from the Large Early Galaxy Astrophysics Census (LEGA-C) survey with structural parameters measured from deep Hubble Space Telescope imaging, we obtain a sample of 1419 massive ($\log(M_*/M_\odot) >10.5$) galaxies that span a wide range in morphology, star formation activity and environment, and therefore is representative of the massive galaxy population at $z\sim0.8$. We find that quiescent and star-forming galaxies occupy the parameter space of the $g$-band FP differently and thus have different distributions in the dynamical mass-to-light ratio ($M_{\rm dyn}/L_g$), largely owing to differences in the stellar age and recent star formation history, and, to a lesser extent, the effects of dust attenuation. In contrast, we show that both star-forming and quiescent galaxies lie on the same mass FP at $z\sim 0.8$, with a comparable level of intrinsic scatter about the plane. We examine the variation in $M_{\rm dyn}/M_*$ through the thickness of the mass FP, finding no significant residual correlations with stellar population properties, Sérsic index, or galaxy overdensity. Our results suggest that, at fixed size and velocity dispersion, the variations in $M_{\rm dyn}/L_g$ of massive galaxies reflect an approximately equal contribution of variations in $M_*/L_g$, and variations in the dark matter fraction or initial mass function. △ Less

Submitted 23 March, 2021; originally announced March 2021.

Comments: 32 pages, 19 figures (including appendices); accepted for publication in ApJ

A low [CII]/[NII] ratio in the center of a massive galaxy at z=3.7: witnessing the transition to quiescence at high-redshift?

Authors: C. Schreiber, K. Glazebrook, C. Papovich, T. Diaz-Santos, A. Verma, D. Elbaz, G. G. Kacprzak, T. Nanayakkara, P. Oesch, M. Pannella, L. Spitler, C. Straatman, K. -V. Tran, T. Wang

Abstract: Understanding the process of quenching is one of the major open questions in galaxy evolution, and crucial insights may be obtained by studying quenched galaxies at high redshifts, at epochs when the Universe and the galaxies were younger and simpler to model. However, establishing the degree of quiescence in high redshift galaxies is a challenging task. One notable example is Hyde, a recently dis… ▽ More Understanding the process of quenching is one of the major open questions in galaxy evolution, and crucial insights may be obtained by studying quenched galaxies at high redshifts, at epochs when the Universe and the galaxies were younger and simpler to model. However, establishing the degree of quiescence in high redshift galaxies is a challenging task. One notable example is Hyde, a recently discovered galaxy at z=3.709. As compact (r~0.5 kpc) and massive (M*~1e11 Msun) as its quenched neighbor Jekyll, it is also extremely obscured yet only moderately luminous in the sub-millimeter. Panchromatic modeling suggested it could be the first galaxy found in transition to quenching at z>3, however the data were also consistent with a broad range of star-formation activity, including moderate SFR in the lower scatter of the galaxy main-sequence (MS). Here, we describe ALMA observations of the [CII] 157um and [NII] 205um far-infrared emission lines. The [CII] emission within the half-light radius is dominated by ionized gas, while the outskirts are dominated by PDRs or neutral gas. This suggests that the ionization in the center is not primarily powered by on-going star formation, and could come instead from remnant stellar populations formed in an older burst, or from a moderate AGN. Accounting for this information in the multi-wavelength modeling provides a tighter constraint on the star formation rate of SFR=$50^{+24}_{-18}$ Msun/yr. This rules out fully quenched solutions, and favors SFRs more than factor of two lower than expected for a galaxy on the MS, confirming the nature of Hyde as a transition galaxy. Theses results suggest that quenching happens from inside-out, and starts before the galaxy expels or consumes all its gas reservoirs. Similar observations of a larger sample would determine whether this is an isolated case or the norm for quenching at high-redshift. [abriged] △ Less

Submitted 27 November, 2020; originally announced November 2020.

Comments: Accepted for publication in A&A. 14 pages, 9 figures

Journal ref: A&A 646, A68 (2021)

Tightly coupled morpho-kinematic evolution for massive star-forming and quiescent galaxies across 7 Gyr of cosmic time

Authors: Anna de Graaff, Rachel Bezanson, Marijn Franx, Arjen van der Wel, Eric F. Bell, Francesco D'Eugenio, Bradford Holden, Michael V. Maseda, Adam Muzzin, Camilla Pacifici, Jesse van de Sande, David Sobral, Caroline M. S. Straatman, Po-Feng Wu

Abstract: We use the Fundamental Plane (FP) to measure the redshift evolution of the dynamical mass-to-light ratio ($M_{\mathrm{dyn}}/L$) and the dynamical-to-stellar mass ratio ($M_{\mathrm{dyn}}/M_*$). Although conventionally used to study the properties of early-type galaxies, we here obtain stellar kinematic measurements from the Large Early Galaxy Astrophysics Census (LEGA-C) Survey for a sample of… ▽ More We use the Fundamental Plane (FP) to measure the redshift evolution of the dynamical mass-to-light ratio ($M_{\mathrm{dyn}}/L$) and the dynamical-to-stellar mass ratio ($M_{\mathrm{dyn}}/M_*$). Although conventionally used to study the properties of early-type galaxies, we here obtain stellar kinematic measurements from the Large Early Galaxy Astrophysics Census (LEGA-C) Survey for a sample of $\sim1400$ massive ($\log( M_*/M_\odot) >10.5$) galaxies at $0.6<z<1.0$ that span a wide range in star formation activity. In line with previous studies, we find a strong evolution in $M_{\mathrm{dyn}}/L_g$ with redshift. In contrast, we find only a weak dependence of the mean value of $M_{\mathrm{dyn}}/M_*$ on the specific star formation rate, and a redshift evolution that likely is explained by systematics. Therefore, we demonstrate that star-forming and quiescent galaxies lie on the same, stable mass FP across $0<z<1$, and that the decrease in $M_{\mathrm{dyn}}/L_g$ toward high redshift can be attributed entirely to evolution of the stellar populations. Moreover, we show that the growth of galaxies in size and mass is constrained to occur within the mass FP. Our results imply either minimal structural evolution in massive galaxies since $z\sim1$, or a tight coupling in the evolution of their morphological and dynamical properties, and establish the mass FP as a tool for studying galaxy evolution with low impact from progenitor bias. △ Less

Submitted 23 October, 2020; originally announced October 2020.

Comments: 7 pages, 4 figures; accepted for publication in ApJL

Consistent dynamical and stellar masses with potential light IMF in massive quiescent galaxies at $3 < z < 4$ using velocity dispersions measurements with MOSFIRE

Authors: James Esdaile, Karl Glazebrook, Ivo Labbe, Edward Taylor, Corentin Schreiber, Themiya Nanayakkara, Glenn G. Kacprzak, Pascal A. Oesch, Kim-Vy H. Tran, Casey Papovich, Lee Spitler, Caroline M. S. Straatman

Abstract: We present the velocity dispersion measurements of four massive $\sim10^{11}M_\odot$ quiescent galaxies at $3.2 < z < 3.7$ based on deep H and K$-$band spectra using the Keck/MOSFIRE near-infrared spectrograph. We find high velocity dispersions of order $σ_e\sim250$ km/s based on strong Balmer absorption lines and combine these with size measurements based on HST/WFC3 F160W imaging to infer dynami… ▽ More We present the velocity dispersion measurements of four massive $\sim10^{11}M_\odot$ quiescent galaxies at $3.2 < z < 3.7$ based on deep H and K$-$band spectra using the Keck/MOSFIRE near-infrared spectrograph. We find high velocity dispersions of order $σ_e\sim250$ km/s based on strong Balmer absorption lines and combine these with size measurements based on HST/WFC3 F160W imaging to infer dynamical masses. The velocity dispersion are broadly consistent with the high stellar masses and small sizes. Together with evidence for quiescent stellar populations, the spectra confirm the existence of a population of massive galaxies that formed rapidly and quenched in the early universe $z>4$. Investigating the evolution at constant velocity dispersion between $z\sim3.5$ and $z\sim2$, we find a large increase in effective radius $0.35\pm0.12$ dex and in dynamical-to-stellar mass ratio $<$log(M$_{\textrm{dyn}}$/M*)$>$ of 0.33$\pm0.08$ dex, with low expected contribution from dark matter. The dynamical masses for our $z\sim3.5$ sample are consistent with the stellar masses for a Chabrier initial mass function (IMF), with the ratio $<$log(M$_{\textrm{dyn}}$/M$^*_{\textrm{Ch}})>$ = -0.13$\pm$0.10 dex suggesting an IMF lighter than Salpeter may be common for massive quiescent galaxies at $z>3$. This is surprising in light of the Salpeter or heavier IMFs found for high velocity dispersion galaxies at $z\sim2$ and cores of present-day ellipticals, which these galaxies are thought to evolve into. Future imaging and spectroscopic observations with resolved kinematics using the upcoming James Webb Space Telescope could rule out potential systematics from rotation, and confirm these results. △ Less

Submitted 30 March, 2021; v1 submitted 19 October, 2020; originally announced October 2020.

Comments: 11 pages, 3 figures. Accepted to ApJ Letters

Dust Attenuation Curves at z $\sim$ 0.8 from LEGA-C: Precise Constraints on the Slope and 2175$Å$ Bump Strength

Authors: Ivana Barisic, Camila Pacifici, Arjen van der Wel, Caroline Straatman, Eric F. Bell, Rachel Bezanson, Gabriel Brammer, Francesco D'Eugenio, Marijn Franx, Josha van Houdt, Michael V. Maseda, Adam Muzzin, David Sobral, Po-Feng Wu

Abstract: We present a novel approach to measure the attenuation curves of 485 individual star-forming galaxies with M$_*$ $>$ 10$^{10}$ M$_{\odot}$ based on deep optical spectra from the VLT/VIMOS LEGA-C survey and multi-band photometry in the COSMOS field. Most importantly, we find that the attenuation curves in the rest-frame $3000-4500$A range are typically almost twice as steep as the Milky Way, LMC, S… ▽ More We present a novel approach to measure the attenuation curves of 485 individual star-forming galaxies with M$_*$ $>$ 10$^{10}$ M$_{\odot}$ based on deep optical spectra from the VLT/VIMOS LEGA-C survey and multi-band photometry in the COSMOS field. Most importantly, we find that the attenuation curves in the rest-frame $3000-4500$A range are typically almost twice as steep as the Milky Way, LMC, SMC, and Calzetti attenuation curves, which is in agreement with recent studies of the integrated light of present-day galaxies. The attenuation at $4500$A and the slope strongly correlate with the galaxy inclination: face-on galaxies show less attenuation and steeper curves compared to edge-on galaxies, suggesting that geometric effects dominate observed variations in attenuation. Our new method produces $2175$A UV bump detections for 260 individual galaxies. Even though obvious correlations between UV bump strength and global galaxy properties are absent, strong UV bumps are most often seen in face-on, lower-mass galaxies (10 $<$ log$_{10}$(M$_*$/M$_{\odot}$) $<$ 10.5) with low overall attenuation. Finally, we produce a typical attenuation curve for star-forming galaxies at $z\sim0.8$; this prescription represents the effect of dust on the integrated spectral energy distributions of high-redshift galaxies more accurately than commonly used attenuation laws. △ Less

Submitted 2 October, 2020; originally announced October 2020.

Comments: 15 pages, 10 figures, 1 table. Accepted for publication in The Astrophysical Journal

Inverse stellar population age gradients of post-starburst galaxies at z=0.8 with LEGA-C

Authors: Francesco D'Eugenio, Arjen van der Wel, Po-Feng Wu, Tania M. Barone, Josha van Houdt, Rachel Bezanson, Caroline M. S. Straatman, Camilla Pacifici, Adam Muzzin, Anna Gallazzi, Vivienne Wild, David Sobral, Eric F. Bell, Stefano Zibetti, Lamiya Mowla, Marijn Franx

Abstract: We use deep, spatially resolved spectroscopy from the LEGA-C Survey to study radial variations in the stellar population of 17 spectroscopically-selected post-starburst (PSB) galaxies. We use spectral fitting to measure two Lick indices, $Hδ_A$ and $Fe4383$, and find that, on average, PSB galaxies have radially decreasing $Hδ_A$ and increasing $Fe4383$ profiles. In contrast, a control sample of qu… ▽ More We use deep, spatially resolved spectroscopy from the LEGA-C Survey to study radial variations in the stellar population of 17 spectroscopically-selected post-starburst (PSB) galaxies. We use spectral fitting to measure two Lick indices, $Hδ_A$ and $Fe4383$, and find that, on average, PSB galaxies have radially decreasing $Hδ_A$ and increasing $Fe4383$ profiles. In contrast, a control sample of quiescent, non-PSB galaxies in the same mass range shows outwardly increasing $Hδ_A$ and decreasing $Fe4383$. The observed gradients are weak ($\approx-0.2$ Å/$R_e$), mainly due to seeing convolution. A two-SSP model suggests intrinsic gradients are as strong as observed in local PSB galaxies ($\approx -0.8$ Å$/R_e$). We interpret these results in terms of inside-out growth (for the bulk of the quiescent population) vs star formation occurring last in the centre (for PSB galaxies). At $z\approx0.8$, central starbursts are often the result of gas-rich mergers, as evidenced by the high fraction of PSB galaxies with disturbed morphologies and tidal features (40%). Our results provide additional evidence for multiple paths to quiescence: a standard path, associated with inside-out disc formation and with gradually decreasing star-formation activity, without fundamental structural transformation, and a fast path, associated with centrally-concentrated starbursts, leaving an inverse age gradient and smaller half-light radius. △ Less

Submitted 1 July, 2020; originally announced July 2020.

Comments: 18 pages, 9 figures

A giant galaxy in the young Universe with a massive ring

Authors: Tiantian Yuan, Ahmed Elagali, Ivo Labbe, Glenn G. Kacprzak, Claudia del P. Lagos, Leo Y. Alcorn, Jonathan H. Cohn, Kim-Vy H. Tran, Karl Glazebrook, Brent A. Groves, Kenneth C. Freeman, Lee R. Spitler, Caroline M. S. Straatman, Deanne B. Fisher, Sarah M. Sweet

Abstract: In the local (redshift z~0) Universe, collisional ring galaxies make up only ~0.01% of galaxies and are formed by head-on galactic collisions that trigger radially propagating density waves. These striking systems provide key snapshots for dissecting galactic disks and are studied extensively in the local Universe. However, not much is known about distant (z>0.1) collisional rings. Here we present… ▽ More In the local (redshift z~0) Universe, collisional ring galaxies make up only ~0.01% of galaxies and are formed by head-on galactic collisions that trigger radially propagating density waves. These striking systems provide key snapshots for dissecting galactic disks and are studied extensively in the local Universe. However, not much is known about distant (z>0.1) collisional rings. Here we present a detailed study of a ring galaxy at a look-back time of 10.8 Gyr (z=2.19). Compared with our Milky Way, this galaxy has a similar stellar mass, but has a stellar half-light radius that is 1.5-2.2 times larger and is forming stars 50 times faster. The large, diffuse stellar light outside the star-forming ring, combined with a radial velocity on the ring and an intruder galaxy nearby, provides evidence for this galaxy hosting a collisional ring. If the ring is secularly evolved, the implied large bar in a giant disk would be inconsistent with the current understanding of the earliest formation of barred spirals. Contrary to previous predictions, this work suggests that massive collisional rings were as rare 11 Gyr ago as they are today. Our discovery offers a unique pathway for studying density waves in young galaxies, as well as constraining the cosmic evolution of spiral disks and galaxy groups. △ Less

Submitted 24 May, 2020; originally announced May 2020.

Comments: Author's version for the main article (10 pages). The Supplementary Information (22 pages) and a combined pdf are provided here http://astronomy.swin.edu.au/~tyuan/paper Published version available online http://dx.doi.org/10.1038/s41550-020-1102-7

Journal ref: Nature Astronomy 2020

MOSEL: Strong [OIII]5007 Å Emitting Galaxies at (3<z<4) from the ZFOURGE Survey

Authors: Kim-Vy H. Tran, Ben Forrest, Leo Y. Alcorn, Tiantian Yuan, Themiya Nanayakkara, Jonathan Cohn, Michael Cowley, Karl Glazebrook, Anshu Gupta, Glenn G. Kacprzak, Lisa Kewley, Ivo Labbe, Casey Papovich, Lee Spitler, Caroline M. S. Straatman, Adam Tomczak

Abstract: To understand how strong emission line galaxies (ELGs) contribute to the overall growth of galaxies and star formation history of the universe, we target Strong ELGs (SELGs) from the ZFOURGE imaging survey that have blended (Hb+[OIII]) rest-frame equivalent widths of >230A and 2.5<zphot<4.0. Using Keck/MOSFIRE, we measure 49 redshifts for galaxies brighter than Ks=25 mag as part of our Multi-Objec… ▽ More To understand how strong emission line galaxies (ELGs) contribute to the overall growth of galaxies and star formation history of the universe, we target Strong ELGs (SELGs) from the ZFOURGE imaging survey that have blended (Hb+[OIII]) rest-frame equivalent widths of >230A and 2.5<zphot<4.0. Using Keck/MOSFIRE, we measure 49 redshifts for galaxies brighter than Ks=25 mag as part of our Multi-Object Spectroscopic Emission Line (MOSEL) survey. Our spectroscopic success rate is ~53% and zphot uncertainty is sigma_z= [Delta(z)/(1+z)]=0.0135. We confirm 31 ELGs at 3<zspec<3.8 and show that Strong ELGs have spectroscopic rest-frame [OIII]5007A equivalent widths of 100-500A and tend to be lower mass systems [log(Mstar/Msun)~8.2-9.6] compared to more typical star-forming galaxies. The Strong ELGs lie ~0.9 dex above the star-forming main-sequence at z~3.5 and have high inferred gas fractions of fgas~>60%, i.e. the inferred gas masses can easily fuel a starburst to double stellar masses within ~10-100 Myr. Combined with recent results using ZFOURGE, our analysis indicates that 1) strong [OIII]5007A emission signals an early episode of intense stellar growth in low mass (Mstar<0.1M*) galaxies and 2) many, if not most, galaxies at z>3 go through this starburst phase. If true, low-mass galaxies with strong [OIII]5007A emission (EW_rest>200A) may be an increasingly important source of ionizing UV radiation at z>3. △ Less

Submitted 27 April, 2020; originally announced April 2020.

Comments: Astrophysical Journal, in press

MOSEL Survey: Tracking the Growth of Massive Galaxies at 2<z<4 using Kinematics and the IllustrisTNG Simulation

Authors: Anshu Gupta, Kim-Vy Tran, Jonathan Cohn, Leo Y. Alcorn, Tiantian Yuan, Vicente Rodriguez-Gomez, Anishya Harshan, Ben Forrest, Lisa J. Kewley, Karl Glazebrook, Caroline M. Straatman, Glenn G. Kacprzak, Themiya Nanayakkara, Ivo Labbé, Casey Papovich, Michael Cowley

Abstract: We use K-band spectroscopic data from the Multi-Object Spectroscopic Emission Line (MOSEL) survey to analyze the kinematic properties of galaxies at z>3. Our sample consists of 34 galaxies at 3.0<zspec<3.8 between 9.0<log(M_star)<11.0. We find that galaxies with log(M_star) > 10.2 at z > 3 have 56 +/- 21 km/s lower integrated velocity dispersion compared to galaxies at z ~ 2 of similar stellar mas… ▽ More We use K-band spectroscopic data from the Multi-Object Spectroscopic Emission Line (MOSEL) survey to analyze the kinematic properties of galaxies at z>3. Our sample consists of 34 galaxies at 3.0<zspec<3.8 between 9.0<log(M_star)<11.0. We find that galaxies with log(M_star) > 10.2 at z > 3 have 56 +/- 21 km/s lower integrated velocity dispersion compared to galaxies at z ~ 2 of similar stellar mass. Massive galaxies at z > 3 have either a flat or declining star formation history (SFH), whereas similar stellar mass galaxies at z~2.0 exhibit a slight peak in the past 500 Myrs. Comparing with the IllustrisTNG cosmological simulation, we find that (i) the dynamical mass of massive galaxies in simulations (log(M_star) > 10.0) increases by ~0.1 dex at a fixed stellar mass between z=2.0-3.0, and (ii) dynamical mass growth is coupled with a rapid rise in the ex-situ stellar mass fraction (stars accreted from other galaxies) for massive galaxies at z < 3.5. We speculate that the rising contribution of ex-situ stellar mass to the total stellar mass growth of massive galaxies is driving the higher integrated velocity dispersion and rising SFHs of massive galaxies at z~2.0 compared to galaxies of similar stellar masses at z > 3. △ Less

Submitted 3 March, 2020; originally announced March 2020.

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

Stellar Kinematics and Environment at z~0.8 in the LEGA-C Survey: Massive, Slow-Rotators are Built First in Overdense Environments

Authors: Justin Cole, Rachel Bezanson, Arjen van der Wel, Eric Bell, Francesco D'Eugenio, Marijn Franx, Anna Gallazzi, Josha van Houdt, Adam Muzzin, Camilla Pacifici, Jesse van de Sande, David Sobral, Caroline Straatman, Po-Feng Wu

Abstract: In this letter, we investigate the impact of environment on integrated and spatially-resolved stellar kinematics of a sample of massive, quiescent galaxies at intermediate redshift ($0.6<z<1.0$). For this analysis, we combine photometric and spectroscopic parameters from the UltraVISTA and Large Early Galaxy Astrophysics Census (LEGA-C) surveys in the COSMOS field and environmental measurements. W… ▽ More In this letter, we investigate the impact of environment on integrated and spatially-resolved stellar kinematics of a sample of massive, quiescent galaxies at intermediate redshift ($0.6<z<1.0$). For this analysis, we combine photometric and spectroscopic parameters from the UltraVISTA and Large Early Galaxy Astrophysics Census (LEGA-C) surveys in the COSMOS field and environmental measurements. We analyze the trends with overdensity (1+$δ$) on the rotational support of quiescent galaxies and find no universal trends at either fixed mass or fixed stellar velocity dispersion. This is consistent with previous studies of the local Universe; rotational support of massive galaxies depends primarily on stellar mass. We highlight two populations of massive galaxies ($\log M_\star/M_\odot\geq11$) that deviate from the average mass relation. First, the most massive galaxies in the most under-dense regions ($(1+δ)\leq1$) exhibit elevated rotational support. Similarly, at the highest masses ($\log M_\star/M_\odot\geq11.25$) the range in rotational support is significant in all but the densest regions. This corresponds to an increasing slow-rotator fraction such that only galaxies in the densest environments ($(1+δ)\geq3.5$) are primarily (90$\pm$10\%) slow-rotators.This effect is not seen at fixed velocity dispersion, suggesting minor merging as the driving mechanism: only in the densest regions have the most massive galaxies experienced significant minor merging, building stellar mass and diminishing rotation without significantly affecting the central stellar velocity dispersion. In the local Universe, most massive galaxies are slow-rotators, regardless of environment, suggesting minor merging occurs at later cosmic times $(z\lesssim0.6)$ in all but the most dense environments. △ Less

Submitted 8 January, 2020; originally announced January 2020.

Comments: 10 pages, 5 figures, submitted to ApJ Letters

Reconstructing the observed ionizing photon production efficiency at z~2 using stellar population models

Authors: Themiya Nanayakkara, Jarle Brinchmann, Karl Glazebrook, Rychard Bouwens, Lisa Kewley, Kim-Vy Tran, Michael Cowley, Deanne Fisher, Glenn G. Kacprzak, Ivo Labbe, Caroline Straatman

Abstract: The ionizing photon production efficiency, $ξ_{ion}$, is a critical parameter that provides a number of physical constraints to the nature of the early Universe, including the contribution of galaxies to the timely completion of the reionization of the Universe. Here we use KECK/MOSFIRE and ZFOURGE multi-band photometric data to explore the $ξ_{ion}$ of a population of galaxies at $z\sim2$ with… ▽ More The ionizing photon production efficiency, $ξ_{ion}$, is a critical parameter that provides a number of physical constraints to the nature of the early Universe, including the contribution of galaxies to the timely completion of the reionization of the Universe. Here we use KECK/MOSFIRE and ZFOURGE multi-band photometric data to explore the $ξ_{ion}$ of a population of galaxies at $z\sim2$ with $log_{10}(M_*/M_\odot)\sim9.0-11.5$. Our 130 \Halpha\ detections show a median $log_{10}(ξ_{ion}[Hz/erg])$ of $24.8\pm0.5$ when dust corrected using a Calzetti et al. (2000) dust prescription. Our values are typical of mass/magnitude selected $ξ_{ion}$ values observed in the $z\sim2$ Universe. Using BPASSv2.2.1 and Starburst99 stellar population models with simple parametric star-formation-histories (SFH), we find that even with models that account for effects of stellar evolution with binaries/stellar rotation, model galaxies at $log_{10}(ξ_{ion}[Hz/erg])\lesssim25.0$ have low H$α$ equivalent widths (EW) and redder colors compared to our $z\sim2$ observed sample. We find that introducing star-bursts to the SFHs resolve the tension with the models, however, due to the rapid time evolution of $ξ_{ion}$, H$α$ EWs, and rest-frame optical colors, our Monte Carlo simulations of star-bursts show that random distribution of star-bursts in evolutionary time of galaxies are unlikely to explain the observed distribution. Thus, either our observed sample is specially selected based on their past SFH or stellar models require additional mechanisms to reproduce the observed high UV luminosity of galaxies for a given production rate of hydrogen ionizing photons. △ Less

Submitted 8 January, 2020; originally announced January 2020.

Comments: Accepted for publication in the ApJ, 8 figures, 15 pages

The colors and sizes of recently quenched galaxies: a result of compact starburst before quenching

Authors: Po-Feng Wu, Arjen van der Wel, Rachel Bezanson, Anna Gallazzi, Camilla Pacifici, Caroline M. S. Straatman, Ivana Barisic, Eric F. Bell, Priscilla Chauke, Francesco D'Eugenio, Marijn Franx, Adam Muzzin, David Sobral, Josha van Houdt

Abstract: We analyze the colors and sizes of 32 quiescent (UVJ-selected) galaxies with strong Balmer absorption ($\mbox{EW}(Hδ) \geq 4$Å) at $z\sim0.8$ drawn from DR2 of the LEGA-C survey to test the hypothesis that these galaxies experienced compact, central starbursts before quenching. These recently quenched galaxies, usually referred to as post-starburst galaxies, span a wide range of colors and we find… ▽ More We analyze the colors and sizes of 32 quiescent (UVJ-selected) galaxies with strong Balmer absorption ($\mbox{EW}(Hδ) \geq 4$Å) at $z\sim0.8$ drawn from DR2 of the LEGA-C survey to test the hypothesis that these galaxies experienced compact, central starbursts before quenching. These recently quenched galaxies, usually referred to as post-starburst galaxies, span a wide range of colors and we find a clear correlation between color and half-light radius, such that bluer galaxies are smaller. We build simple toy models to explain this correlation: a normal star-forming disk plus a central, compact starburst component. Bursts with exponential decay timescale of $\sim$~100 Myr that produce $\sim10\%$ to more than 100\% of the pre-existing masses can reproduce the observed correlation. More significant bursts also produce bluer and smaller descendants. Our findings imply that when galaxies shut down star formation rapidly, they generally had experienced compact, starburst events and that the large, observed spread in sizes and colors mostly reflects a variety of burst strengths. Recently quenched galaxies should have younger stellar ages in the centers; multi-wavelength data with high spatial resolution are required to reveal the age gradient. Highly dissipative processes should be responsible for this type of formation history. While determining the mechanisms for individual galaxies is challenging, some recently quenched galaxies show signs of gravitational interactions, suggesting that mergers are likely an important mechanism in triggering the rapid shut-down of star-formation activities at $z\sim0.8$. △ Less

Submitted 8 December, 2019; originally announced December 2019.

Comments: 16 pages, 8 figures, 1 table, accepted by ApJ

A Tale of Two Clusters: An Analysis of Gas-Phase Metallicity and Nebular Gas Conditions in Proto-cluster Galaxies at z~2

Authors: Leo Y. Alcorn, Anshu Gupta, Kim-Vy Tran, Glenn G. Kacprzak, Tiantian Yuan, Jonathan Cohn, Ben Forrest, Karl Glazebrook, Anishya Harshan, Lisa J. Kewley, Ivo Labbé, Themiya Nanayakkara, Casey Papovich, Lee R. Spitler, Caroline M. S. Straatman

Abstract: The ZFIRE survey has spectroscopically confirmed two proto-clusters using the MOSFIRE instrument on Keck 1: one at z=2.095 in COSMOS and another at z=1.62 in UDS. Here we use an updated ZFIRE dataset to derive the properties of ionized gas regions of proto-cluster galaxies by extracting fluxes from emission lines H$β$ 4861Å, [OIII] 5007Å, H$α$ 6563Å, [NII] 6585Å, and [SII] 6716,6731Å. We measure g… ▽ More The ZFIRE survey has spectroscopically confirmed two proto-clusters using the MOSFIRE instrument on Keck 1: one at z=2.095 in COSMOS and another at z=1.62 in UDS. Here we use an updated ZFIRE dataset to derive the properties of ionized gas regions of proto-cluster galaxies by extracting fluxes from emission lines H$β$ 4861Å, [OIII] 5007Å, H$α$ 6563Å, [NII] 6585Å, and [SII] 6716,6731Å. We measure gas-phase metallicity of members in both proto-clusters using two indicators, including a strong-line indicator relatively independent of ionization parameter and electron density. Proto-cluster and field galaxies in both UDS and COSMOS lie on the same Mass-Metallicity Relation with both metallicity indicators. We compare our results to recent IllustrisTNG results, which reports no significant gas-phase metallicity offset between proto-cluster and field galaxies until z=1.5. This is in agreement with our observed metallicities, where no offset is measured between proto-cluster and field populations. We measure tentative evidence from stacked spectra that UDS high mass proto-cluster and field galaxies have differing [OIII]/H$β$ ratios, however these results are dependent on the sample size of the high mass stacks. △ Less

Submitted 19 August, 2019; originally announced August 2019.

Comments: 18 pages, 6 figures, accepted by the Astrophysical Journal

Rejuvenation in $z\sim0.8$ quiescent galaxies in LEGA-C

Authors: Priscilla Chauke, Arjen van der Wel, Camilla Pacifici, Rachel Bezanson, Po-Feng Wu, Anna Gallazzi, Caroline Straatman, Marijn Franx, Ivana Barišić, Eric F. Bell, Josha van Houdt, Michael V. Maseda, Adam Muzzin, David Sobral, Justin Spilker

Abstract: We use reconstructed star-formation histories (SFHs) of quiescent galaxies at $z=0.6-1$ in the LEGA-C survey to identify secondary star-formation episodes that, after an initial period of quiescence, moved the galaxies back to the star-forming main sequence (blue cloud). $16\pm3$\% of the $z\sim0.8$ quiescent population has experienced such rejuvenation events in the redshift range $0.7<z<1.5$ aft… ▽ More We use reconstructed star-formation histories (SFHs) of quiescent galaxies at $z=0.6-1$ in the LEGA-C survey to identify secondary star-formation episodes that, after an initial period of quiescence, moved the galaxies back to the star-forming main sequence (blue cloud). $16\pm3$\% of the $z\sim0.8$ quiescent population has experienced such rejuvenation events in the redshift range $0.7<z<1.5$ after reaching quiescence at some earlier time. On average, these galaxies first became quiescent at $z=1.2$, and those that rejuvenated, remained quiescent for $\sim1$Gyr before their secondary SF episode which lasted $\sim0.7$Gyr. The stellar mass attributed to rejuvenation is on average 10\% of the galaxy stellar mass, with rare instances of an increase of more than a factor 2. Overall, rejuvenation events only contribute $\sim2$\% of the total stellar mass in $z\sim0.8$ quiescent galaxies and we conclude that rejuvenation is not an important evolutionary channel when considering the growth of the red sequence. However, our results complicate the interpretation of galaxy demographics in color space: the galaxies with rejuvenation events tend to lie in the so-called `green valley', yet their progenitors were quiescent at $z\sim2$. △ Less

Submitted 18 May, 2019; originally announced May 2019.

Comments: 14 pages, 8 figures, accepted by ApJ

1D Kinematics from stars and ionized gas at $z\sim0.8$ from the LEGA-C spectroscopic survey of massive galaxies

Authors: Rachel Bezanson, Arjen van der Wel, Caroline Straatman, Camilla Pacifici, Po-Feng Wu, Ivana Barišić, Eric F. Bell, Charlie Conroy, Francesco D'Eugenio, Marijn Franx, Anna Gallazzi, Josha van Houdt, Michael V. Maseda, Adam Muzzin, Jesse van de Sande, David Sobral, Justin Spilker

Abstract: We present a comparison of the observed, spatially integrated stellar and ionized gas velocity dispersions of $\sim1000$ massive ($\log M_{\star}/M_{\odot}\gtrsim\,10.3$) galaxies in the Large Early Galaxy Astrophysics Census (LEGA-C) survey at $0.6\lesssim\,z\lesssim1.0$. The high $S/N\sim20{\rmÅ^{-1}}$ afforded by 20 hour VLT/VIMOS spectra allows for joint modeling of the stellar continuum and e… ▽ More We present a comparison of the observed, spatially integrated stellar and ionized gas velocity dispersions of $\sim1000$ massive ($\log M_{\star}/M_{\odot}\gtrsim\,10.3$) galaxies in the Large Early Galaxy Astrophysics Census (LEGA-C) survey at $0.6\lesssim\,z\lesssim1.0$. The high $S/N\sim20{\rmÅ^{-1}}$ afforded by 20 hour VLT/VIMOS spectra allows for joint modeling of the stellar continuum and emission lines in all galaxies, spanning the full range of galaxy colors and morphologies. These observed integrated velocity dispersions (denoted as $σ'_{g, int}$ and $σ'_{\star, int}$) are related to the intrinsic velocity dispersions of ionized gas or stars, but also include rotational motions through beam smearing and spectral extraction. We find good average agreement between observed velocity dispersions, with $\langle\log(σ'_{g, int}/σ'_{\star, int})\rangle=-0.003$. This result does not depend strongly on stellar population, structural properties, or alignment with respect to the slit. However, in all regimes we find significant scatter between $σ'_{g, int}$ and $σ'_{\star, int}$, with an overall scatter of 0.13 dex of which 0.05 dex is due to observational uncertainties. For an individual galaxy, the scatter between $σ'_{g, int}$ and $σ'_{\star, int}$ translates to an additional uncertainty of $\sim0.24\rm{dex}$ on dynamical mass derived from $σ'_{g, int}$, on top of measurement errors and uncertainties from Virial constant or size estimates. We measure the $z\sim0.8$ stellar mass Faber-Jackson relation and demonstrate that emission line widths can be used to measure scaling relations. However, these relations will exhibit increased scatter and slopes that are artificially steepened by selecting on subsets of galaxies with progressively brighter emission lines. △ Less

Submitted 19 November, 2018; originally announced November 2018.

Comments: 9 pages, 5 figures, Accepted for publication in ApJ Letters

The Large Early Galaxy Astrophysics Census (LEGA-C) Data Release II: dynamical and stellar population properties of z ~< 1 galaxies in the COSMOS field

Authors: Caroline M. S. Straatman, Arjen van der Wel, Rachel Bezanson, Camilla Pacifici, Anna Gallazzi, Po-Feng Wu, Kai Noeske, Ivana Barisic, Eric F. Bell, Gabriel B. Brammer, Joao Calhau, Priscilla Chauke, Marijn Franx, Josha van Houdt, Ivo Labbe, Michael V. Maseda, Juan C. Munoz-Mateos, Adam Muzzin, Jesse van de Sande, David Sobral, Justin S. Spilker

Abstract: We present the second data release of the Large Early Galaxy Astrophysics Census (LEGA-C), an ESO 130-night public spectroscopic survey conducted with VIMOS on the Very Large Telescope. We release 1988 spectra with typical continuum S / N ~= 20 /Angstrom of galaxies at 0.6 ~< z ~< 1.0, each observed for ~20 hours and fully reduced with a custom-built pipeline. We also release a catalog with spectr… ▽ More We present the second data release of the Large Early Galaxy Astrophysics Census (LEGA-C), an ESO 130-night public spectroscopic survey conducted with VIMOS on the Very Large Telescope. We release 1988 spectra with typical continuum S / N ~= 20 /Angstrom of galaxies at 0.6 ~< z ~< 1.0, each observed for ~20 hours and fully reduced with a custom-built pipeline. We also release a catalog with spectroscopic redshifts, emission line fluxes, Lick/IDS indices, and observed stellar and gas velocity dispersions that are spatially integrated quantities including both rotational motions and genuine dispersion. To illustrate the new parameter space in the intermediate redshift regime probed by LEGA-C we explore relationships between dynamical and stellar population properties. The star-forming galaxies typically have observed stellar velocity dispersions of ~150 km/s and strong Hdelta absorption (Hd_A ~ 5 Angstrom), while passive galaxies have higher observed stellar velocity dispersions ~200 km/s and weak Hdelta absortion (Hd_A ~ 0 Angstrom). Strong [O III]5007 / Hbeta ratios tend to occur mostly for galaxies with weak Hd_A or galaxies with higher observed velocity dispersion. Beyond these broad trends, we find a large diversity of possible combinations of rest-frame colors, absorption line strengths and emission line detections, illustrating the utility of spectroscopic measurements to more accurately understand galaxy evolution. By making the spectra and value-added catalogs publicly available we encourage the community to take advantage of this very substantial investment in telescope time provided by ESO. △ Less

Submitted 21 September, 2018; originally announced September 2018.

Comments: accepted for publication by ApJ; LEGA-C data-release http://www.mpia.de/home/legac/

Fast and slow paths to quiescence: ages and sizes of 400 quiescent galaxies from the LEGA-C survey

Authors: Po-Feng Wu, Arjen van der Wel, Rachel Bezanson, Anna Gallazzi, Camilla Pacifici, Caroline M. S. Straatman, Ivana Barisic, Eric F. Bell, Priscilla Chauke, Josha van Houdt, Marijn Franx, Adam Muzzin, David Sobral, Vivienne Wild

Abstract: We analyze stellar age indicators (D$_n$4000 and EW(H$δ$)) and sizes of 467 quiescent galaxies with $M_\ast \geq 10^{10} M_\odot$ at $z\sim0.7$ drawn from DR2 of the LEGA-C survey. Interpreting index variations in terms of equivalent single stellar population age, we find that the median stellar population is younger for larger galaxies at fixed stellar mass. The effect is significant, yet small;… ▽ More We analyze stellar age indicators (D$_n$4000 and EW(H$δ$)) and sizes of 467 quiescent galaxies with $M_\ast \geq 10^{10} M_\odot$ at $z\sim0.7$ drawn from DR2 of the LEGA-C survey. Interpreting index variations in terms of equivalent single stellar population age, we find that the median stellar population is younger for larger galaxies at fixed stellar mass. The effect is significant, yet small; the ages of the larger and the smaller subsets differ by only $<500$ Myr, much less than the age variation among individual galaxies ($\sim1.5$ Gyr). At the same time, quiescent galaxies with the strongest H$δ$ absorption --- those experienced recent and rapid quenching events --- tend to be smaller than the average. These co-existing trends unify seemingly contradictory results in the literature; the complex correlations between size and age indicators revealed by our large sample of galaxies with high-quality spectra suggest that there are multiple evolutionary pathways to quiescence. Regardless of the specific physical mechanisms responsible for the cessation of star formation in massive galaxies, the large scatter in D$_n$4000 and EW(H$δ$) immediately implies that galaxies follow a large variety in evolutionary pathways. On the one hand, we see evidence for a process that slowly shuts off star-formation and transforms star-forming galaxies to quiescent galaxies without necessarily changing their structures. On the other hand, there is likely a mechanism that rapidly quenches galaxies, an event that coincides with dramatic structural changes, producing small post-starburst galaxies. △ Less

Submitted 11 October, 2018; v1 submitted 4 September, 2018; originally announced September 2018.

Comments: Main part 16 pages, 13 figures. Appendix: 8 figures. Accepted by ApJ

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

Near infrared spectroscopy and star-formation histories of 3<z<4 quiescent galaxies

Authors: C. Schreiber, K. Glazebrook, T. Nanayakkara, G. G. Kacprzak, I. Labbe, P. Oesch, T. Yuan, K. -V. Tran, C. Papovich, L. Spitler, C. Straatman

Abstract: We present Keck-MOSFIRE H and K spectra for a sample of 24 candidate quiescent galaxies (QGs) at 3<z<4, identified from UVJ colors and photometric redshifts in the ZFOURGE and 3DHST surveys. We obtain spectroscopic redshifts for half of the sample, using absorption or emission lines, and confirm the high accuracy of the photometric redshifts with a median error of 1.2%. Two galaxies turn out to be… ▽ More We present Keck-MOSFIRE H and K spectra for a sample of 24 candidate quiescent galaxies (QGs) at 3<z<4, identified from UVJ colors and photometric redshifts in the ZFOURGE and 3DHST surveys. We obtain spectroscopic redshifts for half of the sample, using absorption or emission lines, and confirm the high accuracy of the photometric redshifts with a median error of 1.2%. Two galaxies turn out to be dusty objects at lower redshifts (z<2.5), and are the only two detected in the sub-mm with ALMA. High equivalent-width [OIII] was observed in two galaxies, contributing up to 30% of the K-band flux and mimicking the colors of an old stellar population. This implies a failure rate of only 20% for the UVJ selection at these redshifts. Balmer absorption was identified in 4 of the brighest galaxies, confirming the absence of OB stars. Modeling all QGs with a wide range of star-formation histories, we find sSFR a factor of 10 below the main sequence (MS) for all but one galaxy, and less than 0.01 Gyr$^{-1}$ for half of the sample. This is consistent with the H$β$ and [OII] luminosities, and the ALMA non-detections. We then find that these QGs have quenched on average 300 Myr before observation, between z=3.5 and 5, and that they formed at z~5.5 with a mean SFR~300 Msun/yr. Considering an alternative selection of QGs based solely on the sSFR from SED modeling, we find that galaxies a factor 10 below the MS are 40% more numerous than UVJ-quiescent galaxies, implying that the UVJ selection is pure but incomplete. Current models fail at reproducing our observations and underestimate either the number density of QGs by more than an order of magnitude or the duration of their quiescence by a factor two. Overall, these results confirm the existence of an unexpected population of QGs at z>3, and offer the first insights on their formation history. [abridged] △ Less

Submitted 6 July, 2018; originally announced July 2018.

Comments: 30 pages (+ appendix), 18 figures, accepted for publication in A&A

Journal ref: A&A 618, A85 (2018)

Molecular Gas Contents and Scaling Relations for Massive Passive Galaxies at Intermediate Redshifts from the LEGA-C Survey

Authors: Justin Spilker, Rachel Bezanson, Ivana Barisic, Eric Bell, Claudia del P. Lagos, Michael Maseda, Adam Muzzin, Camilla Pacifici, David Sobral, Caroline Straatman, Arjen van der Wel, Pieter van Dokkum, Benjamin Weiner, Katherine Whitaker, Christina C. Williams, Po-Feng Wu

Abstract: A decade of study has established that the molecular gas properties of star-forming galaxies follow coherent scaling relations out to z~3, suggesting remarkable regularity of the interplay between molecular gas, star formation, and stellar growth. Passive galaxies, however, are expected to be gas-poor and therefore faint, and thus little is known about molecular gas in passive galaxies beyond the… ▽ More A decade of study has established that the molecular gas properties of star-forming galaxies follow coherent scaling relations out to z~3, suggesting remarkable regularity of the interplay between molecular gas, star formation, and stellar growth. Passive galaxies, however, are expected to be gas-poor and therefore faint, and thus little is known about molecular gas in passive galaxies beyond the local universe. Here we present deep Atacama Large Millimeter/submillimeter Array (ALMA) observations of CO(2-1) emission in 8 massive (Mstar ~ 10^11 Msol) galaxies at z~0.7 selected to lie a factor of 3-10 below the star-forming sequence at this redshift, drawn from the Large Early Galaxy Astrophysics Census (LEGA-C) survey. We significantly detect half the sample, finding molecular gas fractions <~0.1. We show that the molecular and stellar rotational axes are broadly consistent, arguing that the molecular gas was not accreted after the galaxies became quiescent. We find that scaling relations extrapolated from the star-forming population over-predict both the gas fraction and gas depletion time for passive objects, suggesting the existence of either a break or large increase in scatter in these relations at low specific star formation rate. Finally, we show that the gas fractions of the passive galaxies we have observed at intermediate redshifts are naturally consistent with evolution into local massive early-type galaxies by continued low-level star formation, with no need for further gas accretion or dynamical stabilization of the gas reservoirs in the intervening 6 billion years. △ Less

Submitted 7 May, 2018; originally announced May 2018.

Comments: 23 pages, 9 figures, accepted to ApJ. Key results in Figures 7, 8, and 9

SFHs OF $Z\sim1$ Galaxies in LEGA-C

Authors: Priscilla Chauke, Arjen van der Wel, Camilla Pacifici, Rachel Bezanson, Po-Feng Wu, Anna Gallazzi, Kai Noeske, Caroline Straatman, Juan-Carlos Muños-Mateos, Marijn Franx, Ivana Barišić, Eric F. Bell, Gabriel B. Brammer, Joao Calhau, Josha van Houdt, Ivo Labbé, Michael V. Maseda, Adam Muzzin, Hans-Walter Rix, David Sobral

Abstract: Using high resolution spectra from the VLT LEGA-C program, we reconstruct the star formation histories (SFHs) of 607 galaxies at redshifts $z = 0.6-1.0$ and stellar masses $\gtrsim10^{10}$M$_{\odot}$ using a custom full spectrum fitting algorithm that incorporates the emcee and FSPS packages. We show that the mass-weighted age of a galaxy correlates strongly with stellar velocity dispersion (… ▽ More Using high resolution spectra from the VLT LEGA-C program, we reconstruct the star formation histories (SFHs) of 607 galaxies at redshifts $z = 0.6-1.0$ and stellar masses $\gtrsim10^{10}$M$_{\odot}$ using a custom full spectrum fitting algorithm that incorporates the emcee and FSPS packages. We show that the mass-weighted age of a galaxy correlates strongly with stellar velocity dispersion ($σ_*$) and ongoing star-formation (SF) activity, with the stellar content in higher-$σ_*$ galaxies having formed earlier and faster. The SFHs of quiescent galaxies are generally consistent with passive evolution since their main SF epoch, but a minority show clear evidence of a rejuvenation event in their recent past. The mean age of stars in galaxies that are star-forming is generally significantly younger, with SF peaking after $z<1.5$ for almost all star-forming galaxies in the sample: many of these still have either constant or rising SFRs on timescales $>100$Myrs. This indicates that $z>2$ progenitors of $z\sim1$ star-forming galaxies are generally far less massive. Finally, despite considerable variance in the individual SFHs, we show that the current SF activity of massive galaxies ($>$L$_*$) at $z\sim1$ correlates with SF levels at least $3$Gyrs prior: SFHs retain `memory' on a large fraction of the Hubble time. Our results illustrate a novel approach to resolve the formation phase of galaxies, and in identifying their individual evolutionary paths, connects progenitors and descendants across cosmic time. This is uniquely enabled by the high-quality continuum spectroscopy provided by the LEGA-C survey. △ Less

Submitted 8 May, 2018; v1 submitted 7 May, 2018; originally announced May 2018.

Comments: 18 pages, 9 figures, accepted by ApJ

ZFIRE: 3D Modeling of Rotation, Dispersion, and Angular Momentum of Star-Forming Galaxies at z~2

Authors: Leo Y. Alcorn, Kim-Vy Tran, Karl Glazebrook, Caroline M. Straatman, Michael Cowley, Ben Forrest, Glenn G. Kacprzak, Lisa J. Kewley, Ivo Labbé, Themiya Nanayakkara, Lee R. Spitler, Adam Tomczak, Tiantian Yuan

Abstract: We perform a kinematic and morphological analysis of 44 star-forming galaxies at $z\sim2$ in the COSMOS legacy field using near-infrared spectroscopy from Keck/MOSFIRE and F160W imaging from CANDELS/3D-HST as part of the ZFIRE survey. Our sample consists of cluster and field galaxies from $2.0 < z < 2.5$ with K band multi-object slit spectroscopic measurements of their H$α$ emission lines. H$α$ ro… ▽ More We perform a kinematic and morphological analysis of 44 star-forming galaxies at $z\sim2$ in the COSMOS legacy field using near-infrared spectroscopy from Keck/MOSFIRE and F160W imaging from CANDELS/3D-HST as part of the ZFIRE survey. Our sample consists of cluster and field galaxies from $2.0 < z < 2.5$ with K band multi-object slit spectroscopic measurements of their H$α$ emission lines. H$α$ rotational velocities and gas velocity dispersions are measured using the Heidelberg Emission Line Algorithm (HELA), which compares directly to simulated 3D data-cubes. Using a suite of simulated emission lines, we determine that HELA reliably recovers input S$_{0.5}$ and angular momentum at small offsets, but $V_{2.2}/σ_g$ values are offset and highly scattered. We examine the role of regular and irregular morphology in the stellar mass kinematic scaling relations, deriving the kinematic measurement S$_{0.5}$, and finding $\log(S_{0.5}) = (0.38\pm0.07)\log(M/M_{\odot}-10) + (2.04\pm0.03)$ with no significant offset between morphological populations and similar levels of scatter ($\sim0.16$ dex). Additionally, we identify a correlation between M$_{\star}$ and $V_{2.2}/σ_g$ for the total sample, showing an increasing level of rotation dominance with increasing M$_{\star}$, and a high level of scatter for both regular and irregular galaxies. We estimate the specific angular momenta ($j_{disk}$) of these galaxies and find a slope of $0.36\pm0.12$, shallower than predicted without mass-dependent disk growth, but this result is possibly due to measurement uncertainty at M$_{\star}$ $<$ 9.5. However, through a K-S test we find irregular galaxies to have marginally higher $j_{disk}$ values than regular galaxies, and high scatter at low masses in both populations. △ Less

Submitted 10 April, 2018; originally announced April 2018.

Comments: Accepted to ApJ

Spatially Resolved Stellar Kinematics from LEGA-C: Increased Rotational Support in z~0.8 Quiescent Galaxies

Authors: Rachel Bezanson, Arjen van der Wel, Camilla Pacifici, Kai Noeske, Ivana Barišić, Eric F. Bell, Gabriel B. Brammer, Joao Calhau, Priscilla Chauke, Pieter van Dokkum, Marijn Franx, Anna Gallazzi, Josha van Houdt, Ivo Labbé, Michael V. Maseda, Juan Carlos Muños-Mateos, Adam Muzzin, Jesse van de Sande, David Sobral, Caroline Straatman, Po-Feng Wu

Abstract: We present stellar rotation curves and velocity dispersion profiles for 104 quiescent galaxies at $z=0.6-1$ from the Large Early Galaxy Astrophysics Census (LEGA-C) spectroscopic survey. Rotation is typically probed across 10-20kpc, or to an average of 2.7${\rm R_e}$. Combined with central stellar velocity dispersions ($σ_0$) this provides the first determination of the dynamical state of a sample… ▽ More We present stellar rotation curves and velocity dispersion profiles for 104 quiescent galaxies at $z=0.6-1$ from the Large Early Galaxy Astrophysics Census (LEGA-C) spectroscopic survey. Rotation is typically probed across 10-20kpc, or to an average of 2.7${\rm R_e}$. Combined with central stellar velocity dispersions ($σ_0$) this provides the first determination of the dynamical state of a sample selected by a lack of star formation activity at large lookback time. The most massive galaxies ($M_{\star}>2\times10^{11}\,M_{\odot}$) generally show no or little rotation measured at 5kpc ($|V_5|/σ_0<0.2$ in 8 of 10 cases), while ${\sim}64\%$ of less massive galaxies show significant rotation. This is reminiscent of local fast- and slow-rotating ellipticals and implies that low- and high-redshift quiescent galaxies have qualitatively similar dynamical structures. We compare $|V_5|/σ_0$ distributions at $z\sim0.8$ and the present day by re-binning and smoothing the kinematic maps of 91 low-redshift quiescent galaxies from the CALIFA survey and find evidence for a decrease in rotational support since $z\sim1$. This result is especially strong when galaxies are compared at fixed velocity dispersion; if velocity dispersion does not evolve for individual galaxies then the rotational velocity at 5kpc was an average of ${94\pm22\%}$ higher in $z\sim0.8$ quiescent galaxies than today. Considering that the number of quiescent galaxies grows with time and that new additions to the population descend from rotationally-supported star-forming galaxies, our results imply that quiescent galaxies must lose angular momentum between $z\sim1$ and the present, presumably through dissipationless merging, and/or that the mechanism that transforms star-forming galaxies also reduces their rotational support. △ Less

Submitted 6 April, 2018; originally announced April 2018.

Comments: 30 pages, 18 figures, Accepted for Publication in ApJ

Stellar Populations of over one thousand $z\sim0.8$ Galaxies from LEGA-C: Ages and Star Formation Histories from D$_n$4000 and H$δ$

Authors: Po-Feng Wu, Arjen van der Wel, Anna Gallazzi, Rachel Bezanson, Camilla Pacifici, Caroline Straatman, Marijn Franx, Ivana Barišić, Eric F. Bell, Gabriel B. Brammer, Joao Calhau, Priscilla Chauke, Josha van Houdt, Michael V. Maseda, Adam Muzzin, Hans-Walter Rix, David Sobral, Justin Spilker, Jesse van de Sande, Pieter van Dokkum, Vivienne Wild

Abstract: Drawing from the LEGA-C dataset, we present the spectroscopic view of the stellar population across a large volume- and mass-selected sample of galaxies at large lookback time. We measure the 4000Å break (D$_n$4000) and Balmer absorption line strengths (probed by H$δ$) from 1019 high-quality spectra of $z=0.6 - 1.0$ galaxies with $M_\ast = 2 \times 10^{10} M_\odot - 3 \times 10^{11} M_\odot$. Our… ▽ More Drawing from the LEGA-C dataset, we present the spectroscopic view of the stellar population across a large volume- and mass-selected sample of galaxies at large lookback time. We measure the 4000Å break (D$_n$4000) and Balmer absorption line strengths (probed by H$δ$) from 1019 high-quality spectra of $z=0.6 - 1.0$ galaxies with $M_\ast = 2 \times 10^{10} M_\odot - 3 \times 10^{11} M_\odot$. Our analysis serves as a first illustration of the power of high-resolution, high-S/N continuum spectroscopy at intermediate redshifts as a qualitatively new tool to constrain galaxy formation models. The observed D$_n$4000-EW(H$δ$) distribution of our sample overlaps with the distribution traced by present-day galaxies, but $z\sim 0.8$ galaxies populate that locus in a fundamentally different manner. While old galaxies dominate the present-day population at all stellar masses $> 2\times10^{10} M_\odot$, we see a bimodal D$_n$4000-EW(H$δ$) distribution at $z\sim0.8$, implying a bimodal light-weighted age distribution. The light-weighted age depends strongly on stellar mass, with the most massive galaxies $>1\times10^{11}M_\odot$ being almost all older than 2 Gyr. At the same time we estimate that galaxies in this high mass range are only $\sim3$ Gyr younger than their $z\sim0.1$ counterparts, at odd with pure passive evolution given a difference in lookback time of $>5$ Gyr; younger galaxies must grow to $>10^{11}M_\odot$ in the meantime, and/or small amounts of young stars must keep the light-weighted ages young. Star-forming galaxies at $z\sim0.8$ have stronger H$δ$ absorption than present-day galaxies with the same D$_n$4000, implying larger short-term variations in star-formation activity. △ Less

Submitted 19 February, 2018; originally announced February 2018.

Comments: 17 pages, 9 figures, accepted by 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

Decoupled Black Hole Accretion and Quenching: The Relationship Between BHAR, SFR, and Quenching in Milky Way and Andromeda-mass Progenitors Since z = 2.5

Authors: Michael J. Cowley, Lee R. Spitler, Ryan F. Quadri, Andy D. Goulding, Casey Papovich, Kim-Vy H. Tran, Ivo Labbe, Leo Alcorn, Rebecca J. Allen, Ben Forrest, Karl Glazebrook, Glenn G. Kacprzak, Glenn Morrison, Themiya Nanayakkara, Caroline M. S. Straatman, Adam R. Tomczak

Abstract: We investigate the relationship between the black hole accretion rate (BHAR) and star-formation rate (SFR) for Milky Way (MW) and Andromeda (M31)-mass progenitors from z = 0.2 - 2.5. We source galaxies from the Ks-band selected ZFOURGE survey, which includes multi-wavelenth data spanning 0.3 - 160um. We use decomposition software to split the observed SEDs of our galaxies into their active galacti… ▽ More We investigate the relationship between the black hole accretion rate (BHAR) and star-formation rate (SFR) for Milky Way (MW) and Andromeda (M31)-mass progenitors from z = 0.2 - 2.5. We source galaxies from the Ks-band selected ZFOURGE survey, which includes multi-wavelenth data spanning 0.3 - 160um. We use decomposition software to split the observed SEDs of our galaxies into their active galactic nuclei (AGN) and star-forming components, which allows us to estimate BHARs and SFRs from the infrared (IR). We perform tests to check the robustness of these estimates, including a comparison to BHARs and SFRs derived from X-ray stacking and far-IR analysis, respectively. We find as the progenit- ors evolve, their relative black hole-galaxy growth (i.e. their BHAR/SFR ratio) increases from low to high redshift. The MW-mass progenitors exhibit a log-log slope of 0.64 +/- 0.11, while the M31-mass progenitors are 0.39 +/- 0.08. This result contrasts with previous studies that find an almost flat slope when adopting X-ray/AGN-selected or mass-limited samples and is likely due to their use of a broad mixture of galaxies with different evolutionary histories. Our use of progenitor-matched samples highlights the potential importance of carefully selecting progenitors when searching for evolutionary relationships between BHAR/SFRs. Additionally, our finding that BHAR/SFR ratios do not track the rate at which progenitors quench casts doubts over the idea that the suppression of star-formation is predominantly driven by luminous AGN feedback (i.e. high BHARs). △ Less

Submitted 11 October, 2017; originally announced October 2017.

Comments: 7 pages, 5 figures, Accepted for publication in MNRAS

Jekyll & Hyde: quiescence and extreme obscuration in a pair of massive galaxies 1.5 Gyr after the Big Bang

Authors: Corentin Schreiber, Ivo Labbé, Karl Glazebrook, Georgios Bekiaris, Casey Papovich, Tiago Costa, David Elbaz, Glenn G. Kacprzak, Themiya Nanayakkara, Pascal Oesch, Maurilio Pannella, Lee Spitler, Caroline Straatman, Kim-Vy Tran, Tao Wang

Abstract: We obtained ALMA spectroscopy and imaging to investigate the origin of the unexpected sub-mm emission toward the most distant quiescent galaxy known to date, ZF-COSMOS-20115 at z=3.717. We show here that this sub-mm emission is produced by another massive, compact and extremely obscured galaxy, located only 3.1 kpc away from the quiescent galaxy. We dub the quiescent and dusty galaxies Jekyll and… ▽ More We obtained ALMA spectroscopy and imaging to investigate the origin of the unexpected sub-mm emission toward the most distant quiescent galaxy known to date, ZF-COSMOS-20115 at z=3.717. We show here that this sub-mm emission is produced by another massive, compact and extremely obscured galaxy, located only 3.1 kpc away from the quiescent galaxy. We dub the quiescent and dusty galaxies Jekyll and Hyde, respectively. No dust emission is detected at the location of Jekyll, implying SFR < 13 Msun/yr, which is the most stringent upper limit ever obtained for a quiescent galaxy at these redshifts. The two sources are confirmed to lie at the same redshift thanks to the detection of [CII]158 in Hyde, which provides one the few robust redshifts for an "H-dropout" galaxy. The line has a rotating-disk velocity profile blueshifted from Jekyll by 549+/-60 km/s, demonstrating that it is produced by another galaxy. Careful de-blending of the Spitzer imaging confirms the existence of Hyde, and its non-detection with Hubble requires extreme attenuation by dust. Modeling the photometry of both galaxies shows that Jekyll has fully quenched >200 Myr prior to observation and still presents a challenge for models, while Hyde only harbors moderate star-formation (SFR<120 Msun/yr) and is located at least a factor 1.4 below the z~4 main sequence. Hyde could also have stopped forming stars <200 Myr before being observed, which would be consistent with its hight compactness similar to z~4 quiescent galaxies and its low [CII]/FIR ratio, but significant SF cannot be ruled out. Finally, we show that Hyde hosts a dense reservoir of gas comparable to that of extreme starbursts, suggesting that its SFR was reduced without expelling the gas outside of the galaxy. We argue that Jekyll and Hyde can be seen as two stages of quenching, and provide a unique laboratory to study this poorly understood phenomenon. [abridged] △ Less

Submitted 8 December, 2017; v1 submitted 11 September, 2017; originally announced September 2017.

Comments: 22 pages, 16 figures, accepted for publication in A&A

Journal ref: A&A 611, A22 (2018)

Stellar Dynamics and Star-Formation Histories of z $\sim$ 1 Radio-loud Galaxies

Authors: Ivana Barišić, Arjen van der Wel, Rachel Bezanson, Camilla Pacifici, Kai Noeske, Juan C. Muñoz - Mateos, Marijn Franx, Vernesa Smolčić, Eric F. Bell, Gabriel Brammer, João Calhau, Priscilla Chauké, Pieter G. van Dokkum, Josha van Houdt, Anna Gallazzi, Ivo Labbé, Michael V. Maseda, Adam Muzzin, David Sobral, Caroline Straatman, Po-Feng Wu

Abstract: We investigate the stellar kinematics and stellar populations of 58 radio-loud galaxies of intermediate luminosities (L$_{3 GHz}$ $>$ 10$^{23}$ W Hz$^{-1}$ ) at 0.6 < z < 1. This sample is constructed by cross-matching galaxies from the deep VLT/VIMOS LEGA-C spectroscopic survey with the VLA 3 GHz dataset. The LEGA-C continuum spectra reveal for the first time stellar velocity dispersions and age… ▽ More We investigate the stellar kinematics and stellar populations of 58 radio-loud galaxies of intermediate luminosities (L$_{3 GHz}$ $>$ 10$^{23}$ W Hz$^{-1}$ ) at 0.6 < z < 1. This sample is constructed by cross-matching galaxies from the deep VLT/VIMOS LEGA-C spectroscopic survey with the VLA 3 GHz dataset. The LEGA-C continuum spectra reveal for the first time stellar velocity dispersions and age indicators of z $\sim$ 1 radio galaxies. We find that $z\sim 1$ radio-loud AGN occur exclusively in predominantly old galaxies with high velocity dispersions: $σ_*>$ 175 km s$^{-1}$, corresponding to black hole masses in excess of $10^8$ M$_{\odot}$. Furthermore, we confirm that at a fixed stellar mass the fraction of radio-loud AGN at z $\sim$ 1 is 5 - 10 times higher than in the local universe, suggesting that quiescent, massive galaxies at z $\sim$ 1 switch on as radio AGN on average once every Gyr. Our results strengthen the existing evidence for a link between high black-hole masses, radio loudness and quiescence at z $\sim$ 1. △ Less

Submitted 28 August, 2017; originally announced August 2017.

Comments: 10 pages, 8 figures, 1 table. Accepted for publication in The Astrophysical Journal

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

ZFIRE: Using H$α$ equivalent widths to investigate the in situ initial mass function at z~2

Authors: Themiya Nanayakkara, Karl Glazebrook, Glenn G. Kacprzak, Tiantian Yuan, David B. Fisher, Kim-Vy Tran, Lisa Kewley, Lee Spitler, Leo Alcorn, Michael Cowley, Ivo Labbe, Caroline Straatman, Adam Tomczak

Abstract: We use the ZFIRE survey (http://zfire.swinburne.edu.au) to investigate the high mass slope of the initial mass function (IMF) for a mass-complete (log10(M$_*$/M$_\odot$)~9.3) sample of 102 star-forming galaxies at z~2 using their H$α$ equivalent widths (H$α$-EW) and rest-frame optical colours. We compare dust-corrected H$α$-EW distributions with predictions of star-formation histories (SFH) from P… ▽ More We use the ZFIRE survey (http://zfire.swinburne.edu.au) to investigate the high mass slope of the initial mass function (IMF) for a mass-complete (log10(M$_*$/M$_\odot$)~9.3) sample of 102 star-forming galaxies at z~2 using their H$α$ equivalent widths (H$α$-EW) and rest-frame optical colours. We compare dust-corrected H$α$-EW distributions with predictions of star-formation histories (SFH) from PEGASE.2 and Starburst99 synthetic stellar population models. We find an excess of high H$α$-EW galaxies that are up to 0.3--0.5 dex above the model-predicted Salpeter IMF locus and the H$α$-EW distribution is much broader (10--500 Å) than can easily be explained by a simple monotonic SFH with a standard Salpeter-slope IMF. Though this discrepancy is somewhat alleviated when it is assumed that there is no relative attenuation difference between stars and nebular lines, the result is robust against observational biases, and no single IMF (i.e. non-Salpeter slope) can reproduce the data. We show using both spectral stacking and Monte Carlo simulations that starbursts cannot explain the EW distribution. We investigate other physical mechanisms including models with variations in stellar rotation, binary star evolution, metallicity, and the IMF upper-mass cutoff. IMF variations and/or highly rotating extreme metal poor stars (Z~0.1Z$_\odot$) with binary interactions are the most plausible explanations for our data. If the IMF varies, then the highest H$α$-EWs would require very shallow slopes ($Γ$>-1.0) with no one slope able to reproduce the data. Thus, the IMF would have to vary stochastically. We conclude that the stellar populations at z~2 show distinct differences from local populations and there is no simple physical model to explain the large variation in H$α$-EWs at z~2. △ Less

Submitted 3 April, 2017; v1 submitted 13 March, 2017; originally announced March 2017.

Comments: Accepted to MNRAS. 43 pages, 27 Figures. Survey website: http://zfire.swinburne.edu.au/

Discovery of Extreme [OIII]+H$β$ Emitting Galaxies Tracing an Overdensity at z~3.5 in CDF-South

Authors: Ben Forrest, Kim-Vy H. Tran, Adam Broussard, Rebecca J. Allen, Miranda Apfel, Michael J. Cowley, Karl Glazebrook, Glenn G. Kacprzak, Ivo Labbé, Themiya Nanayakkara, Casey Papovich, Ryan F. Quadri, Lee R. Spitler, Caroline M. S. Straatman, Adam Tomczak

Abstract: Using deep multi-wavelength photometry of galaxies from ZFOURGE, we group galaxies at $2.5<z<4.0$ by the shape of their spectral energy distributions (SEDs). We identify a population of galaxies with excess emission in the $K_s$-band, which corresponds to [OIII]+H$β$ emission at $2.95<z<3.65$. This population includes 78% of the bluest galaxies with UV slopes steeper than $β= -2$. We de-redshift a… ▽ More Using deep multi-wavelength photometry of galaxies from ZFOURGE, we group galaxies at $2.5<z<4.0$ by the shape of their spectral energy distributions (SEDs). We identify a population of galaxies with excess emission in the $K_s$-band, which corresponds to [OIII]+H$β$ emission at $2.95<z<3.65$. This population includes 78% of the bluest galaxies with UV slopes steeper than $β= -2$. We de-redshift and scale this photometry to build two composite SEDs, enabling us to measure equivalent widths of these Extreme [OIII]+H$β$ Emission Line Galaxies (EELGs) at $z\sim3.5$. We identify 60 galaxies that comprise a composite SED with [OIII]+H$β$ rest-frame equivalent width of $803\pm228$Å and another 218 galaxies in a composite SED with equivalent width of $230\pm90$Å. These EELGs are analogous to the `green peas' found in the SDSS, and are thought to be undergoing their first burst of star formation due to their blue colors ($β< -1.6$), young ages ($\log(\rm{age}/yr)\sim7.2$), and low dust attenuation values. Their strong nebular emission lines and compact sizes (typically $\sim1.4$ kpc) are consistent with the properties of the star-forming galaxies possibly responsible for reionizing the universe at $z>6$. Many of the EELGs also exhibit Lyman-$α$ emission. Additionally, we find that many of these sources are clustered in an overdensity in the Chandra Deep Field South, with five spectroscopically confirmed members at $z=3.474 \pm 0.004$. The spatial distribution and photometric redshifts of the ZFOURGE population further confirm the overdensity highlighted by the EELGs. △ Less

Submitted 10 March, 2017; originally announced March 2017.

Comments: Accepted for publication in ApJL. 8 pages, 4 figures

ZFIRE: The Evolution of the Stellar Mass Tully-Fisher Relation to Redshift 2.0 < Z < 2.5 with MOSFIRE

Authors: Caroline M. S. Straatman, Karl Glazebrook, Glenn G. Kacprzak, Ivo Labbé, Themiya Nanayakkara, Leo Alcorn, Michael Cowley, Lisa J. Kewley, Lee R. Spitler, Kim-Vy H. Tran, Tiantian Yuan

Abstract: Using observations made with MOSFIRE on Keck I as part of the ZFIRE survey, we present the stellar mass Tully-Fisher relation at 2.0 < z < 2.5. The sample was drawn from a stellar mass limited, Ks-band selected catalog from ZFOURGE over the CANDELS area in the COSMOS field. We model the shear of the Halpha emission line to derive rotational velocities at 2.2X the scale radius of an exponential dis… ▽ More Using observations made with MOSFIRE on Keck I as part of the ZFIRE survey, we present the stellar mass Tully-Fisher relation at 2.0 < z < 2.5. The sample was drawn from a stellar mass limited, Ks-band selected catalog from ZFOURGE over the CANDELS area in the COSMOS field. We model the shear of the Halpha emission line to derive rotational velocities at 2.2X the scale radius of an exponential disk (V2.2). We correct for the blurring effect of a two-dimensional PSF and the fact that the MOSFIRE PSF is better approximated by a Moffat than a Gaussian, which is more typically assumed for natural seeing. We find for the Tully-Fisher relation at 2.0 < z < 2.5 that logV2.2 =(2.18 +/- 0.051)+(0.193 +/- 0.108)(logM/Msun - 10) and infer an evolution of the zeropoint of Delta M/Msun = -0.25 +/- 0.16 dex or Delta M/Msun = -0.39 +/- 0.21 dex compared to z = 0 when adopting a fixed slope of 0.29 or 1/4.5, respectively. We also derive the alternative kinematic estimator S0.5, with a best-fit relation logS0.5 =(2.06 +/- 0.032)+(0.211 +/- 0.086)(logM/Msun - 10), and infer an evolution of Delta M/Msun= -0.45 +/- 0.13 dex compared to z < 1.2 if we adopt a fixed slope. We investigate and review various systematics, ranging from PSF effects, projection effects, systematics related to stellar mass derivation, selection biases and slope. We find that discrepancies between the various literature values are reduced when taking these into account. Our observations correspond well with the gradual evolution predicted by semi-analytic models. △ Less

Submitted 28 February, 2017; originally announced March 2017.

Comments: 21 pages, 14 figures, 1 appendix. Accepted for publication by Apj, February 28, 2017

A massive, quiescent galaxy at redshift of z=3.717

Authors: Karl Glazebrook, Corentin Schreiber, Ivo Labbé, Themiya Nanayakkara, Glenn G. Kacprzak, Pascal A. Oesch, Casey Papovich, Lee R Spitler, Caroline M. S. Straatman, Kim-Vy H. Tran, Tiantian Yuan

Abstract: In the early Universe finding massive galaxies that have stopped forming stars present an observational challenge as their rest-frame ultraviolet emission is negligible and they can only be reliably identified by extremely deep near-infrared surveys. These have revealed the presence of massive, quiescent early-type galaxies appearing in the universe as early as z$\sim$2, an epoch 3 Gyr after the B… ▽ More In the early Universe finding massive galaxies that have stopped forming stars present an observational challenge as their rest-frame ultraviolet emission is negligible and they can only be reliably identified by extremely deep near-infrared surveys. These have revealed the presence of massive, quiescent early-type galaxies appearing in the universe as early as z$\sim$2, an epoch 3 Gyr after the Big Bang. Their age and formation processes have now been explained by an improved generation of galaxy formation models where they form rapidly at z$\sim$3-4, consistent with the typical masses and ages derived from their observations. Deeper surveys have now reported evidence for populations of massive, quiescent galaxies at even higher redshifts and earlier times, however the evidence for their existence, and redshift, has relied entirely on coarsely sampled photometry. These early massive, quiescent galaxies are not predicted by the latest generation of theoretical models. Here, we report the spectroscopic confirmation of one of these galaxies at redshift z=3.717 with a stellar mass of 1.7$\times$10$^{11}$ M$_\odot$ whose absorption line spectrum shows no current star-formation and which has a derived age of nearly half the age of the Universe at this redshift. The observations demonstrates that the galaxy must have quickly formed the majority of its stars within the first billion years of cosmic history in an extreme and short starburst. This ancestral event is similar to those starting to be found by sub-mm wavelength surveys pointing to a possible connection between these two populations. Early formation of such massive systems is likely to require significant revisions to our picture of early galaxy assembly. △ Less

Submitted 23 October, 2017; v1 submitted 6 February, 2017; originally announced February 2017.

Comments: 6 pages, 7 figures. This is the final preprint corresponding closely to the published version. Uploaded 6 months after publication in accordance with Nature policy

Journal ref: Nature 2017, vol. 544, 71-74

The Size Evolution of Star-forming Galaxies Since z~7 Using ZFOURGE

Authors: Rebecca J. Allen, Glenn G. Kacprzak, Karl Glazebrook, Ivo Labbe, Kim-Vy H. Tran, Lee R. Spitler, Michael Cowley, Themiya Nanayakkara, Casey Papovich, Ryan Quadri, Caroline M. S. Straatman, Vithal Tilvi, Pieter van Dokkum

Abstract: For the first time, we present the size evolution of a mass-complete (log(M*/Msol)>10) sample of star-forming galaxies over redshifts z=1-7, selected from the FourStar Galaxy Evolution Survey (ZFOURGE). Observed H-band sizes are measured from the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS) Hubble Space Telescope (HST)/F160W imaging. Distributions of individual galaxy m… ▽ More For the first time, we present the size evolution of a mass-complete (log(M*/Msol)>10) sample of star-forming galaxies over redshifts z=1-7, selected from the FourStar Galaxy Evolution Survey (ZFOURGE). Observed H-band sizes are measured from the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS) Hubble Space Telescope (HST)/F160W imaging. Distributions of individual galaxy masses and sizes illustrate that a clear mass-size relation exists up to z~7. At z~7, we find that the average galaxy size from the mass-size relation is more compact at a fixed mass of log(M*/Msol)=10.1, with r_1/2,maj=1.02+/-0.29 kpc, than at lower redshifts. This is consistent with our results from stacking the same CANDELS HST/F160W imaging, when we correct for galaxy position angle alignment. We find that the size evolution of star-forming galaxies is well fit by a power law of the form r_e = 7.07(1 + z)^-0.89 kpc, which is consistent with previous works for normal star-formers at 1<z<4. In order to compare our slope with those derived Lyman break galaxy studies, we correct for different IMFs and methodology and find a slope of -0.97+/-0.02, which is shallower than that reported for the evolution of Lyman break galaxies at z>4 (r_e\propto(1 +z)^-1.2+/-0.06). Therefore, we conclude the Lyman break galaxies likely represent a subset of highly star-forming galaxies that exhibit rapid size growth at z>4. △ Less

Submitted 15 December, 2016; originally announced December 2016.

Comments: 7 pages, 4 figures, accepted to Apj letters

ZFIRE: Similar Stellar Growth in H$α$-emitting Cluster and Field Galaxies at z~2

Authors: Kim-Vy H. Tran, Leo Y. Alcorn, Glenn G. Kacprzak, Themiya Nanayakkara, Caroline Straatman, Tiantian Yuan, Michael Cowley, Romeel Davé, Karl Glazebrook, Lisa J. Kewley, Ivo Labbé, Davidé Martizzi, Casey Papovich, Ryan Quadri, Lee R. Spitler, Adam Tomczak

Abstract: We compare galaxy scaling relations as a function of environment at $z\sim2$ with our ZFIRE survey where we have measured H$α$ fluxes for 90 star-forming galaxies selected from a mass-limited [$\log(M_{\star}/M_{\odot})>9$] sample based on ZFOURGE. The cluster galaxies (37) are part of a confirmed system at z=2.095 and the field galaxies (53) are at $1.9<z<2.4$; all are in the COSMOS legacy field.… ▽ More We compare galaxy scaling relations as a function of environment at $z\sim2$ with our ZFIRE survey where we have measured H$α$ fluxes for 90 star-forming galaxies selected from a mass-limited [$\log(M_{\star}/M_{\odot})>9$] sample based on ZFOURGE. The cluster galaxies (37) are part of a confirmed system at z=2.095 and the field galaxies (53) are at $1.9<z<2.4$; all are in the COSMOS legacy field. There is no statistical difference between H$α$-emitting cluster and field populations when comparing their star formation rate (SFR), stellar mass ($M_{\star}$), galaxy size ($r_{eff}$), SFR surface density [$Σ$(H$α_{star}$)], and stellar age distributions. The only difference is that at fixed stellar mass, the H$α$-emitting cluster galaxies are $\log(r_{eff})\sim0.1$ larger than in the field. Approximately 19% of the H$α$-emitters in the cluster and 26% in the field are IR-luminous ($L_{IR}>2\times10^{11} L_{\odot}$). Because the LIRGs in our combined sample are $\sim5$ times more massive than the low-IR galaxies, their radii are $\sim70$% larger. To track stellar growth, we separate galaxies into those that lie above, on, and below the H$α$ star-forming main sequence (SFMS) using $Δ$SFR$(M_{\star})=\pm0.2$ dex. Galaxies above the SFMS (starbursts) tend to have higher H$α$ SFR surface densities and younger light-weighted stellar ages compared to galaxies below the SFMS. Our results indicate that starbursts (+SFMS) in the cluster and field at $z\sim2$ are growing their stellar cores. Lastly, we compare to the (SFR-$M_{\star}$) relation from RHAPSODY cluster simulations and find the predicted slope is nominally consistent with the observations. However, the predicted cluster SFRs tend to be too low by a factor of $\sim2$ which seems to be a common problem for simulations across environment. △ Less

Submitted 5 December, 2016; originally announced December 2016.

Comments: ApJ in press; full version of Table 1 available from ApJ and upon request. Survey websites are http://zfire.swinburne.edu.au and http://zfourge.tamu.edu

Large Molecular Gas Reservoirs in Ancestors of Milky Way-Mass Galaxies 9 Billion Years Ago

Authors: Casey Papovich, Ivo Labbé, Karl Glazebrook, Ryan Quadri, Georgios Bekiaris, Mark Dickinson, Steven Finkelstein, David Fisher, Hanae Inami, Rachael Livermore, Lee Spitler, Caroline Straatman, Kim-Vy Tran

Abstract: The gas accretion and star-formation histories of galaxies like the Milky Way remain an outstanding problem in astrophysics. Observations show that 8 billion years ago, the progenitors to Milky Way-mass galaxies were forming stars 30 times faster than today and predicted to be rich in molecular gas, in contrast with low present-day gas fractions ($<$10%). Here we show detections of molecular gas f… ▽ More The gas accretion and star-formation histories of galaxies like the Milky Way remain an outstanding problem in astrophysics. Observations show that 8 billion years ago, the progenitors to Milky Way-mass galaxies were forming stars 30 times faster than today and predicted to be rich in molecular gas, in contrast with low present-day gas fractions ($<$10%). Here we show detections of molecular gas from the CO(J=3-2) emission (rest-frame 345.8 GHz) in galaxies at redshifts z=1.2-1.3, selected to have the stellar mass and star-formation rate of the progenitors of today's Milky Way-mass galaxies. The CO emission reveals large molecular gas masses, comparable to or exceeding the galaxy stellar masses, and implying most of the baryons are in cold gas, not stars. The galaxies' total luminosities from star formation and CO luminosities yield long gas-consumption timescales. Compared to local spiral galaxies, the star-formation efficiency, estimated from the ratio of total IR luminosity to CO emission,} has remained nearly constant since redshift z=1.2, despite the order of magnitude decrease in gas fraction, consistent with results for other galaxies at this epoch. Therefore the physical processes that determine the rate at which gas cools to form stars in distant galaxies appear to be similar to that in local galaxies. △ Less

Submitted 17 October, 2016; originally announced October 2016.

Comments: To appear in Nature Astronomy, 9 pages

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/