Researcher profile

Steven L. Finkelstein

Steven L. Finkelstein contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
41works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

41 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&#39; $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&#39; 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&#34;$), 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 &#34;OIR-dark&#34; 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 &#34;unVO87&#34; diagram for 461 galaxies and $\frac{[OIII]}{Hb}$ vs. $\frac{[NeIII]}{[OII]}$ as an &#34;OHNO&#34; 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&#39; 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&#39;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&#39;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&#39;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 &#34;non-parametric&#34; 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 &#34;dry&#34; 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&#39; galaxy formation physics to `top-level&#39; 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{&#39;}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&#39;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&#39;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&#39;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 &#34;over-consumption&#34; 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&#39;&#39;), 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Å).