We investigate the ultraviolet (UV) spectral properties of faint Lyman-$\alpha$ emitters (LAEs) in the redshift range 2.9<z<4.6 and provide material to prepare future observations of the faint Universe. We use data from the MUSE Hubble Ultra Deep Survey to construct mean rest-frame spectra of continuum-faint (median M$_{UV}$ of -18 and down to M$_{UV}$ of -16), low stellar mass (median value of $10^{8.4}$ and down to $10^{7}M_{\odot}$) LAEs at redshift z>3. We compute various averaged spectra of LAEs sub-sampled on the basis of their observational (e.g., Ly$\alpha$ strength, UV magnitude and spectral slope) and physical (e.g., stellar mass and star-formation rate) properties. We search for UV spectral features other than Ly$\alpha$, such as higher-ionization nebular emission lines and absorption features. We successfully observe the OIII]1666 and CIII]909 collisionally excited emission lines and the HeII1640 recombination feature, as well as the resonant CIV1550 doublet either in emission or P-Cygni. We compare the observed spectral properties of the different mean spectra and find the emission lines to vary with the observational and physical properties of the LAEs. In particular, the mean spectra of LAEs with larger Ly$\alpha$ equivalent widths, fainter UV magnitudes, bluer UV spectral slopes and lower stellar masses show the strongest nebular emission. The line ratios of these lines are similar to those measured in the spectra of local metal-poor galaxies, while their equivalent widths are weaker compared to the handful of extreme values detected in individual spectra of z>2 galaxies. This suggests that weak UV features are likely ubiquitous in high z, low-mass and faint LAEs. We publicly release the stacked spectra as they can serve as empirical templates for the design of future observations, such as those with the James Webb Space Telescope and the Extremely Large Telescope.

We present initial results from the Cosmic Ultraviolet Baryon Survey (CUBS). CUBS is designed to map diffuse baryonic structures at redshift z<~1 using absorption-line spectroscopy of 15 UV-bright QSOs with matching deep galaxy survey data. CUBS QSOs are selected based on their NUV brightness to avoid biases against the presence of intervening Lyman Limit Systems (LLSs) at zabs<1. We report five new LLSs of log N(HI)/cm^-2 >~ 17.2 over a total redshift survey pathlength of dz=9.3, and a number density of n(z)=0.43 (-0.18, +0.26). Considering all absorbers with log N(HI)/cm^-2 > 16.5 leads to n(z)=1.08 (-0.25, +0.31) at z<1. All LLSs exhibit a multi-component structure and associated metal transitions from multiple ionization states such as CII, CIII, MgII, SiII, SiIII, and OVI absorption. Differential chemical enrichment levels as well as ionization states are directly observed across individual components in three LLSs. We present deep galaxy survey data obtained using the VLT-MUSE integral field spectrograph and the Magellan Telescopes, reaching sensitivities necessary for detecting galaxies fainter than 0.1L* at d<~300 physical kpc (pkpc) in all five fields. A diverse range of galaxy properties is seen around these LLSs, from a low-mass dwarf galaxy pair, a co-rotating gaseous halo/disk, a star-forming galaxy, a massive quiescent galaxy, to a galaxy group. The closest galaxies have projected distances ranging from d=15 to 72 pkpc and intrinsic luminosities from ~0.01L* to ~3L*. Our study shows that LLSs originate in a variety of galaxy environments and trace gaseous structures with a broad range of metallicities.

A prediction of the classic active galactic nuclei (AGN) unification model is the presence of ionisation cones with different orientations depending on the AGN type. Confirmations of this model exist for present times, but it is less clear in the early Universe. Here, we use the morphology of giant Ly$\alpha$ nebulae around AGNs at redshift z$\sim$3 to probe AGN emission and therefore the validity of the AGN unification model at this redshift. We compare the spatial morphology of 19 nebulae previously found around type I AGNs with a new sample of 4 Ly$\alpha$ nebulae detected around type II AGNs. Using two independent techniques, we find that nebulae around type II AGNs are more asymmetric than around type I, at least at radial distances $r>30$~physical kpc (pkpc) from the ionizing source. We conclude that the type I and type II AGNs in our sample show evidence of different surrounding ionising geometries. This suggests that the classical AGN unification model is also valid for high-redshift sources. Finally, we discuss how the lack of asymmetry in the inner parts (r$\lesssim$30 pkpc) and the associated high values of the HeII to Ly$\alpha$ ratios in these regions could indicate additional sources of (hard) ionizing radiation originating within or in proximity of the AGN host galaxies. This work demonstrates that the morphologies of giant Ly$\alpha$ nebulae can be used to understand and study the geometry of high redshift AGNs on circum-nuclear scales and it lays the foundation for future studies using much larger statistical samples.

We present an IGM HI tomography map in a survey volume of $16 \times 19 \times 131 \ h^{-3} {\rm comoving \ Mpc}^{3}$ (cMpc$^3$) centered at MAMMOTH-1 nebula and three neighbouring quasars at $z=2.3$. MAMMOTH-1 nebula is an enormous Ly$\alpha$ nebula (ELAN), hosted by a type-II quasar dubbed MAMMOTH1-QSO, that extends over $1\ h^{-1}$ cMpc with not fully clear physical origin. Here we investigate the HI-gas distribution around MAMMOTH1-QSO with the ELAN and three neighbouring type-I quasars, making the IGM HI tomography map with a spatial resolution of $2.6\ h^{-1}$ cMpc. Our HI tomography map is reconstructed with HI Ly$\alpha$ forest absorption of bright background objects at $z=2.4-2.9$: one eBOSS quasar and 16 Keck/LRIS galaxy spectra. We estimate the radial profile of HI flux overdensity for MAMMOTH1-QSO, and find that MAMMOTH1-QSO resides in a volume with significantly weak HI absorption. This suggests that MAMMOTH1-QSO has a proximity zone where quasar illuminates and photo-ionizes the surrounding HI gas and suppresses HI absorption, and that the ELAN is probably a photo-ionized cloud embedded in the cosmic web. The HI radial profile of MAMMOTH1-QSO is very similar to those of three neighbouring type-I quasars at $z=2.3$, which is compatible with the AGN unification model. We compare the distributions of the HI absorption and star-forming galaxies in our survey volume, and identify a spatial offset between density peaks of star-forming galaxies and HI gas. This segregation may suggest anisotropic UV background radiation created by star-forming galaxy density fluctuations.

While low-luminosity galaxies dominate number counts at all redshifts, their contribution to cosmic Reionization is poorly understood due to a lack of knowledge of their physical properties. We isolate a sample of 35 z~4-5 continuum-faint Lyman-alpha emitters from deep VLT/MUSE spectroscopy and directly measure their Halpha emission using stacked Spitzer/IRAC Ch. 1 photometry. Based on Hubble Space Telescope imaging, we determine that the average UV continuum magnitude is fainter than -16 (~0.01 L_star), implying a median Lyman-alpha equivalent width of 249 Angstroms. By combining the Halpha measurement with the UV magnitude we determine the ionizing photon production efficiency, xi_ion, a first for such faint galaxies. The measurement of log (xi_ion [Hz/erg]) = 26.28 (+0.28; -0.40) is in excess of literature measurements of both continuum- and emission line-selected samples, implying a more efficient production of ionizing photons in these lower-luminosity, Lyman-alpha-selected systems. We conclude that this elevated efficiency can be explained by stellar populations with metallicities between 4e-4 and 0.008, with light-weighted ages less than 3 Myr.

We present spatially resolved maps of six individually-detected Lyman alpha haloes (LAHs) as well as a first statistical analysis of the Lyman alpha (Lya) spectral signature in the circum-galactic medium of high-redshift star-forming galaxies using MUSE. Our resolved spectroscopic analysis of the LAHs reveals significant intrahalo variations of the Lya line profile. Using a three-dimensional two-component model for the Lya emission, we measure the full width at half maximum (FWHM), the peak velocity shift and the asymmetry of the Lya line in the core and in the halo of 19 galaxies. We find that the Lya line shape is statistically different in the halo compared to the core for ~40% of our galaxies. Similarly to object-by-object based studies and a recent resolved study using lensing, we find a correlation between the peak velocity shift and the width of the Lya line both at the interstellar and circum-galactic scales. While there is a lack of correlation between the spectral properties and the spatial scale lengths of our LAHs, we find a correlation between the width of the line in the LAH and the halo flux fraction. Interestingly, UV bright galaxies show broader, more redshifted and less asymmetric Lya lines in their haloes. The most significant correlation found is for the FWHM of the line and the UV continuum slope of the galaxy, suggesting that the redder galaxies have broader Lya lines. The generally broad and red line shapes found in the halo component suggests that the Lya haloes are powered either by scattering processes through an outflowing medium, fluorescent emission from outflowing cold clumps of gas, or a mix of both. Considering the large diversity of the Lya line profiles observed in our sample and the lack of strong correlation, the interpretation of our results is still broadly open and underlines the need for realistic spatially resolved models of the LAHs.

We present a study of the galaxy environment of 9 strong HI+CIV absorption line systems ($16.2<{\rm log}(N({\rm HI}))<21.2$) spanning a wide range in metallicity at $z\sim4-5$, using MUSE integral field and X-Shooter spectroscopic data collected in a $z\approx 5.26$ quasar field. We identify galaxies within a 250 kpc and $\pm1000$ km s$^{-1}$ window for 6 out of the 9 absorption systems, with 2 of the absorption line systems showing multiple associated galaxies within the MUSE field of view. The space density of Ly$\alpha$ emitting galaxies (LAEs) around the HI and CIV systems is $\approx10-20$ times the average sky density of LAEs given the flux limit of our survey, showing a clear correlation between the absorption and galaxy populations. Further, we find that the strongest CIV systems in our sample are those that are most closely aligned with galaxies in velocity space, i.e. within velocities of $\pm500$ km s$^{-1}$. The two most metal poor systems lie in the most dense galaxy environments, implying we are potentially tracing gas that is infalling for the first time into star-forming groups at high redshift. Finally, we detect an extended Ly$\alpha$ nebula around the $z\approx 5.26$ quasar, which extends up to $\approx50$ kpc at the surface brightness limit of $3.8 \times 10^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$. After scaling for surface brightness dimming, we find that this nebula is centrally brighter, having a steeper radial profile than the average for nebulae studied at $z\sim3$ and is consistent with the mild redshift evolution seen from $z\approx 2$.

During the last years, Ly$\alpha$ nebulae have been routinely detected around high-z, radio-quiet quasars (RQQs) thanks to the advent of sensitive integral field spectrographs. Constraining the physical properties of the Ly$\alpha$ nebulae is crucial for a full understanding of the circum-galactic medium (CGM), which is a venue of feeding and feedback processes. The most luminous quasars are privileged test-beds to study these processes, given their large ionizing fluxes and dense CGM environments in which they are expected to be embedded. We aim at characterizing the rest-frame UV emission lines in the CGM around a hyper-luminous, broad emission line, RQQ at z~3.6, that exhibits powerful outflows at both nuclear and host galaxy scales. We analyze VLT/MUSE observations of the quasar J1538+08 and perform a search for extended UV emission lines to characterize its morphology, emissivity, kinematics and metal content. We report the discovery of a very luminous ($\sim2 \times10^{44}~erg~s^{-1}$), giant Ly$\alpha$ nebula and a likely associated extended CIV nebula. The Ly$\alpha$ nebula emission exhibits moderate blueshift compared with the quasar systemic redshift and large average velocity dispersion ($\sigma_{v}$ ~700 $km~s^{-1}$) across the nebula, while the CIV nebula shows $\sigma_{v}$~$350~km~s^{-1}$. The Ly$\alpha$ line profile exhibits a significant asymmetry towards negative velocity values at 20-30 kpc south of the quasar and is well parameterized by two Gaussian components: a narrow ($\sigma$~$470~km~s^{-1}$) systemic one plus a broad ($\sigma$~1200 $km~s^{-1}$), blueshifted (~1500 $km~s^{-1}$) one. Our analysis of the MUSE observation of J1538+08 reveals metal-enriched CGM around this hyper-luminous quasar and our detection of blueshifted emission in the line profile of the Ly$\alpha$ nebula suggests that powerful nuclear outflows can propagate through the CGM over tens of kpc.

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 the design, methods, and first results of the MUSE Analysis of Gas around Galaxies (MAGG) survey, a large programme on the Multi Unit Spectroscopic Explorer (MUSE) instrument at the Very Large Telescope (VLT) which targets 28 z > 3.2 quasars to investigate the connection between optically-thick gas and galaxies at z~3-4. MAGG maps the environment of 52 strong absorption line systems at z > 3, providing the first statistical sample of galaxies associated with gas-rich structures in the early Universe. In this paper, we study the galaxy population around a very metal poor gas cloud at z~3.5 towards the quasar J124957.23-015928.8. We detect three Lyman alpha emitters within <200km/s of the cloud redshift, at projected separations <185 kpc (physical). The presence of star-forming galaxies near a very metal-poor cloud indicates that metal enrichment is still spatially inhomogeneous at this redshift. Based on its very low metallicity and the presence of nearby galaxies, we propose that the most likely scenario for this LLS is that it lies within a filament which may be accreting onto a nearby galaxy. Taken together with the small number of other LLSs studied with MUSE, the observations to date show a range of different environments near strong absorption systems. The full MAGG survey will significantly expand this sample and enable a statistical analysis of the link between gas and galaxies to pin down the origin of these diverse environments at z~3-4.

Lyman-\alpha (Ly\alpha) emission lines are typically found to be redshifted with respect to the systemic redshifts of galaxies, likely due to resonant scattering of Ly\alpha photons. Here we measure the average velocity offset for a sample of 96 z ~ 3.3 Ly\alpha emitters (LAEs) with a median Ly\alpha flux (luminosity) of ~10^{-17} erg cm^{-2} s^{-1} (~10^{42} erg s^{-1}) and a median unobscured star formation rate (SFR) of ~1.3 Msun/yr, detected by the Multi-Unit Spectroscopic Explorer as part of our MUSEQuBES circumgalactic medium (CGM) survey. By postulating that the stacked CGM absorption profiles of these LAEs, probed by 8 background quasars, must be centered on the systemic redshift, we measure an average velocity offset, Voffset = 169 \pm 10 km/s, between the Ly\alpha emission peak and the systemic redshift. The observed Voffset is lower by factors of ~1.4 and ~2.6 compared to the velocity offsets measured for narrow-band selected LAEs and Lyman break galaxies, respectively, which probe galaxies with higher masses and SFRs. Consistent with earlier studies based on direct measurements for individual objects, we find that the Voffset is correlated with the full width at half-maximum of the red peak of the Ly\alpha line, and anti-correlated with the rest-frame equivalent width. Moreover, we find that Voffset is correlated with SFR with a sub-linear scaling relation, Voffset ~ SFR^(0.16 \pm 0.03). Adopting the mass scaling for main sequence galaxies, such a relation suggests that Voffset scales with the circular velocity of the dark matter halos hosting the LAEs.

Cosmological simulations predict the Universe contains a network of intergalactic gas filaments, within which galaxies form and evolve. However, the faintness of any emission from these filaments has limited tests of this prediction. We report the detection of rest-frame ultraviolet Lyman-alpha radiation from multiple filaments extending more than one megaparsec between galaxies within the SSA 22 proto-cluster at a redshift of 3.1. Intense star formation and supermassive black-hole activity is occurring within the galaxies embedded in these structures, which are the likely sources of the elevated ionizing radiation powering the observed Lyman-alpha emission. Our observations map the gas in filamentary structures of the type thought to fuel the growth of galaxies and black holes in massive proto-clusters.

Motivated by the recent discovery of the near-ubiquity of Ly$\alpha$ emission around $z \gtrsim 3$ QSOs, we performed a systematic study of QSO circumgalactic Ly$\alpha$ emission at $z\approx2$, utilizing the unique capability of the Keck Cosmic Web Imager (KCWI) -- a new wide-field, blue sensitive integral-field spectrograph (IFU). In this paper, we present KCWI observations on a sample of 16 ultraluminous Type-I QSOs at $z=2.1-2.3$ with ionizing luminosities of $L_{\rm{\nu_{\rm{LL}}}}=10^{31.1-32.3}$ erg s$^{-1}$ Hz$^{-1}$. We found that 14 out of 16 QSOs are associated with Ly$\alpha$ nebulae with projected linear-sizes larger than 50 physical kpc (pkpc). Among them, four nebulae have enormous Ly$\alpha$ emission with the Ly$\alpha$ surface brightness $SB_{\rm{Ly\alpha}} >10^{-17}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ on the $>100$ kpc scale, extending beyond the field of view of KCWI. Our KCWI observations reveal that most $z\approx2$ QSO nebulae have a more irregular morphology compared to those at $z\gtrsim3$. In turn, we measure that the circularly-averaged surface brightness (SB) at $z\approx2$ is 0.4 dex fainter than the redshift-corrected, median SB at $z\gtrsim 3$. The Ly$\alpha$ SB profile (SB$_{\rm{Ly\alpha}}$) of QSOs at $z\approx 2$ can be described by a power law of SB$_{\rm{Ly\alpha},z\approx2.3}= 3.7\times10^{-17}\times(r/40)^{-1.8}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$, with the slope similar to that at $z\gtrsim3$. The observed lower redshift-corrected, circularly-averaged SB may be mainly due to the lower covering factor of cool gas clouds in massive halos at $z\approx2$.

We present the goals, design, and first results of the MUSE Ultra Deep Field (MUDF) survey, a large programme using the Multi Unit Spectroscopic Explorer (MUSE) instrument at the ESO Very Large Telescope. The MUDF survey is collecting ~ 150 hours on-source of integral field optical spectroscopy in a 1.5 x 1.2 square arcmin region which hosts several astrophysical structures along the line of sight, including two bright z ~ 3.2 quasars with close separation (~ 500 kpc). Following the description of the data reduction procedures, we present the analysis of the galaxy environment and gaseous properties of seven groups detected at redshifts 0.5 < z < 1.5, spanning a large dynamic range in halo mass, log(Mh/Msun) ~ 11 - 13.5. For four of the groups, we find associated MgII absorbers tracing cool gas in high-resolution spectroscopy of the two quasars, including one case of correlated absorption in both sightlines at distance ~ 480 kpc. The absorption strength associated with the groups is higher than what has been reported for more isolated galaxies of comparable mass and impact parameters. We do not find evidence for widespread cool gas giving rise to strong absorption within these groups. Combining these results with the distribution of neutral and ionised gas seen in emission in lower-redshift groups, we conclude that gravitational interactions in the group environment strip gas from the galaxy haloes into the intragroup medium, boosting the cross section of cool gas and leading to the high fraction of strong MgII absorbers that we detect.

The well-known quasar SDSS J095253.83+011421.9 (J0952+0114) at z=3.02 has one of the most peculiar spectra discovered so far, showing the presence of narrow Ly$\alpha$ and broad metal emission lines. Although recent studies have suggested that a Proximate Damped Ly$\alpha$ system (PDLA) causes this peculiar spectrum, the origin of the gas associated with the PDLA is unknown. Here we report the results of MUSE observations that reveal a new giant ($\approx$ 100 physical kpc) Lyman $\alpha$ nebula. The detailed analysis of the Ly$\alpha$ velocity, velocity dispersion, and surface brightness profiles suggests that the J0952+0114 Ly$\alpha$ nebula shares similar properties of other QSO nebulae previously detected with MUSE, implying that the PDLA in J0952+0144 is covering only a small fraction of the QSO emission solid angle. We also detected bright and spectrally narrow CIV$\lambda$1550 and HeII$\lambda$1640 extended emission around J0952+0114 with velocity centroids similar to the peak of the extended and central narrow Ly$\alpha$ emission. The presence of a peculiarly bright, unresolved, and relatively broad HeII$\lambda$1640 emission in the central region at exactly the same PDLA redshift hints at the possibility that the PDLA originates in a clumpy outflow with a bulk velocity of about 500 km/s. The smaller velocity dispersion of the large scale Ly$\alpha$ emission suggests that the high-speed outflow is confined to the central region. Lastly, the derived spatially resolved HeII/Ly$\alpha$ and CIV/Ly$\alpha$ maps show a positive gradient with the distance to the QSO hinting at a non-homogeneous ionization parameter distribution.

Theory suggests that there are two primary modes of accretion through which dark matter halos acquire the gas to form and fuel galaxies, hot and cold mode accretion. In cold mode accretion, gas streams along cosmic web filaments to the center of the halo, allowing for the efficient delivery of star-forming fuel. Recently, two QSO-illuminated HI Lyman alpha (Ly{\alpha}) emitting objects were reported to have properties of cold, rotating structures (Martin et al. 2015, Martin et al. 2016). However, the spatial and spectral resolution available was insufficient to constrain radial flows associated with connecting filaments. With the Keck Cosmic Web Imager (KCWI) we now have eight times the spatial resolution, permitting the detection of these in-spiraling flows. In order to detect these inflows, we introduce a suite of models which incorporate zonal radial flows, demonstrate their performance on a numerical simulation that exhibits coldflow accretion, and show that they are an excellent match to KCWI velocity maps of two Ly{\alpha} emitters observed around high-redshift quasars. These Multi-Filament Inflow models kinematically isolate zones of radial inflow that correspond to extended filamentary emission. The derived gas flux and inflow path is sufficient to fuel the inferred central galaxy star formation rate and angular momentum. Thus, our kinematic emission maps provide strong evidence for the inflow of gas from the cosmic web building galaxies at the peak of star formation.

We present results from a survey of galaxies in the fields of six z>3 Damped Lyman alpha systems (DLAs) using the Multi Unit Spectroscopic Explorer (MUSE) at the Very Large Telescope (VLT). We report a high detection rate of up to ~80% of galaxies within 1000 km/s from DLAs and with impact parameters between 25 and 280 kpc. In particular, we discovered 5 high-confidence Lyman alpha emitters associated with three DLAs, plus up to 9 additional detections across five of the six fields. The majority of the detections are at relatively large impact parameters (>50 kpc) with two detections being plausible host galaxies. Among our detections, we report four galaxies associated with the most metal-poor DLA in our sample (Z/Z_sun = -2.33), which trace an overdense structure resembling a filament. By comparing our detections with predictions from the Evolution and Assembly of GaLaxies and their Environments (EAGLE) cosmological simulations and a semi-analytic model designed to reproduce the observed bias of DLAs at z>2, we conclude that our observations are consistent with a scenario in which a significant fraction of DLAs trace the neutral regions within halos with a characteristic mass of 10^11-10^12 M_sun, in agreement with the inference made from the large-scale clustering of DLAs. We finally show how larger surveys targeting ~25 absorbers have the potential of constraining the characteristic masses of halos hosting high-redshift DLAs with sufficient accuracy to discriminate between different models.

We investigate the baryon content of the circumgalactic medium (CGM) within the virial radius of $M_h \sim 10^{12} \; M_\odot$ haloes at z ~ 3, by modelling the surface brightness profile of the giant Ly$\alpha$ nebulae recently discovered by MUSE around bright QSOs at this redshift. We initially assume fluorescent emission from cold photo-ionized gas confined by the pressure of a hot halo. Acceptable CGM baryon fractions (equal or smaller than the cosmological value) require that the cold gas occupies $\lesssim$ 1% of the volume, but is about as massive as the hot gas. CGM baryon fractions as low as 30% of the cosmic value, as predicted by some strongly ejective feedback models at this redshift, are not easy to reconcile with observations, under our assumptions, unless both the QSO-hosting haloes at $z\sim3$ are more massive than recent BOSS estimates based on clustering and the photo-ionized gas is colder than expected in a standard QSO ionizing radiation field. We also consider the option that the emission is dominated by photons scattered from the QSO broad line region. In this scenario, a very stringent lower limit to the baryon fraction can be obtained under the extreme assumption of optically thin scattering. We infer in this case a baryon fraction of at least 70% of the cosmic value, for fiducial parameters. Lower values require halo masses or gas temperatures different than expected, or that some mechanism keeps the cold gas systematically over-pressured with respect to the ambient medium.

Cosmic filaments play a crucial role in galaxy evolution transporting gas from the intergalactic medium into galaxies. However, little is known about the efficiency of this process and whether the gas is accreted in a homogenous or clumpy way. Recent observations suggest the presence of broad gas density distributions in the circumgalactic medium which could be related to the accretion of filaments. By means of high-resolution hydrodynamical simulations, we explore here the evolution of cold accreting filaments flowing through the hot circumgalactic medium (CGM) of high-z galaxies. In particular, we examine the nonlinear effects of Kelvin-Helmholtz instability (KHI) on the development of broad gas density distributions and on the formation of cold, dense clumps. We explore a large parameter space in filament and perturbation properties, such as, filament Mach number, initial perturbation wavelength, and thickness of the interface between the filament and the halo. We find that the time averaged density distribution of the cold gas is qualitatively consistent with a skewed log-normal probability distribution function (PDF) plus an additional component in form of a high density tail for high Mach-numbers. Our results suggest a tight correlation between the accreting velocity and the maximum densities developing in the filament which is consistent with the variance-Mach number relation for turbulence. Therefore, cosmological accretion could be a viable mechanism to produce turbulence and broad gas density distributions within the CGM.

We present first results from Multi Unit Spectroscopic Explorer (MUSE) observations at the Very Large Telescope in the MUSE Ultra Deep Field (MUDF), a $\approx 1.2\times 1.4$ arcmin$^2$ region for which we are collecting $\approx$200 hours of integral field spectroscopy. The $\approx 40$-hour observation completed to date reveals the presence of a group of three Ly$\alpha$ nebulae associated with a bright quasar pair at $z\simeq3.23$ with projected separation of $\approx 500\rm~kpc$. Two of the nebulae are physically associated with the quasars which are likely powering the Ly$\alpha$ emission, and extend for $\gtrsim 100~\rm kpc$ at a surface brightness level of $\approx 6\times 10^{-19}~\rm erg~s^{-1}~cm^{-2}~arcsec^{-2}$. A third smaller ($\approx$35 kpc) nebula lies at a velocity offset of $\approx 1550$ km s$^{-1}$. Despite their clustered nature, the two large nebulae have properties similar to those observed in isolated quasars and exhibit no sharp decline in flux at the current depth, suggesting an even more extended distribution of gas around the quasars. We interpret the shape and the alignment of the two brighter nebulae as suggestive of the presence of an extended structure connecting the two quasar host galaxies, as seen for massive galaxies forming within gas-rich filaments in cosmological simulations.

We present Keck/KCWI integral field spectrograph observations of the complex system surrounding SMM J02399$-$0136 (a lensed $z=2.8$ sub-mm galaxy), including an associated Ly$\rm \alpha$ nebula, a dust-obscured, broad-absorption-line quasar, and neighboring galaxies. At a 3$\sigma$ surface brightness contour of 1.6$\times$ 10$^{-17}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$, the Ly$\rm \alpha$ nebula extends over 17 arcsec ($\gtrsim 140$ physical kpc) and has a total Ly$\rm \alpha$ luminosity of $2.5 \times 10^{44} \, \rm erg \, s^{-1}$ (uncorrected for lensing). The nebula exhibits a kinematic shear of $\sim$ 1000 km s$^{-1}$ over 100 pkpc with lowest velocities east of SMM J02399$-$0136 and increasing to the southwest. We also discover a bright, Ly$\rm \alpha$ emitter, separated spatially and kinematically from the nebula, at a projected separation of $\approx$60 kpc from the quasar. This source has no clear central counterpart in deep $Hubble \, Space \, Telescope$ imaging, giving an intrinsic Ly$\rm \alpha$ rest-frame equivalent width greater than 312 $ \rm \mathring{A}$ (5-$\sigma$). We argue that this 'dark cloud' is illuminated by the quasar with a UV flux that is orders of magnitude brighter than the emission along our sightline. This result confirms statistical inferences that luminous quasars at $z>2$ emit UV radiation anisotropically. Future KCWI observations of other lines, e.g. Ly$\rm \beta$, He$\rm _{II}$, C$\rm _{IV}$, etc, and with polarimetry will further reveal the origin of the Ly$\rm \alpha$ nebula and nature of the dark cloud.

Deep optical spectroscopic surveys of galaxies provide us a unique opportunity to investigate rest-frame ultra-violet (UV) emission line properties of galaxies at ${z \sim 2-4.5}$. Here we combine VLT/MUSE Guaranteed Time Observations of the Hubble Deep Field South, Ultra Deep Field, COSMOS, and several quasar fields with other publicly available data from VLT/VIMOS and VLT/FORS2 to construct a catalogue of HeII${\lambda}$1640 emitters at ${z\sim2}$. The deepest areas of our MUSE pointings reach a ${3\sigma}$ line flux limit of 3.1${\times}$10-19 erg/ s/ cm$^2$. After discarding broad line active galactic nuclei we find 13 HeII${\lambda}$1640 detections from MUSE with a median MUV = $-20.1$ and 21 tentative HeII${\lambda}$1640 detections from other public surveys. Excluding Ly${\alpha}$, all except two galaxies in our sample show at least one other rest-UV emission line, with CIII]${\lambda}$1907,${\lambda}$1909 being the most prominent. We use multi-wavelength data available in the Hubble legacy fields to derive basic galaxy properties of our sample via spectral energy distribution fitting techniques. Taking advantage of the high quality spectra obtained by MUSE (${\sim10 - 30}$h of exposure time per pointing), we use photo-ionisation models to study the rest-UV emission line diagnostics of the HeII${\lambda}$1640 emitters. Line ratios of our sample can be reproduced by moderately sub-solar photo-ionisation models, however, we find that including effects of binary stars lead to degeneracies in most free parameters. Even after considering extra ionising photons produced by extreme sub-solar metallicity binary stellar models, photo-ionisation models are unable to reproduce rest-frame HeII${\lambda}$1640 equivalent widths (${\sim}$ 0.2 - 10 A), thus additional mechanisms are necessary in models to match the observed HeII${\lambda}$1640 properties.

This paper presents a study of the galactic environment of a chemically-pristine (<0.6% solar metallicity) Lyman Limit system (LLS) discovered along the sightline toward QSO SDSSJ135726.27+043541.4 (zQSO=1.233) at projected distance d=126 physical kpc (pkpc) from a luminous red galaxy (LRG) at z=0.33. Combining deep Hubble Space Telescope images, MUSE integral field spectroscopic data, and wide-field redshift survey data has enabled an unprecedented, ultra-deep view of the environment around this LRG-LLS pair. A total of 12 galaxies, including the LRG, are found at d<~400 pkpc and line-of-sight velocity dv<600 km/s of the LLS, with intrinsic luminosity ranging from 0.001L* to 2L* and a corresponding stellar mass range of Mstar=10^{7-11} Msun. All 12 galaxies contribute to a total mass of Mstar=1.6e11 Msun with ~80% contained in the LRG. The line-of-sight velocity dispersion of these galaxies is found to be {\sigma}_group=230 km/s with the center of mass at d_group=118 pkpc and line-of-sight velocity offset of {\Delta}v_group=181 km/s from the LLS. Three of these are located at d<~100 pkpc from the LLS, and they are all faint with intrinsic luminosity <0.02 L* and gas phase metallicity of ~10% solar in their interstellar medium. The disparity in the chemical enrichment level between the LLS and the group members suggests that the LLS originates in infalling intergalactic medium and that parts of the intergalactic gas near old and massive galaxies can still remain chemically pristine through the not too distant past.

Enormous Ly$\alpha$ Nebulae (ELANe), unique tracers of galaxy density peaks, are predicted to lie at the nodes and intersections of cosmic filamentary structures. Previous successful searches for ELANe have focused on wide-field narrowband surveys, or have targeted known sources such as ultraluminous quasi-stellar-objects (QSOs) or radio galaxies. Utilizing groups of coherently strong Ly$\alpha$ absorptions (CoSLAs), we have developed a new method to identify high-redshift galaxy overdensities and have identified an extremely massive overdensity, BOSS1441, at $z=2-3$ (Cai et al. 2016a). In its density peak, we discover an ELAN that is associated with a relatively faint continuum. To date, this object has the highest diffuse Ly$\alpha$ nebular luminosity of $L_{\rm{nebula}}=5.1\pm0.1\times10^{44}$ erg s$^{-1}$. Above the 2$\sigma$ surface brightness limit of SB$_{\rm{Ly\alpha}}= 4.8\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$, this nebula has an end-to-end spatial extent of 442 kpc. This radio-quiet source also has extended \civ\ $\lambda1549$ and \heii\ $\lambda1640$ emission on $\gtrsim30$ kpc scales. Note that the Ly$\alpha$, \heii\ and \civ\ emission all have double-peaked line profiles. Each velocity component has a full-width-half-maximum (FWHM) of $\approx700 - 1000$ km s$^{-1}$. We argue that this Ly$\alpha$ nebula could be powered by shocks due to an AGN-driven outflow or/and photoionization by a strongly obscured source.

With a projected size of about 450 kpc at z~2.3, the Slug Ly-alpha nebula is a rare laboratory to study, in emission, the properties of the intergalactic gas in the Cosmic Web. Since its discovery, the Slug has been the subject of several spectroscopic follow-ups to constrain the properties of the emitting gas. Here we report the results of a deep MUSE integral-field spectroscopic search for non-resonant, extended HeII1640 and metal emission. Extended HeII radiation is detected on scales of about 100 kpc, but only in some regions associated with the bright Ly-alpha emission and a continuum-detected source, implying large and abrupt variations in the line ratios across adjacent regions in projected space. The recent detection of associated H-alpha emission and similar abrupt variations in the Ly-alpha kinematics, strongly suggest that the HeII/Ly-alpha gradient is due to large variations in the physical distances between the associated quasar and these regions. This implies that the overall length of the emitting structure could extend to physical Mpc scales and be mostly oriented along our line of sight. At the same time, the relatively low HeII/Ly-alpha values suggest that the emitting gas has a broad density distribution that - if expressed in terms of a lognormal - implies dispersions as high as those expected in the interstellar medium of galaxies. These results strengthen the possibility that the density distribution of intergalactic gas at high-redshift is extremely clumpy and multiphase on scales below our current observational spatial resolution of a few physical kpc.

We report the discovery of six spatially extended (10-100 kpc) line-emitting nebulae in the z=0.57 galaxy group hosting PKS0405-123, one of the most luminous quasars at z<1. The discovery is enabled by the Multi Unit Spectroscopic Explorer (MUSE) and provides tantalizing evidence connecting large-scale gas streams with nuclear activity on scales of <10 proper kpc (pkpc). One of the nebulae exhibits a narrow, filamentary morphology extending over 50 pkpc toward the quasar with narrow internal velocity dispersion (50 km/s) and is not associated with any detected galaxies, consistent with a cool intragroup medium (IGrM) filament. Two of the nebulae are 10 pkpc North and South of the quasar with tidal arm like morphologies. These two nebulae, along with a continuum emitting arm extending 60 pkpc from the quasar are signatures of interactions which are expected to redistribute angular momentum in the host interstellar medium (ISM) to facilitate star formation and quasar fueling in the nucleus. The three remaining nebulae are among the largest and most luminous [O III] emitting `blobs' known (1400-2400 pkpc^2) and correspond both kinematically and morphologically with interacting galaxy pairs in the quasar host group, consistent with arising from stripped ISM rather than large-scale quasar outflows. The presence of these large- and small-scale nebulae in the vicinity of a luminous quasar bears significantly on the effect of large-scale environment on galaxy and black hole fueling, providing a natural explanation for the previously known correlation between quasar luminosity and cool circumgalactic medium (CGM).

Galaxies are surrounded by large reservoirs of gas, mostly hydrogen, fed by inflows from the intergalactic medium and by outflows due to galactic winds. Absorption-line measurements along the sightlines to bright and rare background quasars indicate that this circumgalactic medium pervades far beyond the extent of starlight in galaxies, but very little is known about the spatial distribution of this gas. A new window into circumgalactic environments was recently opened with the discovery of ubiquitous extended Lyman-alpha emission from hydrogen around high-redshift galaxies, facilitated by the extraordinary sensitivity of the MUSE instrument at the ESO Very Large Telescope. Due to the faintness of this emission, such measurements were previously limited to especially favourable systems or to massive statistical averaging. Here we demonstrate that low surface brightness Lyman-alpha emission surrounding faint galaxies at redshifts between 3 and 6 adds up to a projected sky coverage of nearly 100%. The corresponding rate of incidence (the mean number of Lyman-alpha emitters penetrated by any arbitrary line of sight) is well above unity and similar to the incidence rate of high column density absorbers frequently detected in the spectra of distant quasars. This similarity suggests that most circumgalactic atomic hydrogen at these redshifts has now been detected also in emission.

Using an ultra-deep blind survey with the MUSE integral field spectrograph on the ESO Very Large Telescope, we obtain spectroscopic redshifts to a depth never explored before: galaxies with observed magnitudes m > 30 - 32. Specifically, we detect objects via Lyman-alpha emission at 2.9 < z < 6.7 without individual continuum counterparts in areas covered by the deepest optical/near-infrared imaging taken by the Hubble Space Telescope, the Hubble Ultra Deep Field. In total, we find 102 such objects in 9 square arcminutes at these redshifts. Detailed stacking analyses confirm the Lyman-alpha emission as well as the 1216 Angstrom-breaks and faint UV continua (M_UV ~ -15). This makes them the faintest spectroscopically-confirmed objects at these redshifts, similar to the sources believed to reionize the universe. A simple model for the expected fraction of detected/undetected Lyman-alpha emitters as a function of luminosity is consistent with these objects being the high-equivalent width tail of the normal Lyman-alpha-emitter population at these redshifts.

Motivated by the recent discovery of rare Enormous Lyman-Alpha Nebulae (ELAN) around z~2 quasars, we have initiated a long-term observational campaign with the MUSE instrument to directly uncover the astrophysics of the gas around quasars. We present here the first 61 targets of our effort under the acronym QSO MUSEUM (Quasar Snapshot Observations with MUse: Search for Extended Ultraviolet eMission). These quasars have a median redshfit of z=3.17, absolute $i$ magnitude in the range $-29.67\leq M_i(z=2)\leq-27.03$, and different levels of radio-loudness. This sample unveils diverse specimens of Ly$\alpha$ nebulosities extending for tens of kiloparsecs around these quasars (on average out to a maximum projected distance of 80 kpc) above a surface brightness SB$>8.8\times10^{-19}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ ($2\sigma$). The bulk of the extended Ly$\alpha$ emission is within R<50 kpc, and is characterized by relatively quiescent kinematics, with average velocity dispersions of $\langle \sigma_{\rm Ly\alpha}\rangle < 400$~km~s$^{-1}$. Therefore, the motions within all these Ly$\alpha$ nebulosities have amplitudes consistent with gravitational motions expected in dark matter halos hosting quasars at these redshifts, possibly reflecting the complexity in propagating a fast wind on large scales. Our current data suggest a combination of photoionization and resonant scattering as powering mechanisms of the Ly$\alpha$ emission. We discover the first $z\sim3$ ELAN, which confirms a very low probability ($\sim1\%$) of occurrence of such extreme systems at these cosmic epochs. Finally, we discuss the redshift evolution currently seen in extended Ly$\alpha$ emission around radio-quiet quasars from $z\sim3$ to $z\sim2$, concluding that it is possibly linked to a decrease of cool gas mass within the quasars' CGM from $z\sim3$ to $z\sim2$, and thus to the balance of cool vs hot media.

In a recent paper (Fumagalli et al. 2017) we reported on the detection of a diffuse H$\alpha$ glow in the outskirts of the nearby, edge-on disc galaxy UGC 7321 observed with the Multi Unit Spectroscopic Explorer (MUSE) at the ESO Very Large Telescope. By interpreting the H$\alpha$ emission as fluorescence arising from hydrogen ionised by an external (i.e., extragalactic) radiation field, we estimated the UV background (UVB) intensity in terms of HI ionisation rate (per ion) at $z\simeq 0$ to be in the range $\Gamma_{\rm HI}\sim 6-8\times 10^{-14}$ s$^{-1}$. In the present work, by performing radiative transfer calculations over a large set of models of the gaseous disc of UGC 7321, we refine our estimate and through an MCMC analysis derive a value for the photoionisaton rate of $\Gamma_{\rm HI}=7.27^{+2.93}_{-2.90}\times 10^{-14}$ s$^{-1}$. In particular, our analysis demonstrates that this value is robust against large variations in the galaxy model and that the uncertainties are mainly driven by the errors associated with the observed H$\alpha$ surface brightness. Our measurement is consistent with several recent determinations of the same quantity by a completely independent technique (i.e., flux decrement analysis of the Ly$\alpha$-forest), and support the notion that the low redshift UVB is largely dominated by active galactic nuclei (AGNs), possibly with no need of further contribution from star forming galaxies.

We present a study of the [OIII]5007/[OII]3727 (O32) ratios of star-forming galaxies drawn from MUSE data spanning a redshift range 0.28<z<0.85. Recently discovered Lyman continuum (LyC) emitters have extremely high oxygen line ratios: O32>4. Here we aim to understand the properties and the occurrences of galaxies with such high line ratios. Combining data from several MUSE GTO programmes, we select a population of star-forming galaxies with bright emission lines, from which we draw 406 galaxies for our analysis based on their position in the z-dependent star formation rate (SFR) - stellar mass (M*) plane. Out of this sample 15 are identified as extreme oxygen emitters based on their O32 ratios (3.7%) and 104 galaxies have O32>1 (26%). Our analysis shows no significant correlation between M*, SFR, and the distance from the SFR-M* relation with O32. We find a decrease in the fraction of galaxies with O32>1 with increasing M*, however, this is most likely a result of the relationship between O32 and metallicity, rather than between O32 and M*. We draw a comparison sample of local analogues with <z>~0.03 from SDSS, and find similar incidence rates for this sample. In order to investigate the evolution in the fraction of high O32 emitters with redshift, we bin the sample into three redshift subsamples of equal number, but find no evidence for a dependence on redshift. Furthermore, we compare the observed line ratios with those predicted by nebular models with no LyC escape and find that most of the extreme oxygen emitters can be reproduced by low metallicity models. The remaining galaxies are likely LyC emitter candidates. Finally, based on a comparison between electron temperature estimates from the [OIII4363]/[OIII]5007 ratio of the extreme oxygen emitters and nebular models, we argue that the galaxies with the most extreme O32 ratios have young light-weighted ages.

The Slug Nebula is one of the largest and most luminous Lyman-alpha (LyA) nebulae discovered to date, extending over 450 kiloparsecs (kpc) around the bright quasar UM287 at z=2.283. Characterized by high surface brightnesses over intergalactic scales, its LyA emission may either trace high-density ionized gas ("clumps") or large column densities of neutral material. To distinguish between these two possibilities, information from a non-resonant line such as Halpha is crucial. Therefore, we analyzed a deep MOSFIRE observation of one of the brightest LyA emitting regions in the Slug Nebula with the goal of detecting associated Halpha emission. We also obtained a deep, moderate resolution LyA spectrum of the nearby brightest region of the Slug. We detected an Halpha flux of F_(Halpha)= 2.62 +/- 0.47 x 10^-17 erg/cm^2/s (SB_(Halpha)=2.70 +/- 0.48 x 10^-18 erg/cm^2/s/sq") at the expected spatial and spectral location. Combining the Halpha detection with its corresponding LyA flux (determined from the narrow-band imaging) we calculate a flux ratio of F_(LyA_/F_(Halpha)= 5.5 +/- 1.1. The presence of a skyline at the location of the Halpha emission decreases the signal to noise ratio of the detection and our ability to put stringent constraints on the Halpha kinematics. Our measurements argue for the origin of the LyA emission being recombination radiation, suggesting the presence of high-density ionized gas. Finally, our high-resolution spectroscopic study of the LyA emission does not show evidence of a rotating disk pattern and suggest a more complex origin for at least some parts of the Slug Nebula.

Post starburst E+A galaxies show indications of a powerful starburst that was quenched abruptly. Their disturbed, bulge-dominated morphologies suggest that they are merger remnants. The more massive E+A galaxies are suggested to be quenched by AGN feedback, yet little is known about AGN-driven winds in this short-lived phase. We present spatially-resolved IFU spectroscopy by the Keck Cosmic Web Imager of SDSS J003443.68+251020.9, at z=0.118. The system consists of two galaxies, the larger of which is a post starburst E+A galaxy hosting an AGN. Our modelling suggests a 400 Myrs starburst, with a peak star formation rate of 120 Msun/yr. The observations reveal stationary and outflowing gas, photoionized by the central AGN. We detect gas outflows to a distance of 17 kpc from the central galaxy, far beyond the region of the stars (about 3 kpc), inside a conic structure with an opening angle of 70 degrees. We construct self-consistent photoionization and dynamical models for the different gas components and show that the gas outside the galaxy forms a continuous flow, with a mass outflow rate of about 24 Msun/yr. The gas mass in the flow, roughly $10^9$ Msun, is larger than the total gas mass within the galaxy, some of which is outflowing too. The continuity of the flow puts a lower limit of 60 Myrs on the duration of the AGN feedback. Such AGN are capable of removing, in a single episode, most of the gas from their host galaxies and expelling enriched material into the surrounding CGM.

The Lyman alpha (lya) line of Hydrogen is a prominent feature in the spectra of star-forming galaxies, usually redshifted by a few hundreds of km/s compared to the systemic redshift. This large offset hampers follow-up surveys, galaxy pair statistics and correlations with quasar absorption lines when only lya is available. We propose diagnostics that can be used to recover the systemic redshift directly from the properties of the lya line profile. We use spectroscopic observations of Lyman-Alpha Emitters (LAEs) for which a precise measurement of the systemic redshift is available. Our sample contains 13 sources detected between z~3 and z~6 as part of various Multi Unit Spectroscopic Explorer (MUSE) Guaranteed Time Observations (GTO). We also include a compilation of spectroscopic lya data from the literature spanning a wide redshift range (z~0-8). First, restricting our analysis to double-peaked lya spectra, we find a tight correlation between the velocity offset of the red peak with respect to the systemic redshift, Vpeak, and the separation of the peaks. Secondly, we find a correlation between Vpeak and the full width at half maximum of the lya line. Fitting formulas, to estimate systemic redshifts of galaxies with an accuracy of +-100 km/s when only the lya emission line is available, are given for the two methods.

Enormous Ly$\alpha$ nebulae (ELANe) represent the extrema of Ly$\alpha$ nebulosities. They have detected extents of $>200$ kpc in Ly$\alpha$ and Ly$\alpha$ luminosities $>10^{44}$ erg s$^{-1}$. The ELAN population is an ideal laboratory to study the interactions between galaxies and the intergalactic/circumgalactic medium (IGM/CGM) given their brightness and sizes. The current sample size of ELANe is still very small, and the few $z\approx2$ ELANe discovered to date are all associated with local overdensities of active galactic nuclei (AGNs). Inspired by these results, we have initiated a survey of ELANe associated with QSO pairs using the Palomar and Keck Cosmic Web Imagers (PCWI/KCWI). In this letter, we present our first result: the discovery of ELAN0101+0201 associated with a QSO pair at $z=2.45$. Our PCWI discovery data shows that, above a 2-$\sigma$ surface brightness of $1.2\times10^{-17}$ \sbunit, the end-to-end size of ELAN0101+0201 is $\gtrsim 232$ kpc. We have conducted follow-up observations using KCWI, resolving multiple Ly$\alpha$ emitting sources within the rectangular field-of-view of $\approx 130\times165$ projected kpc$^2$, and obtaining their emission line profiles at high signal-to-noise ratios. Combining both KCWI and PCWI, our observations confirm that ELAN0101+0201 resides in an extremely overdense environment. Our observations further support that a large amount of cool ($T\sim10^4$K) gas could exist in massive halos (M$\gtrsim10^{13}$M$_\odot$) at $z\approx2$. Future observations on a larger sample of similar systems will provide statistics of how cool gas is distributed in massive overdensities at high-redshift and strongly constrain the evolution of the intracluster medium (ICM).

In this work we report deep MUSE observations of a Broad Absorption Line (BAL) quasar at z ~ 5, revealing a Lyman alpha nebula with a maximum projected linear size of ~ 60 kpc around the quasar (down to our 2-sigma SB limit per layer of ~ 9e-19 erg/s/cm^2/arcsec^2 for a 1 arcsec^2 aperture). After correcting for the cosmological surface brightness dimming, we find that our nebula, at z ~ 5, has an intrinsically less extended Lyman alpha emission than nebulae at lower redshift. However, such a discrepancy is greatly reduced when referring to comoving distances, which take into account the cosmological growth of dark matter (DM) haloes, suggesting a positive correlation between the size of Lyman alpha nebulae and the sizes of DM haloes/structures around quasars. Differently from the typical nebulae around radio-quiet non-BAL quasars, in the inner regions (~ 10 kpc) of the circumgalactic medium (CGM) of our source, the velocity dispersion of the Lyman alpha emission is very high (FWHM > 1000 km/s), suggesting that in our case we may be probing outflowing material associated with the quasar.

We report the discovery of a 10^4 kpc^2 gaseous structure detected in [OII] in an over-dense region of the COSMOS-Gr30 galaxy group at z~0.725 thanks to deep MUSE Guaranteed Time Observations. We estimate the total amount of diffuse ionised gas to be of the order of (~5+-3)x10^10 Msun and explore its physical properties to understand its origin and the source(s) of the ionisation. The MUSE data allow the identification of a dozen of group members embedded in this structure from emission and absorption lines. We extracted spectra from small apertures defined for both the diffuse ionised gas and the galaxies. We investigated the kinematics and ionisation properties of the various galaxies and extended gas regions thanks to line diagnostics (R23, O32 and [OIII]/H\beta) available within the MUSE wavelength range. We compared these diagnostics to photo-ionisation models and shock models. The structure is divided in two kinematically distinct sub-structures. The most extended sub-structure of ionised gas is likely rotating around a massive galaxy and displays filamentary patterns linking some galaxies. The second sub-structure links another massive galaxy hosting an Active Galactic Nucleus to a low mass galaxy but also extends orthogonally to the AGN host disk over ~35 kpc. This extent is likely ionised by the AGN itself. The location of small diffuse regions in the R23 vs. O32 diagram is compatible with photo-ionisation. However, the location of three of these regions in this diagram (low O32, high R23) can also be explained by shocks, which is supported by their large velocity dispersions. One edge-on galaxy shares the same properties and may be a source of shocks. Whatever the hypothesis, the extended gas seems to be non primordial. We favour a scenario where the gas has been extracted from galaxies by tidal forces and AGN triggered by interactions between at least the two sub-structures.

We present rest-frame Ly$\alpha$ equivalent widths (EW) of 417 Ly$\alpha$ emitters (LAEs) detected with Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT) at $2.9 < z < 6.6$ in the Hubble Ultra Deep Field. Based on the deep MUSE spectroscopy and ancillary Hubble Space Telescope (HST) photometry data, we carefully measured EW values taking into account extended Ly$\alpha$ emission and UV continuum slopes ($\beta$). Our LAEs reach unprecedented depths, both in Ly$\alpha$ luminosities and UV absolute magnitudes, from log($L_{\rm Ly\alpha}$/erg s$^{-1}$) $\sim$41.0 to 43.0 and from Muv $\sim$ -16 to -21 (0.01-1.0 $L^{*}_{\rm z=3}$). The EW values span the range of $\sim$ 5 to 240 \AA\ or larger, and their distribution can be well fitted by an exponential law $N = N_{\rm 0}$ exp($-$EW/$w_{\rm 0}$). Owing to the high dynamic range in Muv, we find that the scale factor, $w_{\rm 0}$, depends on Muv in the sense that including fainter Muv objects increases $w_{\rm 0}$, i.e., the Ando effect. The results indicate that selection functions affect the EW scale factor. Taking these effects into account, we find that our $w_{\rm 0}$ values are consistent with those in the literature within $1\sigma$ uncertainties at $2.9 < z < 6.6$ at a given threshold of Muv and $L_{\rm Ly\alpha}$. Interestingly, we find 12 objects with EW $>200$ \AA\ above $1\sigma$ uncertainties. Two of these 12 LAEs show signatures of merger or AGN activity: the weak CIV $\lambda 1549$ emission line. For the remaining 10 very large EW LAEs, we find that the EW values can be reproduced by young stellar ages ($< 100$ Myr) and low metallicities ($\lesssim 0.02$ $Z_{\rm \odot}$). Otherwise, at least part of the Ly$\alpha$ emission in these LAEs needs to arise from anisotropic radiative transfer effects, fluorescence by hidden AGN or quasi-stellar object activity, or gravitational cooling.

We present a measurement of the fraction of Lyman $\alpha$ (Ly$\alpha$) emitters ($X_{\rm{Ly} \alpha}$) amongst HST continuum-selected galaxies at $3<z<6$ with the Multi-Unit Spectroscopic Explorer (MUSE) on the VLT. Making use of the first 24 MUSE-Wide pointings in GOODS-South, each having an integration time of 1 hour, we detect 100 Ly$\alpha$ emitters and find $X_{\rm{Ly} \alpha}\gtrsim0.5$ for most of the redshift range covered, with 29 per cent of the Ly$\alpha$ sample exhibiting rest equivalent widths (rest-EWs) $\leq$ 15\AA. Adopting a range of rest-EW cuts (0 - 75\AA), we find no evidence of a dependence of $X_{\rm{Ly} \alpha}$ on either redshift or UV luminosity.

We report the detection of extended Lyman-alpha (Lya) haloes around 145 individual star-forming galaxies at redshifts 3<z<6 in the Hubble Ultra Deep Field observed with the Multi-Unit Spectroscopic Explorer at ESO-VLT. Our sample consists of continuum-faint (-15> M_{UV}> -22) Lya emitters (LAEs). Using a 2D, two-component decomposition of Lya emission assuming circular exponential distributions, we measure scale lengths and luminosities of Lya haloes. We find that 80% of our objects having reliable Lya halo measurements show Lya emission that is significantly more extended than the UV continuum detected by HST (by a factor ~4 to >20). The median exponential scale length of the Lya haloes in our sample is ~4.5 kpc. By comparing the maximal detected extent of the Lya emission with the predicted dark matter halo virial radii of simulated galaxies, we show that the detected Lya emission of our selected sample of LAEs probes a significant portion of the cold circum-galactic medium (CGM) of these galaxies (>50% in average). This result shows that there must be significant HI reservoirs in the CGM and reinforces the idea that Lya haloes are ubiquitous around high-redshift Lya emitting galaxies. Our characterization of the Lya haloes indicates that the majority of the Lya flux comes from the halo (~65%) and that their scale lengths seem to be linked to the UV properties of the galaxies. We do not observe a significant Lya halo size evolution with redshift. We also find that the Lya lines cover a large range of full width at half maximum (FWHM) from 118 to 512 km/s. While the FWHM does not seem to be correlated to the Lya scale length, most compact Lya haloes and those that are not detected with high significance tend to have narrower Lya profiles. Finally, we investigate the origin of the extended Lya emission but we conclude that our data do not allow us to disentangle the possible processes.

We present the MUSE Hubble Ultra Deep Survey, a mosaic of nine MUSE fields covering 90\% of the entire HUDF region with a 10-hour deep exposure time, plus a deeper 31-hour exposure in a single 1.15 arcmin2 field. The improved observing strategy and advanced data reduction results in datacubes with sub-arcsecond spatial resolution (0.65 arcsec at 7000 A) and accurate astrometry (0.07 arcsec rms). We compare the broadband photometric properties of the datacubes to HST photometry, finding a good agreement in zeropoint up to mAB=28 but with an increasing scatter for faint objects. We have investigated the noise properties and developed an empirical way to account for the impact of the correlation introduced by the 3D drizzle interpolation. The achieved 3 sigma emission line detection limit for a point source is 1.5 and 3.1 10-19 erg.s-1.cm-2 for the single ultra-deep datacube and the mosaic, respectively. We extracted 6288 sources using an optimal extraction scheme that takes the published HST source locations as prior. In parallel, we performed a blind search of emission line galaxies using an original method based on advanced test statistics and filter matching. The blind search results in 1251 emission line galaxy candidates in the mosaic and 306 in the ultradeep datacube, including 72 sources without HST counterparts (mAB>31). In addition 88 sources missed in the HST catalog but with clear HST counterparts were identified. This data set is the deepest spectroscopic survey ever performed. In just over 100 hours of integration time, it provides nearly an order of magnitude more spectroscopic redshifts compared to the data that has been accumulated on the UDF over the past decade. The depth and high quality of these datacubes enables new and detailed studies of the physical properties of the galaxy population and their environments over a large redshift range.

In an effort to search for Ly$\alpha$ emission from circum- and intergalactic gas on scales of hundreds of kpc around $z\sim3$ quasars, and thus characterise the physical properties of the gas in emission, we have initiated an extensive fast-survey with the Multi Unit Spectroscopic Explorer (MUSE): Quasar Snapshot Observations with MUse: Search for Extended Ultraviolet eMission (QSO MUSEUM). In this work, we report the discovery of an enormous Ly$\alpha$ nebula (ELAN) around the quasar SDSS~J102009.99+104002.7 at $z=3.164$, which we followed-up with deeper MUSE observations. This ELAN spans $\sim297$ projected kpc, has an average Ly$\alpha$ surface brightness ${\rm SB}_{\rm Ly\alpha}\sim 6.04\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ (within the $2\sigma$ isophote), and is associated with an additional four, previously unknown embedded sources: two Ly$\alpha$ emitters and two faint active galactic nuclei (one Type-1 and one Type-2 quasar). By mapping at high significance the line-of-sight velocity in the entirety of the observed structure, we unveiled a large-scale coherent rotation-like pattern spanning $\sim300$ km s$^{-1}$ with a velocity dispersion of $<270$ km s$^{-1}$, which we interpret as a signature of the inspiraling accretion of substructures within the quasar's host halo. Future multiwavelength data will complement our MUSE observations, and are definitely needed to fully characterise such a complex system. None the less, our observations reveal the potential of new sensitive integral-field spectrographs to characterise the dynamical state of diffuse gas on large scales in the young Universe, and thereby witness the assembly of galaxies.

We utilize the Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT) to search for extended Lyman-Alpha emission around the z~6.6 QSO J0305-3150. After carefully subtracting the point-spread-function, we reach a nominal 5-sigma surface brightness limit of SB = 1.9x10$^{-18}$ erg/s/cm$^2$/arcsec$^2$ over a 1 arcsec$^2$ aperture, collapsing 5 wavelength slices centered at the expected location of the redshifted Lyman-Alpha emission (i.e. at 9256 Ang.). Current data suggest the presence (5-sigma, accounting for systematics) of a Lyman-Alpha nebula that extends for 9 kpc around the QSO. This emission is displaced and redshifted by 155 km/s with respect to the location of the QSO host galaxy traced by the [CII] emission line. The total luminosity is L = 3.0x10$^{42}$ erg/s. Our analysis suggests that this emission is unlikely to rise from optically thick clouds illuminated by the ionizing radiation of the QSO. It is more plausible that the Lyman-Alpha emission is due to fluorescence of the highly ionized optically thin gas. This scenario implies a high hydrogen volume density of n$_H$ ~ 6 cm$^{-3}$. In addition, we detect a Lyman-Alpha emitter (LAE) in the immediate vicinity of the QSO: i.e., with a projected separation of 12.5 kpc and a line-of-sight velocity difference of 560 km/s. The luminosity of the LAE is L = 2.1x10$^{42}$ erg/s and its inferred star-formation-rate is SFR ~ 1.3 M$_\odot$/yr. The probability of finding such a close LAE is one order of magnitude above the expectations based on the QSO-galaxy cross-correlation function. This discovery is in agreement with a scenario where dissipative interactions favour the rapid build-up of super-massive black holes at early Cosmic times.

Recent theoretical models suggest that the early phase of galaxy formation could involve an epoch when galaxies are gas-rich but inefficient at forming stars: a "dark galaxy" phase. Here, we report the results of our MUSE (Multi Unit Spectroscopic Explorer) survey for dark galaxies fluorescently illuminated by quasars at $z>3$. Compared to previous studies which are based on deep narrow-band (NB) imaging, our integral field survey provides a nearly uniform sensitivity coverage over a large volume in redshift space around the quasars as well as full spectral information at each location. Thanks to these unique features, we are able to build control samples at large redshift distances from the quasars using the same data taken under the same conditions. By comparing the rest-frame equivalent width (EW$_{0}$) distributions of the Ly$\alpha$ sources detected in proximity to the quasars and in control samples, we detect a clear correlation between the locations of high EW$_{0}$ objects and the quasars. This correlation is not seen in other properties such as Ly$\alpha$ luminosities or volume overdensities, suggesting the possible fluorescent nature of at least some of these objects. Among these, we find 6 sources without continuum counterparts and EW$_{0}$ limits larger than $240\,\mathrm{\AA}$ that are the best candidates for dark galaxies in our survey at $z>3.5$. The volume densities and properties, including inferred gas masses and star formation efficiencies, of these dark galaxy candidates are similar to previously detected candidates at $z\approx2.4$ in NB surveys. Moreover, if the most distant of these are fluorescently illuminated by the quasar, our results also provide a lower limit of $t=60$ Myr on the quasar lifetime.

We present new observations acquired with the Multi Unit Spectroscopic Explorer instrument on the Very Large Telescope in a quasar field that hosts a high column-density damped Ly{\alpha} absorber (DLA) at z~3.25. We detect Ly{\alpha} emission from a nebula at the redshift of the DLA with line luminosity (27+/-1)x1e41 erg/s, which extends over 37+/-1 kpc above a surface brightness limit of 6x1e-19 erg/s/cm2/arcsec2 at a projected distance of 30.5+/-0.5 kpc from the quasar sightline. Two clumps lie inside this nebula, both with Ly{\alpha} rest-frame equivalent width > 50 A and with relative line-of-sight velocities aligned with two main absorption components seen in the DLA spectrum. In addition, we identify a compact galaxy at a projected distance of 19.1+/-0.5 kpc from the quasar sightline. The galaxy spectrum is noisy but consistent with that of a star-forming galaxy at the DLA redshift. We argue that the Ly{\alpha} nebula is ionized by radiation from star formation inside the two clumps, or by radiation from the compact galaxy. In either case, these data imply the presence of a structure with size >>50 kpc inside which galaxies are assembling, a picture consistent with galaxy formation in groups and filaments as predicted by cosmological simulations such as the EAGLE simulations.

We present a clustering analysis of a sample of 238 Ly{$\alpha$}-emitters at redshift 3<z<6 from the MUSE-Wide survey. This survey mosaics extragalactic legacy fields with 1h MUSE pointings to detect statistically relevant samples of emission line galaxies. We analysed the first year observations from MUSE-Wide making use of the clustering signal in the line-of-sight direction. This method relies on comparing pair-counts at close redshifts for a fixed transverse distance and thus exploits the full potential of the redshift range covered by our sample. A clear clustering signal with a correlation length of r0 = 2.9(+1.0/-1.1) Mpc (comoving) is detected. Whilst this result is based on only about a quarter of the full survey size, it already shows the immense potential of MUSE for efficiently observing and studying the clustering of Ly{$\alpha$}-emitters.

Cosmological simulations suggest that most of the matter in the Universe is distributed along filaments connecting galaxies. Illuminated by the cosmic UV background (UVB), these structures are expected to glow in fluorescent Lyman alpha emission with a Surface Brightness (SB) that is well below current observational limits for individual detections. Here, we perform a stacking analysis of the deepest MUSE/VLT data using three-dimensional regions (subcubes) with orientations determined by the position of neighbouring Lyman alpha galaxies (LAEs) at 3<z<4. Our method should increase the probability of detecting filamentary Lyman alpha emission, provided that these structures are Lyman Limit Systems (LLSs). By stacking 390 oriented subcubes we reach a 2 sigma sensitivity level of SB ~ 0.44e-20 erg/s/cm^2/arcsec^2 in an aperture of 1 arcsec^2 x 6.25 Angstrom, which is three times below the expected fluorescent Lyman alpha signal from the Haardt-Madau 2012 (HM12) UVB at z~3.5. No detectable emission is found on intergalactic scales, implying that at least two thirds of our subcubes do not contain oriented LLSs for a HM12 UVB. On the other hand, significant emission is detected in the circum-galactic medium (CGM) of galaxies in the direction of the neighbours. The signal is stronger for galaxies with a larger number of neighbours and appears to be independent of any other galaxy properties such as luminosity, redshift and neighbour distance. We estimate that preferentially oriented satellite galaxies cannot contribute significantly to this signal, suggesting instead that gas densities in the CGM are typically larger in the direction of neighbouring galaxies on cosmological scales.

The radio-quiet quasar Q2059-360 at redshift $z=3.08$ is known to be close to a small Lyman $\alpha$ blob (LAB) and to be absorbed by a proximate damped Ly$\alpha$ (PDLA) system. Here, we present the Multi Unit Spectroscopic Explorer (MUSE) integral field spectroscopy follow-up of this quasi-stellar object (QSO). Our primary goal is to characterize this LAB in detail by mapping it both spatially and spectrally using the Ly$\alpha$ line, and by looking for high-ionization lines to constrain the emission mechanism. Combining the high sensitivity of the MUSE integral field spectrograph mounted on the Yepun telescope at ESO-VLT with the natural coronagraph provided by the PDLA, we map the LAB down to the QSO position, after robust subtraction of QSO light in the spectral domain. In addition to confirming earlier results for the small bright component of the LAB, we unveil a faint filamentary emission protruding to the south over about 80 pkpc (physical kpc); this results in a total size of about 120 pkpc. We derive the velocity field of the LAB (assuming no transfer effects) and map the Ly$\alpha$ line width. Upper limits are set to the flux of the N V $\lambda 1238-1242$, C IV $\lambda 1548-1551$, He II $\lambda 1640$, and C III] $\lambda 1548-1551$ lines. We have discovered two probable Ly$\alpha$ emitters at the same redshift as the LAB and at projected distances of 265 kpc and 207 kpc from the QSO; their Ly$\alpha$ luminosities might well be enhanced by the QSO radiation. We also find an emission line galaxy at $z=0.33$ near the line of sight to the QSO. This LAB shares the same general characteristics as the 17 others surrounding radio-quiet QSOs presented previously. However, there are indications that it may be centered on the PDLA galaxy rather than on the QSO.

We present new MUSE observations of a galaxy group probed by a background quasar. The quasar sightline passes between multiple $z=0.28$ galaxies, whilst showing at the same redshift low ionised metal line species, including Ca II, Mg I, Mg II and Fe II. Based on the galaxy redshifts measured from the MUSE data, we estimate the galaxies to be part of a small galaxy group with a halo mass of $\approx6\times10^{12}$ M$_{\odot}$. We use the MUSE data to reveal the two dimensional dynamical properties of the gas and stars in the group galaxies, and relate these to the absorber kinematics. With these data we consider a number of scenarios for the nature of the gas probed by the sightline absorbers: a co-rotating gas halo associated with a single galaxy within the group; outflowing material from a single group member powered by recent star-formation; and cool dense gas associated with an intra-group medium. We find that the dynamics, galaxy impact parameters, star-formation rates, and the absorber strength suggest the cool gas can not be clearly associated with any single galaxy within the group. Instead we find that the observations are consistent with a superposition of cool gas clouds originating with the observed galaxies as they fall into the group potential, and are now likely in the process of forming the intra-group medium.

We report the detection of extended Halpha emission from the tip of the HI disk of the nearby edge-on galaxy UGC 7321, observed with the Multi Unit Spectroscopic Explorer (MUSE) instrument at the Very Large Telescope. The Halpha surface brightness fades rapidly where the HI column density drops below N(HI) = 10^19 cm^-2 , consistent with fluorescence arising at the ionisation front from gas that is photoionized by the extragalactic ultraviolet background (UVB). The surface brightness measured at this location is (1.2 +/- 0.5)x10^-19 erg/s/cm^2/arcsec^2, where the error is mostly systematic and results from the proximity of the signal to the edge of the MUSE field of view, and from the presence of a sky line next to the redshifted Halpha wavelength. By combining the Halpha and the HI 21 cm maps with a radiative transfer calculation of an exponential disk illuminated by the UVB, we derive a value for the HI photoionization rate of Gamma ~ (6-8)x10^-14 1/s . This value is consistent with transmission statistics of the Lyalpha forest and with recent models of a UVB which is dominated by quasars.

Emission signatures from galactic winds provide an opportunity to directly map the outflowing gas, but this is traditionally challenging because of the low surface brightness. Using deep observations (27 hours) of the Hubble Deep Field South from the Multi Unit Spectroscopic Explorer (MUSE) instrument, we identify signatures of an outflow in both emission and absorption from a spatially resolved galaxy at z = 1.29 with a stellar mass M* = 8 x 10^9 Msun, star formation rate SFR = 77 Msun/yr, and star formation rate surface brightness 1.6 Msun/kpc^2 within the [OII] half-light radius R_1/2,[OII] = 2.76 +- 0.17 kpc. From a component of the strong resonant MgII and FeII absorptions at -350 km/s, we infer a mass outflow rate that is comparable to the star formation rate. We detect non-resonant FeII* emission, at lambda 2626, 2612, 2396, and 2365, at 1.2-2.4-1.5-2.7 x 10^-18 egs s-1 cm-2 respectively. These flux ratios are consistent with the expectations for optically thick gas. By combining the four non-resonant FeII* emission lines, we spatially map the FeII* emission from an individual galaxy for the first time. The FeII* emission has an elliptical morphology that is roughly aligned with the galaxy minor kinematic axis, and its integrated half-light radius R_1/2,FeII* = 4.1 +- 0.4 kpc is 50% larger than the stellar continuum (R_1/2,* = 2.34 +- 0.17 kpc) or the [OII] nebular line. Moreover, the FeII* emission shows a blue wing extending up to -400 km/s, which is more pronounced along the galaxy minor kinematic axis and reveals a C-shaped pattern in a p-v diagram along that axis. These features are consistent with a bi-conical outflow.

Several decades of observations and discoveries have shown that high-redshift AGN and massive galaxies are often surrounded by giant Lyman-alpha nebulae extending in some cases up to 500 kpc in size. In this review, I discuss the properties of the such nebulae discovered at z>2 and their connection with gas flows in and around the galaxies and their halos. In particular, I show how current observations are used to constrain the physical properties and origin of the emitting gas in terms of the Lyman-alpha photon production processes and kinematical signatures. These studies suggest that recombination radiation is the most viable scenario to explain the observed Lyman-alpha luminosities and Surface Brightness for the large majority of the nebulae and imply that a significant amount of dense, ionized and cold clumps should be present within and around the halos of massive galaxies. Spectroscopic studies suggest that, among the giant Lyman-alpha nebulae, the one associated with radio-loud AGN should have kinematics dominated by strong, ionized outflows within at least the inner 30-50 kpc. Radio-quiet nebulae instead present more quiescent kinematics compatible with stationary situation and, in some cases, suggestive of rotating structures. However, definitive evidences for accretion onto galaxies of the gas associated with the giant Lyman-alpha emission are not unambiguously detected yet. Deep surveys currently ongoing using other bright, non-resonant lines such as Hydrogen H-alpha and HeII1640 will be crucial to search for clearer signatures of cosmological gas accretion onto galaxies and AGN.

The total contribution of diffuse halo gas to the galaxy baryon budget strongly depends on its dominant ionization state. In this paper, we address the physical conditions in the highly ionized circumgalactic medium (CGM) traced by OVI absorption lines observed in COS-Halos spectra. We analyze the observed ionic column densities, absorption-line widths and relative velocities, along with the ratios of NV/OVI for 39 fitted Voigt profile components of OVI. We compare these quantities with the predictions given by a wide range of ionization models. Photoionization models that include only extragalactic UV background radiation are ruled out; conservatively, the upper limits to NV/OVI and measurements of N$_{\rm OVI}$ imply unphysically large path lengths > 100 kpc. Furthermore, very broad OVI absorption (b > 40 km s$^{-1}$) is a defining characteristic of the CGM of star-forming L* galaxies. We highlight two possible origins for the bulk of the observed OVI: (1) highly structured gas clouds photoionized primarily by local high energy sources or (2) gas radiatively cooling on large scales behind a supersonic wind. Approximately 20% of circumgalactic OVI does not align with any low-ionization state gas within $\pm$50 km s$^{-1}$ and is found only in halos with M$_{\rm halo}$ < 10$^{12}$ M$_{\odot}$. We suggest that this type of unmatched OVI absorption traces the hot corona itself at a characteristic temperature of 10$^{5.5}$ K. We discuss the implications of these very distinct physical origins for the dynamical state, gas cooling rates, and total baryonic content of L* gaseous halos.

The physical properties of galactic winds are one of the keys to understand galaxy formation and evolution. These properties can be constrained thanks to background quasar lines of sight (LOS) passing near star-forming galaxies (SFGs). We present the first results of the MusE GAs FLOw and Wind (MEGAFLOW) survey obtained of 2 quasar fields which have 8 MgII absorbers of which 3 have rest-equivalent width greater than 0.8 \AA. With the new Multi Unit Spectroscopic Explorer (MUSE) spectrograph on the Very Large Telescope (VLT), we detect 6 (75$\%$) MgII host galaxy candidates withing a radius of 30 arcsec from the quasar LOS. Out of these 6 galaxy--quasar pairs, from geometrical arguments, one is likely probing galactic outflows, two are classified as "ambiguous", two are likely probing extended gaseous disks and one pair seems to be a merger. We focus on the wind$-$pair and constrain the outflow using a high resolution quasar spectra from Ultraviolet and Visual Echelle Spectrograph (UVES). Assuming the metal absorption to be due to gas flowing out of the detected galaxy through a cone along the minor axis, we find outflow velocities of the order of $\approx$ 150 km/s (i.e. smaller than the escape velocity) with a loading factor, $\eta =\dot M_{\rm out}/$SFR, of $\approx$ 0.7. We see evidence for an open conical flow, with a low-density inner core. In the future, MUSE will provide us with about 80 multiple galaxy$-$quasar pairs in two dozen fields.

We present the first estimate of the Ly{\alpha} luminosity function using blind spectroscopy from the Multi Unit Spectroscopic Explorer, MUSE, in the Hubble Deep Field South. Using automatic source-detection software, we assemble a homogeneously-detected sample of 59 Ly{\alpha} emitters covering a flux range of -18.0 < log10 (F) < -16.3 (erg s^-1 cm^-2), corresponding to luminosities of 41.4 < log10 (L) < 42.8 (erg s^-1). As recent studies have shown, Ly{\alpha} fluxes can be underestimated by a factor of two or more via traditional methods, and so we undertake a careful assessment of each object's Ly{\alpha} flux using a curve-of-growth analysis to account for extended emission. We describe our self-consistent method for determining the completeness of the sample, and present an estimate of the global Ly{\alpha} luminosity function between redshifts 2.91 < z < 6.64 using the 1/Vmax estimator. We find the luminosity function is higher than many number densities reported in the literature by a factor of 2 - 3, although our result is consistent at the 1{\sigma} level with most of these studies. Our observed luminosity function is also in good agreement with predictions from semi-analytic models, and shows no evidence for strong evolution between the high- and low-redshift halves of the data. We demonstrate that one's approach to Ly{\alpha} flux estimation does alter the observed luminosity function, and caution that accurate flux assessments will be crucial in measurements of the faint end slope. This is a pilot study for the Ly{\alpha} luminosity function in the MUSE deep-fields, to be built on with data from the Hubble Ultra Deep Field which will increase the size of our sample by almost a factor of 10.

QSO near-zones are an important probe of the the ionization state of the IGM at z ~ 6-7, at the end of reionization. We present here high-resolution cosmological 3D radiative transfer simulations of QSO environments for a wide range of host halo masses, 10$^{10-12.5}$ M_sun. Our simulated near-zones reproduce both the overall decrease of observed near-zone sizes at 6 < z < 7 and their scatter. The observable near-zone properties in our simulations depend only very weakly on the mass of the host halo. The size of the H II region expanding into the IGM is generally limited by (super-)Lyman Limit systems loosely associated with (low-mass) dark matter haloes. This leads to a strong dependence of near-zone size on direction and drives the large observed scatter. In the simulation centred on our most massive host halo, many sightlines show strong red damping wings even for initial volume averaged neutral hydrogen fractions as low as ~ 10$^{-3}$. For QSO lifetimes long enough to allow growth of the central supermassive black hole while optically bright, we can reproduce the observed near-zone of ULAS J1120+0641 only with an IGM that is initially neutral. Our results suggest that larger samples of z > 7 QSOs will provide important constraints on the evolution of the neutral hydrogen fraction and thus on how late reionization ends.

We report on the search for galaxies in the proximity of two very metal-poor gas clouds at z~3 towards the quasar Q0956+122. With a 5-hour MUSE integration in a ~500x500 kpc^2 region centred at the quasar position, we achieve a >80% complete spectroscopic survey of continuum-detected galaxies with m<25 mag and Ly{\alpha} emitters with luminosity L>3e41 erg/s. We do not identify galaxies at the redshift of a z~3.2 Lyman limit system (LLS) with log Z/Zsun = -3.35 +/- 0.05, placing this gas cloud in the intergalactic medium or circumgalactic medium of a galaxy below our sensitivity limits. Conversely, we detect five Ly{\alpha} emitters at the redshift of a pristine z~3.1 LLS with log Z/Zsun < -3.8, while ~0.4 sources were expected given the z~3 Ly{\alpha} luminosity function. Both this high detection rate and the fact that at least three emitters appear aligned in projection with the LLS suggest that this pristine cloud is tracing a gas filament that is feeding one or multiple galaxies. Our observations uncover two different environments for metal-poor LLSs, implying a complex link between these absorbers and galaxy halos, which ongoing MUSE surveys will soon explore in detail. Moreover, in agreement with recent MUSE observations, we detected a ~90 kpc Ly{\alpha} nebula at the quasar redshift and three Ly{\alpha} emitters reminiscent of a "dark galaxy" population.

Direct Ly $\alpha$ imaging of intergalactic gas at $z\sim2$ has recently revealed giant cosmological structures around quasars, e.g. the Slug Nebula (Cantalupo et al. 2014). Despite their high luminosity, the detection rate of such systems in narrow-band and spectroscopic surveys is less than 10%, possibly encoding crucial information on the distribution of gas around quasars and the quasar emission properties. In this study, we use the MUSE integral-field instrument to perform a blind survey for giant Ly $\alpha$ nebulae around 17 bright radio-quiet quasars at $3<z<4$ that does not suffer from most of the limitations of previous surveys. After data reduction and analysis performed with specifically developed tools, we found that each quasar is surrounded by giant Ly $\alpha$ nebulae with projected sizes larger than 100 physical kpc and, in some cases, extending up to 320 kpc. The circularly averaged surface brightness profiles of the nebulae appear very similar to each other despite their different morphologies and are consistent with power laws with slopes $\approx-1.8$. The similarity between the properties of all these nebulae and the Slug Nebula suggests a similar origin for all systems and that a large fraction of gas around bright quasars could be in a relatively "cold" (T$\sim$10$^4$K) and dense phase. In addition, our results imply that such gas is ubiquitous within at least 50 kpc from bright quasars at $3<z<4$ independently of the quasar emission opening angle, or extending up to 200 kpc for quasar isotropic emission.

A toy model is developed to understand how the spatial distribution of fluorescent emitters in the vicinity of bright quasars could be affected by the geometry of the quasar bi-conical radiation field and by its lifetime. The model is then applied to the distribution of high equivalent width Lyman $\alpha$ emitters (with rest-frame equivalent widths above 100 A, threshold used in e.g. Trainor & Steidel, 2013) identified in a deep narrow-band 36x36 arcmin$^2$ image centered on the luminous quasar Q0420-388. These emitters are found to the edge of the field and show some evidence of an azimuthal asymmetry on the sky of the type expected if the quasar is radiating in a bipolar cone. If these sources are being fluorescently illuminated by the quasar, the two most distant objects require a lifetime of at least 15 Myr for an opening angle of 60 degrees or more, increasing to more than 40 Myr if the opening angle is reduced to a minimum 30 degrees. However, few of the other expected signatures of boosted fluorescence are seen at the current survey limits, e.g. a fall off in Lyman $\alpha$ brightness, or equivalent width, with distance. Furthermore, to have most of the Lyman $\alpha$ emission of the two distant sources to be fluorescently boosted would require the quasar to have been significantly brighter in the past. This suggests that these particular sources may not be fluorescent, invalidating the above lifetime constraints. This would cast doubt on the use of this relatively low equivalent width threshold and thus on the lifetime analysis in Trainor & Steidel (2013).

In the context of the FLASHLIGHT survey, we obtained deep narrow band images of 15 $z\sim2$ quasars with GMOS on Gemini-South in an effort to measure Ly$\alpha$ emission from circum- and inter-galactic gas on scales of hundreds of kpc from the central quasar. We do not detect bright giant Ly$\alpha$ nebulae (SB~10$^{-17}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ at distances >50 kpc) around any of our sources, although we routinely ($\simeq47$%) detect smaller scale <50 kpc Ly$\alpha$ emission at this SB level emerging from either the extended narrow emission line regions powered by the quasars or by star-formation in their host galaxies. We stack our 15 deep images to study the average extended Ly$\alpha$ surface brightness profile around $z\sim2$ quasars, carefully PSF-subtracting the unresolved emission component and paying close attention to sources of systematic error. Our analysis, which achieves an unprecedented depth, reveals a surface brightness of SB$_{\rm Ly\alpha}\sim10^{-19}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ at $\sim200$ kpc, with a $2.3\sigma$ detection of Ly$\alpha$ emission at SB$_{\rm Ly\alpha}=(5.5\pm3.1)\times10^{-20}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ within an annulus spanning 50 kpc <R< 500 kpc from the quasars. Assuming this Ly$\alpha$ emission is powered by fluorescence from highly ionized gas illuminated by the bright central quasar, we deduce an average volume density of $n_{\rm H}=0.6\times10^{-2}$ cm$^{-3}$ on these large scales. Our results are in broad agreement with the densities suggested by cosmological hydrodynamical simulations of massive ($M\simeq10^{12.5}M_\odot$) quasar hosts, however they indicate that the typical quasars at these redshifts are surrounded by gas that is a factor of ~100 times less dense than the (~1 cm$^{-3}$) gas responsible for the giant bright Ly$\alpha$ nebulae around quasars recently discovered by our group.

We use a background quasar to detect the presence of circum-galactic gas around a $z=0.91$ low-mass star forming galaxy. Data from the new Multi Unit Spectroscopic Explorer (MUSE) on the VLT show that the host galaxy has a dust-corrected star-formation rate (SFR) of 4.7$\pm$0.2 Msun/yr, with no companion down to 0.22 Msun/yr (5 $\sigma$) within 240 kpc (30"). Using a high-resolution spectrum (UVES) of the background quasar, which is fortuitously aligned with the galaxy major axis (with an azimuth angle $\alpha$ of only $15^\circ$), we find, in the gas kinematics traced by low-ionization lines, distinct signatures consistent with those expected for a "cold flow disk" extending at least 12 kpc ($3\times R_{1/2}$). We estimate the mass accretion rate $\dot M_{\rm in}$ to be at least two to three times larger than the SFR, using the geometric constraints from the IFU data and the HI column density of $\log N_{\rm HI} \simeq 20.4$ obtained from a {\it HST}/COS NUV spectrum. From a detailed analysis of the low-ionization lines (e.g. ZnII, CrII, TiII, MnII, SiII), the accreting material appears to be enriched to about 0.4 $Z_\odot$ (albeit with large uncertainties: $\log Z/Z_\odot=-0.4~\pm~0.4$), which is comparable to the galaxy metallicity ($12+\log \rm O/H=8.7\pm0.2$), implying a large recycling fraction from past outflows. Blue-shifted MgII and FeII absorptions in the galaxy spectrum from the MUSE data reveal the presence of an outflow. The MgII and FeII doublet ratios indicate emission infilling due to scattering processes, but the MUSE data do not show any signs of fluorescent FeII* emission.

We present the results obtained with MUSE on the core of the lensing cluster A1689. Integral-field observations with MUSE provide a unique view of the central region, allowing us to conduct a complete census on both cluster galaxies and lensed background sources, identified based on their spectral features without preselection. We investigate the multiple-image configuration for all known sources in the field. Previous to our survey, 28 different lensed galaxies displaying 46 multiple images were known in the MUSE field of view, most of them based on photometric redshifts and lensing considerations. Among them, we spectroscopically confirm 12 images based on their emission-lines, corresponding to 7 different lensed galaxies between z = 0.95 and 5.0. In addition, 14 new galaxies have been spectroscopically identified in this area, with redshifts ranging between 0.8 and 6.2. All background sources within the MUSE field of view correspond to multiple-imaged systems lensed by A1689. 17 sources in total are found at z > 3 based on their Lyman-alpha emission, with Lyman-alpha luminosities ranging between 40.5 < log(Ly{\alpha}) < 42.5 after correction for magnification. This sample is particularly sensitive to the slope of the LF toward the faintest-end. The density of sources obtained in this survey is consistent with a steep value of {\alpha} < -1.5, although this result still needs further investigation. These results illustrate the efficiency of MUSE in the characterization of lensing clusters on one hand, and the study of faint and distant populations of galaxies on the other hand. In particular, our current survey of lensing clusters should provide a unique census of sources responsible for the reionization in a representative volume at z ~ 4-7.

We report the detection of extended Ly alpha emission around individual star-forming galaxies at redshifts z = 3-6 in an ultradeep exposure of the Hubble Deep Field South obtained with MUSE on the ESO-VLT. The data reach a limiting surface brightness (1sigma) of ~1 x 10^-19 erg s^-1 cm^-2 arcsec^-2 in azimuthally averaged radial profiles, an order of magnitude improvement over previous narrowband imaging. Our sample consists of 26 spectroscopically confirmed Ly alpha-emitting, but mostly continuum-faint (m_AB >~ 27) galaxies. In most objects the Ly alpha emission is considerably more extended than the UV continuum light. While 5 of the faintest galaxies in the sample show no significantly detected Ly alpha haloes, the derived upper limits suggest that this is just due to insufficient S/N. Ly alpha haloes therefore appear to be (nearly) ubiquitous even for low-mass (~10^8-10^9 M_sun) star-forming galaxies at z>3. We decompose the Ly alpha emission of each object into a compact `continuum-like' and an extended halo component, and infer sizes and luminosities of the haloes. The extended Ly alpha emission approximately follows an exponential surface brightness distribution with a scale length of a few kpc. While these haloes are thus quite modest in terms of their absolute sizes, they are larger by a factor of 5-15 than the corresponding rest-frame UV continuum sources as seen by HST. They are also much more extended, by a factor ~5, than Ly alpha haloes around low-redshift star-forming galaxies. Between ~40% and >90% of the observed Ly alpha flux comes from the extended halo component, with no obvious correlation of this fraction with either the absolute or the relative size of the Ly alpha halo. Our observations provide direct insights into the spatial distribution of at least partly neutral gas residing in the circumgalactic medium of low to intermediate mass galaxies at z > 3.

Spatially resolved studies of high redshift galaxies, an essential insight into galaxy formation processes, have been mostly limited to stacking or unusually bright objects. We present here the study of a typical (L$^{*}$, M$_\star$ = 6 $\times 10^9$ $M_\odot$) young lensed galaxy at $z=3.5$, observed with MUSE, for which we obtain 2D resolved spatial information of Ly$\alpha$ and, for the first time, of CIII] emission. The exceptional signal-to-noise of the data reveals UV emission and absorption lines rarely seen at these redshifts, allowing us to derive important physical properties (T$_e\sim$15600 K, n$_e\sim$300 cm$^{-3}$, covering fraction f$_c\sim0.4$) using multiple diagnostics. Inferred stellar and gas-phase metallicities point towards a low metallicity object (Z$_{\mathrm{stellar}}$ = $\sim$ 0.07 Z$_\odot$ and Z$_{\mathrm{ISM}}$ $<$ 0.16 Z$_\odot$). The Ly$\alpha$ emission extends over $\sim$10 kpc across the galaxy and presents a very uniform spectral profile, showing only a small velocity shift which is unrelated to the intrinsic kinematics of the nebular emission. The Ly$\alpha$ extension is $\sim$4 times larger than the continuum emission, and makes this object comparable to low-mass LAEs at low redshift, and more compact than the Lyman-break galaxies and Ly$\alpha$ emitters usually studied at high redshift. We model the Ly$\alpha$ line and surface brightness profile using a radiative transfer code in an expanding gas shell, finding that this model provides a good description of both observables.

All galaxies once passed through a hyperluminous quasar phase powered by accretion onto a supermassive black hole. But because these episodes are brief, quasars are rare objects typically separated by cosmological distances. In a survey for Lyman-alpha emission at redshift z ~ 2, we discovered a physical association of four quasars embedded in a giant nebula. Located within a substantial overdensity of galaxies, this system is probably the progenitor of a massive galaxy cluster. The chance probability of finding a quadruple quasar is estimated to be ~10^-7, implying a physical connection between Lyman-alpha nebulae and the locations of rare protoclusters. Our findings imply that the most massive structures in the distant universe have a tremendous supply (~ 10^11 solar masses) of cool dense (volume density ~1 cm^-3) gas, which is in conflict with current cosmological simulations.

The recent discovery by Cantalupo et al. (2014) of the largest (~500 kpc) and luminous Ly-alpha nebula associated with the quasar UM287 (z=2.279) poses a great challenge to our current understanding of the astrophysics of the halos hosting massive z~2 galaxies. Either an enormous reservoir of cool gas is required $M\simeq10^{12}$ $M_{\odot}$, exceeding the expected baryonic mass available, or one must invoke extreme gas clumping factors not present in high-resolution cosmological simulations. However, observations of Ly-alpha emission alone cannot distinguish between these two scenarios. We have obtained the deepest ever spectroscopic integrations in the HeII and CIV lines with the goal of detecting extended line emission, but detect neither line to a 3$\sigma$ limiting SB $\simeq10^{-18}$ erg/s/cm$^2$/arcsec$^2$. We construct models of the expected emission spectrum in the highly probable scenario that the nebula is powered by photoionization from the central hyper-luminous quasar. The non-detection of HeII implies that the nebular emission arises from a mass $M_{\rm c}\lesssim6.4\times10^{10}$ $M_{\odot}$ of cool gas on ~200 kpc scales, distributed in a population of remarkably dense ($n_{\rm H}\gtrsim3$ cm$^{-3}$) and compact ($R\lesssim20$ pc) clouds, which would clearly be unresolved by current cosmological simulations. Given the large gas motions suggested by the Ly-alpha line ($v\simeq$ 500 km/s), it is unclear how these clouds survive without being disrupted by hydrodynamic instabilities. Our study serves as a benchmark for future deep integrations with current and planned wide-field IFU such as MUSE, KCWI, and KMOS. Our work suggest that a $\simeq$ 10 hr exposure would likely detect ~10 rest-frame UV/optical emission lines, opening up the possibility of conducting detailed photoionization modeling to infer the physical state of gas in the CGM.

We observed the Hubble Deep Field South with the new panoramic integral field spectrograph MUSE that we built and just commissioned at the VLT. The data cube resulting from 27 hours of integration covers one arcmin^2 field of view at an unprecedented depth with a 1 sigma emission line surface brightness limit of 1x$10^{-19}$ erg/s/cm$^2$/arcsec$^2$ and contains ~90,000 spectra. We present the combined and calibrated data cube, and we perform a first-pass analysis of the sources detected in the HDF-S imaging. We measured the redshifts of 189 sources up to a magnitude F814W = 29.5, increasing by more than an order of magnitude the number of known spectroscopic redshifts in this field. We also discovered 26 Lya emitting galaxies which are not detected in the HST WFPC2 deep broad band images. The intermediate spectral resolution of 2.3{\AA} allows us to separate resolved asymmetric Lya emitters, [O II] emitters, and C III] emitters and the large instantaneous wavelength range of 4500{\AA} helps to identify single emission lines. We also show how the three dimensional information of MUSE helps to resolve sources which are confused at ground-based image quality. Overall, secure identifications are provided for 83% of the 227 emission line sources detected in the MUSE data cube and for 32% of the 586 sources identified in the HST catalog of Casertano et al 2000. The overall redshift distribution is fairly flat to z=6.3, with a reduction between z=1.5 to 2.9, in the well-known redshift desert. The field of view of MUSE also allowed us to detect 17 groups within the field. We checked that the number counts of [O II] and Ly-a emitters are roughly consistent with predictions from the literature. Using two examples we demonstrate that MUSE is able to provide exquisite spatially resolved spectroscopic information on intermediate redshift galaxies present in the field.

We use hydrodynamic simulations post-processed with the radiative-transfer code RADAMESH to assess recent claims that the low HeII opacity observed in z>3 quasar spectra may be incompatible with models of HeII reionization driven by the observed population of active galactic nuclei (AGNs). In particular, building upon our previous work, we consider an early population of sources and start the radiative-transfer calculation at redshifts z>=5. Our model faithfully reproduces the emissivity of optically selected AGNs as inferred from measurements of their luminosity function. We find that HeII reionization is very extended in redshift ({\Delta} z>=2) and highly spatially inhomogeneous. In fact, mock spectra extracted from the simulations show a large variability in the evolution of the HeII effective optical depth within chunks of size {\Delta} z=0.04. Regions with low opacity ({\tau}^{eff}_{HeII}<3) can be found at high redshift, in agreement with the most recent observations of UV-transmitting quasars. At the highest redshift currently probed by observations (z~3.4), our updated model predicts a much lower HeII effective optical depth than previous simulations in the literature relieving most of the tension with the current data, that, however, still persists at about the (Gaussian) 1{\sigma} to 2{\sigma} level. Given the very small number of observed lines of sight, our analysis indicates that current data cannot rule out a purely AGN-driven scenario with high statistical significance.

We present optical and near-infrared imaging of the field of the z=7.0842 quasar ULAS J112001.48+064124.3 taken with the Hubble Space Telescope. We use these data to search for galaxies that may be physically associated with the quasar, using the Lyman break technique, and find three such objects, although the detection of one in Spitzer Space Telescope imaging strongly suggests it lies at z~2. This is consistent with the field luminosity function and indicates that there is no excess of >L* galaxies within 1 Mpc of the quasar. A detection of the quasar shortward of the Ly-alpha line is consistent with the previously observed evolution of the intergalactic medium at z>5.5.

Simulations of structure formation in the Universe predict that galaxies are embedded in a "cosmic web", where the majority of baryons reside as rarefied and highly ionized gas. This material has been studied for decades in absorption against background sources, but the sparseness of these inherently one-dimensional probes preclude direct constraints on the three-dimensional morphology of the underlying web. Here we report observations of a cosmic web filament in Lyman-alpha emission, discovered during a survey for cosmic gas fluorescently "illuminated" by bright quasars at z=2.3. With a projected size of approximately 460 physical kpc, the Lyman-alpha emission surrounding the radio-quiet quasar UM287 extends well beyond the virial radius of any plausible associated dark matter halo. The estimated cold gas mass of the nebula from the observed emission is at least ten times larger than what is typically found by cosmological simulations, suggesting that a population of intergalactic gas clumps with sub-kpc sizes may be missing within current numerical models.

We present a first study of the effect of local photoionising radiation on gas cooling in smoothed particle hydrodynamics simulations of galaxy formation. We explore the combined effect of ionising radiation from young and old stellar populations. The method computes the effect of multiple radiative sources using the same tree algorithm used for gravity, so it is computationally efficient and well resolved. The method foregoes calculating absorption and scattering in favour of a constant escape fraction for young stars to keep the calculation efficient enough to simulate the entire evolution of a galaxy in a cosmological context to the present day. This allows us to quantify the effect of the local photoionisation feedback through the whole history of a galaxy`s formation. The simulation of a Milky Way like galaxy using the local photoionisation model forms ~ 40 % less stars than a simulation that only includes a standard uniform background UV field. The local photoionisation model decreases star formation by increasing the cooling time of the gas in the halo and increasing the equilibrium temperature of dense gas in the disc. Coupling the local radiation field to gas cooling from the halo provides a preventive feedback mechanism which keeps the central disc light and produces slowly rising rotation curves without resorting to extreme feedback mechanisms. These preliminary results indicate that the effect of local photoionising sources is significant and should not be ignored in models of galaxy formation.

We use a set of AMR hydrodynamic simulations post-processed with the radiative-transfer code RADAMESH to study how inhomogeneous HeII reionization affects the intergalactic medium (IGM). We propagate radiation from active galactic nuclei (AGNs) considering two scenarios for the time evolution of the ionizing sources. We find that HeII reionization takes place in a very inhomogeneous fashion, through the production of well separated bubbles of the ionized phase that subsequently percolate. Overall, the reionization process is extended in time and lasts for a redshift interval Delta z>1. At fixed gas density, the temperature distribution is bimodal during the early phases of HeII reionization and cannot be described by a simple effective equation of state. When HeII reionization is complete, the IGM is characterized by a polytropic equation of state with index gamma~1.20. This relation is appreciably flatter than at the onset of the reionization process (gamma=1.56) and also presents a much wider dispersion around the mean. We extract HI and HeII Ly-alpha absorption spectra from the simulations and fit Voigt profiles to them. We find that the regions where helium is doubly ionized are characterized by different probability density functions of the curvature and of the Doppler b parameters of the HI Ly-alpha forest as a consequence of the bimodal temperature distribution during the early phases of HeII reionization. The column-density ratio in HeII and HI shows a strong spatial variability. Its probability density function rapidly evolves with time reflecting the increasing volume fraction in which ionizing radiation is harder due to the AGN contribution. Finally we show that the number density of the flux-transmission windows per unit redshift and the mean size of the dark gaps in the HeII spectra have the potential to distinguish between different reionization scenarios. (abridged)

With close pairs of quasars at different redshifts, a background quasar sightline can be used to study a foreground quasar's environment in absorption. We use a sample of 650 projected quasar pairs to study the HI Lya absorption transverse to luminous, z~2 quasars at proper separations of 30kpc < R < 1Mpc. In contrast to measurements along the line-of-sight, regions transverse to quasars exhibit enhanced HI Lya absorption and a larger variance than the ambient intergalactic medium, with increasing absorption and variance toward smaller scales. Analysis of composite spectra reveals excess absorption characterized by a Lya equivalent width profile W = 2.3A (R/100kpc)^-0.46. We also observe a high (~60%) covering factor of strong, optically thick HI absorbers (HI column log NHI > 17.3) at separations R<200kpc, which decreases to ~20% at R~1Mpc, but still represents a significant excess over the cosmic average. This excess of optically thick absorption can be described by a quasar-absorber cross-correlation function xi_QA(r) = (r/r_0)^gamma with a large correlation length r_0 = 12.5+2.7-1.4 Mpc/h (comoving) and gamma = 1.68+0.14-0.30. The HI absorption measured around quasars exceeds that of any previously studied population, consistent with quasars being hosted by massive dark matter halos Mhalo~10^12.5 Msun at z~2.5. The environments of these massive halos are highly biased towards producing optically thick gas, and may even dominate the cosmic abundance of Lyman limit systems and hence the intergalactic opacity to ionizing photons at z~2.5. The anisotropic absorption around quasars implies the transverse direction is much less likely to be illuminated by ionizing radiation than the line-of-sight, which we interpret in terms of the same obscuration effects frequently invoked in unified models of active galactic nuclei.

A deep narrow-band survey for Ly-alpha emission carried out on the VLT-FORS2 has revealed 98 Ly-alpha candidates down to a flux limit of 4.e-18 erg/s/cm^2 in a volume of 5500 comoving Mpc^3 at z=2.4 centered on the hyperluminous quasar HE0109-3518. The properties of the detected sources in terms of their i) equivalent width distribution, ii) luminosity function, and iii) the average luminosity versus projected distance from the quasar, all suggest that a large fraction of these objects have been fluorescently "illuminated" by HE0109-3518. This conclusion is supported by comparison with detailed radiative transfer simulations of the effects of the quasar illumination. 18 objects have a rest-frame Equivalent Width (EW0) larger than 240A, the expected limit for Ly-alpha emission powered by Population II star formation and 12 sources among these do not have any continuum counterpart in a deep V-band imaging of the same field. For these, a stacking analysis indicates EW0>800A, effectively ruling out Ly-alpha powered by internal star formation. These sources are thus the best candidates so far for proto-galactic clouds or "dark" galaxies at high-redshift, whose existence has recently been suggested by several theoretical studies. Assuming they are mostly ionized by the quasar radiation, we estimate that their gas masses would be about 10^9 Msun implying that their star formation efficiencies (SFE) are less than 1.e-11 yr^-1 several times below the SFE of the most gas-rich dwarf galaxies locally, and five hundred times lower than typical massive star-forming galaxies at z~2. We have also discovered extended, filamentary gas, also likely illuminated by the quasar, around some of the brightest continuum-detected sources with EW0>240A. This emission is compatible with the expectations for circum-galactic cold streams but other origins, including tidal stripping, are also possible.

We present a new three-dimensional radiative transfer (RT) code, RADAMESH, based on a ray-tracing, photon-conserving and adaptive (in space and time) scheme. RADAMESH uses a novel Monte Carlo approach to sample the radiation field within the computational domain on a "cell-by-cell" basis. Thanks to this algorithm, the computational efforts are now focused where actually needed, i.e. within the Ionization-fronts (I-fronts). This results in an increased accuracy level and, at the same time, a huge gain in computational speed with respect to a "classical" Monte Carlo RT, especially when combined with an Adaptive Mesh Refinement (AMR) scheme. Among several new features, RADAMESH is able to adaptively refine the computational mesh in correspondence of the I-fronts, allowing to fully resolve them within large, cosmological boxes. We follow the propagation of ionizing radiation from an arbitrary number of sources and from the recombination radiation produced by H and He. The chemical state of six species (HI, HII, HeI, HeII, HeIII, e) and gas temperatures are computed with a time-dependent, non-equilibrium chemistry solver. We present several validating tests of the code, including the standard tests from the RT Code Comparison Project and a new set of tests aimed at substantiating the new characteristics of RADAMESH. Using our AMR scheme, we show that properly resolving the I-front of a bright quasar during Reionization produces a large increase of the predicted gas temperature within the whole HII region. Also, we discuss how H and He recombination radiation is able to substantially change the ionization state of both species (for the classical Stroemgren sphere test) with respect to the widely used "on-the-spot" approximation.

Current models of galaxy formation lack an efficient and physically constrained mechanism to regulate star formation (SF) in low and intermediate mass galaxies. We argue that the missing ingredient could be the effect of photoionization by local sources on the gas cooling. We show that the soft X-ray and EUV flux generated by SF is able to efficiently remove the main coolants (e.g., HeII, OV and FeIX) from the halo gas via direct photoionization. As a consequence, the cooling and accretion time of the gas surrounding star-forming galaxies may increase by one or two orders of magnitude. For a given halo mass and redshift, the effect is directly related to the value of the star formation rate (SFR). Our results suggest the existence of a critical SFR above which "cold" mode accretion is stopped, even for haloes with virial masses well below the critical shock-heating mass suggested by previous studies.The evolution of the critical SFR with redshift, for a given halo mass, resembles the respective steep evolution of the observed SFR for z<1. This suggests that photoionization by local sources would be able to regulate gas accretion and star formation, without the need for additional, strong feedback processes.

We present a new method to directly map the neutral-hydrogen distribution during the reionization epoch and to constrain the emission properties of the highest-redshift quasars (QSOs). As a tracer of HI, we propose to use the Ly-alpha radiation produced by quasar ionization fronts (I-fronts) that expand in the partially ionized intergalactic medium (IGM) before reionization is complete. These Ly-alpha photons are mainly generated by collisional excitations of hydrogen atoms in the boundary of the rapidly expanding HII region. The observable signal is produced by the part of the I-front that lies behind the QSO with respect to the observer. Combining two radiative transfer models (one for the QSO ionizing radiation and one for the Ly-alpha photons), we estimate the expected Ly-alpha spectral shape and surface brightness (SB) for a large number of configurations where we varied both the properties of the ionizing QSO and of the surrounding medium. We find that the expected signal is observable as a single (broad) emission line with a characteristic width of 100-200 km/s. The expected SB produced at redshift z~6.5 within a fully neutral region (at mean density) by a typical QSO I-front lies in the range $10^{-21}-10^{-20}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ and decreases proportionally to $(1+z)^2$ for a given QSO age. QSOs with harder spectra may produce a significantly brighter emission at early phases. The signal may cover up to a few hundred square arcmin on the sky and should be already detectable with current facilities by means of moderate/high resolution spectroscopy. The detection of this Ly-alpha emission can shed new light on the reionization history, the age and the emission properties of the highest-redshift QSOs. (abridged)

We report the results of a survey for fluorescent Ly-alpha emission carried out in the field surrounding the z=3.1 quasar QSO0420-388 using the FORS2 instrument on the VLT. We first review the properties expected for fluorescent Ly-alpha emitters, compared with those of other non-fluorescent Ly-alpha emitters. Our observational search detected 13 Ly-alpha sources sparsely sampling a volume of ~14000 comoving Mpc^3 around the quasar. The properties of these in terms of i) the line equivalent width, ii) the line profile and iii) the value of the surface brightness related to the distance from the quasar, all suggest that several of these may be plausibly fluorescent. Moreover, their number is in good agreement with the expectation from theoretical models. One of the best candidates for fluorescence is sufficiently far behind QSO0420-388 that it would imply that the quasar has been active for (at least) ~60 Myrs. Further studies on such objects will give information about proto-galactic clouds and on the radiative history (and beaming) of the high-redshift quasars.

We combine a high-resolution hydro-simulation of the LambdaCDM cosmology with two radiative transfer schemes (for continuum and line radiation) to predict the properties, spectra and spatial distribution of fluorescent Ly-alpha emission at z~3. We focus on line radiation produced by recombinations in the dense intergalactic medium ionized by UV photons. In particular, we consider both a uniform background and the case where gas clouds are illuminated by a nearby quasar. We find that the emission from optically thick regions is substantially less than predicted from the widely used static, plane-parallel model. The effects induced by a realistic velocity field and by the complex geometric structure of the emitting regions are discussed in detail. We make predictions for the expected brightness and size distributions of the fluorescent sources.Our results account for recent null detections and can be used to plan new observational campaigns both in the field (to measure the intensity of the diffuse UV background) and in the proximity of bright quasars (to understand the origin of high colum-density absorbers).