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.Comment: 11 pages, 10 figures. Accepted for publication in Astronomy and Astrophysic

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 (WFC3) 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. The high sensitivity of these grism observations, with >1σ detections of emission lines to f > 2.5 × 10^(–18) erg s^(–1) cm^(–2), means that the galaxies in the sample are typically ~7 times less massive (median M_* = 10^(9.5) M_☉) than previously studied z ~ 2 emission-line galaxies. Despite their lower mass, the galaxies have [O III]/Hβ 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 [O III] emission line is more spatially concentrated than the Hβ emission line with 98.1% confidence. We additionally stack the X-ray data (all sources are individually undetected), and find that the average L_([O III])/L_(0.5-10keV) ratio is intermediate between typical z ~ 0 obscured active galaxies and star-forming galaxies. Together the compactness of the stacked [O III] 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

High-redshift galaxies and black holes in the eyes of JWST: a population synthesis model from infrared to X-rays

The first billion years of the Universe is a pivotal time: stars, black holes (BHs) and galaxies form and assemble, sowing the seeds of galaxies as we know them today. Detecting, identifying and understand- ing the first galaxies and BHs is one of the current observational and theoretical challenges in galaxy formation. In this paper we present a population synthesis model aimed at galaxies, BHs and Active Galactic Nuclei (AGNs) at high redshift. The model builds a population based on empirical relations. Galaxies are characterized by a spectral energy distribution determined by age and metallicity, and AGNs by a spectral energy distribution determined by BH mass and accretion rate. We validate the model against observational constraints, and then predict properties of galaxies and AGN in other wavelength and/or luminosity ranges, estimating the contamination of stellar populations (normal stars and high-mass X-ray binaries) for AGN searches from the infrared to X-rays, and vice-versa for galaxy searches. For high-redshift galaxies, with stellar ages < 1 Gyr, we find that disentangling stellar and AGN emission is challenging at restframe UV/optical wavelengths, while high-mass X-ray binaries become more important sources of confusion in X-rays. We propose a color-color selection in JWST bands to separate AGN vs star-dominated galaxies in photometric observations. We also esti- mate the AGN contribution, with respect to massive, hot, metal-poor stars, at driving high ionization lines, such as C IV and He II. Finally, we test the influence of the minimum BH mass and occupa- tion fraction of BHs in low mass galaxies on the restframe UV/near-IR and X-ray AGN luminosity function.Comment: Accepted for publication in the Astrophysical Journa

Spitzer Infrared Properties of Lyman α Emitters

We present results on the infrared properties of Lyman α emitters. At z=0.3 we demonstrate that our sample of 50 GALEX-identified emitters possess the same extinction properties as similarly ultraviolet-bright galaxies at the same redshift. The AGN fraction for these z=0.3 sources is low, 10-20%, and they show no correlation between Lyα flux and total infrared luminosity. At z=2.4-3.1 we discuss a sample of Lyman α blobs. Roughly two thirds of these blobs have infrared counterparts with infrared and sub-mm colors indicating they are dominated by star formation. Of the six IRS spectra taken of blob counterparts, four show PAH features. Their equivalent widths indicate two are strongly star formation dominated, while the other two are mixed sources with a significant contribution from an AGN

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.Comment: 17 pages. Accepted for publicatio

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).Comment: 4 pages, 2 figures, accepted for publication in Astronomy and Astrophysic

The escape of Lyman photons from a young starburst: the case of Haro 11

(Abridged) Lyman-alpha (Lya) is a dominant probe of the galaxy population at high-z. However, interpretation of data drawn from Lya alone hinges on the Lya escape fraction which, due to the complex radiative transport, may vary greatly. Here we map the Lya emission from local starburst Haro 11, a Lya emitter and the only known candidate for low-z Lyman continuum emission (LyC). To aid in the interpretation we perform a detailed multi-wavelength analysis and model the stellar population, dust distribution, ionising photon budget, and star-cluster population. We use archival X-ray observations to further constrain properties of the starburst and estimate the HI column density. The Lya morphology is found to be strongly decoupled from stellar and nebular (H-alpha) morphologies. General surface photometry finds only very slight correlation between Lya and H-halpha, E(B-V), and stellar age. Only around the central Lya-bright cluster do we find the Lya/Ha ratio at values predicted by recombination theory. The total Lya escape fraction is found to be just 3%. We compute that ~90% of the Lya photons that escape do so after undergoing multiple resonance scattering events, masking their point of origin. This leads to a largely symmetric distribution and, by increasing the distance that photons must travel to escape, decreases the escape probability significantly. While dust must ultimately be responsible for the destruction of Lya, it plays little role in governing the observed morphology, which is regulated more by ISM kinematics and geometry. We find tentative evidence for local Lya equivalent width in the immediate vicinity of star-clusters being a function of cluster age, consistent with hydrodynamic studies. We estimate the ionising photon production and further constrain the escape fraction at 900 AA to <~9% .Comment: In press for MNRAS. 18 pages, 9 figures. Version with full resolution images to be found at http://www.astro.su.se/~matthew/english/papers/hayes_lya_haro11_hires.pd

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.Comment: 18 pages, 14 figures, Accepted by Ap

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_*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$.Comment: Accepted to MNRA

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^{-\alpha} with \alpha = 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.Comment: Accepted for publication by ApJ Letters. 10 pages, 2 figure