We have discovered a 300kpc-wide giant Lya nebula centered on the massive galaxy group RO-1001 at z=2.91 in the COSMOS field. Keck Cosmic Web Imager observations reveal three cold gas filaments converging into the center of the potential well of its ~4x10^{13}Msun dark matter halo, hosting 1200Msun/yr of star formation as probed by ALMA and NOEMA observations. The inwards deceleration of the filament and the prevalence of blueshifted components in the Lya spectra support a scenario of gas accretion. The upper limits on AGN activity and overall energetics favor gravity as the primary Lya powering source and infall as the main source of gas flows to the system. Although interpretational difficulties remain, with outflows and likely also recombinations still playing a role, this finding provides arguably an ideal environment to quantitatively test models of cold gas accretion and galaxy feeding inside a massive halo at high redshift.

Lyman-$\alpha$ (Ly$\alpha$) spectra provide insights into the small-scale structure and kinematics of neutral hydrogen (HI) within galaxies as well as the ionization state of the intergalactic medium (IGM). The former defines the intrinsic spectrum of a galaxy, which is modified by the latter. These two effects are degenerate. Using the IllustrisTNG100 simulation, we study the impact of the IGM on Ly$\alpha$ spectral shapes between $z\sim 0$ and $5$. We compute the distribution of the expected Ly$\alpha$ peaks and of the peak asymmetry for different intrinsic spectra, redshifts and large-scale environments. We find that the commonly used averaged transmission curves give an incorrect perception of the observed spectral properties. We show that the distributions of peak counts and asymmetry can lift the degeneracy between the intrinsic spectrum and the IGM absorption. For example, we expect a significant number of triple peaked Ly$\alpha$ spectra (up to 30% at $z\sim 3$) if the galaxies' HI distribution become more porous at higher redshift as predicted by cosmological simulations. We provide a public catalog of transmission curves to be used in future simulations and observations to allow a more realistic IGM treatment.

We present the observations of Lyman continuum (LyC) emission in the afterglow spectra of GRB 191004B at $z=3.5055$, together with those of the other two previously known LyC-emitting long gamma-ray bursts (LGRB) (GRB 050908 at $z=3.3467$, and GRB 060607A at $z=3.0749$), to determine their LyC escape fraction and compare their properties. From the afterglow spectrum of GRB 191004B we determine a neutral hydrogen column density at the LGRB redshift of $\log(N_{\rm HI}/cm^{-2})= 17.2 \pm 0.15$, and negligible extinction ($A_{\rm V}=0.03 \pm 0.02$ mag). The only metal absorption lines detected are CIV and SiIV. In contrast to GRB 050908 and GRB 060607A, the host galaxy of GRB 191004B displays significant Ly$\alpha$ emission. From its Ly$\alpha$ emission and the non-detection of Balmer emission lines we constrain its star-formation rate (SFR) to $1 \leq$ SFR $\leq 4.7$ M$_{\odot}\ yr^{-1}$. We fit the Ly$\alpha$ emission with a shell model and find parameters values consistent with the observed ones. The absolute LyC escape fractions we find for GRB 191004B, GRB 050908 and GRB 060607A are of $0.35^{+0.10}_{-0.11}$, $0.08^{+0.05}_{-0.04}$ and $0.20^{+0.05}_{-0.05}$, respectively. We compare the LyC escape fraction of LGRBs to the values of other LyC emitters found from the literature, showing that LGRB afterglows can be powerful tools to study LyC escape for faint high-redshift star-forming galaxies. Indeed we could push LyC leakage studies to much higher absolute magnitudes. The host galaxies of the three LGRB presented here have all $M_{\rm 1600} > -19.5$ mag, with the GRB 060607A host at $M_{\rm 1600} > -16$ mag. LGRB hosts may therefore be particularly suitable for exploring the ionizing escape fraction in galaxies that are too faint or distant for conventional techniques. Furthermore the time investment is very small compared to galaxy studies. [Abridged]

We demonstrate the redshift-evolution of the spectral profile of HI Lyman-alpha (Ly{\alpha}) emission from star-forming galaxies. In this first study we pay special attention to the contribution of blueshifted emission. At redshift z=2.9-6.6, we compile spectra of a sample of 229 Ly{\alpha}-selected galaxies identified with the Multi-Unit Spectroscopic Explorer at the Very Large Telescope, while at low-z (< 0.44) we use a sample of 74 ultraviolet-selected galaxies observed with the Cosmic Origin Spectrograph onboard the Hubble Space Telescope. At low-z, where absorption from the intergalactic medium (IGM) is negligible, we show that the relative contribution of blueshifted emission (blue/red) increases rapidly with increasing equivalent width ($W_{Ly{\alpha}}$). This result does not, however, emerge at z=3-4, and we use bootstrap simulations to demonstrate the blue emission to be suppressed by stochastic IGM absorption. Our main result is that the blue-peak contribution evolves rapidly downwards with increasing redshift: the blue/red flux ratio is $\approx 30$% at $z\approx 0$, but drops to 15% at $z\approx 3$, and to below 3% by $z\approx 6$. Applying further simulations of the IGM absorption to the unabsorbed COS spectrum, we demonstrate that this decrease in the blue-wing contribution can be entirely attributed to the thickening of intervening Ly{\alpha} absorbing systems, with no need for additional HI opacity from local structure, companion galaxies, or cosmic infall. We discuss our results in light of the numerical radiative transfer simulations, the evolving total Ly{\alpha} output of galaxies, and the utility of resolved Ly{\alpha} spectra in the reionization epoch.

The decline in abundance of Lyman-$\alpha$ (Ly$\alpha$) emitting galaxies at $z \gtrsim 6$ is a powerful and commonly used probe to constrain the progress of cosmic reionization. We use the CoDaII simulation, which is a radiation hydrodynamic simulation featuring a box of $\sim 94$ comoving Mpc side length, to compute the Ly$\alpha$ transmission properties of the intergalactic medium (IGM) at $z\sim 5.8$ to $7$. Our results mainly confirm previous studies, i.e., we find a declining Ly$\alpha$ transmission with redshift and a large sightline-to-sightline variation. However, motivated by the recent discovery of blue Ly$\alpha$ peaks at high redshift, we also analyze the IGM transmission on the blue side, which shows a rapid decline at $z\gtrsim 6$ of the blue transmission. This low transmission can be attributed not only to the presence of neutral regions but also to the residual neutral hydrogen within ionized regions, for which a density even as low as $n_{\rm HI}\sim 10^{-9}\,\mathrm{cm}^{-3}$ (sometimes combined with kinematic effects) leads to a significantly reduced visibility. Still, we find that $\sim 5\%$ of sightlines towards $M_{\mathrm{1600AB}}\sim -21$ galaxies at $z\sim 7$ are transparent enough to allow a transmission of a blue Ly$\alpha$ peak. We discuss our results in the context of the interpretation of observations.

Identifying and characterising reionised bubbles enables us to track both their size distribution, which depends on the primary ionising sources, and the relationship between reionisation and galaxy evolution. We demonstrate that spectrally resolved $z\gtrsim6$ Lyman-alpha (Ly$\alpha$) emission can constrain properties of reionised regions. Specifically, the distant from a source to a neutral region sets the minimum observable Ly$\alpha$ velocity offset from systemic. Detection of flux on the blue side of the Ly$\alpha$ resonance implies the source resides in a large, sufficiently ionised region that photons can escape without significant resonant absorption, and thus constrains both the sizes of and the residual neutral fractions within ionised bubbles. We estimate the extent of the region around galaxies which is optically thin to blue Ly$\alpha$ photons, analogous to quasar proximity zones, as a function of the source's ionising photon output and surrounding gas density. This optically thin region is typically $\lesssim 0.3$ pMpc in radius (allowing transmission of flux $\gtrsim -250$ km s$^{-1}$), $\lesssim 10$% of the distance to the neutral region. In a proof-of-concept, we demonstrate the $z\approx6.6$ galaxy COLA1 - with a blue Ly$\alpha$ peak - likely resides in an ionised region $>0.7$ pMpc, with residual neutral fraction $<10^{-5.5}$. To ionise its own proximity zone we infer COLA1 has a steep UV slope ($\beta < -2.27$), a high ionising photon escape fraction ($f_{\mathrm{esc}}>0.48$) and moderate line-of-sight gas density ($<2.61\times$ the cosmic mean).

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.

Quasar proximity zones at $z>5.5$ correspond to over-dense and over-ionized environments. Galaxies found inside proximity zones can therefore display features which would otherwise be masked by absorption in the IGM. We demonstrate the utility of this quasar-galaxy synergy by reporting the discovery of the first three `proximate Lyman-$\alpha$ emitters' (LAEs) within the proximity zone of quasar J0836 at $z=5.802$ (\textit{Aerith A, B} and \textit{C}). \textit{Aerith A}, located behind the quasar with an impact parameter $D_\perp = 278$ pkpc, provides the first detection of a Lyman-$\alpha$ transverse proximity effect. We model the transmission and show it constrains the onset of J0836's quasar phase to $0.2 \text{Myr}<t<20\text{Myr}$ in the past. The second object, \textit{Aerith B} at a distance $D=750$ pkpc from the quasar, displays a bright, broad double-peaked \lal emission line. Based on relations calibrated at $z\leq3$, the peak separation implies a low ionizing $f_{\text{esc}} \lesssim 1\%$, the most direct such constraint on a reionization-era galaxy. We fit the Ly-$\alpha$ line with an outflowing shell model, finding a completely typical central density $\text{log N}_{\text{HI}}/\text{cm}^{-2} = 19.3_{-0.2}^{+0.8}$, outflow velocity $v=16_{-11}^{+4}$ km s$^{-1}$, and gas temperature $\text{log} T/\text{K} = 3.8_{-0.7}^{+0.8}$ compared to $2<z<3$ analogue LAEs. Finally, we detect an emission line at $\lambda=8177$ \AA\ in object \textit{Aerith C} which, if it is \lal at $z=5.726$, would correspond closely with the end of the quasar's proximity zone ($\Delta z<0.02$ from the boundary) and suggests the quasar influences the IGM up to $\sim85$ cMpc away, making it the largest quasar proximity zone. Via the analyses conducted here, we illustrate how proximate LAEs offer unique insight into the ionizing properties of both quasars and galaxies during the epoch of reionization.

Cold $T \sim 10^{4}$K gas morphology could span a spectrum ranging from large discrete clouds to a fine `mist' in a hot medium. This has myriad implications, including dynamics and survival, radiative transfer, and resolution requirements for cosmological simulations. Here, we use 3D hydrodynamic simulations to study the pressure-driven fragmentation of cooling gas. This is a complex, multi-stage process, with an initial Rayleigh-Taylor unstable contraction phase which seeds perturbations, followed by a rapid, violent expansion leading to the dispersion of small cold gas `droplets' in the vicinity of the gas cloud. Finally, due to turbulent motions, and cooling, these droplets may coagulate. Our results show that a gas cloud `shatters' if it is sufficiently perturbed out of pressure balance ($\delta P/P\sim 1$), and has a large final overdensity $\chi_{\mathrm{f}}\gtrsim 300$, with only a weak dependence on the cloud size. Otherwise, the droplets reassemble back into larger pieces. We discuss our results in the context of thermal instability, and clouds embedded in a shock heated environment.

The observed properties of the Lyman-$\alpha$ (Ly$\alpha$) emission line are a powerful probe of neutral gas in and around galaxies. We present spatially resolved Ly$\alpha$ spectroscopy with VLT/MUSE targeting VR7, a UV-luminous galaxy at $z=6.532$ with moderate Ly$\alpha$ equivalent width (EW$_0\approx38$ {\AA}). These data are combined with deep resolved [CII]$_{\rm 158 \mu m}$ spectroscopy obtained with ALMA and UV imaging from {\it HST} and we also detect UV continuum with MUSE. Ly$\alpha$ emission is clearly detected with S/N $\approx40$ and FWHM of 374 km s$^{-1}$. Ly$\alpha$ and [CII] are similarly extended beyond the UV, with effective radius r$_{\rm eff} = 2.1\pm0.2$ kpc for a single exponential model or r$_{\rm eff, Ly\alpha, halo} = 3.45^{+1.08}_{-0.87}$ kpc when measured jointly with the UV continuum. The Ly$\alpha$ profile is broader and redshifted with respect to the [CII] line (by 213 km s$^{-1}$), but there are spatial variations that are qualitatively similar in both lines and coincide with resolved UV components. This suggests that the emission originates from two components with plausibly different HI column densities. We place VR7 in the context of other galaxies at similar and lower redshift. The Ly$\alpha$ halo scale length is similar at different redshifts and velocity shifts with respect to the systemic are typically smaller. Overall, we find little indications of a more neutral vicinity at higher redshift. This means that the local ($\sim 10$ kpc) neutral gas conditions that determine the observed Ly$\alpha$ properties in VR7 resemble the conditions in post-re-ionisation galaxies.

We present spatially-resolved echelle spectroscopy of an intervening MgII-FeII-MgI absorption-line system detected at $z_{\rm abs}=0.73379$ toward the giant gravitational arc PSZ1 G311.65-18.48. The absorbing gas is associated to an inclined disk-like star-forming galaxy, whose major axis is aligned with the two arc-segments reported here. We probe in absorption the galaxy's extended disk continuously, at $\approx 3$ kpc sampling, from its inner region out to $15\times$ the optical radius. We detect strong ($W_0^{2796}>0.3$ \r{A}) coherent absorption along $13$ independent positions at impact parameters $D=0$--$29$ kpc on one side of the galaxy, and no absorption at $D=28$--$57$ kpc on the opposite side (all de-lensed distances at $z_{\rm abs}$). We show that: (1) the gas distribution is anisotropic; (2) $W_0^{2796}$, $W_0^{2600}$, $W_0^{2852}$, and the ratio $W_0^{2600}\!/W_0^{2796}$, all anti-correlate with $D$; (3) the $W_0^{2796}$-$D$ relation is not cuspy and exhibits significantly less scatter than the quasar-absorber statistics; (4) the absorbing gas is co-rotating with the galaxy out to $D \lesssim 20$ kpc, resembling a `flat' rotation curve, but at $D\gtrsim 20$ kpc velocities decline below the expectations from a 3D disk-model extrapolated from the nebular [OII] emission. These signatures constitute unambiguous evidence for rotating extra-planar diffuse gas, possibly also undergoing enriched accretion at its edge. Arguably, we are witnessing some of the long-sought processes of the baryon cycle in a single distant galaxy expected to be representative of such phenomena.

The existence of fast moving, cold gas ubiquitously observed in galactic winds is theoretically puzzling, since the destruction time of cold gas is much smaller than its acceleration time. In previous work, we showed that cold gas can accelerate to wind speeds and grow in mass if the radiative cooling time of mixed gas is shorter than the cloud destruction time. Here, we study this process in much more detail, and find remarkably robust cloud acceleration and growth in a wide variety of scenarios. Radiative cooling, rather than the Kelvin-Helmholtz instability, enables self-sustaining entrainment of hot gas onto the cloud via cooling-induced pressure gradients. Indeed, growth peaks when the cloud is almost co-moving. The entrainment velocity is of order the cold gas sound speed, and growth is accompanied by cloud pulsations. Growth is also robust to the background wind and initial cloud geometry. In an adiabatic Chevalier-Clegg type wind, for instance, the mass growth rate is constant. Although growth rates are similar with magnetic fields, cloud morphology changes dramatically, with low density, magnetically supported filaments which have a small mass fraction but dominate by volume. This could bias absorption line observations. Cloud growth from entraining and cooling hot gas can potentially account for the cold gas content of the CGM. It can also fuel star formation in the disk as cold gas recycled in a galactic fountain accretes and cools halo gas. We speculate that galaxy-scale simulations should converge in cold gas mass once cloud column densities of ${\rm N} \sim 10^{18} \ {\rm cm^{-2}}$ are resolved.

During the epoch of reionisation, neutral gas in the early Universe was ionized by hard ultraviolet radiation emitted by young stars in the first galaxies. To do so, ionizing ultraviolet photons must escape from the host galaxy. We present Hubble Space Telescope observations of the gravitationally lensed galaxy PSZ1-ARC G311.6602-18.4624, revealing bright, multiply-imaged ionizing photon escape from a compact star-forming region through a narrow channel in an optically thick gas. The gravitational lensing magnification shows how ionizing photons escape this galaxy, contributing to the re-ionization of the Universe. The multiple sight lines to the source probe absorption by intergalactic neutral hydrogen on scales of no more than a few hundred, perhaps even less than ten, parsec.

We examine Lyman continuum (LyC) leakage through H II regions regulated by turbulence and radiative feedback in a giant molecular cloud in the context of fully-coupled radiation hydrodynamics (RHD). The physical relations of the LyC escape with H I covering fraction, kinematics, spectral hardness, and the emergent Lyman-$\alpha$ (Ly$\alpha$) line profiles are studied using a series of RHD turbulence simulations performed with RAMSES-RT. The turbulence-regulated mechanism allows ionizing photons to leak out at early times before the onset of supernova feedback. The LyC photons escape through turbulence-generated low column density channels which are evacuated efficiently by radiative feedback via photoheating-induced shocks across the D-type ionization fronts. Ly$\alpha$ photons funnel through the photoionized channels along the paths of LyC escape, resulting in a diverse Ly$\alpha$ spectral morphology including narrow double-peaked profiles. The Ly$\alpha$ peak separation is controlled by the residual H I column density of the channels and the line asymmetry correlates with the porosity and multiphase structure of the H II region. This mechanism through the turbulent H II regions can naturally reproduce the observed Ly$\alpha$ spectral characteristics of some of LyC-leaking galaxies. This RHD turbulence-origin provides an appealing hypothesis to explain high LyC leakage from very young ($\sim3$ Myr) star-forming galaxies found in the local Universe without need of extreme galactic outflows nor supernova feedback. We discuss the implications of the turbulent H II regions on other nebular emission lines and a possible observational test with the Magellanic System and local blue compact dwarf galaxies as analogs of reionization-era systems.

Searches for modified gravity in the large-scale structure try to detect the enhanced amplitude of density fluctuations caused by the fifth force present in many of these theories. Neutrinos, on the other hand, suppress structure growth below their free-streaming length. Both effects take place on comparable scales, and uncertainty in the neutrino mass leads to a degeneracy with modified gravity parameters for probes that are measuring the amplitude of the matter power spectrum. We explore the possibility to break the degeneracy between modified gravity and neutrino effects in the growth of structures by considering kinematic information related to either the growth rate on large scales or the virial velocities inside of collapsed structures. In order to study the degeneracy up to fully non-linear scales, we employ a suite of $N$-body simulations including both $f(R)$ modified gravity and massive neutrinos. Our results indicate that velocity information provides an excellent tool to distinguish massive neutrinos from modified gravity. Models with different values of neutrino masses and modified gravity parameters possessing a comparable matter power spectrum at a given time have different growth rates. This leaves imprints in the velocity divergence, which is therefore better suited than the amplitude of density fluctuations to tell the models apart. In such models with a power spectrum comparable to $\Lambda$CDM today, the growth rate is strictly enhanced. We also find the velocity dispersion of virialised clusters to be well suited to constrain deviations from general relativity without being affected by the uncertainty in the sum of neutrino masses.

Multiple modifications of general relativity (GR) have been proposed in the literature in order to understand the nature of the accelerated expansion of the Universe. However, thus far all the predictions of GR have been confirmed with constantly increasing accuracy. In this work, we study the imprints of a particular class of models -- "screened" modified gravity theories -- on cosmic filaments. We have utilized the $N$-body code ISIS/RAMSES to simulate the symmetron model and the Hu-Sawicky $f(R)$ model, and we post-process the output with DisPerSE to identify the filaments of the cosmic web. We investigated how the global properties of the filaments -- such as their lengths, masses, and thicknesses -- as well as their radial density and speed profiles change under different gravity theories. We find that filaments are, on average, shorter and denser in modified gravity models compared to in $\Lambda$CDM. We also find that the speed profiles of the filaments are enhanced, consistent with theoretical expectations. Overall, our results suggest that cosmic filaments can be an effective complementary probe of screened modified gravity theories on Mpc scales.

Distant luminous Lyman-alpha emitters (LAEs) are excellent targets for spectroscopic observations of galaxies in the epoch of reionisation (EoR). We present deep high-resolution (R=5000) VLT/X-SHOOTER observations, along with an extensive collection of photometric data of `COLA1', a proposed double peaked LAE at z=6.6 (Hu et al. 2016). We rule out that COLA1's emission line is an [OII] doublet at z=1.475 on the basis of i) the asymmetric red line-profile and flux ratio of the peaks (blue/red=$0.31\pm0.03$) and ii) an unphysical [OII]/Ha ratio ([OII]/Ha > 22). We show that COLA1's observed B-band flux is explained by a faint extended foreground LAE, for which we detect Lya and [OIII] at z=2.142. We thus conclude that COLA1 is a real double-peaked LAE at z=6.593, the first discovered at z>6, confirming the result from Hu et al. (2016). COLA1 is UV luminous (M$_{1500}=-21.6\pm0.3$), has a high equivalent width (EW$_{0}$~120 \AA) and very compact Lya emission (r$_{50} = 0.3$ kpc). Relatively weak inferred Hb+[OIII] line-emission from Spitzer/IRAC indicates an extremely low metallicity of Z<1/20 Z$_{\odot}$ or reduced strength of nebular lines due to high escape of ionising photons. The small Lya peak separation of $220\pm20$ km/s implies a low HI column density and an ionising photon escape fraction of ~15-30 %, providing the first direct evidence that such galaxies contribute actively to the reionisation of the Universe at z>6. Based on simple estimates, we find that COLA1 could have provided just enough photons to reionise its own ~0.3 pMpc bubble, allowing the blue Lya line to be observed. However, we also discuss alternative scenarios explaining the detected double peaked nature of COLA1. Our results show that future high-resolution observations of statistical samples of double peaked LAEs at z>5 are a promising probe of the occurrence of ionised regions around galaxies in the EoR.

Both absorption and emission line studies show that cold gas around galaxies is commonly outflowing at speeds of several hundred km$\,\textrm{s}^{-1}$. This observational fact poses a severe challenge to our theoretical models of galaxy evolution since most feedback mechanisms (e.g., supernovae feedback) accelerate hot gas, and the timescale it takes to accelerate a blob of cold gas via a hot wind is much larger than the time it takes to destroy the blob. We revisit this long-standing problem using three-dimensional hydrodynamical simulations with radiative cooling. Our results confirm previous findings, that cooling is often not efficient enough to prevent the destruction of cold gas. However, we also identify regions of parameter space where the cooling efficiency of the mixed, `warm' gas is sufficiently large to contribute new comoving cold gas which can significantly exceed the original cold gas mass. This happens whenever, $t_{\mathrm{cool, mix}}/t_{\mathrm{cc}} < 1$, where $t_{\mathrm{cool,mix}}$ is the cooling time of the mixed warm gas and $t_{\mathrm{cc}}$ is the cloud-crushing time. This criterion is always satisfied for a large enough cloud. Cooling `focuses' stripped material onto the tail where mixing takes place and new cold gas forms. A sufficiently large simulation domain is crucial to capturing this behavior.

Recent magneto-hydrodynamic simulations of the star-forming interstellar medium (ISM) with parsec scale resolution indicate that relativistic cosmic rays support the launching of galactic outflows on scales of a few kpc. If these fundamental constituents of the ISM are injected at the sites of supernova (SN) explosions, the outflows are smoother, colder, and denser than the highly structured, hot-phase driven outflows forming, e.g., by thermal SN energy injection alone. In this Letter we present computations of resonant Lyman-$\alpha$ (Ly$\alpha$) radiation transfer through snapshots of a suite of stratified disk simulations from the SILCC project. For a range of thermal, radiative, and kinetic feedback models only simulations including non-thermal cosmic rays produce Ly$\alpha$ spectra with enhanced red peaks and strong absorption at line center -- similar to observed systems. The absence of cosmic ray feedback leads to spectra incompatible with observations. We attribute this to the smoother neutral gas distribution of cosmic ray supported outflows within a few kpc from the disk midplane.

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.

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 rest-frame ultraviolet and optical spectroscopy of the brightest lensed galaxy yet discovered, at redshift z = 2.4. This source reveals a characteristic, triple-peaked Lyman {\alpha} profile which has been predicted by various theoretical works but to our knowledge has not been unambiguously observed previously. The feature is well fit by a superposition of two components: a double-peak profile emerging from substantial radiative transfer, and a narrow, central component resulting from directly escaping Lyman {\alpha} photons; but is poorly fit by either component alone. We demonstrate that the feature is unlikely to contain contamination from nearby sources, and that the central peak is unaffected by radiative transfer effects apart from very slight absorption. The feature is detected at signal-to-noise ratios exceeding 80 per pixel at line center, and bears strong resemblance to synthetic profiles predicted by numerical models.

We examine the dust geometry and Ly{\alpha} scattering in the galaxies of the Lyman Alpha Reference Sample (LARS), a set of 14 nearby (0.02 < $z$ < 0.2) Ly{\alpha} emitting and starbursting systems with Hubble Space Telescope Ly{\alpha}, H{\alpha}, and H{\beta} imaging. We find that the global dust properties determined by line ratios are consistent with other studies, with some of the LARS galaxies exhibiting clumpy dust media while others of them show significantly lower Ly{\alpha} emission compared to their Balmer decrement. With the LARS imaging, we present Ly{\alpha}/H{\alpha} and H{\alpha}/H{\beta} maps with spatial resolutions as low as $\sim$ 40 pc, and use these data to show that in most galaxies, the dust geometry is best modeled by three distinct regions: a central core where dust acts as a screen, an annulus where dust is distributed in clumps, and an outer envelope where Ly{\alpha} photons only scatter. We show that the dust that affects the escape of Ly{\alpha} is more restricted to the galaxies' central regions, while the larger Ly{\alpha} halos are generated by scattering at large radii. We present an empirical modeling technique to quantify how much Ly{\alpha} scatters in the halo, and find that this "characteristic" scattering distance correlates with the measured size of the Ly{\alpha} halo. We note that there exists a slight anti-correlation between the scattering distance of Ly{\alpha} and global dust properties.

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 compare $237$ Lyman-$\alpha$ (Ly$\alpha$) spectra of the "MUSE-Wide survey" (Herenz et al. 2017) to a suite of radiative transfer simulations consisting of a central luminous source within a concentric, moving shell of neutral gas, and dust. This six parameter shell-model has been used numerously in previous studies, however, on significantly smaller data-sets. We find that the shell-model can reproduce the observed spectral shape very well - better than the also common `Gaussian-minus-Gaussian' model which we also fitted to the dataset. Specifically, we find that $\sim 94\%$ of the fits possess a goodness-of-fit value of $p(\chi^2)>0.1$. The large number of spectra allows us to robustly characterize the shell-model parameter range, and consequently, the spectral shapes typical for realistic spectra. We find that the vast majority of the Ly$\alpha$ spectral shapes require an outflow and only $\sim 5\%$ are well-fitted through an inflowing shell. In addition, we find $\sim 46\%$ of the spectra to be consistent with a neutral hydrogen column density $<10^{17}\,\mathrm{cm}^{-2}$ - suggestive of a non-negligible fraction of continuum leakers in the MUSE-Wide sample. Furthermore, we correlate the spectral against the Ly$\alpha$ halo properties against each other but do not find any strong correlation.

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

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.

We use cosmological hydrodynamical simulations to study the effect of screened modified gravity models on the mass estimates of galaxy clusters. In particular, we focus on two novel aspects: (i) we study modified gravity models in which baryons and dark matter are coupled with different strengths to the scalar field, and, (ii) we put the simulation results into the greater context of a general screened-modified gravity parametrization. We have compared the mass of clusters inferred via lensing versus the mass inferred via kinematical measurements as a probe of violations of the equivalence principle at Mpc scales. We find that estimates of cluster masses via X-ray observations is mainly sensitive to the coupling between the scalar degree of freedom and baryons - while the kinematical mass is mainly sensitive to the coupling to dark matter. Therefore, the relation between the two mass estimates is a probe of a possible non-universal coupling between the scalar field, the standard model fields, and dark matter. Finally, we used observational data of kinetic, thermal and lensing masses to place constraints on deviations from general relativity on cluster scales for a general parametrization of screened modified gravity theories which contains $f(R)$ and Symmetron models. We find that while the kinematic mass can be used to place competitive constraints, using thermal measurements is challenging as a potential non-thermal contribution is degenerate with the imprint of modified gravity.

Several `giant' Lyman-$\alpha$ (Ly$\alpha$) nebulae with extent $\gtrsim 300\,$kpc and observed Ly$\alpha$ luminosity of $\gtrsim 10^{44}\,{\rm erg}\,{\rm s}^{-1}\,{\rm cm}^{-2}\,{\rm arcsec}^{-2}$ have recently been detected, and it has been speculated that their presence hints at a substantial cold gas reservoir in small cool clumps not resolved in modern hydro-dynamical simulations. We use the Illustris simulation to predict the Ly$\alpha$ emission emerging from large halos ($M > 10^{11.5}M_{\odot}$) at $z\sim 2$ and thus test this model. We consider both AGN and star driven ionization, and compare the simulated surface brightness maps, profiles and Ly$\alpha$ spectra to a model where most gas is clumped below the simulation resolution scale. We find that with Illustris no additional clumping is necessary to explain the extents, luminosities and surface brightness profiles of the `giant Ly$\alpha$ nebulae' observed. Furthermore, the maximal extents of the objects show a wide spread for a given luminosity and do not correlate significantly with any halo properties. We also show how the detected size depends strongly on the employed surface brightness cutoff, and predict that further such objects will be found in the near future.

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.

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.

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

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.

Modified gravity models require a screening mechanism to be able to evade the stringent constraints from local gravity experiments and, at the same time, give rise to observable astrophysical and cosmological signatures. Such screened modified gravity models necessarily have dynamics determined by complex nonlinear equations that usually need to be solved on a model-by-model basis to produce predictions. This makes testing them a cumbersome process. In this paper, we investigate whether there is a common signature for all the different models that is suitable to testing them on cluster scales. To do this we propose an observable related to the fifth force, which can be observationally related to the ratio of dynamical-to-lensing mass of a halo, and then show that the predictions for this observable can be rescaled to a near universal form for a large class of modified gravity models. We demonstrate this using the Hu-Sawicki $f(R)$, the Symmetron, the nDGP, and the Dilaton models, as well as unifying parametrizations. The universal form is determined by only three quantities: a strength, a mass, and a width parameter. We also show how these parameters can be derived from a specific theory. This self-similarity in the predictions can hopefully be used to search for signatures of modified gravity on cluster scales in a model-independent way.

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

Screened modified gravity predicts potentially large signatures in the peculiar velocity field that makes it an interesting probe to test gravity on cosmological scales. We investigate the signatures induced by the Symmetron and a Chameleon $f(R)$ model in the peculiar velocity field using $N$-body simulations. By studying fifth force and halo velocity profiles we identify three general categories of effects found in screened modified gravity models: a fully screened regime where we recover $\Lambda$CDM to high precision, an unscreened regime where the fifth force is in full operation, and, a partially screened regime where screening occurs in the inner part of a halo, but the fifth force is active at larger radii. These three regimes can be pointed out very clearly by analyzing the deviation in the maximum cluster velocity. Observationally, the partially screened regime is of particular interest since an uniform increase of the gravitational force - as present in the unscreened regime - is degenerate with the (dynamical) halo mass estimate, and, thus, hard to detect.

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.

Aims. We investigate the gravitational redshift in clusters of galaxies in the symmetron and Hu-Sawicky $f(R)$ models. The characteristic feature of both models is the screening mechanism that hides the fifth force in dense environments recovering general relativity. Methods. We use $N$-body simulations that were run using the code \texttt{Isis}, which includes scalar fields, to analyse the deviation of observables in modified gravity models with respect to $\Lambda$CDM. Results. We find that due to the presence of the screening, the deviation is highly dependent on the halo mass. For instance, the $f(R)$ parameters $|f_{R0}|=10^{-5}$, $n=1$ cause an enhancement of the gravitational signal by up to 50% for halos with masses between $10^{13}\,M_\odot h^{-1}$ and $10^{14}\,M_\odot h^{-1}$. The characteristic mass range where the fifth force is most active varies with the model parameters. The usual assumption is that the presence of a fifth force leads to a deeper potential well and thus, a stronger gravitational redshift. However, we find that in cases in which only the central regions of the halos are screened, there could also be a weaker gravitational redshift.