The gas content of galaxies is a key factor for their growth, starting from star formation and black hole accretion to galaxy mergers. Thus, characterising its properties via observations of tracers like the CO emission line is of big importance in order to understand the bigger picture of galaxy evolution. We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of dust continuum, CO(5-4) and CO(8-7) line emission in the quasar--star-forming companion system SMM J04135+10277 (z=2.84). Earlier low-$J$ CO studies of this system found a huge molecular gas reservoir associated to the companion galaxy, while the quasar appeared gas-poor. Our CO observations revealed that the host galaxy of the quasar is also gas-rich, with an estimated molecular gas mass of $\sim(0.7-2.3)\times10^{10}$ M$_{\odot}$. The CO line profiles of the companion galaxy are broad ($\sim1000$ km s$^{-1}$), and show signs of rotation of a compact, massive system. In contrast to previous far-infrared observations, we resolve the continuum emission and detect both sources, with the companion galaxy dominating the dust continuum and the quasar having a $\sim25\%$ contribution to the total dust emission. By fitting the infrared spectral energy distribution of the sources with \textsc{MR-MOOSE} and empirical templates, the infrared luminosities of the quasar and the companion are in the range of $L_{\rm IR, QSO}\sim(2.1-9.6)\times10^{12}$ L$_{\odot}$ and $L_{\rm IR, Comp.}\sim(2.4-24)\times10^{12}$ L$_{\odot}$, while the estimated star formation rates are $\sim210-960$ M$_{\odot}$ yr$^{-1}$ and $\sim240-2400$ M$_{\odot}$ yr$^{-1}$, respectively. Our results demonstrate that non-detection of low-$J$ CO transition lines in similar sources does not necessarily imply the absence of massive molecular gas reservoir but that the excitation conditions favour the excitation of high-$J$ transitions.

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

We report the results of a pilot study of CO$(4-3)$ emission line of three {\it WISE}-selected hyper-luminous, dust-obscured quasars (QSOs) with sensitive ALMA Band 3 observations. These obscured QSOs with $L_{\rm bol}>10^{14}L_\odot$ are among the most luminous objects in the universe. All three QSO hosts are clearly detected both in continuum and in CO$(4-3)$ emission line. Based on CO$(4-3)$ emission line detection, we derive the molecular gas masses ($\sim 10^{10-11}$ M$_\odot$), suggesting that these QSOs are gas-rich systems. We find that three obscured QSOs in our sample follow the similar $L'_{\rm CO}- L_{\rm FIR}$ relation as unobscured QSOs at high redshifts. We also find the complex velocity structures of CO$(4-3)$ emission line, which provide the possible evidence for gas-rich merger in W0149+2350 and possible molecular outflow in W0220+0137 and W0410$-$0913. Massive molecular outflow can blow away the obscured interstellar medium (ISM) and make obscured QSOs evolve towards the UV/optical bright, unobscured phase. Our result is consistent with the popular AGN feedback scenario involving the co-evolution between the SMBH and host galaxy.

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.

(Abridged) In the last decade several massive molecular gas reservoirs were found <100 kpc distance from active galactic nuclei (AGNs), residing in gas-rich companion galaxies. The study of AGN-gas-rich companion systems opens the opportunity to determine whether the stellar mass of massive local galaxies was formed in their host after a merger event or outside of their host galaxy in a close starbursting companion and later incorporated via mergers. We study the quasar-companion galaxy system of SMM J04135+10277 (z=2.84) and investigate the expected frequency of quasar-starburst galaxy pairs at high redshift using a cosmological galaxy formation model. We use archive data and new APEX ArTeMiS data to construct and model the spectral energy distribution of SMM J04135. We also carry out a comprehensive analysis of the cosmological galaxy formation model GALFORM with the aim of characterising how typical the system of SMM J04135 is and whether quasar-star-forming galaxy pairs may constitute an important stage in galaxy evolution. The companion galaxy of SMM J04135 is a heavily dust-obscured starburst galaxy with a median star formation rate (SFR) of $700\,\rm{M_{\odot}\,yr^{-1}}$, median dust mass of $5.1\times 10^9\,\rm{M_{\odot}}$ and median dust luminosity of $\textrm 9.3 \times 10^{12}\,\rm{L_{\odot}}$. Our simulations, performed at z=2.8, suggest that SMM J04135 is not unique. In fact, at a distance of <100 kpc, 22% of our simulated quasar sample have at least one companion galaxy of a stellar mass $>10^8\, \rm{M_{\odot}}$, and 0.3% have at least one highly star-forming companion ($\rm{SFR}>100\,\rm{M_{\odot}\,yr^{-1}}$). Our results suggest that quasar-gas-rich companion galaxy systems are common phenomena in the early Universe and the high incidence of companions makes the study of such systems crucial to understand the growth and hierarchical build-up of galaxies and black holes.

We compare the kinetic energy and momentum injection rates from intense star formation, bolometric AGN radiation, and radio jets with the kinetic energy and momentum observed in the warm ionized gas in 24 powerful radio galaxies at z~2. These galaxies are amongst our best candidates for being massive galaxies near the end of their active formation period, when intense star formation, quasar activity, and powerful radio jets all co-exist. All galaxies have VLT/SINFONI imaging spectroscopy of the rest-frame optical line emission, showing emission-line regions with large velocity offsets (up to 1500 km/s) and line widths (typically 800-1000 km/s) consistent with very turbulent, often outflowing gas. As part of the HeRGE sample, they also have FIR estimates of the star formation and quasar activity obtained with Herschel/PACS and SPIRE, which enables us to measure the relative energy and momentum release from each of the three main sources of feedback in massive, star-forming AGN host galaxies during their most rapid formation phase. We find that star formation falls short by factors 10-1000 of providing the energy and momentum necessary to power the observed gas kinematics. The obscured quasars in the nuclei of these galaxies provide enough energy and momentum in about half of the sample, however, only if these are transfered to the gas relatively efficiently. We compare with theoretical and observational constraints on the efficiency of the energy and momentum transfer from jet and AGN radiation, which advocates that the radio jet is the main driver of the gas kinematics.

We measure the sizes of redshift ~2 star-forming galaxies by stacking data from the Atacama Large Millimeter/submillimeter Array (ALMA). We use a uv-stacking algorithm in combination with model fitting in the uv-domain and show that this allows for robust measures of the sizes of marginally resolved sources. The analysis is primarily based on the 344 GHz ALMA continuum observations centred on 88 sub-millimeter galaxies in the LABOCA ECDFS Submillimeter Survey (ALESS). We study several samples of galaxies at z~2 with $M_* \sim{} 5\times{}10^{10} $ M$_\odot$ , selected using near-infrared photometry (distant red galaxies, extremely red objects, sBzK-galaxies, and galaxies selected on photometric redshift). We find that the typical sizes of these galaxies are ~0.6 arcsec which corresponds to ~5 kpc at z~2, this agrees well with the median sizes measured in the near-infrared z-band (~0.6 arcsec). We find errors on our size estimates of ~0.1-0.2 arcsec, which agree well with the expected errors for model fitting at the given signal-to-noise ratio. With the uv-coverage of our observations (18-160 m), the size and flux density measurements are sensitive to scales out to 2 . We compare this to a simulated ALMA Cycle 3 dataset with intermediate length baseline coverage, and we find that, using only these baselines, the measured stacked flux density would be an order of magnitude fainter. This highlights the importance of short baselines to recover the full flux density of high-redshift galaxies.

We present a study of the complete ultraviolet to submillimetre spectral energy distributions (SEDs) of twelve 3CR radio galaxy hosts in the redshift range $1.0 < z < 2.5$, which were all detected in the far-infrared by the Herschel Space Observatory. The study employs the new spectro-chemical evolutionary code P\'EGASE.3, in combination with recently published clumpy AGN torus models. We uncover the properties of the massive host galaxy stellar populations, the AGN torus luminosities, and the properties of the recent starbursts, which had earlier been inferred in these objects from their infrared SEDs. The P\'EGASE.3 fitting yields very luminous (up to 10$^{13}$L$_{\odot}$) young stellar populations with ages of several hundred million years in hosts with masses exceeding 10$^{11}$M$_{\odot}$. Dust masses are seen to increase with redshift, and a surprising correlation -- or better upper envelope behaviour -- is found between the AGN torus luminosity and the starburst luminosity, as revealed by their associated dust components. The latter consistently exceeds the former by a constant factor, over a range of one order of magnitude in both quantities.

High-redshift radio galaxies present signs of both star formation and AGN activity, making them ideal candidates to investigate the connection and coevolution of AGN and star formation in the progenitors of present-day massive galaxies. We make use of a sample of 11 powerful radio galaxies spanning 1<z<4 which have complete coverage of their spectral energy distribution (SED) from UV to FIR wavelengths. Using Herschel data, we disentangle the relative contribution of the AGN and star formation by combining the galaxyevolutioncodePEGASE.3 with an AGN torus model. We find that three components are necessary to reproduce the observed SEDs: an evolved and massive stellar component, a submm bright young starburst, and an AGN torus. We find that powerful radio galaxies form at very high-redshift, but experience episodic and important growth at 1<z<4 as the mass of the associated starburst varies from 5 to 50% of the total mass of the system. The properties of star formation differ from source to source, indicating no general trend of the star formation properties in the most infrared luminous high-redshift radio galaxies and no correlation with the AGN bolometric luminosity. Moreover, we find that AGN scattered light have a very limited impact on broad-band SED fitting on our sample. Finally, our analysis also suggests a wide range in origins for the observed star formation,which we partially constrain for some sources.

The far-infrared fine-structure line [CII] at 1900.5\,GHz is known to be one of the brightest cooling lines in local galaxies, and therefore it has been suggested to be an efficient tracer for star-formation in very high-redshift galaxies. However, recent results for galaxies at $z>6$ have yielded numerous non-detections in star-forming galaxies, except for quasars and submillimeter galaxies. We report the results of ALMA observations of two lensed, star-forming galaxies at $z = 6.029$ and $z=6.703$. The galaxy A383-5.1 (star formation rate [SFR] of 3.2 M$_\odot$ yr$^{-1}$ and magnification of $\mu = 11.4\pm1.9$) shows a line detection with $L_{\rm [CII]} = 8.9\times10^{6}$ L$_\odot$, making it the lowest $L_{\rm [CII]}$ detection at $z>6$. For MS0451-H (SFR = 0.4 M$_\odot$ yr$^{-1}$ and $\mu = 100\pm20$) we provide an upper limit of $L_{\rm [CII]} < 3\times10^{5}$ L$_\odot$, which is 1\,dex below the local SFR-$L_{\rm [CII]}$ relations. The results are consistent with predictions for low-metallicity galaxies at $z>6$, however, other effects could also play a role in terms of decreasing $L_{\rm [CII]}$. The detection of A383-5.1 is encouraging and suggests that detections are possible, but much fainter than initially predicted.

We utilize a Bayesian approach to fit the observed mid-IR-to-submm/mm spectral energy distributions (SEDs) of 22 WISE-selected and submm-detected, hyperluminous hot dust-obscured galaxies. By adopting the Torus+GB model, we decompose the observed IR SEDs of Hot DOGs into torus and cold dust components. The main results are: 1) Hot DOGs in our submm-detected sample are hyperluminous, with torus emission dominating the IR energy output. However, cold dust emission is non-negligible, averagely contributing ~24% of total IR luminosity. 2) Compared to QSO and starburst SED templates, the median SED of Hot DOGs shows the highest luminosity ratio between mid-IR and submm at rest-frame, while it is very similar to that of QSOs at 10-50um suggesting that the heating sources of Hot DOGs should be buried AGNs. 3) Hot DOGs have both high dust temperatures ~73K and IR luminosity of cold dust. The T-L relation of Hot DOGs suggests that the increase in IR luminosity for Hot DOGs is mostly due to the increase of the dust temperature, rather than dust mass. Hot DOGs have lower dust masses than those of submillimeter galaxies (SMGs) and QSOs within the similar redshift range. Both high IR luminosity of cold dust and relatively low dust mass in Hot DOGs can be expected by their relatively high dust temperatures. 4) Hot DOGs have high dust covering factors, which deviate the previously proposed trend of the dust covering factor decreasing with increasing bolometric luminosity. Finally, we can reproduce the observed properties in Hot DOGs by employing a physical model of galaxy evolution. The result suggests that Hot DOGs may lie at or close to peaks of both star formation and black hole growth histories, and represent a transit phase during the evolution of massive galaxies, transforming from the dusty starburst dominated phase to the optically bright QSO phase. (abridged)

We present a pilot study of the z=2.923 radio galaxy MRC0943-242, where we for the first time combine information from ALMA and MUSE data cubes. Even with modest integration times, we disentangle an AGN and a starburst dominated set of components. These data reveal a highly complex morphology, as the AGN, starburst, and molecular gas components show up as widely separated sources in dust continuum, optical continuum and CO line emission observations. CO(1-0) and CO(8-7) line emission suggest that there is a molecular gas reservoir offset from both the dust and the optical continuum that is located ~90kpc from the AGN. The UV line emission has a complex structure in emission and absorption. The line emission is mostly due to i) a large scale ionisation cone energised by the AGN, ii) a Ly-alpha emitting bridge of gas between the radio galaxy and a heavily star-forming set of components. Strangely, the ionisation cone has no Ly-alpha emission. We find this is due to an optically thick layer of neutral gas with unity covering fraction spread out over a region of at least ~100kpc from the AGN. Other, less thick absorption components are associated with Ly-alpha emitting gas within a few tens of kpc from the radio galaxy and are connected by a bridge of emission. We speculate that this linear structure of dust, Ly-alpha and CO emission, and the redshifted absorption seen in the circum-nuclear region may represent an accretion flow feeding gas into this massive AGN host galaxy.

The Dragonfly Galaxy (MRC0152-209), at redshift z~2, is one of the most vigorously star-forming radio galaxies in the Universe. What triggered its activity? We present ALMA Cycle 2 observations of cold molecular CO(6-5) gas and dust, which reveal that this is likely a gas-rich triple merger. It consists of a close double nucleus (separation ~4 kpc) and a weak CO-emitter at ~10 kpc distance, all of which have counterparts in HST/NICMOS imagery. The hyper-luminous starburst and powerful radio-AGN were triggered at this precoalescent stage of the merger. The CO(6-5) traces dense molecular gas in the central region, and complements existing CO(1-0) data, which revealed more widespread tidal debris of cold gas. We also find ~10$^{10}$ M(sun) of molecular gas with enhanced excitation at the highest velocities. At least 20-50% of this high-excitation, high-velocity gas shows kinematics that suggests it is being displaced and redistributed within the merger, although with line-of-sight velocities of |v| < 500 km/s, this gas will probably not escape the system. The processes that drive the redistribution of cold gas are likely related to either the gravitational interaction between two kpc-scale discs, or starburst/AGN-driven outflows. We estimate that the rate at which the molecular gas is redistributed is at least ~1200 +- 500 M(sun)/yr, and could perhaps even approach the star formation rate of ~3000 +- 800 M(sun)/yr. The fact that the gas depletion and gas redistribution timescales are similar implies that dynamical processes can be important in the evolution of massive high-z galaxies.

In the low-redshift Universe, the most powerful radio sources are often associated with gas-rich galaxy mergers or interactions. We here present evidence for an advanced, gas-rich (`wet') merger associated with a powerful radio galaxy at a redshift of z~2. This radio galaxy, MRC 0152-209, is the most infrared-luminous high-redshift radio galaxy known in the southern hemisphere. Using the Australia Telescope Compact Array, we obtained high-resolution CO(1-0) data of cold molecular gas, which we complement with HST/WFPC2 imaging and WHT long-slit spectroscopy. We find that, while roughly M(H2) ~ 2 x 10$^{10}$ M$_{\odot}$ of molecular gas coincides with the central host galaxy, another M(H2) ~ 3 x 10$^{10}$ M$_{\odot}$ is spread across a total extent of ~60 kpc. Most of this widespread CO(1-0) appears to follow prominent tidal features visible in the rest-frame near-UV HST/WFPC2 imaging. Ly$\alpha$ emission shows an excess over HeII, but a deficiency over L(IR), which is likely the result of photo-ionisation by enhanced but very obscured star formation that was triggered by the merger. In terms of feedback, the radio source is aligned with widespread CO(1-0) emission, which suggests that there is a physical link between the propagating radio jets and the presence of cold molecular gas on scales of the galaxy's halo. Its optical appearance, combined with the transformational stage at which we witness the evolution of MRC 0152-209, leads us to adopt the name `Dragonfly Galaxy'.

In this paper, we show that the supernova remnant (SNR) masses cumulated from core-collapse supernovae along the star formation history of two powerful z=3.8 radio galaxies 4C41.17 and TN J2007-1316 reach up to > 10^9 Msun, comparable with supermassive black hole (SMBH) masses measured from the SDSS sample at similar redshifts. The SNR mass is measured from the already exploded supernova mass after subtraction of ejecta at the galaxy age where the mass of still luminous stars fits at best the observed spectral energy distribution (SED), continuously extended to the optical-Spitzer-Herschel-submm domains, with the help of the galaxy evolution model P\'egase.3. For the recent and old stellar populations, SNR masses vary on 10^(9 to 10) Msun and the SNR-to-star mass ratio between 1 and 0.1 percent is comparable to the observed low-z SMBH-to-star mass ratio. For the template radio galaxy 4C41.17, SNR and stellar population masses estimated from large aperture (>4arcsec=30kpc) observations are compatible, within one mass order, with the total mass of multiple optical HST (~700pc) structures, associated with VLA radio emissions, both at 0.1 arcsec. Probing the SNR accretion by central black holes is a simple explanation for SMBH growth, requiring physics on star formation, stellar and galaxy dynamics with consequences on various processes (quenching, mergers, negative feedback) and a key to the relation bulge-SMBH.

We present results from a survey of 70 radio galaxies (RGs) at redshifts 1<z<5.2 using the PACS and SPIRE on-board Herschel. Combined with existing mid-IR photometry from Spitzer and observations obtained with LABOCA, the SEDs of galaxies in our sample are continuously covered across 3.6-870um. The total infrared luminosities of these RGs are such that they almost all are either ultra-or hyper-luminous infrared galaxies. We fit the infrared SEDs with a set of empirical templates which represent dust heated (1) by a variety of SB and (2) by a AGN. We find that the SEDs of RGs require the dust to be heated by both AGN and SB, but the luminosities of these two components are not strongly correlated. Assuming empirical relations and simple physical assumptions, we calculate the SFR, the black hole mass accretion rate (MdotBH), and the black hole mass (MBH) for each RG. We find that the host galaxies and their BHs are growing extremely rapidly, having SFR~100-5000 Msun/yr and MdotBH~1-100 Msun/yr. The mean sSFR of RGs at z>2.5 are higher than the sSFR of typical star-forming galaxies over the same redshift range but are similar or perhaps lower than the galaxy population for RGs at z<2.5. By comparing the sSFR and the specific black hole mass accretion rate, we conclude that BHs in radio loud AGN are already, or soon will be, overly massive compared to their host galaxies in terms of expectations from the local MBH-MGal relation. In order to ``catch up'' with the BH, the galaxies require about an order-of magnitude more time to grow in mass, at the observed SFRs, compared to the time the BH is actively accreting. However, during the current cycle of activity, we argue that this catching-up is likely to be difficult due to the short gas depletion times. Finally, we speculate on how the host galaxies might grow sufficiently in stellar mass to ultimately fall onto the local MBH-MGal relation.

We present new results of CO(1-0) spectroscopic observations of 4 SDSS type 2 quasars (QSO2) at z$\sim$0.3, observed with the 30m IRAM telescope. The QSO2 have infrared luminosities in the ULIRG (UltraLuminous Infrared Galaxies) regime. We confirm the CO(1-0) detection in one of our 4 QSO2, SDSS J1543-00, with $L'_{CO}$ and $M_{H_2}$ (1.2$\pm$0.2) $\times$10$^{10}$ K km s$^{-1}$ pc$^2$ and (9.4$\pm$1.4)$\times$10$^9$ M$_{\odot}$, respectively. The CO(1-0) line has $FWHM=$575$\pm$102 km s$^{-1}$. No CO(1-0) emission is detected in SDSS J0903+02, SDSS J1337-01, SDSS J1520-01 above 3 sigma, yielding upper limits on $M(H_2)\sim$ 9.6, 4.3 and 5.1 $\times$10$^9$ M$_{\odot}$ respectively. Together with CO measurements of 9 QSO2 at $z\sim$0.3-1.0 from the ULIRG sample by Combes et al. (2011, 2013), we expand previous studies of the molecular gas content of intermediate $z$ QSO2 into the ULIRG regime. We discuss the location of the 13 ULIRG QSO2 at $z<$1 with available $L'_{CO}$ measurements in the $L'_{CO}$ vs. $z$ and $L'_{CO}$ vs. $L_{FIR}$ diagrams, in comparison with other QSO1 and ULIRG star forming samples.

We present a CO(1-0) survey for cold molecular gas in a representative sample of 13 high-z radio galaxies (HzRGs) at 1.4<z<2.8, using the Australia Telescope Compact Array. We detect CO(1-0) emission associated with five sources: MRC 0114-211, MRC 0152-209, MRC 0156-252, MRC 1138-262 and MRC 2048-272. The CO(1-0) luminosities are in the range $L'_{\rm CO} \sim (5 - 9) \times 10^{10}$ K km/s pc$^{2}$. For MRC 0152-209 and MRC 1138-262 part of the CO(1-0) emission coincides with the radio galaxy, while part is spread on scales of tens of kpc and likely associated with galaxy mergers. The molecular gas mass derived for these two systems is M$_{\rm H2} \sim 6 \times 10^{10}\, {\rm M}_{\odot}$ (M$_{\rm H2}$/$L'_{\rm CO}$=0.8). For the remaining three CO-detected sources, the CO(1-0) emission is located in the halo (~50-kpc) environment. These three HzRGs are among the fainter far-IR emitters in our sample, suggesting that similar reservoirs of cold molecular halo gas may have been missed in earlier studies due to pre-selection of IR-bright sources. In all three cases the CO(1-0) is aligned along the radio axis and found beyond the brightest radio hot-spot, in a region devoid of 4.5$\mu$m emission in Spitzer imaging. The CO(1-0) profiles are broad, with velocity widths of ~ 1000 - 3600 km/s. We discuss several possible scenarios to explain these halo reservoirs of CO(1-0). Following these results, we complement our CO(1-0) study with detections of extended CO from the literature and find at marginal statistical significance (95% level) that CO in HzRGs is preferentially aligned towards the radio jet axis. For the eight sources in which we do not detect CO(1-0), we set realistic upper limits of $L'_{\rm CO} \sim 3-4 \times 10^{10}$ K km/s pc$^{2}$. Our survey reveals a CO(1-0) detection rate of 38%, allowing us to compare the CO(1-0) content of HzRGs with that of other types of high-z galaxies.

We present results of CO(1-0) spectroscopic observations of 10 SDSS type 2 quasars (QSO2) at z~0.2-0.3 observed with the 30m IRAM radiotelescope and the Australia Telescope Compact Array. We report 5 new confirmed CO(1-0) detections and 1 tentative detection. They have L'(CO)~several x 1e9 K km s-1 pc^2, while upper limits for the non detections are L'(CO)<~ several x 1e9 K km s-1 pc^2. This study increases the total number of QSO2 with CO measurements at z<~1 to 20, with a 50% detection rate.

We have characterized the amount, spatial distribution and kinematics of the molecular gas in the merging, double nucleus type 2 quasar SDSS J0025-10 at z=0.30 using the CO(1-0) transition, based on data obtained with ATCA. This is one of the scarce examples of quasar host galaxies where the CO emission has been resolved spatially at any redshift. We infer a molecular gas mass M(H2) = (6 +/- 1) x 1e9 Msun, which is distributed in two main reservoirs separated by ~9 kpc. ~60% of the gas is in the central region, associated with the QSO nucleus and/or the intermediate region between the two nuclei. The other 40% is associated with the northern tidal tail and is therefore unsettled. With its high infrared luminosity L(IR) = (1.1 +/- 0.3) x 1e12 Lsun, SDSS J0025-10 is an analogue of local luminous LIRGs and ULIRGs. On the other hand, the clear evidence for an ongoing major merger of two gas rich progenitors, the high L(IR) dominated by a starburst, the massive reservoir of molecular gas with a large fraction still unsettled, and the quasar activity are all properties consistent with a transition phase in the (U)LIRG-optical QSO evolutionary scenario. We propose that we are observing the system during a particular transient phase, prior to more advanced mergers where the nuclei have already coalesced. We argue that a fraction of the molecular gas reservoir is associated with a tidal dwarf galaxy identified in the optical HST image at the tip of the northern tidal tail. The formation of such structures is predicted by simulations of colliding galaxies.

The new evolutionary code Pegase.3 is used for the spectral synthesis of continuous optical-Spitzer-Herschel-submm spectral energy distributions (SEDs) of two z=3.8 radio galaxies 4C 41.17 and TN J2007-1316. Both galaxies are selected from the Herschel Radio Galaxy Evolution Project for their faint AGN contribution and because they show evidence of a large stellar contribution to their bolometric luminosity. Pegase.3 coherently models the reprocessing of the stellar luminosity to dust emission, allowing to build UV to IR-submm libraries by types. These templates allow to predict SEDs at redshifts z in the observer's frame by using z=0 templates corrected for cosmology and evolution. Our best fits are a sum of two distinct evolving stellar populations: i) a very massive starburst observed ~30 Myrs after formation from optical and far-IR peaks ii) a significantly older (~1Gyr old) massive early-type population revealed from mid-IR Spitzer data. The AGN emission fits a roughly thermal model. These results confirm that many of the stellar populations in high-z radio galaxies were already formed by massive (a few thousands of billion solar masses) starbursts at z=4, confirming the previous K-z diagram interpretation. Gas-rich mergers and/or jet-cloud interactions are favored for triggering the intense star formation, possibly in relation with massive super black holes of these radio galaxies. These multiple stellar populations may be a generic feature of the luminous infrared radio galaxy population as a whole and strongly constrain mass accumulation and galaxy formation models.

We present Herschel observations at 70, 160, 250, 350 and 500 micron of the environment of the radio galaxy 4C+41.17 at z = 3.792. About 65% of the extracted sources are securely identified with mid-IR sources observed with the Spitzer Space Telescope at 3.6, 4.5, 5.8, 8 and 24 micron. We derive simple photometric redshifts, also including existing 850 micron and 1200 micron data, using templates of AGN, starburst-dominated systems and evolved stellar populations. We find that most of the Herschel sources are foreground to the radio galaxy and therefore do not belong to a structure associated with 4C+41.17. We do, however, find that the SED of the closest (~ 25" offset) source to the radio galaxy is fully consistent with being at the same redshift as 4C+41.17. We show that finding such a bright source that close to the radio galaxy at the same redshift is a very unlikely event, making the environment of 4C+41.17 a special case. We demonstrate that multi-wavelength data, in particular on the Rayleigh-Jeans side of the spectral energy distribution, allow us to confirm or rule out the presence of protocluster candidates that were previously selected by single wavelength data sets.

We examine the relative orientation of radio jets and dusty tori surrounding the AGN in powerful radio galaxies at z > 1. The radio core dominance R = P(20 GHz) /P(500 MHz) serves as an orientation indicator, measuring the ratio between the anisotropic Doppler-beamed core extended core emission and the isotropic lobe emission. Assuming a fixed cylindrical geometry for the hot, dusty torus, we derive its inclination i by fitting optically-thick radiative transfer models to spectral energy distributions obtained with the Spitzer Space Telescope. We find a highly significant anti-correlation (p < 0.0001) between R and i in our sample of 35 type 2 AGN combined with a sample of 18 z = 1 3CR sources containing both type 1 and 2 AGN. This analysis provides observational evidence both for the Unified scheme of AGN and for the common assumption that radio jets are in general perpendicular to the plane of the torus. The use of inclinations derived from mid-infrared photometry breaks several degeneracies which have been problematic in earlier analyses. We illustrate this by deriving the core Lorentz factor Gamma from the R-i anti-correlation, finding Gamma > 1.3.

We present a detailed study of the infrared spectral energy distribution of the high-redshift radio galaxy MRC 1138-26 at z = 2.156, also known as the Spiderweb Galaxy. By combining photometry from Spitzer, Herschel and LABOCA we fit the rest-frame 5-300 um emission using a two component, starburst and active galactic nucleus (AGN), model. The total infrared (8 - 1000 um) luminosity of this galaxy is (1.97+/-0.28)x10^13 Lsun with (1.17+/-0.27) and (0.79+/-0.09)x10^13 Lsun due to the AGN and starburst components respectively. The high derived AGN accretion rate of \sim20% Eddington, and the measured star formation rate (SFR) of 1390pm150 Msun/yr, suggest that this massive system is in a special phase of rapid central black hole and host galaxy growth, likely caused by a gas rich merger in a dense environment. The accretion rate is sufficient to power both the jets and the previously observed large outflow. The high SFR and strong outflow suggest this galaxy could potentially exhaust its fuel for stellar growth in a few tens of Myr, although the likely merger of the radio galaxy with nearby satellites suggest bursts of star formation may recur again on time scales of several hundreds of Myr. The age of the radio lobes implies the jet started after the current burst of star formation, and therefore we are possibly witnessing the transition from a merger-induced starburst phase to a radio-loud AGN phase. We also note tentative evidence for [CII]158um emission. This paper marks the first results from the Herschel Galaxy Evolution Project (Project HeRGE), a systematic study of the evolutionary state of 71 high redshift, 1 < z < 5.2, radio galaxies.