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.

The objective of this work is to study how active galactic nuclei (AGN) influence star formation in host galaxies. We present a detailed investigation of the star-formation history and conditions of a $z=2.57$ massive radio galaxy based on VLT/X-SHOOTER and ALMA observations. The deep rest-frame ultraviolet spectrum contains photospheric absorption lines and wind features indicating the presence of OB-type stars. The most significantly detected photospheric features are used to characterize the recent star formation: neither instantaneous nor continuous star-formation history is consistent with the relative strength of the Si II $\lambda$1485 and S V $\lambda$1502 absorption. Rather, at least two bursts of star formation took place in the recent past, at $6^{+1}_{-2}$ Myr and $\gtrsim20$ Myr ago, respectively. We deduce a molecular H$_{2}$ gas mass of $(3.9\pm1.0)\times10^{10}$ M$_{\odot}$ based on ALMA observations of the [C I] $^3$P$_{2}$-$^3$P$_{1}$ emission. The molecular gas mass is only 13 % of its stellar mass. Combined with its high star-formation rate of ($1020^{+190}_{-170}$) Myr, this implies a high star-formation efficiency of $(26\pm8$) Gyr$^{-1}$ and a short depletion time of $(38\pm12)$ Myr. We attribute the efficient star formation to compressive gas motions in order to explain the modest velocity dispersions ($\leqslant$ 55 km s$^{-1}$) of the photospheric lines and of the star-forming gas traced by [C I]. Because of the likely very young age of the radio source, our findings suggest that vigorous star formation consumes much of the gas and works in concert with the AGN to remove any residual molecular gas, and eventually quenching star formation in massive galaxies.

We present 0.3" (band 6) and 1.5" (band 3) ALMA observations of the (sub)millimeter dust continuum emission for 25 radio galaxies at 1<z<5.2. Our survey reaches a rms flux density of ~50$\mu$Jy in band 6 and ~20$\mu$Jy in band 3. This is an order of magnitude deeper than single-dish 850 $\mu$m observations, and reaches fluxes where synchrotron and thermal dust emission are expected to be of the same order of magnitude. Combining our sensitive ALMA observations with radio data from ATCA, VLA, and IR photometry from Herschel and Spitzer, we have disentangled the synchrotron and thermal dust emission. We determine the star-formation rates (SFR) and AGN IR luminosities using our newly developed spectral energy distribution fitting code MrMoose. We find that synchrotron emission contributes substantially at ~1 mm. Through our sensitive flux limits and accounting for a contribution from synchrotron emission in the mm, we revise downward the median SFR by a factor of 7 compared to previous estimates based solely on Herschel and Spitzer data. The hosts of these radio-loud AGN appear predominantly below the main sequence of star-forming galaxies, indicating that the star formation in many of the host galaxies has been quenched. Future growth of the host galaxies without substantial black hole mass growth will be needed to bring these objects on the local relation between the supermassive black holes and their host galaxies. Given the mismatch in the timescales of any star formation that took place in the host galaxies and lifetime of the AGN, we hypothesize that a key role is played by star formation in depleting the gas before the action of the powerful radio jets quickly drives out the remaining gas. This positive feedback loop of efficient star formation rapidly consuming the gas coupled to the action of the radio jets in removing the residual gas is how massive galaxies are rapidly quenched.

We present a multiwavelength study of the AGN-host starburst galaxy PKS 0529-549 at z~2.6. We use (1) new ALMA observations of the dust continuum and of the [CI] 370 um line, tracing molecular gas, (2) SINFONI spectroscopy of the [OIII] 5007 Ang line, tracing ionized gas, and (3) ATCA radio continuum images, tracing synchrotron emission. Both [CI] and [OIII] show regular velocity gradients, but their systemic velocities and position angles differ by ~300 km/s and ~30 degrees, respectively. The [CI] is consistent with a rotating disc, aligned with the dust and stellar continuum, while the [OIII] likely traces an outflow, aligned with two AGN-driven radio lobes. We model the [CI] cube using 3D disc models, which give best-fit rotation velocities V~310 km/s and velocity dispersions sigma<30 km/s. Hence, the [CI] disc has V/sigma>10 and is not particularly turbulent, similar to local galaxy discs. The dynamical mass (~10^11 Msun) is comparable to the baryonic mass within the errors. This suggests that baryons dominate the inner galaxy dynamics, similar to massive galaxies at z=0. Remarkably, PKS 0529-549 lies on the local baryonic Tully-Fisher relation, indicating that at least some massive galaxies are already in place and kinematically relaxed at z~2.6. This work highlights the potential of the [CI] line to trace galaxy dynamics at high z, as well as the importance of multiwavelength data to interpret gas kinematics.

We present the public release of MrMoose, a fitting procedure that is able to perform multi-wavelength and multi-object spectral energy distribution (SED) fitting in a Bayesian framework. This procedure is able to handle a large variety of cases, from an isolated source to blended multi-component sources from an heterogeneous dataset (i.e. a range of observation sensitivities and spectral/spatial resolutions). Furthermore, MrMoose handles upper-limits during the fitting process in a continuous way allowing models to be gradually less probable as upper limits are approached. The aim is to propose a simple-to-use, yet highly-versatile fitting tool fro handling increasing source complexity when combining multi-wavelength datasets with fully customisable filter/model databases. The complete control of the user is one advantage, which avoids the traditional problems related to the "black box" effect, where parameter or model tunings are impossible and can lead to overfitting and/or over-interpretation of the results. Also, while a basic knowledge of Python and statistics is required, the code aims to be sufficiently user-friendly for non-experts. We demonstrate the procedure on three cases: two artificially-generated datasets and a previous result from the literature. In particular, the most complex case (inspired by a real source, combining Herschel, ALMA and VLA data) in the context of extragalactic SED fitting, makes MrMoose a particularly-attractive SED fitting tool when dealing with partially blended sources, without the need for data deconvolution.

Using MUSE on the ESO-VLT, we obtained a 4 hour exposure of the z=3.12 radio galaxy MRC0316-257. We detect features down to ~10^-19 erg/s/cm^2/arcsec^2 with the highest surface brightness regions reaching more than a factor of 100 higher. We find Ly-alpha emission out to ~250 kpc in projection from the active galactic nucleus (AGN). The emission shows arc-like morphologies arising at 150-250 kpc from the nucleus in projection with the connected filamentary structures reaching down into the circum-nuclear region. The most distant arc is offset by 700 km/s relative to circum-nuclear HeII 1640 emission, which we assume to be at the systemic velocity. As we probe emission closer to the nucleus, the filamentary emission narrows in projection on the sky, the relative velocity decreases to ~250 km/s, and line full-width at half maximum range from 300-700 km/s. From UV line ratios, the emission on scales of 10s of kpc from the nucleus along a wide angle in the direction of the radio jets is clearly excited by the radio jets and ionizing radiation of the AGN. Assuming ionization equilibrium, the more extended emission outside of the axis of the jet direction would require 100% or more illumination to explain the observed surface brightness. High speed (>300 km/s) shocks into rare gas would provide sufficiently high surface brightness. We discuss the possibility that the arcs of Ly-alpha emission represent accretion shocks and the filamentary emission represent gas flows into the halo, and compare our results with gas accretion simulations.