The reionization of hydrogen is closely linked to the first structures in the universe, so understanding the timeline of reionization promises to shed light on the nature of these early objects. In particular, transmission of Lyman alpha (Ly$\alpha$) from galaxies through the intergalactic medium (IGM) is sensitive to neutral hydrogen in the IGM, so can be used to probe the reionization timeline. In this work, we implement an improved model of the galaxy UV luminosity to dark matter halo mass relation to infer the volume-averaged fraction of neutral hydrogen in the IGM from Ly$\alpha$ observations. Many models assume that UV-bright galaxies are hosted by massive dark matter haloes in overdense regions of the IGM, so reside in relatively large ionized regions. However, observations and N-body simulations indicate that scatter in the UV luminosity-halo mass relation is expected. Here, we model the scatter (though we assume the IGM topology is unaffected) and assess the impact on Ly$\alpha$ visibility during reionization. We show that UV luminosity-halo mass scatter reduces Ly$\alpha$ visibility compared to models without scatter, and that this is most significant for UV-bright galaxies. We then use our model with scatter to infer the neutral fraction, $\overline{x}_{\mathrm{HI}}$, at $z \sim 7$ using a sample of Lyman-break galaxies in legacy fields. We infer $\overline{x}_{\mathrm{HI}} = 0.55_{-0.13}^{+0.11}$ with scatter, compared to $\overline{x}_{\mathrm{HI}} = 0.59_{-0.14}^{+0.12}$ without scatter, a very slight decrease and consistent within the uncertainties. Finally, we place our results in the context of other constraints on the reionization timeline and discuss implications for future high-redshift galaxy studies.

Lyman-$\alpha$ blobs (LABs) are spatially extended nebulae of emission in the Ly$\alpha$ line of hydrogen, seen at high redshifts$^{1,2}$, and most commonly found in the dense environment of star-forming galaxies$^{3,4}$. The origin of Ly$\alpha$ emission in the LABs is still unclear and under debate$^{5}$. Proposed powering sources generally fall into two categories: (1) photoionization, galactic super-winds/outflows, resonant scattering of Ly$\alpha$ photons from starbursts or active galactic nuclei (AGNs)$^{6,7,8,9,10}$ and (2) cooling radiation from cold streams of gas accreting onto galaxies$^{12}$. Here we analyze the gas kinematics within a LAB providing rare observational evidence for infalling gas. This is consistent with the release of gravitational accretion energy as cold streams radiate Ly$\alpha$ photons. It also provides direct evidence for possible cold streams feeding the central galaxies. The infalling gas is not important by mass but hints at more than one mechanism to explain the origin of the extended Ly$\alpha$ emission around young galaxies. It is also possible that the infalling gas may represent material falling back to the galaxy from where it originated, forming a galactic fountain.

We discovered a strongly lensed (\mu >40) Lya emission at z=6.629 (S/N~18) in the MUSE Deep Lensed Field (MDLF) targeting the Hubble Frontier Field galaxy cluster MACS~J0416. Dedicated lensing simulations imply that the Lya emitting region necessarily crosses the caustic. The arc-like shape of the Lya extends 3 arcsec on the observed plane and is the result of two merged multiple images, each one with a de-lensed Lya luminosity L<~2.8 x 10^(40) erg/s arising from a confined region (< 150 pc effective radius). A spatially unresolved HST counterpart is barely detected at S/N~2 after stacking the near-infrared bands, corresponding to an observed(intrinsic) magnitude m_(1500)>~30.8(>~35.0). The inferred rest-frame Lya equivalent width is EWo > 1120 A if the IGM transmission is T(IGM)<0.5. The low luminosities and the extremely large Lya EWo match the case of a Population~III star complex made of several dozens stars (~ 10^4 Msun) which irradiate a HII region crossing the caustic. While the Lya and stellar continuum are among the faintest ever observed at this redshift, the continuum and the Lya emissions could be affected by differential magnification, possibly biasing the EWo estimate. The aforementioned tentative HST detection tend to favor a large EWo, making such a faint Pop~III candidate a key target for the James Webb Space Telescope and Extremely Large Telescopes.

Recent near-IR power-spectra and panchromatic Extragalactic Background Light measurements provide upper limits on the near-IR surface brightness (SB>31 mag/arcsec^2) that may come from Pop III stars and accretion disks around resulting stellar-mass black holes (BHs) in the epoch of First Light (z=7-17). Physical parameters for zero metallicity Pop III stars at z>7 can be estimated from MESA stellar evolution models through helium-depletion, and for BH accretion disks from quasar microlensing results and multicolor accretion models. Second-generation stars can form at higher multiplicity, so that BH accretion disks may be fed by Roche-lobe overflow from lower-mass companions in their AGB stage. The near-IR SB constraints can be used to calculate the number of caustic transits behind lensing clusters that JWST and the 25~39 m ground-based telescopes may detect for both Pop III stars and stellar mass BH accretion disks. Because Pop III stars and stellar mass BH accretion disks have sizes of a few x 10^-11 arcsec at z>7, typical caustic magnifications can be mu=10^4~10^5, with rise times of hours and decline times of < 1 year for cluster transverse velocities of v<1000 km/s. Microlensing by intracluster medium objects can modify transit magnifications, and lengthen visibility times. Depending on BH masses, accretion-disk radii and feeding efficiencies, stellar-mass BH accretion-disk caustic transits could outnumber those from Pop III stars. To observe Pop III caustic transits directly may require monitoring 3~30 lensing clusters to AB< 29 mag over a decade or more. Such a program must be started with JWST in Cycle 1, and -- depending on the role of microlensing in the Intra Cluster Light -- should be continued for decades with the GMT and TMT, where JWST and the ground-based telescopes each will play a unique and strongly complementary role.

Galactic outflows are critical to our understanding of galaxy formation and evolution. However the details of the underlying feedback process remain unclear. We compare Ly$\alpha$ observations of the circumgalactic medium (CGM) of Lyman Break Galaxies (LBGs) with mock observations of their simulated CGM. We use cosmological hydrodynamical `zoom-in' simulations of an LBG which contains strong, momentum-driven galactic outflows. Simulation snapshots at $z=2.2$ and $z=2.65$ are used, corresponding to the available observational data. The simulation is post-processed with the radiative transfer code \textsc{crash} to account for the impact of ionising photons on hydrogen gas surrounding the simulated LBG. We generate mock absorption line maps for comparison with data derived from observed close galaxy-galaxy pairs. We perform calculations of Ly$\alpha$ photons scattering through the CGM with our newly developed Monte-Carlo code \textsc{slaf}, and compare to observations of diffuse Ly$\alpha$ halos around LBGs. Our fiducial galactic outflow model comes closer to reproducing currently observed characteristics of the CGM in Ly$\alpha$ than a reference inefficient feedback model used for comparison. Nevertheless, our fiducial model still struggles to reproduce the observed data of the inner CGM (at impact parameter $b<30$kpc). Our results suggest that galactic outflows affect Ly$\alpha$ absorption and emission around galaxies mostly at impact parameters $b<50$ kpc, while cold accretion flows dominate at larger distances. We discuss the implications of this result, and underline the potential constraining power of CGM observations - in emission and absorption - on galactic outflow models.

We present a new perturbative approach to "constrained Ly$\alpha$ radiative transfer'" (RT) through the circum- and inter-galactic medium (CGM and IGM). We constrain the HI content and kinematics of the CGM and IGM in a physically motivated model, using the galaxy-Ly$\alpha$ forest clustering data from spectroscopic galaxy surveys in QSO fields at $z\sim2-3$. This enables us to quantify the impact of the CGM/IGM on Ly$\alpha$ emission in an observationally constrained, realistic cosmological environment. Our model predicts that the CGM and IGM at these redshifts transmit $\approx80~\%$ of Ly$\alpha$ photons after having escaped from galaxies. This implies that while the inter-stellar medium primarily regulates Ly$\alpha$ escape, the CGM has a non-negligible impact on the observed Ly$\alpha$ line properties and the inferred Ly$\alpha$ escape fraction, even at $z\sim 2-3$. Ly$\alpha$ scattering in the CGM and IGM further introduces an environmental dependence in the (apparent) Ly$\alpha$ escape fraction, and the observed population of Ly$\alpha$ emitting galaxies: the CGM/IGM more strongly suppresses direct Ly$\alpha$ emission from galaxies in overdense regions in the Universe, and redistributes this emission into brighter Ly$\alpha$ haloes. The resulting mean surface brightness profile of the Ly$\alpha$ haloes is generally found to be a power-law $\propto r^{-2.4}$. Although our model still contains arbitrariness, our results demonstrate how (integral field) spectroscopic surveys of galaxies in QSO fields constrain circumgalactic Ly$\alpha$ RT, and we discuss the potential of these models for studying CGM physics and cosmology.

This chapter describes the prospects for detecting the first black holes in our Universe, with a specific focus on instruments/telescopes that will come online within the next decade, including e.g. the SKA, WFIRST, EUCLID, JWST, and large ground-based facilities such as E-ELTs, GMT and TMT. This chapter focusses on: (1) the indirect detectability of the first generations of stellar mass black holes through their imprint on various Cosmic Background Radiation fields including the 21-cm and X-ray backgrounds; (2) direct detectability of line and continuum emission by more massive black holes (M_BH ~ 1e4-1e6 M_sun) that formed via channels other than `ordinary' stellar evolution.

We present a detailed study of a unusually bright, lensed galaxy at $z=5.424$ discovered within the CFHTLS imaging survey. With an observed flux of $i_{\rm AB}=23.0$, J141446.82+544631.9 is one of the brightest galaxies known at $z>5$. It is characterized by strong Ly$\alpha$ emission, reaching a peak in (observed) flux density of $>10^{-16}~{\rm erg}~{\rm s}^{-1}~{\rm cm}^{-2}~{\rm \AA}^{-1}$. A deep optical spectrum from the LBT places strong constraints on NV and CIV emission, disfavouring an AGN source for the emission. However, a detection of the NIV]~$\lambda$1486 emission line indicates a hard ionizing continuum, possibly from hot, massive stars. Resolved imaging from HST deblends the galaxy from a foreground interloper; these observations include narrowband imaging of the Ly$\alpha$ emission, which is marginally resolved on $\sim$few~kpc scales and has EW$_0~\sim$ 260\AA. The Ly$\alpha$ emission extends over ~2000~${\rm km}~{\rm s}^{-1}$ and is broadly consistent with expanding shell models. SED fitting that includes Spitzer/IRAC photometry suggests a complex star formation history that include both a recent burst and an evolved population. J1414+5446 lies 30" from the centre of a known lensing cluster in the CFHTLS; combined with the foreground contribution this leads to a highly uncertain estimate for the lensing magnification in the range $5 \leq \mu \leq 25$. Because of its unusual brightness J1414+5446 affords unique opportunities for detailed study of an individual galaxy near the epoch of reionization and a preview of what can be expected from upcoming wide-area surveys that will yield hundreds of similar objects.

Recent detections of Lyman alpha (Ly$\alpha$) emission from $z>7.5$ galaxies were somewhat unexpected given a dearth of previous non-detections in this era when the intergalactic medium (IGM) is still highly neutral. But these detections were from UV bright galaxies, which preferentially live in overdensities which reionize early, and have significantly Doppler-shifted Ly$\alpha$ line profiles emerging from their interstellar media (ISM), making them less affected by the global IGM state. Using a combination of reionization simulations and empirical ISM models we show, as a result of these two effects, UV bright galaxies in overdensities have $>2\times$ higher transmission through the $z\sim7$ IGM than typical field galaxies, and this boosted transmission is enhanced as the neutral fraction increases. The boosted transmission is not sufficient to explain the observed high Ly$\alpha$ fraction of $M_\mathrm{UV} \lesssim -22$ galaxies (Stark et al. 2017), suggesting Ly$\alpha$ emitted by these galaxies must be stronger than expected due to enhanced production and/or selection effects. Despite the bias of UV bright galaxies to reside in overdensities we show Ly$\alpha$ observations of such galaxies can accurately measure the global neutral hydrogen fraction, particularly when Ly$\alpha$ from UV faint galaxies is extinguished, making them ideal candidates for spectroscopic follow-up into the cosmic Dark Ages.

Lyman-$\alpha$ (Ly$\alpha$) is a powerful astrophysical probe. Not only is it ubiquitous at high redshifts, it is also a resonant line, making Ly$\alpha$ photons scatter. This scattering process depends on the physical conditions of the gas through which Ly$\alpha$ propagates, and these conditions are imprinted on observables such as the Ly$\alpha$ spectrum and its surface brightness profile. In this work, we focus on a less-used observable capable of probing any scattering process: polarization. We implement the density matrix formalism of polarization into the Monte Carlo radiative transfer code tlac. This allows us to treat it as a quantum mechanical process where single photons develop and lose polarization from scatterings in arbitrary gas geometries. We explore static and expanding ellipsoids, biconical outflows, and clumpy multiphase media. We find that photons become increasingly polarized as they scatter and diffuse into the wings of the line profiles, making scattered Ly$\alpha$ polarized in general. The degree and orientation of Ly$\alpha$ polarization depends on the kinematics and distribution of the scattering HI gas. We find that it generally probes spatial or velocity space asymmetries and aligns itself tangentially to the emission source. We show that the mentioned observables, when studied separately, can leave similar signatures for different source models. We conclude by revealing how a joint analysis of the Ly$\alpha$ spectra, surface brightness profiles, and polarization can break these degeneracies and help us extract unique physical information on galaxies and their environments from their strongest, most prominent emission line.

We present a new flexible Bayesian framework for directly inferring the fraction of neutral hydrogen in the intergalactic medium (IGM) during the Epoch of Reionization (EoR, z~6-10) from detections and non-detections of Lyman Alpha (Ly$\alpha$) emission from Lyman break galaxies (LBGs). Our framework combines sophisticated reionization simulations with empirical models of the interstellar medium (ISM) radiative transfer effects on Ly$\alpha$. We assert that the Ly$\alpha$ line profile emerging from the ISM has an important impact on the resulting transmission of photons through the IGM, and that these line profiles depend on galaxy properties. We model this effect by considering the peak velocity offset of Ly$\alpha$ lines from host galaxies' systemic redshifts, which are empirically correlated with UV luminosity and redshift (or halo mass at fixed redshift). We use our framework on the sample of LBGs presented in Pentericci et al. (2014) and infer a global neutral fraction at z~7 of $\overline{x}_\mathrm{HI} = 0.59_{-0.15}^{+0.11}$, consistent with other robust probes of the EoR and confirming reionization is on-going ~700 Myr after the Big Bang. We show that using the full distribution of Ly$\alpha$ equivalent width detections and upper limits from LBGs places tighter constraints on the evolving IGM than the standard Ly$\alpha$ emitter fraction, and that larger samples are within reach of deep spectroscopic surveys of gravitationally lensed fields and JWST NIRSpec.

We summarize panchromatic Extragalactic Background Light data to place upper limits on the integrated near-infrared surface brightness (SB) that may come from Population III stars and possible accretion disks around their stellar-mass black holes (BHs) in the epoch of First Light, broadly taken from z$\simeq$7-17. Theoretical predictions and recent near-infrared power-spectra provide tighter constraints on their sky-signal. We outline the physical properties of zero metallicity Population III stars from MESA stellar evolution models through helium-depletion and of BH accretion disks at z$\gtrsim$7. We assume that second-generation non-zero metallicity stars can form at higher multiplicity, so that BH accretion disks may be fed by Roche-lobe overflow from lower-mass companions. We use these near-infrared SB constraints to calculate the number of caustic transits behind lensing clusters that the James Webb Space Telescope and the next generation ground-based telescopes may observe for both Population III stars and their BH accretion disks. Typical caustic magnifications can be $\mu$$\simeq$10$^4$-10$^5$, with rise times of hours and decline times of $\lesssim$1 year for cluster transverse velocities of $v_{T}$$\lesssim$1000 km s$^{-1}$. Microlensing by intracluster medium objects can modify transit magnifications, but lengthen visibility times. Depending on BH masses, accretion-disk radii and feeding efficiencies, stellar-mass BH accretion-disk caustic transits could outnumber those from Population III stars. To observe Population III caustic transits directly may require to monitor 3-30 lensing clusters to AB$\lesssim$29 mag over a decade.

We present observations of escaping Lyman Continuum (LyC) radiation from 34 massive star-forming galaxies and 12 weak AGN with reliably measured spectroscopic redshifts at $z$$\simeq$2.3-4.1. We analyzed Hubble Space Telescope Wide Field Camera 3 (WFC3) mosaics of the Early Release Science field in three UVIS filters to sample the rest-frame LyC over this redshift range. With our best current assessment of the WFC3 systematics, we provide 1$\sigma$ upper limits for the average LyC emission of galaxies at $z$=2.35, 2.75, and 3.60 to $\sim$28.5, 28.1, and 30.7 mag in image stacks of 11-15 galaxies in the WFC3/UVIS F225W, F275W, and F336W, respectively. The LyC flux of weak AGN at $z$=2.62 and 3.32 are detected at 28.3 and 27.4 mag with SNRs of $\sim$2.7 and 2.5 in F275W and F336W for stacks of 7 and 3 AGN, respectively, while AGN at $z$=2.37 are constrained to $\gtrsim$27.9 mag at 1$\sigma$ in a stack of 2 AGN. The stacked AGN LyC light profiles are flatter than their corresponding non-ionizing UV continuum profiles out to radii of r$\lesssim$0."9, which may indicate a radial dependence of porosity in the ISM. With synthetic stellar SEDs fit to UV continuum measurements longwards of Ly$\alpha$ and IGM transmission models, we constrain the absolute LyC escape fractions to $f_{\rm esc}^{\rm abs}$$\simeq$$22^{+44}_{-22}$% at $z$=2.35 and $\lesssim$55% at $z$=2.75 and 3.60, respectively. All available data for galaxies, including published work, suggests a more sudden increase of $f_{\rm esc}$ with redshift at $z$$\simeq$2. Dust accumulating in (massive) galaxies over cosmic time correlates with increased HI column density, which may lead to reducing $f_{\rm esc}$ more suddenly at $z$$\lesssim$2. This may suggest that star-forming galaxies collectively contributed to maintaining cosmic reionization at redshifts $z$$\gtrsim$2-4, while AGN likely dominated reionization at $z$$\lesssim$2.

Star-forming Compact Dwarf Galaxies (CDGs) resemble the expected pristine conditions of the first galaxies in the Universe and are the best systems to test models on primordial galaxy formation and evolution. Here we report on one of such CDGs, Tololo 1214-277, which presents a broad, single peaked, highly symmetric Ly$\alpha$ emission line that had evaded theoretical interpretation so far. In this paper we reproduce for the first time these line features with two different physically motivated kinematic models: an interstellar medium composed by outflowing clumps with random motions and an homogeneous gaseous sphere undergoing solid body rotation. The multiphase model requires a clump velocity dispersion of $54.3\pm 0.6$ km s$^{-1}$ with outflows of $54.3\pm 5.1$ km s$^{-1}$, while the bulk rotation velocity is constrained to be $348^{+75}_{-48}$ km s$^{-1}$. We argue that the results from the multiphase model provide a correct interpretation of the data. In that case the clump velocity dispersion implies a dynamical mass of $2\times 10^{9}$ M$_{\odot}$, ten times its baryonic mass. If future kinematic maps of Tololo 1214-277 confirm the velocities suggested by the multiphase model, it would provide additional support to expect such kinematic state in primordial galaxies, opening the opportunity to use the models and methods presented in this paper to constrain the physics of star formation and feedback in the early generation of Ly-$\alpha$ emitting galaxies.

We studied Lyman-$\alpha$ (Ly$\alpha$) escape in a statistical sample of 43 Green Peas with HST/COS Ly$\alpha$ spectra. Green Peas are nearby star-forming galaxies with strong [OIII]$\lambda$5007 emission lines. Our sample is four times larger than the previous sample and covers a much more complete range of Green Pea properties. We found that about 2/3 of Green Peas are strong Ly$\alpha$ line emitters with rest-frame Ly$\alpha$ equivalent width $>20$ \AA. The Ly$\alpha$ profiles of Green Peas are diverse. The Ly$\alpha$ escape fraction, defined as the ratio of observed Ly$\alpha$ flux to intrinsic Ly$\alpha$ flux, shows anti-correlations with a few Ly$\alpha$ kinematic features -- both the blue peak and red peak velocities, the peak separations, and FWHM of the red portion of the Ly$\alpha$ profile. Using properties measured from SDSS optical spectra, we found many correlations -- Ly$\alpha$ escape fraction generally increases at lower dust reddening, lower metallicity, lower stellar mass, and higher [OIII]/[OII] ratio. We fit their Ly$\alpha$ profiles with the HI shell radiative transfer model and found Ly$\alpha$ escape fraction anti-correlates with the best-fit $N_{HI}$. Finally, we fit an empirical linear relation to predict Ly$\alpha$ escape fraction from the dust extinction and Ly$\alpha$ red peak velocity. The standard deviation of this relation is about 0.3 dex. This relation can be used to isolate the effect of IGM scatterings from Ly$\alpha$ escape and to probe the IGM optical depth along the line of sight of each $z>7$ Ly$\alpha$ emission line galaxy in the JWST era.

We present a new method to quantify the value of the escape fraction of ionizing photons, and the existence of ultra-faint galaxies clustered around brighter objects during the epoch of cosmic reionization, using the diffuse Ly$\alpha$, continuum and H$\alpha$ emission observed around galaxies at $z\sim6$. We model the surface brightness profiles of the diffuse halos considering the fluorescent emission powered by ionizing photons escaping from the central galaxies, and the nebular emission from satellite star-forming sources, by extending the formalisms developed in Mas-Ribas & Dijkstra (2016) and Mas-Ribas et al. (2017). The comparison between our predicted profiles and Ly$\alpha$ observations at $z=5.7$ and $z=6.6$ favors a low ionizing escape fraction, $f_{\rm esc}^{\rm ion}\sim5\%$, for galaxies in the range $-19\gtrsim M_{\rm UV} \gtrsim -21.5$. However, uncertainties and possible systematics in the observations do not allow for firm conclusions. We predict H$\alpha$ and rest-frame visible continuum observations with JWST, and show that JWST will be able to detect extended (a few tens of kpc) fluorescent H$\alpha$ emission powered by ionizing photons escaping from a bright, $L\gtrsim 5L^*$, galaxy. Such observations can differentiate fluorescent emission from nebular emission by satellite sources. We discuss how observations and stacking of several objects may provide unique constraints on the escape fraction for faint galaxies and/or the abundance of ultra-faint radiation sources.

We study the spectrophotometric properties of a highly magnified (\mu~40-70) pair of stellar systems identified at z=3.2222 behind the Hubble Frontier Field galaxy cluster MACS~J0416. Five multiple images (out of six) have been spectroscopically confirmed by means of VLT/MUSE and VLT/X-Shooter observations. Each image includes two faint (m_uv~30.6), young (<100 Myr), low-mass (<10^7 Msun), low-metallicity (12+Log(O/H)~7.7, or 1/10 solar) and compact (30 pc effective radius) stellar systems separated by ~300pc, after correcting for lensing amplification. We measured several rest-frame ultraviolet and optical narrow (\sigma_v <~ 25 km/s) high-ionization lines. These features may be the signature of very hot (T>50000 K) stars within dense stellar clusters, whose dynamical mass is likely dominated by the stellar component. Remarkably, the ultraviolet metal lines are not accompanied by Lya emission (e.g., CIV / Lya > 15), despite the fact that the Lya line flux is expected to be 150 times brighter (inferred from the Hbeta flux). A spatially-offset, strongly-magnified (\mu>50) Lya emission with a spatial extent <~7.6 kpc^2 is instead identified 2 kpc away from the system. The origin of such a faint emission can be the result of fluorescent Lya induced by a transverse leakage of ionizing radiation emerging from the stellar systems and/or can be associated to an underlying and barely detected object (with m_uv > 34 de-lensed). This is the first confirmed metal-line emitter at such low-luminosity and redshift without Lya emission, suggesting that, at least in some cases, a non-uniform covering factor of the neutral gas might hamper the Lya detection.

Lyman alpha halos are observed ubiquitously around star-forming galaxies at high redshift, but their origin is still a matter of debate. We demonstrate that the emission from faint unresolved satellite sources, $M_{\rm UV} \gtrsim -17$, clustered around the central galaxies may play a major role in generating spatially extended Ly$\alpha$, continuum (${\rm UV + VIS}$) and H$\alpha$ halos. We apply the analytic formalism developed in Mas-Ribas & Dijkstra (2016) to model the halos around Lyman Alpha Emitters (LAEs) at $z=3.1$, for several different satellite clustering prescriptions. In general, our UV and Ly$\alpha$ surface brightness profiles match the observations well at $20\lesssim r \lesssim 40$ physical kpc from the centers of LAEs. We discuss how our profiles depend on various model assumptions and how these can be tested and constrained with future H$\alpha$ observations by the James Webb Space Telescope (JWST). Our analysis shows how spatially extended halos constrain (i) the presence of otherwise undetectable satellite sources, (ii) the integrated, volumetric production rates of Ly$\alpha$ and LyC photons, and (iii) their population-averaged escape fractions. These quantities are all directly relevant for understanding galaxy formation and evolution and, for high enough redshifts, cosmic reionization.

Motivated by observational and theoretical work which both suggest very small scale ($\lesssim 1\,$pc) structure in the circum-galactic medium of galaxies and in other environments, we study Lyman-$\alpha$ (Ly$\alpha$) radiative transfer in an extremely clumpy medium with many "clouds" of neutral gas along the line of sight. While previous studies have typically considered radiative transfer through sightlines intercepting $\lesssim 10$ clumps, we explore the limit of a very large number of clumps per sightline (up to $f_{\mathrm{c}} \sim 1000$). Our main finding is that, for covering factors greater than some critical threshold, a multiphase medium behaves similar to a homogeneous medium in terms of the emergent Ly$\alpha$ spectrum. The value of this threshold depends on both the clump column density and on the movement of the clumps. We estimate this threshold analytically and compare our findings to radiative transfer simulations with a range of covering factors, clump column densities, radii, and motions. Our results suggest that (i) the success in fitting observed Ly$\alpha$ spectra using homogeneous "shell models" (and the corresponding failure of multiphase models) hints towards the presence of very small-scale structure in neutral gas, in agreement within a number of other observations; and (ii) the recurrent problems of reproducing realistic line profiles from hydrodynamical simulations may be due to their inability to resolve small-scale structure, which causes simulations to underestimate the effective covering factor of neutral gas clouds.

Lecture notes for 8 lectures on the `Physics of Lyman alpha Radiative Transfer', given at the 46th Saas-Fee winter school held in Les Diablerets, Switzerland on March 13-19 2016. These lectures aimed at offering basic insights into Lyman alpha (Lya) radiative processes including emission processes and Lya radiative transfer, and highlighting some of the physics associated with these processes. The notes include derivations in greater detail than what was discussed during the lectures.

We assess the impact of trapped Lyman alpha cooling radiation on the formation of direct collapse black holes (DCBHs). We apply a one-zone chemical and thermal evolution model, accounting for the photodetachment of H- ions, precursors to the key coolant H2, by Lyman alpha photons produced during the collapse of a cloud of primordial gas in an atomic cooling halo at high redshift. We find that photodetachment of H- by trapped Lyman alpha photons may lower the level of the H2-dissociating background radiation field required for DCBH formation substantially, dropping the critical flux by up to a factor of a few. This translates into a potentially large increase in the expected number density of DCBHs in the early Universe, and supports the view that DCBHs may be the seeds for the BHs residing in the centers of a significant fraction of galaxies today. We find that detachment of H- by Lyman alpha has the strongest impact on the critical flux for the relatively high background radiation temperatures expected to characterize the emission from young, hot stars in the early Universe. This lends support to the DCBH origin of the highest redshift quasars.

We present the results of ALMA spectroscopic follow-up of a $z=6.765$ Lyman-$\alpha$ emitting galaxy behind the cluster RXJ1347-1145. We report the detection of [CII]158$\mu$m line fully consistent with the Lyman-$\alpha$ redshift and with the peak of the optical emission. Given the magnification of $\mu=5.0 \pm 0.3$ the intrinsic (corrected for lensing) luminosity of the [CII] line is $L_{[CII]} =1.4^{+0.2}_{-0.3} \times 10^7L_{\odot}$, which is ${\sim}5$ times fainter than other detections of $z\sim 7$ galaxies. The result indicates that low $L_{[CII]}$ in $z\sim 7$ galaxies compared to the local counterparts might be caused by their low metallicities and/or feedback. The small velocity off-set ($\Delta v = 20_{-40}^{+140} \rm km/s$) between the Lyman-$\alpha$ and [CII] line is unusual, and may be indicative of ionizing photons escaping.

Lyman-Alpha (Ly$\alpha$) is the strongest emission line in the Universe and is frequently used to detect and study the most distant galaxies. Because Lya is a resonant line, photons typically scatter prior to escaping; this scattering process complicates the interpretation of Ly$\alpha$ spectra, but also encodes a wealth of information about the structure and kinematics of neutral gas in the galaxy. Modeling the Ly$\alpha$ line therefore allows us to study tiny-scale features of the gas, even in the most distant galaxies. Curiously, observed Ly$\alpha$ spectra can be modeled successfully with very simple, homogeneous geometries (such as an expanding, spherical shell), whereas more realistic, multiphase geometries often fail to reproduce the observed spectra. This seems paradoxical since the gas in galaxies is known to be multiphase. In this Letter, we show that spectra emerging from extremely clumpy geometries with many clouds along the line of sight converge to the predictions from simplified, homogeneous models. We suggest that this resolves the apparent discrepancy, and may provide a way to study the gas structure in galaxies on scales far smaller than can be probed in either cosmological simulations or direct (i.e., spatially-resolved) observations.

We revisit calculations of nebular hydrogen Lya and HeII 1640A line strengths for population III galaxies, undergoing continuous and bursts of star formation. We focus on initial mass functions (IMFs) motivated by recent theoretical studies, which generally span a lower range of stellar masses than earlier works. We also account for case-B departures and the stochastic sampling of the IMF. In agreement with previous works, we find that departures from case-B can enhance the Lya flux by a factor of a few, but we argue that this enhancement is driven mainly by collisional excitation and ionization, and not due to photoionization from the n = 2 state of atomic hydrogen. The increased sensitivity of the Lya flux to the high-energy end of the galaxy spectrum makes it more subject to stochastic sampling of the IMF. The latter introduces a dispersion in the predicted nebular line fluxes around the deterministic value by as much as a factor of ~4. In contrast, the stochastic sampling of the IMF has less impact on the emerging Lyman Werner (LW) photon flux. When case-B departures and stochasticity effects are combined, nebular line emission from population III galaxies can be up to one order of magnitude brighter than predicted by 'standard' calculations that do not include these effects. This enhances the prospects for detection with future facilities such as JWST and large, groundbased telescopes.

In spite of their conjectured importance for the Epoch of Reionization, the properties of low-mass galaxies are currently still under large debate. In this article, we study the stellar and gaseous properties of faint, low-mass galaxies at z>3. We observed the Frontier Fields cluster Abell S1063 with MUSE over a 2 arcmin^2 field, and combined integral-field spectroscopy with gravitational lensing to perform a blind search for intrinsically faint Lya emitters (LAEs). We determined in total the redshift of 172 galaxies of which 14 are lensed LAEs at z=3-6.1. We increased the number of spectroscopically-confirmed multiple-image families from 6 to 17 and updated our gravitational-lensing model accordingly. The lensing-corrected Lya luminosities are with L(Lya) <= 10^41.5 erg/s among the lowest for spectroscopically confirmed LAEs at any redshift. We used expanding gaseous shell models to fit the Lya line profile, and find low column densities and expansion velocities. This is to our knowledge the first time that gaseous properties of such faint galaxies at z>=3 are reported. We performed SED modelling to broadband photometry from the U-band through the infrared to determine the stellar properties of these LAEs. The stellar masses are very low (10^{6-8} Msun), and are accompanied by very young ages of 1-100 Myr. The very high specific star formation rates (~100 Gyr^-1) are characteristic of starburst galaxies. The UV-continuum slopes beta are low in our sample, with beta<-2 for all galaxies with M_star<10^8 Msun. We conclude that our low-mass galaxies at 3<z<6 are forming stars at higher rates when correcting for stellar mass effects than seen locally or in more massive galaxies. The young stellar populations with high star-formation rates and low HI column densities lead to continuum slopes and LyC-escape fractions expected for a scenario where low mass galaxies reionise the Universe.

We present new Atacama Large Millimeter/Submillimeter Array (ALMA) 850um continuum observations of the original Lyman-alpha Blob (LAB) in the SSA22 field at z=3.1 (SSA22-LAB01). The ALMA map resolves the previously identified submillimeter source into three components with total flux density S_850 = 1.68+/-0.06 mJy, corresponding to a star formation rate of ~150 M_sun/yr. The submillimeter sources are associated with several faint (m~27 mag) rest-frame ultraviolet sources identified in Hubble Space Telescope Imaging Spectrograph (STIS) clear filter imaging (~5850A). One of these companions is spectroscopically confirmed with Keck MOSFIRE to lie within 20 projected kpc and 250 km/s of one of the ALMA components. We postulate that some of these STIS sources represent a population of low-mass star-forming satellites surrounding the central submillimeter sources, potentially contributing to their growth and activity through accretion. Using a high resolution cosmological zoom simulation of a 10^13 M_sun halo at z=3, including stellar, dust and Ly-alpha radiative transfer, we can model the ALMA+STIS observations and demonstrate that Ly-alpha photons escaping from the central submillimeter sources are expected to resonantly scatter in neutral hydrogen, the majority of which is predicted to be associated with halo substructure. We show how this process gives rise to extended Ly-alpha emission with similar surface brightness and morphology to observed giant LABs.

The escape of ionizing Lyman Continuum (LyC) photons requires the existence of low-N_HI sightlines, which also promote escape of Lyman-Alpha (Lya). We use a suite of 2500 Lya Monte-Carlo radiative transfer simulations through models of dusty, clumpy interstellar (`multiphase') media from Gronke & Dijkstra (2016), and compare the escape fractions of Lya [f_esc(Lya)] and LyC radiation [f_esc(LyC)]. We find that f_esc(LyC) and f_esc(Lya) are correlated: galaxies with a low f_esc(Lya) consistently have a low f_esc(LyC), while galaxies with a high f_esc(Lya) exhibit a large dispersion in f_esc(LyC). We argue that there is increasing observational evidence that Lya escapes more easily from UV-faint galaxies. The correlation between f_esc(LyC) and f_esc(Lya) then implies that UV-faint galaxies contribute more to the ionizing background than implied by the faint-end slope of the UV-luminosity function. In multiphase gases, the ionizing escape fraction is most strongly affected by the cloud covering factor, f_cl, which implies that f_esc(LyC) is closely connected to the observed Lya spectral line shape. Specifically, LyC emitting galaxies typically having narrower, more symmetric line profiles. This prediction is qualitatively similar to that for `shell models'.

We report on the discovery of extended Ly-alpha nebulae at z~3.3 in the Hubble Ultra Deep Field (HUDF, ~ 40 kpc X 80 kpc) and behind the Hubble Frontier Fields galaxy cluster MACSJ0416 (~ 40kpc), spatially associated with groups of star-forming galaxies. VLT/MUSE integral field spectroscopy reveals a complex structure with a spatially-varying double peaked Ly-alpha emission. Overall, the spectral profiles of the two Ly-alpha nebulae are remarkably similar, both showing a prominent blue emission, more intense and slightly broader than the red peak. From the first nebula, located in the HUDF, no X-ray emission has been detected, disfavoring the possible presence of AGNs. Spectroscopic redshifts have been derived for 11 galaxies within two arcsec from the nebula and spanning the redshift range 1.037<z<5.97. The second nebula, behind MACSJ0416, shows three aligned star-forming galaxies plausibly associated to the emitting gas. In both systems, the associated galaxies reveal possible intense rest-frame-optical nebular emissions lines [OIII]4959-5007+Hbeta with equivalent widths as high as 1500A rest-frame and star formation rates ranging from a few to tens of solar masses per year. A possible scenario is that of a group of young, star-forming galaxies sources of escaping ionising radiation that induce Ly-alpha fluorescence, therefore revealing the kinematics of the surrounding gas. Also Ly-alpha powered by star-formation and/or cooling radiation may resemble the double peaked spectral properties and the morphology observed here. If the intense blue emission is associated with inflowing gas, then we may be witnessing an early phase of galaxy or a proto-cluster (or group) formation.

We report the detection of Ly$\alpha$ emission at $\sim9538$\AA{} in the Keck/DEIMOS and \HST WFC3 G102 grism data from a triply-imaged galaxy at $z=6.846\pm0.001$ behind galaxy cluster MACS J2129.4$-$0741. Combining the emission line wavelength with broadband photometry, line ratio upper limits, and lens modeling, we rule out the scenario that this emission line is \oii at $z=1.57$. After accounting for magnification, we calculate the weighted average of the intrinsic Ly$\alpha$ luminosity to be $\sim1.3\times10^{42}~\mathrm{erg}~\mathrm{s}^{-1}$ and Ly$\alpha$ equivalent width to be $74\pm15$\AA{}. Its intrinsic UV absolute magnitude at 1600\AA{} is $-18.6\pm0.2$ mag and stellar mass $(1.5\pm0.3)\times10^{7}~M_{\odot}$, making it one of the faintest (intrinsic $L_{UV}\sim0.14~L_{UV}^*$) galaxies with Ly$\alpha$ detection at $z\sim7$ to date. Its stellar mass is in the typical range for the galaxies thought to dominate the reionization photon budget at $z\gtrsim7$; the inferred Ly$\alpha$ escape fraction is high ($\gtrsim 10$\%), which could be common for sub-$L^*$ $z\gtrsim7$ galaxies with Ly$\alpha$ emission. This galaxy offers a glimpse of the galaxy population that is thought to drive reionization, and it shows that gravitational lensing is an important avenue to probe the sub-$L^*$ galaxy population.

We perform Lyman-$\alpha$ (Ly$\alpha$) Monte-Carlo radiative transfer calculations on a suite of $2500$ models of multiphase, outflowing media, which are characterized by $14$ parameters. We focus on the Ly$\alpha$ spectra emerging from these media, and investigate which properties are dominant in shaping the emerging Ly$\alpha$ profile. Multiphase models give rise to a wide variety of emerging spectra, including single, double and triple peaked spectra. We find that the dominant parameters in shaping the spectra include (i) the cloud covering factor, $f_c$, in agreement with earlier studies, and (ii) the temperature and number density of residual HI in the hot ionized medium. We attempt to reproduce spectra emerging from multiphase models with `shell models' which are commonly used to fit observed Ly$\alpha$ spectra, and investigate the connection between shell-model parameters and the physical parameters of the clumpy media. In shell models, the neutral hydrogen content of the shell is one of the key parameters controlling Ly$\alpha$ radiative transfer. Because Ly$\alpha$ spectra emerging from multi-phase media depend much less on the neutral hydrogen content of the clumps, the shell model parameters such as HI column density (but also shell velocity and dust content) are generally not well matched to the associated physical parameters of clumpy media.

We present VLT/X-Shooter and MUSE spectroscopy of an faint F814W=28.60+/-0.33 (Muv=-17.0), low mass (~<10^7 Msun) and compact (Reff=62pc) freshly star-forming galaxy at z=3.1169 magnified (16x) by the Hubble Frontier Fields galaxy cluster Abell S1063. Gravitational lensing allows for a significant jump toward low-luminosity regimes, in moderately high resolution spectroscopy (R=lambda/dlambda ~ 3000-7400). We measured CIV1548,1550, HeII1640, OIII]1661,1666, CIII]1907,1909, Hbeta, [OIII]4959,5007, emission lines with FWHM< 50 km/s and (de-lensed) fluxes spanning the interval 1.0x10^-19 - 2.0x10^-18 erg/s/cm2 at S/N=4-30. The double peaked Lya emission with Delta_v(red-blue) = 280(+/-7)km/s and de-lensed fluxes 2.4_(blue)|8.5_(red)x10^-18 erg/s/cm2 (S/N=38_(blue)|110_(red)) indicate a low column density of neutral hydrogen gas consistent with a highly ionized interstellar medium as also inferred from the large [OIII]5007/[OII]3727>10 ratio. We detect CIV1548,1550 resonant doublet in emission, each component with FWHM ~< 45 km/s, and redshifted by +51(+/-10)km/s relative to the systemic redshift. We interpret this as nebular emission tracing an expanding optically-thin interstellar medium. Both CIV1548,1550 and HeII1640 suggest the presence of hot and massive stars (with a possible faint AGN). The ultraviolet slope is remarkably blue, beta =-2.95 +/- 0.20 (F_lambda=lambda^beta), consistent with a dust-free and young ~<20 Myr galaxy. Line ratios suggest an oxygen abundance 12+log(O/H)<7.8. We are witnessing an early episode of star-formation in which a relatively low NHI and negligible dust attenuation might favor a leakage of ionizing radiation. This galaxy currently represents a unique low-luminosity reference object for future studies of the reionization epoch with JWST.

We quantify the contribution of Ly$\alpha$ fluorescence to observed spatially extended Ly$\alpha$ halos around Ly$\alpha$ emitters (LAE) at redshift ${\rm z=3.1}$. The key physical quantities that describe the fluorescent signal include (${\it i}$) the distribution of cold gas in the circum-galactic medium (CGM); we explore simple analytical models and fitting functions to recent hydrodynamical simulations; (${\it ii}$) local variations in the ionizing background due to ionizing sources that cluster around the central galaxy. We account for clustering by boosting the observationally inferred volumetric production rate of ionizing photons, $\epsilon_{\rm LyC}$, by a factor of $1+\xi_{\rm LyC}(r)$, in which $\xi_{\rm LyC}(r)$ quantifies the clustering of ionizing sources around the central galaxy. We compute $\xi_{\rm LyC}(r)$ by assigning an 'effective' bias parameter to the ionizing sources. This novel approach allows us to quantify our ignorance of the population of ionizing sources in a simple parametrized form. We find a maximum enhancement in the local ionizing background in the range $50-200$ at $r \sim 10$ physical kpc. For spatially uncorrelated ionizing sources and fluorescing clouds we find that fluorescence can contribute up to $\sim 50-60\%$ of the observed spatially extended Ly$\alpha$ emission. We briefly discuss how future observations can shed light on the nature of Ly$\alpha$ halos around star forming galaxies.

We analyze archival Ly$\alpha$ spectra of 12 "Green Pea" galaxies observed with the Hubble Space Telescope, model their Ly$\alpha$ profiles with radiative transfer models, and explore the dependence of Ly$\alpha$ escape fraction on various properties. Green Pea galaxies are nearby compact starburst galaxies with [OIII]$\lambda$5007 equivalent widths of hundreds of \AA. All 12 Green Pea galaxies in our sample show Ly$\alpha$ lines in emission, with a Ly$\alpha$ equivalent width distribution similar to high redshift Ly$\alpha$ emitters. Combining the optical and UV spectra of Green Pea galaxies, we estimate their Ly$\alpha$ escape fractions and find correlations between Ly$\alpha$ escape fraction and kinematic features of Ly$\alpha$ profiles. The escape fraction of Ly$\alpha$ in these galaxies ranges from 1.4% to 67%. We also find that the Ly$\alpha$ escape fraction depends strongly on metallicity and moderately on dust extinction. We compare their high-quality Ly$\alpha$ profiles with single HI shell radiative transfer models and find that the Ly$\alpha$ escape fraction anti-correlates with the derived HI column densities. Single shell models fit most Ly$\alpha$ profiles well, but not the ones with highest escape fractions of Ly$\alpha$. Our results suggest that low HI column density and low metallicity are essential for Ly$\alpha$ escape, and make a galaxy a Ly$\alpha$ emitter.

`Direct collapse black holes' (DCBHs) provide possible seeds for supermassive black holes that exist at redshifts as high as z~7. We study Lyman Alpha (Lya) radiative transfer through simplified representations of the DCBH-scenario. We find that gravitational heating of the collapsing cloud gives rise to a Lya cooling luminosity of up to ~ 1e38(M_gas/1e6 Msun)^2 erg/s. The Lya production rate can be significantly larger during the final stages of collapse, but collisional deexcitation efficiently suppresses the emerging Lya flux. Photoionization by a central source boosts the Lya luminosity to L~1e43(M_BH/1e6 M_sun) erg/s during specific evolutionary stages of the cloud, where M_BH denotes the mass of the black hole powering this source. We predict that the width and velocity off-set of the Lya spectral line range from a few tens to few thousands km/s, depending sensitively on the evolutionary state of the cloud. We also compare our predictions to observations of CR7 (Sobral et al. 2015), a luminous Lya emitter at z~7, which is potentially associated with a DCBH. If CR7 is powered by a black hole, then its Lya flux alone requires that M_BH> 1e7 M_sun, which exceeds the mass of DCBHs when they first form. The observed width of the Lya spectrum favors the presence of only a low column density of hydrogen, log [N_HI/cm^-2]~19-20. The shape of the Lya spectrum indicates that this gas is outflowing. These requirements imply that if CR7 harbors a DCBH, then the physical conditions that enabled its formation have been mostly erased, which is in agreement with theoretical expectations.

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-$\alpha$ "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$\alpha$ 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-$\alpha$ 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$\alpha$ 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.

The direct collapse black hole (DCBH) scenario describes the isothermal collapse of a pristine gas cloud directly into a massive, M_BH=10^4-10^6 M_sun black hole. In this paper we show that large HI column densities of primordial gas at T~10^4 K with low molecular abundance - which represent key aspects of the DCBH scenario - provide optimal conditions for pumping of the 2p-level of atomic hydrogen by trapped Lyman alpha (Lya) photons. This Lya pumping mechanism gives rise to inverted level population of the 2s_1/2-2p_3/2 transition, and therefore to stimulated fine structure emission at 3.04 cm (rest-frame). We show that simplified models of the DCBH scenario amplify the CMB by up to a factor of 10^5, above which the maser saturates. Hyperfine splitting of the 3-cm transition gives rise to a characteristic broad (FWHM ~ tens of MHz in the observers frame) asymmetric line profile. This signal subtends an angular scale of ~ 1-10 mas, which translates to a flux of ~ 0.3-3 microJy, which is detectable with ultra-deep surveys being planned with SKA1-MID. While challenging, as the signal is visible for a fraction of the collapse time of the cloud, the matching required physical conditions imply that a detection of the redshifted 3-cm emission line would provide direct evidence for the DCBH scenario.

[abbreviated] We present a census of Ly\alpha\ emission at $z\gtrsim7$ utilizing deep near infrared HST grism spectroscopy from the first six completed clusters of the Grism Lens-Amplified Survey from Space (GLASS). In 24/159 photometrically selected galaxies we detect emission lines consistent with Ly\alpha\ in the GLASS spectra. Based on the distribution of signal-to-noise ratios and on simulations we expect the completeness and the purity of the sample to be 40-100% and 60-90%, respectively. For the objects without detected emission lines we show that the observed (not corrected for lensing magnification) 1$\sigma$ flux limits reaches $5\times10^{-18}$erg/s/cm$^{2}$ per position angle over the full wavelength range of GLASS (0.8-1.7$\mu$m). Based on the conditional probability of Ly\alpha\ emission measured from the ground at $z\sim7$ we would have expected 12-18 Ly\alpha\ emitters. This is consistent with the number of detections, within the uncertainties, confirming the drop in Ly\alpha\ emission with respect to $z\sim6$. These candidates include a promising source at $z=8.1$. The spatial extent of Ly\alpha\ in a deep stack of the most convincing Ly\alpha\ emitters with $\langle z\rangle=7.2$ is consistent with that of the rest-frame UV continuum. Extended Ly$\alpha$ emission, if present, has a surface brightness below our detection limit, consistent with the properties of lower redshift comparison samples. From the stack we estimate upper limits on rest-frame UV emission line ratios and find $f_\textrm{CIV} / f_\textrm{Ly${\alpha}$} \lesssim 0.32$ and $f_\textrm{CIII]} / f_\textrm{Ly$\alpha$} \lesssim 0.23$ in good agreement with other values published in the literature.

Neutral diffuse intergalactic gas that existed during the Epoch of Reionization (EoR) suppresses Lyman Alpha (Lya) flux emitted by background galaxies. In this chapter I summarise the increasing observational support for the claim that Lya photons emitted by galaxies at z>6 are suppressed by intervening HI gas. I describe key physical processes that affect Lya transfer during the EoR. I argue that in spite of the uncertainties associated with this complex multiscale problem, the data on Lya emitting galaxies at z=0-6 strongly suggests that the observed reduction in Lya flux from galaxies at z>6 is due to additional intervening HI gas. The main question is what fraction of this additional HI gas is in the diffuse neutral IGM. I summarise how future surveys on existing and upcoming instruments are expected to reduce existing observational uncertainties enormously. With these improved data we will likely be able to nail down reionization with Lya emitting galaxies.

The visibility of LyA emitting galaxies during the Epoch of Reionization is controlled by both diffuse HI patches in large-scale bubble morphology and small-scale absorbers. To investigate the impact on LyA photons, we apply a novel combination of analytic and numerical calculations to three scenarios: (i) the `bubble' model, where only diffuse HI outside ionized bubbles is present; (ii) the `web' model, where HI exists only in overdense self-shielded gas; and (iii) the more realistic 'web-bubble' model, which contains both. Our analysis confirms that there is a degeneracy between the ionization structure of the intergalactic medium (IGM) and the HI fraction inferred from LyA surveys, as the three models suppress LyA flux equally with very different HI fractions. We argue that a joint analysis of the LyA luminosity function and the rest-frame equivalent width distribution/LyA fraction can break this degeneracy and provide constraints on the reionization history and its topology. We further show that constraints can improve if we consider the full shape of the M_UV-dependent redshift evolution of the LyA fraction of Lyman break galaxies. Contrary to conventional wisdom, we find that (i) a drop of LyA fraction larger for UV-faint than for UV-bright galaxies can be reproduced with web and web-bubble models and therefore does not provide exclusive evidence of patchy reionization, and (ii) the IGM-transmission PDF is unimodal for bubble models and bimodal in web models. We further highlight the importance of galaxy-absorber cross-correlation. Comparing our models to observations, the neutral fraction at z~7 is likely to be of order of tens of per cent when interpreted with bubble or web-bubble models. Alternatively, we obtain a conservative lower limit ~1% in the web models, if we allow for a drop in the photoionization rate by a factor of ~100 from the post-reionized universe. [abridged]

Outflows promote the escape of Lyman-$\alpha$ (Ly$\alpha$) photons from dusty interstellar media. The process of radiative transfer through interstellar outflows is often modelled by a spherically symmetric, geometrically thin shell of gas that scatters photons emitted by a central Ly$\alpha$ source. Despite its simplified geometry, this `shell model' has been surprisingly successful at reproducing observed Ly$\alpha$ line shapes. In this paper we perform automated line fitting on a set of noisy simulated shell model spectra, in order to determine whether degeneracies exist between the different shell model parameters. While there are some significant degeneracies, we find that most parameters are accurately recovered, especially the HI column density ($N_{\rm HI}$) and outflow velocity ($v_{\rm exp}$). This work represents an important first step in determining how the shell model parameters relate to the actual physical properties of Ly$\alpha$ sources. To aid further exploration of the parameter space, we have made our simulated model spectra available through an interactive online tool.

We explore the impact of Fermi-like acceleration of Lyman-alpha (Ly{\alpha}) photons across shock fronts on the observed Ly{\alpha} spectral line shape. We first confirm the result of Neufeld & McKee (1988) that this mechanism gives rise to extended blue wings which may have been observed in some radio galaxies. Our Monte-Carlo radiative transfer calculations further show that in a minor modification of the shell-model, in which we add an additional static shell of hydrogen, this process can naturally explain the small blue bumps observed in a subset of Ly{\alpha} emitting galaxies, which have been difficult to explain with conventional shell-models. Blue bumps can be produced with an additional column density of static hydrogen as small as $N_{HI}^{static} \ll N_{HI}^{shell}$, and typically occur at roughly the outflow velocity of the shell. In our model the spectra of so-called 'blue-bump objects' might reflect an evolutionary stage in which the outflows regulating the escape of Ly{\alpha} photons are still engulfed within a static interstellar medium.

We give an overview of the Grism Lens Amplified Survey from Space (GLASS), a large Hubble Space Telescope program aimed at obtaining grism spectroscopy of the fields of ten massive clusters of galaxies at redshift z=0.308-0.686, including the Hubble Frontier Fields (HFF). The Wide Field Camera 3 yields near infrared spectra of the cluster cores, covering the wavelength range 0.81-1.69mum through grisms G102 and G141, while the Advanced Camera for Surveys in parallel mode provides G800L spectra of the infall regions of the clusters. The WFC3 spectra are taken at two almost orthogonal position angles in order to minimize the effects of confusion. After summarizing the scientific drivers of GLASS, we describe the sample selection as well as the observing strategy and data processing pipeline. We then utilize MACSJ0717.5+3745, a HFF cluster and the first one observed by GLASS, to illustrate the data quality and the high-level data products. Each spectrum brighter than H_AB=23 is visually inspected by at least two co-authors and a redshift is measured when sufficient information is present in the spectra. Furthermore, we conducted a thorough search for emission lines through all the GLASS WFC3 spectra with the aim of measuring redshifts for sources with continuum fainter than H_AB=23. We provide a catalog of 139 emission-line based spectroscopic redshifts for extragalactic sources, including three new redshifts of multiple image systems (one probable, two tentative). In addition to the data itself we also release software tools that are helpful to navigate the data.

We present new results from our search for z~7 galaxies from deep spectroscopic observations of candidate z-dropouts in the CANDELS fields. Despite the extremely low flux limits achieved by our sensitive observations, only 2 galaxies have robust redshift identifications, one from its Lyalpha emission line at z=6.65, the other from its Lyman-break, i.e. the continuum discontinuity at the Lyalpha wavelength consistent with a redshift 6.42, but with no emission line. In addition, for 23 galaxies we present deep limits in the Lyalpha EW derived from the non detections in ultra-deep observations. Using this new data as well as previous samples, we assemble a total of 68 candidate z~7 galaxies with deep spectroscopic observations, of which 12 have a line detection. With this much enlarged sample we can place solid constraints on the declining fraction of Ly$\alpha$ emission in z~7 Lyman break galaxies compared to z~6, both for bright and faint galaxies. Applying a simple analytical model, we show that the present data favor a patchy reionization process rather than a smooth one.

We predict Lyman-$\alpha$ (Ly$\alpha$) luminosity functions (LFs) of Ly$\alpha$-selected galaxies (Ly$\alpha$ emitters, or LAEs) at $z=3-6$ using the phenomenological model of Dijkstra & Wyithe (2012). This model combines observed UV-LFs of Lyman-break galaxies (LBGs, or drop out galaxies), with constraints on their distribution of Ly$\alpha$ line strengths as a function of UV-luminosity and redshift. Our analysis shows that while Ly$\alpha$ LFs of LAEs are generally not Schechter functions, these provide a good description over the luminosity range of $\log_{10}( L_{\alpha}/{\rm erg}\,{\rm s}^{-1})=41-44$. Motivated by this result, we predict Schechter function parameters at $z=3-6$. Our analysis further shows that (i) the faint end slope of the Ly$\alpha$ LF is steeper than that of the UV-LF of Lyman-break galaxies, (with a median $\alpha_{Ly\alpha} < -2.0$ at $z\gtrsim 4$), and (ii) a turn-over in the Ly$\alpha$ LF of LAEs at Ly$\alpha$ luminosities $10^{40}$ erg s$^{-1}<L_{\alpha}\lesssim 10^{41}$ erg s$^{-1}$ may signal a flattening of UV-LF of Lyman-break galaxies at $-12>M_{\rm UV}>-14$. We discuss the implications of these results - which can be tested directly with upcoming surveys - for the Epoch of Reionization.

We present new Hubble Space Telescope Cosmic Origins Spectrograph far-ultraviolet (far-UV) spectroscopy and Keck Echellete optical spectroscopy of 11 ultraluminous infrared galaxies (ULIRGs), a rare population of local galaxies experiencing massive gas inflows, extreme starbursts, and prominent outflows. We detect H Lyman alpha emission from 8 ULIRGs and the companion to IRAS09583+4714. In contrast to the P Cygni profiles often seen in galaxy spectra, the H Lyman alpha profiles exhibit prominent, blueshifted emission out to Doppler shifts exceeding -1000 km/s in three HII-dominated and two AGN-dominated ULIRGs. To better understand the role of resonance scattering in shaping the H Lyman alpha line profiles, we directly compare them to non-resonant emission lines in optical spectra. We find that the line wings are already present in the intrinsic nebular spectra, and scattering merely enhances the wings relative to the line core. The H Lyman alpha attenuation (as measured in the COS aperture) ranges from that of the far-UV continuum to over 100 times more. A simple radiative transfer model suggests the H Lyman alpha photons escape through cavities which have low column densities of neutral hydrogen and become optically thin to the Lyman continuum in the most advanced mergers. We show that the properties of the highly blueshifted line wings on the H Lyman alpha and optical emission-line profiles are consistent with emission from clumps of gas condensing out of a fast, hot wind. The luminosity of the H Lyman alpha emission increases non-linearly with the ULIRG bolometric luminosity and represents about 0.1 to 1% of the radiative cooling from the hot winds in the HII-dominated ULIRGs.

We present results of radiative transfer calculations to measure the impact of gas bulk rotation on the morphology of the Lyman $\alpha$ emission line in distant galaxies. We model a galaxy as a sphere with an homogeneous mixture of dust and hydrogen at a constant temperature. These spheres undergo solid-body rotation with maximum velocities in the range 0-300 \kms and neutral hydrogen optical depths in the range $\tau_{\rm H}=10^{5}-10^{7}$. We consider two types of source distributions in the sphere: central and homogeneous. Our main result is that rotation introduces a dependence of the line morphology with viewing angle and rotational velocity. Observations with a line of sight parallel to the rotation axis yield line morphologies similar to the static case. For lines of sight perpendicular to the rotation axis both the intensity at the line center and the line width increase with rotational velocity. Along the same line of sight, the line becomes single peaked at rotational velocities close to half the line width in the static case. Notably, we find that rotation does not induce any spatial anisotropy in the integrated line flux, the escape fraction or the average number of scatterings. This is because Lyman {\alpha} scattering through a rotating solid-body proceeds identical as in the static case. The only difference is the doppler shift from the different regions in the sphere that move with respect to the observer. This allows us to derive an analytic approximation for the viewing-angle dependence of the emerging spectrum, as a function of rotational velocity.

We quantify the directional dependence of the escape fraction of Lyman-$\alpha$ (Ly$\alpha$) and non-ionising UV-continuum photons from a multiphase medium, and investigate whether there exist directional enhancements in the Ly$\alpha$ equivalent width (EW). Our multiphase medium consists of spherically symmetric distributions of cold, dusty clumps embedded within a hot dust-free medium. We focus on three models from the analysis presented by Laursen et al. (2013). We find that for a Ly$\alpha$ and UV-continuum point source, it is possible to find an EW boost $b(\theta,\phi) > 5 \bar{b}$ in a few per cent of sight lines, where $\bar{b}$ denotes the boost averaged over all photons. For spatially extended sources this directional dependence vanishes quickly when the size of the UV emitting region exceeds the mean distance between cold dusty clumps. Our analysis suggests that directional EW boosting can occur, and that this is mostly driven by reduced escape fractions of UV photons (which gives rise to UV-continuum `shadows'), and less due to an enhanced Ly$\alpha$ escape fraction (or beaming), in certain directions.

The Epoch of Reionization (EoR) represents a milestone in the evolution of our Universe. Star-forming galaxies that existed during the EoR likely emitted a significant fraction (~5-40%) of their bolometric luminosity as Lyman Alpha (Lya) line emission. However, neutral intergalactic gas that existed during the EoR was opaque to Lya emission that escaped from galaxies during this epoch, which makes it difficult to observe. The neutral intergalactic medium (IGM) may thus reveal itself by suppressing the Lya flux from background galaxies. Interestingly, a `sudden' reduction in the observed Lya flux has now been observed in galaxies at z >6. This review contains a detailed summary of Lya radiative processes: I describe (i) the main Lya emission processes, including collisional-excitation & recombination (and derive the famous factor `0.68'), and (ii) basic radiative transfer concepts, including e.g. partially coherent scattering, frequency diffusion, resonant versus wing scattering, optically thick versus 'extremely' optically thick (static/outflowing/collapsing) media, and multiphase media. Following this review, I derive expressions for the Gunn-Peterson optical depth of the IGM during (inhomogeneous) reionization and post-reionization. I then describe why current observations appear to require a very rapid evolution of volume-averaged neutral fraction of hydrogen in the context of realistic inhomogeneous reionization models, and discuss uncertainties in this interpretation. Finally, I describe how existing & futures surveys and instruments can help reduce these uncertainties, and allow us to fully exploit Lya emitting galaxies as a probe of the EoR.

The large cross-section of the Lyman alpha (Lya) line makes it a sensitive probe of the ionization state of the intergalactic medium (IGM). Here we present the most complete study to date of the IGM Lya opacity, and its application to the redshift evolution of the 'Lya fraction', i.e. the fraction of color-selected galaxies with a detectable Lya emission line. We use a tiered approach, which combines large-scale semi-numeric simulations of reionization with moderate-scale hydrodynamic simulations of the ionized IGM. This allows us to simultaneously account for evolution in both: (i) the opacity from an incomplete (patchy) reionization, parameterized by the filling factor of ionized regions, Q_HII; and (ii) the opacity from self-shielded systems in the ionized IGM, parameterized by the average photo-ionization rate inside HII regions, \Gamma. In contrast to recent empirical models, attenuation from patchy reionization has a unimodal distribution along different sightlines, while attenuation from self-shielded systems is more bimodal. We quantify the average IGM transmission in our (Q_HII, \Gamma) parameter space, which can easily be used to interpret new data sets. Using current observations, we predict that the Lya fraction cannot drop by more than a factor of ~2 with IGM attenuation alone, even for HII filling factors as low as Q_HII>0.1. Larger changes in the Lya fraction could result from a co-evolution with galaxy properties. Marginalizing over \Gamma, we find that current observations constrain Q_HII < 0.6 at z=7 [68% confidence level (C.L.)]. However, all of our parameter space is consistent with observations at 95% C.L., highlighting the need for larger observational samples at z >= 6.

The nature of the seeds of high-redshift supermassive black holes (SMBHs) is a key question in cosmology. Direct collapse black holes (DCBH) that form in pristine, atomic-line cooling halos, illuminated by a Lyman-Werner (LW) UV flux exceeding a critical threshold J_crit, represent an attractive possibility. We investigate when and where these conditions are met during cosmic evolution. For the LW intensity, J_LW, we account for departures from the background value in close proximity to star forming galaxies. For the pristine halo fraction, we account for both (i) supernova driven outflows, and (ii) the inherent pollution from progenitor halos. We estimate the abundance of DCBH formation sites, n_DCBH(z), and find that it increases with cosmic time from n_DCBH(z=20) ~ 1e-12 -1e-7 cMpc^-3 to n_DCBH(z=10) ~ 1e-10 - 1e-5 cMpc^-3. Our analysis shows the possible importance of galactic winds, which can suppress the predicted n_DCBH by several orders of magnitude, and cause DCBH formation to preferentially occur around the UV-brightest (M_UV ~ -22 to -20) star forming galaxies. Our analysis further highlights the dependence of these predictions on (i) the escape fraction of LW photons, (ii) J_crit, and (iii) the galactic outflow prescription.

The rapid decline in the number of strong Lyman Alpha (Lya) emitting galaxies at z > 6 provides evidence for neutral hydrogen in the IGM, but is difficult to explain with plausible models for reionization. We demonstrate that the observed reduction in Lya flux from galaxies at z > 6 can be explained by evolution in the escape fraction of ionizing photons, f_esc. We find that the median observed drop in the fraction of galaxies showing strong Lya emission, as well as the observed evolution of the Lya luminosity function both follow from a small increase in f_esc of Delta f_esc ~ 0.1 from f_esc ~ 0.6 at z ~ 6. This high escape fraction may be at odds with current constraints on the ionising photon escape fraction, which favor smaller values of f_esc < 20%. However, models that invoke a redshift evolution of f_ esc that is consistent with these constraints can suppress the z~7 Lya flux to the observed level, if they also include a small evolution in global neutral fraction of Delta x_HI ~ 0.2. Thus, an evolving escape fraction of ionising photons can be a plausible part of the explanation for evolution in the Lya emission of high redshift galaxies. More generally, our analysis also shows that the drop in the Lya fraction is quantitatively consistent with the observed evolution in the Lya luminosity functions of Lya Emitters.

The Grism Lens-Amplified Survey from Space (GLASS) is a Hubble Space Telescope (HST) Large Program, which will obtain 140 orbits of grism spectroscopy of the core and infall regions of 10 galaxy clusters, selected to be among the very best cosmic telescopes. Extensive HST imaging is available from many sources including the CLASH and Frontier Fields programs. We introduce the survey by analyzing spectra of faint multiply-imaged galaxies and $z\gtrsim6$ galaxy candidates obtained from the first seven orbits out of fourteen targeting the core of the Frontier Fields cluster MACS0717.5+3745. Using the G102 and G141 grisms to cover the wavelength range 0.8-1.7$\mu$m, we confirm 4 strongly lensed systems by detecting emission lines in each of the images. For the 9 $z\gtrsim6$ galaxy candidates clear from contamination, we do not detect any emission lines down to a seven-orbit 1$\sigma$ noise level of $\sim$5$\times$10$^{-18}$erg s$^{-1}$cm$^{-2}$. Taking lensing magnification into account, our flux sensitivity reaches $\sim$0.2-5$\times$10$^{-18}$erg s$^{-1}$cm$^{-2}$. These limits over an uninterrupted wavelength range rule out the possibility that the high-$z$ galaxy candidates are instead strong line emitters at lower redshift. These results show that by means of careful modeling of the background - and with the assistance of lensing magnification - interesting flux limits can be reached for large numbers of objects, avoiding pre-selection and the wavelength restrictions inherent to ground-based multi-slit spectroscopy. These observations confirm the power of slitless HST spectroscopy even in fields as crowded as a cluster core.

We investigate the radiative transfer of Lyman alpha photons through simplified anisotropic gas distributions, which represent physically motivated extensions of the popular 'shell-models'. Our study is motivated by the notion that (i) shell models do not always reproduce observed Lyman alpha spectral line profiles, (ii) (typical) shell models do not allow for the escape of ionizing photons, and (iii) the observation & expectation that winds are more complex, anisotropic phenomena. We examine the influence of inclination on the Lyman alpha spectra, relative fluxes and escape fractions. We find the flux to be enhanced/suppressed by factors up to a few depending on the parameter range of the models, corresponding to a boost in equivalent width of the same amplitude if we neglect dust. In general, lower mean optical depths tend to reduce the impact of anisotropies as is expected. We find a correlation between an observed peak in the -- occasionally triple-peaked -- spectrum at the systemic velocity and the existence of a low optical depth cavity along the line of sight. This can be of importance in the search for ionizing photons leaking from high- redshift galaxies since these photons will also be able to escape through the cavity.

We derive empirical constraints on the volume averaged `effective' escape fraction of Lyman Alpha (Lya) photons from star forming galaxies as a function of redshift, by comparing star formation functions inferred directly from observations, to observed Lya luminosity functions. Our analysis shows that the effective escape fraction increases from f_esc^eff ~ 1-5 % at z=0, to f_esc^eff ~ 10 % at z=3-4, and to f_esc^eff=30-50 % at z=6. Our constraint at z=6 lies above predictions by models that do not include winds, and therefore hints at the importance of winds in the Lya transfer process (even) at this redshift. We can reproduce Lya luminosity functions with an f_esc^eff that does not depend on the galaxies star formation rates (SFR) over up to ~2 orders of magnitude in Lya luminosity. It is possible to reproduce the luminosity functions with an f_esc^eff that decreases with SFR - which appears favored by observations of drop-out galaxies - in models which include a large scatter (~ 1.0 dex) in f_esc^eff, and/or in which star forming galaxies only have a non-zero f_esc^eff for a fraction of their life-time or a fraction of sightlines. We provide a fitting formula that summarizes our findings.

We have used the SINFONI near-infrared integral field unit on the VLT to resolve the optical emission line structure of one of the brightest (L~1e44 erg/s) and nearest (z=2.38) of all Lya blobs (LABs). The target, known in the literature as object 'B1' (Francis et al. 1996), lies at a redshift where the main optical emission lines are accessible in the observed near-infrared. We detect luminous [OIII]4959,5007A and Ha emission with a spatial extent of at least 32x40 kpc (4"x5"). The dominant optical emission line component shows relatively broad lines (600-800 km/s, FWHM) and line ratios consistent with AGN-photoionization. The new evidence for AGN photoionization, combined with previously detected CIV and luminous, warm infrared emission, suggest that B1 is the site of a hidden quasar. This is confirmed by the fact that [OII] is relatively weak compared to [OIII] (extinction-corrected [OIII]/[OII] of about 3.8), which is indicative of a high, Seyfert-like ionization parameter. From the [OIII] luminosity we infer a bolometric AGN luminosity of ~3e46 erg/s, and further conclude that the obscured AGN may be Compton-thick given existing X-ray limits. The large line widths observed are consistent with clouds moving within the narrow line region of a luminous QSO. The AGN scenario is capable of producing sufficient ionizing photons to power the Lya, even in the presence of dust. By performing a census of similar objects in the literature, we find that virtually all luminous LABs harbor obscured quasars. Based on simple duty-cycle arguments, we conclude that AGN are the main drivers of the Lya in LABs rather than the gravitational heating and subsequent cooling suggested by cold stream models. We also conclude that the empirical relation between LABs and overdense environments at high redshift must be due to a more fundamental correlation between AGN (or massive galaxies) and environment.

Observations of distant galaxies play a key role in improving our understanding of the Epoch of Reionization (EoR). The observed Ly-alpha emission line strength - quantified by its restframe equivalent width (EW) - provides a valuable diagnostic of stellar populations and dust in galaxies during and after the EoR. In this paper we quantify the effects of star formation stochasticity on the predicted Ly-alpha EW in dwarf galaxies, using the publicly available code SLUG ('Stochastically Light Up Galaxies'). We compute the number of hydrogen ionizing photons, as well as flux in the Far UV for a set of models with star formation rates (SFR) in the range 10-3-1 Msol/yr. From these fluxes we compute the luminosity, L-alpha, and the EW of the Ly-alpha line. We find that stochasticity alone induces a broad distribution in L-alpha and EW at a fixed SFR, and that the widths of these distributions decrease with increasing SFR. We parameterize the EW probability density function (PDF) as an SFR-dependent double power law. We find that it is possible to have EW as low as ~EW0/4 and as high as ~3 times the EW0, where EW0 denotes the expected EW in the absence of stochasticity. We argue that stochasticity may therefore be important when linking drop-out and narrow-band selected galaxies, when identifying population III galaxies, and that it may help to explain the large EW (EW > 100 - 200 A) observed for a fraction of Ly-alpha- selected galaxies. Finally, we show that stochasticity can also affect the inferred escape fraction of ionizing photons from galaxies. In particular, we argue that stochasticity may simultaneously explain the observed anomalous ratios of the Lyman continuum flux density to the (non-ionizing) UV continuum density in so-called Lyman-Bump galaxies at z = 3.1, as well as the absence of such objects among a sample of z = 1.3 drop-out galaxies.

We perform a systematic study of how the inhomogeneities in the Inter-Galactic Medium (IGM) affect the observability of Lyman Alpha Emitters (LAEs) around the Epoch of Reionization. We focus on the IGM close to the galaxies as the detailed ionization distribution and velocity fields of this region could significantly influence the scattering of Ly-alpha photons off neutral H atoms as they traverse the IGM after escaping from the galaxy. We simulate the surface brightness (SB) maps and spectra of more than 100 LAEs at z=7.7 as seen by an observer at z=0. To achieve this, we extract the source properties of galaxies and their surrounding IGM from cosmological simulations of box sizes 5-30 Mpc/h and follow the coupled radiative transfer of ionizing and Ly-alpha radiation through the IGM using CRASH-alpha. We find that the simulated SB profiles are extended and their detailed structure is affected by inhomogeneities in the IGM, especially at high neutral fractions. The detectability of LAEs and the fraction of the flux observed depend heavily on the shape of the SB profile and the SB threshold (SB_th) of the observational campaign. Only ultradeep observations (e.g. SB_th ~ 10^-23 ergs/s/cm^2/arcsec^2) would be able to obtain the true underlying mass-luminosity relation and luminosity functions of LAEs. The details of our results depend on whether Ly-alpha photons are significantly shifted in the galaxy to longer wavelengths, the mean ionization fraction in the IGM and the clustering of ionizing sources. These effects can lead to an easier escape of Ly-alpha photons with less scattering in the IGM and a concentrated SB profile similar to the one of a point source. Finally, we show that the SB profiles are steeper at high ionization fraction for the same LAE sample which can potentially be observed from the stacked profile of a large number of LAEs.

We present constrained radiative transfer calculations of Lyman Alpha (Lya) photons propagating through clumpy, dusty, large scale outflows, and explore whether we can quantitatively explain the Lya halos that have been observed around Lyman Break Galaxies. We construct phenomenological models of large-scale outflows which consist of cold clumps that are in pressure equilibrium with a constant-velocity hot wind. First we consider models in which the cold clumps are distributed symmetrically around the galaxy, and in which the clumps undergo a continuous acceleration in its 'circumgalactic' medium (CGM). We constrain the properties of the cold clumps (radius, velocity, HI column density, & number density) by matching the observed Lya absorption strength of the CGM in the spectra of background galaxies. We then insert a Lya source in the center of this clumpy outflow, which consists of 1e5-1e6 clumps, and compute observable properties of the scattered Lya photons. In these models, the scattered radiation forms halos that are significantly more concentrated than is observed. In order to simultaneously reproduce the observed Lya absorption line strengths and the Lya halos, we require - preferably bipolar - outflows in which the clumps decelerate after their initial acceleration. This deceleration is predicted naturally in 'momentum-driven' wind models of clumpy outflows. In models that simultaneously fit the absorption and emission line data, the predicted linear polarization is ~30-40% at a surface brightness contour of S=1e-18 erg/s/cm^2/arcsec^2. Our work illustrates clearly that Lya emission line halos around starforming galaxies provide valuable constraints on the cold gas distribution & kinematics in their circumgalactic medium, and that these constraints complement those obtained from absorption line studies alone.

We study the visibility of the Lyman Alpha (Lya) emission line during the Epoch of Reionization (EoR). Combining galactic outflow models with large-scale semi-numeric simulations of reionization, we quantify the probability distribution function (PDF) of the fraction of Lya photons transmitted through the intergalactic medium (IGM), T_IGM. Our study focusses on galaxies populating dark matter halos with masses of M_halo=1e10 M_sun at z=8.6, which is inspired by the recent reported discovery of a galaxy at z=8.6 with strong Lya line emission. For reasonable model assumptions, we find that winds cause T_IGM>10% [50%], for the majority of galaxies, even when the Universe is ~80% [60%] neutral by volume. Thus, the observed strong Lya emission from the reported z=8.6 galaxy is consistent with a highly neutral IGM. We also investigate the implications of the recent tentative evidence for a observed decrease in the `LAE fraction' among drop-out galaxies between z=6 and z=7. If confirmed, we show that a rapid evolution in x_HI will be required to explain this observation via the effects of reionization.

Observations of local star forming galaxies have revealed a correlation between the rate at which galaxies form stars and their X-Ray luminosity. We combine this correlation with the most recent observational constraints on the integrated star formation rate density, and find that star forming galaxies account for 5-20% of the total soft and hard X-ray backgrounds, where the precise number depends on the energy band and the assumed average X-ray spectral energy distribution of the galaxies below ~20 keV. If we combine the L_X-SFR relation with recently derived star formation rate function, then we find that star forming galaxies whose X-ray flux falls well (more than a factor of 10) below the detection thresholds of the Chandra Deep Fields, can fully account for the unresolved soft X-ray background, which corresponds to ~6% of its total. Motivated by this result, we put limits on the allowed redshift evolution of the parameter c_X \equiv L_X/SFR, and/or its evolution towards lower and higher star formation rates. If we parametrize the redshift evolution of c_X ~ (1+z)^b, then we find that b \leq 1.3 (95% CL). On the other hand, the observed X-ray luminosity functions (XLFs) of star forming galaxies indicate that c_X may be increasing towards higher redshifts and/or higher star formation rates at levels that are consistent with the X-ray background, but possibly at odds with the locally observed L_X-SFR relation.

Lyman Alpha Emitters (LAEs) are galaxies that have been selected on the basis of a strong Ly{\alpha} emission line in their spectra. Observational campaigns have dramatically increased the sample of known LAEs, which now extends out to z=7. These discoveries have motivated numerous theoretical studies on the subject, which usually define LAEs in their models based on sharp Ly{\alpha} luminosity and equivalent width (EW) cuts. While broadly representative, this procedure does not mimic the selection from observational programs in detail, which instead use cuts in various colour-spaces. We investigate what implications this disjoint may have for studies that aim to model LAEs. We construct an empirical model for the number density of star forming galaxies as a function of their UV and Ly{\alpha} luminosity, utilising measured constraints on the luminosity functions (LFs) of drop-out galaxies, and their luminosity dependent probability distribution function of Ly{\alpha} EW. In particular, we investigate whether the LAE LFs can be reproduced by defining LAEs using a (z-dependent) Ly{\alpha} luminosity and EW threshold. While we are able to reproduce the observed distribution of Ly{\alpha} EW among LAEs out to restframe EW 200 A, we find that our formalism over-predicts both the UV and Ly{\alpha} LFs of LAEs by a factor of 2-3, and is inconsistent with observations at the ~95% level. This tension is partially resolved if we assume the Ly{\alpha} EW-distribution of drop-out galaxies to be truncated at restframe EW>150 A. However the overprediction indicates that modeling LAEs with simple REW and luminosity cuts does not accurately mimic observed selection criteria, and can lead to uncertainties in the predicted number density of LAEs. On the other hand, the predicted z-evolution is not affected. We apply our formalism to drop-out galaxies at z>6, and predict the LFs of LAEs at z=7-9.

We present the final results from our ultra-deep spectroscopic campaign with FORS2 at the ESO/VLT for the confirmation of z~7 "z--band dropout" candidates selected from our VLT/Hawk-I imaging survey over three independent fields. In particular we report on two newly discovered galaxies at redshift ~6.7 in the NTT deep field: both galaxies show a Ly-alpha emission line with rest-frame EWs of the order 15-20 A and luminosities of 2-4 X 10^{42} erg/s. We also present the results of ultra-deep observations of a sample of i-dropout galaxies, from which we set a solid upper limit on the fraction of interlopers. Out of the 20 z-dropouts observed we confirm 5 galaxies at 6.7 < z < 7.1. This is systematically below the expectations drawn on the basis of lower redshift observations: in particular there is a significant lack of objects with intermediate Ly-alpha EWs (between 20 and 55 A). We conclude that the trend for the fraction of Ly-alpha emission in LBGs that is constantly increasing from z~3 to z~6 is most probably reversed from z~6 to z~7. Explaining the observed rapid change in the LAE fraction among the drop-out population with reionization requires a fast evolution of the neutral fraction of hydrogen in the Universe. Assuming that the Universe is completely ionized at z=6 and adopting the semi-analytical models of Dijkstra et al. (2011), we find that our data require a change of the neutral hydrogen fraction of the order Delta chi_{HI} ~ 0.6 in a time Delta z ~ 1, provided that the escape fraction does not increase dramatically over the same redshift interval.

Narrowband observations have detected no Ly-alpha emission within a 70 pMpc^3 volume centered on the z = 2.168 quasar PKS 0424-131. This is in contrast to surveys of Ly-alpha emitters in the field at similar redshifts and flux limits, which indicate that tens of sources should be visible within the same volume. The observed difference indicates that the quasar environment has a significant influence on the observed density of Ly-alpha emitters. To quantify this effect we have constructed a semi-analytic model to simulate the effect of a luminous quasar on nearby Ly-alpha emitters. We find the null detection around PKS 0424-131 implies that the minimum isothermal temperature of Ly-alpha emitter host halos is greater than 3.4 x 10^6 K (68% level), corresponding to a virial mass of ~1.2 x 10^12 M_solar. This indicates that the intense UV emission of the quasar may be suppressing the star formation in nearby galaxies. Our study illustrates that low redshift quasar environments may serve as a surrogate for studying the radiative suppression of galaxy formation during the epoch of reionization.

We show that the dependence of Ly-alpha absorption on environment leads to significant non-gravitational features in the redshift space power-spectrum of Ly-alpha selected galaxies. We derive a physically motivated fitting formula that can be included in clustering analyses, and use this to discuss the predicted features in the Ly-alpha galaxy power-spectrum based on detailed models in which Ly-alpha absorption is influenced by gas infall and/or by strong galactic outflows. We show that power-spectrum measurements could be used to study the astrophysics of the galaxy-IGM connection, and to measure the properties of outflows from star-forming galaxies. Applying the modified redshift space power-spectrum to a Ly-alpha survey with parameters corresponding to the planned Hobby-Eberly Telescope Dark Energy Experiment (HETDEX), we find that the dependence of observed Ly-alpha flux on velocity gradient and ionising background may compromise the ability of Ly-alpha selected galaxy redshift surveys to constrain cosmology using information from the full power-spectrum. This is because the effects of fluctuating ionizing background and velocity gradients effect the shape of the observed power-spectrum in ways that are similar to the shape of the primordial power-spectrum and redshift space distortions respectively. We use the Alcock-Paczynski test to show that without prior knowledge of the details of Ly-alpha absorption in the IGM, the precision of line-of-sight and transverse distance measurements for HETDEX will be ~1.3-1.7%, decreased by a factor of ~1.5-2 relative to the best case precision of ~0.8% available in a traditional galaxy redshift survey. We specify the precision with which modelling of Ly-alpha radiative transfer must be understood in order for HETDEX to achieve distance measurements that are better than 1%.

The next generation of telescopes aim to directly observe the first generation of galaxies that initiated the reionization process in our Universe. The Lyman Alpha (Lya) emission line is robustly predicted to be the most prominent intrinsic spectral feature of these galaxies, making it an ideal target to search for and study high-redshift galaxies. I briefly discuss why Lya emitting galaxies (LAEs) provide a good probe of the intergalactic medium (IGM) during the Epoch of Reionization (EoR). I argue that if we wish to fully exploit LAEs as a probe of the EoR, it is important to understand what drives their observed redshift evolution after reionization is completed. One important uncertainty in interpreting existing LAE observations relates to the impact of the ionized IGM on Lya photons emitted by galaxies, which is strongly connected to the effects of scattering through HI in galactic outflows. Scattering through galactic outflows can also modify the Lya spectral line shape such that > 5% of the emitted Lya radiation is transmitted directly to the observer, even through a fully neutral IGM. Finally, I discuss what is required --observationally and theoretically-- to resolve the uncertainties that affect existing interpretations of data on LAEs.

It is well established that between 380000 and 1 billion years after the Big Bang the Inter Galactic Medium (IGM) underwent a "phase transformation" from cold and fully neutral to warm (~10^4 K) and ionized. Whether this phase transformation was fully driven and completed by photoionization by young hot stars is a question of topical interest in cosmology. AIMS. We propose here that besides the ultraviolet radiation from massive stars, feedback from accreting black holes in high-mass X-ray binaries (BH-HMXBs) was an additional, important source of heating and reionization of the IGM in regions of low gas density at large distances from star-forming galaxies. METHODS. We use current theoretical models on the formation and evolution of primitive massive stars of low metallicity, and the observations of compact stellar remnants in the near and distant universe, to infer that a significant fraction of the first generations of massive stars end up as BH-HMXBs. The total number of energetic ionizing photons from an accreting stellar black hole in an HMXB is comparable to the total number of ionizing photons of its progenitor star. However, the X-ray photons emitted by the accreting black hole are capable of producing several secondary ionizations and the ionizing power of the resulting black hole could be greater than that of its progenitor. Feedback by the large populations of BH-HMXBs heats the IGM to temperatures of ~10^4 K and maintains it ionized on large distance scales. BH-HMXBs determine the early thermal history of the universe and mantain it as ionized over large volumes of space in regions of low density. This has a direct impact on the properties of the faintest galaxies at high redshifts, the smallest dwarf galaxies in the local universe, and on the existing and future surveys at radio wavelengths of atomic hydrogen in the early universe.

Recent studies have shown that galaxies accrete most of their baryons via the cold mode, from streams with temperatures T~10^4-10^5 K. At these temperatures, the streams should radiate primarily in the Lya line and have therefore been proposed as a model to power the extended, high-redshift objects known as Lya blobs and other high-redshift Lya sources. We introduce a new Lya radiative transfer code, aRT, and apply it to cosmological hydrodynamical simulations. We address physical and numerical issues that are critical to making accurate predictions for the cooling luminosity, but that have been mostly neglected or treated simplistically so far. We highlight the importance of self-shielding and of properly treating sub-resolution models in simulations. Most existing simulations do not self-consistently incorporate these effects, which can lead to order-of-magnitude errors in the predicted cooling luminosity. Using a combination of post-processing ionizing radiative transfer and re-simulation techniques, we develop an approximation to the consistent evolution of the self-shielded gas. We quantify the dependence of the Lya cooling luminosity on halo mass at z=3 for the simplified problem of pure gas accretion. While cooling in massive halos (without additional energy input from star formation and AGN) is in principle sufficient to produce L_alpha~10^43-10^44 erg s^-1 blobs, this requires including energy released in gas of density sufficient to form stars, but which is kept 100% gaseous in our optimistic estimates. Excluding emission from such dense gas yields lower luminosities by up to one to two orders of magnitude at high masses, making it difficult to explain the observed Lya blobs with pure cooling. Resonant scattering produces diffuse Lya halos, even for centrally concentrated emission, and broad double peaked line profiles. [Abridged]

The next generation of telescopes aim to directly observe the first generation of galaxies that initiated the reionization process in our Universe. The Lyman Alpha (Lya) emission line is robustly predicted to be the most prominent intrinsic spectral feature of these galaxies, making it an ideal target to search for and study high redshift galaxies. Unfortunately the large Gunn-Peterson optical depth of the surrounding neutral intergalactic medium (IGM) is thought to render this line extremely difficult to detect prior to reionization. In this paper we demonstrate that the radiative transfer effects in the interstellar medium (ISM), which cause Lya flux to emerge from galaxies at frequencies where the Gunn-Peterson optical depth is reduced, can substantially enhance the prospects for detection of the Lya line at high redshift. In particular, scattering off outflows of interstellar HI gas can modify the Lya spectral line shape such that >5% of the emitted Lya radiation is transmitted directly to the observer, even through a fully neutral IGM. It may therefore be possible to directly observe `strong' Lya emission lines (EW > 50 Angstrom rest frame) from the highest redshift galaxies that reside in the smallest HII `bubbles' early in the reionization era with JWST. In addition, we show that outflows can boost the fraction of Lya radiation that is transmitted through the IGM during the latter stages of reionization, and even post-reionization. Coupled with the fact that the first generation of galaxies are thought to have very large intrinsic equivalent Lya equivalent widths, our results suggest that the search for galaxies in their redshifted Lya emission line can be competitive with the drop-out technique out to the highest redshifts that can be probed in the JWST era.

The equivalent width (EW) of the Lyman Alpha (Lya) line is directly related to the ratio of star formation rates determined from Lya flux and UV flux density [SFR(Lya)/SFR(UV)]. We use published data --in the literature EW and SFR(Lya)/SFR(UV) are treated as independent quantities-- to show that the predicted relation holds for the vast majority of observed Lya emitting galaxies (LAEs). We show that the relation between EW and SFR(Lya)/SFR(UV) applies irrespective of a galaxy's `true' underlying star formation rate, and that its only source of scatter is the variation in the spectral slope of the UV continuum between individual galaxies. The derived relation, when combined with the observed EW distribution, implies that the ratio SFR(UV)/SFR(Lya) is described well by a log-normal distribution with a standard deviation of ~0.3-0.35. This result is useful when modelling the statistical properties of LAEs. We further discuss why the relation between EW and SFR(Lya)/SFR(UV) may help identifying galaxies with unusual stellar populations.

Recent hydrodynamic simulations of galaxy formation reveal streams of cold (T ~ 1e4 K) gas flowing into the centers of dark matter halos as massive as 1e12-1e13.5 M_sun at redshifts z~1-3. In this paper we show that if > 20% of the gravitational binding energy of the gas is radiated away, then the simulated cold flows are spatially extended Lyman Alpha (Lya) sources with luminosities, Lya line widths, and number densities that are comparable to those of observed Lya blobs. Furthermore, the filamentary structure of the cold flows can explain the wide range of observed Lya blob morphologies. Since the most massive halos form in dense environments, the association of Lya blobs with overdense regions arise naturally. We argue that Lya blobs - even those which are clearly associated with starburst galaxies or quasars - provide direct observational support for the cold accretion mode of galaxies. We discuss various testable predictions of this association.

Active Galactic Nuclei (AGN) produce a dominant fraction (~80%) of the Soft X-ray background (SXB) at photon energies 0.5<E<2 keV. If dust pervaded throughout the intergalactic medium, its scattering opacity would have produced diffuse X-ray halos around AGN. Taking account of known galaxies and galaxy clusters, only a fraction F_halo <10% of the SXB can be in the form of diffuse X-ray halos around AGN. We therefore limit the intergalactic opacity to optical/infrared photons from large dust grains (with radii in the range a=0.2-2.0 mum) to a level tau_GD<0.15(F_halo/10%) to a redshift z~1. Our results are only weakly dependent on the grain size distribution or the redshift evolution of the intergalactic dust. Stacking X-ray images of AGN can be used to improve our constraints and diminish the importance of dust as a source of systematic uncertainty for future supernova surveys which aim to improve the precision on measuring the redshift evolution of the dark energy equation-of-state.

Scattering of Lyman Alpha (hereafter Lya) photons by neutral hydrogen gas in a single outflowing 'supershell' around star forming regions often explains the shape and offset of the observed Lya emission line from galaxies. We compute the radiation pressure that is exerted by this scattered Lya radiation on the outflowing material. We show that for reasonable physical parameters, Lya radiation pressure alone can accelerate supershells to velocities in the range v_sh=200-400 km/s. These supershells possibly escape from the gravitational potential well of their host galaxies and contribute to the enrichment of the intergalactic medium. We compute the physical properties of expanding supershells that are likely to be present in a sample of known high-redshift (z=2.7-5.0) galaxies, under the assumption that they are driven predominantly by Lya radiation pressure. We predict ranges of radii r_sh=0.1-10 kpc, ages t_sh=1-100 Myr, and energies E_sh=1e53-1e55 ergs, which are in reasonable agreement with the properties of local galactic supershells. Furthermore, we find that the radius, r_sh, of a Lya-driven supershell of constant mass depends uniquely on the intrinsic Lya luminosity of the galaxy, L_alpha, the HI column density of the supershell, N_HI, and the shell speed, v_sh, through the scaling relation r_sh ~ L_alpha/(N_HI v_sh^2). We derive mass outflow rates in supershells that reach ~10-100% of the star formation rates of their host galaxies.

A soft ultraviolet radiation field, 10.2 eV < E <13.6 eV, that permeates neutral intergalactic gas during the Epoch of Reionization (EoR) excites the 2p (directly) and 2s (indirectly) states of atomic hydrogen. Because the 2s state is metastable, the lifetime of atoms in this level is relatively long, which may cause the 2s state to be overpopulated relative to the 2p state. It has recently been proposed that for this reason, neutral intergalactic atomic hydrogen gas may be detected in absorption in its 3-cm fine-structure line (2s_1/2 -> 2p_3/2) against the Cosmic Microwave Background out to very high redshifts. In particular, the optical depth in the fine-structure line through neutral intergalactic gas surrounding bright quasars during the EoR may reach tau~1e-5. The resulting surface brightness temperature of tens of micro K (in absorption) may be detectable with existing radio telescopes. Motivated by this exciting proposal, we perform a detailed analysis of the transfer of Lyman beta,gamma,delta,... radiation, and re-analyze the detectability of the fine-structure line in neutral intergalactic gas surrounding high-redshift quasars. We find that proper radiative transfer modeling causes the fine-structure absorption signature to be reduced tremendously to tau< 1e-10. We therefore conclude that neutral intergalactic gas during the EoR cannot reveal its presence in the 3-cm fine-structure line to existing radio telescopes.

The earliest generation of stars and black holes must have established an early 'Lyman-Werner' background (LWB) at high redshift, prior to the epoch of reionization. Because of the long mean free path of photons with energies E<13.6 eV, the LWB was nearly uniform. However, some variation in the LWB is expected due to the discrete nature of the sources, and their highly clustered spatial distribution. In this paper, we compute the probability distribution function (PDF) of the LW flux that irradiates dark matter (DM) halos collapsing at high-redshift (z~10). Our model accounts for (i) the clustering of DM halos, (ii) Poisson fluctuations in the number of corresponding star forming galaxies, and (iii) scatter in the LW luminosity produced by halos of a given mass (calibrated using local observations). We find that > 99% of the DM halos are illuminated by a LW flux within a factor of 2 of the global mean value. However, a small fraction, ~1e-8 to 1e-6, of DM halos with virial temperatures above 1e4 K have a close luminous neighbor within < 10 kpc, and are exposed to a LW flux exceeding the global mean by a factor of > 20, or to J_(21,LW)> 1e3 (in units of 1e-21 erg/s/Hz/sr/cm^2). This large LW flux can photo--dissociate H_2 molecules in the gas collapsing due to atomic cooling in these halos, and prevent its further cooling and fragmentation. Such close halo pairs therefore provide possible sites in which primordial gas clouds collapse directly into massive black holes (M_BH~ 1e4 - 1e6 M_sun), and subsequently grow into supermassive (M_BH > 1e9 M_sun) black holes by z~6.

The collapse of baryons into the center of a host dark matter halo is accompanied by radiation that may be detectable as compact (< 10 kpc) UV-continuum and Lyman Alpha (hereafter Lya) emission with Lya luminosities as high as ~1e42-1e43 erg/s in halos of mass M=1e11-10e12 solar masses. We show that the observed equivalent width (EW) of the Lya line emitted by these cooling clouds is EW 400 Angstrom (restframe). These luminosities and EWs are comparable to those detected in narrowband surveys for redshifted Lya emission. The rest-frame ultraviolet of Lya emitting cooling clouds radiation may be dominated by two-photon transitions from 2s->1s. The resulting spectrum can distinguish cooling clouds from a broad class of young star forming galaxies.

We present accurate Monte-Carlo calculations of Lyman Alpha (hereafter Lya) radiation pressure in a range of models which represent galaxies during various epochs of our Universe. We show that the radiation force that Lya photons exert on hydrogen gas in the neutral intergalactic medium (IGM), that surrounds minihalos that host the first stars, may exceed gravity by orders of magnitude and drive supersonic winds. Lya radiation pressure may also dominate over gravity in the neutral IGM that surrounds the HII regions produced by the first galaxies. However, the radiation force is likely too weak to result in supersonic outflows in this case. Furthermore, we show that Lya radiation pressure may drive outflows in the interstellar medium of star forming galaxies that reach hundreds of km/s. This mechanism could also operate at lower redshifts z<6, and may have already been indirectly detected in the spectral line shape of observed Lya emission lines.

The high-redshift Universe contains luminous Lyman Alpha (hereafter Lya) emitting sources such as galaxies and quasars. The emitted Lya radiation is often scattered by surrounding neutral hydrogen atoms. We show that the scattered Lya radiation obtains a high level of polarization for a wide range of likely environments of high-redshift galaxies. For example, the back-scattered Lya flux observed from galaxies surrounded by a superwind-driven outflow may reach a fractional polarization as high as ~40%. Equal levels of polarization may be observed from neutral collapsing protogalaxies. Resonant scattering in the diffuse intergalactic medium typically results in a lower polarization amplitude (<7%), which depends on the flux of the ionizing background. Spectral polarimetry can differentiate between Lya scattering off infalling gas and outflowing gas; for an outflow the polarization should increase towards longer wavelengths while for infall the opposite is true. Our numerical results suggest that Lya polarimetry is feasible with existing instruments, and may provide a new diagnostic of the distribution and kinematics of neutral hydrogen around high-redshift galaxies. Moreover, polarimetry may help suppress infrared lines originating in the Earth's atmosphere, and thus improve the sensitivity of ground-based observations to high-redshift Lya emitting galaxies outside the currently available redshift windows.

The ionizing ultraviolet background (UVB) during reionization can suppress the gas content of low-mass galaxies, even those capable of efficient atomic cooling, and thus lead to an extended reionization epoch. In this work, we explore the importance of negative UV radiative feedback on Tvir > 10^4 K halos during the middle and late stages of reionization. We do not try to self-consistently model reionization; instead, we explore a large parameter space in an attempt to draw general, robust conclusions. We do this using a tiered approach. Using 1-D hydrodynamical simulations, we model the collapse of gas onto halos of various masses under UVBs of various intensities. We then generate realistic, parametrized maps of the inhomogeneous UVB, using large-scale semi-numeric simulations. By combining these results, we find that under all reasonably conservative scenarios, UV feedback on atomically-cooled halos is not strong enough to notably delay the bulk of reionization. Such a delay is only likely if ionizing efficiencies of z > 10 sources are much higher (~ two orders of magnitude) than z ~ 6 data seem to imply. We also find that feedback is very strongly dependent on halo mass. Our results suggest that the natural time-scale for the bulk of reionization is the growth of the global collapsed fraction contained in Tvir > 10^4 K halos. Finally, our results underscore the importance of taking into account extended dynamical ranges when modeling reionization.

We perform Monte-Carlo calculations of the radiative transfer of Lyman Alpha (hereafter Lya) photons emitted by a source embedded in a neutral collapsing gas cloud. This represents a young galaxy or quasar during the early stages of the epoch of reionisation (EoR). After computing the Lya spectrum as function of radius and time, we find that the Lya color temperature may be negative in large volumes surrounding the central source. Motivated by this result, we explore the prospects for a population inversion in the hyperfine levels of atomic hydrogen via the Wouthuysen-Field (WF) effect. The reason for this exploration is clear: if 21-cm masers exist during the EoR, they could greatly boost the expected 21-cm flux from this epoch. We find that population inversion is unlikely to occur for several reasons: (1) the required Lya luminosities of the central source exceed ~10e45 erg/s. The radiation pressure exerted by such a large Lya flux likely halts the collapse of the cloud; (2) When quantum corrections to the WF-coupling strength are applied, the required Lya luminosities are (even) larger by orders of magnitude; (3) A relatively low flux of Lya photons that is produced via other channels (x-ray heating, collisional excitation of hydrogen) prevents the Lya color temperature from becoming negative.

A significant fraction of Lyman Alpha (Lya) emitting galaxies (LAEs) at z> 5.7 have rest-frame equivalent widths (EW) greater than ~100 Angstrom. However only a small fraction of the Lya flux produced by a galaxy is transmitted through the IGM, which implies intrinsic Lya EWs that are in excess of the maximum allowed for a population-II stellar population having a Salpeter mass function. In this paper we study characteristics of the sources powering Lya emission in high redshift galaxies. We propose a simple model for Lya emitters in which galaxies undergo a burst of very massive star formation that results in a large intrinsic EW, followed by a phase of population-II star formation with a lower EW. We confront this model with a range of high redshift observations and find that the model is able to simultaneously describe the following eight properties of the high redshift galaxy population with plausible values for parameters like the efficiency and duration of star formation: i-iv) the UV and Lya luminosity functions of LAEs at z=5.7 and 6.5, v-vi) the mean and variance of the EW distribution of Lya selected galaxies at z=5.7, vii) the EW distribution of i-drop galaxies at z~6, and viii) the observed correlation of stellar age with EW. Our modeling suggests that the observed anomalously large intrinsic equivalent widths require a burst of very massive star formation lasting no more than a few to ten percent of the galaxies star forming lifetime. This very massive star formation may indicate the presence of population-III star formation in a few per cent of i-drop galaxies, and in about half of the Lya selected galaxies.

Recent observations imply that the observed number counts of Lya Emitters (LAEs) evolved significantly between z=5.7 and z=6.5. It has been suggested that this was due to a rapid evolution in the ionisation state, and hence transmission of the IGM which caused Lya flux from z=6.5 galaxies to be more strongly suppressed. In this paper we consider the joint evolution of the Lya and UV luminosity functions (LFs) and show that the IGM transmission evolved between z=6.5 and z=5.7 by a factor 1.1 <R < 1.8 (95% CL). This result is insensitive to the underlying model of the Lya LF (as well as cosmic variance). Using a model for IGM transmission, we find that the evolution of the mean IGM density through cosmic expansion alone may result in a value for the ratio of transmissions as high as R=1.3. Thus, the existing LFs do not provide evidence for overlap. Furthermore, the constraint R<1.8 suggests that the Universe at z=6.5 was more than half ionised by volume, i.e. x_i,V>0.5.

We calculate the impact of the intergalactic medium (IGM) on the observed Lyman alpha lines (hereafter Lya) emitted by galaxies in an ionised IGM at z>4. Our model accounts for gas clumping in the IGM and for the fact that high-redshift galaxies reside in overdense regions, which causes the velocity field of the IGM to depart from the Hubble flow. The observed shape of the Lya line varies widely, with dependence on the intrinsic width and systemic velocity of the line, a galaxies star formation rate and the local extra-galactic UV-background. For large star formation rates and levels of the UV-background, absorption in the IGM does not result in a Lya line that is asymmetric as is common among known high-redshift Lya emitters. For models in which the lines do show the observed strong asymmetry, the IGM typically transmits only 10-30% of the Lya flux. The increase in the ionising background that accompanied the completion of reionisation barely increased the IGM transmission, which suggests that LAEs of comparable luminosity should not appear to be significantly dimmer prior to overlap. In this light, we briefly discuss the potential of Lya emitters as a probe into the epoch of reionisation.

The morphological and spectroscopic properties of the giant (>60 kpc) Ly-alpha nebulae associated with three radio galaxies at z~2.5 (MRC 1558-003, MRC 2025-218 and MRC 0140-257) have been investigated using integral field spectroscopic data obtained with VIMOS on VLT. The morphologies are varied. The nebula of one source has a centrally peaked, rounded appearance. In the other two objects, it consists of two spatial components. The three nebulae are aligned with the radio axis within <30 deg. The total Ly-alpha luminosities are in the range (0.3-3.4) x 1e44 erg s-1. The Ly-alpha spectral profile shows strong variation through the nebulae, with FWHM values in the range ~400-1500 km s-1 and velocity shifts V~120-600 km s-1. We present an infall model which can explain successfully most Ly-alpha morphological and spectroscopic properties of the nebula associated with MRC 1558-003. This adds further support to our previous conclusion that the _quiescent_ giant nebulae associated with this and other high redshift powerful radio galaxies are in infall. A problem for this model is the difficulty to reproduce the large Ly-alpha FWHM values. We have discovered a giant (~85 kpc) Ly-alpha nebula associated with the radio galaxy MRC 0140-257 at z=2.64. It shows strikingly relaxed kinematics (FWHM< 300 km s-1 and V<120 km s-1), unique among high z (>2) radio galaxies.

Recent surveys have detected Lya emission from z=4.5-6.5 at luminosities as low as 10^41 erg/s. There is good evidence that low numbers of AGN are among observed faint Lya emitters. Combining these observations with an empirical relation between the intrinsic Lya and B-band luminosities of AGN, we obtain an upper limit on the number density of AGN with absolute magnitudes M_B=[-16,-19] at z=4.5-6.5. These AGN are up to two orders of magnitude fainter than those discovered in the Chandra Deep Field, resulting in the faintest observational constraints to date at these redshifts. At z=4.5, the powerlaw slope of the very faint end of the luminosity function of AGN is shallower than the slope observed at lower redshifts, beta <1.6, at the 98% confidence level. In fact, we find marginal evidence that the luminosity function rises with luminosity, corresponding to a powerlaw slope beta <0, at magnitudes fainter than M_B~-20 (75% confidence level). These results suggest either that accretion onto lower mass black holes is less efficient than onto their more massive counterparts, or that the number of black holes powering AGN with M_B >-20 is lower than expected from the M_BH-sigma relation by one-two orders of magnitude. Extrapolating from reverberation-mapping studies suggests that these black holes would have M_BH=10^6-10^7 Msun. To facilitate the identification of AGN among observed Lya emitters, we derive observational properties of faint AGN in the Lya line, as well as in the X-ray and optical bands.

We model the spectra and surface brightness distributions for the Lyman alpha (Lya) radiation expected from protogalaxies that are caught in the early stages of their assembly. We use the results of a companion paper to characterize the radiation emerging from spherically collapsing gas clouds. We then modify the intrinsic spectra to incorporate the effect of subsequent resonant scattering in the intergalactic medium (IGM). Using these models, we interpret a number of recent observations of extended Lya blobs (LABs) at high redshift. We suggest, based on the angular size, energetics, as well as the relatively shallow surface brightness profiles, and double-peaked spectra, that several of these LABs may be associated with collapsing protogalaxies. We suggest two follow-up observations to diagnose the presence of gas infall. High S/N spectra of LABs should reveal a preferential flattening of the surface brightness profile at the red side of the line. Complementary imaging of the blobs at redshifted Balmer alpha wavelengths should reveal the intrinsic Lya emissivity and allow its separation from radiative transfer effects. We show that Lya scattering by infalling gas can reproduce the observed spectrum of Steidel et al's LAB2 as accurately as a recently proposed outflow model. Finally, we find similar evidence for infall in the spectra of point-like Lyman alpha emitters. The presence of scattering by the infalling gas implies that the intrinsic Lya luminosities, and derived quantities, such as the star-formation rate, in these objects may have been underestimated by about an order of magnitude.

We present Monte Carlo calculations of Lyman alpha (Lya) radiative transfer through collapsing gas clouds, representing proto-galaxies that are caught in the process of their assembly. Such galaxies produce Lya flux over an extended solid angle from a spatially extended Lya emissivity and/or from scattering effects. We study the effect of the gas distribution and kinematics, and of the Lya emissivity profile, on the emergent spectrum and surface brightness distribution. The emergent Lya spectrum is typically double-peaked and asymmetric. In practice, the blue peak is significantly enhanced and the red peak, in most cases, will be undetectable. The resulting effective blueshift, combined with scattering in the intergalactic medium, will render extended Lya emission from collapsing protogalaxies difficult to detect beyond redshift z=4. The surface brightness distribution is typically flat, and a strong wavelength dependence of its slope (with preferential flattening at the red side of the line) would be a robust indication that Lya photons are being generated (rather than just scattered) in a spatially extended region around the galaxy. We also find that for self-ionized clouds whose effective Lya optical depth is less than 10^3, infall and outflow models can produce nearly identical spectra and surface brightness distributions, and are difficult to distinguish from one another. The presence of deuterium with a cosmic abundance may produce a narrow but detectable dip in the spectra of systems with moderate hydrogen column densities, in the range 10^18-10^20 cm^-2. Finally, we present a new analytic solution for the emerging Lya spectrum in the limiting case of a static uniform sphere, extending previous solutions for static plane-parallel slabs.

We explore the possibility of resolving an image of a damped Lyman alpha (DLA) system in absorption against an extended, diffuse background X-ray source. Typical columns of neutral hydrogen in DLAs are high enough to block out up to ~30% of the soft X-ray flux at an observed photon energy of 0.5 keV, and we find that ~ 1% of the area of extended X-ray sources at z > 1 have their 0.5 keV flux reduced by at least 20%. We discuss the observability of such absorption and find that < 2 arcsecond resolution, and > 300 photons per angular resolution element are required in the 0.3-8 keV band for its detection, and in order to distinguish it from intrinsic surface brightness fluctuations. For the surface brightness of the currently known high-redshift extended X-ray sources, this requires an integration time of a few Msec on Chandra. The detection will be within the reach of a routine observation with a next generation X-ray telescope such as XEUS or Generation X.

We discuss currently available observational constraints on the reionization history of the intergalactic medium (IGM), and the extent to which accreting black holes (BHs) can help explain these observations. We show new evidence, based on the combined statistics of Lyman alpha and beta absorption in quasar spectra, that the IGM contains a significant amount of neutral hydrogen, and is experiencing rapid ionization around redshift z=6. However, we argue that quasar BHs, even faint ones that are below the detection thresholds of existing optical surveys, are unlikely to drive the evolution of the neutral fraction around this epoch, because they would over-produce the present-day soft X-ray background. On the other hand, the seeds of the quasar BHs around z=6 likely appeared at much earlier epochs (around z=20), and produced hard ionizing radiation by accretion. These early BHs are promising candidates to account for the high redshift ionization (around z=15) implied by the recent cosmic microwave anisotropy data from WMAP. Using a model for the growth of BHs by accretion and mergers in a hierarchical cosmology, we suggest that the early growth of quasars must include a super-Eddington growth phase, and that, although not yet optically identified, the FIRST radio survey may have already detected several thousand BHs at z>6 with Mbh>10^8 Msun.

A population of black holes (BHs) at high redshifts (z>6) that contributes significantly to the ionization of the intergalactic medium (IGM) would be accompanied by the copious production of hard (>10keV) X-ray photons. The resulting hard X-ray background would redshift and be observed as a present-day soft X-ray background (SXB). Under the hypothesis that BHs are the main producers of reionizing photons in the high-redshift universe, we calculate their contribution to the present-day SXB. Our results, when compared to the unresolved component of the SXB in the range 0.5-2 keV, suggest that accreting BHs (be it luminous quasars or their lower-mass ``miniquasar'' counterparts) did not dominate reionization. Distant miniquasars that produce enough X-rays to only partially ionize the IGM to a level of at most 50% are still allowed, but could be severely constrained by improved determinations of the unresolved component of the SXB.

Absorption between the rest-frame wavelengths of 973 and 1026 Angstroms in quasar spectra arises from two sources (apart from occasional metals): one is due to Lyman-alpha (Lya) absorption by materials at a low redshift, and the other from Lyman-beta (Lyb) at a higher redshift. These two sources of absorption are to a good approximation uncorrelated because of their wide physical separation. Therefore, the two-point correlation of absorption in this region of quasar spectra neatly factorizes into two pieces: the Lyb correlation at high z, and the Lya correlation at low z. The latter can be independently measured from quasar spectra at lower redshifts using current techniques. A simple division then offers a way to statistically separate out the Lyb two-point correlation from the Lya correlation. Several applications of this technique are discussed. First, since the Lyb absorption cross-section is lower than Lya by about a factor of 5, the Lyb forest is a better probe of the intergalactic medium (IGM) at higher redshifts where Lya absorption is often saturated. Second, for the same reason, the Lyb forest allows a better measurement of the equation of state of the IGM at higher overdensities, yielding stronger constraints on its slope when used in conjunction with the Lya forest. Third, models of the Lya forest based on gravitational instability make unique predictions for the Lyb forest, which can be tested against observations. We briefly point out that feedback processes that affect higher density regions but leave low density structure intact may be better constrained by the Lyb forest.

The cosmic ultraviolet (UV) ionizing background impacts the formation of dwarf galaxies in the low-redshift universe (z=3) by suppressing gas infall into galactic halos with circular velocities up to v(circ)=75 km/s. Using a one-dimensional, spherically symmetric hydrodynamics code (Thoul & Weinberg 1995), we examine the effect of an ionizing background on low-mass galaxies forming at high redshifts (z>10). We find that the importance of photoionization feedback is greatly reduced, because (1) at high redshift, dwarf-galaxy sized objects can self-shield against the ionizing background, (2) collisional cooling processes at high redshift are more efficient, (3) the amplitude of the ionizing background at high redshift is lower, and (4) the ionizing radiation turns on when the perturbation that will become the dwarf galaxy has already grown to a substantial overdensity. We find that because of these reasons, gas can collect inside halos with circular velocities as low as v(circ)=10 km/s. This result has important implications for the reionization history of the universe.

We present results of two 21-cm HI surveys performed with the Australia Telescope Compact Array in the nearby Centaurus A and Sculptor galaxy groups. These surveys are sensitive to compact HI clouds and galaxies with HI masses as low as 3E+06 Msun, and are therefore among the most sensitive extragalactic HI surveys to date. The surveys consist of sparsely spaced pointings that sample approximately 2% of the groups' area on the sky. We detected previously known group members, but we found no new HI clouds or galaxies down to the sensitivity limit of the surveys. If the HI mass function had a faint end slope of alpha = 1.5 below M_{HI} = 10^{7.5} Msun in these groups, we would have expected ~3 new objects. Cold dark matter theories of galaxy formation predict the existence of a large number low mass DM sub-halos that might appear as tiny satellites in galaxy groups. Our results support and extend similar conclusions derived from previous HI surveys that a HI rich population of these satellites does not exist.