Source author record

Hakim Atek

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

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

36 published item(s)

preprint2026arXiv

BEACON: JWST NIRCam Pure-parallel Imaging Survey. II. Physical Properties of $z=7-14$ Galaxies

We present photometric properties of 161 galaxy candidates at $z=7-14$ selected from the second data release (DR2) of BEACON, a JWST Cycle 2 pure-parallel NIRCam imaging program. Carefully selected from 36 independent pointings (corresponding to $\sim350$\,arcmin$^2$ sky coverage), and hence with reduced cosmic variance, our galaxy candidates provide an unbiased sample for investigating galaxy properties over a wide range of environments. We measure the physical properties, including UV continuum slope ($β_{\rm UV}$), stellar mass ($M_*$), star formation rate (SFR), and sizes. Our highest redshift galaxy candidate at $z=13.71\pm0.15$ has a remarkably bright UV luminosity of $M_{\rm UV}=-21.19\pm0.08$, making it the brightest galaxy at $z>12$ if spectroscopically confirmed. With an extremely blue UV slope, compact morphology, and high star formation rate surface density ($Σ_{\rm SFR}$), this candidate may have extremely low metallicity, high ionizing photon escape fraction, or contributions from an AGN. Among our multiple independent sightlines, we identify three fields of galaxy number overdensity with $>3σ$ significance. The properties of galaxies in various environments do not exhibit significant differences, implying either that accelerated galaxy evolution in overdense regions is not yet widespread at $z>7$, or that the current constraints are limited by sample size. Our simulations indicate that increasing the sample by an order of magnitude would allow such environmental trends to be robustly confirmed or ruled out, underscoring the importance of future pure-parallel observations.

preprint2026arXiv

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

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

preprint2023arXiv

Constraining the physical properties of the first lensed $z\sim10-16$ galaxy candidates with JWST

The first deep-field observations of the JWST have immediately yielded a surprisingly large number of very high redshift candidates, pushing the frontier of observability well beyond $z\gtrsim10$. We here present a detailed SED-fitting analysis of the 15 gravitationally lensed $z\sim10-16$ galaxy candidates detected behind the galaxy cluster SMACS J0723.3-7327 in Atek et al. (2022) using the BEAGLE tool. Our analysis makes use of dynamical considerations to place limits on the ages of these galaxies and of all three published SL models of the cluster to account for lensing systematics. We find these galaxies to have relatively low stellar masses $M_{\star}\sim10^7-10^8\,\mathrm{M}_{\odot}$ and young ages $t_{\mathrm{age}}\sim10-100$\,Myr. Due to their very blue UV-slopes, down to $β\sim-3$, all of the galaxies in our sample have extremely low dust attenuations $A_V\lesssim0.02$. Placing the measured parameters into relation, we find a very shallow $M_{\star}-M_{\mathrm{UV}}$-slope and high sSFRs above the main sequence of star-formation with no significant redshift-evolution in either relation. This is in agreement with the bright UV luminosities measured for these objects and indicates that we are naturally selecting galaxies that are currently undergoing a star-bursting episode at the time they are observed. Finally, we discuss the robustness of our high-redshift galaxy sample regarding low-redshift interlopers and conclude that low-redshift solutions can safely be ruled out for roughly half of the sample, including the highest-redshift galaxies at $z\sim12-16$. These objects represent compelling targets for spectroscopic follow-up observations with JWST and ALMA.

preprint2023arXiv

UNCOVER: Illuminating the Early Universe -- JWST/NIRSpec Confirmation of $z > 12$ Galaxies

Observations of high-redshift galaxies provide a critical direct test to the theories of early galaxy formation, yet to date, only three have been spectroscopically confirmed at $z>12$. Due to strong gravitational lensing over a wide area, the galaxy cluster field A2744 is ideal for searching for the earliest galaxies. Here we present JWST/NIRSpec observations of two galaxies: a robust detection at $z_{\rm spec} = 12.393^{+0.004}_{-0.001}$, and a plausible candidate at $z_{\rm spec} = 13.079^{+0.013}_{-0.001}$. The galaxies are discovered in JWST/NIRCam imaging and their distances are inferred with JWST/NIRSpec spectroscopy, all from the JWST Cycle 1 UNCOVER Treasury survey. Detailed stellar population modeling using JWST NIRCam and NIRSpec data corroborates the primeval characteristics of these galaxies: low mass ($\sim 10^8~{\rm M_\odot}$), young, rapidly-assembling, metal-poor, and star-forming. Interestingly, both galaxies are spatially resolved, having lensing-corrected rest-UV effective radii on the order of 300-400 pc, which are notably larger than other spectroscopically confirmed systems at similar redshifts. The observed dynamic range of $z \gtrsim 10$ sizes spans over 1 order of magnitude, implying a significant scatter in the size-mass relation at early times. Deep into the epoch of reionization, these discoveries elucidate the emergence of the first galaxies.

preprint2022arXiv

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

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

preprint2022arXiv

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

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

preprint2022arXiv

The Star Formation Burstiness and Ionizing Efficiency of Low-mass Galaxies

We investigate the burstiness of star formation and the ionizing efficiency of a large sample of galaxies at $0.7 < z < 1.5$ using HST grism spectroscopy and deep ultraviolet (UV) imaging in the GOODS-N and GOODS-S fields. The star formation history (SFH) in these strong emission line low-mass galaxies indicates an elevated star formation rate (SFR) based on the H$α$ emission line at a given stellar mass when compared to the standard main sequence. Moreover, when comparing the H$α$ and UV SFR indicators, we find that an excess in SFR(H$α$) compared to SFR(UV) is preferentially observed in lower-mass galaxies below $10^{9}$ M$\odot$, which are also the highest-EW galaxies. These findings suggest that the burstiness parameters of these strong emission line galaxies may differ from those inferred from hydrodynamical simulations and previous observations. For instance, a larger burstiness duty cycle would explain the observed SFR(H$α$) excess. We also estimate the ionizing photon production efficiency $ξ_{ion}$, finding a median value of Log($ξ_{ion}$/erg$^{-1}$ Hz)$=24.80 \pm 0.26$ when adopting a Galactic dust correction for H$α$ and an SMC one for the stellar component. We observe an increase of $ξ_{ion}$ with redshift, further confirming similar results at higher redshifts. We also find that $ξ_{ion}$ is strongly correlated with EW(H$α$), which provides an approach for deriving $ξ_{ion}$ in early galaxies. Lower-mass, lower-luminosity galaxies have a higher $ξ_{ion}$. Overall, these results provide further support for faint galaxies playing a major role in the reionization of the Universe.

preprint2021arXiv

How robustly can we constrain the low-mass end of the $z\sim6-7$ stellar mass function? -- The limits of lensing models and stellar population assumptions in the Hubble Frontier Fields

We present new measurements of the very low-mass end of the galaxy stellar mass function (GSMF) at $z\sim6-7$ computed from a rest-frame ultraviolet selected sample of dropout galaxies. These galaxies lie behind the six Hubble Frontier Fields clusters and are all gravitationally magnified. Using deep Spitzer/IRAC and Hubble Space Telescope imaging, we derive stellar masses by fitting galaxy spectral energy distributions and explore the impact of different model assumptions and parameter degeneracies on the resulting GSMF. Our sample probes stellar masses down to $M_{\star}>10^{6}\,\text{M}_{\odot}$ and we find the $z\sim6-7$ GSMF to be best parametrized by a modified Schechter function which allows for a turnover at very low masses. Using a Monte-Carlo Markov Chain analysis of the GSMF, including accurate treatment of lensing uncertainties, we obtain a relatively steep low-mass end slope $α\simeq-1.96_{-0.08}^{+0.09}$ and a turnover at $\log(M_T/\text{M}_{\odot})\simeq7.10_{-0.56}^{+0.17}$ with a curvature of $β\simeq1.00_{-0.73}^{+0.87}$ for our minimum assumption model with constant star-formation history (SFH) and low dust attenuation, $A_V\leq0.2$. We find that the $z\sim6-7$ GSMF, in particular its very low-mass end, is significantly affected by the assumed functional form of the star formation history and the degeneracy between stellar mass and dust attenuation. For example, the low-mass end slope ranges from $α\simeq-1.82_{-0.07}^{+0.08}$ for an exponentially rising SFH to $α\simeq-2.34_{-0.10}^{+0.11}$ when allowing $A_V$ of up to 3.25. Future observations at longer wavelengths and higher angular resolution with the James Webb Space Telescope are required to break these degeneracies and to robustly constrain the stellar mass of galaxies on the extreme low-mass end of the GSMF.

preprint2021arXiv

Spectroscopically Identified Emission Line Galaxy Pairs in the WISP survey

We identify a sample of spectroscopically measured emission line galaxy (ELG) pairs up to z=1.6 from the WFC3 Infrared Spectroscopic Parallels (WISP) survey. WISP obtained slitless, near-infrared grism spectroscopy along with direct imaging in the J and H bands by observing in the pure-parallel mode with the Wide Field Camera Three (WFC3) on the Hubble Space Telescope (HST). From our search of 419 WISP fields covering an area of ~0.5 deg$^{2}$, we find 413 ELG pair systems, mostly Halpha emitters. We then derive reliable star formation rates (SFRs) based on the attenuation-corrected Halpha fluxes. Compared to isolated galaxies, we find an average SFR enhancement of 40%-65%, which is stronger for major pairs and pairs with smaller velocity separations (Delta_v < 300 km/s). Based on the stacked spectra from various subsamples, we study the trends of emission line ratios in pairs, and find a general consistency with enhanced lower-ionization lines. We study the pair fraction among ELGs, and find a marginally significant increase with redshift $f \propto (1+z)^α$, where the power-law index α=0.58$\pm$0.17 from $z\sim$0.2 to $z\sim$1.6. The fraction of Active galactic Nuclei (AGNs), is found to be the same in the ELG pairs as compared to isolated ELGs.

preprint2020arXiv

Reionization with galaxies and active galactic nuclei

In this work we investigate the properties of the sources that reionized the intergalactic medium (IGM) in the high-redshift Universe. Using a semi-analytical model aimed at reproducing galaxies and black holes in the first 1.5 Gyr of the Universe, we revisit the relative role of star formation and black hole accretion in producing ionizing photons that can escape into the IGM. Both star formation and black hole accretion are regulated by supernova feedback, resulting in black hole accretion being stunted in low-mass halos. We explore a wide range of combinations for the escape fraction of ionizing photons (redshift-dependent, constant and scaling with stellar mass) from both star formation ($\langle f_{\rm esc}^{\rm sf} \rangle$) and AGN ($f_{\rm esc}^{\rm bh}$) to find: (i) the ionizing budget is dominated by stellar radiation from low stellar mass ($M_*<10^9 {\rm M_\odot}$ ) galaxies at $z>6$ with the AGN contribution (driven by $M_{bh}>10^6 {\rm M_\odot}$ black holes in $M_* > 10^9 {\rm M_\odot}$ galaxies) dominating at lower redshifts; (ii) AGN only contribute $10-25\%$ to the cumulative ionizing emissivity by $z=4$ for the models that match the observed reionization constraints; (iii) if the stellar mass dependence of $\langle f_{\rm esc}^{\rm sf} \rangle$ is shallower than $f_{\rm esc}^{\rm bh}$, at $z<7$ a transition stellar mass exists above which AGN dominate the escaping ionizing photon production rate; (iv) the transition stellar mass decreases with decreasing redshift. While AGN dominate the escaping emissivity above the knee of the stellar mass function at $z \sim 6.8$, they take-over at stellar masses that are a tenth of the knee mass by $z=4$.

preprint2020arXiv

The BUFFALO HST Survey

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

preprint2015arXiv

Are Ultra-faint Galaxies at z=6-8 Responsible for Cosmic Reionization ? Combined Constraints from the Hubble Frontier Fields Clusters and Parallels

We use deep Hubble Space Telescope imaging of the Frontier Fields to accurately measure the galaxy rest-frame ultraviolet luminosity function (UV LF) in the redshift range $z \sim 6-8$. We combine observations in three lensing clusters A2744, MACS0416, MACS0717 and their associated parallels fields to select high-redshift dropout candidates. We use the latest lensing models to estimate the flux magnification and the effective survey volume in combination with completeness simulations performed in the source plane. We report the detection of 227 galaxy candidates at $z=6-7$ and 25 candidates at $z \sim 8$. While the total survey area is about 4 arcmin$^{2}$ in each parallel field, it drops to about 0.6 to 1 arcmin$^{2}$ in the cluster core fields because of the strong lensing. We compute the UV luminosity function at $z \sim 7$ using the combined galaxy sample and perform Monte Carlo simulations to determine the best fit Schechter parameters. We are able to reliably constrain the LF down to an absolute magnitude of $M_{UV}=-15.25$, which corresponds to 0.005$L^{\star}$. More importantly, we find that the faint-end slope remains steep down to this magnitude limit with $α=-2.04_{-0.17}^{+0.13}$. Our results confirm the most recent results in deep blank fields but extend the LF measurements more than two magnitudes deeper. The UV LF at $z \sim 8$ is not very well constrained below $M_{UV}=-18$ due to the small number statistics and incompleteness uncertainties. To assess the contribution of galaxies to cosmic reionization we derive the UV luminosity density at $z\sim7$ by integrating the UV LF down to an observationally constrained limit of $M_{UV} = -15$. We show that our determination of Log($ρ_{UV}$)=$26.2\pm0.13$ (erg s$^{-1}$ Hz$^{-1}$ Mpc$^{-3}$) can be sufficient to maintain the IGM ionized.

preprint2015arXiv

Hubble Frontier Fields : A High-Precision Strong-Lensing Mass Model of the Massive Galaxy Cluster Abell 2744 using 180 Multiple Images

We present a high-precision mass model of galaxy cluster Abell 2744, based on a strong-gravitational-lensing analysis of the \emph{Hubble Space Telescope Frontier Fields} (HFF) imaging data, which now include both \emph{Advanced Camera for Surveys} and \emph{Wide-Field Camera 3} observations to the final depth. Taking advantage of the unprecedented depth of the visible and near-infrared data, we identify 34 new multiply imaged galaxies, bringing the total to 61, comprising 181 individual lensed images. In the process, we correct previous erroneous identifications and positions of multiple systems in the northern part of the cluster core. With the \textsc{Lenstool} software and the new sets of multiple images, we model the cluster using two cluster-scale dark matter halos plus galaxy-scale halos for the cluster members. Our best-fit model predicts image positions with an \emph{RMS} error of 0.69$\arcsec$, which constitutes an improvement by almost a factor of two over previous parametric models of this cluster. We measure the total projected mass inside a 200~kpc aperture as ($2.162\pm 0.005$)$\times 10^{14}M_{\odot}$, thus reaching 1\% level precision for the second time, following the recent HFF measurement of MACSJ0416.1-2403. Importantly, the higher quality of the mass model translates into an overall improvement by a factor of 4 of the derived magnification factor. % for the high-redshift lensed background galaxies. Together with our previous HFF gravitational lensing analysis, this work demonstrates that the HFF data enables high-precision mass measurements for massive galaxy clusters and the derivation of robust magnification maps to probe the early Universe.

preprint2015arXiv

The Lyman Alpha Reference Sample: V. The impact of neutral ISM kinematics and geometry on Lyman Alpha escape

We present high-resolution far-UV spectroscopy of the 14 galaxies of the Lyman Alpha Reference Sample; a sample of strongly star-forming galaxies at low redshifts ($0.028 < z < 0.18$). We compare the derived properties to global properties derived from multi band imaging and 21 cm HI interferometry and single dish observations, as well as archival optical SDSS spectra. Besides the Lyman $α$ line, the spectra contain a number of metal absorption features allowing us to probe the kinematics of the neutral ISM and evaluate the optical depth and and covering fraction of the neutral medium as a function of line-of-sight velocity. Furthermore, we show how this, in combination with precise determination of systemic velocity and good Ly$α$ spectra, can be used to distinguish a model in which separate clumps together fully cover the background source, from the "picket fence" model named by Heckman et al. (2011). We find that no one single effect dominates in governing Ly$α$ radiative transfer and escape. Ly$α$ escape in our sample coincides with a maximum velocity-binned covering fraction of $\lesssim 0.9$ and bulk outflow velocities of $\gtrsim 50$ km s$^{-1}$, although a number of galaxies show these characteristics and yet little or no Ly$α$ escape. We find that Ly$α$ peak velocities, where available, are not consistent with a strong backscattered component, but rather with a simpler model of an intrinsic emission line overlaid by a blueshifted absorption profile from the outflowing wind. Finally, we find a strong anticorrelation between H$α$ equivalent width and maximum velocity-binned covering factor, and propose a heuristic explanatory model.

preprint2015arXiv

The Lyman alpha reference sample. VII. Spatially resolved H$α$ kinematics

We present integral field spectroscopic observations with the Potsdam Multi Aperture Spectrophotometer of all 14 galaxies in the $z\sim 0.1$ Lyman Alpha Reference Sample (LARS). We produce 2D line of sight velocity maps and velocity dispersion maps from the Balmer $α$ (H$α$) emission in our data cubes. These maps trace the spectral and spatial properties of the LARS galaxies' intrinsic Ly$α$ radiation field. We show our kinematic maps spatially registered onto the Hubble Space Telescope H$α$ and Lyman $α$ (Ly$α$) images. Only for individual galaxies a causal connection between spatially resolved H$α$ kinematics and Ly$α$ photometry can be conjectured. However, no general trend can be established for the whole sample. Furthermore, we compute non-parametric global kinematical statistics -- intrinsic velocity dispersion $σ_0$, shearing velocity $v_\mathrm{shear}$, and the $v_\mathrm{shear}/σ_0$ ratio -- from our kinematic maps. In general LARS galaxies are characterised by high intrinsic velocity dispersions (54\,km\,s$^{-1}$ median) and low shearing velocities (65\,km\,s$^{-1}$ median). $v_\mathrm{shear}/σ_0$ values range from 0.5 to 3.2 with an average of 1.5. Noteworthy, five galaxies of the sample are dispersion dominated systems with $v_\mathrm{shear}/σ_0 <1$ and are thus kinematically similar to turbulent star forming galaxies seen at high redshift. When linking our kinematical statistics to the global LARS Ly$α$ properties, we find that dispersion dominated systems show higher Ly$α$ equivalent widths and higher Ly$α$ escape fractions than systems with $v_\mathrm{shear}/σ_0 > 1$. Our result indicates that turbulence in actively star-forming systems is causally connected to interstellar medium conditions that favour an escape of Ly$α$ radiation.

preprint2015arXiv

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

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

preprint2014arXiv

Hubble Frontier Fields : A High Precision Strong Lensing Analysis of Galaxy Cluster MACSJ0416.1-2403 using ~200 Multiple Images

We present a high-precision mass model of the galaxy cluster MACSJ0416.1-2403, based on a strong-gravitational-lensing analysis of the recently acquired Hubble Space Telescope Frontier Fields (HFF) imaging data. Taking advantage of the unprecedented depth provided by HST/ACS observations in three passbands, we identify 51 new multiply imaged galaxies, quadrupling the previous census and bringing the grand total to 68, comprising 194 individual lensed images. Having selected a subset of the 57 most securely identified multiply imaged galaxies, we use the Lenstool software package to constrain a lens model comprised of two cluster-scale dark-matter halos and 98 galaxy-scale halos. Our best-fit model predicts image positions with an $RMS$ error of 0.68'', which constitutes an improvement of almost a factor of two over previous, pre-HFF models of this cluster. We find the total projected mass inside a 200~kpc aperture to be $(1.60\pm0.01)\times 10^{14}\ M_\odot$, a measurement that offers a three-fold improvement in precision, reaching the percent level for the first time in any cluster. Finally, we quantify the increase in precision of the derived gravitational magnification of high-redshift galaxies and find an improvement by a factor of $\sim$2.5 in the statistical uncertainty. Our findings impressively confirm that HFF imaging has indeed opened the domain of high-precision mass measurements for massive clusters of galaxies.

preprint2014arXiv

Hubble Frontier Fields: The Geometry and Dynamics of the Massive Galaxy Cluster Merger MACSJ0416.1-2403

We use a joint optical/X-ray analysis to constrain the geometry and history of the ongoing merging event in the massive galaxy cluster MACSJ0416.1-2403 (z=0.397). Our investigation of cluster substructure rests primarily on a combined strong- and weak-lensing mass reconstruction based on the deep, high-resolution images obtained for the Hubble Frontier Fields initiative. To reveal the system's dynamics, we complement this lensing analysis with a study of the intra-cluster gas using shallow Chandra data, and a three-dimensional model of the distribution and motions of cluster galaxies derived from over 100 spectroscopic redshifts. The multi-scale grid model obtained from our combined lensing analysis extends the high-precision strong-lensing mass reconstruction recently performed to cluster-centric distances of almost 1 Mpc. Our analysis detects the two well known mass concentrations in the cluster core. A pronounced offset between collisional and collisionless matter is only observed for the SW cluster component, while excellent alignment is found for the NE cluster. Both the lensing analysis and the distribution of cluster light strongly suggest the presence of a third massive structure, almost 2 arcmin SW of the cluster centre. Since no X-ray emission is detected in this region, we conclude that this structure is non-virialised and speculate that it might be part of a large-scale filament almost aligned with our line of sight. Combining all evidence from the distribution of dark and luminous matter, we propose two alternative scenarios for the trajectories of the components of MACSJ0416.1-2403. Upcoming deep X-ray observations that allow the detection of shock fronts, cold cores, and sloshing gas (all key diagnostics for studies of cluster collisions) will allow us to test, and distinguish between these two scenarios.

preprint2014arXiv

Hubble Space Telescope Grism Spectroscopy of Extreme Starbursts Across Cosmic Time: The Role of Dwarf Galaxies in the Star Formation History of the Universe

Near infrared slitless spectroscopy with the Wide Field Camera 3, onboard the Hubble Space Telescope, offers a unique opportunity to study low-mass galaxy populations at high-redshift ($z\sim$1-2). While most high$-z$ surveys are biased towards massive galaxies, we are able to select sources via their emission lines that have very-faint continua. We investigate the star formation rate (SFR)-stellar mass ($M_{\star}$) relation for about 1000 emission-line galaxies identified over a wide redshift range of $0.3 \lesssim z \lesssim 2.3$. We use the H$_α$ emission as an accurate SFR indicator and correct the broadband photometry for the strong nebular contribution to derive accurate stellar masses down to $M_{\star} \sim 10^{7} M_{\odot}$. We focus here on a subsample of galaxies that show extremely strong emission lines (EELGs) with rest-frame equivalent widths ranging from 200 to 1500 Å. This population consists of outliers to the normal SFR-$M_{\star}$ sequence with much higher specific SFRs ($> 10$ Gyr$^{-1}$). While on-sequence galaxies follow a continuous star formation process, EELGs are thought to be caught during an extreme burst of star formation that can double their stellar mass in less than $100$ Myr. The contribution of starbursts to the total star formation density appears to be larger than what has been reported for more massive galaxies in previous studies. In the complete mass range $8.2 <$ log($M_{\star}/M_{\odot}$) $< 10$ and a SFR lower completeness limit of about 2 $M_{\odot}$ yr$^{-1}$ (10 $M_{\odot}$ yr$^{-1}$) at $z\sim1$ ($z \sim 2$), we find that starbursts having EW$_{rest}$(H$_α$)$>$ 300, 200, and 100 A contribute up to $\sim13$, 18, and 34 %, respectively, to the total SFR of emission-line selected sample at $z\sim1-2$. The comparison with samples of massive galaxies shows an increase in the contribution of starbursts towards lower masses.

preprint2014arXiv

New Constraints on the Faint-end of the UV Luminosity Function at z~7-8 using the Gravitational Lensing of the Hubble Frontier Fields Cluster A2744

Exploiting the power of gravitational lensing, the Hubble Frontier Fields (HFF) program aims at observing six massive galaxy clusters to explore the distant Universe far beyond the depth limits of blank field surveys. Using the complete Hubble Space Telescope observations of the first HFF cluster Abell 2744, we report the detection of 50 galaxy candidates at $z \sim 7$ and eight candidates at $z \sim 8$ in a total survey area of 0.96 arcmin$^{2}$ in the source plane. Three of these galaxies are multiply-imaged by the lensing cluster. Using an updated model of the mass distribution in the cluster we were able to calculate the magnification factor and the effective survey volume for each galaxy in order to compute the ultraviolet galaxy luminosity function at both redshifts 7 and 8. Our new measurements extend the $z \sim 7$ UV LF down to an absolute magnitude of $M_{UV} \sim -15.5$. We find a characteristic magnitude of $M^{\star}_{UV}=-20.63^{+0.69}_{-0.56}$ mag and a faint-end slope $α= -1.88^{+0.17}_{-0.20}$ close to previous determinations in blank fields. We show here for the first time that this slope remains steep down to very faint luminosities of 0.01$L^{\star}$. Although prone to large uncertainties, our results at $z \sim 8$ also seem to confirm a steep faint-end slope below 0.1$L^{\star}$. The HFF program is therefore providing an extremely efficient way to study the faintest galaxy populations at $z > 7$ that would otherwise be inaccessible with current instrumentation. The full sample of six galaxy clusters will provide yet better constraints on the build-up of galaxies at early epochs and their contribution to cosmic reionization.

preprint2014arXiv

Physical Properties of Emission-Line Galaxies at z ~ 2 from Near-Infrared Spectroscopy with Magellan FIRE

We present results from near-infrared spectroscopy of 26 emission-line galaxies at z ~ 2 obtained with the FIRE spectrometer on the Magellan Baade telescope. The sample was selected from the WISP survey, which uses the near-infrared grism of the Hubble Space Telescope Wide Field Camera 3 to detect emission-line galaxies over 0.3 < z < 2.3. Our FIRE follow-up spectroscopy (R~5000) over 1.0-2.5 micron permits detailed measurements of physical properties of the z~2 emission-line galaxies. Dust-corrected star formation rates for the sample range from ~5-100 M_sun yr-1. We derive a median metallicity for the sample of ~0.45 Z_sun, and the estimated stellar masses range from ~10^8.5 - 10^9.5 M_sun. The average ionization parameters measured for the sample are typically much higher than what is found for local star-forming galaxies. We derive composite spectra from the FIRE sample, from which we infer typical nebular electron densities of ~100-400 cm^-3. Based on the location of the galaxies and composite spectra on BPT diagrams, we do not find evidence for significant AGN activity in the sample. Most of the galaxies as well as the composites are offset in the BPT diagram toward higher [O III]/H-beta at a given [N II]/H-alpha, in agreement with other observations of z > 1 star-forming galaxies, but composite spectra derived from the sample do not show an appreciable offset from the local star-forming sequence on the [O III]/H-beta versus [S II]/H-alpha diagram. We infer a high nitrogen-to-oxygen abundance ratio from the composite spectrum, which may contribute to the offset of the high-redshift galaxies from the local star-forming sequence in the [O III]/H-beta versus [N II]/H-alpha diagram. We speculate that the elevated nitrogen abundance could result from substantial numbers of Wolf-Rayet stars in starbursting galaxies at z~2. (Abridged)

preprint2014arXiv

Probing the z>6 Universe with the first Hubble Frontier Fields cluster Abell 2744

The Hubble Frontier Fields (HFF) program combines the capabilities of the Hubble Space Telescope (HST) with the gravitational lensing of massive galaxy clusters to probe the distant Universe to an unprecedented depth. Here, we present the results of the first combined HST and Spitzer observations of the cluster Abell 2744. We combine the full near-infrared data with ancillary optical images to search for gravitationally lensed high-redshift (z > 6) galaxies. We report the detection of 15 I814-dropout candidates at z ~ 6-7 and one Y105-dropout at z ~ 8 in a total survey area of 1.43 arcmin^2 in the source plane. The predictions of our lens model allow us to also identify five multiply-imaged systems lying at redshifts between z ~ 6 and z ~ 8. Thanks to constraints from the mass distribution in the cluster, we were able to estimate the effective survey volume corrected for completeness and magnification effects. This was in turn used to estimate the rest-frame ultraviolet luminosity function (LF) at z ~ 6-8. Our LF results are generally in agreement with the most recent blank field estimates, confirming the feasibility of surveys through lensing clusters. Although based on a shallower observations than what will be achieved in the final dataset including the full ACS observations, the LF presented here extends down to Muv ~ -18.5 at z ~ 7 with one identified object at Muv ~ -15 thanks to the highly-magnified survey areas. This early study forecasts the power of using massive galaxy clusters as cosmic telescopes and its complementarity to blank fields.

preprint2014arXiv

The Lyman-alpha Reference Sample: I. Survey outline and first results for Markarian 259

The Lyman-alpha reference sample (LARS) is a program with the Hubble Space Telescope (HST) that provides a sample of local universe laboratory galaxies in which to study the astrophysics of the visibility and strength of the Lyman-alpha (Lya) line of hydrogen. This article presents an overview of the survey, its selection function and HST imaging observations. The sample was selected from the GALEX+SDSS catalogue at z=0.028-0.19, in order to allow Lya to be captured with combinations of long pass filters in the Solar Blind Channel (SBC) of HST/ACS. In addition, LARS utilises Halpha and Hbeta narrow, and U, B, i broad-band imaging with ACS and WFC3. In order to study galaxies in which large numbers of Lya photons are produced we demanded an Halpha equivalent width > 100Å. The sample of 14 galaxies covers far UV (FUV) luminosities that overlaps with those of high-z Lya emitters and Lyman Break Galaxies (LBGs). We present the reduction steps used to obtain the Lya images, including our LARS eXtraction software (LaXs) which utilises pixel-by-pixel spectral synthesis fitting to subtract the continuum at Lya. We demonstrate that the use of SBC long pass filters dramatically increase the signal to noise compared to the nominal Lya filter in SBC. To exemplify, we present results for Mrk 259. This irregular galaxy shows extended but strongly asymmetric Lya emission. Spectroscopy from HST/COS show a moderate outflow in the neutral ISM and an asymmetric Lya profile. Radiative transfer modeling is able to reproduce the Lya line profile, confirming the presence of an outflow. The integrated photometry provides a Lya luminosity of 1.3E42 erg/s an equivalent width W(Lya) =45Å and a FUV absolute magnitude -19.2. Mrk 259 would hence be detectable in high-z Lya and LBG surveys. The total Lya escape fraction is 12%, similar to that observed at z>4. (shortened)

preprint2013arXiv

Low Masses and High Redshifts: The Evolution of the Mass-Metallicity Relation

We present the first robust measurement of the high redshift mass-metallicity (MZ) relation at 10^{8}< M/M_{\sun} < 10^{10}, obtained by stacking spectra of 83 emission-line galaxies with secure redshifts between 1.3 < z < 2.3. For these redshifts, infrared grism spectroscopy with the Hubble Space Telescope Wide Field Camera 3 is sensitive to the R23 metallicity diagnostic: ([OII]3726,3729 + [OIII] 4959,5007)/Hβ. Using spectra stacked in four mass quartiles, we find a MZ relation that declines significantly with decreasing mass, extending from 12+log(O/H) = 8.8 at M=10^{9.8} M_{\sun} to 12+log(O/H)= 8.2 at M=10^{8.2} M_{\sun}. After correcting for systematic offsets between metallicity indicators, we compare our MZ relation to measurements from the stacked spectra of galaxies with M>10^{9.5} M_{\sun} and z~2.3. Within the statistical uncertainties, our MZ relation agrees with the z~2.3 result, particularly since our somewhat higher metallicities (by around 0.1 dex) are qualitatively consistent with the lower mean redshift z=1.76 of our sample. For the masses probed by our data, the MZ relation shows a steep slope which is suggestive of feedback from energy-driven winds, and a cosmological downsizing evolution where high mass galaxies reach the local MZ relation at earlier times. In addition, we show that our sample falls on an extrapolation of the star-forming main sequence (the SFR-M_{*} relation) at this redshift. This result indicates that grism emission-line selected samples do not have preferentially high SFRs. Finally, we report no evidence for evolution of the mass-metallicity-SFR plane; our stack-averaged measurements show excellent agreement with the local relation.

preprint2013arXiv

On the influence of physical galaxy properties on Lyman-alpha escape in star-forming galaxies

[abridged] Among the different observational techniques used to select high-redshift galaxies, the hydrogen recombination line Lyman-alpha (Lya) is of particular interest as it gives access to the measurement of cosmological quantities such as the star formation rate of distant galaxy populations. However, the interpretation of this line and the calibration of such observables is still subject to serious uncertainties. Therefore, it important to understand under what conditions the Lya line can be used as a reliable star formation diagnostic tool. We use a sample of 24 Lya emitters at z ~ 0.3 with an optical spectroscopic follow-up to calculate the Lya escape fraction and its dependency upon different physical properties. We also examine the reliability of Lya as a star formation rate indicator. We combine these observations with a compilation of Lya emitters selected at z = 0 - 0.3 to assemble a larger sample. The Lya escape fraction depends clearly on the dust extinction following the relation fesc(Lya) = C(Lya) x 10^(-0.4 E(B-V) k(Lya)), but with a shallower slope than previously reported, with k(Lya) ~ 6.67 and C(Lya) = 0.22. However, the correlation does not follow the expected curve for a simple dust attenuation. We explore the various mechanisms than lead to fesc(Lya) values above the continuum extinction curve, i.e. to an enhancement of the Lya output. We also observe that the strength of Lya and the escape fraction appear unrelated to the galaxy metallicity. Regarding the reliability of Lya as a star formation rate (SFR) indicator, we show that the deviation of SFR(Lya) from the true SFR (as traced by the UV continuum) is a function of the observed SFR(UV), which can be seen as the decrease of fesc(Lya) with increasing UV luminosity. Moreover, we observe a redshift-dependence of this relationship revealing the underlying evolution of fesc(Lya) with redshift.

preprint2013arXiv

Predicting Future Space Near-IR Grism Surveys using the WFC3 Infrared Spectroscopic Parallels Survey

We present near-infrared emission line counts and luminosity functions from the HST WFC3 Infrared Spectroscopic Parallels (WISP) program for 29 fields (0.037 deg^2) observed using both the G102 and G141 grisms. Altogether we identify 1048 emission line galaxies with observed equivalent widths greater than 40 Angstroms, 467 of which have multiple detected emission lines. The WISP survey is sensitive to fainter flux levels (3-5x10^-17 ergs/s/cm^2) than the future space near-infrared grism missions aimed at baryonic acoustic oscillation cosmology (1-4x10^-16 ergs/s/cm^2), allowing us to probe the fainter emission line galaxies that the shallower future surveys may miss. Cumulative number counts of 0.7<z<1.5 galaxies reach 10,000 deg^-2 above an H-alpha flux of 2x10^-16 ergs/s/cm^2. H-alpha-emitting galaxies with comparable [OIII] flux are roughly 5 times less common than galaxies with just H-alpha emission at those flux levels. Galaxies with low H-alpha/[OIII] ratios are very rare at the brighter fluxes that future near-infrared grism surveys will probe; our survey finds no galaxies with H-alpha/[OIII] < 0.95 that have H-alpha flux greater than 3x10^-16 ergs/s/cm^2. Our H-alpha luminosity function contains a comparable number density of faint line emitters to that found by the NICMOS near-infrared grism surveys, but significantly fewer (factors of 3-4 less) high luminosity emitters. We also find that our high redshift (z=0.9-1.5) counts are in agreement with the high redshift (z=1.47) narrow band H-alpha survey of HiZELS (Sobral et al. 2013), while our lower redshift luminosity function (z=0.3-0.9) falls slightly below their z=0.84 result. The evolution in both the H-alpha luminosity function from z=0.3--1.5 and the [OIII] luminosity function from z=0.7-2.3 is almost entirely in the L* parameter, which steadily increases with redshift over those ranges.

preprint2013arXiv

The Lyman alpha Reference Sample: Extended Lyman alpha Halos Produced at Low Dust Content

We report on new imaging observations of the Lyman alpha emission line (Lya), performed with the Hubble Space Telescope, that comprise the backbone of the Lyman alpha Reference Sample (LARS). We present images of 14 starburst galaxies at redshifts 0.028 < z < 0.18 in continuum-subtracted Lya, Halpha, and the far ultraviolet continuum. We show that Lya is emitted on scales that systematically exceed those of the massive stellar population and recombination nebulae: as measured by the Petrosian 20 percent radius, RP20, Lya radii are larger than those of Halpha by factors ranging from 1 to 3.6, with an average of 2.4. The average ratio of Lya-to-FUV radii is 2.9. This suggests that much of the Lya light is pushed to large radii by resonance scattering. Defining the "Relative Petrosian Extension" of Lya compared to Halpha, ξ_ext = RP20_Lya / RP20_Ha, we find ξ_ext to be uncorrelated with total Lya luminosity. However ξ_ext is strongly correlated with quantities that scale with dust content, in the sense that a low dust abundance is a necessary requirement (although not the only one) in order to spread Lya photons throughout the interstellar medium and drive a large extended Lya halo.

preprint2013arXiv

The Lyman alpha Reference Sample: II. HST imaging results, integrated properties and trends

We report upon new results regarding the Lya output of galaxies, derived from the Lyman alpha Reference Sample (LARS), focusing on Hubble Space Telescope imaging. For 14 galaxies we present intensity images in Lya, Halpha, and UV, and maps of Halpha/Hbeta, Lya equivalent width (EW), and Lya/Halpha. We present Lya and UV light profiles and show they are well-fitted by Sérsic profiles, but Lya profiles show indices systematically lower than those of the UV (n approx 1-2 instead of >~4). This reveals a general lack of the central concentration in Lya that is ubiquitous in the UV. Photometric growth curves increase more slowly for Lya than the FUV, showing that small apertures may underestimate the EW. For most galaxies, however, flux and EW curves flatten by radii ~10 kpc, suggesting that if placed at high-z, only a few of our galaxies would suffer from large flux losses. We compute global properties of the sample in large apertures, and show total luminosities to be independent of all other quantities. Normalized Lya throughput, however, shows significant correlations: escape is found to be higher in galaxies of lower star formation rate, dust content, mass, and several quantities that suggest harder ionizing continuum and lower metallicity. Eight galaxies could be selected as high-z Lya emitters, based upon their luminosity and EW. We discuss the results in the context of high-z Lya and UV samples. A few galaxies have EWs above 50 AA, and one shows f_escLya of 80%; such objects have not previously been reported at low-z.

preprint2013arXiv

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

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

preprint2012arXiv

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

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

preprint2012arXiv

Discovery of Three Distant, Cold Brown Dwarfs in the WFC3 Infrared Spectroscopic Parallels Survey

We present the discovery of three late type (>T4) brown dwarfs, including a probable Y dwarf, in the WFC3 Infrared Spectroscopic Parallels (WISP) Survey. We use the G141 grism spectra to determine the spectral types of the dwarfs and derive distance estimates based on a comparison with nearby T dwarfs with known parallaxes. These are the most distant spectroscopically confirmed T/Y dwarfs, with the farthest at an estimated distance of ~400 pc. We compare the number of cold dwarfs found in the WISP survey with simulations of the brown dwarf mass function. The number found is generally consistent with an initial stellar mass function dN/dM \propto M^{-α} with α= 0.0--0.5, although the identification of a Y dwarf is somewhat surprising and may be indicative of either a flatter absolute magnitude/spectral type relation than previously reported or an upturn in the number of very late type brown dwarfs in the observed volume.

preprint2011arXiv

On the redshift-evolution of the Lyman-alpha escape fraction and the dust content of galaxies

The Ly-alpha emission line has been proven a powerful tool by which to study evolving galaxies at the highest redshifts. However, in order to use Lya as a physical probe of galaxies, it becomes vital to know the Lya escape fraction (fescLya). Unfortunately, due to the resonant nature of Lya, fescLya may vary unpredictably and requires empirical measurement. Here we compile Lya luminosity functions between redshift z=0 and 8 and, combined with H-alpha and ultraviolet data, assess how fescLya evolves with redshift. We find a strong upwards evolution in fescLya over the range z=0.3-6, which is well-fit by the power-law fescLya \propto (1+z)^ξ, with ξ=(2.57_-0.12^+0.19). This predicts that fescLya should reach unity at z=11.1. By comparing fescLya and E(B-V) in individual galaxies we derive an empirical relationship between fescLya and E(B-V), which includes resonance scattering and can explain the redshift evolution of fescLya between z=0 and 6 purely as a function of the evolution in the dust content of galaxies. Beyond z~6.5, fescLya drops more substantially; an effect attributed to either ionizing photon leakage, or an increase in the neutral gas fraction of the intergalactic medium. While distinguishing between those two scenarios may be extremely challenging, by framing the problem this way we remove the uncertainty of the halo mass from Lya-based tests of reionization. We finally derive a new method by which to estimate the dust content of galaxies based purely upon the observed Lya and UV LFs. These data are characterized by an exponential with an e-folding redshift of ~3.5.

preprint2010arXiv

Escape of about five per cent of Lyman-alpha photons from high-redshift star-forming galaxies

The Lyman-alpha (Lya) emission line is the primary observational signature of star-forming galaxies at the highest redshifts, and has enabled the compilation of large samples of galaxies with which to study cosmic evolution. The resonant nature of the line, however, means that Lya photons scatter in the neutral interstellar medium of their host galaxies, and their sensitivity to absorption by interstellar dust may therefore be enhanced greatly. This implies that the Lya luminosity may be significantly reduced, or even completely suppressed. Hitherto, no unbiased empirical test of the escaping fraction (f_esc) of Lya photons has been performed at high redshifts. Here we report that the average fesc from star-forming galaxies at redshift z = 2.2 is just 5 per cent by performing a blind narrowband survey in Lya and Ha. This implies that numerous conclusions based on Lya-selected samples will require upwards revision by an order of magnitude and we provide a benchmark for this revision. We demonstrate that almost 90 per cent of star-forming galaxies emit insufficient Lya to be detected by standard selection criteria. Both samples show an anti-correlation of fesc with dust content, and we show that Lya- and Ha-selection recovers populations that differ substantially in dust content and fesc.

preprint2009arXiv

Empirical Estimate of Lyman-alpha Escape Fraction in a Statistical Sample of Lyman-alpha Emitters

The Lyman-alpha (Lya) recombination line is a fundamental tool for galaxy evolution studies and modern observational cosmology. However, subsequent interpretations are still prone to a number of uncertainties. Besides numerical efforts, empirical data are urgently needed for a better understanding of Lya escape process. We empirically estimate the Lyman-alpha escape fraction fesc(Lya) in a statistically significant sample of z ~ 0 - 0.3 galaxies in order to calibrate high-redshift Lyman-alpha observations. An optical spectroscopic follow-up of a sub-sample of 24 Lyman-alpha emitters (LAEs) detected by GALEX at z ~ 0.2-0.3, combined with a UV-optical sample of local starbursts, both with matched apertures, allow us to quantify the dust extinction through Balmer lines, and to estimate the Lyman-alpha escape fraction from the Halpha flux corrected for extinction in the framework of the recombination theory. The global escape fraction of Lyman-alpha radiation spans nearly the entire range of values, from 0.5 to 100 %, and fesc(Lya) clearly decreases with increasing nebular dust extinction E(B-V). Several objects show fesc(Lya) greater than fesc(continuum) which may be an observational evidence for clumpy ISM geometry or for an aspherical ISM. Selection biases and aperture size effects may still prevail between z ~ 0.2-0.3 LAEs and local starbursts, which may explain the difference observed for fesc(Lya).

preprint2009arXiv

On The Origin of Lyman-alpha Absorption in Nearby Starbursts and Implications for Other Galaxies

(Abridged) Despite the privileged position that Lyman-alpha (Lya) emission line holds in the exploration of the distant universe and modern observational cosmology, the origin of the observed diversity of lya profiles remains to be thoroughly explained. Observations of nearby star forming galaxies bring their batch of apparent contradictions between Lya emission and their physical parameters, and call for a detailed understanding of the physical processes at work. IZw 18, one of the most metal-poor galaxies known is of particular interest in this context. We use a 3D Lya radiation transfer code to model Hubble Space Telescope (HST) observations of IZw 18 and to fit its Lya spectrum. Different geometrical configurations of the source and the neutral gas are explored. The integrated Lya profile of NW region of IZw 18 is reproduced using the observed small amount of dust (E(B-V) ~ 0.05) and a spherical HI shell with N(HI) = 6.5 x 10^(21) cm^(-2). Such a high column density makes it possible to transform a strong Lya emission (EW(Lya) = 60 A) into a damped absorption even with a small extinction. When a slab geometry is applied and a given line of sight is chosen, the Lya profile can be successfully reproduced with no dust at all and N(HI) = 3 x 10^(21) cm^(-2). The spatial variations of the profile shape are naturally explained by radiation transfer effects. In the case of outflowing Inter Stellar Medium (ISM), as commonly observed in Lyman Break Galaxies (LBGs), a high N(H) and dust content are required to observe Lya in absorption. For nearly static neutral gas as observed in IZw 18 and other low luminosity galaxies only a small amount of dust is required provided a sufficiently high N(H) covers the galaxy.

preprint2009arXiv

The Lyman alpha morphology of local starburst galaxies: release of calibrated images

We present reduced and calibrated high resolution Lyman-alpha (Lya) images for a sample of 6 local star forming galaxies. Targets were selected to represent a range in luminosity and metallicity and to include both known Lya emitters and non-emitters. Far ultraviolet imaging was carried out with the Solar Blind Channel of the ACS on HST in the F122M (Lya on-line) and F140LP (continuum) filters. The resulting Lya images are the product of careful modeling of both the stellar and nebular continua, facilitated by supporting HST imaging in Ha and 5 continuum bands, combined with Starburst99 evolutionary synthesis models, and prescriptions for dust extinction on the continuum. In all, the resulting morphologies in Lya, Ha, and UV-continuum are qualitatively very different and we show that the bulk of Lya emerges in a diffuse component resulting from resonant scattering events. Lya escape fractions, computed from integrated Ha luminosities and recombination theory, are found never to exceed 14%. Internal dust extinction is estimated in each pixel and used to correct Lya fluxes. However, the extinction corrections are far too small (factors from 2.6 to infinity) to reconcile the global Lya luminosities with recombination theory. Surprisingly, when comparing the global equivalent widths of Lya and Ha, the two quantities appear anti-correlated, which may be due to the evolution of mechanical feedback. This calls for caution in the interpretation of Lya observations. The images presented have a physical resolution 3 orders of magnitude better than attainable at high-z from the ground with current instrumentation and our images may therefore serve as useful templates for comparing with observations and modeling of primeval galaxy formation. We therefore provide the reduced Lya, Ha, and continuum images to the community.