Researcher profile

Daniel P. Stark

Daniel P. Stark contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
23works
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

23 published item(s)

preprint2026arXiv

REBELS-MOSFIRE: Weak CIII] Emission is Typical Among Extremely UV-bright, Massive Galaxies at $z\sim7$

We present Keck/MOSFIRE H-band spectroscopic measurements covering the [CIII]1907, CIII]1909 doublet for a sample of 8 z~7 spectroscopically-confirmed star-forming galaxies drawn from the Reionization Era Bright Emission Line Survey (REBELS). This REBELS-MOSFIRE sample is notable for its bright median UV luminosity (Muv=-22.5 AB) and large median stellar mass (log(Mstar/Msun)=9.2). Although three sources show tentative evidence of a CIII] detection, we obtain no confident detections for any of the 8 REBELS-MOSFIRE sources. The median [CIII]1907+CIII]1909 3-sigma upper limit in equivalent width (EW) for the REBELS-MOSFIRE sample is 6.5 AA, and a stack of their H-band MOSFIRE spectra yields a non-detection with an associated 3-sigma upper limit of 2.6 AA. These upper limits fall significantly below the CIII] EW measured in a composite spectrum of representative z~7 star-forming galaxies, as well as those measured for notable early star-forming galaxies such as GN-z11, GHZ2, GS-z12, and RXCJ2248-ID. The lack of strong CIII] emission can be understood within the context of the stellar populations of the REBELS galaxies, as well as the ionization conditions and gas-phase metallicity implied by rest-frame optical spectroscopic properties ([OIII]+Hb EWs, and [OIII]5007/[OII]3727 and [NeIII]3869/[OII]3727 line ratios). The REBELS-MOSFIRE sample represents the higher-mass, higher-metallicity, lower-excitation tail of the z~7 galaxy population, whose ionizing properties must be fully characterized to constrain the role of star-forming galaxies during cosmic reionization.

preprint2025arXiv

Constraints on the $z\sim6-13$ intergalactic medium from JWST spectroscopy of Lyman-alpha damping wings in galaxies

JWST provides a unique dataset to study reionization&#39;s earliest stages, promising insights into the first galaxies. Many JWST/NIRSpec prism spectra reveal smooth Lyman-alpha breaks in z>5 galaxies, implying damping wing scattering by neutral hydrogen. We investigate what current prism spectra imply about the intergalactic medium (IGM), and how best to use NIRSpec spectra to recover IGM properties. We use a sample of 99 z~5.5-13 galaxies with high S/N prism spectra in the public archive, including 12 at z>10. We analyse these spectra using damping wing sightlines from inhomogeneous reionizing IGM simulations, mapping between the distance of a source from neutral IGM and the mean IGM neutral fraction. We marginalise over absorption by local HI around galaxies, and Lyman-alpha emission. We observe a decline in the median and variance of flux around the Lyman-alpha break with increasing redshift, consistent with an increasingly neutral IGM, as ionized regions become smaller and rarer. At $z\gtrsim9$ the spectra become consistent with an almost fully neutral IGM. We find S/N>15 per pixel is required to robustly estimate IGM properties from prism spectra. We fit a sub-sample of high S/N spectra and infer mean neutral fractions $\overline{x}_\mathrm{HI}=0.33^{+0.18}_{-0.27}, 0.64^{+0.17}_{-0.23}$ ($>0.70$ excluding GNz11) at $z \approx 6.5, 9.3$. We also investigate local HI absorption, finding median column density $\log_{10}N_\mathrm{HI}\approx10^{20.8}$ cm$^{-2}$, comparable to $z\sim3$ Lyman-break galaxies, with no significant redshift evolution. We find galaxies showing the strongest absorption are more likely to be in close associations (<500 pkpc), implying enhanced absorption in massive dark matter halos. Future deep prism and grating spectroscopy of z>9 sources will provide tighter constraints on the earliest stages of reionization, key for understanding the onset of star formation.

preprint2024arXiv

Inside the bubble: exploring the environments of reionisation-era Lyman-$α$ emitting galaxies with JADES and FRESCO

We present a study of the environments of 17 Lyman-$α$ (Ly$α$) emitting galaxies (LAEs) in the reionisation era ($5.8 < z < 8$) identified by JWST/NIRSpec as part of the JWST Advanced Deep Extragalactic Survey (JADES). Unless situated in sufficiently (re)ionised regions, Ly$α$ emission from these galaxies would be strongly absorbed by neutral gas in the intergalactic medium (IGM). We conservatively estimate sizes of the ionised regions required to reconcile the relatively low Ly$α$ velocity offsets ($Δv_\text{Ly$α$}<300\,\mathrm{km\,s^{-1}}$) with moderately high Ly$α$ escape fractions ($f_\mathrm{esc,\,Lyα}>5\%$) observed in our sample of LAEs, suggesting the presence of ionised hydrogen along the line of sight towards at least eight out of 17 LAEs. We find minimum physical `bubble&#39; sizes of the order of $R_\text{ion}\sim0.1$-$1\,\mathrm{pMpc}$ are required in a patchy reionisation scenario where ionised bubbles containing the LAEs are embedded in a fully neutral IGM. Around half of the LAEs in our sample are found to coincide with large-scale galaxy overdensities seen in FRESCO at $z \sim 5.8$-$5.9$ and $z\sim7.3$, suggesting Ly$α$ transmission is strongly enhanced in such overdense regions, and underlining the importance of LAEs as tracers of the first large-scale ionised bubbles. Considering only spectroscopically confirmed galaxies, we find our sample of UV-faint LAEs ($M_\text{UV}\gtrsim-20\,\mathrm{mag}$) and their direct neighbours are generally not able to produce the required ionised regions based on the Ly$α$ transmission properties, suggesting lower-luminosity sources likely play an important role in carving out these bubbles. These observations demonstrate the combined power of JWST multi-object and slitless spectroscopy in acquiring a unique view of the early Universe during cosmic reionisation via the most distant LAEs.

preprint2023arXiv

Star formation histories of UV-luminous galaxies at $z \simeq 6.8$: implications for stellar mass assembly at early cosmic times

The variety of star formation histories (SFHs) of $z\gtrsim6$ galaxies provides important insights into early star formation, but has been difficult to systematically quantify. Some observations suggest that many $z\sim6-9$ galaxies are dominated by $\gtrsim200$ Myr stellar populations, implying significant star formation at $z\gtrsim9$, while others find that most reionization era galaxies are $\lesssim10$ Myr, consistent with little $z\gtrsim9$ star formation. Here, we quantify the distribution of ages of UV-bright ($-22.5\lesssim M_{UV}\lesssim-21$) galaxies colour-selected to lie at $z\simeq6.6-6.9$, an ideal redshift range to systematically study the SFHs of reionization era galaxies with ground-based observatories and Spitzer. We infer galaxy properties with two SED modelling codes and compare results, finding that stellar masses are largely insensitive to the model, but the inferred ages can vary by an order of magnitude. We infer a distribution of ages assuming a simple, parametric SFH model, finding a median age of $\sim30-70$ Myr depending on SED model. We quantify the fractions of $\leq10$ Myr and $\geq250$ Myr galaxies, finding that these systems comprise $\sim15-30$ percent and $\sim20-25$ percent of the population, respectively. With a flexible SFH model, the shapes of the SFHs are consistent with those implied by the simple model (e.g. young galaxies have rapidly rising SFHs). However, stellar masses can differ significantly, with those of young systems sometimes being more than an order of magnitude larger with the flexible SFH. We quantify the implications of these results for $z\gtrsim9$ stellar mass assembly and discuss improvements expected from JWST.

preprint2022arXiv

A double-peaked Lyman-$α$ emitter with a strong blue peak multiply imaged by the galaxy cluster RXC J0018.5+1626

We report the discovery of a double-peaked Lyman-$α$ (Ly$α$) emitter (LAE) at $z=3.2177\pm0.0001$ in VLT/MUSE data. The galaxy is strongly lensed by the galaxy cluster RXC~J0018.5+1626 recently observed in the RELICS survey, and the double-peaked Ly$α$ emission is clearly detected in the two counter images in the MUSE field-of-view. We measure a relatively high Ly$α$ rest-frame equivalent width (EW) of $\mathrm{EW}_{\mathrm{Ly}α,0}=(63\pm2)\,\mathring{\mathrm{A}}$. Additional near-infrared (NIR) spectroscopy allows us to measure the H$β$, [OIII]$\lambda4959\,\mathring{\mathrm{A}}$ and [OIII]$\lambda5007\,\mathring{\mathrm{A}}$ emission lines, which show moderate rest-frame EWs of the order of a few $\sim10-100\,\mathring{\mathrm{A}}$, an [OIII]$\lambda5007\,\mathring{\mathrm{A}}$/H$β$ ratio of $4.8\pm0.7$, and a lower limit on the [OIII]/[OII] ratio of $>5.6$. The galaxy has very blue UV-continuum slopes of $β_{\mathrm{FUV}}=-2.23\pm0.06$ and $β_{\mathrm{NUV}}=-3.0\pm0.2$, and is magnified by factors $μ\sim7-10$ in each of the two images, thus enabling a view into a low-mass ($M_{\star}\simeq10^{7.5}\,\mathrm{M}_{\odot}$) high-redshift galaxy analog. Notably, the blue peak of the Ly$α$ profile is significantly stronger than the red peak, which suggests an inflow of matter and possibly very low HI column densities in its circumgalactic gas. Combined with the high lensing magnification and image multiplicity, these properties make this galaxy a prime candidate for follow-up observations to search for LyC emission and constrain the LyC photon escape fraction.

preprint2022arXiv

A Panchromatic Study of Massive Stars in the Extremely Metal-Poor Local Group Dwarf Galaxy Leo A

We characterize massive stars (M>8 M_sun) in the nearby (D~0.8 Mpc) extremely metal-poor (Z~5% Z_sun) galaxy Leo A using Hubble Space Telescope ultra-violet (UV), optical, and near-infrared (NIR) imaging along with Keck/LRIS and MMT/Binospec optical spectroscopy for 18 main sequence OB stars. We find that: (a) 12 of our 18 stars show emission lines, despite not being associated with an H II region, suggestive of stellar activity (e.g., mass loss, accretion, binary star interaction), which is consistent with previous predictions of enhanced activity at low metallicity; (b) 6 are Be stars, which are the first to be spectroscopically studied at such low metallicity -- these Be stars have unusual panchromatic SEDs; (c) for stars well-fit by the TLUSTY non-local thermodynamic equilibrium (non-LTE) models, the photometric and spectroscopic values of T_eff and log(g) agree to within ~0.01 dex and ~0.18 dex, respectively, indicating that NUV/optical/NIR imaging can be used to reliably characterize massive (M ~ 8-30 M_sun) main sequence star properties relative to optical spectroscopy; (d) the properties of the most massive stars in H II regions are consistent with constraints from previous nebular emission line studies; and (e) 13 stars with M>8 M_sun are >40 pc from a known star cluster or H II region. Our sample comprises ~50% of all known massive stars at Z < 10% Z_sun with derived stellar parameters, high-quality optical spectra, and panchromatic photometry.

preprint2022arXiv

CLASSY III: The Properties of Starburst-Driven Warm Ionized Outflows

We report the results of analyses of galactic outflows in a sample of 45 low-redshift starburst galaxies in the COS Legacy Archive Spectroscopic SurveY (CLASSY), augmented by five additional similar starbursts with COS data. The outflows are traced by blueshifted absorption-lines of metals spanning a wide range of ionization potential. The high quality and broad spectral coverage of CLASSY data enable us to disentangle the absorption due to the static ISM from that due to outflows. We further use different line multiplets and doublets to determine the covering fraction, column density, and ionization state as a function of velocity for each outflow. We measure the outflow&#39;s mean velocity and velocity width, and find that both correlate in a highly significant way with the star-formation rate, galaxy mass, and circular velocity over ranges of four orders-of-magnitude for the first two properties. We also estimate outflow rates of metals, mass, momentum, and kinetic energy. We find that, at most, only about 20% of silicon created and ejected by supernovae in the starburst is carried in the warm phase we observe. The outflows&#39; mass-loading factor increases steeply and inversely with both circular and outflow velocity (log-log slope $\sim$ -1.6), and reaches $\sim 10$ for dwarf galaxies. We find that the outflows typically carry about 10 to 100% of the momentum injected by massive stars and about 1 to 20% of the kinetic energy. We show that these results place interesting constraints on, and new insights into, models and simulations of galactic winds.

preprint2022arXiv

Radio and Far-IR Emission Associated with a Massive Star-forming Galaxy Candidate at z$\simeq$6.8: A Radio-Loud AGN in the Reionization Era?

We report the identification of radio (0.144-3 GHz), mid-IR, far-IR, and sub-mm (24-850$μ$m) emission at the position of one of 41 UV-bright (M$_\mathrm{UV}^{}\lesssim-21.25$) $z\simeq6.6-6.9$ Lyman-break galaxy candidates in the 1.5 deg$^2$ COSMOS field. This source, COS-87259, exhibits a sharp flux discontinuity (factor $>$3) between two narrow/intermediate bands at 9450 and 9700 Angstroms and is undetected in all nine bands blueward of 9600 Angstroms, as expected from a Lyman-alpha break at $z\simeq6.8$. The full multi-wavelength (X-ray through radio) data of COS-87529 can be self-consistently explained by a very massive (M$_{\ast}=10^{10.8}$ M$_{\odot}$) and extremely red (rest-UV slope $β=-0.59$) $z\simeq6.8$ galaxy with hyperluminous infrared emission (L$_{\mathrm{IR}}=10^{13.6}$ L$_{\odot}$) powered by both an intense burst of highly-obscured star formation (SFR$\approx$1800 M$_{\odot}$ yr$^{-1}$) and an obscured ($τ_{\mathrm{9.7μm}}=7.7\pm2.5$) radio-loud (L$_{\mathrm{1.4\ GHz}}\approx10^{25.4}$ W Hz$^{-1}$) AGN. The radio emission is compact (1.04$\pm$0.12 arcsec) and exhibits an ultra-steep spectrum between 1.32-3 GHz ($α=-1.57^{+0.22}_{-0.21}$) that flattens at lower frequencies ($α=-0.86^{+0.22}_{-0.16}$ between 0.144-1.32 GHz), consistent with known $z>4$ radio galaxies. We also demonstrate that COS-87259 may reside in a significant (11$\times$) galaxy overdensity at $z\simeq6.6-6.9$, as common for systems hosting radio-loud AGN. Nonetheless, a spectroscopic redshift will ultimately be required to establish the true nature of COS-87259 as we cannot yet completely rule out low-redshift solutions. If confirmed to lie at $z\simeq6.8$, the properties of COS-87259 would be consistent with a picture wherein AGN and highly-obscured star formation activity are fairly common among very massive (M$_{\ast}>10^{10}$ M$_{\odot}$) reionization-era galaxies.

preprint2022arXiv

Stellar populations and star formation histories of the most extreme [OIII] emitters at $z=1.3-3.7$

As the James Webb Space Telescope approaches scientific operation, there is much interest in exploring the redshift range beyond that accessible with Hubble Space Telescope imaging. Currently, the only means to gauge the presence of such early galaxies is to age-date the stellar population of systems in the reionisation era. As a significant fraction of $z\simeq7-8$ galaxies are inferred from Spitzer photometry to have extremely intense [OIII] emission lines, it is commonly believed these are genuinely young systems that formed at redshifts $z<10$, consistent with a claimed rapid rise in the star formation density at that time. Here we study a spectroscopically-confirmed sample of extreme [OIII] emitters at $z=1.3-3.7$, using both dynamical masses estimated from [OIII] line widths and rest-frame UV to near-infrared photometry to illustrate the dangers of assuming such systems are genuinely young. For the most extreme of our intermediate redshift line emitters, we find dynamical masses $10-100$ times that associated with a young stellar population mass, which are difficult to explain solely by the presence of additional dark matter or gaseous reservoirs. Adopting nonparametric star formation histories, we show how the near-infrared photometry of a subset of our sample reveals an underlying old ($>100$ Myr) population whose stellar mass is $\simeq40$ times that associated with the starburst responsible for the extreme line emission. Without adequate rest-frame near-infrared photometry we argue it may be premature to conclude that extreme line emitters in the reionisation era are low mass systems that formed at redshifts below $z\simeq10$.

preprint2022arXiv

Strong Ly$α$ Emission in an Overdense Region at $z=6.8$: A Very Large ($R\sim3$ physical Mpc) Ionized Bubble in COSMOS?

Our understanding of reionization has advanced considerably over the past decade, with several results now demonstrating that the IGM transitioned from substantially neutral at $z=7$ to largely reionized at $z=6$. However, little remains known about the sizes of ionized bubbles at $z\gtrsim7$ as well as the galaxy overdensities which drive their growth. Fortunately, rest-UV spectroscopic observations offer a pathway towards characterizing these ionized bubbles thanks to the resonant nature of Lyman-alpha photons. In a previous work, we presented Ly$α$ detections from three closely-separated Lyman-break galaxies at $z\simeq6.8$, suggesting the presence of a large ($R>1$ physical Mpc) ionized bubble in the 1.5 deg$^2$ COSMOS field. Here, we present new deep Ly$α$ spectra of ten UV-bright ($\mathrm{M}_{\mathrm{UV}}^{} \leq -20.4$) $z\simeq6.6-6.9$ galaxies in the surrounding area, enabling us to better characterize this potential ionized bubble. We confidently detect (S/N$>$7) Ly$α$ emission at $z=6.701-6.882$ in nine of ten observed galaxies, revealing that the large-scale volume spanned by these sources (characteristic radius $R = 3.2$ physical Mpc) traces a strong galaxy overdensity ($N/\langle N\rangle \gtrsim 3$). Our data additionally confirm that the Ly$α$ emission of UV-bright galaxies in this volume is significantly enhanced, with 40% (4/10) showing strong Ly$α$ emission (equivalent width$>$25 $\mathrm{\mathring{A}}$) compared to the 8$-$9% found on average at $z\sim7$. The median Ly$α$ equivalent width of our observed galaxies is also $\approx$2$\times$ that typical at $z\sim7$, consistent with expectations if a very large ($R\sim3$ physical Mpc) ionized bubble is allowing the Ly$α$ photons to cosmologically redshift far into the damping wing before encountering HI.

preprint2022arXiv

The ALMA REBELS Survey: Specific Star-Formation Rates in the Reionization Era

We present specific star-formation rates for 40 UV-bright galaxies at $z\sim7-8$ observed as part of the Reionization Era Bright Emission Line Survey (REBELS) ALMA large program. The sSFRs are derived using improved measures of SFR and stellar masses, made possible by measurements of far-infrared (FIR) continuum emission and [CII]-based spectroscopic redshifts. For each source in the sample, we derive stellar masses from SED fitting and total SFRs from calibrations of the UV and FIR emission. The median sSFR is $18_{-5}^{+7}$ Gyr$^{-1}$, significantly larger than literature measurements lacking constraints in the FIR. The increase in sSFR reflects the larger obscured SFRs we derive from the dust continuum relative to that implied by the UV+optical SED. We suggest that such differences may reflect spatial variations in dust across these luminous galaxies, with the component dominating the FIR distinct from that dominating the UV. We demonstrate that the inferred stellar masses (and hence sSFRs) are strongly-dependent on the assumed star formation history (SFH) in reionization-era galaxies. When large sSFR galaxies are modeled with non-parametric SFHs, the derived stellar masses can increase by an order of magnitude relative to constant star formation models, owing to the presence of a significant old stellar population that is outshined by the recent burst. The [CII] line widths in the largest sSFR systems are often very broad, suggesting dynamical masses that are easily able to accommodate the dominant old stellar population suggested by non-parametric models. Regardless of these systematic uncertainties in the derived parameters, we find that the sSFR increases rapidly toward higher redshifts for massive galaxies ($9.6<\log(\rm M_*/M_{\odot})<9.8$), with a power law that goes as $(1+z)^{1.7\pm0.3}$, broadly consistent with expectations from the evolving baryon accretion rates.

preprint2021arXiv

Deep Realistic Extragalactic Model (DREaM) Galaxy Catalogs: Predictions for a Roman Ultra-Deep Field

In the next decade, deep galaxy surveys from telescopes such as the James Webb Space Telescope and Roman Space Telescope will provide transformational data sets that will greatly enhance the understanding of galaxy formation during the epoch of reionization (EoR). In this work, we present the Deep Realistic Extragalactic Model (DREaM) for creating synthetic galaxy catalogs. Our model combines dark matter simulations, subhalo abundance matching and empirical models, and includes galaxy positions, morphologies, and spectral energy distributions (SEDs). The resulting synthetic catalog extends to redshifts $z \sim 12 $, and galaxy masses $\log_{10}(M/M_{\odot}) = 5 $ covering an area of $1\, {\rm deg}^2$ on the sky. We use DREaM to explore the science returns of a $1\, {\rm deg}^2$ \emph{Roman} UDF, and to provide a resource for optimizing ultra-deep survey designs. We find that a \emph{Roman} UDF to $\sim 30\,m_{\rm AB}$ will potentially detect more than $10^6$ $M_{\rm UV}<-17$ galaxies, with more than $10^4$ at redshifts $z>7$, offering an unparalleled dataset for constraining galaxy properties during the EoR. Our synthetic catalogs and simulated images are made publicly available to provide the community with a tool to prepare for upcoming data.

preprint2021arXiv

Direct Constraints on the Extremely Metal-Poor Massive Stars Underlying Nebular C IV Emission from Ultra-Deep HST/COS Ultraviolet Spectroscopy

Metal-poor nearby galaxies hosting massive stars have a fundamental role to play in our understanding of both high-redshift galaxies and low metallicity stellar populations. But while much attention has been focused on their bright nebular gas emission, the massive stars that power it remain challenging to constrain. Here we present exceptionally deep Hubble Space Telescope ultraviolet spectra targeting six galaxies that power strong nebular C IV emission approaching that encountered at $z>6$. We find that the strength and spectral profile of the nebular C IV in these new spectra follow a sequence evocative of resonant scattering models, indicating that the hot circumgalactic medium likely plays a key role in regulating C IV escape locally. We constrain the metallicity of the massive stars in each galaxy by fitting the forest of photospheric absorption lines, reporting measurements driven by iron that lie uniformly below 10% solar. Comparison with the gas-phase oxygen abundances reveals evidence for enhancement in O/Fe above solar across the sample, robust to assumptions about the absolute gas-phase metallicity scale. This supports the idea that these local systems are more chemically-similar to their primordial high-redshift counterparts than to the bulk of nearby galaxies. Finally, we find significant tension between the strong stellar wind profiles observed and our population synthesis models constrained by the photospheric forest in our highest-quality spectra. This reinforces the need for caution in interpreting wind lines in isolation at high-redshift, but also suggests a unique path towards validating fundamental massive star physics at extremely low metallicity with integrated ultraviolet spectra.

preprint2021arXiv

The [OIII]$\lambda5007$ equivalent width distribution at z $\sim2$: The redshift evolution of the extreme emission line galaxies

We determine the [OIII]$\lambda5007$ equivalent width (EW) distribution of $1.700<\rm{z}<2.274$ rest-frame UV-selected (M$_{\rm{UV}}<-19$) star-forming galaxies in the GOODS North and South fields. We make use of deep HDUV broadband photometry catalogues for selection and 3D-HST WFC3/IR grism spectra for measurement of line properties. The [OIII]$\lambda5007$ EW distribution allows us to measure the abundance of extreme emission line galaxies (EELGs) within this population. We model a log-normal distribution to the [OIII]$\lambda5007$ rest-frame equivalent widths of galaxies in our sample, with location parameter $μ=4.24\pm0.07$ and variance parameter $σ= 1.33\pm0.06$. This EW distribution has a mean [OIII]$\lambda5007$ EW of 168$\pm1Å$. The fractions of $\rm{z}\sim2$ rest-UV-selected galaxies with [OIII]$\lambda5007$ EWs greater than $500, 750$ and $1000Å$ are measured to be $6.8^{+1.0}_{-0.9}\%$, $3.6^{+0.7}_{-0.6}\%$, and $2.2^{+0.5}_{-0.4}\%$ respectively. The EELG fractions do not vary strongly with UV luminosity in the range ($-21.6<M_{\rm{UV}}<-19.0$) considered in this paper, consistent with findings at higher redshifts. We compare our results to $\rm{z}\sim5$ and $\rm{z}\sim7$ studies where candidate EELGs have been discovered through Spitzer/IRAC colours, and we identify rapid evolution with redshift in the fraction of star-forming galaxies observed in an extreme emission line phase (a rise by a factor $\sim10$ between $\rm{z}\sim2$ and $\rm{z}\sim7$). This evolution is consistent with an increased incidence of strong bursts in the galaxy population of the reionisation era. While this population makes a sub-dominant contribution of the ionising emissivity at $\rm{z}\simeq2$, EELGs are likely to dominate the ionising output in the reionisation era.

preprint2020arXiv

Effects of self-consistent rest-ultraviolet colours in semi-empirical galaxy formation models

Connecting the observed rest-ultraviolet (UV) luminosities of high-$z$ galaxies to their intrinsic luminosities (and thus star formation rates) requires correcting for the presence of dust. We bypass a common dust-correction approach that uses empirical relationships between infrared (IR) emission and UV colours, and instead augment a semi-empirical model for galaxy formation with a simple -- but self-consistent -- dust model and use it to jointly fit high-$z$ rest-UV luminosity functions (LFs) and colour-magnitude relations ($M_{\mathrm{UV}}$-$β$). In doing so, we find that UV colours evolve with redshift (at fixed UV magnitude), as suggested by observations, even in cases without underlying evolution in dust production, destruction, absorption, or geometry. The observed evolution in our model arises due to the reduction in the mean stellar age and rise in specific star formation rates with increasing $z$. The UV extinction, $A_{\mathrm{UV}}$, evolves similarly with redshift, though we find a systematically shallower relation between $A_{\mathrm{UV}}$ and $M_{\mathrm{UV}}$ than that predicted by IRX-$β$ relationships derived from $z \sim 3$ galaxy samples. Finally, assuming that high $1600 Å$ transmission ($\gtrsim 0.6$) is a reliable LAE indicator, modest scatter in the effective dust surface density of galaxies can explain the evolution both in $M_{\mathrm{UV}}$-$β$ and LAE fractions. These predictions are readily testable by deep surveys with the James Webb Space Telescope.

preprint2020arXiv

Exploring Chemical Homogeneity in Dwarf Galaxies: A VLT-MUSE study of JKB 18

Deciphering the distribution of metals throughout galaxies is fundamental in our understanding of galaxy evolution. Nearby, low-metallicity, star-forming dwarf galaxies in particular can offer detailed insight into the metal-dependent processes that may have occurred within galaxies in the early Universe. Here we present VLT/MUSE observations of one such system, JKB 18, a blue diffuse dwarf galaxy with a metallicity of only 12+log(O/H)=7.6$\pm$0.2 (~0.08 Zsol). Using high spatial-resolution integral field spectroscopy of the entire system, we calculate chemical abundances for individual HII regions using the direct method and derive oxygen abundance maps using strong-line metallicity diagnostics. With large-scale dispersions in O/H, N/H and N/O of ~0.5-0.6 dex and regions harbouring chemical abundances outside this 1$σ$ distribution, we deem JKB 18 to be chemically inhomogeneous. We explore this finding in the context of other chemically inhomogeneous dwarf galaxies and conclude that neither the accretion of metal-poor gas, short mixing timescales, or self-enrichment from Wolf-Rayet stars are accountable. Using a galaxy-scale, multi-phase, hydrodynamical simulation of a low-mass dwarf galaxy, we find that chemical inhomogeneities of this level may be attributable to the removal of gas via supernovae and the specific timing of the observations with respect to star-formation activity. This study not only draws attention to the fact that dwarf galaxies can be chemically inhomogeneous, but also that the methods used in the assessment of this characteristic can be subject to bias.

preprint2020arXiv

High-mass X-ray binaries in nearby metal-poor galaxies: on the contribution to nebular HeII emission

Despite significant progress both observationally and theoretically, the origin of high-ionization nebular HeII emission in galaxies dominated by stellar photoionization remains unclear. Accretion-powered radiation from high-mass X-ray binaries (HMXBs) is still one of the leading proposed explanations for the missing $\mathrm{He^+}$-ionizing photons, but this scenario has yet to be conclusively tested. In this paper, we present nebular line predictions from a grid of photoionization models with input SEDs containing the joint contribution of both stellar atmospheres and a multi-color disk model for HMXBs. This grid demonstrates that HMXBs are inefficient producers of the photons necessary to power HeII, and can only boost this line substantially in galaxies with HMXB populations large enough to power X-ray luminosities of $10^{42}$ erg/s per unit star formation rate (SFR). To test this, we assemble a sample of eleven low-redshift star-forming galaxies with high-quality constraints on both X-ray emission from Chandra and HeII emission from deep optical spectra, including new observations with the MMT. These data reveal that the HMXB populations of these nearby systems are insufficient to account for the observed HeII strengths, with typical X-ray luminosities or upper limits thereon of only $10^{40}$-$10^{41}$ erg/s per SFR. This indicates that HMXBs are not the dominant source of $\mathrm{He^+}$ ionization in these metal-poor star-forming galaxies. We suggest that the solution may instead reside in revisions to stellar wind predictions, softer X-ray sources, or very hot products of binary evolution at low metallicity.

preprint2020arXiv

RELICS: A Very Large ($θ_{E}\sim40&#34;$) Cluster Lens -- RXC J0032.1+1808

Extensive surveys with the \textit{Hubble Space Telescope} (HST) over the past decade, targeting some of the most massive clusters in the sky, have uncovered dozens of galaxy-cluster strong lenses. The massive cluster strong-lens scale is typically $θ_{E}\sim10\arcsec$ to $\sim30-35\arcsec$, with only a handful of clusters known with Einstein radii $θ_{E}\sim40\arcsec$ or above (for $z_{source}=2$, nominally). Here we report another very large cluster lens, RXC J0032.1+1808 ($z=0.3956$), the second richest cluster in the redMapper cluster catalog and the 85th most massive cluster in the Planck Sunyaev-Zel&#39;dovich catalog. With our Light-Traces-Mass and fully parametric (dPIEeNFW) approaches, we construct strong lensing models based on 18 multiple images of 5 background galaxies newly identified in the \textit{Hubble} data mainly from the \textit{Reionization Lensing Cluster Survey} (RELICS), in addition to a known sextuply imaged system in this cluster. Furthermore, we compare these models to Lenstool and GLAFIC models that were produced independently as part of the RELICS program. All models reveal a large effective Einstein radius of $θ_{E}\simeq40\arcsec$ ($z_{source}=2$), owing to the obvious concentration of substructures near the cluster center. Although RXC J0032.1+1808 has a very large critical area and high lensing strength, only three magnified high-redshift candidates are found within the field targeted by RELICS. Nevertheless, we expect many more high-redshift candidates will be seen in wider and deeper observations with \textit{Hubble} or \emph{JWST}. Finally, the comparison between several algorithms demonstrates that the total error budget is largely dominated by systematic uncertainties.

preprint2020arXiv

Searching for z > 6.5 Analogs Near the Peak of Cosmic Star Formation

Strong [OIII]$λλ$4959,5007+H$β$ emission appears to be typical in star-forming galaxies at z > 6.5. As likely contributors to cosmic reionization, these galaxies and the physical conditions within them are of great interest. At z > 6.5, where Ly$α$ is greatly attenuated by the intergalactic medium, rest-UV metal emission lines provide an alternative measure of redshift and also constraints on the physical properties of star-forming regions and massive stars. We present the first statistical sample of rest-UV line measurements in z $\sim$ 2 galaxies selected as analogs of those in the reionization era based on [OIII]$λλ$4959,5007 EW or rest-frame U-B color. Our sample is drawn from the 3D-HST Survey and spans the redshift range 1.36 $\leqslant$ z $\leqslant$ 2.49. We find that the median Ly$α$ and CIII]$λλ$1907,1909 EWs of our sample are significantly greater than those of z $\sim$ 2 UV-continuum-selected star-forming galaxies. Measurements from both individual and composite spectra indicate a monotonic, positive correlation between CIII] and [OIII], while a lack of trend is observed between Ly$α$ and [OIII] at [OIII] EW < 1000$\unicode{x212B}$. At higher [OIII] EW, extreme Ly$α$ emission starts to emerge. Using stacked spectra, we find that Ly$α$ and CIII] are significantly enhanced in galaxies with lower metallicity. Two objects in our sample appear comparable to z > 6.5 galaxies with exceptionally strong rest-UV metal line emission. These objects have significant CIV$λλ$1548,1550, HeII$λ$1640, and OIII]$λλ$1661,1665 emission in addition to intense Ly$α$ or CIII]. Detailed characterization of these lower-redshift analogs provides unique insights into the physical conditions in z > 6.5 star-forming regions, motivating future observations of reionization-era analogs at lower redshifts.

preprint2020arXiv

The Ionizing Photon Production Efficiency ($ξ_{ion}$) Of Lensed Dwarf Galaxies At $z \sim 2 $

We measure the ionizing photon production efficiency ($ξ_{ion}$) of low-mass galaxies ($10^{7.8}$-$10^{9.8}$ $M_{\odot}$) at $1.4<z<2.7$, allowing us to better understand the contribution of dwarf galaxies to the ionizing background and cosmic reionization. We target galaxies that are magnified by the strong lensing galaxy clusters Abell 1689, MACS J0717, and MACS J1149. We utilize Keck/MOSFIRE spectra to measure optical nebular emission line fluxes and HST imaging to measure the rest-UV and rest-optical photometry. We present two methods of stacking. First, we take the average of the log(L$_{Hα}$ /L$_{UV}$) of galaxies in our sample to determine the typical log($ξ_{ion}$). Second, we take the logarithm of the total L$_{Hα}$ over the total L$_{UV}$. We prefer the latter as it provides the total ionizing UV luminosity density of galaxies when multiplied by the non-ionizing UV luminosity density from the UV luminosity function. log($ξ_{ion}$) calculated from the second method is $\sim$ 0.2 dex higher than the first method. We do not find any strong dependence between log($ξ_{ion}$) and stellar mass, M$_{UV}$ or UV spectral slope ($β$). We report a value of log($ξ_{ion}$) $\sim25.47\pm 0.09$ for our UV-complete sample ($-22<M_{UV}<-17.3$) and $\sim25.37\pm0.11$ for our mass-complete sample ($7.8<\log(M_*)<9.8)$. These values are consistent with measurements of more massive, more luminous galaxies in other high-redshift studies that use the same stacking technique. Our log($ξ_{ion}$) is $0.2-0.3$ dex higher than low-redshift galaxies of similar mass, indicating an evolution in the stellar properties, possibly due to metallicity, age, or the prevalence of binary stars. We also find a correlation between log($ξ_{ion}$) and the equivalent widths of H$α$ and [OIII]$λ$5007 fluxes, confirming that these equivalent widths can be used to estimate $ξ_{ion}$.

preprint2019arXiv

Clustering with JWST: Constraining Galaxy Host Halo Masses, Satellite Quenching Efficiencies, and Merger Rates at z=4-10

Galaxy clustering measurements can be used to constrain many aspects of galaxy evolution, including galaxy host halo masses, satellite quenching efficiencies, and merger rates. We simulate JWST galaxy clustering measurements at z$\sim$4$-$10 by utilizing mock galaxy samples produced by an empirical model, the UniverseMachine. We also adopt the survey footprints and typical depths of the planned joint NIRCam and NIRSpec Guaranteed Time Observation program planned for Cycle 1 to generate realistic JWST survey realizations and to model high-redshift galaxy selection completeness. We find that galaxy clustering will be measured with $\gtrsim$5$σ$ significance at z$\sim$4$-$10. Halo mass precisions resulting from Cycle 1 angular clustering measurements will be $\sim$0.2 dex for faint (-18 $\gtrsim$ $\mathit{M}_{\mathrm{UV}}^{ }$ $\gtrsim$ -19) galaxies at z$\sim$4$-$10 as well as $\sim$0.3 dex for bright ($\mathit{M}_{\mathrm{UV}}^{ }$ $\sim$ -20) galaxies at z$\sim$4$-$7. Dedicated spectroscopic follow-up over $\sim$150 arcmin$^2$ would improve these precisions by $\sim$0.1 dex by removing chance projections and low-redshift contaminants. Future JWST observations will therefore provide the first constraints on the stellar-halo mass relation in the epoch of reionization and substantially clarify how this relation evolves at z$>$4. We also find that $\sim$1000 individual satellites will be identifiable at z$\sim$4$-$8 with JWST, enabling strong tests of satellite quenching evolution beyond currently available data (z$\lesssim$2). Finally, we find that JWST observations can measure the evolution of galaxy major merger pair fractions at z$\sim$4$-$8 with $\sim$0.1$-$0.2 dex uncertainties. Such measurements would help determine the relative role of mergers to the build-up of stellar mass into the epoch of reionization.

preprint2019arXiv

The Detection of [O III]$λ$4363 in a Lensed, Dwarf Galaxy at $z$ = 2.59: Testing Metallicity Indicators and Scaling Relations at High Redshift and Low Mass

We present Keck/MOSFIRE and Keck/LRIS spectroscopy of A1689-217, a lensed (magnification $\sim7.9$), star-forming (SFR $\sim$ 16 $\rm M_\odot\ yr^{-1}$), dwarf (log($M_\ast/M_\odot$) = $8.07-8.59$), Ly$α$-emitter ($\rm EW_0\sim138$ A) at $z$ = 2.5918. Dwarf galaxies similar to A1689-217 are common at high redshift and likely responsible for reionization, yet few have been studied with detailed spectroscopy. We report a 4.2$σ$ detection of the electron-temperature-sensitive [O III]$λ$4363 emission line and use this line to directly measure an oxygen abundance of 12+$\log$(O/H) = 8.06 $\pm$ 0.12 ($\sim1/4\ \rm Z_\odot$). A1689-217 is the lowest mass galaxy at $z>$ 2 with an [O III]$λ$4363 detection. Using the rest-optical emission lines, we measure A1689-217&#39;s other nebular conditions including electron temperature ($T_e$([O III]) $\sim$ 14,000 K), electron density ($n_e\sim$ 220 $\rm cm^{-3}$) and reddening ($E(B-V)\sim$ 0.39). We study relations between strong-line ratios and direct metallicities with A1689-217 and other galaxies with [O III]$λ$4363 detections at $z\sim0-3.1$, showing that the locally-calibrated, oxygen-based, strong-line relations are consistent from $z\sim0-3.1$. We also show additional evidence that the $\rm O_{32}$ vs. $\rm R_{23}$ excitation diagram can be utilized as a redshift-invariant, direct-metallicity-based, oxygen abundance diagnostic out to $z\sim$ 3.1. From this excitation diagram and the strong-line ratio $-$ metallicity plots, we observe that the ionization parameter at fixed O/H is consistent with no redshift evolution. Although A1689-217 is metal-rich for its $M_\ast$ and SFR, we find it to be consistent within the large scatter of the low-mass end of the Fundamental Metallicity Relation.

preprint2019arXiv

The Lyman Continuum Escape Survey -- II: Ionizing Radiation as a Function of the [OIII]/[OII] Line Ratio

We discuss the rest-frame optical emission line spectra of a large (~50) sample of z=3.1 Lyman alpha emitters (LAEs) whose physical properties suggest such sources are promising analogs of galaxies in the reionization era. Reliable Lyman continuum escape fractions have now been determined for a large sample of such LAEs from the Lyman Continuum Escape Survey (LACES) undertaken via deep HST imaging in the SSA22 survey area reported in Fletcher et al. (2019). Using new measures of [OII] emission secured from Keck MOSFIRE spectra we re-examine, for a larger sample, earlier claims that Lyman continuum leakages may correlate with the nebular emission line ratio [OIII]/[OII] as expected for density-bound HII regions. We find that a large [OIII]/[OII] line ratio is indeed a necessary condition for Lyman continuum leakage, strengthening earlier claims made using smaller samples at various redshifts. However, not all LAEs with large [OIII]/[OII] line ratios are leakers and leaking radiation appears not to be associated with differences in other spectral diagnostics. This suggests the detection of leaking radiation is modulated by an additional property, possibly the viewing angle for porous HII regions. We discuss our new results in the context of the striking bimodality of LAE leakers and non-leakers found in the LACES program and the implications for the sources of cosmic reionization.