Source author record

Daniel P. Stark

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

57works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

57 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'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' 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'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' 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"$) 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'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'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.

preprint2016arXiv

A Keck Survey of Gravitationally-Lensed Star-Forming Galaxies: High Spatial Resolution Studies of Kinematics and Metallicity Gradients

We discuss spatially resolved emission line spectroscopy secured for a total sample of 15 gravitationally lensed star-forming galaxies at a mean redshift of $z\simeq2$ based on Keck laser-assisted adaptive optics observations undertaken with the recently-improved OSIRIS integral field unit (IFU) spectrograph. By exploiting gravitationally lensed sources drawn primarily from the CASSOWARY survey, we sample these sub-L$^{\ast}$ galaxies with source-plane resolutions of a few hundred parsecs ensuring well-sampled 2-D velocity data and resolved variations in the gas-phase metallicity. Such high spatial resolution data offers a critical check on the structural properties of larger samples derived with coarser sampling using multiple-IFU instruments. We demonstrate how kinematic complexities essential to understanding the maturity of an early star-forming galaxy can often only be revealed with better sampled data. Although we include four sources from our earlier work, the present study provides a more representative sample unbiased with respect to emission line strength. Contrary to earlier suggestions, our data indicates a more diverse range of kinematic and metal gradient behavior inconsistent with a simple picture of well-ordered rotation developing concurrently with established steep metal gradients in all but merging systems. Comparing our observations with the predictions of hydrodynamical simulations suggests that gas and metals have been mixed by outflows or other strong feedback processes, flattening the metal gradients in early star-forming galaxies.

preprint2016arXiv

Absorption Line Spectroscopy of Gravitationally-Lensed Galaxies: Further Constraints on the Escape Fraction of Ionizing Photons at High Redshift

The fraction of ionizing photons escaping from high-redshift star-forming galaxies remains a key obstacle in evaluating whether galaxies were the primary agents of cosmic reionization. We previously proposed using the covering fraction of low-ionization gas, measured via deep absorption line spectroscopy, as a proxy. We now present a significant update, sampling seven gravitationally-lensed sources at $4<z<5$. We show that the absorbing gas in our sources is spatially inhomogeneous with a median covering fraction of 66\%. Correcting for reddening according to a dust-in-cloud model, this implies an estimated absolute escape fraction of $\simeq19\pm6$\%. With possible biases and uncertainties, collectively we find the average escape fraction could be reduced to no less than 11\%, excluding the effect of spatial variations. For one of our lensed sources, we have sufficient signal/noise to demonstrate the presence of such spatial variations and scatter in its dependence on the Ly$α$ equivalent width, consistent with recent simulations. If this source is typical, our lower limit to the escape fraction could be reduced by a further factor $\simeq$2. Across our sample, we find a modest anti-correlation between the inferred escape fraction and the local star formation rate, consistent with a time delay between a burst and leaking Lyman continuum photons. Our analysis demonstrates considerable variations in the escape fraction consistent with being governed by the small-scale behavior of star-forming regions, whose activities fluctuate over short timescales. This supports the suggestion that the escape fraction may increase toward the reionization era when star formation becomes more energetic and burst-like.

preprint2016arXiv

Lyman-alpha and CIII] Emission in z=7-9 Galaxies: Accelerated Reionization Around Luminous Star Forming Systems?

We discuss new Keck/MOSFIRE spectroscopic observations of four luminous galaxies at z~7-9 selected to have intense optical line emission by Roberts-Borsani et al. (2016). Previous follow-up has revealed Lyman-alpha in two of the four galaxies. Our new MOSFIRE observations confirm that Lyman-alpha is present in the entire sample. We detect Lyman-alpha emission in COS-zs7-1, confirming its redshift as z=7.154, and we detect Lyman-alpha in EGS-zs8-2 at z=7.477, verifying a tentative detection presented in an earlier study. The ubiquity of Lyman-alpha in this sample is puzzling given that the IGM is likely significantly neutral over 7<z<9. To investigate this result in more detail, we have initiated a campaign to target UV metal emission in the four Lyman-alpha emitters as a probe of both the radiation field and the velocity offset of Lyman-alpha. Here we present the detection of intense CIII] emission in EGS-zs8-1, a galaxy from this sample previously shown to have Lyman-alpha at z=7.73. Photoionization models indicate that an intense radiation field and low metallicity are required to reproduce the intense CIII] and optical line emission. We argue that this extreme radiation field is likely to affect the local environment, increasing the transmission of Lyman-alpha through the galaxy. Moreover, the centroid of CIII] indicates that Lyman-alpha is redshifted from the systemic value by 340 km/s. This velocity offset is larger than that seen in less luminous systems, providing an additional explanation for the transmission of Lyman-alpha emission through the IGM. Since the transmission is further enhanced by the likelihood that such systems are also situated in the densest regions with the largest ionized bubbles, the visibility of Lyman-alpha at z>7 is expected to be strongly luminosity-dependent, with the most effective transmission occurring in systems with intense star formation.

preprint2016arXiv

The Evolution Of The Faint End Of The UV Luminosity Function During The Peak Epoch Of Star Formation (1<z<3)

[Abridged] We present a robust measurement of the rest-frame UV luminosity function (LF) and its evolution during the peak epoch of cosmic star formation at 1<z<3. We use our deep near ultraviolet imaging from WFC3/UVIS on the Hubble Space Telescope (HST) and existing ACS/WFC and WFC3/IR imaging of three lensing galaxy clusters, Abell 2744 and MACSJ0717 from the Hubble Frontier Field survey and Abell 1689. We use photometric redshifts to identify 780 ultra-faint galaxies with $M_{UV}$<-12.5 AB mag at 1<z<3. From these samples, we identified 5 new, faint, multiply imaged systems in A1689. We compute the rest-frame UV LF and find the best-fit faint-end slopes of $α=-1.56\pm0.04$, $α=-1.72\pm0.04$ and $α=-1.94\pm0.06$ at 1.0<z<1.6, 1.6<z<2.2 and 2.2<z<3.0, respectively. Our results demonstrate that the UV LF becomes steeper from z\sim1.3 to z\sim2.6 with no sign of a turnover down to $M_{UV}=-14$ AB mag. We further derive the UV LFs using the Lyman break "dropout" selection and confirm the robustness of our conclusions against different selection methodologies. Because the sample sizes are so large, and extend to such faint luminosities, the statistical uncertainties are quite small, and systematic uncertainties (due to the assumed size distribution, for example), likely dominate. If we restrict our analysis to galaxies and volumes above > 50% completeness in order to minimize these systematics, we still find that the faint-end slope is steep and getting steeper with redshift, though with slightly shallower (less negative) values ($α=-1.55\pm0.06$, $-1.69\pm0.07$ and $-1.79\pm0.08$ for $z\sim1.3$, 1.9 and 2.6, respectively). Finally, we conclude that the faint star-forming galaxies with UV magnitudes of $-18.5<M_{UV}<-12.5$ covered in this study, produce the majority (55%-60%) of the unobscured UV luminosity density at 1<z<3.

preprint2015arXiv

Consequences of bursty star formation on galaxy observables at high redshifts

The star formation histories (SFHs) of dwarf galaxies are thought to be \emph{bursty}, with large -- order of magnitude -- changes in the star formation rate on timescales similar to O-star lifetimes. As a result, the standard interpretations of many galaxy observables (which assume a slowly varying SFH) are often incorrect. Here, we use the SFHs from hydro-dynamical simulations to investigate the effects of bursty SFHs on sample selection and interpretation of observables and make predictions to confirm such SFHs in future surveys. First, because dwarf galaxies' star formation rates change rapidly, the mass-to-light ratio is also changing rapidly in both the ionizing continuum and, to a lesser extent, the non-ionizing UV continuum. Therefore, flux limited surveys are highly biased toward selecting galaxies in the \emph{burst} phase and very deep observations are required to detect all dwarf galaxies at a given stellar mass. Second, we show that a $\log_{10}[νL_ν(1500{\rm Å})/L_{{\rm H}α}]>2.5$ implies a very recent quenching of star formation and can be used as evidence of stellar feedback regulating star formation. Third, we show that the ionizing continuum can be significantly higher than when assuming a constant SFH, which can affect the interpretation of nebular emission line equivalent widths and direct ionizing continuum detections. Finally, we show that a star formation rate estimate based on continuum measurements only (and not on nebular tracers such as the hydrogen Balmer lines) will not trace the rapid changes in star formation and will give the false impression of a star-forming main sequence with low dispersion.

preprint2015arXiv

Lyman-alpha Emission from a Luminous z=8.68 Galaxy: Implications for Galaxies as Tracers of Cosmic Reionization

We report the discovery of Lyman-alpha emission (Ly$α$) in the bright galaxy EGSY-2008532660 (hereafter EGSY8p7) using the MOSFIRE spectrograph at the Keck Observatory. First reported by Roberts-Borsani et al. (2015), it was selected for spectroscopic observations because of its photometric redshift ($z_{phot}=8.57^{+0.22}_{-0.43}$), apparent brightness (H$_{160}=25.26\pm0.09$) and red Spitzer/IRAC [3.6]-[4.5] color indicative of contamination by strong oxygen emission in the [4.5] band. With a total integration of $\sim$4.3 hours, our data reveal an emission line at $\simeq$11776 Å which we argue is likely Ly$α$ at a redshift $z_{spec}=8.683^{+0.001}_{-0.004}$, in good agreement with the photometric estimate. The line was detected independently on two nights using different slit orientations and its detection significance is $\sim7.5σ$. An overlapping skyline contributes significantly to the uncertainty on the total line flux although the significance of the detected line is robust to a variety of skyline-masking procedures. By direct addition and a Gaussian fit, we estimate a 95\% confidence range of 1.0--2.5$\times10^{-17}$ erg s$^{-1}$ cm$^{-2}$, corresponding to a rest-frame equivalent width of 17--42 Å. EGSY8p7 is the most distant galaxy confirmed spectroscopically to date, and the third luminous source in the EGS field beyond $z_{phot}\gtrsim7.5$ with detectable Ly$α$ emission viewed at a time when the intergalactic medium is believed to be fairly neutral. Although the reionization process was probably patchy, we discuss whether luminous sources with prominent IRAC color excesses may harbor harder ionizing spectra than the dominant fainter population thereby creating earlier ionized bubbles. Further spectroscopic follow-up of such bright sources promises important insight into the early formation of galaxies.

preprint2015arXiv

Spectroscopic detection of CIV in a galaxy at z=7.045: Implications for the ionizing spectra of reionization-era galaxies

We present Keck/MOSFIRE observations of UV metal lines in four bright gravitationally-lensed z~6-8 galaxies behind the cluster Abell 1703. The spectrum of A1703-zd6, a highly-magnified star forming galaxy with a Lyman-alpha redshift of z=7.045, reveals a confident detection of the nebular CIV emission line (unresolved with FWHM < 125 km/s). UV metal lines are not detected in the three other galaxies. At z~2-3, nebular CIV emission is observed in just 1% of UV-selected galaxies. The presence of strong CIV emission in one of the small sample of galaxies targeted in this paper may indicate hard ionizing spectra are more common at z~7. The total estimated equivalent width of the CIV doublet (38 A) and CIV/Lyman-alpha flux ratio (0.3) are comparable to measurements of narrow-lined AGNs. Photoionization models show that the nebular CIV line can also be reproduced by a young stellar population, with very hot metal poor stars dominating the photon flux responsible for triply ionizing carbon. Regardless of the origin of the CIV, we show that the ionizing spectrum of A1703-zd6 is different from that of typical galaxies at z~2, producing more H ionizing photons per unit 1500A luminosity and a larger flux density at 30-50 eV. If such extreme radiation fields are typical in UV-selected systems at z>7, it would indicate that reionization-era galaxies are more efficient ionizing agents than previously thought. Alternatively, we suggest that the small sample of Lyman-alpha emitters at z>7 may trace a rare population with intense radiation fields capable of ionizing their surrounding hydrogen distribution. Additional constraints on high ionization emission lines in galaxies with and without Lyman-alpha detections will help clarify whether hard ionizing spectra are common in the reionization era.

preprint2015arXiv

Uncovering Blue Diffuse Dwarf Galaxies

Extremely metal poor (XMP) galaxies are known to be very rare, despite the large numbers of low-mass galaxies predicted by the local galaxy luminosity function. This paper presents a sub-sample of galaxies that were selected via a morphology-based search on SDSS images with the aim of finding these elusive XMP galaxies. By using the recently discovered extremely metal-poor galaxy, Leo P, as a guide, we obtained a collection of faint, blue systems, each with isolated HII regions embedded in a diffuse continuum, that have remained optically undetected until now. Here we show the first results from optical spectroscopic follow-up observations of 12 of ~100 of these blue, diffuse dwarf (BDD) galaxies yielded by our search algorithm. Oxygen abundances were obtained via the direct method for eight galaxies, and found to be in the range 7.45<12+log(O/H)<8.0, with two galaxies being classified as XMPs. All BDDs were found to currently have a young star-forming population (<10 Myr) and relatively high ionisation parameters of their HII regions. Despite their low luminosities (-11< M_B < -18) and low surface brightnesses (~23-25 mag arcsec^-2), the galaxies were found to be actively star-forming, with current star-formation rates between 0.0003 and 0.078 Msol/yr. From our current subsample, BDD galaxies appear to be a population of non-quiescent dwarf irregular (dIrr) galaxies, or the diffuse counterparts to blue compact galaxies (BCDs) and as such may bridge the gap between these two populations. Our search algorithm demonstrates that morphology-based searches are successful in uncovering more diffuse metal-poor star-forming galaxies, which traditional emission-line based searches overlook.

preprint2014arXiv

Accounting for Cosmic Variance in Studies of Gravitationally-Lensed High-Redshift Galaxies in the Hubble Frontier Field Clusters

Strong gravitational lensing provides a powerful means for studying faint galaxies in the distant universe. By magnifying the apparent brightness of background sources, massive clusters enable the detection of galaxies fainter than the usual sensitivity limit for blank fields. However, this gain in effective sensitivity comes at the cost of a reduced survey volume and, in this {\it Letter}, we demonstrate there is an associated increase in the cosmic variance uncertainty. As an example, we show that the cosmic variance uncertainty of the high redshift population viewed through the Hubble Space Telescope Frontier Field cluster Abell 2744 increases from ~35% at redshift z~7 to >~65% at z~10. Previous studies of high redshift galaxies identified in the Frontier Fields have underestimated the cosmic variance uncertainty that will affect the ultimate constraints on both the faint end slope of the high-redshift luminosity function and the cosmic star formation rate density, key goals of the Frontier Field program.

preprint2014arXiv

Line Emitting Galaxies Beyond a Redshift of 7: An Improved Method for Estimating the Evolving Neutrality of the Intergalactic Medium

The redshift-dependent fraction of color-selected galaxies revealing Lyman alpha emission has become the most valuable constraint on the evolving neutrality of the early intergalactic medium. However, in addition to resonant scattering by neutral gas, the visibility of Lyman alpha is also dependent on the intrinsic properties of the host galaxy, including its stellar population, dust content and the nature of outflowing gas. Taking advantage of significant progress we have made in determining the line emitting properties of $z \simeq 4-6$ galaxies, we propose an improved method, based on using the measured slopes of the rest-frame ultraviolet continua of galaxies, to interpret the growing body of near-infrared spectra of $z>7$ galaxies in order to take into account these host galaxy dependencies. In a first application of our new method, we demonstrate its potential via a new spectroscopic survey of $7<z<8$ galaxies undertaken with the Keck MOSFIRE spectrograph. Together with earlier published data our data provides improved estimates of the evolving visibility of Lyman alpha, particularly at redshift $z\simeq 8$. As a byproduct, we also present a new line emitting galaxy at a redshift $z=7.62$ which supersedes an earlier redshift record. We discuss the improving constraints on the evolving neutral fraction over $6<z<8$ and the implications for cosmic reionization.

preprint2014arXiv

Spectroscopic detections of CIII]1909 at z~6-7: A new probe of early star forming galaxies and cosmic reionisation

Deep spectroscopic observations of z~6.5 galaxies have revealed a marked decline with increasing redshift in the detectability of Lyman-alpha emission. While this may offer valuable insight into the end of the reionisation process, it presents a fundamental challenge to the detailed spectroscopic study of the many hundreds of photometrically-selected distant sources now being found via deep HST imaging, and particularly those bright sources viewed through foreground lensing clusters. In this paper we demonstrate the validity of a new way forward via the convincing detection of an alternative diagnostic line, CIII]1909, seen in spectroscopic exposures of two star forming galaxies at z=6.029 and 7.213. The former detection is based on a 3.5 hour X-shooter spectrum of a bright (J=25.2) gravitationally-lensed galaxy behind the cluster Abell 383. The latter detection is based on a 4.2 hour MOSFIRE spectra of one of the most distant spectroscopically confirmed galaxies, GN-108036, with J=25.2. Both targets were chosen for their continuum brightness and previously-known redshift (based on Lyman-alpha), ensuring that any CIII] emission would be located in a favorable portion of the near-infrared sky spectrum. We compare our CIII] and Lyman-alpha equivalent widths in the context of those found at z~2 from earlier work and discuss the motivation for using lines other than Lyman-alpha to study galaxies in the reionisation era.

preprint2014arXiv

Ultraviolet Emission Lines in Young Low Mass Galaxies at z~2: Physical Properties and Implications for Studies at z>7

We present deep spectroscopy of 17 very low mass (M* ~ 2.0x10^6 Msun to 1.4x10^9 Msun) and low luminosity (M_UV ~ -13.7 to -19.9) gravitationally lensed galaxies in the redshift range z~1.5-3.0. Deep rest-frame ultraviolet spectra reveal large equivalent width emission from numerous lines (NIV], OIII], CIV, Si III], CIII]) which are rarely seen in individual spectra of more massive star forming galaxies. CIII] is detected in 16 of 17 low mass star forming systems with rest-frame equivalent widths as large as 13.5 Angstroms. Nebular CIV emission is present in the most extreme CIII] emitters, requiring an ionizing source capable of producing a substantial component of photons with energies in excess of 47.9 eV. Photoionization models support a picture whereby the large equivalent widths are driven by the increased electron temperature and enhanced ionizing output arising from metal poor gas and stars, young stellar populations, and large ionization parameters. The young ages implied by the emission lines and continuum SEDs indicate that the extreme line emitters in our sample are in the midst of a significant upturn in their star formation activity. The low stellar masses, blue UV colors, and large sSFRs of our sample are similar to those of typical z>6 galaxies. Given the strong attenuation of Ly-alpha in z>6 galaxies we suggest that CIII] is likely to provide our best probe of early star forming galaxies with ground-based spectrographs and one of the most efficient means of confirming z>10 galaxies with the James Webb Space Telescope.

preprint2013arXiv

Contamination of Broad-Band Photometry by Nebular Emission in High Redshift Galaxies: Investigations with Keck's MOSFIRE Near-Infrared Spectrograph

Earlier work has raised the potential importance of nebular emission in the derivation of the physical characteristics of high redshift Lyman break galaxies. Within certain redshift ranges, and especially at z ~ 6-7, such lines may be strong enough to reduce estimates of the stellar masses and ages of galaxies compared those derived assuming broad-band photometry represents stellar light alone. To test this hypothesis at the highest redshifts where such lines can be probed with ground-based facilities, we examine the near-infrared spectra of a representative sample of 20 3.0 < z < 3.8 Lyman break galaxies using the newly-commissioned MOSFIRE near-infrared spectrograph at the Keck I telescope. We use this data to derive the rest-frame equivalent widths (EW) of [O III] emission and show that these are comparable to estimates derived using the SED fitting technique introduced for sources of known redshift by Stark et al (2013). Although our current sample is modest, its [O III] EW distribution is consistent with that inferred for Hα based on SED fitting of Stark et al's larger sample of 3.8 < z < 5 galaxies. For a subset of survey galaxies, we use the combination of optical and near-infrared spectroscopy to quantify kinematics of outflows in z ~ 3.5 star-forming galaxies, and discuss the implications for reionization measurements. The trends we uncover underline the dangers of relying purely on broad-band photometry to estimate the physical properties of high redshift galaxies and emphasize the important role of diagnostic spectroscopy.

preprint2013arXiv

Deep LBT/LUCI Spectroscopy of a Lyman-alpha Emitter Candidate at z ~ 7.7

We present deep spectroscopic observations of a Lyman-alpha emitter (LAE) candidate at z ~ 7.7 using the infrared spectrograph LUCI on the 2 x 8.4m Large Binocular Telescope (LBT). The candidate is the brightest among the four z ~ 7.7 LAE candidates found in a narrow-band imaging survey by Krug et al. 2012. Our spectroscopic data include a total of 7.5 hours of integration with LBT/LUCI and are deep enough to significantly (3.2-4.9 sigma) detect the Lyman-alpha emission line of this candidate, based on its Lyman-alpha flux 1.2 x 10^{-17} erg s^{-1} cm^{-2} estimated from the narrow-band photometry. However, we do not find any convincing signal at the expected position of its Lyman-alpha emission line, suggesting that this source is not an LAE at z ~ 7.7. The non-detection in this work, together with the previous studies of z ~ 7.7 LAEs, puts a strong constraint on the bright-end Lyman-alpha luminosity function (LF) at z ~ 7.7. We find a rapid evolution of the Lyman-alpha LF from z ~ 6.5 to 7.7: the upper limit of the z ~ 7.7 LF is more than 5 times lower than the z ~ 6.5 LF at the bright end (f > 1.0 x 10^{-17} erg s^{-1} cm^{-2}, or L > 6.9 x 10^{42} erg s^{-1}). This is likely caused by an increasing neutral fraction in the IGM that substantially attenuates Lyman-alpha emission at z ~ 7.7.

preprint2013arXiv

Keck Spectroscopy of Gravitationally Lensed z=4 Galaxies: Improved Constraints on the Escape Fraction of Ionizing Photons

The fraction of ionizing photons that escape from young star-forming galaxies is one of the largest uncertainties in determining the role of galaxies in cosmic reionization. Yet traditional techniques for measuring this fraction are inapplicable at the redshifts of interest due to foreground screening by the Lyman alpha forest. In an earlier study, we demonstrated a reduction in the equivalent width of low-ionization absorption lines in composite spectra of Lyman break galaxies at z=4 compared to similar measures at z=3. This might imply a lower covering fraction of neutral gas and hence an increase with redshift in the escape fraction of ionizing photons. However, our spectral resolution was inadequate to differentiate between several alternative explanations, including changes with redshift in the outflow kinematics. Here we present higher quality spectra of 3 gravitationally lensed Lyman break galaxies at z=4 with a spectral resolution sufficient to break this degeneracy of interpretation. We present a method for deriving the covering fraction of low-ionization gas as a function of outflow velocity and compare the results with similar quality data taken for galaxies at lower redshift. We find a significant trend of lower covering fractions of low-ionization gas for galaxies with strong \Lya emission. In combination with the demographic trends of \Lya emission with redshift from our earlier work, our results provide new evidence for a reduction in the average H I covering fraction, and hence an increase in the escape fraction of ionizing radiation from Lyman break galaxies, with redshift.

preprint2013arXiv

New Constraints on Cosmic Reionization from the 2012 Hubble Ultra Deep Field Campaign

Understanding cosmic reionization requires the identification and characterization of early sources of hydrogen-ionizing photons. The 2012 Hubble Ultra Deep Field (UDF12) campaign has acquired the deepest infrared images with the Wide Field Camera 3 aboard Hubble Space Telescope and, for the first time, systematically explored the galaxy population deep into the era when cosmic microwave background (CMB) data indicates reionization was underway. The UDF12 campaign thus provides the best constraints to date on the abundance, luminosity distribution, and spectral properties of early star-forming galaxies. We synthesize the new UDF12 results with the most recent constraints from CMB observations to infer redshift-dependent ultraviolet (UV) luminosity densities, reionization histories, and electron scattering optical depth evolution consistent with the available data. Under reasonable assumptions about the escape fraction of hydrogen ionizing photons and the intergalactic medium clumping factor, we find that to fully reionize the universe by redshift z~6 the population of star-forming galaxies at redshifts z~7-9 likely must extend in luminosity below the UDF12 limits to absolute UV magnitudes of M_UV\sim -13 or fainter. Moreover, low levels of star formation extending to redshifts z~15-25, as suggested by the normal UV colors of z\simeq7-8 galaxies and the smooth decline in abundance with redshift observed by UDF12 to z\simeq10, are additionally likely required to reproduce the optical depth to electron scattering inferred from CMB observations.

preprint2013arXiv

The CASSOWARY spectroscopy survey: A new sample of gravitationally lensed galaxies in SDSS

Bright gravitationally lensed galaxies provide our most detailed view of galaxies at high redshift. Yet as a result of the small number of ultra-bright z~2 lensed systems with confirmed redshifts, most detailed spectroscopic studies have been limited in their scope. With the goal of increasing the number of bright lensed galaxies available for detailed follow-up, we have undertaken a spectroscopic campaign targeting wide separation (>3 arcsec) galaxy-galaxy lens candidates within the Sloan Digital Sky Survey (SDSS). Building on the earlier efforts of our CASSOWARY survey, we target a large sample of candidate galaxy-galaxy lens systems in SDSS using a well-established search algorithm which identifies blue arc-like structures situated around luminous red galaxies. In this paper, we present a new redshift catalog containing 25 lensed sources in SDSS confirmed through spectroscopic follow-up of candidate galaxy-galaxy lens systems. Included in this new sample are two of the brightest galaxies (r=19.6 and 19.7) galaxies known at z~2, a low metallicity (12 + log (O/H)~8.0) extreme nebular line emitting galaxy at z=1.43, and numerous systems for which detailed follow-up will be possible. The source redshifts span 0.9<z<2.5 (median redshift of 1.9), and their optical magnitudes are in the range 19.6<r<22.3. We present a brief source-by-source discussion of the spectroscopic properties extracted from our confirmatory spectra and discuss some initial science results. With more than 50 gravitationally lensed z>1 galaxies now confirmed within SDSS, it will soon be possible for the first time to develop generalized conclusions from detailed spectroscopic studies of the brightest lensed systems at high redshift.

preprint2013arXiv

Ultra-faint Ultraviolet Galaxies at z~2 Behind the Lensing Cluster Abell 1689: the Luminosity Function, Dust Extinction and Star Formation Rate Density

We identified the z~2 Lyman break galaxies using deep HST ultraviolet (F275W/F336W) imaging of Abell 1689. Because of the imaging depth and the large magnification provided by the cluster, we detect galaxies 100x fainter (-19.5< M_1500 <-13) than previous surveys at this redshift. We are able to calculate the intrinsic sensitivity of the observations as a function of source plane position, allowing determinations of effective volume as a function of luminosity. We fit the faint-end slope of the luminosity function to be alpha = -1.74 +/-0.08, consistent with the values obtained for 2.5 < z < 6. There is no turnover in the luminosity function down to MUV = -13. The trend of increasingly redder UV spectral slopes with luminosity at higher redshifts is observed in our sample, but with redder slopes at all luminosities and average reddening of < E(B - V) >= 0.15. We assume the stars in these galaxies are metal poor (0.2Z_sun) compared to their brighter counterparts (Z_sun), resulting in bluer assumed intrinsic UV slopes and larger derived dust extinction. The total UV luminosity density at z ~ 2 is 4.31x10^26 erg/s/Hz/Mpc^3, more than 70% of which is emitted by galaxies in the luminosity range of our sample. We determine the star formation rate density at z ~ 2 (assuming constant dust extinction correction of 4.2 over all luminosities and a Kroupa IMF) of 0.148 M/yr/Mpc^3, significantly higher than previous determinations because of the additional population of fainter galaxies and the larger dust correction factors.[abridged]

preprint2012arXiv

Evolution of the Sizes of Galaxies over 7<z<12 Revealed by the 2012 Hubble Ultra Deep Field Campaign

We analyze the redshift- and luminosity-dependent sizes of dropout galaxy candidates in the redshift range z~7-12 using deep images from the UDF12 campaign, data which offers two distinct advantages over that used in earlier work. Firstly, we utilize the increased S/N ratio offered by the UDF12 imaging to provide improved size measurements for known galaxies at z=6.5-8 in the HUDF. Specifically, we stack the new deep F140W image with the existing F125W data in order to provide improved measurements of the half-light radii of z-dropouts. Similarly we stack this image with the new deep UDF12 F160W image to obtain new size measurements for a sample of Y-dropouts. Secondly, because the UDF12 data have allowed the construction of the first robust galaxy sample in the HUDF at z>8, we have been able to extend the measurement of average galaxy size out to significantly higher redshifts. Restricting our size measurements to sources which are now detected at >15sigma, we confirm earlier indications that the average half-light radii of z~7-12 galaxies are extremely small, 0.3-0.4 kpc, comparable to the sizes of giant molecular associations in local star-forming galaxies. We also confirm that there is a clear trend of decreasing half-light radius with increasing redshift, and provide the first evidence that this trend continues beyond z~8. Modeling the evolution of the average half-light radius as a power-law (1+z)^s, we obtain a best-fit index of s=-1.28+/-0.13 over the redshift range z~4-12, mid-way between the physically expected evolution for baryons embedded in dark halos of constant mass (s=-1) and constant velocity (s=-1.5). A clear size-luminosity relation, such as that found at lower redshift, is also evident in both our z- and Y-dropout sample. This relation can be interpreted in terms of a constant surface density of star formation over a range in luminosity of 0.05-1.0L*_z=3.(abridged)

preprint2012arXiv

Keck Spectroscopy of 3<z<7 Faint Lyman Break Galaxies: The Importance of Nebular Emission in Understanding the Specific Star Formation Rate and Stellar Mass Density

The physical properties inferred from the SEDs of z>3 galaxies have been influential in shaping our understanding of early galaxy formation and the role galaxies may play in cosmic reionization. Of particular importance is the stellar mass density at early times which represents the integral of earlier star formation. An important puzzle arising from the measurements so far reported is that the specific star formation rates (sSFR) evolve far less rapidly than expected in most theoretical models. Yet the observations underpinning these results remain very uncertain, owing in part to the possible contamination of rest-optical broadband light from strong nebular emission lines. To quantify the contribution of nebular emission to broad-band fluxes, we investigate the SEDs of 92 spectroscopically-confirmed galaxies in the redshift range 3.8<z<5.0 chosen because the H-alpha line lies within the Spitzer/IRAC 3.6 um filter. We demonstrate that the 3.6 um flux is systematically in excess of that expected from stellar continuum, which we derive by fitting the SED with population synthesis models. No such excess is seen in a control sample at 3.1<z<3.6 in which there is no nebular contamination in the IRAC filters. From the distribution of our 3.6 um flux excesses, we derive an H-alpha equivalent width (EW) distribution. The mean rest-frame H-alpha EW we infer at 3.8<z<5.0 (270 A) indicates that nebular emission contributes at least 30% of the 3.6 um flux. Via our empirically-derived EW distribution we correct the available stellar mass densities and show that the sSFR evolves more rapidly at z>4 than previously thought, supporting up to a 5x increase between z~2 and 7. Such a trend is much closer to theoretical expectations. Given our findings, we discuss the prospects for verifying quantitatively the nebular emission line strengths prior to the launch of the James Webb Space Telescope.

preprint2012arXiv

The UV Luminosity Function of star-forming galaxies via dropout selection at redshifts z ~ 7 and 8 from the 2012 Ultra Deep Field campaign

We present a catalog of high redshift star-forming galaxies selected to lie within the redshift range z ~ 7-8 using the Ultra Deep Field 2012 (UDF12), the deepest near-infrared (near-IR) exposures yet taken with the Hubble Space Telescope. As a result of the increased near-infrared exposure time compared to previous HST imaging in this field, we probe 0.65 (0.25) mag fainter in absolute UV magnitude, at z ~ 7 (8), which increases confidence in a measurement of the faint end slope of the galaxy luminosity function. Through a 0.7 mag deeper limit in the key F105W filter that encompasses or lies just longward of the Lyman break, we also achieve a much-refined color-color selection that balances high redshift completeness and a low expected contamination fraction. We improve the number of drop-out selected UDF sources to 47 at z ~ 7 and 27 at z ~ 8. Incorporating brighter archival and ground-based samples, we measure the z ~ 7 UV luminosity function to an absolute magnitude limit of M_UV = -17 and find a faint end Schechter slope of α= -1.87+/- 0.18. Using a similar color-color selection at z ~ 8 that takes account of our newly-added imaging in the F140W filter, and incorporating archival data from the HIPPIES and BoRG campaigns, we provide a robust estimate of the faint end slope at z ~ 8, α= -1.94 +/- 0.23. We briefly discuss our results in the context of earlier work and that derived using the same UDF12 data but with an independent photometric redshift technique (McLure et al 2012).

preprint2011arXiv

Assessing the redshift evolution of massive black holes and their hosts

Motivated by recent observational results that focus on high redshift black holes, we explore the effect of scatter and observational biases on the ability to recover the intrinsic properties of the black hole population at high redshift. We find that scatter and selection biases can hide the intrinsic correlations between black holes and their hosts, with 'observable' subsamples of the whole population suggesting, on average, positive evolution even when the underlying population is characterized by no- or negative evolution. We create theoretical mass functions of black holes convolving the mass function of dark matter halos with standard relationships linking black holes with their hosts. Under these assumptions, we find that the local MBH - sigma correlation is unable to fit the z = 6 black hole mass function proposed by Willott et al. (2010), overestimating the number density of all but the most massive black holes. Positive evolution or including scatter in the MBH - sigma correlation makes the discrepancy worse, as it further increases the number density of observable black holes. We notice that if the MBH - sigma correlation at z = 6 is steeper than today, then the mass function becomes shallower. This helps reproducing the mass function of z = 6 black holes proposed by Willott et al. (2010). Alternatively, it is possible that very few halos (of order 1/1000) host an active massive black hole at z = 6, or that most AGN are obscured, hindering their detection in optical surveys. Current measurements of the high redshift black hole mass function might be underestimating the density of low mass black holes if the active fraction or luminosity are a function of host or black hole mass. Finally, we discuss physical scenarios that can possibly lead to a steeper MBH - sigma relation at high redshift.

preprint2011arXiv

Keck Spectroscopy of Faint 3 < z < 7 Lyman Break Galaxies: III. The Mean Ultraviolet Spectrum at z=4

We present and discuss the mean rest-frame ultraviolet spectrum for a sample of 81 Lyman Break Galaxies (LBGs) selected to be B-band dropouts with a mean redshift of z=3.9 and apparent magnitudes z_AB<26. Most of the individual spectra are drawn from our ongoing survey in the GOODS fields with the Keck DEIMOS spectrograph, and we have augmented our sample with published data taken with FORS2 on the VLT. In general we find similar trends in the spectral diagnostics to those found in the earlier, more extensive survey of LBGs at z=3 undertaken by Shapley et al (2003). Specifically, we find low-ionization absorption lines which trace the presence of neutral outflowing gas are weaker in galaxies with stronger Lyman-alpha emission, bluer UV spectral slopes, lower stellar masses, lower UV luminosities and smaller half-light radii. This is consistent with a physical picture whereby star formation drives outflows of neutral gas which scatters Lyman-alpha and gives rise to strong low-ionization absorption lines, while increasing the stellar mass, size, metallicity, and dust content of galaxies. Typical galaxies are thus expected to have stronger Lyman-alpha emission and weaker low-ionization absorption at earlier times (higher redshifts). Indeed, our mean spectrum at z=4 shows somewhat weaker low-ionization absorption lines than at z=3 and available data at higher redshift indicates a rapid decrease in low-ionization absorption strength with redshift. We argue that the reduced low-ionization absorption is likely caused by a decrease in the covering fraction and/or velocity range of outflowing neutral gas at earlier epochs. Our continuing survey will enable us to extend these diagnostics more reliably to higher redshift and determine the implications for the escape fraction of ionizing photons which governs the role of early galaxies in cosmic reionization. [Abridged]

preprint2011arXiv

LBG properties from z~3 to z~6

We analyse the spectral energy distribution (SED) of U, B, V and i-dropout samples from GOODS-MUSIC and we determine their physical properties, such as stellar age and mass, dust attenuation and star formation rate (SFR). Furthermore, we examine how the strength of Lyα emission can be constrained from broad-band SED fits instead of relying in spectroscopy. We use our SED fitting tool including the effects of nebular emission and we explore different star formation histories (SFHs). We find that SEDs are statistically better fitted with nebular emission and exponentially decreasing star formation. Considering this result, stellar mass and star formation rate (SFR) estimations modify the specific SFR (SFR/Mstar) - redshift relation, in comparison to previous studies. Finally, our inferred Lyα properties are in good agreement with the available spectroscopic observations.

preprint2011arXiv

On Lya emission in z~3-6 UV-selected galaxies

Determining Lya properties of distant galaxies is of great interest for various astrophysical studies. We examine how the strength of Lya emission can be constrained from broad-band SED fits instead of relying on spectroscopy. We use our SED fitting tool including the effects of nebular emission, considering in particular Lya emission as a free parameter, and we demonstrate our method with simulations of mock galaxies. Using this tool we analyse a large sample of U, B, V, and i dropout galaxies with multi-band photometry. We find significant trends of the fraction of galaxies with Lya emission increasing both with redshift z and towards fainter magnitude (at fixed z), and similar trends for the Lya equivalent width. Our inferred Lya properties are in good agreement with the available spectroscopic observations and other data. These results demonstrate that the strength of Lya emission in distant star-forming galaxies can be inferred quantitatively from broad-band SED fits, at least statistically for sufficiently large samples with a good photometric coverage.

preprint2011arXiv

The Emission Line Properties of Gravitationally-lensed 1.5 < z < 5 Galaxies

We present and analyse near-infrared spectroscopy for a sample of 28 gravitationally- lensed star-forming galaxies in the redshift range 1.5 < z < 5, observed mostly with the Keck II telescope. With typical magnifications of ~1.5-4 magnitudes, our survey provides a valuable census of star formation rates, gas-phase metallicities and dynamical masses for a representative sample of low luminosity galaxies seen at a formative period in cosmic history. We find less evolution in the mass-metallicity relation compared to earlier work that focused on more luminous systems with z - 2-3, especially in the low mass (- 10^9 Msol) where our sample is - 0.25 dex more metal-rich. We interpret this offset as a result of the lower star formation rates (typically a factor of -10 lower) for a given stellar mass in our sub-luminous systems. Taking this effect into account, we conclude our objects are consistent with a fundamental metallicity relation recently proposed from unlensed observations.

preprint2010arXiv

Early star-forming galaxies and the reionization of the Universe

Star forming galaxies represent a valuable tracer of cosmic history. Recent observational progress with Hubble Space Telescope has led to the discovery and study of the earliest-known galaxies corresponding to a period when the Universe was only ~800 million years old. Intense ultraviolet radiation from these early galaxies probably induced a major event in cosmic history: the reionization of intergalactic hydrogen. New techniques are being developed to understand the properties of these most distant galaxies and determine their influence on the evolution of the universe.

preprint2010arXiv

Keck Spectroscopy of Faint 3<z<7 Lyman Break Galaxies: - I. New constraints on cosmic reionisation from the luminosity and redshift-dependent fraction of Lyman-alpha emission

We present results from a new Keck spectroscopic survey of UV-faint LBGs in the redshift range 3<z<7. Combined with earlier Keck and published ESO VLT data, our sample contains more than 600 dropouts, offering new insight into the nature of sub-L* sources typical of those likely to dominate the cosmic reionisation process. Here we use this sample to characterise the fraction of strong Lya emitters within the continuum-selected dropouts. By quantifying how the "Lya fraction" varies with redshift, we seek to constrain changes in Lya transmission associated with reionisation. In order to distinguish the effects of reionisation from other factors which affect the Lya fraction (e.g. dust, ISM kinematics), we study the luminosity and redshift-dependence of the Lya fraction over 3<z<6, when the IGM is known to be ionised. These results reveal that low luminosity galaxies show strong Lya emission much more frequently than luminous systems, and that at fixed luminosity, the prevalence of strong Lya emission increases moderately with redshift over 3 < z < 6. Based on the correlation between blue UV slopes and strong Lya emitting galaxies in our dataset, we argue that the Lya fraction trends are governed by redshift and luminosity-dependent variations in the dust obscuration, with likely additional contributions from trends in the kinematics and covering fraction of neutral hydrogen. We find a tentative decrease in the Lya fraction at z~7 based on the limited IR spectra for candidate z~7 lensed LBGs, a result which, if confirmed with future surveys, would suggest an increase in the neutral fraction by this epoch. Given the supply of z and Y-drops now available from Hubble WFC3/IR surveys, we show it will soon be possible to significantly improve estimates of the Lya fraction using optical and near-IR spectrographs, thereby extending the study conducted in this paper to 7<z<8.

preprint2010arXiv

Keck Spectroscopy of Faint 3<z<7 Lyman Break Galaxies:- II. A High Fraction of Line Emitters at Redshift Six

As Lyman-alpha photons are scattered by neutral hydrogen, a change with redshift in the Lyman-alpha equivalent width distribution of distant galaxies offers a promising probe of the degree of ionization in the intergalactic medium and hence when cosmic reionization ended. This simple test is complicated by the fact that Lyman-alpha emission can also be affected by the evolving astrophysical details of the host galaxies. In the first paper in this series, we demonstrated both a luminosity and redshift dependent trend in the fraction of Lyman-alpha emitters seen within color-selected Lyman-break galaxies (LBGs) over the range 3<z<6; lower luminosity galaxies and those at higher redshift show an increased likelihood of strong emission. Here we present the results from much deeper 12.5 hour exposures with the Keck DEIMOS spectrograph focused primarily on LBGs at z~6 which enable us to confirm the redshift dependence of line emission more robustly and to higher redshift than was hitherto possible. We find 54+/-11% of faint z~6 LBGs show strong (W_0>25 A) emission, an increase of 1.6x from a similar sample observed at z~4. With a total sample of 74 z~6 LBGs, we determine the luminosity-dependent Lyman-alpha equivalent width distribution. Assuming continuity in these trends to the new population of z~7 sources located with the Hubble WFC3/IR camera, we predict that unless the neutral fraction rises in the intervening 200 Myr, the success rate for spectroscopic confirmation using Lyman-alpha emission should be high.

preprint2010arXiv

Star formation trends in high-redshift galaxy surveys: the elephant or the tail?

Star formation rate and accummulated stellar mass are two fundamental physical quantities that describe the evolutionary state of a forming galaxy. Two recent attempts to determine the relationship between these quantities, by interpreting a sample of star-forming galaxies at redshift of z~4, have led to opposite conclusions. We use a model galaxy population to investigate possible causes for this discrepancy and conclude that minor errors in the conversion from observables to physical quantities can lead to major misrepresentation when applied without awareness of sample selection. We also investigate, in a general way, the physical origin of the correlation between star formation rate and stellar mass within hierarchical galaxy formation theory.

preprint2009arXiv

A Study of Interstellar Gas and Stars in the Gravitationally Lensed Galaxy `The Cosmic Eye' from Rest-Frame Ultraviolet Spectroscopy

We report the results of a study of the rest-frame UV spectrum of the Cosmic Eye, a luminous Lyman break galaxy at z=3.07331 gravitationally lensed by a factor of 25. The spectrum, recorded with the ESI spectrograph on the Keck II telescope, is rich in absorption features from the gas and massive stars in this galaxy. The interstellar absorption lines are resolved into two components of approximately equal strength and each spanning several hundred km/s in velocity. One component has a net blueshift of -70 km/s relative to the stars and H II regions and presumably arises in a galaxy-scale outflow similar to those seen in most star-forming galaxies at z = 2-3. The other is more unusual in showing a mean redshift of +350 km/s relative to the systemic redshift; possible interpretations include a merging clump, or material ejected by a previous star formation episode and now falling back onto the galaxy, or more simply a chance alignment with a foreground galaxy. In the metal absorption lines, both components only partially cover the OB stars against which they are being viewed. We tentatively associate the redshifted component with the strong damped Lyman alpha line, indicative of a column density N(H I) = (3.0 +/- 0.8) x 10(21) atoms/cm2, and propose that it provides the dust `foreground screen' responsible for the low ratio of far-infrared to UV luminosities of the Cosmic Eye. Compared to other well-studied examples of strongly lensed galaxies, we find that the young stellar population of the Cosmic Eye is essentially indistinguishable from those of the Cosmic Horseshoe and MS 1512-cB58, while the interstellar spectra of all three galaxies are markedly different, attesting to the real complexity of the interplay between starbursts and ambient interstellar matter in young galaxies (abridged).

preprint2009arXiv

Detection of PAH and Far-Infrared Emission from the Cosmic Eye: Probing the Dust and Star Formation of Lyman Break Galaxies

We report the results of a Spitzer infrared study of the Cosmic Eye, a strongly lensed, L*_UV Lyman Break Galaxy (LBG) at z=3.074. We obtained Spitzer IRS spectroscopy as well as MIPS 24 and 70 micron photometry. The Eye is detected with high significance at both 24 and 70 microns and, when including a flux limit at 3.5 mm, we estimate an infrared luminosity of L_IR = 8.3 (+4.7-4.4) x10^11 L_sun assuming a magnification of 28+-3. This L_IR is eight times lower than that predicted from the rest-frame UV properties assuming a Calzetti reddening law. This has also been observed in other young LBGs, and indicates that the dust reddening law may be steeper in these galaxies. The mid-IR spectrum shows strong PAH emission at 6.2 and 7.7 microns, with equivalent widths near the maximum values observed in star-forming galaxies at any redshift. The L_PAH-to-L_IR ratio lies close to the relation measured in local starbursts. Therefore, L_PAH or L_MIR may be used to estimate L_IR and thus, star formation rate, of LBGs, whose fluxes at longer wavelengths are typically below current confusion limits. We also report the highest redshift detection of the 3.3 micron PAH emission feature. The PAH ratio, L_6.2/L_3.3=5.1+- 2.7, and the PAH-to-L_IR ratio, L_3.3/L_IR = 8.5 +- 4.7 x10^-4, are both in agreement with measurements in local starbursts and ULIRGs, suggesting that this line may serve as a good proxy for L_PAH or L_IR at z > 3 with the James Webb Space Telescope.

preprint2009arXiv

The Evolutionary History of Lyman Break Galaxies Between Redshift 4 and 6: Observing Successive Generations of Massive Galaxies in Formation

We present new measurements of the evolution in the Lyman break galaxy (LBG) population between z~4 and z~6. By utilizing the extensive multiwavelength datasets available in the GOODS fields, we identify 2443 B, 506 V, and 137 i'-band dropout galaxies likely to be at z~4, 5, and 6. With the goal of understanding the duration of typical star formation episodes in galaxies at z>4, we examine the distribution of stellar masses and ages as a function of cosmic time. We find that at a fixed rest-UV luminosity, the average stellar masses and ages of galaxies do not increase significantly between z~6 and 4. In order to maintain this near equilibrium in the average properties of high redshift LBGs, we argue that there must be a steady flux of young, newly-luminous objects at each successive redshift. When considered along with the short duty cycles inferred from clustering measurements, these results may suggest that galaxies are undergoing star formation episodes lasting only several hundred million years. In contrast to the unchanging relationship between the average stellar mass and rest-UV luminosity, we find that the number density of massive galaxies increases considerably with time over 4<z<6. Given this rapid increase of UV luminous massive galaxies, we explore the possibility that a significant fraction of massive (M*>1e11 Msun) z~2-3 distant red galaxies (DRGs) were in part assembled in an LBG phase at earlier times. Integrating the growth in the stellar mass function of actively forming LBGs over 4<z<6 down to z~2, we find that z>3 LBGs could have contributed significantly to the quiescent DRG population, indicating that the intense star-forming systems probed by current sub-millimeter observations are not the only route toward the assembly of DRGs at z~2.

preprint2007arXiv

A Keck Survey for Gravitationally-Lensed Lyman-alpha Emitters in the Redshift Range 8.5 < z < 10.4: New Constraints on the Contribution of Low Luminosity Sources to Cosmic Reionization

We discuss new observational constraints on the abundance of faint high redshift Lyman-alpha emitters secured from a deep Keck near-infrared spectroscopic survey which utilizes the strong magnification provided by lensing galaxy clusters at intermediate redshift (z=0.2-0.5). In each of 9 clusters we have undertaken a systematic `blind' search for line emission with NIRSPEC in the J-band within carefully-selected regions which offer very high magnifications (>10-50x) for background sources with redshifts z~10. The high magnification enables the detection of emission at unprecedented limits (10^41 - 10^42 erg/s), much fainter than those of conventional narrow band imaging and other spectroscopic searches. As the comoving volumes probed are small, our survey is designed to address the important question of whether intrinsically feeble star forming galaxies could provide the dominant ionizing flux at z~10. Our survey has yielded six promising (>5 sigma) candidate Lyman-alpha emitters which lie between z=8.7 and z=10.2. While the predicted reionization photon budget depends upon a large number of physical assumptions, our first glimpse at the z~10 Universe suggests that low luminosity star-forming galaxies contribute a significant proportion of the UV photons necessary for cosmic reionization. [Abridged]