Source author record

Steven L. Finkelstein

Steven L. Finkelstein 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

97works
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

97 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

Little Red Dots and their Progenitors from Direct Collapse Black Holes

The James Webb Space Telescope (JWST) has discovered a new population of objects, the Little Red Dots (LRDs), characterized by V-shaped spectra indicative of strong breaks around the Balmer limit and compact morphology that gave them their name. A popular explanation is that they are a sub-population of active galactic nuclei/supermassive black holes (AGN/SMBHs) predominantly found in the high-redshift Universe ($z\gtrsim3$). Similarly, direct collapse black holes (DCBHs), theorized to form from collapsing massive, extremely metal-poor gas clouds, have been invoked to explain high-redshift quasars, the most massive AGN sub-population. Here, we employ the semi-analytical code A-SLOTH to produce a population of DCBHs and compare them against observed LRD demographics and properties. Specifically, we compare the DCBH-seeded SMBH population against the standard stellar-remnant seeds and find that DCBH models agree better with observed LRD population statistics and host halo properties. Furthermore, for the most extreme and earliest LRD detections, interpreted to be systems with an AGN but little stellar component, DCBHs are able to reproduce the observed spectral shape and properties under multiple scenarios - high dust attenuation or AGN surrounded by dense gas - that have been proposed to explain the unique shape of LRD spectra. Even when super-Eddington accretion, invoked previously to explain the nature of LRDs, is enforced on stellar remnant seeds, the spectral characteristics of extreme LRDs cannot be reproduced. We emphasize the importance of gas-metallicity observations as an additional dimension besides the widely used SMBH-stellar mass ratios to further constrain the progenitors of LRDs.

preprint2026arXiv

The BoRG-$JWST$ Survey: Abundance and Mass-to-light Ratio of Luminous $z=7-9$ Galaxies from Independent Sight Lines with NIRSpec

We present new results on the rest-frame UV luminosity function (UVLF) and stellar mass-to-light (M/L) ratio of bright (M$_{\rm UV}\lesssim-20$ mag) spectroscopically-confirmed galaxies at $z=7-9$ derived from the BoRG-$JWST$ survey, a unique data set of NIRSpec prism follow up of $HST$-selected sources from random-pointing imaging. By selecting galaxies from over 200 independent sight lines, the survey minimizes cosmic variance ensuring a statistically robust sample of the bright-galaxy population during the epoch of reionization. The data is used to constrain, for the first time, the bright end of the UVLF at $z=7-9$ from spectroscopically-confirmed galaxies over eight independent fields. We find that the bright end of the UVLF is higher than found using imaging over $JWST$ legacy fields, suggesting the latter may be significantly affected by cosmic variance, and thus reducing the tension with recent findings from $JWST$ at $z>10$ and comparable to models invoking little dust attenuation and bursty star formation. Additionally, we use the galaxies' $JWST$ spectra to infer their stellar masses and M/L ratios relative to other $HST$ and $JWST$ studies. We show that the stellar mass scales almost linearly with UV luminosity (M$_* \propto L_{\rm UV}^{0.85\pm0.12}$), albeit with large ($\sim0.5$ dex) intrinsic scatter, consistent with stochastic bursts of star formation in early galaxy formation.

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.

preprint2025arXiv

Simulation-based inference of galaxy properties from JWST pixels

We present an efficient Bayesian SED-fitting framework tailored to multiwavelength pixel photometry from the JWST Advanced Deep Extragalactic Survey (JADES). Our method employs simulation-based inference to enable rapid posterior sampling across galaxy pixels, leveraging the unprecedented spatial resolution, wavelength coverage, and depth provided by the survey. It is trained on synthetic photometry generated from MILES stellar population models, incorporating both parametric and non-parametric SFHs, realistic noise, and JADES-like filter sensitivity thresholds. We validate this amortised inference approach on mock datasets, achieving robust and well-calibrated posterior distributions, with an $R^2$ score of 0.99 for stellar mass. Applying our pipeline to real observations, we derive spatially resolved maps of stellar population properties down to $\mathrm{S/N}_{\rm{pixel}}=5$ (averaged over F277W, F356W, F444W) for 1083 JADES galaxies and ~2 million pixels with spectroscopic redshifts. These maps enable the identification of dusty or starburst regions and offer insights into mass growth and the structural assembly. We assess the outshining phenomenon by comparing pixel-based and integrated stellar mass estimates, finding limited impact only in low-mass galaxies ($<10^8M_{\odot}$) but systematic differences of ~0.20 dex linked to SFH priors. With an average posterior sampling speed of $10^{-4}$ seconds per pixel and a total inference time of ~1 CPU-day for the full dataset, our model offers a scalable solution for extracting high-fidelity stellar population properties from HST+JWST datasets, opening the way for statistical studies at sub-galactic scales.

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

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

Lyα Halos around [O III]-Selected Galaxies in HETDEX

We present extended Lyman-α (Lyα) emission out to 800 kpc of 1034 [O III]-selected galaxies at redshifts 1.9<z<2.35 using the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). The locations and redshifts of the galaxies are taken from the 3D-HST survey. The median-stacked surface brightness profile of Lyα emission of the [O III]-selected galaxies agrees well with that of 968 bright Lyα-emitting galaxies (LAEs) at r>40 kpc from the galaxy centers. The surface brightness in the inner parts (r<10 kpc) around the [O III]-selected galaxies, however, is ten times fainter than that of the LAEs. Our results are consistent with the notion that photons dominating the outer regions of the Lyα halos are not produced in the central galaxies but originate outside of them.

preprint2022arXiv

Maximal X-ray feedback in the pre-reionization universe

X-ray feedback in the pre-reionization Universe provided one of the major energy sources for reionization and the thermal evolution of the early intergalactic medium. However, X-ray sources at high redshift have remained largely inaccessible to observations. One alternative approach to study the overall effect of X-ray feedback in the early Universe is a full cosmological simulation. Toward this goal, in this paper we create an analytic model of X-ray feedback from accretion onto supermassive black holes (SMBHs), to be used as a sub-grid model in future cosmological simulations. Our analytic model provides a relation between the mass of a dark matter halo and the SMBH it hosts, where the efficiency is governed by an energy balance argument between thermal feedback and the confining gravitational potential of the halo. To calibrate the model, we couple the halo-level recipe with the Press-Schechter halo mass function and derive global mass and energy densities. We then compare our model to various observational constraints, such as the resulting soft X-ray and IR cosmic radiation backgrounds, to test our choice of model parameters. We in particular derive model parameters that do not violate any constraints, while providing maximal X-ray feedback prior to reionization. In addition, we consider the contribution of SMBH X-ray sources to reionization and the global 21 cm absorption signal.

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

Semi-analytic forecasts for JWST -- VI. Simulated lightcones and galaxy clustering predictions

In anticipation of the new era of high-redshift exploration marked by the commissioning of the James Webb Space Telescope (JWST), we present two sets of galaxy catalogues that are designed to aid the planning and interpretation of observing programs. We provide a set of 40 wide-field lightcones with footprints spanning approximately ~ 1,000 sq. arcmin, containing galaxies up to z = 10, and a new set of 8 ultra-deep lightcones with 132 sq. arcmin footprints, containing galaxies up to z ~ 12 down to the magnitudes expected to be reached in the deepest JWST surveys. These mock lightcones are extracted from dissipationless N-body simulations and populated with galaxies using the well-established, computationally efficient Santa Cruz semi-analytic model for galaxy formation. We provide a wide range of predicted physical properties, and simulated photometry from NIRCam and many other instruments. We explore the predicted counts and luminosity functions and angular two-point correlation functions for galaxies in these simulated lightcones. We also explore the predicted field-to-field variance using multiple lightcone realizations. We find that these lightcones reproduce the available measurements of observed clustering from 0.2 < z < 7.5 very well. We provide predictions for galaxy clustering at high redshift that may be obtained from future JWST observations. All of the lightcones presented here are made available through a web-based, interactive data release portal https://flathub.flatironinstitute.org/group/sam-forecasts.

preprint2022arXiv

Stellar Populations of Lyman-alpha Emitting Galaxies in the HETDEX Survey I: An Analysis of LAEs in the GOODS-N Field

We present the results of a stellar-population analysis of Lyman-alpha emitting galaxies (LAES) in GOODS-N at 1.9 < z < 3.5 spectroscopically identified by the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). We provide a method for connecting emission-line detections from the blind spectroscopic survey to imaging counterparts, a crucial tool needed as HETDEX builds a massive database of ~1 million Lyman-alpha detections. Using photometric data spanning as many as 11 filters covering 0.4-4.5 microns from the Hubble and Spitzer Space Telescopes, we study the objects' global properties and explore which properties impact the strength of Lyman-alpha emission. We measure a median stellar mass of 0.8 (^+2.9_-0.5) x 10^9 Msol and conclude that the physical properties of HETDEX spectroscopically-selected LAEs are comparable to LAEs selected by previous deep narrow band studies. We find that stellar mass and star formation rate correlate strongly with the Lyman-alpha equivalent width. We then use a known sample of z>7 LAEs to perform a proto-study of predicting Lyman-alpha emission from galaxies in the Epoch of Reionization, finding agreement at the 1-sigma level between prediction and observation for the majority of strong emitters.

preprint2022arXiv

Surface Brightness Profile of Lyman-$α$ Halos out to 320 kpc in HETDEX

We present the median-stacked Lyman-$α$ surface brightness profile of 968 spectroscopically selected Lyman-$α$ emitting galaxies (LAEs) at redshifts $1.9<z<3.5$ in the early data of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). The selected LAEs are high-confidence Lyman-$α$ detections with large signal-to-noise ratios observed with good seeing conditions (point-spread-function full-width-at-half-maximum $<1.4"$), excluding active galactic nuclei (AGN). The Lyman-$α$ luminosities of the LAEs are $10^{42.4}-10^{43}\, \mathrm{erg}\, \mathrm{s}^{-1}$. We detect faint emission in the median-stacked radial profiles at the level of $(3.6\pm 1.3)\times 10^{-20}\,\mathrm{erg}\,\mathrm{s}^{-1}\,\mathrm{cm}^{-2}\,\mathrm{arcsec}^{-2}$ from the surrounding Lyman-$α$ halos out to $r\simeq 160$ kpc (physical). The shape of the median-stacked radial profile is consistent at $r<80\,\mathrm{kpc}$ with that of much fainter LAEs at $3<z<4$ observed with the Multi Unit Spectroscopic Explorer (MUSE), indicating that the median-stacked Lyman-$α$ profiles have similar shapes at redshifts $2<z<4$ and across a factor of $10$ in Lyman-$α$ luminosity. While we agree with the results from the MUSE sample at $r<80\,\mathrm{kpc}$, we extend the profile over a factor of two in radius. At $r>80\,\mathrm{kpc}$, our profile is flatter than the MUSE model. The measured profile agrees at most radii with that of galaxies in the Byrohl et al. (2021) cosmological radiative transfer simulation at $z=3$. This suggests that the surface brightness of a Lyman-$α$ halo at $r\lesssim 100$ kpc is dominated by resonant scattering of Lyman-$α$ photons from star-forming regions in the central galaxy, whereas at $r > 100$ kpc it is dominated by photons from galaxies in surrounding dark matter halos.

preprint2022arXiv

The Active Galactic Nuclei in the Hobby-Eberly Telescope Dark Energy Experiment Survey (HETDEX) I. Sample selection

We present the first Active Galactic Nuclei (AGN) catalog in the Hobby-Eberly Telescope Dark Energy Experiment Survey (HETDEX) observed between January 2017 and June 2020. HETDEX is an ongoing spectroscopic survey with no pre-selection based on magnitudes, colors or morphologies, enabling us to select AGN based on their spectral features. Both luminous quasars and low-luminosity Seyferts are found in our catalog. AGN candidates are selected with at least two significant AGN emission lines, such as the LyA and CIV line pair, or with single broad emission lines (FWHM > 1000 km/s). Each source is further confirmed by visual inspections. This catalog contains 5,322 AGN, covering an effective sky coverage of 30.61 deg^2. A total of 3,733 of these AGN have secure redshifts, and we provide redshift estimates for the remaining 1,589 single broad-line AGN with no cross matched spectral redshifts from SDSS DR14Q. The redshift range of the AGN catalog is 0.25 < z < 4.32, with a median of z = 2.1. The bolometric luminosity range is 10^9-10^14 Lsun with a median of 10^12 Lsun. The median r-band magnitude of the AGN is 21.6 mag, with 34% of the AGN have r > 22.5, and 2.6% reaching the detection limit at r ~ 26 mag of the deepest imaging surveys we searched. We also provide a composite spectrum of the AGN sample covering 700 AA - 4400 AA.

preprint2022arXiv

The Active Galactic Nuclei in the Hobby-Eberly Telescope Dark Energy Experiment Survey (HETDEX) II. Luminosity Function

We present the LyA emission line luminosity function (LF) of the Active Galactic Nuclei (AGN) in the first release of the Hobby-Eberly Telescope Dark Energy Experiment Survey (HETDEX) AGN catalog (Liu et al. 2022, Paper I). The AGN are selected either by emission-line pairs characteristic of AGN or by single broad emission line, free of any photometric pre-selections (magnitude/color/morphology). The sample consists of 2,346 AGN spanning 1.88<z<3.53, covering an effective area of 30.61 deg^2. Approximately 2.6 of the HETDEX AGN are not detected at $>5σ$ confidence at r~26 in the deepest $r$-band images we have searched. The LyA line luminosity ranges from ~10^42.3 to ~10^45.9 erg s^-1. Our LyA LF shows a turnover luminosity with opposite slopes on the bright end and the faint end: The space density is highest at L_LyA^*=10^43.4 erg s^-1. We explore the evolution of the AGN LF over a broader redshift range (0.8<z<3); constructing the rest-frame ultraviolet (UV) LF with the 1450 AA monochromatic luminosity of the power-law component of the continuum ($\rm M_{1450}$) from M_1450~-18 to ~-27.5. We divide the sample into three redshift bins (z~1.5, 2.1, and 2.6). In all three redshift bins, our UV LFs indicate that the space density of AGN is highest at the turnover luminosity M_1450^* with opposite slopes on the bright end and the faint end. The M_1450 LFs in the three redshift bins can be well-fit with a luminosity-evolution-density-evolution (LEDE) model: the turnover luminosity (M_1450^*) increases and the turnover density (Phi^*) decreases with increasing redshift.

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.

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

Characterization of Two 2mm-detected Optically-Obscured Dusty Star-Forming Galaxies

The 2mm Mapping Obscuration to Reionization with ALMA (MORA) Survey was designed to detect high redshift ($z\gtrsim4$), massive, dusty star-forming galaxies (DSFGs). Here we present two, likely high redshift sources, identified in the survey whose physical characteristics are consistent with a class of optical/near-infrared (OIR) invisible DSFGs found elsewhere in the literature. We first perform a rigorous analysis of all available photometric data to fit spectral energy distributions and estimate redshifts before deriving physical properties based on our findings. Our results suggest the two galaxies, called MORA-5 and MORA-9, represent two extremes of the "OIR-dark" class of DSFGs. MORA-5 ($z_{\rm phot}=4.3^{+1.5}_{-1.3}$) is a significantly more active starburst with a star-formation rate of 830$^{+340}_{-190}$M$_\odot$yr$^{-1}$ compared to MORA-9 ($z_{\rm phot}=4.3^{+1.3}_{-1.0}$) whose star-formation rate is a modest 200$^{+250}_{-60}$M$_\odot$yr$^{-1}$. Based on the stellar masses (M$_{\star}\approx10^{10-11}$M$_\odot$), space density ($n\sim(5\pm2)\times10^{-6}$Mpc$^{-3}$, which incorporates two other spectroscopically confirmed OIR-dark DSFGs in the MORA sample at $z=4.6$ and $z=5.9$), and gas depletion timescales ($<1$Gyr) of these sources, we find evidence supporting the theory that OIR-dark DSFGs are the progenitors of recently discovered $3<z<4$ massive quiescent galaxies.

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

First HETDEX Spectroscopic Determinations of Ly$α$ and UV Luminosity Functions at $z=2-3$: Bridging a Gap Between Faint AGN and Bright Galaxies

We present Ly$α$ and ultraviolet-continuum (UV) luminosity functions (LFs) of galaxies and active galactic nuclei (AGN) at $z=2.0-3.5$ determined by the un-targetted optical spectroscopic survey of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). We combine deep Subaru imaging with HETDEX spectra resulting in $11.4$ deg$^2$ of fiber-spectra sky coverage, obtaining $18320$ galaxies spectroscopically identified with Ly$α$ emission, $2126$ of which host type 1 AGN showing broad (FWHM$~>1000$ km s$^{-1}$) Ly$α$ emission lines. We derive the Ly$α$ (UV) LF over 2 orders of magnitude covering bright galaxies and AGN in $\log L_\mathrm{Lyα}/\mathrm{[erg~s^{-1}]}=43.3-45.5$ ($-27<M_\mathrm{UV}<-20$) by the $1/V_\mathrm{max}$ estimator. Our results reveal the bright-end hump of the Ly$α$ LF is composed of type 1 AGN. In conjunction with previous spectroscopic results at the faint end, we measure a slope of the best-fit Schechter function to be $α_\mathrm{Sch}=-1.70^{+0.13}_{-0.14}$, which indicates $α_\mathrm{Sch}$ steepens from $z=2-3$ towards high redshift. Our UV LF agrees well with previous AGN UV LFs, and extends to faint-AGN and bright-galaxy regimes. The number fraction of Ly$α$-emitting objects ($X_\mathrm{LAE}$) increases from $M_\mathrm{UV}^*\sim-21$ to bright magnitude due to the contribution of type 1 AGN, while previous studies claim that $X_\mathrm{Lyα}$ decreases from faint magnitude to $M_\mathrm{UV}^*$, suggesting a valley in the $X_\mathrm{Lyα}-$magnitude relation at $M_\mathrm{UV}^*$. Comparing our UV LF of type 1 AGN at $z=2-3$ with those at $z=0$, we find that the number density of faint ($M_\mathrm{UV}>-21$) type 1 AGN increases from $z\sim2$ to $z\sim0$ as opposed to the evolution of bright ($M_\mathrm{UV}<-21$) type 1 AGN, suggesting the AGN downsizing in the rest-frame UV luminosity.

preprint2021arXiv

Mapping Obscuration to Reionization with ALMA (MORA): 2mm Efficiently Selects the Highest-Redshift Obscured Galaxies

We present the characteristics of 2mm-selected sources from the largest Atacama Large Millimeter and submillimeter Array (ALMA) blank-field contiguous survey conducted to-date, the Mapping Obscuration to Reionization with ALMA (MORA) survey covering 184arcmin$^2$ at 2mm. Twelve of the thirteen detections above 5$σ$ are attributed to emission from galaxies, eleven of which are dominated by cold dust emission. These sources have a median redshift of $\langle z_{\rm 2mm}\rangle=3.6^{+0.4}_{-0.3}$ primarily based on optical/near-infrared (OIR) photometric redshifts with some spectroscopic redshifts, with 77$\pm$11% of sources at $z>3$ and 38$\pm$12% of sources at $z>4$. This implies that 2mm selection is an efficient method for identifying the highest redshift dusty star-forming galaxies (DSFGs). Lower redshift DSFGs ($z<3$) are far more numerous than those at $z>3$ yet likely to drop out at 2mm. MORA shows that DSFGs with star-formation rates in excess of 300M$_\odot$ yr$^{-1}$ and relative rarity of $\sim$10$^{-5}$ Mpc$^{-3}$ contribute $\sim$30% to the integrated star-formation rate density between $3<z<6$. The volume density of 2mm-selected DSFGs is consistent with predictions from some cosmological simulations and is similar to the volume density of their hypothesized descendants: massive, quiescent galaxies at $z>2$. Analysis of MORA sources' spectral energy distributions hint at steeper empirically-measured dust emissivity indices than typical literature studies, with $\langleβ\rangle=2.2^{+0.5}_{-0.4}$. The MORA survey represents an important step in taking census of obscured star-formation in the Universe's first few billion years, but larger area 2mm surveys are needed to more fully characterize this rare population and push to the detection of the Universe's first dusty galaxies.

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

A comprehensive study of H$α$ emitters at $z \sim$ 0.62 in the DAWN survey: the need for deep and wide regions

We present new estimates of the luminosity function (LF) and star formation rate density (SFRD) for an H$α$ selected sample at $z\sim0.62$ from the Deep And Wide Narrow-band (DAWN) survey. Our results are based on a new H$α$ sample in the extended COSMOS region (compared to Coughlin et al. 2018) with the inclusion of flanking fields, resulting in a total area coverage of $\sim$1.5 deg$^2$. A total of 241 H$α$ emitters were selected based on robust selection criteria using spectro-photometric redshifts and broadband color-color classification. We explore the effect of different dust correction prescriptions by calculating the LF and SFRD using a constant dust extinction correction, A{$_{\textrm{H}α}=1$} mag, a luminosity-dependent correction, and a stellar-mass dependent correction. The resulting H$α$ LFs are well fitted using Schechter functions with best-fit parameters: L$^*=10^{42.24}$ erg s$^{-1}$, $ϕ^*=10^{-2.85}$ Mpc$^{-3}$, $α= -1.62$ for constant dust correction, L$^*=10^{42.31}$ erg s$^{-1}$, $ϕ^*=10^{-2.8}$ Mpc$^{-3}$, $α=-1.39$ for luminosity-dependent dust correction, and L$^*=10^{42.36}$ erg s$^{-1}$, $ϕ^*=10^{-2.91}$ Mpc$^{-3}$, $α= -1.48$, for stellar mass-dependent dust correction. The deep and wide nature of the DAWN survey effectively samples H$α$ emitters over a wide range of luminosities, thereby providing better constraints on both the faint and bright end of the LF. Also, the SFRD estimates $ρ_{\textrm{SFR}}=10^{-1.39}$ M$_{\odot}$yr$^{-1}$Mpc$^{-3}$ (constant dust correction), $ρ_{\textrm{SFR}}=10^{-1.47}$ M$_{\odot}$yr$^{-1}$Mpc$^{-3}$ (luminosity-dependent dust correction), and $ρ_{\textrm{SFR}}=10^{-1.49}$ M$_{\odot}$yr$^{-1}$Mpc$^{-3}$ (stellar mass-dependent dust correction) are in good agreement with the evolution of SFRD across redshifts ($0 < z < 2$) seen from previous H$α$ surveys.

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

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

Probing the Bright End of the Rest-Frame Ultraviolet Luminosity Function at z = 8-10 with Hubble Pure-Parallel Imaging

Looking for bright galaxies born in the early universe is fundamental to investigating the Epoch of Reionization, the era when the first stars and galaxies ionized the intergalactic medium. We utilize Hubble Space Telescope pure parallel imaging to select galaxy candidates at a time 500 to 650 million years after the Big Bang, which corresponds to redshifts z ~ 8-10. These data come from the Brightest of Reionizing Galaxies Survey (BoRG) Cycle 22 dataset, which consists of pure-parallel imaging in ~ 90 different lines of sight that sum up to an area of ~ 420 arcmin^2. This survey uses five filters and has the advantage (compared to the Cycle 21 BoRG program) of including imaging in the JH140 band, covering continuous wavelengths from the visible to near-infrared (lambda = 0.35um - 1.7um). This allows us to perform reliable selection of galaxies at z>8 using the photometric redshift technique. We use these galaxy candidates to constrain the bright end of the rest-frame ultraviolet luminosity function in this epoch. These candidates are excellent targets for follow-up observations, particularly with the James Webb Space Telescope.

preprint2020arXiv

Semi-analytic forecasts for JWST -- III. Intrinsic production efficiency of Lyman-continuum radiation

The James Webb Space Telescope (JWST) is expected to enable transformational progress in studying galaxy populations in the very early Universe, during the Epoch of Reionization (EoR). A critical parameter for understanding the sources that reionized the Universe is the Lyman-continuum production efficiency, $ξ_{\rm ion}$, defined as the rate of production of ionizing photons divided by the intrinsic UV luminosity. In this work, we combine self-consistent star formation and chemical enrichment histories predicted by semi-analytic models of galaxy formation with stellar population synthesis (SPS) models to predict the expected dependence of $ξ_{\rm ion}$ on galaxy properties and cosmic epoch from $z=4-10$. We then explore the sensitivity of the production rate of ionizing photons, $\dot{N}_{\rm ion}$, to the choice of SPS model and the treatment of stellar feedback in our galaxy formation model. We compare our results to those of other simulations, constraints from empirical models, and observations. We find that adopting SPS models that include binary stars predict about a factor of two more ionizing radiation than models that only assume single stellar populations. We find that UV-faint, low-mass galaxies have values of $ξ_{\rm ion}$ about 0.25 dex higher than those of more massive galaxies, but find weak evolution with cosmic time, about 0.2 dex from $z\sim 12-4$ at fixed rest-UV luminosity. We provide predictions of $\dot{N}_{\rm ion}$ as a function of $M_{\rm h}$ and a number of other galaxy properties. All results presented in this work are available at https://www.simonsfoundation.org/semi-analytic-forecasts-for-jwst/.

preprint2020arXiv

Semi-analytic forecasts for JWST -- IV. Implications for cosmic reionization and LyC escape fraction

Galaxies forming in low-mass halos are thought to be primarily responsible for reionizing the Universe during the first billion years after the Big Bang. Yet, these halos are extremely inefficient at forming stars in the nearby Universe. In this work, we address this apparent tension, and ask whether a physically motivated model of galaxy formation that reproduces the observed abundance of faint galaxies in the nearby Universe is also consistent with available observational constraints on the reionization history. By interfacing the Santa Cruz semi-analytic model for galaxy formation with an analytic reionization model, we constructed a computationally efficient pipeline that connects `ground-level' galaxy formation physics to `top-level' cosmological-scale observables. Based on photometric properties of the galaxy populations predicted up to $z=15$, we compute the reionization history of intergalactic hydrogen. We quantify the three degenerate quantities that influence the total ionizing photon budget, including the abundance of galaxies, the intrinsic production rate of ionizing photons, and the LyC escape fraction. We explore covariances between these quantities using a Markov chain Monte Carlo method. We find that our locally calibrated model is consistent with all currently available constraints on the reionization history, under reasonable assumptions about the LyC escape fraction. We quantify the fraction of ionizing photons contributed by galaxies of different luminosities and find that the galaxies expected to be detected in JWST NIRCam wide and deep surveys are responsible for producing $\sim 40$-$80\%$ of ionizing photons throughout the EoR. All results presented in this work are available at https://www.simonsfoundation.org/semi-analytic-forecasts-for-jwst/.

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.

preprint2019arXiv

A Catalog of Emission-Line Galaxies from the Faint Infrared Grism Survey: Studying Environmental Influence on Star Formation

We present a catalog of 208 $0.3 < z < 2.1$ Emission Line Galaxies (ELG) selected from 1D slitless spectroscopy obtained using Hubble's WFC3 G102 grism, as part of the Faint Infrared Grism Survey (FIGS). We identify ELG candidates by searching for significant peaks in all continuum-subtracted G102 spectra, and, where possible, confirm candidates by identifying consistent emission lines in other available spectra or with published spectroscopic redshifts. We provide derived emission line fluxes and errors, redshifts, and equivalent widths (EW) for H$α$ $\lambda6563$, [OIII]$λ\lambda4959,5007$, and [OII]$λ\lambda3727$ emission lines, for emission line galaxies down to AB(F105W) $ > 28$ and $> 10^{-17}$ erg cm$^{-2}$ s$^{-1}$ line flux. We use the resulting line catalog to investigate a possible relationship between line emission and a galaxy's environment. We use 7th-nearest-neighbor distances to investigate the typical surroundings of ELGs compared to non-ELGs, and we find that [OIII] emitters are preferentially found at intermediate galaxy densities near galaxy groups. We characterize these ELGs in terms of the galaxy specific star formation rate (SSFR) versus stellar mass, and find no significant influence of environment on that relation. We calculate star formation rates (SFR), and find no dependence of SFR on local galaxy surface density for $0.3 < z < 0.8$ H$α$ emitters and for $0.8<z<1.3$ [OIII] emitters. We find similar rates of close-pair interaction between ELGs and non-ELGs. For galaxy surface densities $Σ\leq 30$ Mpc$^{-2}$, we find no consistent effect of environment on star formation.

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

FIGS: Spectral fitting constraints on the star formation history of massive galaxies since Cosmic Noon

We constrain the stellar population properties of a sample of 52 massive galaxies, with stellar mass log Ms>10.5, over the redshift range 0.5<z<2 by use of observer-frame optical and near-infrared slitless spectra from HST's ACS and WFC3 grisms. The deep exposures (~100 ks) allow us to target individual spectra of massive galaxies to F160W=22.5AB. Our spectral fitting approach uses a set of six base models adapted to the redshift and spectral resolution of each observation, and fits the weights of the base models, including potential dust attenuation, via an MCMC method. Our sample comprises a mixed distribution of quiescent (19) and star-forming galaxies (33). We quantify the width of the age distribution (Dt) that is found to dominate the variance of the retrieved parameters according to Principal Component Analysis. The population parameters follow the expected trend towards older ages with increasing mass, and Dt appears to weakly anti-correlate with stellar mass, suggesting a more efficient star formation at the massive end. As expected, the redshift dependence of the relative stellar age (measured in units of the age of the Universe at the source) in the quiescent sample rejects the hypothesis of a single burst (aka monolithic collapse). Radial colour gradients within each galaxy are also explored, finding a wider scatter in the star-forming subsample, but no conclusive trend with respect to the population parameters.

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.

preprint2019arXiv

Physical Characterization of an Unlensed Dusty Star-Forming Galaxy at $z=5.85$

We present a physical characterization of MMJ100026.36+021527.9 (a.k.a. ``MAMBO-9''), a dusty star-forming galaxy (DSFG) at $z=5.850\pm0.001$. This is the highest redshift unlensed DSFG (and fourth most distant overall) found to-date, and is the first source identified in a new 2mm blank-field map in the COSMOS field. Though identified in prior samples of DSFGs at 850$μ$m-1.2mm with unknown redshift, the detection at 2mm prompted further follow-up as it indicated a much higher probability that the source was likely to sit at $z>4$. Deep observations from the Atacama Large Millimeter and submillimeter Array (ALMA) presented here confirm the redshift through the secure detection of $^{12}$CO($J\!=$6$\rightarrow$5) and p-H$_{2}$O(2$_{1,1}\!\rightarrow$2$_{0,2}$). MAMBO-9 is comprised of a pair of galaxies separated by 6kpc with corresponding star-formation rates of 590M$_\odot$yr$^{-1}$ and 220M$_\odot$yr$^{-1}$ total molecular hydrogen gas mass of (1.7$\pm$0.4)$\times10^{11}$M$_\odot$, dust mass of (1.3$\pm$0.3)$\times10^{9}$M$_\odot$ and stellar mass of (3.2$^{+1.0}_{-1.5}$)$\times10^{9}$M$_\odot$. The total halo mass, (3.3$\pm$0.8)$\times10^{12}$M$_\odot$, is predicted to exceed $>10^{15}$M$_\odot$ by $z=0$. The system is undergoing a merger-driven starburst which will increase the stellar mass of the system tenfold in $τ_{\rm depl}=40-80$Myr, converting its large molecular gas reservoir (gas fraction of 96$^{+1}_{-2}$%) into stars. MAMBO-9 evaded firm spectroscopic identification for a decade, following a pattern that has emerged for some of the highest redshift DSFGs found. And yet, the systematic identification of unlensed DSFGs like MAMBO-9 is key to measuring the global contribution of obscured star-formation to the star-formation rate density at $z>4$, the formation of the first massive galaxies, and the formation of interstellar dust at early times ($<$1Gyr).

preprint2017arXiv

Discovery of a z = 7.452 High Equivalent Width Lyman-α Emitter from the Hubble Space Telescope Faint Infrared Grism Survey

We present the results of an unbiased search for Lyα emission from continuum-selected 6 < z < 8 galaxies. Our dataset consists of 160 orbits of G102 slitless grism spectroscopy obtained with the Hubble Space Telescope (HST) Wide Field Camera 3 as part of the Faint Infrared Grism Survey (FIGS; PI: Malhotra), which obtains deep slitless spectra of all sources in four fields, and was designed to minimize contamination in observations of previously-identified high-redshift galaxy candidates. The FIGS data can potentially spectroscopically confirm the redshifts of galaxies, and as Lyα emission is resonantly scattered by neutral gas, FIGS can also constrain the ionization state of the intergalactic medium (IGM) during the epoch of reionization. These data have sufficient depth to detect Lyα emission in this epoch, as Tilvi et al. (2016) have published the FIGS detection of previously known (Finkelstein et al. 2013) Lyα emission at z = 7.51. The FIGS data use five separate roll-angles of HST to mitigate the contamination by nearby galaxies. We created a method that accounts for and removes the contamination from surrounding galaxies, and also removes any dispersed continuum light from each individual spectrum (Pirzkal et al. 2017). We searched for significant (> 4σ) emission lines using two different automated detection methods, free of any visual inspection biases. Applying these methods on photometrically-selected high-redshift candidates between 6 < z < 8 we find two emission lines, one previously published by Tilvi et al. (2016), and a new line at 1.028μm. We identify this lines as Lyα at z = 7.452 +/- 0.003. This newly spectroscopically confirmed galaxy has the highest Lyα rest-frame equivalent width (EW Lyα) yet published at z > 7 (140.3 +/- 19.0Å).

preprint2016arXiv

Connecting the Dots: Tracking Galaxy Evolution Using Constant Cumulative Number Density at 3<z<7

Using the cosmological smoothed particle hydrodynamical code GADGET-3 we make a realistic assessment of the technique of using constant cumulative number density as a tracer of galaxy evolution at high redshift. We find that over a redshift range of $3\leq z \leq7$ one can on average track the growth of the stellar mass of a population of galaxies selected from the same cumulative number density bin to within $\sim 0.20$ dex. Over the stellar mass range we probe ($10^{10.39}\leq M_s/M_\odot \leq 10^{10.75}$ at $z =$ 3 and $10^{8.48}\leq M_s/M_\odot \leq 10^{9.55}$ at $z =$ 7) one can reduce this bias by selecting galaxies based on an evolving cumulative number density. We find the cumulative number density evolution exhibits a trend towards higher values which can be quantified by simple linear formulations going as $-0.10Δz$ for descendants and $0.12Δz$ for progenitors. Utilizing such an evolving cumulative number density increases the accuracy of descendant/progenitor tracking by a factor of $\sim2$. This result is in excellent agreement, within $0.10$ dex, with abundance matching results over the same redshift range. However, we find that our more realistic cosmological hydrodynamic simulations produce a much larger scatter in descendant/progenitor stellar masses than previous studies, particularly when tracking progenitors. This large scatter makes the application of either the constant cumulative number density or evolving cumulative number density technique limited to average stellar masses of populations only, as the diverse mass assembly histories caused by stochastic physical processes such as gas accretion, mergers, and star formation of individual galaxies will lead to a larger scatter in other physical properties such as metallicity and star-formation rate.

preprint2016arXiv

How Lyman Alpha Emission Depends On Galaxy Stellar Mass

In this work, we show how the stellar mass (M) of galaxies affects the 3<z<4.6 Ly-alpha equivalent width (EW) distribution. To this end, we design a sample of 629 galaxies in the M range 7.6 < logM/Msun < 10.6 from the 3D-HST/CANDELS survey. We perform spectroscopic observations of this sample using the Michigan/Magellan Fiber System, allowing us to measure Ly-alpha fluxes and use 3D-HST/CANDELS ancillary data. In order to study the Ly-alpha EW distribution dependence on M, we split the whole sample in three stellar mass bins. We find that, in all bins, the distribution is best represented by an exponential profile of the form dN(M)/dEW= A(M)exp(-EW/W0(M))/W0(M). Through a Bayesian analysis, we confirm that lower M galaxies have higher Ly-alpha EWs. We also find that the fraction A of galaxies featuring emission and the e-folding scale W0 of the distribution anti- correlate with M, recovering expressions of the forms A(M)= -0.26(.13) logM/Msun+3.01(1.2) and W0(M)= -15.6(3.5) logM/Msun +166(34). These results are crucial for proper interpretation of Ly-alpha emission trends reported in the literature that may be affected by strong M selection biases.

preprint2016arXiv

Keck/MOSFIRE Spectroscopy of z=7-8 Galaxies: Ly$α$ Emission from a Galaxy at z=7.66

We report the results from some of the deepest Keck/Multi-Object Spectrometer For Infra-Red Exploration data yet obtained for candidate $z \gtrsim 7$ galaxies. Our data show one significant line detection with 6.5$σ$ significance in our combined 10 hr of integration which is independently detected on more than one night, thus ruling out the possibility that the detection is spurious. The asymmetric line profile and non-detection in the optical bands strongly imply that the detected line is Ly$α$ emission from a galaxy at $z$(Ly$α)=7.6637 \pm 0.0011$, making it the fourth spectroscopically confirmed galaxy via Ly$α$ at $z>7.5$. This galaxy is bright in the rest-frame ultraviolet (UV; $M_{\rm UV} \sim -21.2$) with a moderately blue UV slope ($β=-2.2^{+0.3}_{-0.2}$), and exhibits a rest-frame Ly$α$ equivalent width of EW(Ly$α$) $\sim 15.6^{+5.9}_{-3.6}$ Å. The non-detection of the 11 other $z \sim$ 7-8 galaxies in our long 10 hr integration, reaching a median 5$σ$ sensitivity of 28 Å in the rest-frame EW(Ly$α$), implies a 1.3$σ$ deviation from the null hypothesis of a non-evolving distribution in the rest-frame EW(Ly$α$) between $3<z<6$ and $z=$ 7-8. Our results are consistent with previous studies finding a decline in Ly$α$ emission at $z>6.5$, which may signal the evolving neutral fraction in the intergalactic medium at the end of the reionization epoch, although our weak evidence suggests the need for a larger statistical sample to allow for a more robust conclusion.

preprint2016arXiv

Lyman-Alpha Emitter Galaxies at z ~ 2.8 in the Extended Chandra Deep Field-South: I. Tracing the Large-Scale Structure via Lyman-Alpha Imaging

We present a narrowband survey with three adjacent filters for z=2.8--2.9 Lyman Alpha Emitter (LAE) galaxies in the Extended Chandra Deep Field South (ECDFS), along with spectroscopic followup. With a complete sample of 96 LAEs in the narrowband NB466, we confirm a large-scale structure at z~ 2.8. Compared to the blank field in NB470 and NB475, the LAE density excess in the NB466 field is ~6.0+/-0.8 times the standard deviation expected at z~2.8, assuming a linear bias of 2. The overdense large scale structure in NB466 can be decomposed into 4 protoclusters, whose overdensities are 4.6 - 6.6. These 4 protoclusters are expected to evolve into a Coma-like cluster at z~ 0. In the meanwhile, we investigate the average star-formation rates derived from Lyα, rest-frame UV and X-ray, the Lyα luminosity functions, the Lyα photon densities and their dependence on the environment. We find that the Lyα photon density in the overdense field (NB466) is ~50\% higher than that in the blank field (NB470+NB475). The 3 brightest LAEs, including a quasar at z=2.81, are all detected in X-ray and in NB466. These three LAE-AGNs contribute an extra 20--30\% Lyα photon density. Furthermore, we find that LAEs in overdense regions are younger and less dusty. We conclude that the structure we found is a significant and rare density peak, and narrowband imaging is an efficient method to detect and study such structures in the high-z universe.

preprint2016arXiv

Observational Searches for Galaxies at z > 6

Although the universe at redshifts greater than six represents only the first one billion years (<10%) of cosmic time, the dense nature of the early universe led to vigorous galaxy formation and evolution activity which we are only now starting to piece together. Technological improvements have, over only the past decade, allowed large samples of galaxies at such high redshifts to be collected, providing a glimpse into the epoch of formation of the first stars and galaxies. A wide variety of observational techniques have led to the discovery of thousands of galaxy candidates at z > 6, with spectroscopically confirmed galaxies out to nearly z = 9. Using these large samples, we have begun to gain a physical insight into the processes inherent in galaxy evolution at early times. In this review, I will discuss i) the selection techniques for finding distant galaxies, including a summary of previous and ongoing ground and space-based searches, and spectroscopic followup efforts, ii) insights into galaxy evolution gleaned from measures such as the rest-frame ultraviolet luminosity function, the stellar mass function, and galaxy star-formation rates, and iii) the effect of galaxies on their surrounding environment, including the chemical enrichment of the universe, and the reionization of the intergalactic medium. Finally, I conclude with prospects for future observational study of the distant universe, using a bevy of new state-of-the-art facilities coming online over the next decade and beyond.

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.

preprint2015arXiv

A Critical Assessment of Stellar Mass Measurement Methods

In this paper we perform a comprehensive study of the main sources of random and systematic errors in stellar mass measurement for galaxies using their Spectral Energy Distributions (SEDs). We use mock galaxy catalogs with simulated multi-waveband photometry (from U-band to mid-infrared) and known redshift, stellar mass, age and extinction for individual galaxies. Given different parameters affecting stellar mass measurement (photometric S/N ratios, SED fitting errors, systematic effects, the inherent degeneracies and correlated errors), we formulated different simulated galaxy catalogs to quantify these effects individually. We studied the sensitivity of stellar mass estimates to the codes/methods used, population synthesis models, star formation histories, nebular emission line contributions, photometric uncertainties, extinction and age. For each simulated galaxy, the difference between the input stellar masses and those estimated using different simulation catalogs, $Δ\log(M)$, was calculated and used to identify the most fundamental parameters affecting stellar masses. We measured different components of the error budget, with the results listed as follows: (1). no significant bias was found among different codes/methods, with all having comparable scatter; (2). A source of error is found to be due to photometric uncertainties and low resolution in age and extinction grids; (3). The median of stellar masses among different methods provides a stable measure of the mass associated with any given galaxy; (4). The deviations in stellar mass strongly correlate with those in age, with a weaker correlation with extinction; (5). the scatter in the stellar masses due to free parameters are quantified, with the sensitivity of the stellar mass to both the population synthesis codes and inclusion of nebular emission lines studied.

preprint2015arXiv

A WFC3 Grism Emission Line Redshift Catalog in the GOODS-South Field

We combine HST/WFC3 imaging and G141 grism observations from the CANDELS and 3D-HST surveys to produce a catalog of grism spectroscopic redshifts for galaxies in the CANDELS/GOODS-South field. The WFC3/G141 grism spectra cover a wavelength range of 1.1<lambda<1.7 microns with a resolving power of R~130 for point sources, thus providing rest-frame optical spectra for galaxies out to z~3.5. The catalog is selected in the H-band (F160W) and includes both galaxies with and without previously published spectroscopic redshifts. Grism spectra are extracted for all H-band detected galaxies with H<24 and a CANDELS photometric redshift z_phot > 0.6. The resulting spectra are visually inspected to identify emission lines and redshifts are determined using cross-correlation with empirical spectral templates. To establish the accuracy of our redshifts, we compare our results against high-quality spectroscopic redshifts from the literature. Using a sample of 411 control galaxies, this analysis yields a precision of sigma_NMAD=0.0028 for the grism-derived redshifts, which is consistent with the accuracy reported by the 3D-HST team. Our final catalog covers an area of 153 square arcmin and contains 1019 redshifts for galaxies in GOODS-S. Roughly 60% (608/1019) of these redshifts are for galaxies with no previously published spectroscopic redshift. These new redshifts span a range of 0.677 < z < 3.456 and have a median redshift of z=1.282. The catalog contains a total of 234 new redshifts for galaxies at z>1.5. In addition, we present 20 galaxy pair candidates identified for the first time using the grism redshifts in our catalog, including four new galaxy pairs at z~2, nearly doubling the number of such pairs previously identified.

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

HST Emission Line Galaxies at z ~ 2: Comparing Physical Properties of Lyman Alpha and Optical Emission Line Selected Galaxies

We compare the physical and morphological properties of z ~ 2 Lyman-alpha emitting galaxies (LAEs) identified in the HETDEX Pilot Survey and narrow band studies with those of z ~ 2 optical emission line galaxies (oELGs) identified via HST WFC3 infrared grism spectroscopy. Both sets of galaxies extend over the same range in stellar mass (7.5 < logM < 10.5), size (0.5 < R < 3.0 kpc), and star-formation rate (~1 < SFR < 100). Remarkably, a comparison of the most commonly used physical and morphological parameters -- stellar mass, half-light radius, UV slope, star formation rate, ellipticity, nearest neighbor distance, star formation surface density, specific star formation rate, [O III] luminosity, and [O III] equivalent width -- reveals no statistically significant differences between the populations. This suggests that the processes and conditions which regulate the escape of Ly-alpha from a z ~ 2 star-forming galaxy do not depend on these quantities. In particular, the lack of dependence on the UV slope suggests that Ly-alpha emission is not being significantly modulated by diffuse dust in the interstellar medium. We develop a simple model of Ly-alpha emission that connects LAEs to all high-redshift star forming galaxies where the escape of Ly-alpha depends on the sightline through the galaxy. Using this model, we find that mean solid angle for Ly-alpha escape is 2.4+/-0.8 steradians; this value is consistent with those calculated from other studies.

preprint2015arXiv

Probing the Physical Properties of z=4.5 Lyman Alpha Emitters with Spitzer

We present the results from a stellar population modeling analysis of a sample of 162 z=4.5, and 14 z=5.7 Lyman alpha emitting galaxies (LAEs) in the Bootes field, using deep Spitzer/IRAC data at 3.6 and 4.5 um from the Spitzer Lyman Alpha Survey, along with Hubble Space Telescope NICMOS and WFC3 imaging at 1.1 and 1.6 um for a subset of the LAEs. This represents one of the largest samples of high-redshift LAEs imaged with Spitzer IRAC. We find that 30/162 (19%) of the z=4.5 LAEs and 9/14 (64%) of the z=5.7 LAEs are detected at >3-sigma in at least one IRAC band. Individual z=4.5 IRAC-detected LAEs have a large range of stellar mass, from 5x10^8 to 10^11 Msol. One-third of the IRAC-detected LAEs have older stellar population ages of 100 Myr - 1 Gyr, while the remainder have ages < 100 Myr. A stacking analysis of IRAC-undetected LAEs shows this population to be primarily low mass (8 -- 20 x 10^8 Msol) and young (64 - 570 Myr). We find a correlation between stellar mass and the dust-corrected ultraviolet-based star-formation rate (SFR) similar to that at lower redshifts, in that higher mass galaxies exhibit higher SFRs. However, the z=4.5 LAE correlation is elevated 4-5 times in SFR compared to continuum-selected galaxies at similar redshifts. The exception is the most massive LAEs which have SFRs similar to galaxies at lower redshifts suggesting that they may represent a different population of galaxies than the traditional lower-mass LAEs, perhaps with a different mechanism promoting Lyman alpha photon escape.

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

Surveying Galaxy Proto-clusters in Emission: A Large-scale Structure at z=2.44 and the Outlook for HETDEX

Galaxy proto-clusters at z >~ 2 provide a direct probe of the rapid mass assembly and galaxy growth of present day massive clusters. Because of the need of precise galaxy redshifts for density mapping and the prevalence of star formation before quenching, nearly all the proto-clusters known to date were confirmed by spectroscopy of galaxies with strong emission lines. Therefore, large emission-line galaxy surveys provide an efficient way to identify proto-clusters directly. Here we report the discovery of a large-scale structure at z = 2.44 in the HETDEX Pilot Survey. On a scale of a few tens of Mpc comoving, this structure shows a complex overdensity of Lya emitters (LAE), which coincides with broad-band selected galaxies in the COSMOS/UltraVISTA photometric and zCOSMOS spectroscopic catalogs, as well as overdensities of intergalactic gas revealed in the Lya absorption maps of Lee et al. (2014). We construct mock LAE catalogs to predict the cosmic evolution of this structure. We find that such an overdensity should have already broken away from the Hubble flow, and part of the structure will collapse to form a galaxy cluster with 10^14.5 +- 0.4 M_sun by z = 0. The structure contains a higher median stellar mass of broad-band selected galaxies, a boost of extended Lya nebulae, and a marginal excess of active galactic nuclei relative to the field, supporting a scenario of accelerated galaxy evolution in cluster progenitors. Based on the correlation between galaxy overdensity and the z = 0 descendant halo mass calibrated in the simulation, we predict that several hundred 1.9 < z < 3.5 proto-clusters with z = 0 mass of > 10^14.5 M_sun will be discovered in the 8.5 Gpc^3 of space surveyed by the Hobby Eberly Telescope Dark Energy Experiment.

preprint2015arXiv

The Case for a James Webb Space Telescope Extragalactic Key Project

The upcoming launch of the James Webb Space Telescope (JWST) in less than three years is certain to bring a revolution in our understanding of many area of astrophysics, with one of the key goals being galaxy evolution. As the first proposals will be due in a little over two years, the time is ripe to take a holistic look at the science goals which the community would wish to accomplish with this observatory. Contrary to our experiences with the Hubble Space Telescope, which has now operated successfully for over two decades due to several timely servicing missions, the lifetime of JWST is finite and relatively short, with a lifetime requirement of five years, and a ten-year goal. Following the discussion session at the "Exploring the Universe with JWST" conference at ESA-ESTEC in October 2015, we highlight in this document the (non-local) extragalactic science goals for JWST. We describe how a concerted community effort could best address these, ensuring that the desired survey can be completed during the JWST mission.

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

UVUDF: Ultraviolet Through Near-infrared Catalog and Photometric Redshifts of Galaxies in the Hubble Ultra Deep Field

We present photometry and derived redshifts from up to eleven bandpasses for 9927 galaxies in the Hubble Ultra Deep field (UDF), covering an observed wavelength range from the near-ultraviolet (NUV) to the near-infrared (NIR) with Hubble Space Telescope observations. Our Wide Field Camera 3 (WFC3)/UV F225W, F275W, and F336W image mosaics from the ultra-violet UDF (UVUDF) imaging campaign are newly calibrated to correct for charge transfer inefficiency, and use new dark calibrations to minimize background gradients and pattern noise. Our NIR WFC3/IR image mosaics combine the imaging from the UDF09 and UDF12 campaigns with CANDELS data to provide NIR coverage for the entire UDF field of view. We use aperture-matched point-spread function corrected photometry to measure photometric redshifts in the UDF, sampling both the Lyman break and Balmer break of galaxies at z~0.8-3.4, and one of the breaks over the rest of the redshift range. Our comparison of these results with a compilation of robust spectroscopic redshifts shows an improvement in the galaxy photometric redshifts by a factor of two in scatter and a factor three in outlier fraction over previous UDF catalogs. The inclusion of the new NUV data is responsible for a factor of two decrease in the outlier fraction compared to redshifts determined from only the optical and NIR data, and improves the scatter at z<0.5 and at z>2. The panchromatic coverage of the UDF from the NUV through the NIR yields robust photometric redshifts of the UDF, with the lowest outlier fraction available.

preprint2014arXiv

Galactic winds and stellar populations in Lyman $α$ emitting galaxies at z ~ 3.1

We present a sample of 33 spectroscopically confirmed z ~ 3.1 Ly$α$-emitting galaxies (LAEs) in the Cosmological Evolution Survey (COSMOS) field. This paper details the narrow-band survey we conducted to detect the LAE sample, the optical spectroscopy we performed to confirm the nature of these LAEs, and a new near-infrared spectroscopic detection of the [O III] 5007 Å line in one of these LAEs. This detection is in addition to two [O III] detections in two z ~ 3.1 LAEs we have reported on previously (McLinden et al 2011). The bulk of the paper then presents detailed constraints on the physical characteristics of the entire LAE sample from spectral energy distribution (SED) fitting. These characteristics include mass, age, star-formation history, dust content, and metallicity. We also detail an approach to account for nebular emission lines in the SED fitting process - wherein our models predict the strength of the [O III] line in an LAE spectrum. We are able to study the success of this prediction because we can compare the model predictions to our actual near-infrared observations both in galaxies that have [O III] detections and those that yielded non-detections. We find a median stellar mass of 6.9 $\times$ 10$^8$ M$_{\odot}$ and a median star formation rate weighted stellar population age of 4.5 $\times$ 10$^6$ yr. In addition to SED fitting, we quantify the velocity offset between the [O III] and Ly$α$ lines in the galaxy with the new [O III] detection, finding that the Ly$α$ line is shifted 52 km s$^{-1}$ redwards of the [O III] line, which defines the systemic velocity of the galaxy.

preprint2014arXiv

Galaxy Zoo: CANDELS Barred Disks and Bar Fractions

The formation of bars in disk galaxies is a tracer of the dynamical maturity of the population. Previous studies have found that the incidence of bars in disks decreases from the local Universe to z ~ 1, and by z > 1 simulations predict that bar features in dynamically mature disks should be extremely rare. Here we report the discovery of strong barred structures in massive disk galaxies at z ~ 1.5 in deep rest-frame optical images from CANDELS. From within a sample of 876 disk galaxies identified by visual classification in Galaxy Zoo, we identify 123 barred galaxies. Selecting a sub-sample within the same region of the evolving galaxy luminosity function (brighter than L*), we find that the bar fraction across the redshift range 0.5< z < 2 (f_bar = 10.7 +6.3 -3.5% after correcting for incompleteness) does not significantly evolve. We discuss the implications of this discovery in the context of existing simulations and our current understanding of the way disk galaxies have evolved over the last 11 billion years.

preprint2014arXiv

Lyα Equivalent Width Distribution of Lyα Emitting Galaxies at Redshift z $\sim$ 4.5

Lyα line EWs provide important clues to the physical nature of high redshift LAEs. However, measuring the Lyα EW distribution of high-z narrowband selected LAEs can be hard because many sources do not have broadband photometry. We investigate the possible biases in measuring the intrinsic Lyα EW distribution for a LAE sample at z $\sim$ 4.5 in the Extended CDFS. Only weak Malmquist-type bias in both the intrinsic Lyα luminosity function and the Lyα EW distribution were found. However, the observed EW distribution is severely biased if one only considers LAEs with detections in the continuum. Taking the broadband non-detections into account requires fitting the distribution of the broadband-to-narrowband ratio, which then gives a larger EW distribution scale length. Assuming an exponential form of the intrinsic Lyα EW distribution dN/dEW = N exp(-EW/W$_0$), we obtain W$_0$ = 167+/-44Å (uncorrected for IGM absorption of Lyα). We discuss the likely range of IGM absorption effects in light of recent measurements of Lyα line profiles and velocity offsets. Our data are consistent with Lyα EW being independent of UV luminosity (i.e., we do not see evidence for the "Ando" effect). Our simulations also imply that broad-band images should be 0.5-1 magnitude deeper than narrowband images for an effective and reasonably complete LAE survey. Comparing with consistent measurements at other redshifts, we see a strong evolution in Lyα EW distribution with redshift which goes as a power-law form of W$_0$ $\prop$(1+z)$^$ξ, with ξ = 1.1$\pm$0.1 (0.6$\pm$0.1) if no IGM corrections are applied to the Lyα line; or ξ = 1.7$\pm$0.1 (1.2$\pm$0.1) after applying a maximal IGM-absorption correction to Lyα line) for an exponential (a gaussian) EW distribution from z = 0.3 to 6.5.

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

Spectral Energy Distribution Fitting of HETDEX Pilot Survey Lyman-alpha Emitters in COSMOS and GOODS-N

We use broadband photometry extending from the rest-frame UV to the near-IR to fit the individual spectral energy distributions (SEDs) of 63 bright (L(Ly-alpha) > 10^43 ergs/s) Ly-alpha emitting galaxies (LAEs) in the redshift range 1.9 < z < 3.6. We find that these LAEs are quite heterogeneous, with stellar masses that span over three orders of magnitude, from 7.5 < log M < 10.5. Moreover, although most LAEs have small amounts of extinction, some high-mass objects have stellar reddenings as large as E(B-V) ~0.4. Interestingly, in dusty objects the optical depths for Ly-alpha and the UV continuum are always similar, indicating that Ly-alpha photons are not undergoing many scatters before escaping their galaxy. In contrast, the ratio of optical depths in low-reddening systems can vary widely, illustrating the diverse nature of the systems. Finally, we show that in the star formation rate (SFR)-log mass diagram, our LAEs fall above the "main-sequence" defined by z ~ 3 continuum selected star-forming galaxies. In this respect, they are similar to sub-mm-selected galaxies, although most LAEs have much lower mass.

preprint2014arXiv

The HETDEX Pilot Survey V: The Physical Origin of Lyman-alpha Emitters Probed by Near-infrared Spectroscopy

We present the results from a VLT/SINFONI and Keck/NIRSPEC near-infrared spectroscopic survey of 16 Lyman-alpha emitters (LAEs) at $z$ = 2.1 - 2.5 in the COSMOS and GOODS-N fields discovered from the HETDEX Pilot Survey. We detect rest-frame optical nebular lines (H$α$ and/or [OIII]$λ$5007) for 10 of the LAEs and measure physical properties, including the star formation rate (SFR), gas-phase metallicity, gas-mass fraction, and Ly$α$ velocity offset. We find that LAEs may lie below the mass-metallicity relation for continuum-selected star-forming galaxies at the same redshift. The LAEs all show velocity shifts of Ly$α$ relative to the systemic redshift ranging between +85 and +296 km s$^{-1}$ with a mean of +180 km s$^{-1}$. This value is smaller than measured for continuum-selected star-forming galaxies at similar redshifts. The Ly$α$ velocity offsets show a moderate correlation with the measured star formation rate (2.5$σ$), but no significant correlations are seen with the SFR surface density, specific SFR, stellar mass, or dynamical mass ($\lesssim$ 1.5$σ$). Exploring the role of dust, kinematics of the interstellar medium (ISM), and geometry on the escape of Ly$α$ photons, we find no signature of selective quenching of resonantly scattered Ly$α$ photons. However, we also find no evidence that a clumpy ISM is enhancing the Ly$α$ equivalent width. Our results suggest that the low metallicity in LAEs may be responsible for yielding an environment with a low neutral hydrogen column density as well as less dust, easing the escape of Ly$α$ photons over that in continuum-selected star-forming galaxies.

preprint2014arXiv

The interstellar medium and feedback in the progenitors of the compact passive galaxies at z~2

Quenched galaxies at z>2 are nearly all very compact relative to z~0, suggesting a physical connection between high stellar density and efficient, rapid cessation of star-formation. We present restframe UV spectra of Lyman-break galaxies (LBGs) at z~3 selected to be candidate progenitors of quenched galaxies at z~2 based on their compact restframe optical sizes and high surface density of star-formation. We compare their UV properties to those of more extended LBGs of similar mass and star formation rate (non-candidates). We find that candidate progenitors have faster ISM gas velocities and higher equivalent widths of interstellar absorption lines, implying larger velocity spread among absorbing clouds. Candidates deviate from the relationship between equivalent widths of Lyman-alpha and interstellar absorption lines in that their Lyman-alpha emission remains strong despite high interstellar absorption, possibly indicating that the neutral HI fraction is patchy such that Lyman-alpha photons can escape. We detect stronger CIV P-Cygni features (emission and absorption) and HeII emission in candidates, indicative of larger populations of metal rich Wolf-Rayet stars compared to non-candidates. The faster bulk motions, broader spread of gas velocity, and Lyman-alpha properties of candidates are consistent with their ISM being subject to more energetic feedback than non-candidates. Together with their larger metallicity (implying more evolved star-formation activity) this leads us to propose, if speculatively, that they are likely to quench sooner than non-candidates, supporting the validity of selection criteria used to identify them as progenitors of z~2 passive galaxies. We propose that massive, compact galaxies undergo more rapid growth of stellar mass content, perhaps because the gas accretion mechanisms are different, and quench sooner than normally-sized LBGs at these early epochs.

preprint2014arXiv

The UV Continuum of z > 1 Star-forming Galaxies in the Hubble Ultraviolet UltraDeep Field

We estimate the UV continuum slope, beta, for 923 galaxies in the range 1 < z < 8 in the Hubble Ultradeep Field (HUDF). These data include 460 galaxies at 1 < z < 2 down to an absolute magnitude M_{UV} = -14 (~0.006 L*_{z=1}; 0.02 L*_{z=0}), comparable to dwarf galaxies in the local universe. We combine deep HST/UVIS photometry in F225W, F275W, F336W wavebands (UVUDF) with recent data from HST/WFC3/IR (HUDF12). Galaxies in the range 1 < z < 2 are significantly bluer than local dwarf galaxies. We find their mean (median) values <beta> = -1.382 (-1.830) +/- 0.002 (random) +/- 0.1 (systematic). We find comparable scatter in beta (standard deviation = 0.43) to local dwarf galaxies and 30% larger scatter than z > 2 galaxies. We study the trends of beta with redshift and absolute magnitude for binned sub-samples and find a modest color-magnitude relation, dbeta/dM = -0.11 +/- 0.01 and no evolution in dbeta/dM with redshift. A modest increase in dust reddening with redshift and luminosity, Delta E(B-V) ~ 0.1, and a comparable increase in the dispersion of dust reddening at z < 2, appears likely to explain the observed trends. At z > 2, we find trends that are consistent with previous works; combining our data with the literature in the range 1 < z < 8, we find a color evolution with redshift, dbeta/dz = -0.09 +/-0.01 for low luminosity (0.05 L*_{z=3}), and dbeta/dz = -0.06 +/-0.01 for medium luminosity (0.25 L*_{z=3}) galaxies.

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

Lyman alpha Luminosity Functions at Redshift z = 4.5

We present a spectroscopically confirmed sample of Lyman alpha emitting galaxies (LAEs) at z ~ 4.5 in the Extended Chandra Deep Field South (ECDFS), which we combine with a sample of z ~ 4.5 LAEs from the Large Area Lyman Alpha (LALA) survey to build a unified Lya luminosity function (LF). We spectroscopically observed 64 candidate LAEs in the ECDFS, confirming 46 objects as z~4.5 LAEs. We did not detect significant flux from neither C_iv 1549Å nor the He_ii 1640Å emission in individual LAE spectra, even with a coadded spectrum. With the coadded line ratio of He_ii to Lya constraining the Population III star formation rate (SFR) to be <0.3% of the total SFR, and <1.25% of the observed SFR (both at the 2-$σ$ level). Only one LAE was detected in both the X-ray and radio, while the other objects remained undetected, even when stacked. The Lya LF in our two deepest narrowband filters in the ECDFS differ at >2$σ$ significance, and the product $L^*Φ^*$ differs by a factor of >3. Similar LF differences have been used to infer evolution in the neutral gas fraction in the intergalactic medium at z>6, yet here the difference is likely due to cosmic variance, given that the two samples are from adjoining line-of-sight volumes. Combining our new sample of LAEs with those from previous LALA narrowband surveys at z = 4.5, we obtain one of the best measured Lya LFs to date of L* = 42.83 $\pm$ 0.06 and $Φ^*$ = -3.48 $\pm$ 0.09. We compare our new LF to others from the literature to study the evolution of the Lya luminosity density at 0 < z < 7. We find tentative evidence for evolution in the product $L^* Φ^*$, which approximately tracks the cosmic SFR density, but since field-to-field and survey-to-survey variations are in some cases as large as the possible evolution, some caution is needed in interpreting this trend.

preprint2013arXiv

The Dynamical Masses, Densities, and Star Formation Scaling Relations of Lyman Alpha Galaxies

We present the first dynamical mass measurements for Lyman alpha galaxies at high redshift, based on velocity dispersion measurements from rest-frame optical emission lines and size measurements from HST imaging, for a sample of nine galaxies drawn from four surveys. These measurements enable us to study the nature of Lyman alpha galaxies in the context of galaxy scaling relations. The resulting dynamical masses range from 1e9 to 1e10 solar masses. We also fit stellar population models to our sample, and use them to plot the Lyman alpha sample on a stellar mass vs. line width relation. Overall, the Lyman alpha galaxies follow well the scaling relation established by observing star forming galaxies at lower redshift (and without regard for Lyman alpha emission), though in 1/3 of the Lyman alpha galaxies, lower-mass fits are also acceptable. In all cases, the dynamical masses agree with established stellarmass-linewidth relation. Using the dynamical masses as an upper limit on gas mass, we show that Lyman alpha galaxies resemble starbursts (rather than "normal" galaxies) in the relation between gas mass surface density and star formation activity, in spite of relatively modest star formation rates. Finally, we examine the mass densities of these galaxies, and show that their future evolution likely requires dissipational ("wet") merging. In short, we find that Lyman alpha galaxies are low mass cousins of larger starbursts.

preprint2013arXiv

The HETDEX Pilot Survey. IV. The Evolution of [O II] Emitting Galaxies from z ~ 0.5 to z ~ 0

We present an analysis of the luminosities and equivalent widths of the 284 z < 0.56 [O II]-emitting galaxies found in the 169 square arcmin pilot survey for the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). By combining emission-line fluxes obtained from the Mitchell spectrograph on the McDonald 2.7-m telescope with deep broadband photometry from archival data, we derive each galaxy's de-reddened [O II] 3727 luminosity and calculate its total star formation rate. We show that over the last ~5 Gyr of cosmic time there has been substantial evolution in the [O II] emission-line luminosity function, with L* decreasing by ~0.6 +/-0.2 dex in the observed function, and by ~0.9 +/-0.2 dex in the de-reddened relation. Accompanying this decline is a significant shift in the distribution of [O II] equivalent widths, with the fraction of high equivalent-width emitters declining dramatically with time. Overall, the data imply that the relative intensity of star formation within galaxies has decreased over the past ~5 Gyr, and that the star formation rate density of the universe has declined by a factor of ~2.5 between z ~ 0.5 and z ~ 0. These observations represent the first [O II]-based star formation rate density measurements in this redshift range, and foreshadow the advancements which will be generated by the main HETDEX survey.

preprint2013arXiv

The Spectrally Resolved Lyman-alpha Emission of Three Lyman-alpha Selected Field Galaxies at z~2.4 from the HETDEX Pilot Survey

We present new results on the spectrally resolved Lyman-alpha (LyA) emission of three LyA emitting field galaxies at z~2.4 with high LyA equivalent width (>100 Angstroms) and LyA luminosity (~10^43 erg/s). At 120 km/s (FWHM) spectral resolution, the prominent double-peaked LyA profile straddles the systemic velocity, where the velocity zero-point is determined from spectroscopy of the galaxies' rest-frame optical nebular emission lines. The average velocity offset from systemic of the stronger redshifted emission component for our sample is 176 km/s while the average total separation between the redshifted and main blueshifted emission components is 380 km/s. These measurements are a factor of ~2 smaller than for UV continuum-selected galaxies that show LyA in emission with lower LyA equivalent width. We compare our LyA spectra to the predicted line profiles of a spherical "expanding shell" LyA radiative transfer grid that models large-scale galaxy outflows. Specifically blueward of the systemic velocity where two galaxies show a weak, highly blueshifted (by ~1000 km/s) tertiary emission peak, the model line profiles are a relatively poor representation of the observed spectra. Since the neutral gas column density has a dominant influence over the shape of the LyA line profile, we caution against equating the observed LyA velocity offset with a physical outflow velocity, especially at lower spectral resolution where the unresolved LyA velocity offset is a convoluted function of several degenerate parameters. Referring to rest-frame ultraviolet and optical Hubble Space Telescope imaging, we find that galaxy-galaxy interactions may play an important role in inducing a starburst that results in copious LyA emission, as well as perturbing the gas distribution and velocity field which have strong influence over the LyA emission line profile.

preprint2013arXiv

To Stack or Not to Stack: Spectral Energy Distribution Properties of Lya-Emitting Galaxies at z=2.1

We use the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS) GOODS-S multi-wavelength catalog to identify counterparts for 20 Lya Emitting (LAE) galaxies at z=2.1. We build several types of stacked Spectral Energy Distributions (SEDs) of these objects. We combine photometry to form average and median flux-stacked SEDs, and postage stamp images to form average and median image-stacked SEDs. We also introduce scaled flux stacks that eliminate the influence of variation in overall brightness. We use the SED fitting code SpeedyMC to constrain the physical properties of individual objects and stacks. Our LAEs at z = 2.1 have stellar masses ranging from 2x10^7 Msun - 8x10^9 Msun (median = 3x10^8 Msun), ages ranging from 4 Myr to 500 Myr (median =100 Myr), and E(B-V) between 0.02 and 0.24 (median = 0.12). We do not observe strong correlations between Lya equivalent width (EW) and stellar mass, age, or E(B-V). The Lya radiative transfer (q) factors of our sample are predominantly close to one and do not correlate strongly with EW or E(B-V), implying that Lya radiative transfer prevents Lya photons from resonantly scattering in dusty regions. The SED parameters of the flux stacks match the average and median values of the individual objects, with the flux-scaled median SED performing best with reduced uncertainties. Median image-stacked SEDs provide a poor representation of the median individual object, and none of the stacking methods captures the large dispersion of LAE properties.

preprint2013arXiv

UVUDF: Ultraviolet Imaging of the Hubble Ultradeep Field with Wide-field Camera 3

We present an overview of a 90-orbit Hubble Space Telescope treasury program to obtain near ultraviolet imaging of the Hubble Ultra Deep Field using the Wide Field Camera 3 UVIS detector with the F225W, F275W, and F336W filters. This survey is designed to: (i) Investigate the episode of peak star formation activity in galaxies at 1<z<2.5; (ii) Probe the evolution of massive galaxies by resolving sub-galactic units (clumps); (iii) Examine the escape fraction of ionizing radiation from galaxies at z~2-3; (iv) Greatly improve the reliability of photometric redshift estimates; and (v) Measure the star formation rate efficiency of neutral atomic-dominated hydrogen gas at z~1-3. In this overview paper, we describe the survey details and data reduction challenges, including both the necessity of specialized calibrations and the effects of charge transfer inefficiency. We provide a stark demonstration of the effects of charge transfer inefficiency on resultant data products, which when uncorrected, result in uncertain photometry, elongation of morphology in the readout direction, and loss of faint sources far from the readout. We agree with the STScI recommendation that future UVIS observations that require very sensitive measurements use the instrument's capability to add background light through a "post-flash". Preliminary results on number counts of UV-selected galaxies and morphology of galaxies at z~1 are presented. We find that the number density of UV dropouts at redshifts 1.7, 2.1, and 2.7 is largely consistent with the number predicted by published luminosity functions. We also confirm that the image mosaics have sufficient sensitivity and resolution to support the analysis of the evolution of star-forming clumps, reaching 28-29th magnitude depth at 5 sigma in a 0.2 arcsecond radius aperture depending on filter and observing epoch.

preprint2012arXiv

A CANDELS WFC3 Grism Study of Emission-Line Galaxies at z~2: A Mix of Nuclear Activity and Low-Metallicity Star Formation

We present Hubble Space Telescope Wide Field Camera 3 slitless grism spectroscopy of 28 emission-line galaxies at z~2, in the GOODS-S region of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). The high sensitivity of these grism observations, with 1-sigma detections of emission lines to f > 2.5x10^{-18} erg/s/cm^2, means that the galaxies in the sample are typically ~7 times less massive (median M_* = 10^{9.5} M_sun) than previously studied z~2 emission-line galaxies. Despite their lower mass, the galaxies have OIII/Hb ratios which are very similar to previously studied z~2 galaxies and much higher than the typical emission-line ratios of local galaxies. The WFC3 grism allows for unique studies of spatial gradients in emission lines, and we stack the two-dimensional spectra of the galaxies for this purpose. In the stacked data the OIII emission line is more spatially concentrated than the Hb emission line with 98.1 confidence. We additionally stack the X-ray data (all sources are individually undetected), and find that the average L(OIII)/L(0.5-10 keV) ratio is intermediate between typical z~0 obscured active galaxies and star-forming galaxies. Together the compactness of the stacked OIII spatial profile and the stacked X-ray data suggest that at least some of these low-mass, low-metallicity galaxies harbor weak active galactic nuclei.

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

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.

preprint2011arXiv

CANDELS: The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) is designed to document the first third of galactic evolution, over the approximate redshift (z) range 8--1.5. It will image >250,000 distant galaxies using three separate cameras on the Hubble Space Telescope, from the mid-ultraviolet to the near-infrared, and will find and measure Type Ia supernovae at z>1.5 to test their accuracy as standardizable candles for cosmology. Five premier multi-wavelength sky regions are selected, each with extensive ancillary data. The use of five widely separated fields mitigates cosmic variance and yields statistically robust and complete samples of galaxies down to a stellar mass of 10^9 M_\odot to z \approx 2, reaching the knee of the ultraviolet luminosity function (UVLF) of galaxies to z \approx 8. The survey covers approximately 800 arcmin^2 and is divided into two parts. The CANDELS/Deep survey (5σ point-source limit H=27.7 mag) covers \sim 125 arcmin^2 within GOODS-N and GOODS-S. The CANDELS/Wide survey includes GOODS and three additional fields (EGS, COSMOS, and UDS) and covers the full area to a 5σ point-source limit of H \gtrsim 27.0 mag. Together with the Hubble Ultra Deep Fields, the strategy creates a three-tiered "wedding cake" approach that has proven efficient for extragalactic surveys. Data from the survey are nonproprietary and are useful for a wide variety of science investigations. In this paper, we describe the basic motivations for the survey, the CANDELS team science goals and the resulting observational requirements, the field selection and geometry, and the observing design. The Hubble data processing and products are described in a companion paper.

preprint2011arXiv

CANDELS: The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey - The Hubble Space Telescope Observations, Imaging Data Products and Mosaics

This paper describes the Hubble Space Telescope imaging data products and data reduction procedures for the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). This survey is designed to document the evolution of galaxies and black holes at $z\sim1.5-8$, and to study Type Ia SNe beyond $z>1.5$. Five premier multi-wavelength sky regions are selected, each with extensive multiwavelength observations. The primary CANDELS data consist of imaging obtained in the Wide Field Camera 3 / infrared channel (WFC3/IR) and UVIS channel, along with the Advanced Camera for Surveys (ACS). The CANDELS/Deep survey covers \sim125 square arcminutes within GOODS-N and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of \sim800 square arcminutes across GOODS and three additional fields (EGS, COSMOS, and UDS). We summarize the observational aspects of the survey as motivated by the scientific goals and present a detailed description of the data reduction procedures and products from the survey. Our data reduction methods utilize the most up to date calibration files and image combination procedures. We have paid special attention to correcting a range of instrumental effects, including CTE degradation for ACS, removal of electronic bias-striping present in ACS data after SM4, and persistence effects and other artifacts in WFC3/IR. For each field, we release mosaics for individual epochs and eventual mosaics containing data from all epochs combined, to facilitate photometric variability studies and the deepest possible photometry. A more detailed overview of the science goals and observational design of the survey are presented in a companion paper.

preprint2011arXiv

First spectroscopic measurements of [OIII] emission from Lyman-alpha selected field galaxies at z ~ 3.1

We present the first spectroscopic measurements of the [OIII] 5007 A line in two z ~ 3.1 Lyman-alpha emitting galaxies (LAEs) using the new near-infrared instrument LUCIFER1 on the 8.4m Large Binocular Telescope (LBT). We also describe the optical imaging and spectroscopic observations used to identify these Lya emitting galaxies. Using the [OIII] line we have measured accurate systemic redshifts for these two galaxies, and discovered a velocity offset between the [OIII] and Ly-alpha lines in both, with the Lya line peaking 342 and 125 km/s redward of the systemic velocity. These velocity offsets imply that there are powerful outflows in high-redshift LAEs. They also ease the transmission of Lya photons through the interstellar medium and intergalactic medium around the galaxies. By measuring these offsets directly, we can refine both Lya-based tests for reionization, and Lya luminosity function measurements where the Lya forest affects the blue wing of the line. Our work also provides the first direct constraints on the strength of the [OIII] line in high-redshift LAEs. We find [OIII] fluxes of 7 and 36 x 10^-17 erg s^-1 cm^-2 in two z ~ 3.1 LAEs. These lines are strong enough to dominate broad-band flux measurements that include the line (in thiscase, K_s band photometry). Spectral energy distribution fits that do not account for the lines would therefore overestimate the 4000 A (and/or Balmer) break strength in such galaxies, and hence also the ages and stellar masses of such high-z galaxies.

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

Sizing up Lyman-alpha and Lyman Break Galaxies

We show that populations of high redshift galaxies grow more luminous as they grow in linear size. This is because the brightness per unit area has a distinct upper limit due to the self-regulating nature of star-formation. As a corollary, we show that the observed increase in characteristic luminosity of Lyman Break Galaxies (L*) with time can be explained by their increase in size, which scales as 1/H(z). In contrast, Lyman-alpha selected galaxies have a characteristic, constant, small size between redshift z=2.25 and 6.5. Coupled with a characteristic surface brightness, this can explain their non-evolving ultraviolet continuum luminosity function. This compact physical size seems to be a critical determining factor in whether a galaxy will show Lyman-alpha emission. We base these conclusions on new size measurements for a sample of about 170 Lyman-alpha selected galaxies with Hubble Space Telescope broad-band imaging, over the redshift range 2.25 < z < 6. We combine these with a similar number of Lyman-alpha selected galaxies with half-light radii from the literature, and derive surface brightnesses for the entire combined sample.

preprint2011arXiv

The HETDEX Pilot Survey. II. The Evolution of the Ly-alpha Escape Fraction from the UV Slope and Luminosity Function of 1.9<z<3.8 LAEs

We study the escape of Ly-alpha photons from Ly-alpha emitting galaxies (LAEs) and the overall galaxy population using a sample of 99 LAEs at 1.9<z<3.8 detected through integral-field spectroscopy of blank fields by the HETDEX Pilot Survey. For 89 LAEs showing counterparts in deep broad-band images we measure the rest-frame UV luminosity and the UV slope, which we use to estimate E(B-V) under the assumption of a constant intrinsic UV slope for LAEs. These two quantities are used to measure the dust-corrected star formation rate (SFR). A comparison between the observed Ly-alpha luminosity and that predicted by the dust-corrected SFR yields the Ly-alpha escape fraction. We also measure the Ly-alpha luminosity function. Integration of the luminosity function provides a measurement of the Ly-alpha luminosity density across our redshift range. We combine our data with that from other surveys at 0.3<z<7.7 to trace the evolution of the Ly-alpha luminosity density. We then compare it to that expected from the star-formation history of the universe in order to characterize the evolution of the Ly-alpha escape fraction of the overall galaxy population [abriged]

preprint2011arXiv

X-ray Constraints on the Lyman-Alpha Escape Fraction

We have coadded X-ray flux of all known Lyman Alpha Emitters (LAEs) in the 4 Msec Chandra Deep Field South (CDF-S) region, to place sensitive upper limits on the average unobscured star-formation rate (SFR_X) in these galaxies. A very small fraction of Lyman-Alpha galaxies in the field are individually detected in the X-rays, implying a low fraction of AGN activity. After excluding the few X-ray detected LAEs, we stack the undetected LAEs located in the 4 Ms CDF-S data and 250 ks Extended CDF-S (ECDFS) data, and compute a 1-σupper limit on SFR_X < 14, 28, 28, 140, 440, 880 M$_{\sun}$ yr$^{-1}$ for LAEs located at z = 2.1, 3.1, 3.2, 4.5, 5.7 and 6.5, respectively. The upper limit of SFR_X in LAEs can be then be compared to SFR$_{Lyα}$ derived from Lyman-Alpha line and thus can constrain on the Lyman-Alpha escape fraction ($f^{Esc}_{Lyα}$). We derive a lower limit on f(Lyman-Alpha Escape) > 14% (84 % confidence level, 1-σlower limit) for LAEs at redshift z ~ 2.1 and z ~ 3.1-3.2. At z > 4, the current LAE samples are not of sufficient size to constrain SFR_X well. By averaging all the LAEs at z> 2, the X-ray non-detection constrains f(Lyman-Alpha Escape) > 17% (84 % confidence level, 1-σlower limit), and rejects f(Lyman-Alpha Escape) < 5.7% at the 99.87% confidence level from 2.1 < z < 6.5.

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

HETDEX pilot survey for emission-line galaxies - I. Survey design, performance, and catalog

We present a catalog of emission-line galaxies selected solely by their emission-line fluxes using a wide-field integral field spectrograph. This work is partially motivated as a pilot survey for the upcoming Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). We describe the observations, reductions, detections, redshift classifications, line fluxes, and counterpart information for 397 emission-line galaxies detected over 169 sq.arcmin with a 3500-5800 Ang. bandpass under 5 Ang. full-width-half-maximum (FWHM) spectral resolution. The survey's best sensitivity for unresolved objects under photometric conditions is between 4-20 E-17 erg/s/sq.cm depending on the wavelength, and Ly-alpha luminosities between 3-6 E42 erg/s are detectable. This survey method complements narrowband and color-selection techniques in the search for high redshift galaxies with its different selection properties and large volume probed. The four survey fields within the COSMOS, GOODS-N, MUNICS, and XMM-LSS areas are rich with existing, complementary data. We find 104 galaxies via their high redshift Ly-alpha emission at 1.9<z<3.8, and the majority of the remainder objects are low redshift [OII]3727 emitters at z<0.56. The classification between low and high redshift objects depends on rest frame equivalent width, as well as other indicators, where available. Based on matches to X-ray catalogs, the active galactic nuclei (AGN) fraction amongst the Ly-alpha emitters (LAEs) is 6%. We also analyze the survey's completeness and contamination properties through simulations. We find five high-z, highly-significant, resolved objects with full-width-half-maximum sizes >44 sq.arcsec which appear to be extended Ly-alpha nebulae. We also find three high-z objects with rest frame Ly-alpha equivalent widths above the level believed to be achievable with normal star formation, EW(rest)>240 Ang.

preprint2010arXiv

Hubble Space Telescope Imaging of Lyman Alpha Emission at z=4.4

We present the highest redshift detections of resolved Lyman alpha emission, using Hubble Space Telescope/ACS F658N narrowband-imaging data taken in parallel with the Wide Field Camera 3 Early Release Science program in the GOODS CDF-S. We detect Lyman alpha emission from three spectroscopically confirmed z = 4.4 Lyman alpha emitting galaxies (LAEs), more than doubling the sample of LAEs with resolved Lyman alpha emission. Comparing the light distribution between the rest-frame ultraviolet continuum and narrowband images, we investigate the escape of Lyman alpha photons at high redshift. While our data do not support a positional offset between the Lyman alpha and rest-frame ultraviolet (UV) continuum emission, the half-light radii in two out of the three galaxies are significantly larger in Lyman alpha than in the rest-frame UV continuum. This result is confirmed when comparing object sizes in a stack of all objects in both bands. Additionally, the narrowband flux detected with HST is significantly less than observed in similar filters from the ground. These results together imply that the Lyman alpha emission is not strictly confined to its indigenous star-forming regions. Rather, the Lyman alpha emission is more extended, with the missing HST flux likely existing in a diffuse outer halo. This suggests that the radiative transfer of Lyman alpha photons in high-redshift LAEs is complicated, with the interstellar-medium geometry and/or outflows playing a significant role in galaxies at these redshifts.

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

Probing Stellar Populations at z ~ 7 - 8

In this proceeding we present the results from a study of very high-redshift galaxies with the newly commissioned Wide Field Camera 3 on the Hubble Space Telescope. With the deepest near-infrared data ever taken, we discovered 31 galaxies at 6.3 < z < 8.6. The rest-frame ultraviolet (UV) colors of these galaxies are extremely blue, showing significant (> 4 sigma) evolution from z ~ 3, over only 1 Gyr of cosmic time. While we cannot yet diagnose the exact cause of the bluer colors, it appears a low dust content is the primary factor. The stellar masses of these galaxies are less than comparably selected galaxies at 3 < z < 6, highlighting evolution in the stellar mass of characteristic (L*) galaxies with redshift. Lastly, the measured rest-UV luminosity density of galaxies in our sample seems sufficient to sustain reionization at z ~ 7 when we account for the likely contribution from galaxies below our magnitude limit.

preprint2010arXiv

Reversal of Fortune: Confirmation of an Increasing Star Formation-Density Relation in a Cluster at z=1.62

We measure the rest-frame colors (dust-corrected), infrared luminosities, star formation rates, and stellar masses of 92 galaxies in a Spitzer-selected cluster at z=1.62. By fitting spectral energy distributions (SEDs) to 10-band photometry (0.4 micron < lambda(obs) <8 micron) and measuring 24 micron fluxes for the 12 spectroscopically confirmed and 80 photometrically selected members, we discover an exceptionally high level of star formation in the cluster core of ~1700 Msun/yr per Mpc^2. The cluster galaxies define a strong blue sequence in (U-V) color and span a range in color. We identify 17 members with L(IR)>10^(11) Lsun, and these IR luminous members follow the same trend of increasing star formation with stellar mass that is observed in the field at z~2. Using rates derived from both the 24 micron imaging and SED fitting, we find that the relative fraction of star-forming members triples from the lowest to highest galaxy density regions, e.g. the IR luminous fraction increases from ~8% at Sigma~10 gal per Mpc^2 to ~25% at Sigma>100 gal per Mpc^2. The observed increase is a reversal of the well-documented trend at z<1 and signals that we have reached the epoch when massive cluster galaxies are still forming a substantial fraction of their stars.

preprint2010arXiv

Searching for the First Galaxies

As some of the first known objects to exist in the Universe, Lyman alpha emitting galaxies (LAEs) naturally draw a lot of interest. First discovered over a decade ago, they have allowed us to probe the early Universe, as their strong emission line compensates for their faint continuum light. While initially thought to be indicative of the first galaxies forming in the Universe, recent studies have shown them to be increasingly complex, as some fraction appear evolved, and many LAEs appear to be dusty, which one would not expect from primordial galaxies. Presently, much interest resides in discovering not only the highest redshift galaxies to constrain theories of reionization, but also pushing closer to home, as previous ground-based studies have only found LAEs at z > 3 due to observational limitations. In this review talk I will cover everything from the first theoretical predictions of LAEs, to their future prospects for study, including the HETDEX survey here in Texas.

preprint2010arXiv

The Luminosity Function of Lyman alpha Emitters at Redshift z=7.7

Lyman alpha (Lya) emission lines should be attenuated in a neutral intergalactic medium (IGM). Therefore the visibility of Lya emitters at high redshifts can serve as a valuable probe of reionization at about the 50% level. We present an imaging search for z=7.7 Lya emitting galaxies using an ultra-narrowband filter (filter width= 9A) on the NEWFIRM imager at the Kitt Peak National Observatory. We found four candidate Lya emitters in a survey volume of 1.4 x 10^4 Mpc^3, with a line flux brighter than 6x10^-18 erg/cm^2/s (5 sigma in 2" aperture). We also performed a detailed Monte-Carlo simulation incorporating the instrumental effects to estimate the expected number of Lya emitters in our survey, and found that we should expect to detect one Lya emitter, assuming a non-evolving Lya luminosity function (LF) between z=6.5 and z=7.7. Even if one of the present candidates is spectroscopically confirmed as a z~8 Lya emitter, it would indicate that there is no significant evolution of the Lya LF from z=3.1 to z~8. While firm conclusions would need both spectroscopic confirmations and larger surveys to boost the number counts of galaxies, we successfully demonstrate the feasibility of sensitive near-infrared (1.06 um) narrow-band searches using custom filters designed to avoid the OH emission lines that make up most of the sky background.

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

The Shining Future of UV Spectral Synthesis

With the coming generation of instruments and telescopes capable of spectroscopy of high redshift galaxies, the spectral synthesis technique in the rest-frame UV and Far-UV range will become one of a few number of tools remaining to study their young stellar populations in detail. The rest-frame UV lines and continuum of high redshift galaxies, observed with visible and infrared telescopes on Earth, can be used for accurate line profile fitting such as PV@1118,1128A, CIII@1176A, and CIV@1550A. These lines are very precise diagnostic tools to estimate ages, metallicities, and masses of stellar populations. Here we discuss the potential for spectral synthesis of rest-frame UV spectra obtained at the Keck telescope. As an example, we study the 8 o'clock arc, a lensed galaxy at z=2.7322. We show that the poor spectral type coverage of the actual UV empirical spectral libraries limits the age and metallicity diagnostic. In order to improve our knowledge of high redshift galaxies using spectral synthesis, UV stellar libraries need to be extended to obtain accurate age, metallicity, and mass estimates likely to be occuring in young stellar populations observed in the early universe.

preprint2008arXiv

Lyman Alpha Galaxies: Primitive, Dusty or Evolved Galaxies?

We present stellar population modeling results for 10 newly discovered Lyman alpha emitting galaxies (LAEs), as well as four previously known LAEs at z ~ 4.5 in the Chandra Deep Field - South. We fit stellar population models to these objects in order to learn specifically if there exists more than one class of LAE. Past observational and theoretical evidence has shown that while many LAEs appear to be young, they may be much older, with Lyman alpha EWs enhanced due to resonant scattering of Lyman alpha photons in a clumpy interstellar medium (ISM). Our results show a large range of stellar population age (3 - 500 Myr), stellar mass (1.6 x 10^8 - 5.0 x 10^10 Msol) and dust extinction (A_1200 = 0.3 - 4.5 mag), broadly consistent with previous studies. With such a large number of individually analyzed objects, we have looked at the distribution of stellar population ages in LAEs for the first time, and we find a very interesting bimodality, in that our objects are either very young (< 15 Myr) or old (> 450 Myr). This bimodality may be caused by dust, and it could explain the Lyman alpha duty cycle which has been proposed in the literature. We find that eight of the young objects are best fit with a clumpy ISM. We find that dust geometry appears to play a large role in shaping the SEDs that we observe, and that it may be a major factor in the observed Lyman alpha equivalent width distribution in high redshift Lyman alpha galaxies, although other factors (i.e. outflows) may be in play. We conclude that 12 out of our 14 LAEs are dusty star-forming galaxies, with the other two LAEs being evolved galaxies.

preprint2007arXiv

The Ages and Masses of Lyman Alpha Galaxies at Redshift $z\sim 4.5$

We examine the stellar populations of a sample of 98 redshift 4.5 Lyman alpha emitting galaxies using their broadband colors derived from deep photometry at the MMT. These galaxies were selected by narrowband excess from the Large Area Lyman Alpha survey. Twenty-two galaxies are detected in two or more of our MMT filters (g', r', i' and z'). By comparing broad and narrowband colors of these galaxies to synthetic colors from stellar population models, we determine their ages and stellar masses. The highest equivalent width objects have an average age of 4 Myr, consistent with ongoing star formation. The lowest EW objects show an age of 40 - 200 Myr, consistent with the expectation that larger numbers of stars are causing low EWs. We found masses ranging from 2e7 solar masses for the youngest objects in the sample to 2e9 solar masses for the oldest. It is possible that dust effects could produce large equivalent widths even in older populations by allowing the Lyman alpha photons to escape, even while the continuum is extinguished, and we present models for this scenario also.