Physical properties and evolutionary state of the Lyman alpha emitting starburst galaxy IRAS 08339+6517

Though Lyman alpha emission (Lya) is one of the most used tracers of massive star formation at high redshift, a correct understanding of radiation transfer effects by neutral gas is required to properly quantify the star formation rate along the history of the Universe. We are embarked in a program to study the properties of the Lya emission (spectral profile, spatial distribution, relation to Balmer lines intensity,...) in several local starburst galaxies. We present here the results obtained for IRAS 08339+6517. Using evolutionary population synthesis models, we have characterized the properties of the starburst (UV continuum, Halpha, total infrared and X-ray emissions, etc.), which transformed 1.4e+8 Mo of gas into stars around 5-6 Myr ago. In addition to the central compact emission blob, we have identified a diffuse Lya emission component smoothly distributed over the whole central area of IRAS 08339+6517. This diffuse emission is spatially decoupled from the UV continuum, the Halpha emission or the Halpha/Hbeta ratio. Both locally and globally, the Lya/Halpha ratio is lower than the Case B predictions, even after reddening correction, with an overall Lya escape fraction of only 4%. We conclude that in IRAS 08339+6517 the resonant scattering of Lya photons by an outflowing shell of neutral gas causes their highly-efficient destruction by dust, which explains the low Lya escape fraction measured. These results stress again the importance of a proper correction of scattering and transfer effects when using Lya to derive the star formation rate in high-redshift galaxies.Comment: Accepted for publication in A&A, 17 pages, 13 figures, 8 tables. If problems with quality of images, see https://cloud.cab.inta-csic.es/public.php?service=files&file=%2Fotih%2Ffiles%2Foti_mas%2Firas%2Firas-v53.ack_referee.pd

Multiwavelength analysis of the Lyman alpha emitting galaxy Haro 2: relation between the diffuse Lyman alpha and soft X-ray emissions

In order to use Lyman alpha (Lya) emission as star formation tracer in cosmological studies, we must understand how the resonant scattering affects the escape fraction of the Lya photons. Thus, high spatial resolution multiwavelength studies of nearby Lya emitters, like Haro 2, are highly needed. For that purpose, we have used Chandra X-ray and HST (UV, optical and NIR) images of Haro 2, and STIS and ground-based spectral images along its major and minor axes, to characterize the Lya emission and the properties of the stellar population. The UV, Ha (Halpha) and FIR luminosities of the Haro 2 nuclear starburst are reproduced using evolutionary synthesis models assuming a young stellar population with ages ~3.5-5.0 Myr, affected by differential interstellar extinctions. The observed X-ray emission is attributed to gas heated by the mechanical energy released by the starburst (soft component) and a Ultra-Luminous X-ray source candidate (hard). Both compact and diffuse Lya components are observed. Whereas Lya is spatially decoupled from Balmer lines emission, Balmer decrement and UV continuum, the diffuse Lya component is spatially correlated with the diffuse soft X-ray emission. Moreover, unlike the compact one, diffuse Lya shows luminosities larger than predicted from Ha, assuming case B recombination and dust extinction as derived from Ha/Hbeta. We propose that, whereas the compact Lya emission is associated to the massive stellar clusters and is affected by outflows and dust extinction, the diffuse Lya originates in gas ionized by the hot plasma responsible for the soft X-ray radiation, as suggested by their spatial correlation and by the measured L(Ha)/LsoftX ratios. Calibration of Lya as star formation rate tracer should therefore include both effects (destruction vs. enhancement) to avoid biases in the study of galaxies at cosmological distances.Comment: Accepted for publication in A&A, 18 pages, 8 figures, 9 tables. If problems with quality of images, see http://www.cab.inta-csic.es/users/otih/haro2-v63.clean.pd

Extremely Small Sizes for Faint z~2-8 Galaxies in the Hubble Frontier Fields: A Key Input For Establishing their Volume Density and UV Emissivity

We provide the first observational constraints on the sizes of the faintest galaxies lensed by the Hubble Frontier Fields (HFF) clusters. Ionizing radiation from faint galaxies likely drives cosmic reionization, and the HFF initiative provides a key opportunity to find such galaxies. Yet, we cannot really assess their ionizing emissivity without a robust measurement of their sizes, since this is key to quantifying both their prevalence and the faint-end slope to the UV luminosity function. Here we provide the first such size constraints with 2 new techniques. The first utilizes the fact that the detectability of highly-magnified galaxies as a function of shear is very dependent on a galaxy's size. Only the most compact galaxies will remain detectable in regions of high shear (vs. a larger detectable size range for low shear), a phenomenon we carefully quantify using simulations. Remarkably, however, no correlation is found between the surface density of faint galaxies and the predicted shear, using 87 faint high-magnification mu>10 z~2-8 galaxies seen behind the first 4 HFF clusters. This can only be the case if such faint (~-15 mag) galaxies have significantly smaller sizes than luminous galaxies. We constrain their half-light radii to be <~30 mas (<160-240 pc). As a 2nd size probe, we rotate and stack 26 faint high-magnification sources along the major shear axis. Less elongation is found than even for objects with an intrinsic half-light radius of 10 mas. Together these results indicate that extremely faint z~2-8 galaxies have near point-source profiles in the HFF dataset (half-light radii conservatively <30 mas and likely 5-10 mas). These results suggest smaller completeness corrections and hence much lower volume densities for faint z~2-8 galaxies and shallower faint-end slopes than have been derived in many recent studies (by factors of ~2-3 and by dalpha>~0.1-0.3).Comment: 19 pages, 15 figures, 3 tables, accepted for publication in Ap

Grid of Lya radiation transfer models for the interpretation of distant galaxies

Lya is a key diagnostic for numerous observations of distant star-forming galaxies. It's interpretation requires, however, detailed radiation transfer models. We provide an extensive grid of 3D radiation transfer models simulating the Lya and UV continuum radiation transfer in the interstellar medium of star-forming galaxies. We have improved our Monte Carlo MCLya code, and have used it to compute a grid of 6240 radiation transfer models for homogeneous spherical shells containing HI and dust surrounding a central source. The simulations cover a wide range of parameter space. We present the detailed predictions from our models including in particular the Lya escape fraction fesc, the continuum attenuation, and detailed Lya line profiles. The Lya escape fraction is shown to depend strongly on dust content, but also on other parameters (HI column density and radial velocity). The predicted line profiles show a great diversity of morphologies ranging from broad absorption lines to emission lines with complex features. The results from our simulations are distributed in electronic format. Our models should be of use for the interpretation of observations from distant galaxies, for other simulations, and should also serve as an important base for comparison for future, more refined, radiation transfer models.Comment: Accepted for publication in Astronomy & Astrophysics. Results from simulations available at http://obswww.unige.ch/sf

The Lyman Alpha Reference Sample. VIII. Characterizing Lyman-Alpha Scattering in Nearby Galaxies

We examine the dust geometry and Ly{\alpha} scattering in the galaxies of the Lyman Alpha Reference Sample (LARS), a set of 14 nearby (0.02 < $z$ < 0.2) Ly{\alpha} emitting and starbursting systems with Hubble Space Telescope Ly{\alpha}, H{\alpha}, and H{\beta} imaging. We find that the global dust properties determined by line ratios are consistent with other studies, with some of the LARS galaxies exhibiting clumpy dust media while others of them show significantly lower Ly{\alpha} emission compared to their Balmer decrement. With the LARS imaging, we present Ly{\alpha}/H{\alpha} and H{\alpha}/H{\beta} maps with spatial resolutions as low as $\sim$ 40 pc, and use these data to show that in most galaxies, the dust geometry is best modeled by three distinct regions: a central core where dust acts as a screen, an annulus where dust is distributed in clumps, and an outer envelope where Ly{\alpha} photons only scatter. We show that the dust that affects the escape of Ly{\alpha} is more restricted to the galaxies' central regions, while the larger Ly{\alpha} halos are generated by scattering at large radii. We present an empirical modeling technique to quantify how much Ly{\alpha} scatters in the halo, and find that this "characteristic" scattering distance correlates with the measured size of the Ly{\alpha} halo. We note that there exists a slight anti-correlation between the scattering distance of Ly{\alpha} and global dust properties.Comment: 32 pages, 51 figures, accepted to Ap

Dust extinction from Balmer decrements of star-forming galaxies at 0.75<z<1.5 with HST/WFC3 spectroscopy from the WISP survey

Spectroscopic observations of Halpha and Hbeta emission lines of 128 star-forming galaxies in the redshift range 0.75<z<1.5 are presented. These data were taken with slitless spectroscopy using the G102 and G141 grisms of the Wide-Field-Camera 3 (WFC3) on board the Hubble Space Telescope as part of the WFC3 Infrared Spectroscopic Parallel (WISP) survey. Interstellar dust extinction is measured from stacked spectra that cover the Balmer decrement (Halpha/Hbeta). We present dust extinction as a function of Halpha luminosity (down to 3 x 10^{41} erg/s), galaxy stellar mass (reaching 4 x 10^{8} Msun), and rest-frame Halpha equivalent width. The faintest galaxies are two times fainter in Halpha luminosity than galaxies previously studied at z~1.5. An evolution is observed where galaxies of the same Halpha luminosity have lower extinction at higher redshifts, whereas no evolution is found within our error bars with stellar mass. The lower Halpha luminosity galaxies in our sample are found to be consistent with no dust extinction. We find an anti-correlation of the [OIII]5007/Halpha flux ratio as a function of luminosity where galaxies with L_{Halpha}<5 x 10^{41} erg/s are brighter in [OIII]5007 than Halpha. This trend is evident even after extinction correction, suggesting that the increased [OIII]5007/Halpha ratio in low luminosity galaxies is likely due to lower metallicity and/or higher ionization parameters.Comment: 11 pages, 9 figures, 2 tables; version addressing the referee comment

The Lyman alpha reference sample. VII. Spatially resolved Hα\alpha 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 $\alpha$ (H$\alpha$) emission in our data cubes. These maps trace the spectral and spatial properties of the LARS galaxies' intrinsic Ly$\alpha$ radiation field. We show our kinematic maps spatially registered onto the Hubble Space Telescope H$\alpha$ and Lyman $\alpha$ (Ly$\alpha$) images. Only for individual galaxies a causal connection between spatially resolved H$\alpha$ kinematics and Ly$\alpha$ 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 $\sigma_0$, shearing velocity $v_\mathrm{shear}$, and the $v_\mathrm{shear}/\sigma_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}/\sigma_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}/\sigma_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$\alpha$ properties, we find that dispersion dominated systems show higher Ly$\alpha$ equivalent widths and higher Ly$\alpha$ escape fractions than systems with $v_\mathrm{shear}/\sigma_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$\alpha$ radiation.Comment: 26 pages, 15 figures, accepted for publication in A&

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, \xi_ext = RP20_Lya / RP20_Ha, we find \xi_ext to be uncorrelated with total Lya luminosity. However \xi_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.Comment: Published in ApJ Letters ~~ 6 pages using emulateapj, 4 figures ~~ Higher-resolution, larger, nicer jpeg versions of Figures 1 and 2 can be found here: http://xayes.org/pub/press_lars.htm

Overdensities of Y-dropout Galaxies from the Brightest-of-Reionizing Galaxies Survey: A Candidate Protocluster at Redshift z~8

Theoretical and numerical modeling of dark-matter halo assembly predicts that the most luminous galaxies at high redshift are surrounded by overdensities of fainter companions. We test this prediction with HST observations acquired by our Brightest of Reionizing Galaxies (BoRG) survey, which identified four very bright z~8 candidates as Y-dropout sources in four of the 23 non-contiguous WFC3 fields observed. We extend here the search for Y-dropouts to fainter luminosities (M_* galaxies with M_AB\sim-20), with detections at >5sigma confidence (compared to >8sigma confidence adopted earlier) identifying 17 new candidates. We demonstrate that there is a correlation between number counts of faint and bright Y-dropouts at >99.84% confidence. Field BoRG58, which contains the best bright z\sim8 candidate (M_AB=-21.3), has the most significant overdensity of faint Y-dropouts. Four new sources are located within 70arcsec (corresponding to 3.1 comoving Mpc at z=8) from the previously known brighter z\sim8 candidate. The overdensity of Y-dropouts in this field has a physical origin to high confidence (p>99.975%), independent of completeness and contamination rate of the Y-dropout selection. We modeled the overdensity by means of cosmological simulations and estimate that the principal dark matter halo has mass M_h\sim(4-7)x10^11Msun (\sim5sigma density peak) and is surrounded by several M_h\sim10^11Msun halos which could host the fainter dropouts. In this scenario, we predict that all halos will eventually merge into a M_h>2x10^14Msun galaxy cluster by z=0. Follow-up observations with ground and space based telescopes are required to secure the z\sim8 nature of the overdensity, discover new members, and measure their precise redshift.Comment: Minor revision: ApJ accepted [17 pages (emulateapj style), 7 figures, 2 tables

The Hubble Legacy Archive ACS Grism Data

A public release of slitless spectra, obtained with ACS/WFC and the G800L grism, is presented. Spectra were automatically extracted in a uniform way from 153 archival fields (or "associations") distributed across the two Galactic caps, covering all observations to 2008. The ACS G800L grism provides a wavelength range of 0.55-1.00 \mu$m, with a dispersion of$40 \ \AA / pixel$and a resolution of$\sim 80\ \AA$for point-like sources. The ACS G800L images and matched direct images were reduced with an automatic pipeline that handles all steps from archive retrieval, alignment and astrometric calibration, direct image combination, catalogue generation, spectral extraction and collection of metadata. The large number of extracted spectra (73,581) demanded automatic methods for quality control and an automated classification algorithm was trained on the visual inspection of several thousand spectra. The final sample of quality controlled spectra includes 47,919 datasets (65of extracted spectra) for$32,149$unique objects, with a median$i_{\rm AB}$-band magnitude of 23.7, reaching 26.5 AB for the faintest objects. Each released dataset contains science-ready 1D and 2D spectra, as well as multi-band image cutouts of corresponding sources and a useful preview page summarising the direct and slitless data, astrometric and photometric parameters. In order to characterize the slitless spectra, emission-line flux and equivalent width sensitivity of the ACS data were compared with public ground-based spectra in the GOODS-South field. An example list of emission line galaxies with two or more identified lines is also included, covering the redshift range$0.2-4.6$.Comment: Accepted for publication in Astronomy and Astrophysics; 29 pages, 16 Figures, 4 Tables in text and 3Tables in Appendi