The discovery of quasars few hundred megayears after the Big Bang represents a major challenge to our understanding of black holes and galaxy formation and evolution. Their luminosity is produced by extreme gas accretion onto black holes, which already reached masses of 10$^9$ M$_\odot$ by z ~ 6. Simultaneously, their host galaxies form hundreds of stars per year, using up gas in the process. To understand which environments are able to sustain the rapid formation of these extreme sources we started a VLT/MUSE effort aimed at characterizing the surroundings of a sample of 5.7 < z < 6.6 quasars dubbed: the Reionization Epoch QUasar InvEstigation with MUSE (REQUIEM) survey. We here present results of our searches for extended Ly-Alpha halos around the first 31 targets observed as part of this program. Reaching 5-sigma surface brightness limits of 0.1-1.1 x 10$^{-17}$ erg/s/cm$^2$/arcsec$^2$ over a 1 arcsec$^2$ aperture, we were able to unveil the presence of 12 Ly-Alpha nebulae, 8 of which are newly discovered. The detected nebulae show a variety of emission properties and morphologies with luminosities ranging from 8 x 10$^{42}$ to 2 x 10$^{44}$ erg/s, FWHMs between 300 and 1700 km/s, sizes < 30 pkpc, and redshifts consistent with those of the quasar host galaxies. As the first statistical and homogeneous investigation of the circum-galactic medium of massive galaxies at the end of the reionization epoch, the REQUIEM survey enables the study of the evolution of the cool gas surrounding quasars in the first 3 Gyr of the Universe. A comparison with the extended Ly-Alpha emission observed around bright (M$_{1450}$ < -25 mag) quasars at intermediate redshift indicates little variations on the properties of the cool gas from z ~ 6 to z ~ 3 followed by a decline in the average surface brightness down to z ~ 2.

New and unique science opportunities in several different fields of astrophysics are offered by conducting spectroscopic studies of the Universe in the ultraviolet (UV), a wavelength regime that is not accessible from the ground. We present some of the scientific challenges that can be addressed with a space-based mission in 2035 - 2050. (1) By detecting the intergalactic medium in emission it will be possible to unveil the cosmic web, whose existence is predicted by current theories of structure formation, but it has not been probed yet. (2) Observations of the neutral gas distribution (by mapping the Lyman-alpha emission) in low-redshift galaxy cluster members will clarify the efficiency with which ram-pressure stripping removes the gas from galaxies and the role of the environment in quenching star formation. (3) By observing statistical samples of supernovae in the UV it will be possible to characterize the progenitor population of core-collapse supernovae, providing the initial conditions for any forward-modeling simulation and allowing the community to progress in the understanding of the explosion mechanism of stars and the final stages of stellar evolution. (4) Targeting populations of accreting white dwarfs in globular clusters it will be possible to constrain the evolution and fate of these stars and investigate the properties of the most compact systems with the shortest orbital periods which are expected to be the brightest low frequency gravitational wave sources. A UV-optimized telescope (wavelength range ~ 90 - 350 nm), equipped with a panoramic integral field spectrograph with a large field of view (FoV ~ 1 x 1 arcmin^2), with medium spectral (R = 4000) and spatial (~ 1" - 3") resolution will allow the community to simultaneously obtain spectral and photometric information of the targets, and tackle the science questions presented in this paper.

We have used the Multi-Unit Spectroscopic Explorer (MUSE) to study the circumgalactic medium (CGM) of a z = 2.92 radio galaxy, MRC 0943-242 by parametrising its emitting and absorbing gas. In both Ly$\alpha$ $\lambda$1216 and He II $\lambda$1640 lines, we observe emission with velocity shifts of $\Delta v \simeq-1000$ km s$^{-1}$ from the systemic redshift of the galaxy. These blueshifted components represent kinematically perturbed gas that is aligned with the radio axis, which we interpret as jet-driven outflows. Three of the four known Ly$\alpha$ absorbers are detected at the same velocity as C IV $\lambda\lambda1548,1551$ and N V $\lambda\lambda1239,1243$ absorbers, proving that the gas is metal enriched more so than previously thought. At the velocity of a strong Ly$\alpha$ absorber with an HI column of $N_{\rm HI}/{\rm cm}^{-2} = 10^{19.2}$ and velocity shift of $\Delta v \simeq -400$ km s$^{-1},$ we also detect Si II $\lambda$1260 and Si II $\lambda$1527 absorption, which suggests that the absorbing gas is ionisation bounded. With the added sensitivity of this MUSE observation, we are more capable of adding constraints to absorber column densities and consequently determining what powers their ionisation. To do this, we obtain photoionisation grid models in \pkg{cloudy} which show that AGN radiation is capable of ionising the gas and producing the observed column densities in a gas of metallicity of Z/Z$_\odot \simeq$ 0.01 with a nitrogen abundance a factor of 10 greater than that of hydrogen. This metal-enriched absorbing gas, which is also spatially extended over a projected distance of $r \gtrsim 60$ kpc, is likely to have undergone chemical enrichment through stellar winds that have swept up metals from the interstellar-medium and deposited them in the outer regions of the galaxy's halo.

We recently discovered that the luminous radio-quiet QSO LBQS0302-0019 at z=3.286 is likely accompanied by an obscured AGN at 20kpc projected distance, which we dubbed Jil. It represents the tightest candidate obscured/unobscured dual AGN system at z>3. To verify the dual AGN scenario we obtained deep $K_s$ band (rest-frame $V$ band) imaging with the VLT/HAWK-I+GRAAL instrument at 0.4" resolution during science verification in January 2018. Indeed, we detect the individual host galaxies of the QSO and Jil with estimated stellar masses of $\log(M_\star/M_{\odot})=11.4\pm0.5$ and $\log(M_\star/M_{\odot})=10.9\pm0.5$, respectively. Near-IR spectra obtained with VLT-KMOS reveal a clear [OIII] $\lambda$5007 line detection at the location of Jil which does not contribute significantly to the $K_s$ band flux. Both observations therefore corroborate the dual AGN scenario. A comparison to Illustris cosmological simulations suggests a parent halo mass of $\log(M_\mathrm{halo}/M_{\odot})=13.2\pm0.5$ for this interacting galaxy system, corresponding to a very massive dark matter halo at that epoch.

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.

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.

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.