Source author record

Casey Papovich

Casey Papovich appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

89works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

89 published item(s)

preprint2026arXiv

From "The Cliff" to "Virgil": Mapping the Spectral Diversity of Little Red Dots with JWST/NIRSpec

One of JWST's most unexpected discoveries is the emergence of "Little Red Dots'' (LRDs): compact sources at $z \gtrsim 3$ with blue rest-frame UV continua, red optical slopes, and broad Balmer emission lines that challenge standard models and suggest a population of early, unusual active galactic nuclei (AGNs). Using a comprehensive photometric selection and public NIRSpec/PRISM spectroscopy across six JWST deep fields, we identify a large sample of 118 LRDs with high-S/N spectra, enabling a population-wide analysis of their UV-optical continuum and emission lines. We find clear correlations between rest-frame color ([0.3-0.9\,$μ$m]) and slopes: bluer LRDs have blue UV slopes ($β_{ν,\mathrm{UV}} \sim 0.3$) and red optical slopes, while redder LRDs exhibit redder UV slopes ($β_{ν,\mathrm{UV}} \sim 1.1$). The continuum shape shows a similar trend: redder LRDs display prominent Balmer breaks and curvature, while bluer LRDs follow power-law-like optical SEDs. From literature compilations, $\sim$60% of known broad-line AGNs satisfy our LRD criteria, and up to 90% of LRDs show broad Balmer lines. Emission-line diagnostics reveal a shift from high H$_α$/H$_β$ and low [OIII]$\lambda5007$/H$_β$ in redder LRDs to the opposite in bluer ones, along with stronger narrow-line equivalent widths, suggesting a transition from AGN- to host-dominated emission. We fit the spectra with a two-component model combining a gas-enshrouded black hole (BH) and a galaxy host. Redder LRDs require higher-luminosity, unreddened BHs and modestly reddened hosts; bluer LRDs require lower-luminosity, reddened BHs and dust-free galaxies. This framework reproduces the diversity in colors and spectral shape by varying BH luminosity, obscuration, and host-to-BH luminosity ratio.

preprint2026arXiv

Galaxies in the Epoch of Reionization Are All Bark and No Bite -- Plenty of Ionizing Photons, Low Escape Fractions

Early results from JWST suggest that epoch-of-reionization (EoR) galaxies produce copious ionizing photons, which, if they escape efficiently, could cause reionization to occur too early. We study this problem using \jwst\ imaging and prism spectroscopy for 412 galaxies at 4.5 < z < 9.0. We fit these data simultaneously with stellar-population and nebular-emission models that include a parameter for the fraction of ionizing photons that escape the galaxy, $f_\mathrm{esc}$. We find that the ionization production efficiency, $ξ_\mathrm{ion}$ = Q(H) / L(UV), increases with redshift and decreasing UV luminosity, but shows significant scatter, $σ( \log ξ_\mathrm{ion})$ = 0.3 dex. The inferred escape fractions averaged over the population are low, ranging from $\langle f_\mathrm{esc} \rangle$ = $2.6\pm 1.4$\% at 6 < z < 9 to $6.5\pm 2.2$\% at 4.5 < z < 6 with weak or no indication of evolution with redshift. This implies that in our models most of the ionizing photons need to be absorbed to account for the nebular emission. We compute the impact of our results on reionization, including the distributions for $ξ_\mathrm{ion}$ and $f_\mathrm{esc}$, and the evolution and uncertainty of the UV luminosity function. Considering galaxies brighter than M(UV) < -16 mag, we would produce an IGM hydrogen-ionized fraction of $x_e = 0.5$ at 5.3 < z < 5.8, possibly too late compared to constraints from from QSO sightlines. Including fainter galaxies, M(UV) < -14 mag, we obtain $x_e = 0.5$ at 6.0 < z < 8.1, fully consistent with QSO and CMB data. This implies that EoR galaxies produce plenty of ionizing photons, but these do not efficiently escape. This may be a result of high gas column densities combined with burstier star-formation histories, which limit the time massive stars are able to clear channels through the gas for ionizing photons to escape.

preprint2026arXiv

VENUS: Two Faint Little Red Dots Separated by $\sim70\,\mathrm{pc}$ Hidden in a Single Lensed Galaxy at $z\sim7$

We report the identification of a pair of faint little red dots (LRDs), dubbed Red Eyes, in a strongly-lensed galaxy at $z\sim7$ behind the PLCKG004.5-10.5 cluster, identified from the JWST Treasury program VENUS. Red Eyes are spatially resolved on the image plane with distinct colors, while the critical curve lies far north of Red Eyes, clearly requiring two different LRDs rather than a single LRD. Red Eyes is an extremely close pair of LRDs separated by $\sim70\,\mathrm{pc}$ in the source plane with a magnification of $μ\sim20$, which consistently explains another counter-image detected to the north-west. Red Eyes is hosted in a typical star-forming galaxy with $M_{\mathrm{UV,int}}\sim -19$, but its own UV emission is very faint ($M_{\mathrm{UV,int}} \gtrsim -16$). Moreover, Red Eyes does not reside at the galaxy center but lies at an offset position of approximately one effective radius $R_{\mathrm{e}}$ away from the galaxy center. If observed without lensing, Red Eyes would appear as a typical star-forming galaxy at $z\sim 7$ with $M_{\mathrm{UV}}\sim -19$, showing no apparent LRD signatures in either morphology or SED. These results suggest that multiple off-center LRDs, similar to Red Eyes, may be commonly hidden in a typical high-$z$ star-forming galaxy. In this case, various plausible scenarios may emerge, one of which is that intermediate-mass black holes (IMBHs) with $M_\mathrm{BH}\sim10^{4\text{--}6}\,M_\odot$ may form in star clusters on a stellar disk and contribute to the growth of the central supermassive black hole via mergers, with some IMBHs detectable as luminous LRDs in a sufficiently active and massive phase.

preprint2023arXiv

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

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

preprint2022arXiv

Bright z~9 Galaxies in Parallel: The Bright End of the Rest-UV Luminosity Function from HST Parallel Programs

The abundance of bright galaxies at z>8 can provide key constraints on models of galaxy formation and evolution, as the predicted abundance varies greatly when different physical prescriptions for gas cooling and star formation are implemented. We present the results of a search for bright z=9-10 galaxies selected from pure-parallel Hubble Space Telescope imaging programs. We include 132 fields observed as part of the Brightest of Reionizing Galaxies survey, the Hubble Infrared Pure Parallel Imaging Extragalactic Survey, and the WFC3 Infrared Spectroscopic Parallel survey. These observations cover a total of 620 sq. arcmin, about 70% of which is also covered with Spitzer Space Telescope infrared imaging. We identify thirteen candidate galaxies in the range 8.3<z<11 with 24.5 < m_H < 26.5 (-22.9 < M_UV < -21.2). This sample capitalizes on the uncorrelated nature of pure parallel observations to overcome cosmic variance and leverages a full multi-wavelength selection process to minimize contamination without sacrificing completeness. We perform detailed completeness and contamination analyses, and present measurements of the bright end of the UV luminosity function using a pseudo-binning technique. We find a number density consistent with results from Finkelstein et al. (2022) and other searches in HST parallel fields. These bright candidates likely reside in overdensities, potentially representing some of the earliest sites of cosmic reionization. These new candidates are excellent targets for follow-up with JWST, and four of them will be observed with the NIRSpec prism in Cycle 1.

preprint2022arXiv

Fitting AGN/galaxy X-ray-to-radio SEDs with CIGALE and improvement of the code

Modern and future surveys effectively provide a panchromatic view for large numbers of extragalactic objects. Consistently modeling these multiwavelength survey data is a critical but challenging task for extragalactic studies. The Code Investigating GALaxy Emission (CIGALE) is an efficient PYTHON code for spectral energy distribution (SED) fitting of galaxies and active galactic nuclei (AGNs). Recently, a major extension of CIGALE (named X-CIGALE) has been developed to account for AGN/galaxy X-ray emission and improve AGN modeling at UV-to-IR wavelengths. Here, we apply X-CIGALE to different samples, including COSMOS spectroscopic type 2 AGNs, CDF-S X-ray detected normal galaxies, SDSS quasars, and COSMOS radio objects. From these tests, we identify several weaknesses of X-CIGALE and improve the code accordingly. These improvements are mainly related to AGN intrinsic X-ray anisotropy, X-ray binary emission, AGN accretion-disk SED shape, and AGN radio emission. These updates improve the fit quality and allow new interpretation of the results, based on which we discuss physical implications. For example, we find that AGN intrinsic X-ray anisotropy is moderate, and can be modeled as $L_X(θ) \propto 1+\cos θ$, where $θ$ is the viewing angle measured from the AGN axis. We merge the new code into the major branch of CIGALE, and publicly release this new version as CIGALE v2022.0 on https://cigale.lam.fr

preprint2022arXiv

Identifying Galaxy Mergers in Simulated CEERS NIRCam Images using Random Forests

Identifying merging galaxies is an important - but difficult - step in galaxy evolution studies. We present random forest classifications of galaxy mergers from simulated JWST images based on various standard morphological parameters. We describe (a) constructing the simulated images from IllustrisTNG and the Santa Cruz SAM, and modifying them to mimic future CEERS observations as well as nearly noiseless observations, (b) measuring morphological parameters from these images, and (c) constructing and training the random forests using the merger history information for the simulated galaxies available from IllustrisTNG. The random forests correctly classify $\sim60\%$ of non-merging and merging galaxies across $0.5 < z < 4.0$. Rest-frame asymmetry parameters appear more important for lower redshift merger classifications, while rest-frame bulge and clump parameters appear more important for higher redshift classifications. Adjusting the classification probability threshold does not improve the performance of the forests. Finally, the shape and slope of the resulting merger fraction and merger rate derived from the random forest classifications match with theoretical Illustris predictions, but are underestimated by a factor of $\sim 0.5$.

preprint2022arXiv

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

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.

preprint2022arXiv

On the Stellar Populations of Galaxies at z=9-11: The Growth of Metals and Stellar Mass at Early Times

We present a detailed stellar population analysis of 11 bright ($H<26.6$) galaxies at $z=9-11$ (three spectroscopically confirmed) to constrain the chemical enrichment and growth of stellar mass of early galaxies. We use the flexible Bayesian spectral energy distribution (SED) fitting code Prospector with a range of star-formation histories (SFHs), a flexible dust attenuation law and a self-consistent modeling of emission lines. This approach allows us to assess how different priors affect our results, and how well we can break degeneracies between dust attenuation, stellar ages, metallicity and emission lines using data which probe only the rest-frame ultraviolet to optical wavelengths. We measure a median observed ultraviolet spectral slope $β=-1.87_{-0.43}^{+0.35}$ for relatively massive star-forming galaxies ($9<\log(M_{\star}/M_{\odot})<10$), consistent with no change from $z=4$ to $z=9-10$ at these stellar masses, implying rapid enrichment. Our SED-fitting results are consistent with a star-forming main sequence with sub-linear slope ($0.7\pm0.2$) and specific star-formation rates of $3-10~\mathrm{Gyr}^{-1}$. However, the stellar ages and SFHs are less well constrained. Using different SFH priors, we cannot distinguish between median mass-weighted ages of $\sim50-150$ Myr, which corresponds to 50\% formation redshifts of $z_{50}\sim10-12$ at $z\sim9$ and is of the order of the dynamical timescales of these systems. Importantly, the models with different SFH priors are able to fit the data equally well. We conclude that the current observational data cannot tightly constrain the mass-buildup timescales of these $z=9-11$ galaxies, with our results consistent with SFHs implying both a shallow and steep increase of the cosmic SFR density with time at $z>10$.

preprint2022arXiv

The Low-Redshift Lyman Continuum Survey I: New, Diverse Local Lyman-Continuum Emitters

The origins of Lyman continuum (LyC) photons responsible for the reionization of the universe are as of yet unknown and highly contested. Detecting LyC photons from the epoch of reionization is not possible due to absorption by the intergalactic medium, which has prompted the development of several indirect diagnostics to infer the rate at which galaxies contribute LyC photons to reionize the universe by studying lower-redshift analogs. We present the Low-redshift Lyman Continuum Survey (LzLCS) comprising measurements made with HST/COS for a z=0.2-0.4 sample of 66 galaxies. After careful processing of the FUV spectra, we obtain a total of 35 Lyman continuum emitters (LCEs) detected with 97.725% confidence, nearly tripling the number of known local LCEs. We estimate escape fractions from the detected LyC flux and upper limits on the undetected LyC flux, finding a range of LyC escape fractions up to 50%. Of the 35 LzLCS LCEs, 12 have LyC escape fractions greater than 5%, more than doubling the number of known local LCEs with cosmologically relevant LyC escape.

preprint2022arXiv

The Low-Redshift Lyman Continuum Survey II: New Insights into LyC Diagnostics

The Lyman continuum (LyC) cannot be observed at the epoch of reionization (z {\gtrsim} 6) due to intergalactic H I absorption. To identify Lyman continuum emitters (LCEs) and infer the fraction of escaping LyC, astronomers have developed various indirect diagnostics of LyC escape. Using measurements of the LyC from the Low-redshift Lyman Continuum Survey (LzLCS), we present the first statistical test of these diagnostics. While optical depth indicators based on Lyα, such as peak velocity separation and equivalent width, perform well, we also find that other diagnostics, such as the [O III]/[O II] flux ratio and star formation rate surface density, predict whether a galaxy is a LCE. The relationship between these galaxy properties and the fraction of escaping LyC flux suggests that LyC escape depends strongly on H I column density, ionization parameter, and stellar feedback. We find LCEs occupy a range of stellar masses, metallicities, star formation histories, and ionization parameters, which may indicate episodic and/or different physical causes of LyC escape.

preprint2022arXiv

Using Multiple Emission Line Ratios to Constrain the Slope of the Dust Attenuation Law

We explore the possibility and practical limitations of using a three-line approach to measure both the slope and normalization of the dust attenuation law in individual galaxies. To do this, we focus on a sample of eleven galaxies with existing ground-based Balmer H$α$ and H$β$ measurements from slit spectra, plus space-based grism constraints on Paschen-$β$. When accounting for observational uncertainties, we show that one galaxy has a well-constrained dust law slope and normalization in the range expected from theoretical arguments; this galaxy therefore provides an example of what may be possible in the future. However, most of the galaxies are best-fit by unusually steep or shallow slopes. We then explore whether additional astrophysical effects or observational biases could explain the elevated Paschen-$β$/H$α$ ratios driving these results. We find that galaxies with high Paschen-$β$/H$α$ ratios may be explained by slightly sub-unity covering fractions ($>$97%). Alternatively, differing slit losses for different lines can have a large impact on the results, emphasizing the importance of measuring all three lines with a consistent spectroscopic aperture. We conclude that while the three-line approach to constraining the shape of the dust attenuation law in individual galaxies is promising, deep observations and a consistent observational strategy will be required to minimize observational biases and to disentangle the astrophysically interesting effect of differing covering fractions. The James Webb Space Telescope will provide more sensitive measurements of Balmer and Paschen lines for galaxies at $z\approx0.3-2$, enabling uniform constraints on the optical-infrared dust attenuation law and its intrinsic variation.

preprint2021arXiv

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

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

preprint2021arXiv

CLEAR: Boosted Ly$α$ Transmission of the Intergalactic Medium in UV bright Galaxies

Reionization is an inhomogeneous process, thought to begin in small ionized bubbles of the intergalactic medium (IGM) around overdense regions of galaxies. Recent Lyman-alpha (Ly$α$) studies during the epoch of reionization show growing evidence that ionized bubbles formed earlier around brighter galaxies, suggesting higher IGM transmission of Ly$α$ from these galaxies. We investigate this problem using IR slitless spectroscopy from the Hubble Space Telescope (HST) Wide-Field Camera 3 (WFC3) G102 grism observations of 148 galaxies selected via photometric redshifts at $6.0<z<8.2$. These galaxies have spectra extracted from the CANDELS Ly$α$ Emission at Reionization (CLEAR) survey. We combine the CLEAR data for 275 galaxies with the Keck/DEIMOS+MOSFIRE dataset from the Texas Spectroscopic Search for Ly$α$ Emission at the End of Reionization Survey. We then constrain the Ly$α$ equivalent-width (EW) distribution at $6.0<z<8.2$, which is described by an exponential form, $dN/d\text{EW}\propto\text{exp(-EW)}/W_0$, with the characteristic $e$-folding scale width ($W_0$). We confirm a significant drop of the Ly$α$ strength (or $W_0$) at $z>6$. Furthermore, we compare the redshift evolution of $W_0$ between galaxies at different UV luminosities. The UV-bright ($M_{\text{UV}}<-21$, or $L_{\text{UV}}>L^{*}$) galaxies show weaker evolution with a decrease of 0.4 ($\pm$0.2) dex in $W_0$ at $z>6$ while UV-faint ($M_{\text{UV}}>-21$, or $L_{\text{UV}}<L^{*}$) galaxies exhibit a significant drop by a factor of 0.7-0.8 ($\pm0.2$) dex in $W_0$ from $z<6$ to $z>6$. Our results add to the accumulating evidence that UV-bright galaxies exhibit boosted Ly$α$ transmission in the IGM, suggesting that reionization completes sooner in regions proximate to galaxies of higher UV luminosity.

preprint2021arXiv

CLEAR: Emission Line Ratios at Cosmic High Noon

We use Hubble Space Telescope WFC3 G102 and G141 grism spectroscopy to measure rest-optical emission-line ratios of 533 galaxies at $z\sim1.5$ in the CANDELS Ly$α$ Emission at Reionization (CLEAR) survey. We compare $\frac{[OIII]}{Hβ}$ vs. $\frac{[SII]}{(Hα+[NII])}$ as an "unVO87" diagram for 461 galaxies and $\frac{[OIII]}{Hb}$ vs. $\frac{[NeIII]}{[OII]}$ as an "OHNO" diagram for 91 galaxies. The unVO87 diagram does not effectively separate active galactic nuclei (AGN) and $[NeV]$ sources from star-forming galaxies, indicating that the unVO87 properties of star-forming galaxies evolve with redshift and overlap with AGN emission-line signatures at $z>1$. The OHNO diagram effectively separates X-ray AGN and $[NeV]$-emitting galaxies from the rest of the population. We find that the $\frac{[OIII]}{Hβ}$ line ratios are significantly anti-correlated with stellar mass and significantly correlated with $\log(L_{Hβ})$, while $\frac{[SII]}{(Hα+[NII])}$ is significantly anti-correlated with $\log(L_{Hβ})$. Comparison with MAPPINGS~V photoionization models indicates that these trends are consistent with lower metallicity and higher ionization in low-mass and high-SFR galaxies. We do not find evidence for redshift evolution of the emission-line ratios outside of the correlations with mass and SFR.Our results suggest that the OHNO diagram of $\frac{[OIII]}{Hb}$ vs. $\frac{[NeIII]}{[OII]}$ will be a useful indicator of AGN content and gas conditions in very high-redshift galaxies to be observed by the James Webb Space Telescope.

preprint2021arXiv

The Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) Survey Design, Reductions, and Detections

We describe the survey design, calibration, commissioning, and emission-line detection algorithms for the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). The goal of HETDEX is to measure the redshifts of over a million Ly$α$ emitting galaxies between 1.88<z<3.52, in a 540 deg^2 area encompassing a co-moving volume of 10.9 Gpc^3. No pre-selection of targets is involved; instead the HETDEX measurements are accomplished via a spectroscopic survey using a suite of wide-field integral field units distributed over the focal plane of the telescope. This survey measures the Hubble expansion parameter and angular diameter distance, with a final expected accuracy of better than 1%. We detail the project's observational strategy, reduction pipeline, source detection, and catalog generation, and present initial results for science verification in the COSMOS, Extended Groth Strip, and GOODS-N fields. We demonstrate that our data reach the required specifications in throughput, astrometric accuracy, flux limit, and object detection, with the end products being a catalog of emission-line sources, their object classifications, and flux-calibrated spectra.

preprint2020arXiv

CLEAR II: Evidence for Early Formation of the Most Compact Quiescent Galaxies at High Redshift

The origin of the correlations between mass, morphology, quenched fraction, and formation history in galaxies is difficult to define, primarily due to the uncertainties in galaxy star-formation histories. Star-formation histories are better constrained for higher redshift galaxies, observed closer to their formation and quenching epochs. Here we use "non-parametric" star-formation histories and a nested sampling method to derive constraints on the formation and quenching timescales of quiescent galaxies at $0.7<z<2.5$. We model deep HST grism spectroscopy and photometry from the CLEAR (CANDELS Lyman$-α$ Emission at Reionization) survey. The galaxy formation redshifts, $z_{50}$ (defined as the point where they had formed 50\% of their stellar mass) range from $z_{50}\sim 2$ (shortly prior to the observed epoch) up to $z_{50} \simeq 5-8$. \editone{We find that early formation redshifts are correlated with high stellar-mass surface densities, $\log Σ_1 / (M_\odot\ \mathrm{kpc}^{-2}) >$10.25, where $Σ_1$ is the stellar mass within 1~pkpc (proper kpc). Quiescent galaxies with the highest stellar-mass surface density, $\logΣ_1 / (M_\odot\ \mathrm{kpc}^{-2}) > 10.25$, } show a \textit{minimum} formation redshift: all such objects in our sample have $z_{50} > 2.9$. Quiescent galaxies with lower surface density, $\log Σ_1 / (M_\odot\ \mathrm{kpc}^{-2}) = 9.5 - 10.25$, show a range of formation epochs ($z_{50} \simeq 1.5 - 8$), implying these galaxies experienced a range of formation and assembly histories. We argue that the surface density threshold $\logΣ_1/(M_\odot\ \mathrm{kpc}^{-2})>10.25$ uniquely identifies galaxies that formed in the first few Gyr after the Big Bang, and we discuss the implications this has for galaxy formation models.

preprint2020arXiv

CLEAR: The Gas-Phase Metallicity Gradients of Star-Forming Galaxies at 0.6 < z < 2.6

We report on the gas-phase metallicity gradients of a sample of 264 star-forming galaxies at 0.6 < z < 2.6, measured through deep near-infrared Hubble Space Telescope slitless spectroscopy. The observations include 12-orbit depth Hubble/WFC3 G102 grism spectra taken as a part of the CANDELS Lya Emission at Reionization (CLEAR) survey, and archival WFC3 G102+G141 grism spectra overlapping the CLEAR footprint. The majority of galaxies (84%) in this sample are consistent with a zero or slightly positive metallicity gradient across the full mass range probed (8.5 < log M_*/M_sun < 10.5). We measure the intrinsic population scatter of the metallicity gradients, and show that it increases with decreasing stellar mass---consistent with previous reports in the literature, but confirmed here with a much larger sample. To understand the physical mechanisms governing this scatter, we search for correlations between the observed gradient and various stellar population properties at fixed mass. However, we find no evidence for a correlation with the galaxy properties we consider---including star-formation rates, sizes, star-formation rate surface densities, and star-formation rates per gravitational potential energy. We use the observed weakness of these correlations to provide material constraints for predicted intrinsic correlations from theoretical models.

preprint2020arXiv

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

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

preprint2020arXiv

JWST/MIRI Simulated Imaging: Insights into Obscured Star-Formation and AGN for Distant Galaxies in Deep Surveys

The JWST MIRI instrument will revolutionize extragalactic astronomy with unprecedented sensitivity and angular resolution in mid-IR. Here, we assess the potential of MIRI photometry to constrain galaxy properties in the Cosmic Evolution Early Release Science (CEERS) survey. We derive estimated MIRI fluxes from the spectral energy distributions (SEDs) of real sources that fall in a planned MIRI pointing. We also obtain MIRI fluxes for hypothetical AGN-galaxy mixed models varying the AGN fractional contribution to the total IR luminosity ($\rm frac_{AGN}$). Based on these model fluxes, we simulate CEERS imaging (3.6-hour exposure) in 6 bands from F770W to F2100W using MIRISIM, and reduce these data using JWST PIPELINE. We perform PSF-matched photometry with TPHOT, and fit the source SEDs with X-CIGALE, simultaneously modeling photometric redshift and other physical properties. Adding the MIRI data, the accuracy of both redshift and $\rm frac_{AGN}$ is generally improved by factors of $\gtrsim 2$ for all sources at $z\lesssim 3$. Notably, for pure-galaxy inputs ($\rm frac_{AGN}=0$), the accuracy of $\rm frac_{AGN}$ is improved by $\sim 100$ times thanks to MIRI. The simulated CEERS MIRI data are slightly more sensitive to AGN detections than the deepest X-ray survey, based on the empirical $L_{\rm X}$-$L_{\rm 6μm}$ relation. Like X-ray observations, MIRI can also be used to constrain the AGN accretion power (accuracy $\approx 0.3$ dex). Our work demonstrates that MIRI will be able to place strong constraints on the mid-IR luminosities from star formation and AGN, and thereby facilitate studies of the galaxy/AGN co-evolution.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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

preprint2020arXiv

Selection of massive evolved galaxies at $3 \leq z \leq 4.5$ in the CANDELS fields

Using the CANDELS photometric catalogs for the HST/ACS and WFC3, we identified massive evolved galaxies at $3 < z < 4.5$, employing three different selection methods. We find the comoving number density of these objects to be $\sim 2 \times 10^{-5}$ and $8 \times 10^{-6}Mpc^{-3}$ after correction for completeness for two redshift bins centered at $z=3.4, 4.7$. We quantify a measure of how much confidence we should have for each candidate galaxy from different selections and what are the conservative error estimates propagated into our selection. Then we compare the evolution of the corresponding number densities and their stellar mass density with numerical simulations, semi-analytical models, and previous observational estimates, which shows slight tension at higher redshifts as the models tend to underestimate the number and mass densities. By estimating the average halo masses of the candidates ($M_h \approx 4.2, 1.9, 1.3 \times 10^{12} M_\odot$ for redshift bins centered at $z=3.4, 4.1, 4.7$), we find them to be consistent with halos that were efficient in turning baryons to stars and were relatively immune to the feedback effects and on the verge of transition into hot-mode accretion. This can suggest the relative cosmological starvation of the cold gas followed by an overconsumption phase in which the galaxy consumes the available cold gas rapidly as one of the possible drivers for the quenching of the massive evolved population at high redshift.

preprint2020arXiv

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

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

preprint2020arXiv

ZFIRE: Measuring Electron Density with [OII] as a function of environment at z = 1.62

The global star formation rates (SFR) of galaxies at fixed stellar masses increase with redshift and are known to vary with environment unto z~2. We explore here whether the changes in the star formation rates can also apply to the electron densities of the inter-stellar medium (ISM) by measuring the [OII] (3727A/3729A) ratio for cluster and field galaxies at z~2. We measure a median electron density of ne = 366+/-84 cm-3 for six galaxies (with 1-sigma scatter = 163 cm-3) in the UDS proton-cluster at z=1.62. We find that the median electron density of galaxies in the UDS photo-cluster environment is three times higher compared to the median electron density of field galaxies (ne = 113+/- 63 cm-3 and 1-sigma scatter = 79 cm-3) at comparable redshifts, stellar mass and SFR. However, we note that a sample of six photo-cluster galaxies is insufficient to reliably measure the electron density in the average porto-cluster environment at z~2. We conclude that the electron density increases with redshift in both cluster and field environments up to z~2 (ne = 30 +/- 1 cm-3 for z ~ 0 to ne =254+/- 76 cm-3 for z~1.5). We find tentative evidence (~2.6 sigma ) for a possible dependence of electron density on environment, but the results require confirmation with larger sample sizes.

preprint2019arXiv

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

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

preprint2019arXiv

HST Imaging of the Ionizing Radiation from a Star-forming Galaxy at z = 3.794

We report on the HST detection of the Lyman-continuum (LyC) radiation emitted by a galaxy at redshift z=3.794, dubbed Ion1 (Vanzella et al. 2012). The LyC from Ion1 is detected at rest-frame wavelength 820$\sim$890 Åwith HST WFC3/UVIS in the F410M band ($m_{410}=27.60\pm0.36$ magnitude (AB), peak SNR = 4.17 in a circular aperture with radius r = 0.12'') and at 700$\sim$830 Åwith the VLT/VIMOS in the U-band ($m_U = 27.84\pm0.19$ magnitude (AB), peak SNR = 6.7 with a r = 0.6'' aperture). A 20-hr VLT/VIMOS spectrum shows low- and high-ionization interstellar metal absorption lines, the P-Cygni profile of CIV and Ly$α$ in absorption. The latter spectral feature differs from what observed in known LyC emitters, which show strong Ly$α$ emission. An HST far-UV color map reveals that the LyC emission escapes from a region of the galaxy that is bluer than the rest, presumably because of lower dust obscuration. The F410M image shows that the centroid of the LyC emission is offset from the centroid of the non-ionizing UV emission by 0.12''$\pm$0.03'', corresponding to 0.85$\pm$0.21 kpc (physical), and that its morphology is likely moderately resolved. These morphological characteristics favor a scenario where the LyC photons produced by massive stars escape from low HI column-density "cavities" in the ISM, possibly carved by stellar winds and/or supernova. We also collect the VIMOS U-band images of a sample of 107 Lyman-break galaxies with spectroscopic redshifts at $3.40<z<3.95$, i.e. sampling the LyC, and stack them with inverse-variance weights. No LyC emission is detected in the stacked image, resulting in a 32.5 magnitude (AB) flux limit (1$σ$) and an upper limit of absolute LyC escape fraction $f_{esc}^{abs} < 0.63\%$. LyC emitters like Ion1 are very likely at the bright-end of the LyC luminosity function.

preprint2019arXiv

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

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

preprint2016arXiv

A New Star-Formation Rate Calibration from Polycyclic Aromatic Hydrocarbon Emission Features and Application to High Redshift Galaxies

We calibrate the integrated luminosity from the polycyclic aromatic hydrocarbon (PAH) features at 6.2\micron, 7.7\micron\ and 11.3\micron\ in galaxies as a measure of the star-formation rate (SFR). These features are strong (containing as much as 5-10\% of the total infrared luminosity) and suffer minimal extinction. Our calibration uses \spitzer\ Infrared Spectrograph (IRS) measurements of 105 galaxies at $0 < z < 0.4$, infrared (IR) luminosities of $10^9 - 10^{12} \lsol$, combined with other well-calibrated SFR indicators. The PAH luminosity correlates linearly with the SFR as measured by the extinction-corrected \ha\ luminosity over the range of luminosities in our calibration sample. The scatter is 0.14 dex comparable to that between SFRs derived from the \paa\ and extinction-corrected \ha\ emission lines, implying the PAH features may be as accurate a SFR indicator as hydrogen recombination lines. The PAH SFR relation depends on gas-phase metallicity, for which we supply an empirical correction for galaxies with $0.2 < \mathrm{Z} \lsim 0.7$~\zsol. We present a case study in advance of the \textit{James Webb Space Telescope} (\jwst), which will be capable of measuring SFRs from PAHs in distant galaxies at the peak of the SFR density in the universe ($z\sim2$) with SFRs as low as $\sim$~10~\sfrunits. We use \spitzer/IRS observations of the PAH features and \paa\ emission plus \ha\ measurements in lensed star-forming galaxies at $1 < z < 3$ to demonstrate the ability of the PAHs to derive accurate SFRs. We also demonstrate that because the PAH features dominate the mid-IR fluxes, broad-band mid-IR photometric measurements from \jwst\ will trace both the SFR and provide a way to exclude galaxies dominated by an AGN.

preprint2016arXiv

Beyond Spheroids and Discs: Classifications of CANDELS Galaxy Structure at 1.4 < z < 2 via Principal Component Analysis

Important but rare and subtle processes driving galaxy morphology and star-formation may be missed by traditional spiral, elliptical, irregular or Sérsic bulge/disk classifications. To overcome this limitation, we use a principal component analysis of non-parametric morphological indicators (concentration, asymmetry, Gini coefficient, $M_{20}$, multi-mode, intensity and deviation) measured at rest-frame $B$-band (corresponding to HST/WFC3 F125W at 1.4 $< z <$ 2) to trace the natural distribution of massive ($>10^{10} M_{\odot}$) galaxy morphologies. Principal component analysis (PCA) quantifies the correlations between these morphological indicators and determines the relative importance of each. The first three principal components (PCs) capture $\sim$75 per cent of the variance inherent to our sample. We interpret the first principal component (PC) as bulge strength, the second PC as dominated by concentration and the third PC as dominated by asymmetry. Both PC1 and PC2 correlate with the visual appearance of a central bulge and predict galaxy quiescence. PC1 is a better predictor of quenching than stellar mass, as as good as other structural indicators (Sérsic-n or compactness). We divide the PCA results into groups using an agglomerative hierarchical clustering method. Unlike Sérsic, this classification scheme separates compact galaxies from larger, smooth proto-elliptical systems, and star-forming disk-dominated clumpy galaxies from star-forming bulge-dominated asymmetric galaxies. Distinguishing between these galaxy structural types in a quantitative manner is an important step towards understanding the connections between morphology, galaxy assembly and star-formation.

preprint2016arXiv

Breaking the Curve with CANDELS: A Bayesian Approach to Reveal the Non-Universality of the Dust-Attenuation Law at High Redshift

Dust attenuation affects nearly all observational aspects of galaxy evolution, yet very little is known about the form of the dust-attenuation law in the distant Universe. Here, we model the spectral energy distributions (SEDs) of galaxies at z = 1.5--3 from CANDELS with rest-frame UV to near-IR imaging under different assumptions about the dust law, and compare the amount of inferred attenuated light with the observed infrared (IR) luminosities. Some individual galaxies show strong Bayesian evidence in preference of one dust law over another, and this preference agrees with their observed location on the plane of infrared excess (IRX, $L_{\text{TIR}}/L_{\text{UV}}$) and UV slope ($β$). We generalize the shape of the dust law with an empirical model, $A_{λ,δ}=E(B-V)\ k_λ (λ/λ_V)^δ$ where $k_λ$ is the dust law of Calzetti et al. (2000), and show that there exists a correlation between the color excess ${E(B-V)}$ and tilt $δ$ with ${δ=(0.62\pm0.05)\log(E(B-V))}$+ ${(0.26~\pm~0.02)}$. Galaxies with high color excess have a shallower, starburst-like law, and those with low color excess have a steeper, SMC-like law. Surprisingly, the galaxies in our sample show no correlation between the shape of the dust law and stellar mass, star-formation rate, or $β$. The change in the dust law with color excess is consistent with a model where attenuation is caused by by scattering, a mixed star-dust geometry, and/or trends with stellar population age, metallicity, and dust grain size. This rest-frame UV-to-near-IR method shows potential to constrain the dust law at even higher ($z>3$) redshifts.

preprint2016arXiv

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

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.

preprint2016arXiv

The Evolution of the Galaxy Stellar Mass Function at z= 4-8: A Steepening Low-mass-end Slope with Increasing Redshift

We present galaxy stellar mass functions (GSMFs) at $z=$ 4-8 from a rest-frame ultraviolet (UV) selected sample of $\sim$4500 galaxies, found via photometric redshifts over an area of $\sim$280 arcmin$^2$ in the CANDELS/GOODS fields and the Hubble Ultra Deep Field. The deepest Spitzer/IRAC data yet-to-date and the relatively large volume allow us to place a better constraint at both the low- and high-mass ends of the GSMFs compared to previous space-based studies from pre-CANDELS observations. Supplemented by a stacking analysis, we find a linear correlation between the rest-frame UV absolute magnitude at 1500 Å ($M_{\rm UV}$) and logarithmic stellar mass ($\log M_*$) that holds for galaxies with $\log(M_*/M_{\odot}) \lesssim 10$. We use simulations to validate our method of measuring the slope of the $\log M_*$-$M_{\rm UV}$ relation, finding that the bias is minimized with a hybrid technique combining photometry of individual bright galaxies with stacked photometry for faint galaxies. The resultant measured slopes do not significantly evolve over $z=$ 4-8, while the normalization of the trend exhibits a weak evolution toward lower masses at higher redshift. We combine the $\log M_*$-$M_{\rm UV}$ distribution with observed rest-frame UV luminosity functions at each redshift to derive the GSMFs, finding that the low-mass-end slope becomes steeper with increasing redshift from $α=-1.55^{+0.08}_{-0.07}$ at $z=4$ to $α=-2.25^{+0.72}_{-0.35}$ at $z=8$. The inferred stellar mass density, when integrated over $M_*=10^8$-$10^{13} M_{\odot}$, increases by a factor of $10^{+30}_{-2}$ between $z=7$ and $z=4$ and is in good agreement with the time integral of the cosmic star formation rate density.

preprint2016arXiv

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

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$.

preprint2016arXiv

The Spitzer-HETDEX Exploratory Large-Area Survey

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

preprint2016arXiv

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

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

preprint2016arXiv

ZFOURGE catalogue of AGN candidates: an enhancement of 160μm-derived star-formation rates in active galaxies to $z$ = 3.2

We investigate active galactic nuclei (AGN) candidates within the FourStar Galaxy Evolution Survey (ZFOURGE) to determine the impact they have on star-formation in their host galaxies. We first identify a population of radio, X-ray, and infrared-selected AGN by cross-matching the deep $K_{s}$-band imaging of ZFOURGE with overlapping multi-wavelength data. From this, we construct a mass-complete (log(M$_{*}$/M$_{\odot}$) $\ge$ 9.75), AGN luminosity limited sample of 235 AGN hosts over z = 0.2 - 3.2. We compare the rest-frame U - V versus V - J (UVJ) colours and specific star-formation rates (sSFRs) of the AGN hosts to a mass-matched control sample of inactive (non-AGN) galaxies. UVJ diagnostics reveal AGN tend to be hosted in a lower fraction of quiescent galaxies and a higher fraction of dusty galaxies than the control sample. Using 160μm Herschel PACS data, we find the mean specific star-formation rate of AGN hosts to be elevated by 0.34$\pm$0.07 dex with respect to the control sample across all redshifts. This offset is primarily driven by infrared-selected AGN, where the mean sSFR is found to be elevated by as much as a factor of ~5. The remaining population, comprised predominantly of X-ray AGN hosts, is found mostly consistent with inactive galaxies, exhibiting only a marginal elevation. We discuss scenarios that may explain these findings and postulate that AGN are less likely to be a dominant mechanism for moderating galaxy growth via quenching than has previously been suggested.

preprint2015arXiv

An Increasing Stellar Baryon Fraction in Bright Galaxies at High Redshift

Recent observations have shown that the characteristic luminosity of the rest-frame ultraviolet (UV) luminosity function does not significantly evolve at 4 < z < 7 and is approximately M*_UV ~ -21. We investigate this apparent non-evolution by examining a sample of 178 bright, M_UV < -21 galaxies at z=4 to 7, analyzing their stellar populations and host halo masses. Including deep Spitzer/IRAC imaging to constrain the rest-frame optical light, we find that M*_UV galaxies at z=4-7 have similar stellar masses of log(M/Msol)=9.6-9.9 and are thus relatively massive for these high redshifts. However, bright galaxies at z=4-7 are less massive and have younger inferred ages than similarly bright galaxies at z=2-3, even though the two populations have similar star formation rates and levels of dust attenuation. We match the abundances of these bright z=4-7 galaxies to halo mass functions from the Bolshoi Lambda-CDM simulation to estimate the halo masses. We find that the typical halo masses in ~M*_UV galaxies decrease from log(M_h/Msol)=11.9 at z=4 to log(M_h/Msol)=11.4 at z=7. Thus, although we are studying galaxies at a similar mass across multiple redshifts, these galaxies live in lower mass halos at higher redshift. The stellar baryon fraction in units of the cosmic mean Omega_b/Omega_m rises from 5.1% at z=4 to 11.7% at z=7; this evolution is significant at the ~3-sigma level. This rise does not agree with simple expectations of how galaxies grow, and implies that some effect, perhaps a diminishing efficiency of feedback, is allowing a higher fraction of available baryons to be converted into stars at high redshifts.

preprint2015arXiv

Dusty galaxies and the degeneracy between their dust distributions and the attenuation formula

Do spatial distributions of dust grains in galaxies have typical forms, as do spatial distributions of stars? We investigate whether or not the distributions resemble uniform foreground screens, as commonly assumed by the high-redshift galaxy community. We use rest-frame infrared, ultraviolet, and H$α$ line luminosities of dust-poor and dusty galaxies at z ~ 0 and z ~ 1 to compare measured H$α$ escape fractions with those predicted by the Calzetti attenuation formula. The predictions, based on UV escape fractions, overestimate the measured H$α$ escape fractions for all samples. The interpretation of this result for dust-poor z ~ 0 galaxies is that regions with ionizing stars have more dust than regions with nonionizing UV-emitting stars. Dust distributions for these galaxies are nonuniform. The interpretation of the overestimates for dusty galaxies at both redshifts is less clear. If the attenuation formula is inapplicable to these galaxies, perhaps the disagreements are unphysical; perhaps dust distributions in these galaxies are uniform. If the attenuation formula does apply, then dusty galaxies have nonuniform dust distributions; the distributions are more uniform than they are in dust-poor galaxies. A broad range of H$α$ escape fractions at a given UV escape fraction for z ~ 1 dusty galaxies, if real, indicates diverse dust morphologies and the implausibility of the screen assumption.

preprint2015arXiv

Quenching and Morphological Transformation in Semi-Analytic Models and CANDELS

We examine the spheroid growth and star formation quenching experienced by galaxies from z~3 to the present by studying the evolution with redshift of the quiescent and spheroid-dominated fractions of galaxies from the CANDELS and GAMA surveys. We compare the observed fractions with predictions from a semi-analytic model which includes prescriptions for bulge growth and AGN feedback due to mergers and disk instabilities. We facilitate direct morphological comparison by converting our model bulge-to-total stellar mass ratios to Sersic indices. We then subdivide our population into the four quadrants of the sSFR-Sersic index plane and study the buildup of each of these subpopulations. We find that the fraction of star forming disks declines steadily, while the fraction of quiescent spheroids builds up over cosmic time. The fractions of star forming spheroids and quiescent disks are both non-negligible, and stay nearly constant over the period we have studied, at about 10% and 15-20% respectively. Our model is qualitatively successful at reproducing the evolution of the two "main" populations (star forming disk-dominated galaxies and quiescent spheroid-dominated galaxies), and approximately reproduces the relative fractions of all four types, but predicts a stronger decline in star forming spheroids, and increase in quiescent disks, than seen in the observations. A model with an additional channel for bulge growth via disk instabilities agrees better overall with the observations than a model in which bulges may grow only through mergers. We study evolutionary tracks of some individual galaxies as they experience morphological transformation and quenching, and examine the importance of different physical drivers of this transformation (major and minor mergers and disk instabilities). We find that complex histories with multiple transformative events are the norm.

preprint2015arXiv

Satellite Quenching and Galactic Conformity at 0.3 < z < 2.5

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

preprint2015arXiv

Star-forming blue ETGs in two newly discovered galaxy overdensities in the HUDF at z=1.84 and 1.9: unveiling the progenitors of passive ETGs in cluster cores

We present the discovery of two galaxy overdensities in the Hubble Space Telescope UDF: a proto-cluster, HUDFJ0332.4-2746.6 at $z = 1.84 \pm 0.01$, and a group, HUDFJ0332.5-2747.3 at $z =1.90 \pm 0.01$. Assuming viralization, the velocity dispersion of HUDFJ0332.4-2746.6 implies a mass of $M_{200}= (2.2 \pm 1.8) \times 10^{14} M_{\odot}$, consistent with the lack of extended X-ray emission. Neither overdensity shows evidence of a red sequence. About $50\%$ of their members show interactions and/or disturbed morphologies, which are signatures of merger remnants or disk instability. Most of their ETGs have blue colors and show recent star-formation. These observations reveal for the first time large fractions of spectroscopically confirmed star-forming blue ETGs in proto-clusters at $z\approx 2$. These star-forming ETGs are most likely among the progenitors of the quiescent population in clusters at more recent epochs. Their mass-size relation is consistent with that of passive ETGs in clusters at $z\sim0.7-1.5$. If these galaxies are the progenitors of cluster ETGs at these lower redshifts, their size would evolve according to a similar mass-size relation. It is noteworthy that quiescent ETGs in clusters at $z=1.8-2$ also do not show any significant size evolution over this redshift range, contrary to field ETGs. The ETG fraction is $\lesssim 50\%$, compared to the typical quiescent ETG fraction of $\approx 80\%$ in cluster cores at $z< 1$. The fraction, masses, and colors of the newly discovered ETGs imply that other cluster ETGs will be formed/accreted at a later time.

preprint2015arXiv

The Differential Size Growth of Field and Cluster Galaxies at z=2.1 Using the ZFOURGE Survey

There is ongoing debate regarding the extent that environment affects galaxy size growth beyond z>1. To investigate the differences in star-forming and quiescent galaxy properties as a function of environment at z=2.1, we create a mass-complete sample of 59 cluster galaxies Spitler et al. (2012) and 478 field galaxies with log(M)>9 using photometric redshifts from the ZFOURGE survey. We compare the mass-size relation of field and cluster galaxies using measured galaxy semi-major axis half-light radii ($r_{1/2,maj}$) from CANDELS HST/F160W imaging. We find consistent mass normalized (log(M)=10.7) sizes for quiescent field galaxies ($r_{1/2,maj}=1.81\pm0.29$ kpc) and quiescent cluster galaxies ($r_{1/2,maj}=2.17\pm0.63$ kpc). The mass normalized size of star-forming cluster galaxies ($r_{1/2,maj}=4.00\pm0.26$ kpc ) is 12% larger (KS test $2.1σ$) than star-forming field galaxies ($r_{1/2,maj}=3.57\pm0.10$ kpc). From the mass-color relation we find that quiescent field galaxies with 9.7<log(M)<10.4 are slightly redder (KS test $3.6σ$) than quiescent cluster galaxies, while cluster and field quiescent galaxies with log(M)>10.4 have consistent colors. We find that star-forming cluster galaxies are on average 20% redder than star-forming field galaxies at all masses. Furthermore, we stack galaxy images to measure average radial color profiles as a function of mass. Negative color gradients are only present for massive star-forming field and cluster galaxies with log(M)>10.4, the remaining galaxy masses and types have flat profiles. Our results suggest given the observed differences in size and color of star-forming field and cluster galaxies, that the environment has begun to influence/accelerate their evolution. However, the lack of differences between field and cluster quiescent galaxies indicates that the environment has not begun to significantly influence their evolution at z~2.

preprint2015arXiv

The Evolution of the Galaxy Rest-Frame Ultraviolet Luminosity Function Over the First Two Billion Years

We present a robust measurement and analysis of the rest-frame ultraviolet (UV) luminosity function at z=4-8. We use deep Hubble Space Telescope imaging over the CANDELS/GOODS fields, the Hubble Ultra Deep Field and the Year 1 Hubble Frontier Field deep parallel observations. These surveys provides an effective volume of 0.6-1.2 x 10^6 Mpc^3 over this epoch, allowing us to perform a robust search for faint (M_UV=-18) and bright (M_UV < -21) galaxies. We select candidate galaxies using a well-tested photometric redshift technique with careful screening of contaminants, finding a sample of 7446 galaxies at 3.5<z<8.5, with >1000 galaxies at z~6-8. We measure the luminosity function using a Markov Chain Monte Carlo analysis to measure robust uncertainties. At the faint end our results agree with previous studies, yet we find a higher abundance of UV-bright galaxies at z>6, with M* ~ -21 at z>5, different than that inferred based on previous trends at lower redshift. At z=8, a single power-law provides an equally good fit to the UV luminosity function, while at z=6 and 7, an exponential cutoff at the bright-end is moderately preferred. We compare to semi-analytical models, and find that the lack of evolution in M* is consistent with models where the impact of dust attenuation on the bright-end of the luminosity function decreases at higher redshift. We measure the evolution of the cosmic star-formation rate density, correcting for dust attenuation, and find that it declines as (1+z)^(-4.3 +/- 0.5) at z>4, consistent with observations at z>9. Our observations are consistent with a reionization history that starts at z>10, completes at z>6, and reaches a midpoint (x_HII = 0.5) at 6.7<z<9.4. Finally, our observations predict that the abundance of bright z=9 galaxies is likely higher than previous constraints, though consistent with recent estimates of bright z~10 galaxies. [abridged]

preprint2015arXiv

The Herschel view of the dominant mode of galaxy growth from z=4 to the present day

We present an analysis of the deepest Herschel images in four major extragalactic fields GOODS-North, GOODS-South, UDS and COSMOS obtained within the GOODS-Herschel and CANDELS-Herschel key programs. The picture provided by 10497 individual far-infrared detections is supplemented by the stacking analysis of a mass-complete sample of 62361 star-forming galaxies from the CANDELS-HST H band-selected catalogs and from two deep ground-based Ks band-selected catalogs in the GOODS-North and the COSMOS-wide fields, in order to obtain one of the most accurate and unbiased understanding to date of the stellar mass growth over the cosmic history. We show, for the first time, that stacking also provides a powerful tool to determine the dispersion of a physical correlation and describe our method called "scatter stacking" that may be easily generalized to other experiments. We demonstrate that galaxies of all masses from z=4 to 0 follow a universal scaling law, the so-called main sequence of star-forming galaxies. We find a universal close-to-linear slope of the logSFR-logM* relation with evidence for a flattening of the main sequence at high masses (log(M*/Msun) > 10.5) that becomes less prominent with increasing redshift and almost vanishes by z~2. This flattening may be due to the parallel stellar growth of quiescent bulges in star-forming galaxies. Within the main sequence, we measure a non varying SFR dispersion of 0.3 dex. The specific SFR (sSFR=SFR/M*) of star-forming galaxies is found to continuously increase from z=0 to 4. Finally we discuss the implications of our findings on the cosmic SFR history and show that more than 2/3 of present-day stars must have formed in a regime dominated by the main sequence mode. As a consequence we conclude that, although omnipresent in the distant Universe, galaxy mergers had little impact in shaping the global star formation history over the last 12.5 Gyr.

preprint2015arXiv

The Relation Between SFR and Stellar Mass for Galaxies at 3.5 $\le z\le$ 6.5 in CANDELS

Distant star-forming galaxies show a correlation between their star formation rates (SFR) and stellar masses, and this has deep implications for galaxy formation. Here, we present a study on the evolution of the slope and scatter of the SFR-stellar mass relation for galaxies at $3.5\leq z\leq 6.5$ using multi-wavelength photometry in GOODS-S from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) and Spitzer Extended Deep Survey. We describe an updated, Bayesian spectral-energy distribution fitting method that incorporates effects of nebular line emission, star formation histories that are constant or rising with time, and different dust attenuation prescriptions (starburst and Small Magellanic Cloud). From $z$=6.5 to $z$=3.5 star-forming galaxies in CANDELS follow a nearly unevolving correlation between stellar mass and SFR that follows SFR $\sim$ $M_\star^a$ with $a = 0.54 \pm 0.16$ at $z\sim 6$ and $0.70 \pm 0.21$ at $z\sim 4$. This evolution requires a star formation history that increases with decreasing redshift (on average, the SFRs of individual galaxies rise with time). The observed scatter in the SFR-stellar mass relation is tight, $σ(\log \mathrm{SFR}/\mathrm{M}_\odot$ yr$^{-1})< 0.3\ - $ 0.4 dex, for galaxies with $\log M_\star/\mathrm{M}_\odot > 9$ dex. Assuming that the SFR is tied to the net gas inflow rate (SFR $\sim$ $\dot{M}_\mathrm{gas}$), then the scatter in the gas inflow rate is also smaller than 0.3$-$0.4 dex for star-forming galaxies in these stellar mass and redshift ranges, at least when averaged over the timescale of star formation. We further show that the implied star formation history of objects selected on the basis of their co-moving number densities is consistent with the evolution in the SFR-stellar mass relation.

preprint2015arXiv

The SFR-M* Relation and Empirical Star-Formation Histories from ZFOURGE at 0.5 < z < 4

We explore star-formation histories (SFHs) of galaxies based on the evolution of the star-formation rate stellar mass relation (SFR-M*). Using data from the FourStar Galaxy Evolution Survey (ZFOURGE) in combination with far-IR imaging from the Spitzer and Herschel observatories we measure the SFR-M* relation at 0.5 < z < 4. Similar to recent works we find that the average infrared SEDs of galaxies are roughly consistent with a single infrared template across a broad range of redshifts and stellar masses, with evidence for only weak deviations. We find that the SFR-M* relation is not consistent with a single power-law of the form SFR ~ M*^a at any redshift; it has a power-law slope of a~1 at low masses, and becomes shallower above a turnover mass (M_0) that ranges from 10^9.5 - 10^10.8 Msol, with evidence that M_0 increases with redshift. We compare our measurements to results from state-of-the-art cosmological simulations, and find general agreement in the slope of the SFR-M* relation albeit with systematic offsets. We use the evolving SFR-M* sequence to generate SFHs, finding that typical SFRs of individual galaxies rise at early times and decline after reaching a peak. This peak occurs earlier for more massive galaxies. We integrate these SFHs to generate mass-growth histories and compare to the implied mass-growth from the evolution of the stellar mass function. We find that these two estimates are in broad qualitative agreement, but that there is room for improvement at a more detailed level. At early times the SFHs suggest mass-growth rates that are as much as 10x higher than inferred from the stellar mass function. However, at later times the SFHs under-predict the inferred evolution, as is expected in the case of additional growth due to mergers.

preprint2015arXiv

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

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

preprint2014arXiv

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

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

preprint2014arXiv

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

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

preprint2014arXiv

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

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

preprint2014arXiv

Kpc-scale Properties of Emission-line Galaxies

We perform a detailed study of the resolved properties of emission-line galaxies at kpc-scale to investigate how small-scale and global properties of galaxies are related. 119 galaxies with high-resolution Keck/DEIMOS spectra are selected to cover a wide range in morphologies over the redshift range 0.2<z<1.3. Using the HST/ACS and HST/WFC3 imaging data taken as a part of the CANDELS project, for each galaxy we perform SED fitting per resolution element, producing resolved rest-frame U-V color, stellar mass, star formation rate, age and extinction maps. We develop a technique to identify blue and red "regions" within individual galaxies, using their rest-frame color maps. As expected, for any given galaxy, the red regions are found to have higher stellar mass surface densities and older ages compared to the blue regions. Furthermore, we quantify the spatial distribution of red and blue regions with respect to both redshift and stellar mass, finding that the stronger concentration of red regions toward the centers of galaxies is not a significant function of either redshift or stellar mass. We find that the "main sequence" of star forming galaxies exists among both red and blue regions inside galaxies, with the median of blue regions forming a tighter relation with a slope of 1.1+/-0.1 and a scatter of ~0.2 dex compared to red regions with a slope of 1.3+/-0.1 and a scatter of ~0.6 dex. The blue regions show higher specific Star Formation Rates (sSFR) than their red counterparts with the sSFR decreasing since z~1, driver primarily by the stellar mass surface densities rather than the SFRs at a giver resolution element.

preprint2014arXiv

Mid-Infrared Selected Quasars I: Virial Black Hole Mass and Eddington Ratios

We provide a catalog of 391 mid-infrared-selected (MIR, 24$μ$m) broad-emission-line (BEL, type 1) quasars in the 22 deg$^2$ SWIRE Lockman Hole field. This quasar sample is selected in the MIR from Spitzer MIPS with $S_{\rm 24} > 400μ$Jy, jointly with an optical magnitude limit of r (AB) $<$ 22.5 for broad line identification. The catalog is based on MMT and SDSS spectroscopy to select BEL quasars, extends the SDSS coverage to fainter magnitudes and lower redshifts, and recovers a more complete quasar population. The MIR-selected quasar sample peaks at $z\sim$1.4, and recovers a significant and constant (20\%) fraction of extended objects with SDSS photometry across magnitudes, which was not included in the SDSS quasar survey dominated by point sources. This sample also recovers a significant population of $z < 3$ quasars at $i > 19.1$. We then investigate the continuum luminosity and line profiles of these MIR quasars, and estimate their virial black hole masses and the Eddington ratios. The SMBH mass shows evidence of downsizing, though the Eddington ratios remain constant at $1 < z < 4$. Compared to point sources in the same redshift range, extended sources at $z < 1$ show systematically lower Eddington ratios. The catalog and spectra are publicly available online.

preprint2014arXiv

Rapid Decline of Lyman-alpha Emission Toward the Reionization Era

The observed deficit of strongly Lyman-alpha emitting galaxies at z>6.5 is attributed to either increasing neutral hydrogen in the intergalactic medium (IGM) and/or to the evolving galaxy properties. To investigate this, we have performed very deep near-IR spectroscopy of z>7 galaxies using MOSFIRE on the Keck-I Telescope. We measure the Lyman-alpha fraction at z~8 (combined photometric redshift peak at z=7.7) using two methods. First, we derived NLyα/Ntot directly using extensive simulations to correct for incompleteness. Second, we used a Bayesian formalism (introduced by Treu et al. 2012) that compares the z>7 galaxy spectra to models of the Lyman-alpha equivalent width (WLyα) distribution at z~6. We explored two simple evolutionary scenarios: smooth evolution where Lyman-alpha is attenuated in all galaxies by a constant factor (perhaps owing to processes from galaxy evolution or a slowly increasing IGM opacity), and patchy evolution where Lyman-alpha is blocked in some fraction of galaxies (perhaps due to the IGM being opaque along only some fraction of sightlines). The Bayesian formalism places stronger constraints compared with the direct method. Combining our data with that in the literature we find that at z~8 the Lyman-alpha fraction has dropped by a factor >3(84% confidence interval) using both the smooth and patchy scenarios compared to the z~6 values. Furthermore, we find a tentative evidence that the data favor the patchy scenario over smooth (with "positive" Bayesian evidence), extending trends observed at z~7 to higher redshift. If this decrease is a result of reionization as predicted by theory, then our data imply the volume averaged neutral hydrogen fraction in the IGM to be >0.3 suggesting that the reionization of the universe is in progress at z~8.

preprint2014arXiv

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

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

preprint2014arXiv

ZFOURGE/CANDELS: On the Evolution of M* Galaxy Progenitors from z=3 to 0.5

Galaxies with stellar masses near M* contain the majority of stellar mass in the universe, and are therefore of special interest in the study of galaxy evolution. The Milky Way (MW) and Andromeda (M31) have present day stellar masses near M*, at 5x10^10 Msol (MW-mass) and 10^11 Msol (M31-mass). We study the typical progenitors of these galaxies using ZFOURGE, a deep medium-band near-IR imaging survey, which is sensitive to the progenitors of these galaxies out to z~3. We use abundance-matching techniques to identify the main progenitors of these galaxies at higher redshifts. We measure the evolution in the stellar mass, rest-frame colors, morphologies, far-IR luminosities, and star-formation rates combining our deep multiwavelength imaging with near-IR HST imaging from CANDELS, and far-IR imaging from GOODS-H and CANDELS-H. The typical MW-mass and M31-mass progenitors passed through the same evolution stages, evolving from blue, star-forming disk galaxies at the earliest stages, to redder dust-obscured IR-luminous galaxies in intermediate stages, and to red, more quiescent galaxies at their latest stages. The progenitors of the MW-mass galaxies reached each evolutionary stage at later times (lower redshifts) and with stellar masses that are a factor of 2-3 lower than the progenitors of the M31-mass galaxies. The process driving this evolution, including the suppression of star-formation in present-day M* galaxies requires an evolving stellar-mass/halo-mass ratio and/or evolving halo-mass threshold for quiescent galaxies. The effective size and star-formation rates imply that the baryonic cold-gas fractions drop as galaxies evolve from high redshift to z~0 and are strongly anticorrelated with an increase in the Sérsic index. Therefore, the growth of galaxy bulges in M* galaxies corresponds to a rapid decline in the galaxy gas fractions and/or a decrease in the star-formation efficiency.

preprint2013arXiv

A Critical Assessment of Photometric Redshift Methods: A CANDELS Investigation

We present results from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) photometric redshift methods investigation. In this investigation, the results from eleven participants, each using a different combination of photometric redshift code, template spectral energy distributions (SEDs) and priors, are used to examine the properties of photometric redshifts applied to deep fields with broad-band multi-wavelength coverage. The photometry used includes U-band through mid-infrared filters and was derived using the TFIT method. Comparing the results, we find that there is no particular code or set of template SEDs that results in significantly better photometric redshifts compared to others. However, we find codes producing the lowest scatter and outlier fraction utilize a training sample to optimize photometric redshifts by adding zero-point offsets, template adjusting or adding extra smoothing errors. These results therefore stress the importance of the training procedure. We find a strong dependence of the photometric redshift accuracy on the signal-to-noise ratio of the photometry. On the other hand, we find a weak dependence of the photometric redshift scatter with redshift and galaxy color. We find that most photometric redshift codes quote redshift errors (e.g., 68% confidence intervals) that are too small compared to that expected from the spectroscopic control sample. We find that all codes show a statistically significant bias in the photometric redshifts. However, the bias is in all cases smaller than the scatter, the latter therefore dominates the errors. Finally, we find that combining results from multiple codes significantly decreases the photometric redshift scatter and outlier fraction. We discuss different ways of combining data to produce accurate photometric redshifts and error estimates.

preprint2013arXiv

CANDELS Observations of the Environmental Dependence of the Color-Mass-Morphology Relation at z = 1.6

We study the environmental dependence of color, stellar mass, and morphology by comparing galaxies in a forming cluster to those in the field at z = 1:6 with Hubble Space Telescope near-infrared imaging in the CANDELS/UDS field. We quantify the morphology of the galaxies using the effective radius, reff, and Sérsic index, n. In both the cluster and field, approximately half of the bulge-dominated galaxies (n > 2) reside on the red sequence of the color-magnitude diagram, and most disk-dominated galaxies (n < 2) have colors expected for star-forming galaxies. There is weak evidence that cluster galaxies have redder rest-frame U - B colors and higher stellar masses compared to the field. Star-forming galaxies in both the cluster and field show no significant differences in their morphologies. In contrast, there is evidence that quiescent galaxies in the cluster have larger median effective radii and smaller Sérsic indices compared to the field with a significance of 2?. These differences are most pronounced for galaxies at clustercentric distances 1 Mpc < Rproj < 1.5 Mpc, which have low Sérsic indices and possibly larger effective radii, more consistent with star-forming galaxies at this epoch and in contrast to other quiescent galaxies. We argue that star-forming galaxies are processed under the influence of the cluster environment at distances greater than the cluster-halo virial radius. Our results are consistent with models where gas accretion onto these galaxies is suppressed from processes associated with the cluster environment.

preprint2013arXiv

Spitzer Spectroscopy of Infrared-Luminous Galaxies: Diagnostics of AGN and Star Formation and Contribution to Total Infrared Luminosity

We use Spitzer IRS spectroscopy to study the nature 65 IR-luminous galaxies at 0.02 < z < 0.6 with F(24micron) > 1.2mJy. The IRS spectra cover wavelengths spanning the PAH features and important atomic diagnostic lines. Our sample corresponds to L(IR) = L(8-1000micron) = 10^10-10^12 Lsun. We divide our galaxies into those with Spitzer IRAC colors indicative of warm dust heated by an AGN (IRAGN) and those whose colors indicate star-formation processes (non-IRAGN). Compared to the non-IRAGN, the IRAGN show smaller PAH emission EWs, which we attribute to an increase in mid-IR continuum from the AGN. We find that in both the IRAGN and non-IRAGN samples, the PAH luminosities correlate strongly with the [Ne II] emission, from which we conclude that the PAH luminosity directly traces the instantaneous SFR in both the IRAGN and non-IRAGN galaxies. We compare the ratio of PAH luminosity to the total IR luminosity and show that for most IRAGN star-formation accounts for 10-50% of the L(IR). We also find no measurable difference between the PAH luminosity ratios of L(11.3)/L(7.7) and L(6.2)/L(7.7) for the IRAGN and non-IRAGN, suggesting that AGN do not significantly excite or destroy PAH molecules on galaxy-wide scales. A small subset of galaxies show excess of [O IV] emission compared to their PAH emission, which indicates the presence of heavily-obscured AGN, including 3 galaxies that are not otherwise selected as IRAGN. The low PAH emission and low [Ne II] emission of the IRAGN and [O IV]-excess objects imply they have low SFRs and their IR luminosity is dominated by processes associated with the AGN.

preprint2012arXiv

A large scale structure traced by [OII] emitters hosting a distant cluster at z=1.62

We present a panoramic narrow-band imaging survey of [OII] emitters in and around the ClG J0218.3-0510 cluster at z=1.62 with Suprime-Cam on Subaru telescope. 352 [OII] emitters were identified on the basis of narrow-band excesses and photometric redshifts. We discovered a huge filamentary structure with some clumps traced by [OII] emitters and found that the ClG J0218.3-0510 cluster is embedded in an even larger super-structure than the one reported previously. 31 [OII] emitters were spectroscopically confirmed with the detection of H-alpha and/or [OIII] emission lines by FMOS observations. In the high density regions such as cluster core and clumps, star-forming [OII] emitters show a high overdensity by a factor of more than 10 compared to the field region. Although the star formation activity is very high even in the cluster core, some massive quiescent galaxies also exits at the same time. Furthermore, the properties of the individual [OII] emitters, such as star formation rates, stellar masses and specific star formation rates, do not show a significant dependence on the local density, either. Such lack of environmental dependence is consistent with our earlier result by Hayashi et al. (2011) on a z=1.5 cluster and its surrounding region. The fact that the star-forming activity of galaxies in the cluster core is as high as that in the field at z~1.6 may suggest that the star-forming galaxies are probably just in a transition phase from a starburst mode to a quiescent mode, and are thus showing comparable level of star formation rates to those in lower density environments. We may be witnessing the start of the reversal of the local SFR--density relation due to the "biased" galaxy formation and evolution in high density regions at high this redshift, beyond which massive galaxies would be forming vigorously in a more biased way in proto-cluster cores.

preprint2012arXiv

A Tale of Dwarfs and Giants: Using a z=1.62 Cluster to Understand How the Red Sequence Grew Over The Last 9.5 Billion Years

We study the red sequence in a cluster of galaxies at z=1.62 and follow its evolution over the intervening 9.5 Gyr to the present day. Using deep YJKs imaging with the HAWK-I instrument on the VLT we identify a tight red sequence and construct its rest-frame i-band luminosity function (LF). There is a marked deficit of faint red galaxies in the cluster that causes a turnover in the LF. We compare the red sequence LF to that for clusters at z<0.8 correcting the luminosities for passive evolution. The shape of the cluster red sequence LF does not evolve between z=1.62 and z=0.6 but at z<0.6 the faint population builds up significantly. Meanwhile, between z=1.62 to 0.6 the inferred total light on the red sequence grows by a factor of about 2 and the bright end of the LF becomes more populated. We construct a simple model for red sequence evolution that grows the red sequence in total luminosity and matches the constant LF shape at z>0.6. In this model the cluster accretes blue galaxies from the field that are then quenched and subsequently allowed to merge. We find that 3--4 mergers among cluster galaxies during the 4 Gyr between z=1.62 and z=0.6 matches the observed luminosity function evolution between the two redshifts. The inferred merger rate is consistent with other studies of this cluster. Our result supports the picture that galaxy merging during the major growth phase of massive clusters is an important process in shaping the red sequence population at all luminosities.

preprint2012arXiv

CANDELS Observations of the Structural Properties and Evolution of Galaxies in a Cluster at z=1.62

We discuss the structural and morphological properties of galaxies in a z=1.62 proto-cluster using near-IR imaging data from Hubble Space Telescope Wide Field Camera 3 data of the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). The cluster galaxies exhibit a clear color-morphology relation: galaxies with colors of quiescent stellar populations generally have morphologies consistent with spheroids, and galaxies with colors consistent with ongoing star formation have disk-like and irregular morphologies. The size distribution of the quiescent cluster galaxies shows a deficit of compact (< 1kpc), massive galaxies compared to CANDELS field galaxies at z=1.6. As a result the cluster quiescent galaxies have larger average effective sizes compared to field galaxies at fixed mass at greater than 90% significance. Combined with data from the literature, the size evolution of quiescent cluster galaxies is relatively slow from z~1.6 to the present, growing as (1+z)^(-0.6+/-0.1). If this result is generalizable, then it implies that physical processes associated with the denser cluster region seems to have caused accelerated size growth in quiescent galaxies prior to z=1.6 and slower subsequent growth at z<1.6 compared to galaxies in the lower density field. The quiescent cluster galaxies at z=1.6 have higher ellipticities compared to lower redshift samples at fixed mass, and their surface-brightness profiles suggest that they contain extended stellar disks. We argue the cluster galaxies require dissipationless (i.e., gas-poor or "dry") mergers to reorganize the disk material and to match the relations for ellipticity, stellar mass, size, and color of early-type galaxies in z<1 clusters.

preprint2012arXiv

CANDELS: The Contribution of the Observed Galaxy Population to Cosmic Reionization

We present measurements of the specific ultraviolet luminosity density from a sample of 483 galaxies at 6<z<8. These galaxies were selected from new deep near-infrared HST imaging from the CANDELS, HUDF09 and ERS programs. In contrast to the majority of previous analyses, which assume that the distribution of galaxy ultraviolet (UV) luminosities follows a Schechter distribution, and that the distribution continues to luminosities far below our observable limit, we investigate the contribution to reionization from galaxies which we can observe, free from these assumptions. We find that the observable population of galaxies can sustain a fully reionized IGM at z=6, if the average ionizing photon escape fraction (f_esc) is ~30%. A number of previous studies have measured UV luminosity densities at these redshifts that vary by 5X, with many concluding that galaxies could not complete reionization by z=6 unless a large population of galaxies fainter than the detection limit were invoked, or extremely high values of f_esc were present. The observed UV luminosity density from our observed galaxy samples at z=7-8 is not sufficient to maintain a fully reionized IGM unless f_esc>50%. Combining our observations with constraints on the emission rate of ionizing photons from Ly-alpha forest observations at z=6, we can constrain f_esc<34% (2-sigma) if the observed galaxies are the only contributors to reionization, or <13% (2-sigma) if the luminosity function extends to M_UV = -13. These escape fractions are sufficient to complete reionization by z=6. These constraints imply that the volume ionized fraction of the IGM becomes less than unity at z>7, consistent with a number of complementary reionization probes. If faint galaxies dominate reionization, future JWST observations will probe deep enough to see them, providing an indirect constraint on the ionizing photon escape fraction [abridged].

preprint2012arXiv

CANDELS: The Evolution of Galaxy Rest-Frame Ultraviolet Colors from z = 8 to 4

We study the evolution of galaxy rest-frame ultraviolet (UV) colors in the epoch 4 < z < 8. We use new wide-field near-infrared data in GOODS-S from the CANDELS, HUDF09 and ERS programs to select galaxies via photometric redshift measurements. Our sample consists of 2812 candidate galaxies at z > 3.5, including 113 at z = 7 to 8. We fit the observed spectral energy distribution to a suite of synthetic stellar population models, and measure the value of the UV spectral slope (beta) from the best-fit model spectrum. The median value of beta evolves significantly from -1.82 (+0.00,-0.04) at z = 4, to -2.37 (+0.26,-0.06) at z = 7. Additionally, we find that faint galaxies at z = 7 have beta = -2.68 (+0.39,-0.24) (~ -2.4 after correcting for observational bias); this is redder than previous claims in the literature, and does not require "exotic" stellar populations to explain their colors. This evolution can be explained by an increase in dust extinction, with the timescale consistent with low-mass AGB stars forming the bulk of the dust. We find no significant (< 2-sigma) correlation between beta and M_UV when measuring M_UV at a consistent rest-frame wavelength of 1500 A. This is particularly true at bright magnitudes, though our results do show evidence for a weak correlation at faint magnitudes when galaxies in the HUDF are considered separately, hinting that dynamic range in sample luminosities may play a role. We do find a strong correlation between beta and the stellar mass at all redshifts, in that more massive galaxies exhibit redder colors. The most massive galaxies in our sample have red colors at each redshift, implying that dust can build up quickly in massive galaxies, and that feedback is likely removing dust from low-mass galaxies at z > 7. Thus the stellar-mass - metallicity relation, previously observed up to z ~ 3, may extend out to z = 7 - 8.

preprint2012arXiv

First Results from Z-FOURGE: Discovery of a Candidate Cluster at z = 2.2 in COSMOS

We report the first results from the Z-FOURGE survey: the discovery of a candidate galaxy cluster at z = 2.2 consisting of two compact overdensities with red galaxies detected at >20σ above the mean surface density. The discovery was made possible by a new deep (Ks < 24.8 AB 5σ) Magellan/FourStar near-IR imaging survey with 5 custom medium-bandwidth filters. The filters pinpoint the location of the Balmer/4000Å break in evolved stellar populations at 1.5 < z < 3.5, yielding significantly more accurate photometric redshifts than possible with broadband imaging alone. The overdensities are within 1' of each other in the COSMOS field and appear to be embedded in a larger structure that contains at least one additional overdensity (\sim10σ). Considering the global properties of the overdensities, the z = 2.2 system appears to be the most distant example of a galaxy cluster with a population of red galaxies. A comparison to a large ΛCDM simulation suggests that the system may consist of merging subclusters, with properties in between those of z > 2 protoclusters with more diffuse distributions of blue galaxies and the lower-redshift galaxy clusters with prominent red sequences. The structure is completely absent in public optical catalogs in COSMOS and only weakly visible in a shallower near-IR survey. The discovery showcases the potential of deep near-IR surveys with medium-band filters to advance the understanding of environment and galaxy evolution at z > 1.5.

preprint2012arXiv

How Do Star-Forming Galaxies at Z>3 Assemble Their Masses?

We investigate how star-forming galaxies typically assemble their masses at high redshift. Using the deep multi-wavelength coverage of the GOODS dataset, we measure stellar mass of a large sample of star-forming galaxies at z~4 and 5, and make a robust determination of stellar mass function (SMF). We report a broad correlation between stellar mass and UV luminosity, such that more UV-luminous galaxies are, on average, more massive. However, the correlation has a substantial intrinsic scatter evidenced by a non-negligible number of UV-faint but massive galaxies. Furthermore, the low-mass end of the SMF does not rise as steeply as the UV luminosity function (alpha_UVLF} -(1.7-1.8), alpha_SMF -(1.3-1.4)) of the same galaxies. In a smooth formation scenario where star formation (SF) is sustained at the observed rates for a long time, these galaxies would have accumulated more mass (by a factor of ~3) than observed and therefore the SMF would mirror more closely that of the UVLF. The relatively shallow slope of the SMF is due to the fact that many of the UV-selected galaxies are not massive enough, and therefore are too faint in their rest-frame optical bands, to be detected in the current observations. Our results favor a episodic formation history in which SF in low-mass galaxies vary significantly over time, a scenario favored by galaxy clustering. Our findings for the UV-faint galaxies are in contrast with those found for more UV-luminous galaxies, which exhibit tighter SFR-M_star correlations. The discrepancy may suggest that galaxies at different luminosities may have different evolutionary paths.

preprint2012arXiv

Luminous and High Stellar Mass Candidate Galaxies at z ~ 8 Discovered in CANDELS

One key goal of the Hubble Space Telescope Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey is to track galaxy evolution back to z ~ 8. Its two-tiered "wide and deep" strategy bridges significant gaps in existing near-infrared surveys. Here we report on z ~ 8 galaxy candidates selected as F105W-band dropouts in one of its deep fields, which covers 50.1 square arcmin to 4 ks depth in each of three near-infrared bands in the Great Observatories Origins Deep Survey southern field. Two of our candidates have J<26.2 mag, and are > 1 mag brighter than any previously known F105W-dropouts. We derive constraints on the bright-end of the rest-frame ultraviolet luminosity function of galaxies at z ~ 8, and show that the number density of such very bright objects is higher than expected from the previous Schechter luminosity function estimates at this redshift. Another two candidates are securely detected in Spitzer Infrared Array Camera images, which are the first such individual detections at z ~ 8. Their derived stellar masses are on the order of a few x 10^9 M_sun, from which we obtain the first measurement of the high-mass end of the galaxy stellar mass function at z ~ 8. The high number density of very luminous and very massive galaxies at z ~ 8, if real, could imply a large stellar-to-halo mass ratio and an efficient conversion of baryons to stars at such an early time.

preprint2012arXiv

The Massive Distant Clusters of WISE Survey: The First Distant Galaxy Cluster Discovered by WISE

We present spectroscopic confirmation of a z=0.99 galaxy cluster discovered using data from the Wide-field Infrared Survey Explorer (WISE). This is the first z~1 cluster candidate from the Massive Distant Clusters of WISE Survey (MaDCoWS) to be confirmed. It was selected as an overdensity of probable z>~1 sources using a combination of WISE and SDSS-DR8 photometric catalogs. Deeper follow-up imaging data from Subaru and WIYN reveal the cluster to be a rich system of galaxies, and multi-object spectroscopic observations from Keck confirm five cluster members at z=0.99. The detection and confirmation of this cluster represents a first step towards constructing a uniformly-selected sample of distant, high-mass galaxy clusters over the full extragalactic sky using WISE data.

preprint2012arXiv

What turns galaxies off? The different morphologies of star-forming and quiescent galaxies since z~2 from CANDELS

We use HST/WFC3 imaging from the CANDELS Multicycle Treasury Survey, in conjunction with the Sloan Digital Sky Survey, to explore the evolution of galactic structure for galaxies with stellar masses >3e10M_sun from z=2.2 to the present epoch, a time span of 10Gyr. We explore the relationship between rest-frame optical color, stellar mass, star formation activity and galaxy structure. We confirm the dramatic increase from z=2.2 to the present day in the number density of non-star-forming galaxies above 3e10M_sun reported by others. We further find that the vast majority of these quiescent systems have concentrated light profiles, as parametrized by the Sersic index, and the population of concentrated galaxies grows similarly rapidly. We examine the joint distribution of star formation activity, Sersic index, stellar mass, inferred velocity dispersion, and stellar surface density. Quiescence correlates poorly with stellar mass at all z<2.2. Quiescence correlates well with Sersic index at all redshifts. Quiescence correlates well with `velocity dispersion' and stellar surface density at z>1.3, and somewhat less well at lower redshifts. Yet, there is significant scatter between quiescence and galaxy structure: while the vast majority of quiescent galaxies have prominent bulges, many of them have significant disks, and a number of bulge-dominated galaxies have significant star formation. Noting the rarity of quiescent galaxies without prominent bulges, we argue that a prominent bulge (and perhaps, by association, a supermassive black hole) is an important condition for quenching star formation on galactic scales over the last 10Gyr, in qualitative agreement with the AGN feedback paradigm.

preprint2011arXiv

AKARI Observation of the North Ecliptic Pole (NEP) Supercluster at z = 0.087: mid-infrared view of transition galaxies

We present the mid-infrared (MIR) properties of galaxies within a supercluster in the North Ecliptic Pole region at z?0.087 observed with the AKARI satellite. We use data from the AKARI NEP-Wide (5.4 deg2) IR survey and the CLusters of galaxies EVoLution studies (CLEVL) mission program. We show that near-IR (3 μm)-mid- IR (11 μm) color can be used as an indicator of the specific star formation rate and the presence of intermediate age stellar populations. From the MIR observations, we find that red-sequence galaxies consist not only of passively evolving red early-type galaxies, but also of 1) "weak-SFG" (disk-dominated star-forming galaxies which have star formation rates lower by \sim 4 \times than blue-cloud galaxies), and 2) "intermediate- MXG" (bulge-dominated galaxies showing stronger MIR dust emission than normal red early-type galaxies). Those two populations can be a set of transition galaxies from blue, star-forming, late-type galaxies evolving into red, quiescent, early-type ones. We find that the weak-SFG are predominant at intermediate masses (1010M\odot < M\star < 1010.5M\odot) and are typically found in local densities similar to the outskirts of galaxy clusters. As much as 40% of the supercluster member galaxies in this mass range can be classified as weak-SFGs, but their proportion decreases to < 10% at larger masses (M\star > 1010.5 M\odot) at any galaxy density. The fraction of the intermediate-MXG among red- sequence galaxies at 1010M\odot < M\star < 1011M\odot also decreases as the density and mass increase. In particular, \sim42% of the red-sequence galaxies with early-type morphologies are classified as intermediate-MXG at intermediate densities. These results suggest that the star formation activity is strongly dependent on the stellar mass, but that the morphological transformation is mainly controlled by the environment.

preprint2011arXiv

Color and Stellar Population Gradients in Passively Evolving Galaxies at z~2 from HST/WFC3 Deep Imaging in the Hubble Ultra Deep Field

We report the detection of color gradients in six massive (stellar mass > 10^{10} M_{sun}) and passively evolving (specific SFR < 10^{-11}/yr) galaxies at redshift 1.3<z<2.5 identified in the HUDF using HST ACS and WFC3/IR images. After matching different PSFs, we obtain color maps and multi-band optical/near-IR photometry (BVizYJH) in concentric annuli, from the smallest resolved radial (~1.7 kpc) up to several times the H-band effective radius. We find that the inner regions of these galaxies have redder rest-frame UV--optical colors than the outer parts. The slopes of the color gradients mildly depend on the overall dust obscuration and rest-frame (U-V) color, with more obscured or redder galaxies having steeper color gradients. The z~2 color gradients are also steeper than those of local early-types. The gradient of a single parameter (age, extinction or metallicity) cannot fully explain the observed color gradients. Fitting spatially resolved HST seven-band photometry to stellar population synthesis models, we find that, regardless of assumptions for metallicity gradient, the redder inner regions of the galaxies have slightly higher dust obscuration than the bluer outer regions, although the magnitude depends on the assumed extinction law. The derived age gradient depends on the assumptions for metallicity gradient. We discuss the implications of a number of assumptions for metallicity gradient on the formation and evolution of these galaxies. We find that the evolution of the mass--size relationship from z~2 to z~0 cannot be driven by in--situ extended star formation, implying that accretion or merger is mostly responsible for the evolution. The lack of a correlation between color gradient and stellar mass argues against the metallicity gradient predicted by the monolithic collapse, which would require significant major mergers to evolve into the one observed at z~0. (Abridged)

preprint2011arXiv

Extended Schmidt Law: Role Of Existing Stars In Current Star Formation

We propose an "extended Schmidt law" with explicit dependence of the star formation efficiency (SFE=SFR/Mgas) on the stellar mass surface density. This relation has a power-law index of 0.48+-0.04 and an 1-sigma observed scatter on the SFE of 0.4 dex, which holds over 5 orders of magnitude in the stellar density for individual global galaxies including various types especially the low-surface-brightness (LSB) galaxies that deviate significantly from the Kennicutt-Schmidt law. When applying it to regions at sub-kpc resolution of a sample of 12 spiral galaxies, the extended Schmidt law not only holds for LSB regions but also shows significantly smaller scatters both within and across galaxies compared to the Kennicutt-Schmidt law. We argue that this new relation points to the role of existing stars in regulating the SFE, thus encoding better the star formation physics. Comparison with physical models of star formation recipes shows that the extended Schmidt law can be reproduced by some models including gas free-fall in a stellar-gravitational potential and pressure-supported star formation. By implementing this new law into the analytic model of gas accretion in Lambda CDM, we show that it can re-produce the observed main sequence of star-forming galaxies (a relation between the SFR and stellar mass) from z=0 up to z=2.

preprint2011arXiv

Probing the Star Formation History and Initial Mass Function of the z~2.5 Lensed Galaxy SMM J163554.2+661225 with Herschel

We present the analysis of Herschel SPIRE far-infrared (FIR) observations of the z = 2.515 lensed galaxy SMM J163554.2+661225. Combining new 250, 350, and 500 micron observations with existing data, we make an improved fit to the FIR spectral energy distribution (SED) of this galaxy. We find a total infrared (IR) luminosity of L(8--1000 micron) = 6.9 +/- 0.6x10^11 Lsol; a factor of 3 more precise over previous L_IR estimates for this galaxy, and one of the most accurate measurements for any galaxy at these redshifts. This FIR luminosity implies an unlensed star formation rate (SFR) for this galaxy of 119 +/- 10 Msol per yr, which is a factor of 1.9 +/- 0.35 lower than the SFR derived from the nebular Pa-alpha emission line (a 2.5-sigma discrepancy). Both SFR indicators assume identical Salpeter initial mass functions (IMF) with slope Gamma=2.35 over a mass range of 0.1 - 100 Msol, thus this discrepancy suggests that more ionizing photons may be necessary to account for the higher Pa-alpha-derived SFR. We examine a number of scenarios and find that the observations can be explained with a varying star formation history (SFH) due to an increasing star formation rate (SFR), paired with a slight flattening of the IMF. If the SFR is constant in time, then larger changes need to be made to the IMF by either increasing the upper-mass cutoff to ~ 200 Msol, or a flattening of the IMF slope to 1.9 +/- 0.15, or a combination of the two. These scenarios result in up to double the number of stars with masses above 20 Msol, which produce the requisite increase in ionizing photons over a Salpeter IMF with a constant SFH.

preprint2011arXiv

The mid-IR Luminosity Function at z<0.3 from 5MUSES: Understanding the Star-formation/AGN Balance from a Spectroscopic View

We present rest-frame 15 and 24 um luminosity functions and the corresponding star-forming luminosity functions at z<0.3 derived from the 5MUSES sample. Spectroscopic redshifts have been obtained for ~98% of the objects and the median redshift is ~0.12. The 5-35 um IRS spectra allow us to estimate accurately the luminosities and build the luminosity functions. Using a combination of starburst and quasar templates, we quantify the star-formation and AGN contributions in the mid-IR SED. We then compute the star-formation luminosity functions at 15 um and 24 um, and compare with the total 15 um and 24 um luminosity functions. When we remove the contribution of AGN, the bright end of the luminosity function exhibits a strong decline, consistent with the exponential cutoff of a Schechter function. Integrating the differential luminosity function, we find that the fractional contribution by star formation to the energy density is 58% at 15 um and 78% at 24 um, while it goes up to ~86% when we extrapolate our mid-IR results to the total IR luminosity density. We confirm that the active galactic nuclei play more important roles energetically at high luminosities. Finally, we compare our results with work at z~0.7 and confirm that evolution on both luminosity and density is required to explain the difference in the LFs at different redshifts.

preprint2010arXiv

Dust Extinction and Metallicities of Star-Forming Lyman Alpha Emitting Galaxies at Low Redshift

We present the results of an optical spectroscopic study of 12 GALEX-discovered star-forming Lyman alpha emitting galaxies (LAEs) at z ~ 0.3. We measure the emission line fluxes from these galaxies by fitting their observed spectra to stellar population models in order to correct for underlying stellar absorption. We revisit earlier stellar population model fitting results, finding that excluding now-known AGNs lowers the typical stellar population age and stellar mass of this sample to ~ 300 Myr and ~ 4 x 10^9 Msol, respectively. We calculate their dust extinction using the Balmer decrement, and find a typical visual attenuation of Av ~ 0.3 mag, similar to that seen in many high-redshift LAEs. Comparing the ratio of Lyalpha/Halpha and the Lyman alpha equivalent widths to the measured dust extinction, we find that the ISMs in these objects appear to be neither enhancing nor seriously attenuating the Lyman alpha equivalent widths, as would be the case in a quasi-clumpy ISM. Lastly, we perform a detailed analysis of the gas-phase metallicities of these galaxies, and we find that most galaxies in our sample have Z < 0.4 Zsol. We find that at a fixed stellar mass, these low-redshift LAE analogs are offset by ~ 0.6 dex lower in metallicity from the general galaxy population at similar redshifts based on the local mass-metallicity relationship. This implies that galaxies with Lyman alpha in emission may be systematically more metal poor than star-forming galaxies at the same stellar mass and redshift, similar to preliminary results at z ~ 2.

preprint2010arXiv

Infrared Luminosities and Aromatic Features in the 24um Flux Limited Sample of 5MUSES

We study a 24\,$μ$m selected sample of 330 galaxies observed with the Infrared Spectrograph for the 5\,mJy Unbiased Spitzer Extragalactic Survey. We estimate accurate total infrared luminosities by combining mid-IR spectroscopy and mid-to-far infrared photometry, and by utilizing new empirical spectral templates from {\em Spitzer} data. The infrared luminosities of this sample range mostly from 10$^9$L$_\odot$ to $10^{13.5}$L$_\odot$, with 83% in the range 10$^{10}$L$_\odot$$<$L$_{\rm IR}$$<10^{12}$L$_\odot$. The redshifts range from 0.008 to 4.27, with a median of 0.144. The equivalent widths of the 6.2\,$μ$m aromatic feature have a bimodal distribution. We use the 6.2\,$μ$m PAH EW to classify our objects as SB-dominated (44%), SB-AGN composite (22%), and AGN-dominated (34%). The high EW objects (SB-dominated) tend to have steeper mid-IR to far-IR spectral slopes and lower L$_{\rm IR}$ and redshifts. The low EW objects (AGN-dominated) tend to have less steep spectral slopes and higher L$_{\rm IR}$ and redshifts. This dichotomy leads to a gross correlation between EW and slope, which does not hold within either group. AGN dominated sources tend to have lower log(L$_{\rm PAH 7.7μm}$/L$_{\rm PAH 11.3μm}$) ratios than star-forming galaxies, possibly due to preferential destruction of the smaller aromatics by the AGN. The log(L$_{\rm PAH 7.7μm}$/L$_{\rm PAH 11.3μm}$) ratios for star-forming galaxies are lower in our sample than the ratios measured from the nuclear spectra of nearby normal galaxies, most probably indicating a difference in the ionization state or grain size distribution between the nuclear regions and the entire galaxy. Finally, we provide a calibration relating the monochromatic 5.8, 8, 14 and 24um continuum or Aromatic Feature luminosity to L$_{\rm IR}$ for different types of objects.

preprint2010arXiv

On the Stellar Populations and Evolution of Star-Forming Galaxies at 6.3 < z < 8.6

We study the physical characteristics of galaxies at 6.3 < z < 8.6, selected from deep near-infrared imaging with the Wide Field Camera 3 (WFC3) on board the Hubble Space Telescope. Accounting for the photometric scatter using simulations, galaxies at z ~ 7 have bluer UV colors compared to typical local starburst galaxies at > 4 sigma confidence. Although these colors necessitate young ages (<100 Myr), low or zero dust attenuation, and low metallicities, these are explicable by normal (albeit unreddened) stellar populations, with no evidence for near-zero metallicities and/or top-heavy initial mass functions. The age of the Universe at these redshifts limits the amount of stellar mass in late-type populations, and the WFC3 photometry implies galaxy stellar masses ~ 10^8 - 10^9 Msol for Salpeter initial mass functions to a limiting magnitude of M_1500 ~ -18. The masses of ``characteristic'' (L*) z > 7 galaxies are smaller than those of L* Lyman break galaxies (LBGs) at lower redshifts, and are comparable to less evolved galaxies selected on the basis of their Lyman alpha emission at 3 < z < 6, implying that the 6.3 < z < 8.6 galaxies are the progenitors of more evolved galaxies at lower redshifts. We estimate that Lyman alpha emission is able to contribute to the observed WFC3 colors of galaxies at these redshifts, with an estimated typical line flux of ~ 10^-18 erg s^-1 cm^-2, roughly a factor of four below currently planned surveys. The integrated UV specific luminosity for the detected galaxies at z ~ 7 and z ~ 8 is within factors of a few of that required to reionize the IGM assuming low clumping factors, implying that in order to reionize the Universe galaxies at these redshifts have a high ( ~ 50%) escape fraction of Lyman continuum photons, possibly substantiated by the very blue colors of this population.

preprint2010arXiv

Spitzer IRS 16 micron Observations of the GOODS Fields

We present Spitzer 16 micron imaging of the Great Observatories Origins Deep Survey (GOODS) fields. We survey 150 square arcminutes in each of the two GOODS fields (North and South), to an average 3 sigma depth of 40 and 65 micro-Jy respectively. We detect about 1300 sources in both fields combined. We validate the photometry using the 3-24 micron spectral energy distribution of stars in the fields compared to Spitzer spectroscopic templates. Comparison with ISOCAM and AKARI observations in the same fields show reasonable agreement, though the uncertainties are large. We provide a catalog of photometry, with sources cross correlated with available Spitzer, Chandra, and HST data. Galaxy number counts show good agreement with previous results from ISOCAM and AKARI, with improved uncertainties. We examine the 16 to 24 micron flux ratio and find that for most sources it lies within the expected locus for starbursts and infrared luminous galaxies. A color cut of S_{16}/S_{24}>1.4 selects mostly sources which lie at 1.1<z<1.6, where the 24 micron passband contains both the redshifted 9.7 micron silicate absorption and the minimum between PAH emission peaks. We measure the integrated galaxy light of 16 micron sources, and find a lower limit on the galaxy contribution to the extragalactic background light at this wavelength to be 2.2\pm 0.2$ nW m^{-2} sr^{-1}.

preprint2010arXiv

The Evolution of the Star Formation Rate of Galaxies at 0.0 < z < 1.2

We present the 24 micron rest-frame luminosity function (LF) of star-forming galaxies in the redshift range 0.0 < z < 0.6 constructed from 4047 spectroscopic redshifts from the AGN and Galaxy Evolution Survey of 24 micron selected sources in the Bootes field of the NOAO Deep Wide-Field Survey. This sample provides the best available combination of large area (9 deg^2), depth, and statistically complete spectroscopic observations, allowing us to probe the evolution of the 24 micron LF of galaxies at low and intermediate redshifts while minimizing the effects of cosmic variance. In order to use the observed 24 micron luminosity as a tracer for star formation, active galactic nuclei (AGNs) that could contribute significantly at 24 micron are identified and excluded from our star-forming galaxy sample based on their mid-IR spectral energy distributions or the detection of X-ray emission. The evolution of the 24 micron LF of star-forming galaxies for redshifts of z < 0.65 is consistent with a pure luminosity evolution where the characteristic 24 micron luminosity evolves as (1+z)^(3.8+/-0.3). We extend our evolutionary study to encompass 0.0 < z < 1.2 by combining our data with that of the Far-Infrared Deep Extragalactic Legacy Survey. Over this entire redshift range the evolution of the characteristic 24 micron luminosity is described by a slightly shallower power law of (1+z)^(3.4+/-0.2). We find a local star formation rate density of (1.09+/-0.21) x 10^-2 Msun/yr/Mpc^-3, and that it evolves as (1+z)^(3.5+/-0.2) over 0.0 < z < 1.2. These estimates are in good agreement with the rates using optical and UV fluxes corrected for the effects of intrinsic extinction in the observed sources. This agreement confirms that star formation at z <~ 1.2 is robustly traced by 24 micron observations and that it largely occurs in obscured regions of galaxies. (Abridged)

preprint2010arXiv

The Rising Star-Formation Histories of Distant Galaxies and Implications for Gas Accretion with Time

Distant galaxies show correlations between their current star-formation rates (SFRs) and stellar masses, implying that their star-formation histories (SFHs) are highly similar. Moreover, observations show that the UV luminosities and stellar masses grow from z=8 to 3, implying that the SFRs increase with time. We compare the cosmologically averaged evolution in galaxies at 3 < z < 8 at constant comoving number density, n = 2 x 10^-4 Mpc^-3. This allows us to study the evolution of stellar mass and star formation in the galaxy predecessors and descendants in ways not possible using galaxies selected at constant stellar mass or SFR, quantities that evolve strongly in time. We show that the average SFH of these galaxies increase smoothly from z=8 to 3 as SFR ~ t^alpha with alpha = 1.7 +/- 0.2. This conflicts with assumptions that the SFR is either constant or declines exponentially in time. We show that the stellar mass growth in these galaxies is consistent with this derived SFH. This provides evidence that the slope of the high-mass end of the IMF is approximately Salpeter unless the duty cycle of star formation is much less than unity. We argue that these relations follow from gas accretion (either through accretion or delivered by mergers) coupled with galaxy disk growth under the assumption that the SFR depends on the local gas surface density. This predicts that gas fractions decrease from z=8 to 3 on average as f_gas ~ (1+z)^0.9 for galaxies with this number density. The implied galaxy gas accretion rates at z > 4 are as fast and may even exceed the SFR: this is the "gas accretion epoch". At z < 4 the SFR overtakes the implied gas accretion rate, indicating a period where galaxies consume gas faster than it is acquired. At z < 3, galaxies with this number density depart from these relations implying that star formation and gas accretion are slowed at later times.

preprint2009arXiv

History of Galaxy Interactions and their Impact on Star Formation over the Last 7 Gyr from GEMS

We perform a comprehensive estimate of the frequency of galaxy mergers and their impact on star formation over z~0.24--0.80 (lookback time T_b~3--7 Gyr) using 3698 (M*>=1e9 Msun) galaxies with GEMS HST, COMBO-17, and Spitzer data. Our results are: (1) Among 790 high mass (M*>=2.5e10 Msun) galaxies, the visually-based merger fraction over z~0.24--0.80, ranges from 9%+-5% to 8%+-2%. Lower limits on the major and minor merger fractions over this interval range from 1.1% to 3.5%, and 3.6% to 7.5%, respectively. This is the first approximate empirical estimate of the frequency of minor mergers at z<1. For a visibility timescale of ~0.5 Gyr, it follows that over T_b~3--7 Gyr, ~68% of high mass systems have undergone a merger of mass ratio >1/10, with ~16%, 45%, and 7% of these corresponding respectively to major, minor, and ambiguous `major or minor' mergers. The mean merger rate is a few x 1e-4 Gyr-1 Mpc-3. (2) We compare the empirical merger fraction and rate for high mass galaxies to a suite of Lambda CDM-based models: halo occupation distribution models, semi-analytic models, and hydrodynamic SPH simulations. We find qualitative agreement between observations and models such that the (major+minor) merger fraction or rate from different models bracket the observations, and show a factor of five dispersion. Near-future improvements can now start to rule out certain merger scenarios. (3) Among ~3698 M*>=1e9 Msun galaxies, we find that the mean SFR of visibly merging systems is only modestly enhanced compared to non-interacting galaxies over z~0.24--0.80. Visibly merging systems only account for less than 30% of the cosmic SFR density over T_b~3--7 Gyr. This suggests that the behavior of the cosmic SFR density over the last 7 Gyr is predominantly shaped by non-interacting galaxies.

preprint2009arXiv

Less than 10 percent of star formation in z=0.6 massive galaxies is triggered by major interactions

Both observations and simulations show that major tidal interactions or mergers between gas-rich galaxies can lead to intense bursts of starformation. Yet, the average enhancement in star formation rate (SFR) in major mergers and the contribution of such events to the cosmic SFR are not well estimated. Here we use photometric redshifts, stellar masses and UV SFRs from COMBO-17, 24 micron SFRs from Spitzer and morphologies from two deep HST cosmological survey fields (ECDFS/GEMS and A901/STAGES) to study the enhancement in SFR as a function of projected galaxy separation. We apply two-point projected correlation function techniques, which we augment with morphologically-selected very close pairs (separation <2 arcsec) and merger remnants from the HST imaging. Our analysis confirms that the most intensely star-forming systems are indeed interacting or merging. Yet, for massive (M* > 10^10 Msun) star-forming galaxies at 0.4<z<0.8, we find that the SFRs of galaxies undergoing a major interaction (mass ratios <1:4 and separations < 40 kpc) are only 1.80 +/- 0.30 times higher than the SFRs of non-interacting galaxies when averaged over all interactions and all stages of the interaction, in good agreement with other observational works. We demonstrate that these results imply that <10% of star formation at 0.4 < z < 0.8 is triggered directly by major mergers and interactions; these events are not important factors in the build-up of stellar mass since z=1.

preprint2009arXiv

Optically-passive spirals: The missing link in gradual star formation suppression upon cluster infall

Galaxies migrate from the blue cloud to the red sequence when their star formation is quenched. Here, we report on galaxies quenched by environmental effects and not by mergers or strong AGN as often invoked: They form stars at a reduced rate which is optically even less conspicuous, and manifest a transition population of blue spirals evolving into S0 galaxies. These 'optically passive' or 'red spirals' are found in large numbers in the STAGES project (and by Galaxy Zoo) in the infall region of clusters and groups.

preprint2009arXiv

Paschen-alpha Emission in the Gravitationally Lensed Galaxy SMM J163554.2+661225

We report the detection of the Paschen-alpha emission line in the z=2.515 galaxy SMM J163554.2+661225 using Spitzer spectroscopy. SMM J163554.2+661225 is a sub-millimeter-selected infrared (IR)-luminous galaxy maintaining a high star-formation rate (SFR), with no evidence of an AGN from optical or infrared spectroscopy, nor X-ray emission. This galaxy is lensed gravitationally by the cluster Abell 2218, making it accessible to Spitzer spectroscopy. Correcting for nebular extinction derived from the H-alpha and Pa-alpha lines, the dust-corrected luminosity is L(Pa-alpha) = (2.57+/-0.43) x 10^43 erg s^-1, which corresponds to an ionization rate, Q = (1.6+/-0.3) x 10^55 photons s^-1. The instantaneous SFR is 171+/-28 solar masses per year, assuming a Salpeter-like initial mass function. The total IR luminosity derived using 70, 450, and 850 micron data is L(IR) = (5-10) x 10^11 solar luminosities, corrected for gravitational lensing. This corresponds to a SFR=90-180 solar masses per year, where the upper range is consistent with that derived from the Paschen-alpha luminosity. While the L(8 micron) / L(Pa-alpha) ratio is consistent with the extrapolated relation observed in local galaxies and star-forming regions, the rest-frame 24 micron luminosity is significantly lower with respect to local galaxies of comparable Paschen-alpha luminosity. Thus, SMM J163554.2+661225 arguably lacks a warmer dust component (T ~ 70 K), which is associated with deeply embedded star formation, and which contrasts with local galaxies with comparable SFRs. Rather, the starburst is consistent with star-forming local galaxies with intrinsic luminosities, L(IR) ~ 10^10 solar luminosities, but "scaled-up" by a factor of 10-100.

preprint2007arXiv

The Warm Spitzer Mission: Prospects for Studies of the Distant Universe

IRAC excels at detecting distant objects. Due to a combination of the shapes of the spectral energy distributions of galaxies and the low background achieved from space, IRAC reaches greater depth in comparable exposure time at 3.6 and 4.5 micron than any ground- or space-based facility currently can at 2.2 micron. Furthermore, the longer wavelengths probed by IRAC enable studies of the rest-frame optical and near-infrared light of galaxies and AGN to much higher redshift than is possible from the ground. This white paper explores the merits of different survey strategies for studying the distant universe during the warm mission. A three-tiered approach serves a wide range of science goals and uses the spacecraft effectively: 1) an ultra-deep survey of ~0.04 square degrees to a depth of ~250 hrs (in conjunction with an HST/WFC3 program), to study the Universe at 7<z<14; 2) a survey of ~2 square degrees to the GOODS depth of 20 hrs, to identify luminous galaxies at z>6 and characterize the relation between the build-up of dark matter halos and their constituent galaxies at 2<z<6, and 3) a 500 square degree survey to the SWIRE depth of 120 s, to systematically study large scale structure at 1<z<2 and characterize high redshift AGN. One or more of these programs could conceivably be implemented by the SSC, following the example of the Hubble Deep Field campaigns. As priorities in this field continuously shift it is also crucial that a fraction of the exposure time remains unassigned, thus enabling science that will reflect the frontiers of 2010 and beyond rather than those of 2007.

preprint2006arXiv

A MMT/Hectospec Redshift Survey of 24 Micron Sources in the Spitzer First Look Survey

We present a spectroscopic survey using the MMT/Hectospec fiber spectrograph of 24 micron sources selected with the Spitzer Space Telescope in the Spitzer First Look Survey. We report 1296 new redshifts for 24 micron sources, including 599 with f(24micron) > 1 mJy. Combined with 291 additional redshifts for sources from the Sloan Digital Sky Survey (SDSS), our observing program was highly efficient and is ~90% complete for i' < 21 mag and f(24micron) > 1 mJy, and is 35% complete for i' < 20.5 mag and 0.3 mJy < f(24micron) < 1 mJy. Our Hectospec survey includes 1078 and 168 objects spectroscopically classified as galaxies and QSOs, respectively. Combining the Hectospec and SDSS samples, we find 24 micron-selected galaxies to z < 0.98 and QSOs to z < 3.6, with mean redshifts of <z(gal)> = 0.27 and <z(QSO)> =1.1. As part of this publication, we include the redshift catalogs and the reduced spectra; these are also available online (http://mips.as.arizona.edu/~papovich/fls) and through the NASA/IPAC Infrared Science Archive (http://irsa.ipac.caltech.edu).

preprint2006arXiv

The Star Formation History and Stellar Assembly of High Redshift Galaxies

I discuss current observational constraints on the star-formation and stellar-assembly histories of galaxies at high redshifts. The data on massive galaxies at z<1 implies that their stellar populations formed at z>2, and that their morphological configuration was in place soon thereafter. Spitzer Space Telescope 24 micron observations indicate that a substantial fraction of massive galaxies at z ~ 1.5-3 have high IR luminosities, suggesting they are rapidly forming stars, accreting material onto supermassive black holes, or both. I compare how observations of these IR-active phases in the histories of massive galaxies constrain current galaxy-formation models.

preprint2002arXiv

The Evolution of the Global Stellar Mass Density at 0 < z < 3

The build-up of stellar mass in galaxies is the consequence of their past star formation and merging histories. Here we report measurements of rest-frame optical light and calculations of stellar mass at high redshift based on an infrared-selected sample of galaxies from the Hubble Deep Field North. The bright envelope of rest-frame B-band galaxy luminosities is similar from 0<z<3, and the co-moving luminosity density is constant to within a factor of 3 over that redshift range. However, galaxies at higher redshifts are bluer, and stellar population modeling indicates that they had significantly lower mass-to-light ratios than those of present-day L* galaxies. This leads to a global stellar mass density, Omega_*(z), which rises with time from z=3 to the present. This measurement essentially traces the integral of the cosmic star formation history that has been the subject of previous investigations. 50-75% of the present-day stellar mass density had formed by z ~ 1, but at z ~= 2.7 we find only 3-14% of today's stars were present. This increase in Omega_* with time is broadly consistent with observations of the evolving global star formation rate once dust extinction is taken into account, but is steeper at 1<z<3 than predicted by some recent semi-analytic models of galaxy formation. The observations appear to be inconsistent with scenarios in which the bulk of stars in present-day galactic spheroids formed at z>>2.