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 Sersic 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.

I Zw 18 as morphological paradigm for rapidly assembling high-z galaxies

IZw18, ever since regarded as the prototypical blue compact dwarf (BCD) galaxy, is, quite ironically, the most atypical BCD known. This is because its large exponential low-surface brightness envelope is not due to an old stellar host but entirely due to extended nebular emission (ne) (Papaderos et al. 2002; P02). We study IZw18 and IZw18C down to an unprecedently faint surface brightness level using HST ACS data. We argue that the properties of IZw18C can be consistently accounted for by propagating star formation over the past ~100 Myr, in combination with stellar diffusion and the associated radial stellar mass filtering effect (P02). As for IZw18, we find that ne extends out to ~16 stellar scale lengths and provides at least 1/3 of the total optical emission. The case of IZw18 suggests caution in studies of distant galaxies in dominant stages of their evolution, rapidly assembling their stellar mass at high specific star formation rates (SSFRs). It calls attention to the fact that ne is not necessarily cospatial with the underlying ionizing and non-ionizing stellar background, neither has to scale with its surface density. The prodigious energetic output during dominant phases of galaxy evolution may result in large exponential ne envelopes, extending much beyond the still compact stellar component, just like in IZw18. Therefore, the morphological paradigm of IZw18, while probably unique in the nearby Universe, may be ubiquitous among high-SSFR galaxies at high redshift. Using IZw18 as reference, we show that extended ne may introduce substantial observational biases and significantly affect fundamental galaxy relations. Among others, we show that the surface brightness profiles of distant morphological analogs to IZw18 may be barely distinguishable from Sersic profiles with an exponent 2<n<5, thus mimicking the profiles of massive galaxy spheroids. (abridged)

The nature of z ~ 2.3 Lyman-alpha emitters

We study the multi-wavelength properties of a set of 171 Ly-alpha emitting candidates at redshift z = 2.25 found in the COSMOS field, with the aim of understanding the underlying stellar populations in the galaxies. We especially seek to understand what the dust contents, ages and stellar masses of the galaxies are, and how they relate to similar properties of Ly-alpha emitters at other redshifts. The candidates here are shown to have different properties from those of Ly-alpha emitters found at higher redshift, by fitting the spectral energy distributions (SEDs) using a Monte-Carlo Markov-Chain technique and including nebular emission in the spectra. The stellar masses, and possibly the dust contents, are higher, with stellar masses in the range log M_* = 8.5 - 11.0 M_sun and A_V = 0.0 - 2.5 mag. Young population ages are well constrained, but the ages of older populations are typically unconstrained. In 15% of the galaxies only a single, young population of stars is observed. We show that the Ly-alpha fluxes of the best fit galaxies are correlated with their dust properties, with higher dust extinction in Ly-alpha faint galaxies. Testing for whether results derived from a light-weighted stack of objects correlate to those found when fitting individual objects we see that stellar masses are robust to stacking, but ages and especially dust extinctions are derived incorrectly from stacks. We conclude that the stellar properties of Ly-alpha emitters at z = 2.25 are different from those at higher redshift and that they are diverse. Ly-alpha selection appears to be tracing systematically different galaxies at different redshifts.

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.

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.

Continuum subtracting Lyman-alpha images: Low redshift studies using the Solar Blind Channel of HST/ACS

[ABRIDGED] We are undertaking an imaging study of local star-forming galaxies in the Lyman-alpha (Lya) emission line using the Solar Blind Channel (SBC) of the ACS onboard the HST. Observations have been obtained in Lya and H-alpha (Ha) and 6 continuum filters between ~1500AA and the I-band. Previously (Hayes et al 2005) we demonstrated that the production of Lya line-only images in the SBC-only data-set is non-trivial and that supporting data is a requirement. We here develop various methods of continuum subtraction and assess their relative performance using a variety of spectral energy distributions (SED) as input. We conclude that simple assumptions about the behavior of the ultraviolet continuum consistently lead to results that are wildly erroneous, and determine that an SED fitting approach is essential. Moreover, a single component stellar or stellar+nebular spectrum is not always sufficient for realistic template SEDs and, in order to successfully recover the input observables, care must be taken to control contributions of nebular gas and underlying stellar populations. Independent measurements of the metallicity must first be obtained, but details of the initial mass function play only a small role. We identify the need to bin together pixels in our data to obtain signal-to-noise (S/N) of ~10 in each band before processing. At S/N=10 we can recover Lya fluxes accurate to ~30% for Lya lines with intrinsic equivalent width [W(Lya)] of 10AA. This accuracy improves to ~10% for W(Lya)=100AA. We describe the image processing applied to the observations presented in Ostlin et al. (2009) and the associated data-release. We also present simulations for an observing strategy for an alternative low-redshift SBC Lya imaging campaign using adjacent combinations of long-pass filters at slightly higher redshift.

Stellar kinematics in blue compact galaxies

In a programme of observations of local luminous blue compact galaxies (BCGs), we are investigating kinematics by using tracers of both stars and ionized gas. Here we summarise our program and present new data on the local Lyman break galaxy analogue Haro 11. From spatially-resolved spectroscopy around the near-infrared Ca II triplet, we find that its stars and ionized gas have similar velocity fields. Our programme so far indicates however that emission line velocities can differ locally by a few tens of km/s from the Ca II values. Comparing our data to simple stellar population models, we assess which stellar population the Ca II triplet traces and its potential beyond the local universe.

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% .