Source author record

Daniel Schaerer

Daniel Schaerer appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

59works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

59 published item(s)

preprint2026arXiv

VENUS: Two Faint Little Red Dots Separated by $\sim70\,\mathrm{pc}$ Hidden in a Single Lensed Galaxy at $z\sim7$

We report the identification of a pair of faint little red dots (LRDs), dubbed Red Eyes, in a strongly-lensed galaxy at $z\sim7$ behind the PLCKG004.5-10.5 cluster, identified from the JWST Treasury program VENUS. Red Eyes are spatially resolved on the image plane with distinct colors, while the critical curve lies far north of Red Eyes, clearly requiring two different LRDs rather than a single LRD. Red Eyes is an extremely close pair of LRDs separated by $\sim70\,\mathrm{pc}$ in the source plane with a magnification of $μ\sim20$, which consistently explains another counter-image detected to the north-west. Red Eyes is hosted in a typical star-forming galaxy with $M_{\mathrm{UV,int}}\sim -19$, but its own UV emission is very faint ($M_{\mathrm{UV,int}} \gtrsim -16$). Moreover, Red Eyes does not reside at the galaxy center but lies at an offset position of approximately one effective radius $R_{\mathrm{e}}$ away from the galaxy center. If observed without lensing, Red Eyes would appear as a typical star-forming galaxy at $z\sim 7$ with $M_{\mathrm{UV}}\sim -19$, showing no apparent LRD signatures in either morphology or SED. These results suggest that multiple off-center LRDs, similar to Red Eyes, may be commonly hidden in a typical high-$z$ star-forming galaxy. In this case, various plausible scenarios may emerge, one of which is that intermediate-mass black holes (IMBHs) with $M_\mathrm{BH}\sim10^{4\text{--}6}\,M_\odot$ may form in star clusters on a stellar disk and contribute to the growth of the central supermassive black hole via mergers, with some IMBHs detectable as luminous LRDs in a sufficiently active and massive phase.

preprint2022arXiv

(Re)Solving Reionization with Lyα: How Bright Lyα Emitters account for the $z\approx2-8$ Cosmic Ionizing Background

The cosmic ionizing emissivity from star-forming galaxies has long been anchored to UV luminosity functions. Here we introduce an emissivity framework based on Ly$α$ emitters (LAEs), which naturally hones in on the subset of galaxies responsible for the ionizing background due to the intimate connection between the production and escape of Ly$α$ and LyC photons. Using constraints on the escape fractions of bright LAEs ($L_{\rm{Lyα}}>0.2 L^{*}$) at $z\approx2$ obtained from resolved Ly$α$ profiles, and arguing for their redshift-invariance, we show that: (i) quasars and LAEs together reproduce the relatively flat emissivity at $z\approx2-6$, which is non-trivial given the strong evolution in both the star-formation density and quasar number density at these epochs and (ii) LAEs produce late and rapid reionization between $z\approx6-9$ under plausible assumptions. Within this framework, the $>10\times$ rise in the UV population-averaged $f_{\rm{esc}}$ between $z\approx3-7$ naturally arises due to the same phenomena that drive the growing Ly$α$ emitter fraction with redshift. Generally, a LAE dominated emissivity yields a peak in the distribution of the ionizing budget with UV luminosity as reported in latest simulations. Using our adopted parameters ($f_{\rm{esc}}=50\%$, $ξ_{\rm{ion}}=10^{25.9}$ Hz erg$^{-1}$ for half the bright LAEs), a highly ionizing minority of galaxies with $M_{\rm UV}<-17$ accounts for the entire ionizing budget from star-forming galaxies. Rapid flashes of LyC from such rare galaxies produce a "disco" ionizing background. We conclude proposing tests to further develop our suggested Ly$α$-anchored formalism.

preprint2022arXiv

Conditions for detecting lensed Population III galaxies in blind surveys with the James Webb Space Telescope, the Roman Space Telescope and Euclid

Dark matter halos that reach the HI-cooling mass without prior star formation or external metal pollution represent potential sites for the formation of small - extremely faint - Population III galaxies at high redshifts. Gravitational lensing may in rare cases boost their fluxes to detectable levels, but to find even a small number of such objects in randomly selected regions of the sky requires very large areas to be surveyed. Because of this, a small, wide-field telescope can in principle offer better detection prospects than a large telescope with a smaller field of view. Here, we derive the minimum comoving number density required to allow gravitational lensing to lift such objects at redshift $z=5-16$ above the detection thresholds of blind surveys carried out with the James Webb space telescope (JWST), the Roman space telescope (RST) and Euclid. We find that the prospects for photometric detections of Pop III galaxies is promising, and that they are better for RST than for JWST and Euclid. However, the Pop III galaxies favoured by current simulations have number densities too low to allow spectroscopic detections based on the strength of the HeII1640 emission line in any of the considered surveys unless very high star formation efficiencies ($ε> 0.1$) are envoked. We argue that targeting individual cluster lenses instead of the wide field surveys considered in this paper results in better spectroscopic detection prospects, while for photometric detection, the wide field surveys perform considerably better.

preprint2022arXiv

Ionization of HeII in star-forming galaxies by X-rays from cluster winds and superbubbles

The nature of the sources powering nebular HeII emission in star-forming galaxies remains debated, and various types of objects have been considered, including Wolf-Rayet stars, X-ray binaries, and Population III stars. Modern X-ray observations show the ubiquitous presence of hot gas filling star-forming galaxies. We use a collisional ionization plasma code to compute the specific HeII ionizing flux produced by hot gas and show that if its temperature is not too high (less than 2.5 MK), then the observed levels of soft diffuse X-ray radiation could explain HeII ionization in galaxies. To gain a physical understanding of this result, we propose a model that combines the hydrodynamics of cluster winds and hot superbubbles with observed populations of young massive clusters in galaxies. We find that in low-metallicity galaxies, the temperature of hot gas is lower and the production rate of HeII ionizing photons is higher compared to high-metallicity galaxies. The reason is that the slower stellar winds of massive stars in lower-metallicity galaxies input less mechanical energy in the ambient medium. Furthermore, we show that ensembles of star clusters up to 10-20 Myr old in galaxies can produce enough soft X-rays to induce nebular HeII emission. We discuss observations of the template low-metallicity galaxy I Zw 18 and suggest that the HeII nebula in this galaxy is powered by a hot superbubble. Finally, appreciating the complex nature of stellar feedback, we suggest that soft X-rays from hot superbubbles are among the dominant sources of HeII ionizing flux in low-metallicity star-forming galaxies.

preprint2022arXiv

The ALPINE-ALMA [CII] survey: dust attenuation curves at z=4.4-5.5

There is now ample evidence that dust is already present in abundance at high z. However, given the faintness of distant galaxies in the optical and the NIR, datasets are still limited and how the dust affects the emerging radiation of galaxies at very high redshift is not yet fully understood. Using the ALPINE survey, our objective is to quantify the dust attenuation properties in galaxies at z=4.4-5.5, and in particular the shape of their attenuation curve. Using the CIGALE code, we model the stellar populations and their interaction with the dust in order to measure some of the physical properties of a subsample of 23 main-sequence ALPINE galaxies. We find that the attenuation curves span a broad range of properties, from curves that are much steeper than the SMC extinction curve, to shallower than the starburst attenuation curve. The shape of the attenuation curves strongly depends on the V-band attenuation. Galaxies with the lowest attenuation also present the steepest curves. The steepness of such curves is probably the consequence of the combination of the intrinsic physical properties of the dust, the relative distribution of stars and dust in the interstellar medium, and the differential reddening. The broad range of attenuation curves found at z~5 shows that no single attenuation curve is appropriate for main sequence galaxies and that assuming a fixed curve can lead to large errors, for instance in the interpretation and use of the IRX-beta diagram, if SED modeling is not feasible. Great caution should be exercised when correcting high redshift galaxies for the presence of dust using the UV slope beta as it can affect the estimation of both SFR and stellar mass even at low V-band attenuation due to the steepness of the attenuation curve. However, when SED modeling can be used, the impact of the choice of the attenuation curve on the SFR and the stellar mass is limited.

preprint2022arXiv

The ALPINE-ALMA [CII] survey: The infrared-radio correlation and AGN fraction of star-forming galaxies at z $\sim$ 4.4-5.9

We present the radio properties of 66 spectroscopically-confirmed normal star-forming galaxies (SFGs) at $4.4<z<5.9$ in the COSMOS field that were [C II] detected in the Atacama Large Millimeter Array (ALMA) Large Program to INvestigate [C II] at Early times (ALPINE). We separate these galaxies ("CII-detected-all") into lower redshift ("CII-detected-lz", $\langle z\rangle=4.5$) and higher redshift ("CII-detected-hz", $\langle z\rangle=5.6$) sub-samples and stack multi-wavelength imaging for each sub-sample from X-ray to radio bands. A radio signal is detected in the stacked 3 GHz image of CII-detected-all and -lz samples at $\gtrsim3σ$. We find that the infrared-radio correlation of our sample, quantified by $q_{\mathrm{TIR}}$, is lower than the local relation for normal SFGs at $\sim$3$σ$ significance level, and is instead broadly consistent with that of bright sub-mm galaxies at $2<z<5$. Neither of these samples show evidence of dominant AGN activity in their stacked Spectral Energy Distributions (SEDs), rest-frame UV spectra, or X-ray images. Although we cannot rule out the possible effect of the assumed spectral index and the applied infrared SED templates as at least partially causing these differences, the lower obscured fraction of star formation than at lower redshift can alleviate the tension between our stacked $q_{\mathrm{TIR}}$s and that of local normal SFGs. It is possible that the dust buildup, which primarily governs the IR emission in addition to older stellar populations, has not had enough time to occur fully in these galaxies, whereas the radio emission can respond on a more rapid timescale. Therefore, we might expect a lower $q_{\mathrm{TIR}}$ to be a general property of high-redshift SFGs.

preprint2022arXiv

The Low-Redshift Lyman Continuum Survey I: New, Diverse Local Lyman-Continuum Emitters

The origins of Lyman continuum (LyC) photons responsible for the reionization of the universe are as of yet unknown and highly contested. Detecting LyC photons from the epoch of reionization is not possible due to absorption by the intergalactic medium, which has prompted the development of several indirect diagnostics to infer the rate at which galaxies contribute LyC photons to reionize the universe by studying lower-redshift analogs. We present the Low-redshift Lyman Continuum Survey (LzLCS) comprising measurements made with HST/COS for a z=0.2-0.4 sample of 66 galaxies. After careful processing of the FUV spectra, we obtain a total of 35 Lyman continuum emitters (LCEs) detected with 97.725% confidence, nearly tripling the number of known local LCEs. We estimate escape fractions from the detected LyC flux and upper limits on the undetected LyC flux, finding a range of LyC escape fractions up to 50%. Of the 35 LzLCS LCEs, 12 have LyC escape fractions greater than 5%, more than doubling the number of known local LCEs with cosmologically relevant LyC escape.

preprint2022arXiv

The Low-Redshift Lyman Continuum Survey II: New Insights into LyC Diagnostics

The Lyman continuum (LyC) cannot be observed at the epoch of reionization (z {\gtrsim} 6) due to intergalactic H I absorption. To identify Lyman continuum emitters (LCEs) and infer the fraction of escaping LyC, astronomers have developed various indirect diagnostics of LyC escape. Using measurements of the LyC from the Low-redshift Lyman Continuum Survey (LzLCS), we present the first statistical test of these diagnostics. While optical depth indicators based on Lyα, such as peak velocity separation and equivalent width, perform well, we also find that other diagnostics, such as the [O III]/[O II] flux ratio and star formation rate surface density, predict whether a galaxy is a LCE. The relationship between these galaxy properties and the fraction of escaping LyC flux suggests that LyC escape depends strongly on H I column density, ionization parameter, and stellar feedback. We find LCEs occupy a range of stellar masses, metallicities, star formation histories, and ionization parameters, which may indicate episodic and/or different physical causes of LyC escape.

preprint2022arXiv

The Low-Redshift Lyman Continuum Survey. Unveiling the ISM properties of low-$z$ Lyman continuum emitters

Combining 66 ultraviolet (UV) spectra and ancillary data from the Low-Redshift Lyman Continuum Survey (LzLCS) and 23 LyC observations by earlier studies, we form a statistical sample of star-forming galaxies at $z \sim 0.3$ to study the role of the cold interstellar medium (ISM) gas in the leakage of ionizing radiation. We first constrain the massive star content (ages and metallicities) and UV attenuation, by fitting the stellar continuum with a combination of simple stellar population models. The models, together with accurate LyC flux measurements, allow to determine the absolute LyC photon escape fraction for each galaxy ($f_{\rm esc}^{\rm abs}$). We measure the equivalent widths and residual fluxes of multiple HI and low-ionization state (LIS) lines, and the geometrical covering fraction adopting the picket-fence model. The $f_{\rm esc}^{\rm abs}$ spans a wide range, with a median (0.16, 0.84 quantiles) of 0.04 (0.02, 0.20), and 50 out of the 89 galaxies detected in the LyC. The HI and LIS line equivalent widths scale with the UV luminosity and attenuation, and inversely with the residual flux of the lines. The HI and LIS residual fluxes are correlated, indicating that the neutral gas is spatially traced by the LIS transitions. We find the observed trends of the absorption lines and the UV attenuation are primarily driven by the covering fraction. The non-uniform gas coverage demonstrates that LyC photons escape through low-column density channels in the ISM. The equivalent widths and residual fluxes of the UV lines strongly correlate with $f_{\rm esc}^{\rm abs}$: strong LyC leakers show weak absorption lines, low UV attenuation, and large Ly$α$ equivalent widths. We finally show that simultaneous UV absorption line and dust attenuation measurements can predict, on average, the escape fraction of galaxies and the method can be applied to galaxies across a wide redshift range.

preprint2022arXiv

Tracing Lyman-alpha and Lyman Continuum Escape in Galaxies with Mg II Emission

Star-forming galaxies are considered the likeliest source of the H I ionizing Lyman Continuum (LyC) photons that reionized the intergalactic medium at high redshifts. However, above z >~ 6, the neutral intergalactic medium prevents direct observations of LyC. Therefore, recent years have seen the development of indirect indicators for LyC that can be calibrated at lower redshifts and applied in the Epoch of Reionization. Emission from Mg II \ly\ly 2796, 2803 doublet has been proposed as a promising LyC proxy. In this paper, we present new Hubble Space Telescope/Cosmic Origins Spectrograph observations for 8 LyC emitter candidates, selected to have strong Mg II emission lines. We securely detect LyC emission in 50% (4/8) galaxies with 2$σ$ significance. This high detection rate suggests that strong Mg II emitters might be more likely to leak LyC than similar galaxies without strong Mg II. Using photoionization models, we constrain the escape fraction of Mg II as ~ 15 -- 60%. We confirm that the escape fraction of Mg II correlates tightly with that of Lyman-alpha (LyA), which we interpret as an indication that the escape fraction of both species is controlled by resonant scattering in the same low column density gas. Furthermore, we show that the combination of the Mg II emission and dust attenuation can be used to estimate the escape fraction of LyC statistically. These findings confirm that Mg II emission can be adopted to estimate the escape fraction of LyA and LyC in local star-forming galaxies and may serve as a useful indirect indicator at the Epoch of Reionization.

preprint2021arXiv

ALMA 1.3 mm Survey of Lensed Submillimeter Galaxies (SMGs) Selected by Herschel: Discovery of Spatially Extended SMGs and Implications

We present an ALMA 1.3 mm (Band 6) continuum survey of lensed submillimeter galaxies (SMGs) at $z=1.0\sim3.2$ with an angular resolution of $\sim0.2$". These galaxies were uncovered by the Herschel Lensing Survey (HLS), and feature exceptionally bright far-infrared continuum emission ($S_\mathrm{peak} \gtrsim 90$ mJy) owing to their lensing magnification. We detect 29 sources in 20 fields of massive galaxy clusters with ALMA. Using both the Spitzer/IRAC (3.6/4.5 $\mathrm{μm}$) and ALMA data, we have successfully modeled the surface brightness profiles of 26 sources in the rest-frame near- and far-infrared. Similar to previous studies, we find the median dust-to-stellar continuum size ratio to be small ($R_\mathrm{e,dust}/R_\mathrm{e,star} = 0.38\pm0.14$) for the observed SMGs, indicating that star formation is centrally concentrated. This is, however, not the case for two spatially extended main-sequence SMGs with a low surface brightness at 1.3 mm ($\lesssim 0.1$ mJy arcsec$^{-2}$), in which the star formation is distributed over the entire galaxy ($R_\mathrm{e,dust}/R_\mathrm{e,star}>1$). As a whole, our SMG sample shows a tight anti-correlation between ($R_\mathrm{e,dust}/R_\mathrm{e,star}$) and far-infrared surface brightness ($Σ_\mathrm{IR}$) over a factor of $\simeq$ 1000 in $Σ_\mathrm{IR}$. This indicates that SMGs with less vigorous star formation (i.e., lower $Σ_\mathrm{IR}$) lack central starburst and are likely to retain a broader spatial distribution of star formation over the whole galaxies (i.e., larger $R_\mathrm{e,dust}/R_\mathrm{e,star}$). The same trend can be reproduced with cosmological simulations as a result of central starburst and potentially subsequent "inside-out" quenching, which likely accounts for the emergence of compact quiescent galaxies at $z\sim2$.

preprint2021arXiv

Radial profiles of lensed $z \sim 1$ galaxies on sub-kiloparsec scales

We study the spatially resolved physical properties of the Cosmic Snake arc in MACS J1206.2-0847 and the arc in Abell 0521 (A521). These are two strongly lensed galaxies at redshifts $z=1.036$ and $z=1.044$. We used observations of the Hubble Space Telescope (HST) and the Atacama Large Millimeter/submillimeter Array (ALMA). The former gives access to the star formation rate (SFR) and stellar mass ($M_\star$), and the latter to the H$_2$ molecular gas mass ($M_{\mathrm{mol}}$). HST and ALMA observations have similar angular resolutions of $0.15^{\prime \prime}-0.2^{\prime \prime}$, which with the help of strong gravitational lensing enable us to reach spatial resolutions down to $\sim 30\,\mathrm{pc}$ and $\sim 100\,\mathrm{pc}$ in these two galaxies, respectively. These resolutions are close to the resolution of observations of nearby galaxies. We study the radial profiles of SFR, $M_\star$, and $M_{\mathrm{mol}}$ surface densities of these high-redshift galaxies and compare the corresponding exponential scale lengths with those of local galaxies. We find that the scale lengths in the Cosmic Snake are about $0.5\,\mathrm{kpc}-1.5\,\mathrm{kpc}$, and they are 3 to 10 times larger in A521. This is a significant difference knowing that the two galaxies have comparable integrated properties. These high-redshift scale lengths are nevertheless comparable to those of local galaxies, which cover a wide distribution. The particularity of our high-redshift radial profiles is the normalisation of the $M_{\mathrm{mol}}$ surface density profiles ($ΣM_{\mathrm{mol}}$), which are offset by up to a factor of 20 with respect to the profiles of $z=0$ counterparts. The SFR surface density profiles are also offset by the same factor as $ΣM_{\mathrm{mol}}$, as expected from the Kennicutt-Schmidt law.

preprint2021arXiv

SDSS-IV MaNGA: an observational evidence of density bounded region in a Lyman-$α$ emitter

Using Integral Field Unit (IFU) spectroscopy, we present here the spatially resolved morphologies of [SII]$λ$6717,6731/H$α$ and [SII]$λ$6717,6731/[OIII]$λ$5007 emission line ratios for the first time in a Blueberry Lyman-$α$ emitter (BBLAE) at z $\sim$ 0.047. Our derived morphologies show that the extreme starburst region of the BBLAE, populated by young ($\leqslant$ 10 Myr), massive Wolf-Rayet stars, is [SII]-deficient, while the rest of the galaxy is [SII]-enhanced. We infer that the extreme starburst region is density-bounded (i.e., optically thin to ionizing photons), and the rest of the galaxy is ionization-bounded $-$ indicating a Blister-type morphology. We find that the previously reported small escape fraction (10%) of Ly$α$ photons is from our identified density-bounded HII region of the BBLAE. This escape fraction is likely constrained by a porous dust distribution. We further report a moderate correlation between [SII]-deficiency and inferred Lyman Continuum (LyC) escape fraction using a sample of confirmed LyC leakers studied in the literature, including the BBLAE studied here. The observed correlation also reveals its dependency on the stellar mass and gas-phase metallicity of the leaky galaxies. Finally, the future scope and implications of our work are discussed in detail.

preprint2020arXiv

AstroSat detection of Lyman continuum emission from a z=1.42 galaxy

One of the outstanding problems of current observational cosmology is to understand the nature of sources that produced the bulk of the ionizing radiation after the Cosmic Dark Age. Direct detection of these reionization sources is practically infeasible at high redshift due to the steep decline of intergalactic medium transmission. However, a number of low-redshift analogs emitting Lyman continuum at~900 Angstrom ~restframe are now detected at $z< 0.4$ and there are detections in the range $ 2.5< z< 3.5$ also. Here, we report the detection of Lyman continuum emission with a high escape fraction (>20%) from a low-mass clumpy galaxy at z=1.42, in the middle of the redshift range where no detection has been made before and near the peak of the Cosmic Star-formation history. The observation was made in the Hubble Extreme Deep field by the wide-field Ultra-Violet-Imaging Telescope onboard AstroSat. This is the first detection of Extreme Ultraviolet radiation from a distant galaxy at a rest-frame wavelength of 600 Angstrom, and it opens up a new window to constrain the shape of the ionization spectrum. Further observations with AstroSat should significantly increase the sample of Lyman continuum leaking galaxies at Cosmic Noon.

preprint2020arXiv

Luminous Lyman-alpha Emitters with Very Blue UV-continuum Slopes at Redshift 5.7 <= z <= 6.6

We study six luminous Lyman-alpha emitters (LAEs) with very blue rest-frame UV continua at $5.7\le z \le 6.6$. These LAEs have previous HST and Spitzer IRAC observations. Combining our newly acquired HST images, we find that their UV-continuum slopes $β$ are in a range of $-3.4\le β\le -2.6$. Unlike previous, tentative detections of $β\simeq -3$ in photometrically selected, low-luminosity galaxies, our LAEs are spectroscopically confirmed and luminous ($M_{\rm UV}<-20$ mag). We model their broadband spectral energy distributions (SEDs), and find that two $β\simeq-2.6\pm0.2$ galaxies can be well fitted with young and dust-free stellar populations. However, it becomes increasingly difficult to fit bluer galaxies. We explore further interpretations by including non-zero LyC escape fraction $f_{\rm esc}$, very low metallicities, and/or AGN contributions. Assuming $f_{\rm esc}\simeq0.2$, we achieve the bluest slopes $β\simeq-2.7$ when nebular emission is considered. This can nearly explain the SEDs of two galaxies with $β\simeq-2.8$ and --2.9 ($σ_β=0.15$). Larger $f_{\rm esc}$ values and very low metallicities are not favored by the strong nebular line emission (evidenced by the IRAC flux) or the observed (IRAC 1 - IRAC 2) color. Finally, we find that the $β\simeq-2.9$ galaxy can potentially be well explained by the combination of a very young population with a high $f_{\rm esc}$ ($\ge0.5$) and an old, dusty population. We are not able to produce two $β\simeq -3.4 \pm0.4$ galaxies. Future deep spectroscopic observations are needed to fully understand these galaxies.

preprint2020arXiv

The ALPINE-ALMA [C II] survey: the luminosity function of serendipitous [C II] line emitters at $z\sim 5$

We present the first [CII] 158 $μ$m luminosity function (LF) at $z\sim 5$ from a sample of serendipitous lines detected in the ALMA Large Program to INvestigate [CII] at Early times (ALPINE). A search performed over the 118 ALPINE pointings revealed several serendipitous lines. Based on their fidelity, we selected 14 lines for the final catalog. According to the redshift of their counterparts, we identified 8 out of 14 detections as [CII] lines at $z\sim 5$, and two as CO transitions at lower redshifts. The remaining 4 lines have an elusive identification in the available catalogs and we considered them as [CII] candidates. We used the 8 confirmed [CII] and the 4 [CII] candidates to build one of the first [CII] LFs at $z\sim 5$. We found that 11 out of these 12 sources have a redshift very similar to that of the ALPINE target in the same pointing, suggesting the presence of overdensities around the targets. Therefore, we split the sample in two (a "clustered" and "field" sub-sample) according to their redshift separation and built two separate LFs. Our estimates suggest that there could be an evolution of the [CII] LF between $z \sim 5$ and $z \sim 0$. By converting the [CII] luminosity to star formation rate we evaluated the cosmic star formation rate density (SFRD) at $z\sim 5$. The clustered sample results in a SFRD $\sim 10$ times higher than previous measurements from UV-selected galaxies. On the other hand, from the field sample (likely representing the average galaxy population) we derived a SFRD $\sim 1.6$ higher compared to current estimates from UV surveys but compatible within the errors. Because of the large uncertainties, observations of larger samples are necessary to better constrain the SFRD at $z\sim 5$. This study represents one of the first efforts aimed at characterizing the demography of [CII] emitters at $z\sim 5$ using a mm-selection of galaxies.

preprint2019arXiv

HST Imaging of the Ionizing Radiation from a Star-forming Galaxy at z = 3.794

We report on the HST detection of the Lyman-continuum (LyC) radiation emitted by a galaxy at redshift z=3.794, dubbed Ion1 (Vanzella et al. 2012). The LyC from Ion1 is detected at rest-frame wavelength 820$\sim$890 Åwith HST WFC3/UVIS in the F410M band ($m_{410}=27.60\pm0.36$ magnitude (AB), peak SNR = 4.17 in a circular aperture with radius r = 0.12'') and at 700$\sim$830 Åwith the VLT/VIMOS in the U-band ($m_U = 27.84\pm0.19$ magnitude (AB), peak SNR = 6.7 with a r = 0.6'' aperture). A 20-hr VLT/VIMOS spectrum shows low- and high-ionization interstellar metal absorption lines, the P-Cygni profile of CIV and Ly$α$ in absorption. The latter spectral feature differs from what observed in known LyC emitters, which show strong Ly$α$ emission. An HST far-UV color map reveals that the LyC emission escapes from a region of the galaxy that is bluer than the rest, presumably because of lower dust obscuration. The F410M image shows that the centroid of the LyC emission is offset from the centroid of the non-ionizing UV emission by 0.12''$\pm$0.03'', corresponding to 0.85$\pm$0.21 kpc (physical), and that its morphology is likely moderately resolved. These morphological characteristics favor a scenario where the LyC photons produced by massive stars escape from low HI column-density "cavities" in the ISM, possibly carved by stellar winds and/or supernova. We also collect the VIMOS U-band images of a sample of 107 Lyman-break galaxies with spectroscopic redshifts at $3.40<z<3.95$, i.e. sampling the LyC, and stack them with inverse-variance weights. No LyC emission is detected in the stacked image, resulting in a 32.5 magnitude (AB) flux limit (1$σ$) and an upper limit of absolute LyC escape fraction $f_{esc}^{abs} < 0.63\%$. LyC emitters like Ion1 are very likely at the bright-end of the LyC luminosity function.

preprint2016arXiv

A Close Comparison between Observed and Modeled Lyα Lines for z ~ 2.2 Lyman Alpha Emitters

We present the results of a Lya profile analysis of 12 Lya emitters (LAEs) at z = 2.2 with high-resolution Lya spectra. We find that all 12 objects have a Lya profile with the main peak redward of the systemic redshift defined by nebular lines, and five have a weak, secondary peak blueward of the systemic redshift (blue bump). The average velocity offset of the red main peak (the blue bump, if any) with respect to the systemic redshift is Delta_v_Lya,r = 174+/- 19 km s-1 (Delta_v_Lya,b = -316+/-45 km s-1), which is smaller than (comparable to) that of Lyman-break galaxies (LBGs). The outflow velocities inferred from metal absorption lines in three individual and one stacked spectra are comparable to those of LBGs. The expanding shell model constructed by Verhamme et al. (2006) reproduces not only the Lya profiles but also other observed quantities including the outflow velocity and the FWHM of nebular lines for the non-blue bump objects. On the other hand, the model predicts too high FWHMs of nebular lines for the blue bump objects, although this discrepancy may disappear if we introduce additional Lya photons produced by gravitational cooling. We show that the small Delta_v_Lya,r values of our sample can be explained by low neutral-hydrogen column densities of log(NHI) = 18.9 cm-2 on average. This value is more than one order of magnitude lower than those of LBGs but is consistent with recent findings that LAEs have high ionization parameters and low Hi gas masses. This result suggests that low NHI values, giving reduced numbers of resonant scattering of Lya photons, are the key to the strong Lya emission of LAEs.

preprint2016arXiv

Ly$α$ Emitters with Very Large Ly$α$ Equivalent Widths, EW$_{\rm 0}$(Ly$α$) $\simeq 200-400$ Å, at $z\sim 2$

We present physical properties of spectroscopically confirmed Ly$α$ emitters (LAEs) with very large rest-frame Ly$α$ equivalent widths EW$_{\rm 0}$(Ly$α$). Although the definition of large EW$_{\rm 0}$(Ly$α$) LAEs is usually difficult due to limited statistical and systematic uncertainties, we identify six LAEs selected from $\sim 3000$ LAEs at $z\sim 2$ with reliable measurements of EW$_{\rm 0}$ (Ly$α$) $\simeq 200-400$ Å given by careful continuum determinations with our deep photometric and spectroscopic data. These large EW$_{\rm 0}$(Ly$α$) LAEs do not have signatures of AGN, but notably small stellar masses of $M_{\rm *} = 10^{7-8}$ $M_{\rm \odot}$ and high specific star-formation rates (star formation rate per unit galaxy stellar mass) of $\sim 100$ Gyr$^{-1}$. These LAEs are characterized by the median values of $L({\rm Lyα})=3.7\times 10^{42}$ erg s$^{-1}$ and $M_{\rm UV}=-18.0$ as well as the blue UV continuum slope of $β= -2.5\pm0.2$ and the low dust extinction $E(B-V)_{\rm *} = 0.02^{+0.04}_{-0.02}$, which indicate a high median Ly$α$ escape fraction of $f_{\rm esc}^{\rm Lyα}=0.68\pm0.30$. This large $f_{\rm esc}^{\rm Lyα}$ value is explained by the low {\sc Hi} column density in the ISM that is consistent with FWHM of the Ly$α$ line, ${\rm FWHM (Lyα)}=212\pm32$ km s$^{-1}$, significantly narrower than those of small EW$_{\rm 0}$(Ly$α$) LAEs. Based on the stellar evolution models, our observational constraints of the large EW$_{\rm 0}$ (Ly$α$), the small $β$, and the rest-frame He{\sc ii} equivalent width imply that at least a half of our large EW$_{\rm 0}$(Ly$α$) LAEs would have young stellar ages of $\lesssim 20$ Myr and very low metallicities of $Z<0.02 Z_\odot$ regardless of the star-formation history.

preprint2015arXiv

Are Ultra-faint Galaxies at z=6-8 Responsible for Cosmic Reionization ? Combined Constraints from the Hubble Frontier Fields Clusters and Parallels

We use deep Hubble Space Telescope imaging of the Frontier Fields to accurately measure the galaxy rest-frame ultraviolet luminosity function (UV LF) in the redshift range $z \sim 6-8$. We combine observations in three lensing clusters A2744, MACS0416, MACS0717 and their associated parallels fields to select high-redshift dropout candidates. We use the latest lensing models to estimate the flux magnification and the effective survey volume in combination with completeness simulations performed in the source plane. We report the detection of 227 galaxy candidates at $z=6-7$ and 25 candidates at $z \sim 8$. While the total survey area is about 4 arcmin$^{2}$ in each parallel field, it drops to about 0.6 to 1 arcmin$^{2}$ in the cluster core fields because of the strong lensing. We compute the UV luminosity function at $z \sim 7$ using the combined galaxy sample and perform Monte Carlo simulations to determine the best fit Schechter parameters. We are able to reliably constrain the LF down to an absolute magnitude of $M_{UV}=-15.25$, which corresponds to 0.005$L^{\star}$. More importantly, we find that the faint-end slope remains steep down to this magnitude limit with $α=-2.04_{-0.17}^{+0.13}$. Our results confirm the most recent results in deep blank fields but extend the LF measurements more than two magnitudes deeper. The UV LF at $z \sim 8$ is not very well constrained below $M_{UV}=-18$ due to the small number statistics and incompleteness uncertainties. To assess the contribution of galaxies to cosmic reionization we derive the UV luminosity density at $z\sim7$ by integrating the UV LF down to an observationally constrained limit of $M_{UV} = -15$. We show that our determination of Log($ρ_{UV}$)=$26.2\pm0.13$ (erg s$^{-1}$ Hz$^{-1}$ Mpc$^{-3}$) can be sufficient to maintain the IGM ionized.

preprint2015arXiv

Evidence for PopIII-like stellar populations in the most luminous Lyman-$α$ emitters at the epoch of re-ionisation: spectroscopic confirmation

Faint Lyman-$α$ (Ly$α$) emitters become increasingly rarer towards the re-ionisation epoch (z~6-7). However, observations from a very large (~5deg$^2$) Ly$α$ survey at z=6.6 (Matthee et al. 2015) show that this is not the case for the most luminous emitters. Here we present follow-up observations of the two most luminous z~6.6 Ly$α$ candidates in the COSMOS field: `MASOSA' and `CR7'. We used X-SHOOTER, SINFONI and FORS2 (VLT), and DEIMOS (Keck), to confirm both candidates beyond any doubt. We find redshifts of z=6.541 and z=6.604 for MASOSA and CR7, respectively. MASOSA has a strong detection in Ly$α$ with a line width of $386\pm30$ km/s (FWHM) and with high EW$_0$ (>200 Å), but it is undetected in the continuum. CR7, with an observed Ly$α$ luminosity of $10^{43.93\pm0.05}$erg/s is the most luminous Ly$α$ emitter ever found at z>6. CR7 reveals a narrow Ly$α$ line with $266\pm15$ km/s FWHM, being detected in the NIR (rest-frame UV, with $β=-2.3\pm0.1$) with an excess in $J$, and also strongly detected in IRAC/Spitzer. We detect a narrow HeII1640$Å$ emission line ($6σ$) which explains the excess seen in the $J$ band photometry (EW$_0$~80 Å). We find no other emission lines from the UV to the NIR in our X-SHOOTER spectra (HeII/OIII]1663A>3 and HeII/CIII]1908A>2.5). We find that CR7 is best explained by a combination of a PopIII-like population which dominates the rest-frame UV and the nebular emission, and a more normal stellar population which dominates the mass. HST/WFC3 observations show that the light is indeed spatially separated between a very blue component, coincident with Ly$α$ and HeII emission, and two red components (~5 kpc away), which dominate the mass. Our findings are consistent with theoretical predictions of a PopIII wave, with PopIII star formation migrating away from the original sites of star formation.

preprint2015arXiv

The Lyman alpha Reference Sample VI: Lyman alpha escape from the edge-on disk galaxy Mrk1486

While numerical simulations suggest that the strength of the Lyman alpha (Lya) line of star-forming disk galaxies strongly depends on the inclination at which they are observed (i.e. from edge-on to face-on, we expect to see a change from an attenuated Lya line to a strong Lya emission line), recent observations with the Hubble space telescope (HST) have highlighted few low-redshift highly inclined (edge-on) disk galaxies that breaks this trend. We aim to understand how a strong Lya emission line is able to escape from one of those inclined disk galaxies, named Mrk1486 (z=0.0338). For that purpose we used a large set of HST imaging and spectroscopic data to investigate both the ISM structure and the dominant source of Lya radiation inside Mrk1486. Moreover, we used a 3D Monte Carlo Lya radiation transfer code to study the radiative transfer of Lya and UV continuum photons inside a 3D geometry of neutral hydrogen (HI) and dust that models the ISM structure at the galaxy center. The analysis of IFU Halpha spectroscopic data of Mrk1486 indicates the presence of two bipolar galactic winds of HI gas above and bellow the disk plane of Mrk1486. Furthermore, comparing different diagnostic diagrams (such as [OIII]5007/Hbeta versus [OI]6300/Halpha) to photo- and shock-ionization models, we find that the Lya production of Mrk1486 is dominated by photoionization inside the galaxy disk. From this perspective, our numerical simulations succeed in reproducing the strength and spectral shape of the observed Lya line of Mrk1486 by assuming a scenario in which the Lya photons are produced inside the disk, travel along the galactic winds and scatter on cool HI materials toward the observer. As bipolar galactic winds are ubiquitous in star-forming disk galaxies, this mechanism may explain the origin of strong Lya emission lines recently observed from highly inclined galaxies at high-redshift.

preprint2015arXiv

The Lyman Alpha Reference Sample: V. The impact of neutral ISM kinematics and geometry on Lyman Alpha escape

We present high-resolution far-UV spectroscopy of the 14 galaxies of the Lyman Alpha Reference Sample; a sample of strongly star-forming galaxies at low redshifts ($0.028 < z < 0.18$). We compare the derived properties to global properties derived from multi band imaging and 21 cm HI interferometry and single dish observations, as well as archival optical SDSS spectra. Besides the Lyman $α$ line, the spectra contain a number of metal absorption features allowing us to probe the kinematics of the neutral ISM and evaluate the optical depth and and covering fraction of the neutral medium as a function of line-of-sight velocity. Furthermore, we show how this, in combination with precise determination of systemic velocity and good Ly$α$ spectra, can be used to distinguish a model in which separate clumps together fully cover the background source, from the "picket fence" model named by Heckman et al. (2011). We find that no one single effect dominates in governing Ly$α$ radiative transfer and escape. Ly$α$ escape in our sample coincides with a maximum velocity-binned covering fraction of $\lesssim 0.9$ and bulk outflow velocities of $\gtrsim 50$ km s$^{-1}$, although a number of galaxies show these characteristics and yet little or no Ly$α$ escape. We find that Ly$α$ peak velocities, where available, are not consistent with a strong backscattered component, but rather with a simpler model of an intrinsic emission line overlaid by a blueshifted absorption profile from the outflowing wind. Finally, we find a strong anticorrelation between H$α$ equivalent width and maximum velocity-binned covering factor, and propose a heuristic explanatory model.

preprint2015arXiv

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

preprint2015arXiv

The Science Case for Multi-Object Spectroscopy on the European ELT

This White Paper presents the scientific motivations for a multi-object spectrograph (MOS) on the European Extremely Large Telescope (E-ELT). The MOS case draws on all fields of contemporary astronomy, from extra-solar planets, to the study of the halo of the Milky Way and its satellites, and from resolved stellar populations in nearby galaxies out to observations of the earliest 'first-light' structures in the partially-reionised Universe. The material presented here results from thorough discussions within the community over the past four years, building on the past competitive studies to agree a common strategy toward realising a MOS capability on the E-ELT. The cases have been distilled to a set of common requirements which will be used to define the MOSAIC instrument, entailing two observational modes ('high multiplex' and 'high definition'). When combined with the unprecedented sensitivity of the E-ELT, MOSAIC will be the world's leading MOS facility. In analysing the requirements we also identify a high-multiplex MOS for the longer-term plans for the E-ELT, with an even greater multiplex (>1000 targets) to enable studies of large-scale structures in the high-redshift Universe. Following the green light for the construction of the E-ELT the MOS community, structured through the MOSAIC consortium, is eager to realise a MOS on the E-ELT as soon as possible. We argue that several of the most compelling cases for ELT science, in highly competitive areas of modern astronomy, demand such a capability. For example, MOS observations in the early stages of E-ELT operations will be essential for follow-up of sources identified by the James Webb Space Telescope (JWST). In particular, multi-object adaptive optics and accurate sky subtraction with fibres have both recently been demonstrated on sky, making fast-track development of MOSAIC feasible.

preprint2014arXiv

New Constraints on the Faint-end of the UV Luminosity Function at z~7-8 using the Gravitational Lensing of the Hubble Frontier Fields Cluster A2744

Exploiting the power of gravitational lensing, the Hubble Frontier Fields (HFF) program aims at observing six massive galaxy clusters to explore the distant Universe far beyond the depth limits of blank field surveys. Using the complete Hubble Space Telescope observations of the first HFF cluster Abell 2744, we report the detection of 50 galaxy candidates at $z \sim 7$ and eight candidates at $z \sim 8$ in a total survey area of 0.96 arcmin$^{2}$ in the source plane. Three of these galaxies are multiply-imaged by the lensing cluster. Using an updated model of the mass distribution in the cluster we were able to calculate the magnification factor and the effective survey volume for each galaxy in order to compute the ultraviolet galaxy luminosity function at both redshifts 7 and 8. Our new measurements extend the $z \sim 7$ UV LF down to an absolute magnitude of $M_{UV} \sim -15.5$. We find a characteristic magnitude of $M^{\star}_{UV}=-20.63^{+0.69}_{-0.56}$ mag and a faint-end slope $α= -1.88^{+0.17}_{-0.20}$ close to previous determinations in blank fields. We show here for the first time that this slope remains steep down to very faint luminosities of 0.01$L^{\star}$. Although prone to large uncertainties, our results at $z \sim 8$ also seem to confirm a steep faint-end slope below 0.1$L^{\star}$. The HFF program is therefore providing an extremely efficient way to study the faintest galaxy populations at $z > 7$ that would otherwise be inaccessible with current instrumentation. The full sample of six galaxy clusters will provide yet better constraints on the build-up of galaxies at early epochs and their contribution to cosmic reionization.

preprint2014arXiv

The Effects of Stellar Rotation. II. A Comprehensive Set of Starburst99 Models

We present a new set of synthesis models for stellar populations obtained with Starburst99, which are based on new stellar evolutionary tracks with rotation. We discuss models with zero rotation velocity and with velocities of 40% of the break-up velocity on the zero-age main-sequence. These values are expected to bracket realistic rotation velocity distributions in stellar populations. The new rotating models for massive stars are more luminous and hotter due to a larger convective core and enhanced surface abundances. This results in pronounced changes in the integrated spectral energy distribution of a population containing massive stars. The changes are most significant at the shortest wavelengths where an increase of the ionizing luminosity by up to a factor of 5 is predicted. We also show that high equivalent widths of recombination lines may not necessarily indicate a very young age but can be achieved at ages as late as 10 Myr. Comparison of these two boundary cases (0 and 40% of the break-up velocity) will allow users to evaluate the effects of rotation and provide guidance for calibrating the stellar evolution models. We also introduce a new theoretical ultraviolet spectral library built from the Potsdam Wolf-Rayet (PoWR) atmospheres. Its purpose is to help identify signatures of Wolf-Rayet stars in the ultraviolet whose strength is sensitive to the particulars of the evolution models. The new models are available for solar and 1/7th solar metallicities. A complete suite of models can be generated on the Starburst99 website (www.stsci.edu/science/starburst99/). The updated Starburst99 package can be retrieved from this website as well.

preprint2014arXiv

The Lyman Alpha Reference Sample: III. Properties of the Neutral ISM from GBT and VLA Observations

We present new H I imaging and spectroscopy of the 14 UV-selected star-forming galaxies in the Lyman Alpha Reference Sample (LARS), aimed for a detailed study of the processes governing the production, propagation, and escape of Ly$α$ photons. New H I spectroscopy, obtained with the 100m Green Bank Telescope (GBT), robustly detects the H I spectral line in 11 of the 14 observed LARS galaxies (although the profiles of two of the galaxies are likely confused by other sources within the GBT beam); the three highest redshift galaxies are not detected at our current sensitivity limits. The GBT profiles are used to derive fundamental H I line properties of the LARS galaxies. We also present new pilot H I spectral line imaging of 5 of the LARS galaxies obtained with the Karl G. Jansky Very Large Array (VLA). This imaging localizes the H I gas and provides a measurement of the total H I mass in each galaxy. In one system, LARS 03 (UGC 8335 or Arp 238), VLA observations reveal an enormous tidal structure that extends over 160 kpc from the main interacting systems and that contains $>$10$^9$ M$_{\odot}$ of H I. We compare various H I properties with global Ly$α$ quantities derived from HST measurements. The measurements of the Ly$α$ escape fraction are coupled with the new direct measurements of H I mass and significantly disturbed H I velocities. Our robustly detected sample reveals that both total H I mass and linewidth are tentatively correlated with key Ly$α$ tracers. Further, on global scales, these data support a complex coupling between Ly$α$ propagation and the H I properties of the surrounding medium.

preprint2014arXiv

The Lyman-alpha Reference Sample: I. Survey outline and first results for Markarian 259

The Lyman-alpha reference sample (LARS) is a program with the Hubble Space Telescope (HST) that provides a sample of local universe laboratory galaxies in which to study the astrophysics of the visibility and strength of the Lyman-alpha (Lya) line of hydrogen. This article presents an overview of the survey, its selection function and HST imaging observations. The sample was selected from the GALEX+SDSS catalogue at z=0.028-0.19, in order to allow Lya to be captured with combinations of long pass filters in the Solar Blind Channel (SBC) of HST/ACS. In addition, LARS utilises Halpha and Hbeta narrow, and U, B, i broad-band imaging with ACS and WFC3. In order to study galaxies in which large numbers of Lya photons are produced we demanded an Halpha equivalent width > 100Å. The sample of 14 galaxies covers far UV (FUV) luminosities that overlaps with those of high-z Lya emitters and Lyman Break Galaxies (LBGs). We present the reduction steps used to obtain the Lya images, including our LARS eXtraction software (LaXs) which utilises pixel-by-pixel spectral synthesis fitting to subtract the continuum at Lya. We demonstrate that the use of SBC long pass filters dramatically increase the signal to noise compared to the nominal Lya filter in SBC. To exemplify, we present results for Mrk 259. This irregular galaxy shows extended but strongly asymmetric Lya emission. Spectroscopy from HST/COS show a moderate outflow in the neutral ISM and an asymmetric Lya profile. Radiative transfer modeling is able to reproduce the Lya line profile, confirming the presence of an outflow. The integrated photometry provides a Lya luminosity of 1.3E42 erg/s an equivalent width W(Lya) =45Å and a FUV absolute magnitude -19.2. Mrk 259 would hence be detectable in high-z Lya and LBG surveys. The total Lya escape fraction is 12%, similar to that observed at z>4. (shortened)

preprint2013arXiv

ELT-MOS White Paper: Science Overview & Requirements

The workhorse instruments of the 8-10m class observatories have become their multi-object spectrographs (MOS), providing comprehensive follow-up to both ground-based and space-borne imaging. With the advent of deeper imaging surveys from, e.g., the HST and VISTA, there are a plethora of spectroscopic targets which are already beyond the sensitivity limits of current facilities. This wealth of targets will grow even more rapidly in the coming years, e.g., after the completion of ALMA, the launch of the JWST and Euclid, and the advent of the LSST. Thus, one of the key requirements underlying plans for the next generation of ground-based telescopes, the Extremely Large Telescopes (ELTs), is for even greater sensitivity for optical and infrared spectroscopy. Here we revisit the scientific motivation for a MOS capability on the European ELT, combining updated elements of science cases advanced from the Phase A instrument studies with new science cases which draw on the latest results and discoveries. These science cases address key questions related to galaxy evolution over cosmic time, from studies of resolved stellar populations in nearby galaxies out to observations of the most distant galaxies, and are used to identify the top-level requirements on an 'E-ELT/MOS'. We argue that several of the most compelling ELT science cases demand MOS observations, in highly competitive areas of modern astronomy. Recent technical studies have demonstrated that important issues related to e.g. sky subtraction and multi-object AO can be solved, making fast- track development of a MOS instrument feasible. To ensure that ESO retains world leadership in exploring the most distant objects in the Universe, galaxy evolution and stellar populations, we are convinced that a MOS should have high priority in the instrumentation plan for the E-ELT.

preprint2013arXiv

Far-Ultraviolet Observations of Outflows from IR-Luminous Galaxies

We obtained medium-resolution ultraviolet (UV) spectra between 1150 and 1450 Angstroms of the four UV-bright, infrared (IR)-luminous starburst galaxies IRAS F08339+6517, NGC 3256, NGC 6090, and NGC 7552 using the Cosmic Origins Spectrograph onboard the Hubble Space Telescope. The selected sightlines towards the starburst nuclei probe the properties of the recently formed massive stars and the physical conditions in the starburst-driven galactic superwinds. Despite being metal-rich and dusty, all four galaxies are strong Lyman-alpha emitters with equivalent widths ranging between 2 and 13 Angstroms. The UV spectra show strong P Cygni-type high-ionization features indicative of stellar winds and blueshifted low-ionization lines formed in the interstellar and circumgalactic medium. We detect outflowing gas with bulk velocities of about 400 km/s and maximum velocities of almost 900 km/s. These are among the highest values found in the local universe and comparable to outflow velocities found in luminous Lyman-break galaxies at intermediate and high redshift. The outflow velocities are unlikely to be high enough to cause escape of material from the galactic gravitational potential. However, the winds are significant for the evolution of the galaxies by transporting heavy elements from the starburst nuclei and enriching the galaxy halos. The derived mass outflow rates of ~100 Msol/yr are comparable to, or even higher than the star-formation rates. The outflows can quench star formation and ultimately regulate the starburst as has been suggested for high-redshift galaxies.

preprint2013arXiv

Lyman-alpha line and continuum radiative transfer in a clumpy interstellar medium

Aims: We aim to study the effects of an inhomogeneous interstellar medium (ISM) on the strength and the shape of the Lyman alpha (Lya) line in starburst galaxies. Methods: Using our 3D Monte Carlo Lya radiation transfer code, we study the radiative transfer of Lya, UV and optical continuum photons in homogeneous and clumpy shells of neutral hydrogen and dust surrounding a central source. Our simulations predict the Lya and continuum escape fraction, the Lya equivalent width EW(Lya), the Lya line profile and their dependence on the gas geometry and the main input physical parameters. Results: The ISM clumpiness is found to have a strong impact on the Lya line radiative transfer, entailing a strong dependence of the emergent features of the Lya line (escape fraction, EW(Lya)) on the ISM morphology. Although a clumpy and dusty ISM appears more transparent to radiation (both line and continuum) compared to an equivalent homogeneous ISM of equal dust optical depth, we find that the Lya photons are, in general, still more attenuated than UV continuum radiation. As a consequence, the observed Lya equivalent width (EWobs(Lya)) is lower than the intrinsic one (EWint(Lya)) for nearly all clumpy ISM configurations considered. There are, however, special conditions under which Lya photons escape more easily than the continuum, resulting in an enhanced EWobs(Lya). The requirement for this to happen is that the ISM is almost static (galactic outflows < 200 km/s), extremely clumpy (with density contrasts >10^7 in HI between clumps and the interclump medium), and very dusty (E(B-V) > 0.30). When these conditions are fulfilled the emergent Lya line profile shows no velocity shift and little asymmetry. Given the asymmetry and velocity shifts generally observed in starburst galaxies with Lya emission, we conclude that clumping is unlikely to significantly enhance their relative Lya/UV transmission.

preprint2013arXiv

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.

preprint2013arXiv

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, ξ_ext = RP20_Lya / RP20_Ha, we find ξ_ext to be uncorrelated with total Lya luminosity. However ξ_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.

preprint2013arXiv

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 Sérsic 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.

preprint2013arXiv

Why the globular cluster NGC 6752 contains no sodium-rich second-generation AGB star

(Abridged) Globular clusters host multiple stellar populations showing different sodium enrichments. These various populations can be observed along the main sequence, red giant and horizontal branch phases. Recently it was shown, however, that at least in the globular cluster NGC 6752, no sodium-rich stars are observed along the early asymptotic giant branch, posing an apparent problem for stellar evolution. We present an explanation for this lack of sodium-rich stars in this region of the colour-magnitude diagram.

preprint2012arXiv

A peculiar galaxy appears at redshift 11: properties of a moderate redshift interloper

Laporte et al. (2011) reported a very high redshift galaxy candidate: a lensed J-band dropout (A2667-J1). J1 has a photometric redshift of z=9.6-12, the probability density function for which permits no low or intermediate z solution. We here report new spectroscopic observations of this galaxy with VLT/XShooter, which show clear [OIII]5007AA, Ly-alpha, H-alpha, and H-beta emission and place the galaxy firmly at z=2.082. The oxygen lines contribute only ~25% to the H-band flux, and do not significantly affect the dropout selection of J1. After correcting the broadband fluxes for line emission, we identify two roughly equally plausible natures for A2667-J1: either it is young heavily reddened starburst, or a maximally old system with a very pronounced 4000AA break, upon which a minor secondary burst of star formation is superimposed. Fits show that to make a 3 sigma detection of this object in the B-band (V-band), imaging of depth AB=30.2 (29.5) would be required - despite the relatively bright NIR magnitude, we would need optical data of equivalent depth to the Hubble Ultra Deep Field to rule out the mid-z solution on purely photometric grounds. Assuming that this stellar population can be scaled to the NIR magnitudes of recent HST/WFC3 IR-selected galaxies, we conclude that infeasibly deep optical data AB~32 would be required for the same level of security. There is a population of galaxies at z~2 with continuum colours alone that mimic those of our z=7-12 candidates.

preprint2012arXiv

Evolutionary synthesis models as a tool and guide towards the first galaxies

We summarize the principles and fundamental ingredients of evolutionary synthesis models, which are stellar evolution, stellar atmospheres, the IMF, star-formation histories, nebular emission, and also attenuation from the ISM and IGM. The chapter focusses in particular on issues of importance for predictions of metal-poor and Population III dominated galaxies. We review recent predictions for the main physical properties and related observables of star-forming galaxies based on up-to-date inputs. The predicted metallicity dependence of these quantities and their physical causes are discussed. The predicted observables include in particular the restframe UV-to-optical domain with continuum emission from stars and the ionized ISM, as well as emission lines from H, He, and metals. Based on these predictions we summarize the main observational signatures (emission line strengths, colors etc.), which can be used to distinguish "normal" stellar populations from very metal-poor objects or even Pop III. Evolutionary synthesis models provide an important and fundamental tool for studies of galaxy formation and evolution, from the nearby Universe back to first galaxies. They are used in many applications to interpret existing observations, to predict and guide future missions/instruments, and to allow direct comparisons between state-of-the-art galaxy simulations and observations.

preprint2012arXiv

Modelling high redshift Lyman-alpha Emitters

We present a new model for high redshift Lyman-Alpha Emitters (LAEs) in the cosmological context which takes into account the resonant scattering of Ly-a photons through expanding gas. The GALICS semi-analytic model provides us with the physical properties of a large sample of high redshift galaxies. We implement a gas outflow model for each galaxy based on simple scaling arguments. The coupling with a library of numerical experiments of Ly-a transfer through expanding or static dusty shells of gas allows us to derive the Ly-a escape fractions and profiles. The predicted distribution of Ly-a photons escape fraction shows that galaxies with a low star formation rate have a f_esc of the order of unity, suggesting that, for those objects, Ly-a may be used to trace the star formation rate assuming a given conversion law. In galaxies forming stars intensely, the escape fraction spans the whole range from 0 to 1. The model is able to get a good match to the UV and Ly-a luminosity function (LF) data at 3 < z < 5. We find that we are in good agreement with both the bright Ly-a data and the faint population observed by Rauch et al. (2008) at z=3. Most of the Ly-a profiles of our LAEs are redshifted by the diffusion in the outflow which suppresses IGM absorption. The bulk of the observed Ly-a equivalent width (EW) distribution is recovered by our model, but we fail to obtain the very large values sometimes detected. Predictions for stellar masses and UV LFs of LAEs show a satisfactory agreement with observational estimates. The UV-brightest galaxies are found to show only low Ly-a EWs in our model, as it is reported by many observations of high redshift LAEs. We interpret this effect as the joint consequence of old stellar populations hosted by UV-bright galaxies, and high HI column densities that we predict for these objects, which quench preferentially resonant Ly-a photons via dust extinction.

preprint2012arXiv

Properties of z~3-6 Lyman break galaxies. II. Testing star formation histories and the SFR-mass relation with ALMA and near-IR spectroscopy

We examine the dependence of derived physical parameters of distant Lyman break galaxies (LBGs) on the assumed star formation histories (SFHs), their implications on the SFR-mass relation, and we propose observational tests to better constrain these quantities. We use our SED-fitting tool including nebular emission to analyze a large sample of LBGs, assuming five different star formation histories, extending our first analysis of this sample (de Barros et al. 2012, paper I). In addition we predict the IR luminosities consistently with the SED fits. Compared to "standard" SED fits assuming constant SFR and neglecting nebular lines, assuming variable SFHs yield systematically lower stellar masses, higher extinction, higher SFR, higher IR luminosities, and a wider range of equivalent widths for optical emission lines. Exponentially declining and delayed SFHs yield basically identical results and generally fit best. Exponentially rising SFHs yield similar masses, but somewhat higher extinction. We find significant deviations between the derived SFR and IR luminosity from the commonly used SFR(IR) or SFR(IR+UV) calibration, due to differences in the SFHs and ages. Models with variable SFHs, favored statistically, yield generally a large scatter in the SFR-mass relation. We show that the true scatter in the SFR-mass relation can be significantly larger than inferred using SFR(UV) and/or SFR(IR). Different SFHs, and hence differences in the derived SFR-mass relation and in the specific star formation rates, can be tested/constrained observationally with future IR observations with ALMA. Measurement of emission lines, such as Halpha and [OII]3727, can also provide useful constraints on the SED models. We conclude that our findings of a large scatter in SFR-mass at high-z and an increase of the specific star formation rate above z>~3 (paper I) can be tested observationally. (abriged)

preprint2012arXiv

The Effects of Stellar Rotation. I. Impact on the Ionizing Spectra and Integrated Properties of Stellar Populations

We present a sample of synthetic massive stellar populations created using the Starburst99 evolutionary synthesis code and new sets of stellar evolutionary tracks, including one set that adopts a detailed treatment of rotation. Using the outputs of the Starburst99 code, we compare the populations' integrated properties, including ionizing radiation fields, bolometric luminosities, and colors. With these comparisons we are able to probe the specific effects of rotation on the properties of a stellar population. We find that a population of rotating stars produces a much harder ionizing radiation field and a higher bolometric luminosity, changes that are primarily attributable to the effects of rotational mixing on the lifetimes, luminosities, effective temperatures, and mass loss rates of massive stars. We consider the implications of the profound effects that rotation can have on a stellar population, and discuss the importance of refining stellar evolutionary models for future work in the study of extragalactic, and particularly high-redshift, stellar populations.

preprint2011arXiv

A new perspective on globular clusters, their initial mass function, and their contribution to the stellar halo and to cosmic reionisation

We examine various implications from a dynamical and chemical model of globular clusters (GCs), which successfully reproduces the observed abundance patterns and the multiple populations of stars in these systems assuming chemical enrichment from fast rotating massive stars. Using the model of Decressin et al. (2007) we determine the ratio between the observed, present-day mass of globular clusters and their initial stellar mass as a function of the stellar initial mass function (IMF). We also compute the mass of low-mass stars ejected, and the amount of hydrogen ionising photons emitted by the proto globular clusters. Typically, we find that the initial masses of GCs must be ~8-10 times (or up to 25 times, if second generation stars also escape from GCs) larger than the present-day stellar mass. The present-day Galactic GC population must then have contributed to approximately 5-8% (10-20%) of the low-mass stars in the Galactic halo. We also show that the detection of second generation stars in the Galactic halo, recently announced by different groups, provides a new constraint on the GC initial mass function (GCIMF). These observations appear to rule out a power-law GCIMF, whereas they are compatible with a log-normal one. Finally, the high initial masses also imply that GCs must have emitted a large amount of ionising photons in the early Universe. Our results reopen the question on the initial mass function of GCs, and reinforce earlier conclusions that old GCs could have represented a significant contribution to reionise the inter-galactic medium at high redshift.

preprint2011arXiv

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.

preprint2011arXiv

LBG properties from z~3 to z~6

We analyse the spectral energy distribution (SED) of U, B, V and i-dropout samples from GOODS-MUSIC and we determine their physical properties, such as stellar age and mass, dust attenuation and star formation rate (SFR). Furthermore, we examine how the strength of Lyα emission can be constrained from broad-band SED fits instead of relying in spectroscopy. We use our SED fitting tool including the effects of nebular emission and we explore different star formation histories (SFHs). We find that SEDs are statistically better fitted with nebular emission and exponentially decreasing star formation. Considering this result, stellar mass and star formation rate (SFR) estimations modify the specific SFR (SFR/Mstar) - redshift relation, in comparison to previous studies. Finally, our inferred Lyα properties are in good agreement with the available spectroscopic observations.

preprint2011arXiv

On Lya emission in z~3-6 UV-selected galaxies

Determining Lya properties of distant galaxies is of great interest for various astrophysical studies. We examine how the strength of Lya emission can be constrained from broad-band SED fits instead of relying on spectroscopy. We use our SED fitting tool including the effects of nebular emission, considering in particular Lya emission as a free parameter, and we demonstrate our method with simulations of mock galaxies. Using this tool we analyse a large sample of U, B, V, and i dropout galaxies with multi-band photometry. We find significant trends of the fraction of galaxies with Lya emission increasing both with redshift z and towards fainter magnitude (at fixed z), and similar trends for the Lya equivalent width. Our inferred Lya properties are in good agreement with the available spectroscopic observations and other data. These results demonstrate that the strength of Lya emission in distant star-forming galaxies can be inferred quantitatively from broad-band SED fits, at least statistically for sufficiently large samples with a good photometric coverage.

preprint2011arXiv

On the redshift-evolution of the Lyman-alpha escape fraction and the dust content of galaxies

The Ly-alpha emission line has been proven a powerful tool by which to study evolving galaxies at the highest redshifts. However, in order to use Lya as a physical probe of galaxies, it becomes vital to know the Lya escape fraction (fescLya). Unfortunately, due to the resonant nature of Lya, fescLya may vary unpredictably and requires empirical measurement. Here we compile Lya luminosity functions between redshift z=0 and 8 and, combined with H-alpha and ultraviolet data, assess how fescLya evolves with redshift. We find a strong upwards evolution in fescLya over the range z=0.3-6, which is well-fit by the power-law fescLya \propto (1+z)^ξ, with ξ=(2.57_-0.12^+0.19). This predicts that fescLya should reach unity at z=11.1. By comparing fescLya and E(B-V) in individual galaxies we derive an empirical relationship between fescLya and E(B-V), which includes resonance scattering and can explain the redshift evolution of fescLya between z=0 and 6 purely as a function of the evolution in the dust content of galaxies. Beyond z~6.5, fescLya drops more substantially; an effect attributed to either ionizing photon leakage, or an increase in the neutral gas fraction of the intergalactic medium. While distinguishing between those two scenarios may be extremely challenging, by framing the problem this way we remove the uncertainty of the halo mass from Lya-based tests of reionization. We finally derive a new method by which to estimate the dust content of galaxies based purely upon the observed Lya and UV LFs. These data are characterized by an exponential with an e-folding redshift of ~3.5.

preprint2011arXiv

Optical dropout galaxies lensed by the cluster A2667

We investigate the nature and the physical properties of z, Y and J-dropout galaxies selected behind the lensing cluster A2667. This field is part of our project aimed at identifying z~7-10 candidates accessible to spectroscopic studies, based on deep photometry with ESO/VLT HAWK-I and FORS2 (zYJH and Ks-band images, AB(3 sigma)~26-27) on a sample of lensing clusters extracted from our multi-wavelength combined surveys with SPITZER, HST, and Herschel. In this paper we focus on the complete Y and J-dropout sample, as well as the bright z-dropouts fulfilling the selection criteria by Capak et al. (2011). 10 candidates are selected within the common field of ~33 arcmin2 (effective area once corrected for contamination and lensing dilution). All of them are detected in H and Ks bands in addition to J and/or IRAC 3.6/4.5, with H(AB)~23.4 to 25.2, and have modest magnification factors. Although best-fit photometric redshifts place all these candidates at high-z, the contamination by low-z interlopers is estimated at 50-75% level based on previous studies, and the comparison with the blank-field WIRCAM Ultra-Deep Survey (WUDS). The same result is obtained when photometric redshifts include a luminosity prior, allowing us to remove half of the original sample as likely z~1.7-3 interlopers with young stellar pulations and strong extinction. Two additional sources among the remaining sample could be identified at low-z based on a detection at 24 microns and on the HST z_850 band. These low-z interlopers are not well described by current templates given the large break, and cannot be easily identified based solely on optical and near-IR photometry. Given the estimated dust extinction and high SFRs, some of them could be also detected in the IR or sub-mm bands. After correction for likely contaminants, the observed counts at z>7.5 seem to be in agreement with an evolving LF. (abridged)

preprint2011arXiv

The importance of nebular emission for SED modeling of distant star-forming galaxies

We highlight and discuss the importance of accounting for nebular emission in the SEDs of high redshift galaxies, as lines and continuum emission can contribute significantly or subtly to broad-band photometry. Physical parameters such as the galaxy age, mass, star-formation rate, dust attenuation and others inferred from SED fits can be affected to different extent by the treatment of nebular emission. We analyse a large sample of Lyman break galaxies from z~3-6, and show some main results illustrating e.g. the importance of nebular emission for determinations of the mass-SFR relation, attenuation and age. We suggest that a fairly large scatter in such relations could be intrinsic. We find that the majority of objects (~60-70%) is better fit with SEDs accounting for nebular emission; the remaining galaxies are found to show relatively weak or no emission lines. Our modeling, and supporting empirical evidence, suggests the existence of two categories of galaxies, "starbursts" and "post-starbursts" (lower SFR and older galaxies) among the LBG population, and relatively short star-formation timescales.

preprint2010arXiv

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.

preprint2010arXiv

On the physical properties of z~6-8 galaxies

We analyse the spectral energy distributions (SEDs) of the most distant galaxies discovered with the Hubble Space telescope and from the COSMOS survey and determine their physical properties, such as stellar age and mass, dust attenuation, and star-formation rate. We use our SED fitting tool including the effects of nebular emission to analyse three samples of z ~6-8 galaxies with observed magnitudes J_AB~23 to 29. Our models cover a wide parameter space. We find that the physical parameters of most galaxies cover a wide range of acceptable values. Stellar ages, in particular, are not strongly constrained, even for objects detected longward of the Balmer break. As already pointed out earlier, the effects of nebular lines significantly affect the age determinations of star-forming galaxies at z ~6-8. We find no need for stellar populations with extreme metallicities or other non-standard assumptions (IMF, escape fraction) to explain the observed properties of faint z-dropout galaxies. Albeit with large uncertainties, our fit results show indications of dust attenuation in some of the z ~6-8 galaxies, which have best-fit values of A_V up to ~1. Furthermore, we find a possible trend of increasing dust attenuation with galaxy mass, and a relatively large scatter in specific star-formation rates, SFR/M*. The physical parameters of very high-z galaxies may be more uncertain than indicated by previous studies. Dust attenuation seems also to be present in some z ~6-8 galaxies, and may be correlated with galaxy mass, as is also the case for SFR.

preprint2010arXiv

Predicted UV properties of very metal-poor starburst galaxies

We study the expected properties of starbursts in order to provide the point of reference for interpretation of high-z galaxy surveys and of very metal-poor galaxies. We concentrate mainly on the UV characteristics such as the ionizing spectra, the UV continuum, the Ly alpha and HeII 1640 A line and two-photon continuum emission. We use evolutionary synthesis models covering metallicities from Pop III to solar and a wide range of IMFs. We also combine the synthetic SEDs with the CLOUDY photoionization code for more accurate predictions of nebular emission, and to study possible departures from case B assumed in the synthesis models. The ionizing fluxes, UV continuum properties, and predicted Ly alpha and HeII 1640 A line strengths are presented for synthesis models covering a wider range of parameter space than our earlier calculations. Strong departures from case B predictions are obtained for Ly alpha and two-photon continuum at low metallicities. At low nebular densities both are shown to be enhanced proportionally to the mean energy carried by the Lyman continuum photons emitted by the ionizing source. Larger Ly alpha equivalent widths are therefore predicted at low metallicity. The HeII 1640 A line can be weaker than case B predicts (in terms of flux as well as the equivalent width) due to its ionization parameter dependence and to the enhanced underlying two-photon continuum. Our results have implications for the interpretation of star-forming metal-poor and/or high redshift galaxies, for galaxies among the Ly alpha emitters (LAE) and Lyman Break galaxy (LBG) populations, and for searches of Population III stars in the distant Universe.

preprint2009arXiv

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

preprint2009arXiv

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.

preprint2009arXiv

The H-alpha luminosity function at redshift 2.2: A new determination using VLT/HAWK-I

We aim to place new, strengthened constraints on the luminosity function (LF) of H-alpha emitting galaxies at redshift z=2.2, and to further constrain the instantaneous star-formation rate density of the universe (rho*). We have used the new HAWK-I instrument at ESO-VLT to obtain extremely deep narrow-band (line; NB2090) and broad-band (continuum; Ks) imaging observations. The target field is in the GOODS-South, providing us with a rich multi-wavelength auxiliary data set, which we utilise for redshift confirmation and to estimate dust content. We use this new data to measure the faint-end slope (alpha) of LF(H-alpha) with unprecedented precision. The data are well fit by a Schechter function and also a single power-law, yielding alpha=(-1.72 +/- 0.20) and (-1.77 +/- 0.21), respectively. Thus we are able to confirm the steepening of alpha from low- to high-z predicted by a number of authors and observed at other wavelengths. We combine our LF data-points with those from a much shallower but wider survey at z=2.2 (Geach et al. 2008), constructing a LF spanning a factor of 50 in luminosity. Re-fitting the Schechter parameters, we obtain log L*=(43.07+/-0.22)erg s^-1 ; log phi*=(-3.45+/-0.52)Mpc^-3 ; alpha=(-1.60+/-0.15). We integrate over LF(Halpha) and apply a correction for dust attenuation to determine the instantaneous cosmic star-formation rate density at z=2 without assuming alpha or extrapolating it from lower-z. Our measurement of rho* is (0.215+/-0.090) Msun yr^-1 Mpc^-3, integrated over a range of 37 <log(LHhalpha / erg s^-1) < 47.

preprint2009arXiv

The impact of nebular emission on the ages of z~6 galaxies

We examine the influence of nebular continuous and line emission in high redshift star forming galaxies on determinations of their age, formation redshift and other properties from SED fits. We include nebular emission consistently with the stellar emission in our SED fitting tool and analyse differentially a sample of 10 z~6 galaxies in the GOODS-S field studied earlier by Eyles et al. (2007). We find that the apparent Balmer/4000 Ang breaks observed in a number of z~6 galaxies detected at >~3.6 micron with IRAC/Spitzer can be mimicked by the presence of strong restframe optical emission lines, implying in particular younger ages than previously thought. Applying these models to the small sample of z~6 galaxies, we find that this effect may lead to a typical downward revision of their stellar ages by a factor ~3. In consequence their average formation redshift may drastically be reduced, and these objects may not have contributed to cosmic reionisation at z>6. Extinction and stellar mass estimates may also be somewhat modified, but to a lesser extent. Careful SED fits including nebular emission and treating properly uncertainties and degeneracies are necessary for more accurate determinations of the physical parameters of high-z galaxies.

preprint2007arXiv

Models for Massive Stellar Populations with Rotation

We present and discuss evolutionary synthesis models for massive stellar populations generated with the Starburst99 code in combination with a new set of stellar evolution models accounting for rotation. The new stellar evolution models were compiled from several data releases of the Geneva group and cover heavy-element abundances ranging from twice solar to one fifth solar. The evolution models were computed for rotation velocities on the zero-age main-sequence of 0 and 300 km/s and with the latest revision of stellar mass-loss rates. Since the mass coverage is incomplete, in particular at non-solar chemical composition, our parameter study is still preliminary and must be viewed as exploratory. Stellar population properties computed with Starburst99 and the new evolution models show some marked differences in comparison with models obtained using earlier tracks. Since individual stars now tend to be more luminous and bluer when on the blue side of the Hertzsprung-Russell diagram, the populations mirror this trend. For instance, increases by factors of two or more are found for the light-to-mass ratios at ultraviolet to near-infrared wavelengths, as well as for the output of hydrogen ionizing photons. If these results are confirmed once the evolution models have matured, recalibrations of certain star-formation and initial mass function indicators will be required.

preprint2006arXiv

Subarcsecond Resolution Mid-Infrared Observations of Super Star Clusters in the Antennae (NGC4038/4039)

In this letter we present ground-based subarcsecond mid-infrared imaging and spectroscopy of young super star clusters in the overlap region of the merging galaxies NGC4038/4039 (the Antennae) obtained with the VLT Imager and Spectrometer for mid-Infrared (VISIR). With its unprecedented spatial resolution VISIR begins to resolve the HII/PDR complexes around the star-forming regions for the first time. In the N-band spectra of two young star clusters unexpectedly low polycyclic aromatic hydrocarbon (PAH) emission is observed, compared to what is seen with the Infrared Space Observatory (ISO) and with the Spitzer Space Telescope. We conclude that a large fraction of the PAH emission cannot directly be associated with the super star clusters, but originate from an extended region of at least 215 pc radius around the clusters. In the distribution of [NeII] 12.81 micron emission a highly obscured cluster is discovered that does not have an optical or near-infrared counterpart.

preprint2004arXiv

Stellar Population in LLAGN.I: Ground-based observations

(Abridge): We present a spectroscopic study of the stellar populations of Low Luminosity AGN (LLAGN). Our main goal is to determine whether the stars who live in the innermost (100 pc-scale) regions of these galaxies are in some way related to the emission line properties, which would imply a link between the stellar population and the ionization mechanism. High signal to noise, ground based long-slit spectra in the 3500--5500 A interval were collected for 60 galaxies.Our main findings are: (1) Few LLAGN have a detectable young (< 10 Myr) starburst component, indicating that very massive stars do not contribute significantly to the optical continuum. In particular, no features due to Wolf-Rayet stars were convincingly detected. (2) High Order Balmer absorption lines of HI (HOBLs), on the other hand, are detected in 40% of LLAGN. These features, which are strongest in 100 Myr--1 Gyr intermediate age stellar populations, are accompanied by diluted metal absorption lines and bluer colors than other objects in the sample. (3) These intermediate age populations are very common (50 %) in LLAGN with relatively weak [OI] emission, but rare (10 %) in LLAGN with stronger [OI]. This is intriguing since LLAGN with weak [OI] have been previously hypothesized to be ``transition objects'' in which both an AGN and young stars contribute to the emission-line excitation. Massive stars, if present, are completely outshone by intermediate age and old stars in the optical. This happens in at least a couple of objects where independent UV spectroscopy detects young starbursts not seen in the optical. (4) Objects with predominantly old stars span the whole range of [OI]/Halpha values, but (5) sources with significant young and/or intermediate age populations are nearly all (90%) weak [OI] emitters.

preprint1996arXiv

Combined stellar structure and atmosphere models for massive stars. III. Spectral evolution and revised ionising fluxes of O3-B0 stars

We provide an extensive set of theoretical spectral energy distributions of massive stars derived from our "combined stellar structure and atmosphere models". The calculations covering the entire main sequence evolution for initial masses M_i=20 - 120 Msun (O3-B0 stars of all luminosity classes) are used for a systematic study of the ionizing fluxes of O and early B stars. We demonstrate the importance of accounting simultaneously for non-LTE effects, line blanketing and stellar winds. The main results from our spectra are the following: (1) The flux in the HeII continuum is increased by 2 to 3 (3 to 6) orders of magnitudes compared to predictions from plane parallel non-LTE (LTE) model atmospheres. (2) The flux in the HeI continuum is increased due to non-LTE effects. However, we find that it is also influenced by wind effects. The combined effect of a mass outflow and line blanketing leads to a flatter energy distribution in the HeI continuum, which confirms the results of Sellmaier et al. (1996) for a wider range of stellar parameters. (3) The flux in the Lyman continuum is also modified due to line blanketing and the presence of a stellar winds, although to a lesser degree than the spectrum at higher energies. We derive revised ionizing fluxes for O3 to B0 stars based on the recent Teff and gravity calibrations of Vacca et al. (1996). The total number of Lyman continuum photons is found to be slightly lower than previous derivations. Due to the increased flux in the HeI continuum the hardness ratio of the HeI to H continuum is increased by 1.6 to 2.5 depending on spectral type and luminosity class. A critical discussion of current model assumptions shows that for stars of types later than B0, having weak stellar winds, reliable predictions of ionizing fluxes are not yet possible.