Using Lyman-alpha to detect galaxies that leak Lyman continuum

We propose to infer ionising continuum leaking properties of galaxies by looking at their Lyman-alpha line profiles. We carry out Lyman-alpha radiation transfer calculations in two models of HII regions which are porous to ionising continuum escape: 1) the so-called "density bounded" media, in which massive stars produce enough ionising photons to keep the surrounding interstellar medium transparent to the ionising continuum, i.e almost totally ionised, and 2) "riddled ionisation-bounded" media, surrounded by neutral interstellar medium, but with holes, i.e. with a covering factor lower than unity. The Lyman-alpha spectra emergent from these configurations have distinctive features: 1) a "classical" asymmetric redshifted profile in the first case, but with a small shift of the maximum of the profile compare to the systemic redshift (Vpeak < 150 km/s); 2) a main peak at the systemic redshift in the second case (Vpeak = 0 km/s), with, as a consequence, a non-zero Lyman-alpha flux bluewards the systemic redshift. Assuming that in a galaxy leaking ionising photons, the Lyman-alpha component emerging from the leaking star cluster(s) dominates the total Lyman-alpha spectrum, the Lyman-alpha shape may be used as a pre-selection tool to detect Lyman continuum (LyC) leaking galaxies, in objects with well determined systemic redshift, and high spectral resolution Lyman-alpha spectra (R >= 4000). The examination of a sample of 10 local starbursts with high resolution HST-COS Lyman-alpha spectra and known in the literature as LyC leakers or leaking candidates, corroborates our predictions. Observations of Lyman-alpha profiles at high resolution should show definite signatures revealing the escape of Lyman continuum photons from star-forming galaxies.Comment: A&A in pres

Confidence limits of evolutionary synthesis models III. On time-integrated quantities

Evolutionary synthesis models are a fundamental tool to interpret the properties of observed stellar systems. In order to achieve a meaningful comparison between models and real data, it is necessary to calibrate the models themselves, i.e. to evaluate the dispersion due to the discreteness of star formation as well as the possible model errors. In this paper we show that linear interpolations in the log M - log t_k plane, that are customary in the evaluation of isochrones in evolutionary synthesis codes, produce unphysical results. We also show that some of the methods used in the calculation of time-integrated quantities (kinetic energy, and total ejected masses of different elements) may produce unrealistic results. We propose alternative solutions to solve both problems. Moreover, we have quantified the expected dispersion of these quantities due to stochastic effects in stellar populations. As a particular result, we show that the dispersion in the 14N/12C ratio increases with time.Comment: 11 pages, 8 figures, accepted by A&

The Massive Stellar Content in the Starburst NGC3049: A Test for Hot-Star Mode

We have obtained high-spatial resolution ultraviolet and optical STIS spectroscopy and imaging of the metal-rich nuclear starburst in NGC3049. The stellar continuum and the absorption line spectrum in the ultraviolet are used to constrain the massive stellar population. The strong, blueshifted stellar lines of CIV and SiIV detected in the UV spectra indicate a metal-rich, compact, massive (1E6 Msol) cluster of age 3--4 Myr emitting the UV-optical continuum. We find strong evidence against a depletion of massive stars in this metal-rich cluster. The derived age and the upper mass-limit cut-off of the initial mass function are also consistent with the detection of Wolf-Rayet (WR) features at optical wavelengths. As a second independentconstraint on the massive stellar content, the nebular emission-line spectrum is modeled with photoionization codes using stellar spectra from evolutionary synthesis models. However, the nebular lines are badly reproduced by 3--4 Myr instantaneous bursts, as required by the UV line spectrum, when unblanketed WR and/or Kurucz stellar atmospheres are used. The corresponding number of photons above 24 and 54 eV in the synthetic models is too high in comparison with values suggested by the observed line ratios. Since the ionizing spectrum in this regime is dominated by emission from WR stars, this discrepancy between observations and models is most likely the result of incorrect assumptions about the WR stars. Thus we conclude that the nebular spectrum of high-metallicity starbursts is poorly reproduced by models for WR dominated populations. However, the new model set of Smith et al. (2002) with blanketed WR and O atmospheres and adjusted WR temperatures predicts a softer far-UV radiation field, providing a better match to the data.Comment: To be published in ApJ, Dec. issue 17 figures, 3 in gif forma

New constraints on dust emission and UV attenuation of z=6.5-7.5 galaxies from millimeter observations

We have targeted two recently discovered Lyman break galaxies (LBGs) to search for dust continuum and [CII] 158 micron line emission. The strongly lensed z~6.8 LBG A1703-zD1 behind the galaxy cluster Abell 1703, and the spectroscopically confirmed z=7.508 LBG z8-GND-5296 in the GOODS-N field have been observed with the Plateau de Bure interferometer (PdBI) at 1.2mm. These observations have been combined with those of three z>6.5 Lya emitters (named HCM6A, Himiko, and IOK-1), for which deep measurements were recently obtained with the PdBI and ALMA. [CII] is undetected in both galaxies, providing a deep upper limit for Abell1703-zD1, comparable to recent ALMA non-detections. Dust continuum emission from Abell1703-zD1 and z8-GND-5296 is not detected with an rms of 0.12 and 0.16 mJy/beam. From these non-detections we derive upper limits on their IR luminosity and star formation rate, dust mass, and UV attenuation. Thanks to strong gravitational lensing the limit for Abell1703-zD1 is probing the sub-LIRG regime ($L_{IR} <8.1 \times 10^{10}$ Lsun) and very low dust masses ($M_d<1.6 \times 10^7$ Msun). We find that all five galaxies are compatible with the Calzetti IRX-$\beta$ relation, their UV attenuation is compatible with several indirect estimates from other methods (the UV slope, extrapolation of the attenuation measured from the IR/UV ratio at lower redshift, and SED fits), and the dust-to-stellar mass ratio is not incompatible with that of galaxies from z=0 to 3. For their stellar mass the high-z galaxies studied here have an attenuation below the one expected from the mean relation of low redshift (z<1.5) galaxies. More and deeper (sub)-mm data are clearly needed to directly determine the UV attenuation and dust content of the dominant population of high-z star-forming galaxies and to establish more firmly their dependence on stellar mass, redshift, and other properties.Comment: 10 pages, 7 figures. Minor revisions. Accepted for publication in A&

Starburst99: Synthesis Models for Galaxies with Active Star Formation

Starburst99 is a comprehensive set of model predictions for spectrophotometric and related properties of galaxies with active star formation. The models are an improved and extended version of the data set previously published by Leitherer & Heckman (1995). We have upgraded our code by implementing the latest set of stellar evolution models of the Geneva group and the model atmosphere grid compiled by Lejeune et al. (1997). Several predictions which were not included in the previous publication are shown here for the first time. The models are presented in a homogeneous way for five metallicities between Z = 0.040 and 0.001 and three choices of the initial mass function. The age coverage is 10^6 to 10^9 yr. We also show the spectral energy distributions which are used to compute colors and other quantities. The full data set is available for retrieval at http://www.stsci.edu/science/starburst99/. This website allows users to run specific models with non-standard parameters as well. We also make the source code available to the community.Comment: 32 pages, LaTeX. All the Figures and the summary Table are located at http://www.stsci.edu/science/starburst99/, ApJ accepte

HST's view of the youngest massive stars in the Magellanic Clouds

Accurate physical parameters of newborn massive stars are essential ingredients to shed light on their formation, which is still an unsolved problem. The rare class of compact H II regions in the Magellanic Clouds (MCs), termed ``high-excitation blobs'' (HEBs), presents a unique opportunity to acquire this information. These objects (~ 4" to 10", ~ 1 to 3 pc, in diameter) harbor the youngest massive stars of the OB association/molecular cloud complexes in the MCs accessible through high-resolution near-IR and optical techniques. We present a brief overview of the results obtained with HST mainly on two HEBs, one in the LMC (N159-5) and the other in the SMC (N81).Comment: 5 pages, to appear in the Proceedings of the 41st ESLAB Symposium "The Impact of HST on European Astronomy", 29 May to 1 June 2007, ESTEC, Noordwijk, Netherlands; eds. Guido De Marchi and Duccio Macchett

Neutral gas in Lyman-alpha emitting galaxies Haro 11 and ESO 338-IG04 measured through sodium absorption

Context. The Lyman alpha emission line of galaxies is an important tool for finding galaxies at high redshift, and thus probe the structure of the early universe. However, the resonance nature of the line and its sensitivity to dust and neutral gas is still not fully understood. Aims. We present measurements of the velocity, covering fraction and optical depth of neutral gas in front of two well known local blue compact galaxies that show Lyman alpha in emission: ESO 338-IG 04 and Haro 11. We thus test observationally the hypothesis that Lyman alpha can escape through neutral gas by being Doppler shifted out of resonance. Methods. We present integral field spectroscopy from the GIRAFFE/Argus spectrograph at VLT/FLAMES in Paranal, Chile. The excellent wavelength resolution allows us to accurately measure the velocity of the ionized and neutral gas through the H-alpha emission and Na D absorption, which traces the ionized medium and cold interstellar gas, respectively. We also present independent measurements with the VLT/X-shooter spectrograph which confirm our results. Results. For ESO 338-IG04, we measure no significant shift of neutral gas. The best fit velocity is -15 (16) km/s. For Haro 11, we see an outflow from knot B at 44 (13) km/s and infalling gas towards knot C with 32 (12) km/s. Based on the relative strength of the Na D absorption lines, we estimate low covering fractions of neutral gas (down to 10%) in all three cases. Conclusions. The Na D absorption likely occurs in dense clumps with higher column densities than where the bulk of the Ly-alpha scattering takes place. Still, we find no strong correlation between outflowing neutral gas and a high Lyman alpha escape fraction. The Lyman alpha photons from these two galaxies are therefore likely escaping due to a low column density and/or covering fraction.Comment: 9 pages, 3 figure