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&

Insights on star formation histories and physical properties of 1.2≤z≲41.2 \leq z \lesssim 4 Herschel-detected galaxies

We test the impact of using variable star forming histories (SFHs) and the use of the IR luminosity (LIR) as a constrain on the physical parameters of high redshift dusty star-forming galaxies. We explore in particular the stellar properties of galaxies in relation with their location on the SFR-M* diagram. We perform SED fitting of the UV-NIR and FIR emissions of a large sample of GOODS-Herschel galaxies, for which rich multi-wavelength observations are available. We test different SFHs and imposing energy conservation in the SED fitting process, to face issues like the age-extinction degeneracy and produce SEDs consistent with observations. Our models work well for the majority of the sample, with the notable exception of the high LIR end, for which we have indications that our simple energy conservation approach cannot hold true. We find trends in the SFHs fitting our sources depending on stellar mass M* and z. Trends also emerge in the characteristic timescales of the SED models depending on the location on the SFR-M* diagram. We show that whilst using the same available observational data, we can produce galaxies less star-forming than usually inferred, if we allow declining SFHs, while properly reproducing their observables. These sources can be post-starbursts undergoing quenching, and their SFRs are potentially overestimated if inferred from their LIR. Fitting without the IR constrain leads to a strong preference for declining SFHs, while its inclusion increases the preference of rising SFHs, more so at high z, in tentative agreement with the cosmic star formation history. Keeping in mind that the sample is biased towards high LIR, the evolution shaped by our model appears as both bursty (initially) and steady-lasting (later on). The global SFH of the sample follows the cosmic SFH with a small scatter, and is compatible with the "downsizing" scenario of galaxy evolution.Comment: 28 pages, 26 figures, one appendix, Accepted for publication in Astronomy & Astrophysic

Grids of stellar models. VIII. From 0.4 to 1.0 Msun at Z=0.020 and Z=0.001, with the MHD equation of state

We present stellar evolutionary models covering the mass range from 0.4 to 1 Msun calculated for metallicities Z=0.020 and 0.001 with the MHD equation of state (Hummer & Mihalas, 1988; Mihalas et al. 1988; D\"appen et al. 1988). A parallel calculation using the OPAL (Rogers et al. 1996) equation of state has been made to demonstrate the adequacy of the MHD equation of state in the range of 1.0 to 0.8 Msun (the lower end of the OPAL tables). Below, down to 0.4 Msun, we have justified the use of the MHD equation of state by theoretical arguments and the findings of Chabrier & Baraffe (1997). We use the radiative opacities by Iglesias & Rogers (1996), completed with the atomic and molecular opacities by Alexander & Fergusson (1994). We follow the evolution from the Hayashi fully convective configuration up to the red giant tip for the most massive stars, and up to an age of 20 Gyr for the less massive ones. We compare our solar-metallicity models with recent models computed by other groups and with observations. The present stellar models complete the set of grids computed with the same up-to-date input physics by the Geneva group [Z=0.020 and 0.001, Schaller et al. (1992), Bernasconi (1996), and Charbonnel et al. (1996); Z=0.008, Schaerer et al. (1992); Z=0.004, Charbonnel et al. (1993); Z=0.040, Schaerer et al. (1993); Z=0.10, Mowlavi et al. (1998); enhanced mass loss rate evolutionary tracks, Meynet et al. (1994)].Comment: Accepted for publication in A&A Supplement Serie

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

26Al yields from rotating Wolf--Rayet star models

We present new $^{26}$Al stellar yields from rotating Wolf--Rayet stellar models which, at solar metallicity, well reproduce the observed properties of the Wolf-Rayet populations. These new yields are enhanced with respect to non--rotating models, even with respect to non--rotating models computed with enhanced mass loss rates. We briefly discuss some implications of the use of these new yields for estimating the global contribution of Wolf-Rayet stars to the quantity of $^{26}$Al now present in the Milky Way.Comment: 6 pages, 2 figures, to appear in New Astronomy Review

Lyman-alpha spectral properties of five newly discovered Lyman continuum emitters

We have recently reported the discovery of five low redshift Lyman continuum (LyC) emitters (LCEs, hereafter) with absolute escape fractions fesc(LyC) ranging from 6 to 13%, higher than previously found, and which more than doubles the number of low redshift LCEs.We use these observations to test theoretical predictions about a link between the characteristics of the Lyman-alpha (Lya) line from galaxies and the escape of ionising photons. We analyse the Lya spectra of eight LCEs of the local Universe observed with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope (our five leakers and three galaxies from the litterature), and compare their strengths and shapes to the theoretical criteria and comparison samples of local galaxies: the Lyman Alpha Reference Survey, Lyman Break Analogs, Green Peas, and the high-redshift strong LyC leaker Ion2. Our LCEs are found to be strong Lya emitters, with high equivalent widths, EW(Lya)> 70 {\AA}, and large Lya escape fractions, fesc(Lya) > 20%. The Lya profiles are all double-peaked with a small peak separation, in agreement with our theoretical expectations. They also have no underlying absorption at the Lya position. All these characteristics are very different from the Lya properties of typical star-forming galaxies of the local Universe. A subset of the comparison samples (2-3 Green Pea galaxies) share these extreme values, indicating that they could also be leaking. We also find a strong correlation between the star formation rate surface density and the escape fraction of ionising photons, indicating that the compactness of star-forming regions plays a role in shaping low column density paths in the interstellar medium of LCEs. The Lya properties of LCEs are peculiar: Lya can be used as a reliable tracer of LyC escape from galaxies, in complement to other indirect diagnostics proposed in the literature.Comment: 11 pages, 10 figures, accepted for publication in A&

The Brightest Lyα\alpha Emitter: Pop III or Black Hole?

CR7 is the brightest $z=6.6 \, {\rm Ly}\alpha$ emitter (LAE) known to date, and spectroscopic follow-up by Sobral et al. (2015) suggests that CR7 might host Population (Pop) III stars. We examine this interpretation using cosmological hydrodynamical simulations. Several simulated galaxies show the same "Pop III wave" pattern observed in CR7. However, to reproduce the extreme CR7 ${\rm Ly}\alpha$/HeII1640 line luminosities ($L_{\rm \alpha/He II}$) a top-heavy IMF and a massive ($>10^{7}{\rm M}_{\odot}$) PopIII burst with age $<2$ Myr are required. Assuming that the observed properties of ${\rm Ly}\alpha$ and HeII emission are typical for Pop III, we predict that in the COSMOS/UDS/SA22 fields, 14 out of the 30 LAEs at $z=6.6$ with $L_{\alpha} >10^{43.3}{\rm erg}\,{\rm s}^{-1}$ should also host Pop III stars producing an observable $L_{\rm He II}>10^{42.7}{\rm erg}\,{\rm s}^{-1}$. As an alternate explanation, we explore the possibility that CR7 is instead powered by accretion onto a Direct Collapse Black Hole (DCBH). Our model predicts $L_{\alpha}$, $L_{\rm He II}$, and X-ray luminosities that are in agreement with the observations. In any case, the observed properties of CR7 indicate that this galaxy is most likely powered by sources formed from pristine gas. We propose that further X-ray observations can distinguish between the two above scenarios.Comment: 6 pages, 4 figure