We present a new implementation for active galactic nucleus (AGN) feedback through small-scale, ultra-fast winds in the moving-mesh hydrodynamic code AREPO. The wind is injected by prescribing mass, momentum and energy fluxes across a spherical boundary centred on a supermassive black hole according to available constraints for accretion disc winds. After sweeping-up a mass equal to their own, small-scale winds thermalise, powering energy-driven outflows with dynamics, structure and cooling properties in excellent agreement with those of analytic wind solutions. Momentum-driven solutions do not easily occur, because the Compton cooling radius is usually much smaller than the free-expansion radius of the small-scale winds. Through various convergence tests, we demonstrate that our implementation yields wind solutions which are well converged down to the typical resolution achieved in cosmological simulations. We test our model in hydrodynamic simulations of isolated Milky Way - mass galaxies. Above a critical AGN luminosity, initially spherical, small-scale winds power bipolar, energy-driven super-winds that break out of the galactic nucleus, flowing at speeds $> 1000 \rm \, km \, s^{-1}$ out to $\sim 10 \, \rm kpc$. These energy-driven outflows result in moderate, but long-term, reduction in star formation, which becomes more pronounced for higher AGN luminosities and faster small-scale winds. Suppression of star formation proceeds through a rapid mode that involves the removal of the highest-density, nuclear gas and through a slower mode that effectively halts halo gas accretion. Our new implementation makes it possible to model AGN-driven winds in a physically meaningful and validated way in simulations of galaxy evolution, the interstellar medium and black hole accretion flows.

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.

Quasars at $z \,=\, 6$ are powered by accretion onto supermassive black holes with masses $M_{\rm BH} \sim 10^9 \rm \, M_{\odot}$. Their rapid assembly requires efficient gas inflow into the galactic nucleus, sustaining black hole accretion at a rate close to the Eddington limit, but also high central star formation rates. Using a set of cosmological 'zoom-in' hydrodynamic simulations performed with the moving mesh code Arepo, we show that $z \,=\, 6$ quasar host galaxies develop extremely tightly bound stellar bulges with peak circular velocities $300$ - $500$ km s$^{-1}$ and half-mass radii $\approx 0.5 \, \rm kpc$. Despite their high binding energy, we find that these compact bulges expand at $z \, < \, 6$, with their half-mass radii reaching $ \approx 5$ kpc by $z \, = \, 3$. The circular velocity drops by factors $\approx 2$ from their initial values to $200$ - $300$ km s$^{-1}$ at $z \, \approx \, 3$ and the stellar profile undergoes a cusp-core transformation. By tracking individual stellar populations, we find that the gradual expansion of the stellar component is mainly driven by fluctuations in the gravitational potential induced by bursty AGN feedback. We also find that galaxy size growth and the development of a cored stellar profile does not occur if AGN feedback is ineffective. Our findings suggest that AGN-driven outflows may have profound implications for the internal structure of massive galaxies, possibly accounting for their size growth, the formation of cored ellipticals as well as for the saturation of the $M_{\rm BH}$ - $\sigma_{\star}$ seen at high velocity dispersions $\sigma_{\star}$.

Using cosmological, radiation-hydrodynamic simulations targeting a rare $\approx \, 2 \times 10^{12} \, \rm M_\odot$ halo at $z \, = \, 6$, we show that the number counts and internal properties of satellite galaxies within the massive halo are sensitively regulated by a combination of local stellar radiative feedback and strong tidal forces. Radiative feedback operates before the first supernova explosions erupt and results in less tightly-bound galaxies. Satellites are therefore more vulnerable to tidal stripping when they accrete onto the main progenitor and are tidally disrupted on a significantly shorter timescale. Consequently, the number of satellites with $M_{\rm \star} > 10^{7} \, \rm M_\odot$ within the parent system's virial radius drops by up to $60 \%$ with respect to an identical simulation performed without stellar radiative feedback. Radiative feedback also impacts the central galaxy, whose effective radius increases by a factor $\lesssim 3$ due to the presence of a more extended and diffuse stellar component. We suggest that the number of satellites in the vicinity of massive high-redshift galaxies is an indication of the strength of stellar radiative feedback and and can be anomalously low in the extreme cosmic environments of high-redshift quasars.

We study molecular outflows in a sample of 45 local galaxies, both star forming and AGN, primarily by using CO data from the ALMA archive and from the literature. For a subsample we also compare the molecular outflow with the ionized and neutral atomic phases. We infer an empirical analytical function relating the outflow rate simultaneously to the SFR, $L_{\rm AGN}$, and galaxy stellar mass; this relation is much tighter than the relations with the individual quantities. The outflow kinetic power shows a larger scatter than in previous, more biased studies, spanning from 0.1 to 5 per cent of $L_{\rm AGN}$, while the momentum rate ranges from 1 to 30 times $L_{\rm AGN}/c$, indicating that these outflows can be both energy-driven, but with a broad range of coupling efficiencies with the ISM, and radiation pressure-driven. For about 10 per cent of the objects the outflow energetics significantly exceed the maximum theoretical values; we interpret these as 'fossil outflows' resulting from activity of a past strong AGN, which has now faded. We estimate that, in the stellar mass range probed here ($>$ 10$^{10}~\rm M_{\odot}$), less than 5 per cent of the outflowing gas escapes the galaxy. The molecular gas depletion time associated with the outflow can be as short as a few million years in powerful AGN, however, the total gas (H$_2$+HI) depletion times are much longer. Altogether, our findings suggest that even AGN-driven outflows might be relatively ineffective in clearing galaxies of their entire gas content, although they are likely capable of clearing and quenching the central region.

We present cosmological radiation-hydrodynamic simulations, performed with the code Ramses-RT, of radiatively-driven outflows in a massive quasar host halo at $z = 6$. Our simulations include both single- and multi-scattered radiation pressure on dust from a quasar and are compared against simulations performed with thermal feedback. For radiation pressure-driving, we show that there is a critical quasar luminosity above which a galactic outflow is launched, set by the equilibrium of gravitational and radiation forces. While this critical luminosity is unrealistically high in the single-scattering limit for plausible black hole masses, it is in line with a $\approx 3 \times 10^9 \, \rm M_\odot$ black hole accreting at its Eddington limit, if infrared (IR) multi-scattering radiation pressure is included. The outflows are fast ($v \, \gtrsim \, 1000 \, \rm km \, s^{-1}$) and strongly mass-loaded with peak mass outflow rates $\approx 10^3 - 10^4 \, \rm M_\odot \, yr^{-1}$, but short-lived ($< 10 \, \rm Myr$). Outflowing material is multi-phase, though predominantly composed of cool gas, forming via a thermal instability in the shocked swept-up component. Radiation pressure- and thermally-driven outflows both affect their host galaxies significantly, but in different, complementary ways. Thermally-driven outflows couple more efficiently to diffuse halo gas, generating more powerful, hotter and more volume-filling outflows. IR radiation, through its ability to penetrate dense gas via diffusion, is more efficient at ejecting gas from the bulge. The combination of gas ejection through outflows with internal pressurisation by trapped IR radiation leads to a complete shut down of star formation in the bulge. We hence argue that radiation pressure-driven feedback may be an important ingredient in regulating star formation in compact starbursts, especially during the quasar's `obscured' phase.

We present radiation-hydrodynamic simulations of radiatively-driven gas shells launched by bright active galactic nuclei (AGN) in isolated dark matter haloes. Our goals are (1) to investigate the ability of AGN radiation pressure on dust to launch galactic outflows and (2) to constrain the efficiency of infrared (IR) multi-scattering in boosting outflow acceleration. Our simulations are performed with the radiation-hydrodynamic code RAMSES-RT and include both single- and multi-scattered radiation pressure from an AGN, radiative cooling and self-gravity. Since outflowing shells always eventually become transparent to the incident radiation field, outflows that sweep up all intervening gas are likely to remain gravitationally bound to their halo even at high AGN luminosities. The expansion of outflowing shells is well described by simple analytic models as long as the shells are mildly optically thick to IR radiation. In this case, an enhancement in the acceleration of shells through IR multi-scattering occurs as predicted, i.e. a force dP/dt = tau_IR L/c is exerted on the gas. For high optical depths tau_IR > 50, however, momentum transfer between outflowing optically thick gas and IR radiation is rapidly suppressed, even if the radiation is efficiently confined. At high tau_IR, the characteristic flow time becomes shorter than the required trapping time of IR radiation such that the momentum flux dP/dt << tau_IR L/c. We argue that while unlikely to unbind massive galactic gaseous haloes, AGN radiation pressure on dust could play an important role in regulating star formation and black hole accretion in the nuclei of massive compact galaxies at high redshift.

It is the twentieth anniversary of the publication of the seminal papers by Magorrian et al. and Silk & Rees which, along with other related work, ignited an explosion of publications connecting active galactic nuclei (AGN)-driven outflows to galaxy evolution. With a surge in observations of AGN outflows, studies are attempting to directly test AGN feedback models using the outflow properties. With a focus on outflows traced by optical and CO emission lines, we discuss significant challenges which greatly complicate this task from both an observational and theoretical perspective. We highlight observational uncertainties involved, and the assumptions required, when deriving kinetic coupling efficiencies (i.e., outflow kinetic power as a fraction of AGN luminosity) from typical observations. Based on recent models we demonstrate that extreme caution should taken when comparing observationally-derived kinetic coupling efficiencies to coupling efficiencies from fiducial feedback models.

We present the results from a $\sim500$ ks Chandra observation of the $z=6.31$ QSO SDSS J1030+0524. This is the deepest X-ray observation to date of a $z\sim6$ QSO. The QSO is detected with a total of 125 net counts in the full ($0.5-7$ keV) band and its spectrum can be modeled by a single power-law model with photon index of $\Gamma = 1.81 \pm 0.18$ and full band flux of $f=3.95\times 10^{-15}$ erg s$^{-1}$ cm$^{-2}$. When compared with the data obtained by XMM-Newton in 2003, our Chandra observation in 2017 shows a harder ($\Delta \Gamma \approx -0.6$) spectrum and a 2.5 times fainter flux. Such a variation, in a timespan of $\sim2$ yrs rest-frame, is unexpected for such a luminous QSO powered by a $> 10^9 \: M_{\odot}$ black hole. The observed source hardening and weakening could be related to an intrinsic variation in the accretion rate. However, the limited photon statistics does not allow us to discriminate between an intrinsic luminosity and spectral change, and an absorption event produced by an intervening gas cloud along the line of sight. We also report the discovery of diffuse X-ray emission that extends for 30"x20" southward the QSO with a signal-to-noise ratio of $\sim$6, hardness ratio of $HR=0.03_{-0.25}^{+0.20}$, and soft band flux of $f_{0.5-2 \: keV}= 1.1_{-0.3}^{+0.3} \times 10^{-15}$ erg s$^{-1}$ cm$^{-2}$, that is not associated to a group or cluster of galaxies. We discuss two possible explanations for the extended emission, which may be either associated with the radio lobe of a nearby, foreground radio galaxy (at $z \approx 1-2$), or ascribed to the feedback from the QSO itself acting on its surrounding environment, as proposed by simulations of early black hole formation.

We obtained ALMA spectroscopy and imaging to investigate the origin of the unexpected sub-mm emission toward the most distant quiescent galaxy known to date, ZF-COSMOS-20115 at z=3.717. We show here that this sub-mm emission is produced by another massive, compact and extremely obscured galaxy, located only 3.1 kpc away from the quiescent galaxy. We dub the quiescent and dusty galaxies Jekyll and Hyde, respectively. No dust emission is detected at the location of Jekyll, implying SFR < 13 Msun/yr, which is the most stringent upper limit ever obtained for a quiescent galaxy at these redshifts. The two sources are confirmed to lie at the same redshift thanks to the detection of [CII]158 in Hyde, which provides one the few robust redshifts for an "H-dropout" galaxy. The line has a rotating-disk velocity profile blueshifted from Jekyll by 549+/-60 km/s, demonstrating that it is produced by another galaxy. Careful de-blending of the Spitzer imaging confirms the existence of Hyde, and its non-detection with Hubble requires extreme attenuation by dust. Modeling the photometry of both galaxies shows that Jekyll has fully quenched >200 Myr prior to observation and still presents a challenge for models, while Hyde only harbors moderate star-formation (SFR<120 Msun/yr) and is located at least a factor 1.4 below the z~4 main sequence. Hyde could also have stopped forming stars <200 Myr before being observed, which would be consistent with its hight compactness similar to z~4 quiescent galaxies and its low [CII]/FIR ratio, but significant SF cannot be ruled out. Finally, we show that Hyde hosts a dense reservoir of gas comparable to that of extreme starbursts, suggesting that its SFR was reduced without expelling the gas outside of the galaxy. We argue that Jekyll and Hyde can be seen as two stages of quenching, and provide a unique laboratory to study this poorly understood phenomenon. [abridged]

We employ cosmological hydrodynamical simulations to investigate models in which the supermassive black holes (BHs) powering luminous z ~ 6 QSOs grow from massive seeds. We simulate at high resolution 18 fields sampling regions with densities ranging from the mean cosmic density all the way to the highest sigma peaks in the Millennium simulation volume. Only in the most massive halos, BHs can grow to masses up to ~ 10^9 Msun by z ~ 6 without invoking super-Eddington accretion. Accretion onto the most massive BHs becomes limited by thermal AGN feedback by z ~ 9-8 with further BH growth proceeding in short Eddington limited bursts. Our modelling suggests that current flux-limited surveys of QSOs at high redshift preferentially detect objects at their peak luminosity and therefore miss a substantial population of QSOs powered by similarly massive BHs but with low accretion rates. To test whether the required host halo masses are consistent with the observed galaxy environment of z ~ 6 QSOs, we produce realistic rest-frame UV images of our simulated galaxies. Without strong stellar feedback, our simulations predict numbers of bright galaxies larger than observed by a factor ten or more. Supernova-driven galactic winds reduce the predicted numbers to a level consistent with observations indicating that stellar feedback was already very efficient at high redshifts. We have further investigated the effect of thermal AGN feedback on the surrounding gas. Our adopted AGN feedback prescription drives mostly energy-driven highly anisotropic outflows with gas speeds of >= 1000 km/s to distances of >= 10 kpc consistent with observations. The spatially extended thermal X-ray emission around bright QSOs powered by these outflows can exceed by large factors the emission expected without AGN feedback and is an important diagnostic of the mechanism whereby AGN feedback energy couples to surrounding gas.

Observations of the emission from spatially extended cold gas around bright high-redshift QSOs reveal surprisingly large velocity widths exceeding 2000 km s^(-1), out to projected distances as large as 30 kpc. The high velocity widths have been interpreted as the signature of powerful AGN-driven outflows. Naively, these findings appear in tension with hydrodynamic models in which AGN-driven outflows are energy-driven and thus very hot with typical temperatures T = 10^6-7 K. Using the moving-mesh code Arepo, we perform 'zoom-in' cosmological simulations of a z = 6 QSO and its environment, following black hole growth and feedback via energy-driven outflows. In the simulations, the QSO host galaxy is surrounded by a clumpy circum-galactic medium pre-enriched with metals due to supernovae-driven galactic outflows. As a result, part of the AGN-driven hot outflowing gas can cool radiatively, leading to large amounts (> 10^9 M_sun) of cold gas comoving with the hot bipolar outflow. This results in velocity widths of spatially extended cold gas similar to those observed. We caution, however, that gas inflows, random motions in the deep potential well of the QSO host galaxy and cooling of supernovae-driven winds contribute significantly to the large velocity width of the cold gas in the simulations, complicating the interpretation of observational data.

We employ hydrodynamical simulations using the moving-mesh code AREPO to investigate the role of energy and momentum input from Active Galactic Nuclei (AGN) in driving large-scale galactic outflows. We start by reproducing analytic solutions for both energy- and momentum-driven outflowing shells in simulations of a spherical isolated dark matter potential with gas in hydrostatic equilibrium and with no radiative cooling. We confirm that for this simplified setup, galactic outflows driven by a momentum input rate of order L_Edd/c can establish an M_BH - sigma relation with slope and normalisation similar to that observed. We show that momentum input at a rate of L_Edd/c is however insufficient to drive efficient outflows once cooling and gas inflows as predicted by cosmological simulations at resolved scales are taken into account. We argue that observed large-scale AGN-driven outflows are instead likely to be energy-driven and show that such outflows can reach momentum fluxes exceeding 10 L_Edd/c within the innermost 10 kpc of the galaxy. The outflows are highly anisotropic, with outflow rates and a velocity structure found to be inadequately described by spherical outflow models. We verify that the hot energy-driven outflowing gas is expected to be strongly affected by metal-line cooling, leading to significant amounts (>10^9 M_sun) of entrained cold gas.