Search

Next-generation radio surveys are expected to detect tens of millions of active galactic nuclei (AGNs) with a median redshift of $z$ ≥ 1. Beyond targeted surveys, the vast majority of these objects will not have spectroscopic redshifts, while photometric redshifts for high-redshift AGNs are of limited quality, and even then require optical and infrared photometry. We propose a new approach to measure the redshifts of lobed radio galaxies based exclusively on radio-frequency imaging and broad-band radio photometry. Specifically, our algorithm uses the lobe flux density, angular size and width, and spectral shape to derive probability density functions for the most likely source redshift based on the Radio AGN in Semi-analytic Environments dynamical model. The full physically based model explains 70 per cent of the variation in the spectroscopic redshifts of a high-redshift (2 < $z$ < 4) sample of radio AGNs, compared to at most 27 per cent for any one of the observed attributes in isolation. We find that upper bounds on the angular size, as expected for unresolved sources, are sufficient to yield accurate redshift measurements at $z$ ≥ 2. The error in the model upon calibration using at least nine sources with known spectroscopic redshifts is

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., Comment: Accepted by MNRAS; 14 pages, 9 figures

The recent detection of the quasi-stellar object (QSO) VIKING J231818.3$-$311346 (hereafter VIK J2318$-$3113) at redshift $z=6.44$ in the Rapid ASKAP Continuum Survey (RACS) uncovered its radio-loud nature, making it one of the most distant known to date in this class. By using data from several radio surveys of the Galaxy and Mass Assembly 23$^\mathrm{h}$ field and from dedicated follow-up, we were able to constrain the radio spectrum of VIK J2318$-$3113 in the observed range $\sim$0.1--10 GHz. At high frequencies (0.888--5.5 GHz in the observed frame) the QSO presents a steep spectrum ($\alpha_{\rm r}$=1.24, with $S_\nu\propto \nu^{-\alpha_{\rm r}}$), while at lower frequencies (0.4--0.888 GHz in the observed frame) it is nearly flat. The overall spectrum can be modelled by either a curved function with a rest-frame turnover around 5 GHz, or with a smoothly varying double power law that is flat below a rest-frame break frequency of about 20 GHz and which significantly steepens above it. Based on the model adopted, we estimated that the radio jets of VIK J2318$-$3113 must be a few hundred years old, in the case of a turnover, or less than few$\times$10$^4$ years, in the case of a break in the spectrum. Having multiple observations at two frequencies (888 MHz and 5.5 GHz), we further investigated the radio variability previously reported for this source. We found that the marginally significant flux density variations are consistent with the expectations from refractive interstellar scintillation, even though relativistic effects related to the orientation of the source may still play a non-negligible role. Further radio and X-ray observations are required to conclusively discern the nature of this variation., Comment: Accepted for publication in A&A

We report the 888 MHz radio detection in the Rapid ASKAP Continuum Survey (RACS) of VIK J2318$-$3113, a $z$=6.44 quasar. Its radio luminosity (1.2 $\times 10^{26}$ W Hz$^{-1}$ at 5 GHz) compared to the optical luminosity (1.8 $\times 10^{24}$ W Hz$^{-1}$ at 4400 A) makes it the most distant radio-loud quasar observed so far, with a radio loudness R$\sim$70 (R$=L_\mathrm{{5GHz}}/L_\mathrm{{4400A}}$). Moreover, the high bolometric luminosity of the source (L$_\mathrm{{bol}}$=7.4 $\times 10^{46}$ erg s$^{-1}$) suggests the presence of a supermassive black hole with a high mass ($\gtrsim$6 $\times$10$^8$ M$_\odot$) at a time when the Universe was younger than a billion years. Combining the new radio data from RACS with previous ASKAP observations at the same frequency, we found that the flux density of the source may have varied by a factor of $\sim$2, which could suggest the presence of a relativistic jet oriented towards the line of sight, that is, a blazar nature. However, currently available radio data do not allow us to firmly characterise the orientation of the source. Further radio and X-ray observations are needed., Comment: 5 pages, 2 figures, 3 tables. Accepted for publication in A&A Letters

We report a detailed CO(1-0) survey of a galaxy protocluster field at $z=2.16$, based on 475 hours of observations with the Australia Telescope Compact Array. We constructed a large mosaic of 13 individual pointings, covering an area of 21 arcmin$^2$ and $\pm6500$ km/s range in velocity. We obtain a robust sample of 46 CO(1-0) detections spanning $z=2.09-2.22$, constituting the largest sample of molecular gas measurements in protoclusters to date. The CO emitters show an overdensity at $z=2.12-2.21$, suggesting a galaxy super-protocluster or a protocluster connected to large-scale filaments with ~120 cMpc size. We find that 90% CO emitters have distances $>0'.5-4'$ to the center galaxy, indicating that small area surveys would miss the majority of gas reservoirs in similar structures. Half of the CO emitters have velocities larger than escape velocities, which appears gravitationally unbound to the cluster core. These unbound sources are barely found within the $R_{200}$ radius around the center, which is consistent with a picture in which the cluster core is collapsed while outer regions are still in formation. Compared to other protoclusters, this structure contains relatively more CO emitters with relatively narrow line width and high luminosity, indicating galaxy mergers. We use these CO emitters to place the first constraint on the CO luminosity function and molecular gas density in an overdense environment. The amplitude of the CO luminosity function is 1.6$\pm$0.5 orders of magnitudes higher than observed for field galaxy samples at $z\sim2$, and one order of magnitude higher than predictions for galaxy protoclusters from semi-analytical SHARK models. We derive a high molecular gas density of $0.6-1.3\times10^{9}$ $M_\odot$ cMpc$^{-3}$ for this structure, consistent with predictions for cold gas density of massive structures from hydro-dynamical DIANOGA simulations., 21 pages + Appendices, 9 figures, 6 tables. Accepted for publication in A&A

We previously reported a putative detection of a radio galaxy at z=10.15, selected from the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey. The redshift of this source, GLEAM J0917-0012, was based upon three weakly detected molecular emission lines observed with the Atacama Large Millimetre Array (ALMA). In order to confirm this result, we conducted deep spectroscopic follow-up observations with ALMA and the Karl Jansky Very Large Array (VLA). The ALMA observations targeted the same CO lines previously reported in Band 3 (84-115GHz) and the VLA targeted the CO(4-3) and [CI(1-0)] lines for an independent confirmation in Q-band (41 and 44GHz). Neither observation detected any emission lines, removing support for our original interpretation. Adding publicly available optical data from the Hyper Suprime-Cam survey, WISE and Herschel Space Observatory in the infrared, as well as, 19 pages, 6 figures, published in PASA

We present the results of a new selection technique to identify powerful ($L_{\rm 500\,MHz}>10^{27}\,$WHz$^{-1}$) radio galaxies towards the end of the Epoch of Reionisation. Our method is based on the selection of bright radio sources showing radio spectral curvature at the lowest frequency ($\sim 100\,$MHz) combined with the traditional faintness in $K-$band for high redshift galaxies. This technique is only possible thanks to the Galactic and Extra-galactic All-sky Murchison wide-field Array (GLEAM) survey which provides us with 20 flux measurements across the $70-230\,$MHz range. For this pilot project, we focus on the GAMA 09 field to demonstrate our technique. We present the results of our follow-up campaign with the Very Large Telescope, Australian Telescope Compact Array and the Atacama Large Millimetre Array (ALMA) to locate the host galaxy and to determine its redshift. Of our four candidate high redshift sources, we find two powerful radio galaxies in the $15.5$ has a high $25-50\%$ success rate., Comment: 19 pages, 9 figures, published in PASA

Both theoretical and observational results suggest that high-redshift radio galaxies (HzRGs) inhabit overdense regions of the universe and might be the progenitors of local, massive galaxies residing in the centre of galaxy clusters. In this paper we present CO(3-2) line observations of the HzRG TXS 0828+193 (z=2.57) and its environment using the Atacama Large Millimeter/submillimeter Array. In contrast to previous observations, we detect CO emission associated with the HzRG and derive a molecular gas mass of $(0.9\pm0.3)\times10^{10}\,\rm M_{\odot}$. Moreover, we confirm the presence of a previously detected off-source CO emitting region (companion #1), and detect three new potential companions. The molecular gas mass of each companion is comparable to that of the HzRG. Companion #1 is aligned with the axis of the radio jet and has stellar emission detected by Spitzer. Thus this source might be a normal star-forming galaxy or alternatively a result of jet-induced star formation. The newly found CO sources do not have counterparts in any other observing band and could be high-density clouds in the halo of TXS 0828+193 and thus potentially linked to the large-scale filamentary structure of the cosmic web., Comment: Accepted by MNRAS; 9 pages, 4 figures

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., Comment: 21 pages, 20 figures, accepted for publication in A&A

International audience; We present the public release of MR-MOOSE, 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 a heterogeneous data set (i.e. a range of observation sensitivities and spectral/spatial resolutions). Furthermore, MR-MOOSE 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 for handling increasing source complexity when combining multi-wavelength data sets with fully customisable filter/model data bases. 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 data sets 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 MR-MOOSE 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., Comment: 3 figures, 2 tables

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, Comment: A&A in press. Corner plots for the fits are available from http://www.eso.org/~cdebreuc/sed_paper/

The HeRGÉ (Herschel Radio Galaxy Evolution) project consists of a sample of 70 radio galaxies in the range 1 < z < 5.2. They benefit from continuous coverage from 3 to 870μm with Spitzer, Herschel and sub-mm ground-based instruments (SCUBA, LABOCA). As a calorimeter, IR is an excellent proxy to estimate the contribution of both AGN and starburst, making of radio galaxies perfect candidates to provide new insights into the relationship between AGN and their host galaxies. The IR SED fitting with empirical templates reveals that radio galaxies are luminous and that their black holes and their host galaxies are not growing simultaneously. Extending the SED to optical/near-IR on a subsample of 12 radio galaxies spanning 1 < z < 4 reveal the necessity of three components to reproduce the observations. Making use of the evolutionary code PEGASE.3 and an AGN torus model, we are able to estimate parameters from the AGN torus, the evolved stellar population and the starburst (SB). They reveal that radio galaxies are massive, evolved, forming the bulk of their mass at very high redshift in a short timescale, but experience episodic, strong SB events, often associated with an AGN activity.

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.

Using the new evolutionary code PEGASE. 3, we undertook an evolutionary spectral synthesis of the optical-IR-submm spectral energy distribution of two distant (z = 3.8) radio galaxies, 4C 41.17 and TN J2007-1316. These two radio galaxies were selected from the HeRGE (Herschel Radio Galaxies Evolution) Project in particular for their faint active galactic nucleus 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 us to build UV to IR-submm spectral energy distribution libraries that can then be used to fit spectral energy distributions in the observer's frame. Our principal conclusion is that a single stellar population is insufficient to fit the spectral energy distribution of either radio galaxy. Our best fits are a sum of two evolving stellar populations - a recent starburst plus an old population - plus the thermal emission from an active galactic nucleus (which provides a good fit to the mid-IR emission). The two stellar components are: (i) a massive (similar or equal to 10(11) M-circle dot) starburst similar or equal to 30 Myr after formation, which is required simultaneously to fit the far-IR Herschel to submm data and the optical data; and (ii) an older massive (similar or equal to 10(11-12) M-circle dot) early-type galaxy population, similar or equal to 1.0 Gyr old, which is required principally to fit the mid-IR Spitzer/IRAC data. A young population alone is insufficient because an evolved giant star population produces a 1-mu m rest-frame peak that is observed in the IRAC photometry. This discovery confirms that many of the stellar populations in high-redshift radio galaxies were formed by massive starbursts in the early Universe. Gas-rich mergers and/or jet-cloud interactions are favoured for triggering the intense star formation necessary to explain the properties of the spectral energy distributions. The discovery of similar characteristics in two distant radio galaxies suggests that multiple stellar populations, one old and one young, may be a generic feature of the luminous infrared radio galaxy population.

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., Accepted for publication in MNRAS

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., 11 pages, 2 figures, accepted for publication in ApJ

(Abridged) In the last decade several massive molecular gas reservoirs were found 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., 10 pages, 5 figures, accepted for publication in A&A

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), Comment: 16 pages, 12 figures. Accepted by ApJ

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��GASE.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��GASE.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., 15 pages, 9 figures, accepted for publication in MNRAS

Distant radio galaxies, hosted by massive ellipticals, follow the galaxy evolution process on an extremely large (0 ≥ z ≥7) time-scale ≥ 1012Gyrs, up to primeval galaxies. The new evolutionary code Pégase.3 predicts on similar time-scales, the coupled stellar and dust emissions of various galaxy types: starbursts and Hubble sequence types. All z=0 templates are fitted on local observations at ages ≃13 Gyrs (except irregulars at 9 Gyrs). The multi-λ spectral energy distributions (SEDs) of two z=3.8 radiogalaxies, including the most recent Herschel data from the HeRGÉ consortium, are interpreted in the observer's frame by Rocca-Volmerange et al. (2012) with Pégase.3. The apparent SEDs are fitted at best with the sum of a young starburst and an older early-type population, an AGN simple model is taken into account. These results favor massive gas-rich mergers at work in evolved galaxies at z≃4. Massive starbursts would be at the origin of galaxy evolution initiated at the earliest epochs (zfor≥10). The possible relation with super massive black holes is still debated.

The Dragonfly Galaxy (MRC 0152-209), at redshift z similar to 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 similar to 4 kpc) and a weak CO-emitter at similar to 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 reveal more widespread tidal debris of cold gas. We also find similar to 10(10) M-circle dot 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 vertical bar v vertical bar < 500 km s(-1), 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. M similar to 1200 +/- 500 M-circle dot yr(-1), and could perhaps even approach the star formation rate of similar to 3000 +/- 800 M-circle dot yr(-1). 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'., Comment: Accepted for publication in MNRAS (12 pages, 5 figures)

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., Comment: 15 pages, 8 figures, accepted for publication in A&A

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., A&A accepted. Comments welcome

We present results from a survey of 70 radio galaxies (RGs) at redshifts 12.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, Comment: 37 pages, accepted to A&A

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, Comment: Accepted for publication in MNRAS; 19 pages, 11 figures

We present results of CO(1-0) spectroscopic observations of 10 Sloan Digital Sky Survey type 2 quasars (QSO2) at z ~ 0.2-0.3 observed with the 30-m IRAM radio telescope and the Australia Telescope Compact Array. We report five new confirmed CO(1-0) detections and one tentative detection. They have L'CO ~ several × 109 K km s-1 pc2, while upper limits for the non-detections are L'CO < 3s = several × 109 K km s-1 pc2. This study increases the total number of QSO2 with CO measurements at z 1 to 20, with a 50 per cent detection rate. The vast majority are at z ~ 0.1-0.4. Assuming a conversion factor α = 0.8 M⊙(K km s-1pc2)-1, the implied molecular gas masses are in the range MH2 4 × 108 to ~5 × 109 M⊙. We compare with samples of type 1 quasars (QSO1), luminous and ultraluminous infrared galaxies. We find no difference in the molecular gas content of QSO1 and QSO2 of a given infrared luminosity, although the QSO2 sample is affected by small number statistics. This result, if confirmed, is consistent with the unification model for quasars. QSO2 fall on the L'COversus z, L'COversus LFIR and λ = LFIR LCO versus LFIR correlations defined by quasars at different z. The location of the QSO2 in these diagrams is discussed in comparison with samples of QSO1, luminous and ultraluminous infrared galaxies, and high-z submm sources. CO(1-0) has full width at half-maximum (FWHM) ~ 180-370 km s-1 when detected, with a variety of kinematic profiles (single or double horned). In general, the CO line is narrower than [O III]5007, as observed in low-z QSO1, with FWHM[O III]/FWHMCO ~ 1-2). This probably reveals different spatial sizes and/or geometry of the ionized and molecular phases and a higher sensitivity of the [O III] emission to non-gravitational motions, such as outflows. Considering the z ~ 0.1-0.4 range, where CO measurements for both QSO1 and QSO2 exist, we find no difference in FWHMCO between them, although this result is tentative. In the unification scenario between QSO1 and QSO2, this suggests that the distribution of CO gas is not related to the obscuring torus., Funding for the SDSS and SDSS-II has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, the U.S. Department of Energy, the National Aeronautics and Space Administration, the Japanese Monbukagakusho, the Max Planck Society and the Higher Education Funding Council for England.

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., 17 pages, accepted for publication in the Monthly Notices of the Royal Astronomical Society

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., 18 pages, 12 figures, 6 tables

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��gase.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., 13 pages, 2 double figures, 1 table, accepted for publication in Astrophysical Journal Letters, in press

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, Comment: 4 pages, 2 figures, Accepted for publication in A&A

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 galaxy evolution code PEGASE.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.