Source author record

Matthew Hayes

Matthew Hayes 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

42works
3topics
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

42 published item(s)

preprint2024arXiv

Lyman Continuum Emission from AGN at 2.3$\lesssim$z$\lesssim$3.7 in the UVCANDELS Fields

We present the results of our search for Lyman continuum (LyC) emitting AGN at redshifts 2.3$\lesssim$z$\lesssim$4.9 from HST WFC3 F275W observations in the UVCANDELS fields. We also include LyC emission from AGN using HST WFC3 F225W, F275W, and F336W found in the ERS and HDUV data. We performed exhaustive queries of the Vizier database to locate AGN with high quality spectroscopic redshifts. In total, we found 51 AGN that met our criteria within the UVCANDELS and ERS footprints. Of these 51, we find 12 AGN had $\geq$4$σ$ detected LyC flux in the WFC3/UVIS images. Using space- and ground-based data from X-ray to radio, we fit the multi-wavelength photometric data of each AGN to a CIGALE SED and correlate various SED parameters to the LyC flux. KS-tests of the SED parameter distributions for the LyC-detected and non-detected AGN showed they are likely not distinct samples. However, we find that X-ray luminosity, star-formation onset age, and disk luminosity show strong correlations relative to their emitted LyC flux. We also find strong correlation of the LyC flux to several dust parameters, i.e., polar and toroidal dust emission, 6 $μm$ luminosity, and anti-correlation with metallicity and $A_{FUV}$. We simulate the LyC escape fraction ($f_{esc}$) using the CIGALE and IGM transmission models for the LyC-detected AGN and find an average $f_{esc}$$\simeq$18%, weighted by uncertainties. We stack the LyC flux of subsamples of AGN according to the wavelength continuum region in which they are detected and find no significant distinctions in their LyC emission, although our $sub-mm\ detected$ F336W sample shows the brightest stacked LyC flux. These findings indicate that LyC-production and -escape in AGN is more complicated than the simple assumption of thermal emission and a 100% escape fraction. Further testing of AGN models with larger samples than presented here is needed.

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

CLASSY III: The Properties of Starburst-Driven Warm Ionized Outflows

We report the results of analyses of galactic outflows in a sample of 45 low-redshift starburst galaxies in the COS Legacy Archive Spectroscopic SurveY (CLASSY), augmented by five additional similar starbursts with COS data. The outflows are traced by blueshifted absorption-lines of metals spanning a wide range of ionization potential. The high quality and broad spectral coverage of CLASSY data enable us to disentangle the absorption due to the static ISM from that due to outflows. We further use different line multiplets and doublets to determine the covering fraction, column density, and ionization state as a function of velocity for each outflow. We measure the outflow's mean velocity and velocity width, and find that both correlate in a highly significant way with the star-formation rate, galaxy mass, and circular velocity over ranges of four orders-of-magnitude for the first two properties. We also estimate outflow rates of metals, mass, momentum, and kinetic energy. We find that, at most, only about 20% of silicon created and ejected by supernovae in the starburst is carried in the warm phase we observe. The outflows' mass-loading factor increases steeply and inversely with both circular and outflow velocity (log-log slope $\sim$ -1.6), and reaches $\sim 10$ for dwarf galaxies. We find that the outflows typically carry about 10 to 100% of the momentum injected by massive stars and about 1 to 20% of the kinetic energy. We show that these results place interesting constraints on, and new insights into, models and simulations of galactic winds.

preprint2022arXiv

CLASSY V: The impact of aperture effects on the inferred nebular properties of local star-forming galaxies

Strong nebular emission lines are an important diagnostic tool for tracing the evolution of star-forming galaxies across cosmic time. However, different observational setups can affect these lines, and the derivation of the physical nebular properties. We analyze 12 local star-forming galaxies from the COS Legacy Spectroscopy SurveY (CLASSY) to assess the impact of using different aperture combinations on the determination of the physical conditions and gas-phase metallicity. We compare optical spectra observed with the SDSS aperture, which has a 3" of diameter similar to COS, to IFU and longslit spectra, including new LBT/MODS observations of five CLASSY galaxies. We calculate the reddening, electron densities and temperatures, metallicities, star formation rates, and equivalent widths (EWs). We find that measurements of the electron densities and temperatures, and metallicity remained roughly constant with aperture size, indicating that the gas conditions are relatively uniform for this sample. However, using the IFU observations of 3 galaxies, we find that the E(B-V) values derived from the Balmer ratios decrease ( by up to 53%) with increasing aperture size. The values change most significantly in the center of the galaxies, and level out near the COS aperture diameter of 2.5". We examine the relative contributions from the gas and stars using the H$α$ and [OIII] $λ$5007 EWs as a function of aperture light fraction, but find little to no variations within a given galaxy. These results imply that the optical spectra provide nebular properties appropriate for the FUV CLASSY spectra, even when narrow 1.0" long-slit observations are used.

preprint2022arXiv

The COS Legacy Archive Spectroscopy SurveY (CLASSY) Treasury Atlas

Far-ultraviolet (FUV; ~1200-2000 angstroms) spectra are fundamental to our understanding of star-forming galaxies, providing a unique window on massive stellar populations, chemical evolution, feedback processes, and reionization. The launch of JWST will soon usher in a new era, pushing the UV spectroscopic frontier to higher redshifts than ever before, however, its success hinges on a comprehensive understanding of the massive star populations and gas conditions that power the observed UV spectral features. This requires a level of detail that is only possible with a combination of ample wavelength coverage, signal-to-noise, spectral-resolution, and sample diversity that has not yet been achieved by any FUV spectral database. We present the COS Legacy Spectroscopic SurveY (CLASSY) treasury and its first high level science product, the CLASSY atlas. CLASSY builds on the HST archive to construct the first high-quality (S/N_1500 >~ 5/resel), high-resolution (R~15,000) FUV spectral database of 45 nearby (0.002 < z < 0.182) star-forming galaxies. The CLASSY atlas, available to the public via the CLASSY website, is the result of optimally extracting and coadding 170 archival+new spectra from 312 orbits of HST observations. The CLASSY sample covers a broad range of properties including stellar mass (6.2 < logM_star(M_sol) < 10.1), star formation rate (-2.0 < log SFR (M_sol/yr) < +1.6), direct gas-phase metallicity (7.0 < 12+log(O/H) < 8.8), ionization (0.5 < O_32 < 38.0), reddening (0.02 < E(B-V < 0.67), and nebular density (10 < n_e (cm^-3) < 1120). CLASSY is biased to UV-bright star-forming galaxies, resulting in a sample that is consistent with z~0 mass-metallicity relationship, but is offset to higher SFRs by roughly 2 dex, similar to z >~2 galaxies. This unique set of properties makes the CLASSY atlas the benchmark training set for star-forming galaxies across cosmic time.

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 Lyman Alpha Spectral Database (LASD)

Lyman $α$ (Ly$α$) emission from star-forming galaxies is an important tool to study a large range of astrophysical questions: it has the potential to carry information about the source galaxy, its nearby circumgalactic medium, and also the surrounding intergalactic medium. Much observational and theoretical work has therefore focused on understanding the details of this emission line. These efforts have been hampered, however, by an absence of spectroscopic reference samples that can be used both as comparisons for observational studies and as critical tests for theoretical work. For this reason we have compiled a large sample of Ly$α$ spectra, at both low and high redshift, and created a publicly available online database, at http://lasd.lyman-alpha.com. The Lyman alpha Spectral Database (LASD) hosts these spectra, as well as large set of spectral and kinematic quantities that have been homogeneously measured for the entire sample. These can easily be viewed online and downloaded in tabular form. The LASD has the capacity for users to easily upload their own Ly$α$ spectra, and all the same spectral measurements will be made, reported, and ingested into the database. We actively invite the community to do so, and the LASD is intended to be a long-term community resource. In this paper we present the design of the database as well as descriptions of the underlying algorithms and the initial Ly$α$ emitter samples that are in the database.

preprint2022arXiv

The molecular gas resolved by ALMA in the low-metallicity dwarf merging galaxy Haro 11

The physical mechanisms for starburst or quenching in less massive ($M_* < 10^{10} M_{\odot}$) galaxies are unclear. The merger is one of the inescapable processes referred to as both starburst and quenching in massive galaxies. However, the effects of the merger on star formation in dwarf galaxies and their evolution results are still uncertain. We aim to explore how to trigger and quench star formation in dwarf galaxies by studying the metal-poor gas-rich dwarf mergers based on the multi-band observations at a spatial resolution of $\sim$ 460 pc. We use the archival data of ALMA (band 3, 8) and VLT/MUSE to map CO($J=$1-0), [CI]($^3$P$_1 - ^3$P$_0$), and H$α$ emission in one of the most extreme starburst merging dwarf galaxies, Haro 11. We find the molecular gas is assembled around the central two star-forming regions. The molecular/ionized gas and stellar components show complex kinematics, indicating that the gas is probably at a combined stage of collision of clouds and feedback from star formation. The peak location and distribution of [CI](1-0) strongly resemble the CO(1-0) emission, meaning that it might trace the same molecular gas as CO in such a dwarf merger starburst galaxy. The enhancement of line ratios ($\sim 0.5$) of [CI]/CO around knot C is probably generated by the dissociation of CO molecules by cosmic rays and far-ultraviolet photons. Globally, Haro 11 and its star-forming regions share similar SFEs as the high-$z$ starburst galaxies or the clumps in nearby (U)LIRGs. Given the high SFE, sSFR, small stellar mass, low metallicity, and deficient HI gas, Haro 11 could be an analog of high-$z$ dwarf starburst and the potential progenitor of the nearby less massive elliptical galaxies. The significantly smaller turbulent pressure and viral parameter will probably trigger the intense starbursts. We also predict that it will quench at $M_* < 8.5 \times 10^9 M_{\odot}$.

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

The Lyman Alpha Reference Sample XII: Morphology of extended Lyman alpha emission in star-forming galaxies

We use Hubble space telescope data of 45 nearby star-forming galaxies to investigate properties of Lyman-alpha (Ly$α$) halos, Ly$α$ morphology, and the star-forming characteristics of galaxies. We study how the morphology of Ly$α$ emission is related to other Ly$α$ observables. Furthermore, we study the interdependencies of Ly$α$ morphological quantities. We studied the spatial extent of Ly$α$ using surface brightness profiles in i) using circular apertures and ii) within faint Ly$α$ isophotes. We measured the average intensity and the size of the regions with a high star formation rate density. The morphology of the galaxies was quantified by computing centroid position, axis ratio, and position angle in the Ly$α$, ultraviolet continuum, and I band maps. We found that galaxies with more extended star-forming regions possess larger Ly$α$ halos. Furthermore, galaxies with more elongated Ly$α$ morphology are also more extended in Ly$α$. Our data suggest that Ly$α$ bright galaxies appear rounder in their Ly$α$ morphology, and there is less of a contribution from their Ly$α$ halo to their overall luminosity. We compared our results with studies at high redshift and found that whilst the Ly$α$ extent in the inner regions of the galaxies in our sample are similar to the high$-z$ Ly$α$ emitters (LAEs), Ly$α$ halos are more extended in high$-z$ LAEs. Our analysis suggests that the Ly$α$ morphology affects the measurement of other observable quantities concerning Ly$α$ emission, and some of the conclusions drawn from high redshift LAEs might be biased towards galaxies with specific Ly$α$ shapes. In particular, faint Ly$α$ emitters have larger Ly$α$ scale lengths and halo fractions. This implies that faint Ly$α$ emitters are harder to detect at high redshift than previously believed.

preprint2020arXiv

Lyman Alpha Reference Sample: X. Predicting Lyman alpha output from starforming galaxies using multivariate regression

Understanding the production and escape of Lyman $α$ (Ly$α$) radiation from star-forming galaxies is a long standing problem in astrophysics. The ability to predict the Ly$α$ luminosity of galaxies would open up new ways of exploring the Epoch of Reionization (EoR), and to estimate Ly$α$ emission from galaxies in cosmological simulations where radiative transfer calculations cannot be done. We apply multivariate regression methods to the Lyman Alpha Reference Sample dataset to obtain a relation between the galaxy properties and the emitted Ly$α$. The derived relation predicts the Ly$α$ luminosity of our galaxy sample to good accuracy, regardless of whether we consider only direct observables (root-mean-square (RMS) dispersion around the relation of $\sim 0.19$ dex) or derived physical quantities (RMS $\sim 0.27$ dex). We confirm the predictive ability on a separate sample of compact star-forming galaxies and find that the prediction works well, but that aperture effects on measured Ly$α$ luminosity may be important, depending on the redshift of the galaxy. We apply statistical feature selection techniques to determine an order of importance of the variables in our dataset, enabling future observations to be optimized for predictive ability. When using physical variables, we are able to determine that the most important predictive parameters are, in order, star formation rate, dust extinction, compactness and the gas covering fraction. We discuss the application of our results in terms of studying the EoR and intensity mapping experiments.

preprint2016arXiv

O VI Emission Imaging of a Galaxy with the Hubble Space Telescope: a Warm Gas Halo Surrounding the Intense Starburst SDSS J115630.63+500822.1

We report results from a new HST study of the OVI 1032,1038Å doublet in emission around intensely star-forming galaxies. The programme aims to characterize the energy balance in starburst galaxies and gas cooling in the difficult-to-map coronal temperature regime of 2-5 x $10^5$K. We present the first resolved image of gas emission in the OVI line. Our target, SDSS J1156+5008, is very compact in the continuum but displays OVI emission to radii of 23 kpc. The surface brightness profile is well fit by an exponential with a scale of 7.5kpc. This is ten times the size of the photoionized gas, and we estimate that 1/6 the total OVI luminosity comes from resonantly scattered continuum radiation. Spectroscopy - which closely resembles a stacked sample of archival spectra - confirms the OVI emission, and determines the column density and outflow velocity from blueshifted absorption. The combination of measurements enables several new calculations with few assumptions. The OVI regions fill only ~$10^{-3}$ of the volume. By comparing the cooling time with the cloud sound-crossing time, the cooling distance with the size, and the pressure in the OVI and nebular gas, we conclude that the OVI-bearing gas cannot have been lifted to the scale height at this temperature, and must be cooling in situ through this coronal temperature regime. The coronal phase contains ~1% of the ionized mass, and its kinetic energy is currently ~1% of the budget set by supernova feedback. However a much larger amount of the gas must have cooled through this phase during the star formation episode. The outflow exceeds the escape velocity and the gas may become unbound, but it will recombine before it escapes and become visible to Lyman (and OI) spectroscopy. The mapping of this gas represents a crucial step in further constraining galaxy formation scenarios and guiding the development of future satellites.

preprint2016arXiv

Spectro-polarimetry confirms central powering of a Ly$α$ nebula at z=3.09

We present a follow-up study to the imaging polarimetry performed by Hayes et. al. 2011 on LAB1 in the SSA22 protocluster region. Arguably the most well-known Lyman-$α$ "blob", this radio-quiet emission-line nebula likely hosts a galaxy which is either undergoing significant star formation or hosts an AGN, or both. We obtain deep, spatially resolved spectro-polarimetry of the Ly$α$ emission and detect integrated linear polarization of $9$-$13\%\pm2$-$3\%$ at a distance of approximately 15 kpc north and south of the peak of the Lyman-$α$ surface brightness with polarization vectors lying tangential to the galactic central source. In these same regions, we also detect a wavelength dependence in the polarization which is low at the center of the Ly$α$ line profile and rises substantially in the wings of the profile. These polarization signatures are easily explained by a weak out-flowing shell model. The spectral dependence of the polarization presented here provide a framework for future observations and interpretations of the southern portion of LAB1 in that any model for this system must be able to reproduce this particular spectral dependence. However, questions still remain for the northern-most spur of LAB1. In this region we detect total linear polarization of between $3$ and $20\%$ at the $5\%$ significance level. Simulations predict that polarization should increase with radius for a symmetric geometry. That the northern spur does not suggests either that this region is not symmetric (which is likely) and exhibits variations in columns density, or that it is kinematically distinct from the rest of LAB1 and powered by another mechanism altogether.

preprint2015arXiv

Kinematics of Haro11 - the miniature Antennae

(abridged) Luminous blue compact galaxies are among the most active galaxies in the local universe in terms of their star formation rate per unit mass. They may be seen as the local analogs of higher redshift Lyman Break Galaxies. Studies of their kinematics is key to understanding what triggers their unusually active star formation In this work we investigate the kinematics of stars and ionised gas in Haro11, one of the most luminous blue compact galaxies in the local universe. Previous works have indicated that many such galaxies may be triggered by galaxy mergers. We have employed Fabry-Perot interferometry, long-slit spectroscopy and Integral Field Unit (IFU) spectroscopy to explore the kinematics of Haro11. We target the near infrared Calcium triplet to derive the stellar velocity field and velocity dispersion. Ionised gas is analysed through emission lines from hydrogen, [OIII] , and [SIII]. When spectral resolution and signal to noise allows we investigate the the line profile in detail and identify multiple velocity components when present. We find that to first order, the velocity field and velocity dispersions derived from stars and ionised gas agree. Hence the complexities reveal real dynamical disturbances providing further evidence for a merger in Haro11. Through decomposition of emission lines we find evidence for kinematically distinct components, for instance a tidal arm behind the galaxy. The ionised gas velocity field can be traced to large galactocentric radii, and shows significant velocity dispersion even far out in the halo. We discuss the origin of the line width, and interpreted as virial motions it indicates a mass of ~1E11 M_sun. Morphologically and kinematically Haro11 shows many resemblances with the famous Antennae galaxies, but is much denser which is the likely explanation for the higher star formation efficiency in Haro11.

preprint2015arXiv

Lyman alpha Emitting Galaxies in the Nearby Universe

The Lya emission line of HI is intrinsically the brightest feature in the spectrum of astrophysical nebulae, making it a very attractive observational tool with which to survey galaxies. Moreover as a UV resonance line, Lya possesses several unique characteristics that make it useful to study the ISM and ionizing stellar population at all cosmic epochs. In this review I present a summary of Lya observations of galaxies in the nearby universe. At UV magnitudes reachable with current facilities, only ~5% of the local galaxy population shows a Lya equivalent width (EW_Lya) that exceeds 20Å. This fraction increases dramatically at higher z, but only in the local universe can we study galaxies in detail and assemble unprecedented multi-wavelength datasets. I discuss many local Lya observations, showing that when galaxies show net Lya emission, they ubiquitously produce large halos of scattered Lya, that dominate the integrated luminosity. We discuss how global EW_Lya and the Lya escape fraction (fescLya) are higher (EW_Lya >~ 20Å and fescLya> 10%) in galaxies that represent the less massive and younger end of the distributions for local objects. This is connected with various properties, such that Lya-emitters have lower metallicities (median value of 12+log(O/H) ~ 8.1) and dust reddening. However, the presence of galactic outflows is also vital to Doppler shift the Lya line out of resonance with the HI, as high EW_Lya is found only among galaxies with winds faster than ~50km/s. The evidence is then assembled into a coherent picture, and the requirement for star formation driven feedback is discussed in the context of an evolutionary sequence where the ISM is accelerated and/or subject to fluid instabilities, which reduce the scattering of Lya. Concluding remarks take the form of perspectives upon the most pressing questions that can be answered by observation.

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 Lyman Continuum Escape Fraction of Low-Mass Star-Forming Galaxies at z~1

To date no direct detection of Lyman continuum emission has been measured for intermediate--redshift z~1 star-forming galaxies . We combine HST grism spectroscopy with GALEX UV and ground--based optical imaging to extend the search for escaping Lyman continuum to a large (~600) sample of z~1 low-mass, moderately star-forming galaxies selected initially on H$α$ emission. The characteristic escape fraction of LyC from SFGs that populate this parameter space remains weakly constrained by previous surveys, but these faint SFGs are assumed to play a significant role in the reionization of neutral hydrogen in the intergalactic medium (IGM) at high redshift (z>6). We do not make an unambiguous detection of escaping LyC radiation from this $z\sim1$ sample, individual non--detections to constrain the absolute Lyman continuum escape fraction, $f_{esc}$<2.1% (3$σ$). We measure upper limits of $f_{esc}$<9.6% from a sample of SFGs selected on high H$α$ equivalent width (EW>200Å), which are thought to be close analogs of high redshift sources of reionization. For reference, we also present an emissivity--weighted escape fraction which is useful as a measurement of the general contribution of the SFGs to the z~1 ionizing UV background. In the discussion, we consider the implications of these intermediate redshift constraints for the reionization of hydrogen in the intergalactic medium at high (z > 6) redshift. If we assume our $z\sim1$ SFGs, for which we measure this emissivity-weighted $f_{esc}$, are analogs to the high redshift sources of reionization, we find is difficult reconcile reionization by faint (M$_{UV}<-13$) SFGs with a low escape fraction ($f_{esc}<$3%), with constraints from independent high redshift observations. If $f_{esc}$ evolves with redshift, reionization by SFGs may be consistent with observations from Planck.

preprint2015arXiv

The super star cluster driven feedback in ESO338-IG04 and Haro 11

The stellar content of young massive star clusters emit large amounts of Lyman continuum photons and inject momentum into the inter stellar medium (ISM) by the strong stellar winds of the most massive stars in the cluster. When the most massive stars explode as supernovae, large amounts of mechanical energy are injected in the ISM. A detailed study of the ISM around these massive cluster provides insights on the effect of cluster feedback. We present high quality integral field spectroscopy taken with VLT/MUSE of two starburst galaxies: ESO 338-IG04 and Haro 11. Both galaxies contain a significant number of super star clusters. The MUSE data provide us with an unprecedented view of the state and kinematics of the ionized gas in the galaxy allowing us to study the effect of stellar feedback on small and large spatial scales. We present our recent results on studying the ISM state of these two galaxies. The data of both galaxies show that the mechanical and ionization feedback of the super star clusters in the galaxy modify the state and kinematics of the ISM substancially by creating highly ionized bubbles around the cluster, making the central part of the galaxy highly ionized. This shows that the HII regions around the individual clusters are density bounded, allowing the ionizing photons to escape and ionize the ISM further out.

preprint2014arXiv

Studying reionization with the next generation of Ly-alpha emitter surveys

We study the prospects for constraining the ionized fraction of the intergalactic medium (IGM) at $z>6$ with the next generation of large Ly$α$ emitter surveys. We make predictions for the upcoming Subaru Hyper Suprime-Cam (HSC) Ly$α$ survey and a hypothetical spectroscopic survey performed with the James Webb Space Telescope (JWST). Considering various scenarios where the observed evolution of the Ly$α$ luminosity function of Ly$α$ emitters at $z>6$ is explained partly by an increasingly neutral IGM and partly by intrinsic galaxy evolution, we show how clustering measurements will be able to distinguish between these scenarios. We find that the HSC survey should be able to detect the additional clustering induced by a neutral IGM if the global IGM neutral fraction is greater than $\sim$20 per cent at $z=6.5$. If measurements of the Ly$α$ equivalent widths (EWs) are also available, neutral fractions as small as 10 per cent may be detectable by looking for correlation between the EW and the local number density of objects. In this case, if it should turn out that the IGM is significantly neutral at $z=6.5$ and the intrinsic EW distribution is relatively narrow, the observed EWs can also be used to construct a map of the locations and approximate sizes of the largest ionized regions. For the JWST survey, the results appear a bit less optimistic. Since such surveys probe a large range of redshifts, the effects of the IGM will be mixed up with any intrinsic galaxy evolution that is present, making it difficult to disentangle the effects. However, we show that a survey with the JWST will have a possibility of observing a large group of galaxies at $z\sim7$, which would be a strong indication of a partially neutral IGM.

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

Escape of Lyman continuum radiation from local galaxies - Detection of leakage from the young starburst Tol 1247-232

The escape fraction of hydrogen ionizing photons (f_esc) from galaxies has been suggested to be evolving with time, but the picture is far from clear. While evidence for significant escape fractions has been found at high redshifts in several studies, the picture looks different in the more nearby universe. Here, we apply a new background subtraction routine on archival data from the Far Ultraviolet Spectroscopic Explorer (FUSE), in order to study local galaxies in search for possible Lyman Continuum (LyC) leakage. In the process, for the first time a stacked spectrum in the LyC is produced for local galaxies. With this small sample, we also make a more tentative approach to look for possible correlations between f_esc and physical parameters such as internal absorption E(B-V)_i, mass, H I mass, specific star formation rate (SSFR), metallicity, and Ly-alpha emission. Eight star forming galaxies with redshifts z > 0.015 from the FUSE archive were re-examined. Also, a sub-sample of an additional four galaxies with lower redshifts were included, for which the escape fraction was estimated from residual flux in the low ionization interstellar C II(1036Å) line. Out of the eight galaxies, only one was found to have significant LyC leakage, Tol 1247-232 (S/N=5.2). This is the second detection of a leaking galaxy in the local universe. We find an absolute escape fraction for Tol 1247-232 of f_esc = 2.4(+0.9/-0.8) %. The stacked sample show an excess in the LyC with f_esc = 1.4(+0.6/-0.5) %, but we note that there might be important selection biases involved. With the small sample, we suggest a possible trend for higher f_esc with lower mass and with enhanced SSFR. None of the galaxies with high values of E(B-V)_i were found to show any sign of leakage.

preprint2013arXiv

Magellan/MMIRS near-infrared multi-object spectroscopy of nebular emission from star forming galaxies at 2<z<3

To investigate the ingredients, which allow star-forming galaxies to present Lyalpha line in emission, we studied the kinematics and gas phase metallicity (Z) of the interstellar medium. We used multi-object NIR spectroscopy with Magellan/MMIRS to study nebular emission from z=2-3 star-forming galaxies discovered in 3 MUSYC fields. We detected emission lines from four active galactic nuclei and 13 high-z star-forming galaxies, including Halpha lines down to a flux of 4.E-17 erg/sec/cm^2. This yielded 7 new redshifts. The most common emission line detected is [OIII]5007, which is sensitive to Z. We were able to measure Z for 2 galaxies and to set upper(lower) limits for another 2(2). The Z values are consistent with 0.3<Z/Zsun<1.2. Comparing the Lyalpha central wavelength with the systemic redshift, we find Delta_v(Lyalpha-[OIII])=70-270 km/sec. High-redshift star-forming galaxies, Lyalpha emitting (LAE) galaxies, and Halpha emitters appear to be located in the low mass, high star-formation rate (SFR) region of the SFR versus stellar mass diagram, confirming that they are experiencing burst episodes of star formation, which are building up their stellar mass. Their Zs are consistent with the relation found for z<2.2 galaxies in the Z versus stellar mass plane. The measured Delta_v(Lyalpha-[OIII]) values imply that outflows of material, driven by star formation, could be present in the z=2-3 LAEs of our sample. Comparing with the literature, we note that galaxies with lower Z than ours are also characterized by similar Delta_v(Lyalpha-[OIII]) velocity offsets. Strong [OIII] is detected in many Lyalpha emitters. Therefore, we propose the Lyalpha/[OIII] flux ratio as a tool for the study of high-z galaxies; while influenced by Z, ionization, and Lyalpha radiative transfer in the ISM, it may be possible to calibrate this ratio to primarily trace one of these effects.

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.

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

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

Project Lyman: Quantifying 11 Gyrs of Metagalactic Ionizing Background Evolution

The timing and duration of the reionization epoch is crucial to the emergence and evolution of structure in the universe. The relative roles that star-forming galaxies, active galactic nuclei and quasars play in contributing to the metagalactic ionizing background across cosmic time remains uncertain. Deep quasar counts provide insights into their role, but the potentially crucial contribution from star-formation is highly uncertain due to our poor understanding of the processes that allow ionizing radiation to escape into the intergalactic medium (IGM). The fraction of ionizing photons that escape from star-forming galaxies is a fundamental free parameter used in models to "fine-tune" the timing and duration of the reionization epoch that occurred somewhere between 13.4 and 12.7 Gyrs ago (redshifts between 12 > z > 6). However, direct observation of Lyman continuum (LyC) photons emitted below the rest frame \ion{H}{1} ionization edge at 912 Å is increasingly improbable at redshifts z > 3, due to the steady increase of intervening Lyman limit systems towards high z. Thus UV and U-band optical bandpasses provide the only hope for direct, up close and in depth, observations of the types of environment that favor LyC escape. By quantifying the evolution over the past 11 billion years (z < 3) of the relationships between LyC escape and local and global parameters ..., we can provide definitive information on the LyC escape fraction that is so crucial to answering the question of, how did the universe come to be ionized? Here we provide estimates of the ionizing continuum flux emitted by "characteristic" (L_{uv}^*) star-forming galaxies as a function of look back time and escape fraction, finding that at z = 1 (7.6 Gyrs ago) L_{uv}^* galaxies with an escape fraction of 1% have a flux of 10^{-19} ergs cm^{-2} s^{-1} Å^{-1}.

preprint2011arXiv

Central Powering of the Largest Lyman-alpha Nebula is Revealed by Polarized Radiation

High-redshift Lyman-alpha blobs are extended, luminous, but rare structures that appear to be associated with the highest peaks in the matter density of the Universe. Their energy output and morphology are similar to powerful radio galaxies, but the source of the luminosity is unclear. Some blobs are associated with ultraviolet or infrared bright galaxies, suggesting an extreme starburst event or accretion onto a central black hole. Another possibility is gas that is shock excited by supernovae. However some blobs are not associated with galaxies, and may instead be heated by gas falling into a dark matter halo. The polarization of the Ly-alpha emission can in principle distinguish between these options, but a previous attempt to detect this signature returned a null detection. Here we report on the detection of polarized Ly-alpha from the blob LAB1. Although the central region shows no measurable polarization, the polarized fraction (P) increases to ~20 per cent at a radius of 45 kpc, forming an almost complete polarized ring. The detection of polarized radiation is inconsistent with the in situ production of Ly-alpha photons, and we conclude that they must have been produced in the galaxies hosted within the nebula, and re-scattered by neutral hydrogen.

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

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

The Massive Star Clusters in the Dwarf Merger ESO 185-IG13: is the Red Excess Ubiquitous in Starbursts?

We have investigated the starburst properties of the luminous blue compact galaxy ESO 185-IG13. The galaxy has been imaged with the high resolution cameras onboard to the Hubble Space Telescope. From the UV to the IR, the data reveal a system shaped by hundreds of young star clusters, and fine structures, like a tidal stream and a shell. The presence of numerous clusters and the perturbed morphology indicate that the galaxy has been involved in a recent merger event. Using previous simulations of shell formation in galaxy mergers we constrain potential progenitors of ESO 185-IG13. The analysis of the star cluster population is used to investigate the properties of the present starburst and to date the final merger event, which has produced hundreds of clusters younger than 100 Myr. We have found a peak of cluster formation only 3.5 Myr old. A large fraction of these clusters will not survive after 10-20 Myr, due to the "infant mortality" caused by gas expulsion. However, this sample of clusters represents an unique chance to investigate the youngest phases of cluster evolution. As already observed in the analog blue compact galaxy Haro 11, a fraction of young clusters are affected by a flux excess at wavelengths longer than 8000 Å. Ages, masses, and extinctions of clusters with this NIR excess are estimated from UV and optical data. We discuss similarities and differences of the observed NIR excess in ESO 185-IG13 clusters with other cases in the literature. The cluster ages and masses are used to distinguish among the potential causes of the excess. We observe, as in Haro 11, that the use of the IR and the (commonly used) I band data results in overestimates of age and mass in clusters affected by the NIR excess. This has important implications for a number of related studies of star clusters.

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 origin of the red excess in very young super star clusters: the case of SBS 0335-052E

The spectral energy distribution analysis of very young unresolved star clusters challenges our understanding of the cluster formation process. Studies of resolved massive clusters in the Milky Way and in the nearby Magellanic Clouds show us that the contribution from photoionized gas is very important during the first Myr of cluster evolution. We present our models which include both a self-consistent treatment of the photoionized gas and the stellar continuum and quantify the impact of such nebular component on the total flux of young unresolved star clusters. A comparison with other available models is considered. The very young star clusters in the SBS 0335-052E dwarf starburst galaxy are used as a test for our models. Due to the low metallicity of the galactic medium our models predict a longer lasted nebular phase which contributes between 10-40% of the total near infrared (NIR) fluxes at around 10 Myr. We propose thus a possible solution for the observed flux excess in the 6 bright super star clusters of SBS 0335-052E. Reines et al. showed that the observed cluster fluxes, in the red-optical and NIR range, sit irreconcilably above the provided stellar continuum models. We find that in the age range estimated from the H_alpha emission we can explain the red excess in all the 6 super star clusters as due to nebular emission, which at cluster ages around 10 Myr still affects the NIR wavebands substantially.

preprint2009arXiv

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

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

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

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 Lyman alpha morphology of local starburst galaxies: release of calibrated images

We present reduced and calibrated high resolution Lyman-alpha (Lya) images for a sample of 6 local star forming galaxies. Targets were selected to represent a range in luminosity and metallicity and to include both known Lya emitters and non-emitters. Far ultraviolet imaging was carried out with the Solar Blind Channel of the ACS on HST in the F122M (Lya on-line) and F140LP (continuum) filters. The resulting Lya images are the product of careful modeling of both the stellar and nebular continua, facilitated by supporting HST imaging in Ha and 5 continuum bands, combined with Starburst99 evolutionary synthesis models, and prescriptions for dust extinction on the continuum. In all, the resulting morphologies in Lya, Ha, and UV-continuum are qualitatively very different and we show that the bulk of Lya emerges in a diffuse component resulting from resonant scattering events. Lya escape fractions, computed from integrated Ha luminosities and recombination theory, are found never to exceed 14%. Internal dust extinction is estimated in each pixel and used to correct Lya fluxes. However, the extinction corrections are far too small (factors from 2.6 to infinity) to reconcile the global Lya luminosities with recombination theory. Surprisingly, when comparing the global equivalent widths of Lya and Ha, the two quantities appear anti-correlated, which may be due to the evolution of mechanical feedback. This calls for caution in the interpretation of Lya observations. The images presented have a physical resolution 3 orders of magnitude better than attainable at high-z from the ground with current instrumentation and our images may therefore serve as useful templates for comparing with observations and modeling of primeval galaxy formation. We therefore provide the reduced Lya, Ha, and continuum images to the community.