Passive early-type galaxies dominate cluster cores at z $\lesssim$1.5. At higher redshift, cluster core galaxies are observed to have still on-going star-formation, fuelled by cold molecular gas. We measure the molecular gas reservoir of the central region around the radio-loud AGN in the cluster CARLA J1103+3449 at z=1.44 with NOEMA. The AGN synchrotron emission dominates the continuum emission at 94.48 GHz, and we measure its flux at the AGN position and at the position of two radio jets. Combining our measurements with published results over the range 4.71 GHz-94.5 GHz, we obtain a flat spectral index $\alpha = 0.14 \pm 0.03$ for the AGN core emission, and a steeper index $\alpha = 1.43 \pm 0.04$ and $\alpha = 1.15 \pm 0.04$ at positions close to the western and eastern lobe, respectively. The total spectral index is $\alpha = 0.92 \pm 0.02$ over the range 73.8 MHz-94.5 GHz. We detect two CO(2-1) emission lines, both blue-shifted with respect to the AGN. Their emission corresponds to two regions, ~17 kpc south-east and ~14 kpc south-west of the AGN, not associated with galaxies. In these two regions, we find a total massive molecular gas reservoir of $M_{gas}$ = 3.9 $\pm$ 0.4 $10^{10} M_{\odot}$, which dominates (~ 60%) the central total molecular gas reservoir. These results can be explained by massive cool gas flows in the center of the cluster. The AGN early-type host is not yet quenched; its star formation rate is consistent with being on the main sequence of star-forming galaxies in the field (SFR~30-140 $M_{\odot}$/yr), and the cluster core molecular gas reservoir is expected to feed the AGN and the host star-formation before quiescence. The other cluster confirmed members show star formation rates at ~2 $\sigma$ below the field main sequence at similar redshifts and do not have molecular gas masses larger than galaxies of similar stellar mass in the field.

We present a multi-wavelength investigation of the radio galaxy population in the galaxy cluster MOO J1506+5137 at $z$=1.09$\pm$0.03, which in previous work we identified as having multiple complex radio sources. The combined dataset used in this work includes data from the Low-Frequency Array Two-metre Sky Survey (LoTSS), NSF's Karl G. Jansky Very Large Array (VLA), the Robert C. Byrd Green Bank Telescope (GBT), the Spitzer Space Telescope, and the Dark Energy Camera Legacy Survey (DECaLS). We find that there are five radio sources which are all located within 500 kpc ($\sim$1$^{\prime}$) of the cluster center and have radio luminosities $P_{\mathrm{1.4GHz}}$ > 1.6$\times$10$^{24}$ W Hz$^{-1}$. The typical host galaxies are among the highest stellar mass galaxies in the cluster. The exceptional radio activity among the massive galaxy population appears to be linked to the dynamical state of the cluster. The galaxy distribution suggests an ongoing merger, with a subgroup found to the northwest of the main cluster. Further, two of the five sources are classified as bent-tail sources with one being a potential wide-angle tail (WAT)/hybrid morphology radio source (HyMoRS) indicating a dynamic environment. The cluster also lies in a region of the mass-richness plane occupied by other merging clusters in the Massive and Distant Clusters of WISE Survey (MaDCoWS). The data suggest that during the merger phase radio activity can be dramatically enhanced, which would contribute to the observed trend of increased radio activity in clusters with increasing redshift.

The properties of galaxy clusters as a function of redshift can be utilized as an important cosmological tool. We present initial results from a program of follow-up observations of the Sunyaev-Zeldovich effect (SZE) in high redshift galaxy clusters detected at infrared wavelengths in the Massive and Distant Clusters of \textit{WISE} Survey (MaDCoWS). Using typical on-source integration times of 3--4 hours per cluster, MUSTANG2 on the Green Bank Telescope was able to measure strong detections of SZE decrements and statistically significant masses on 14 out of 16 targets. On the remaining two, weaker (3.7$\sigma$) detections of the SZE signal and strong upper limits on the masses were obtained. In this paper we present masses and pressure profiles of each target and outline the data analysis used to recover these quantities. Of the clusters with strong detections, three show significantly flatter pressure profiles while, from the MUSTANG2 data, five others show signs of disruption at their cores. However, outside of the cores of the clusters, we were unable to detect significant amounts of asymmetry. Finally, there are indications that the relationship between optical richness used by MaDCoWS and SZE-inferred mass may be significantly flatter than indicated in previous studies.

The Massive and Distant Clusters of WISE Survey (MaDCoWS) provides a catalog of high-redshift ($0.7\lesssim z\lesssim 1.5$) infrared-selected galaxy clusters. However, the verification of the ionized intracluster medium, indicative of a collapsed and nearly virialized system, is made challenging by the high redshifts of the sample members. The main goal of this work is to test the capabilities of the Atacama Compact Array (ACA; also known as the Morita Array) Band 3 observations, centered at about 97.5 GHz, to provide robust validation of cluster detections via the thermal Sunyaev-Zeldovich (SZ) effect. Using a pilot sample that comprises ten MaDCoWS galaxy clusters, accessible to ACA and representative of the median sample richness, we infer the masses of the selected galaxy clusters and respective detection significance by means of a Bayesian analysis of the interferometric data. Our test of the "Verification with the ACA - Localization and Cluster Analysis" (VACA LoCA) program demonstrates that the ACA can robustly confirm the presence of the virialized intracluster medium in galaxy clusters previously identified in full-sky surveys. In particular, we obtain a significant detection of the SZ effect for seven out of the ten VACA LoCA clusters. We note that this result is independent of the assumed pressure profile. However, the limited angular dynamic range of the ACA in Band 3 alone, short observational integration times, and possible contamination from unresolved sources limit the detailed characterization of the cluster properties and the inference of the cluster masses within scales appropriate for the robust calibration of mass-richness scaling relations.

The prevalence of outflow and feedback signatures in active galactic nuclei (AGN) is a major unresolved question which large integral field unit (IFU) surveys now allow to address. In this paper, we present a kinematic analysis of the ionised gas in 2778 galaxies at z $\sim$ 0.05 observed by SDSS-IV MaNGA. Specifically, we measure the kinematics of the [OIII] {\lambda}5007{\AA} emission line in each spatial element and fit multiple Gaussian components to account for possible non-gravitational motions of gas. Comparing the kinematics of the ionised gas between 308 MaNGA-selected AGN that have been previously identified through emission line diagnostics and sources not classified as AGN, we find that while 25% of MaNGA-selected AGN show [OIII] components with emission line widths of $>$ 500 km/s in more than 10% of their spaxels, only 7% of MaNGA non-AGN show a similar signature. Even the AGN that do not show nuclear AGN photoionisation signatures and that were only identified as AGN based on their larger scale photoionisation signatures show similar kinematic characteristics. In addition to obscuration, another possibility is that outflow and mechanical feedback signatures are longer lived than the AGN itself. Our measurements demonstrate that high velocity gas is more prevalent in AGN compared to non-AGN and that outflow and feedback signatures in low-luminosity, low-redshift AGN may so far have been underestimated. We show that higher luminosity MaNGA-selected AGN are able to drive larger scale outflows than lower luminosity AGN. But estimates of the kinetic coupling efficiencies are $<$ 1% and suggest that the feedback signatures probed in this paper are unlikely to have a significant impact on the AGN host galaxies. However, continuous energy injection may still heat a fraction of the cool gas and delay or suppress star formation in individual galaxies even when the AGN is weak.

We present the results from a study with NSF's Karl G. Jansky Very Large Array (VLA) to determine the radio morphologies of extended radio sources and the properties of their host galaxies in 50 massive galaxy clusters at z~1. We find a majority of the radio morphologies to be Fanaroff-Riley (FR) type IIs. By analyzing the infrared counterparts of the radio sources, we find that ~40% of the host galaxies are the candidate brightest cluster galaxy (BCG) and ~83% are consistent with being one of the top six most massive galaxies in the cluster. We investigate the role of environmental factors on the radio-loud AGN population by examining correlations between environmental and radio-galaxy properties. We find that the highest stellar mass hosts ($M_{*} \gtrsim$ 4$\times 10^{11} M_{\odot}$) are confined to the cluster center and host compact jets. There is evidence for an increase in the size of the jets with cluster-centric radius, which may be attributed to the decreased ICM pressure confinement with increasing radius. Besides this correlation, there are no other significant correlations between the properties of the radio-AGN (luminosity, morphology, or size) and environmental properties (cluster richness and location within the cluster). The fact that there are more AGN in the cluster environment than the field at this epoch, combined with the lack of strong correlation between galaxy and environmental properties, argues that the cluster environment fosters radio activity but does not solely drive the evolution of these sources at this redshift.

We present VLT/SINFONI observations of 35 quasars at 2.1 < z < 3.2, the majority of which were selected from the Clusters Around Radio-Loud AGN (CARLA) survey. CARLA quasars have large CIV-based black hole (BH) masses (M(BH) > 10^9 Msun) and powerful radio emission (P(500MHz) > 27.5 W/Hz). We estimate Ha-based M(BH), finding a scatter of 0.35 dex compared to CIV. We evaluate several recipes for correcting CIV-based masses, which reduce the scatter to 0.24 dex. The radio power of the radio-loud quasars is at most weakly correlated with the interconnected quantities Ha-width, L(5100A) and M(BH), suggesting that it is governed by different physical processes. However, we do find a strong inverse correlation between CIV blueshift and radio power linked to higher Eddington ratios and L(5100A). Under standard assumptions, the BH growth time is longer than the cosmic age for many CARLA quasars, suggesting that they must have experienced more efficient growth in the past. If these BHs were growing from seeds since the epoch of reionization, it is possible that they grew at the Eddington limit like the quasars at z ~ 6-7, and then continued to grow at the reduced rates observed until z ~ 2-3. Finally, we study the relation between M(BH) and environment, finding a weak positive correlation between M(BH) and galaxy density measured by CARLA.

We present Hubble Space Telescope 1.4-1.6 micron images of the hosts of ten extremely red quasars (ERQs) and six type 2 quasar candidates at z=2-3. ERQs, whose bolometric luminosities range between 10^47 and 10^48 erg/sec, show spectroscopic signs of powerful ionized winds, whereas type 2 quasar candidates are less luminous and show only mild outflows. After performing careful subtraction of the quasar light, we clearly detect almost all host galaxies. The median rest-frame B-band luminosity of the ERQ hosts in our sample is 10^11.2 L_Sun, or 4 L* at this redshift. Two of the ten hosts of ERQs are in ongoing mergers. The hosts of the type 2 quasar candidates are 0.6 dex less luminous, with 2/6 in likely ongoing mergers. Intriguingly, despite some signs of interaction and presence of low-mass companions, our objects do not show nearly as much major merger activity as do high-redshift radio-loud galaxies and quasars. In the absence of an overt connection to major ongoing gas-rich merger activity, our observations are consistent with a model in which the near-Eddington accretion and strong feedback of ERQs are associated with relatively late stages of mergers resulting in early-type remnants. These results are in some tension with theoretical expectations of galaxy formation models, in which rapid black hole growth occurs within a short time of a major merger. Type 2 quasar candidates are less luminous, so they may instead be powered by internal galactic processes.

We investigate Extremely Red Quasars (ERQs), a remarkable population of heavily-reddened quasars at redshift z$\sim$2 - 3 that might be caught during a short-lived 'blow-out' phase of quasar/galaxy evolution. We perform a near-IR observational campaign using Keck/NIRSPEC, VLT/X-shooter and Gemini/GNIRS to measure rest-frame optical spectra of 28 ERQs with median infrared luminosity $\langle$log L(erg/s)$\rangle$ $\sim$46.2. They exhibit the broadest and most blue-shifted [OIII] $\lambda$4959,5007 emission lines ever reported, with widths (w90) ranging between 2053 and 7227 km/s, and maximum outflow speeds (v98) up to 6702 km/s. ERQs on average have [OIII] outflows velocities about 3 times larger than those of luminosity-matched blue quasar samples. We show that the faster [OIII] outflows in ERQs are strongly correlated with their extreme red colors and not with radio-loudness, larger quasar luminosities, nor higher Eddington ratios. We estimate for these objects that at least 3-5 per cent of their bolometric luminosity is being converted into the kinetic power of the observed wind. Our results reveal that ERQs have the potential to strongly affect the evolution of host galaxies.

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 measurements of the stellar mass fractions ($f_\star$) for a sample of high-redshift ($0.93 \le z \le 1.32$) infrared-selected galaxy clusters from the Massive and Distant Clusters of WISE Survey (MaDCoWS) and compare them to the stellar mass fractions of Sunyaev-Zel'dovich (SZ) effect-selected clusters in a similar mass and redshift range from the South Pole Telescope (SPT)-SZ Survey. We do not find a significant difference in mean $f_\star$ between the two selection methods, though we do find an unexpectedly large range in $f_\star$ for the SZ-selected clusters. In addition, we measure the luminosity function of the MaDCoWS clusters and find $m^*= 19.41\pm0.07$, similar to other studies of clusters at or near our redshift range. Finally, we present SZ detections and masses for seven MaDCoWS clusters and new spectroscopic redshifts for five MaDCoWS clusters. One of these new clusters, MOO J1521+0452 at $z=1.31$, is the most distant MaDCoWS cluster confirmed to date.

We investigate the effects of Active Galactic Nuclei (AGN) on the gas kinematics of their host galaxies, using MaNGA data for a sample of 62 AGN hosts and 109 control galaxies (inactive galaxies). We compare orientation of the line of nodes (kinematic Position Angle - PA) measured from the gas and stellar velocity fields for the two samples. We found that AGN hosts and control galaxies display similar kinematic PA offsets between gas and stars. However, we note that AGN have larger fractional velocity dispersion $\sigma$ differences between gas and stars [$\sigma_{frac}=(\sigma_{\rm gas}-\sigma_{stars})/\sigma_{\rm stars}$] when compared to their controls, as obtained from the velocity dispersion values of the central (nuclear) pixel (2.5" diameter). The AGN have a median value of $\sigma_{\rm frac}$ of $<\sigma_{frac}>_{\rm AGN}=0.04$, while the the median value for the control galaxies is $<\sigma_{frac}>_{\rm CTR}=-0.23$. 75% of the AGN show $\sigma_{frac}>-0.13$, while 75% of the normal galaxies show $\sigma_{\rm frac}<-0.04$, thus we suggest that the parameter $\sigma_{\rm frac}$ can be used as an indicative of AGN activity. We find a correlation between the [OIII]$\lambda$5007 luminosity and $\sigma_{frac}$ for our sample. Our main conclusion is that the AGN already observed with MaNGA are not powerful enough to produce important outflows at galactic scales, but at 1-2 kpc scales, AGN feedback signatures are always present on their host galaxies.

We present the Massive and Distant Clusters of WISE Survey (MaDCoWS), a search for galaxy clusters at 0.7<z<1.5 based upon data from the Wide-field Infrared Survey Explorer (WISE) mission. MaDCoWS is the first cluster survey capable of discovering massive clusters at these redshifts over the full extragalactic sky. The search is divided into two regions -- the region of the extragalactic sky covered by Pan-STARRS (Dec>-30 degrees) and the remainder of the southern extragalactic sky at Dec<-30 degrees for which shallower optical data from SuperCOSMOS Sky Survey are available. In this paper we describe the search algorithm, characterize the sample, and present the first MaDCoWS data release -- catalogs of the 2433 highest amplitude detections in the WISE--Pan-STARRS region and the 250 highest amplitude detections in the WISE--SuperCOSMOS region. A total of 1723 of the detections from the WISE--Pan-STARRS sample have also been observed with the Spitzer Space Telescope, providing photometric redshifts and richnesses, and an additional 64 detections within the WISE--SuperCOSMOS region also have photometric redshifts and richnesses. Spectroscopic redshifts for 38 MaDCoWS clusters with IRAC photometry demonstrate that the photometric redshifts have an uncertainty of $\sigma_z/(1+z)\sim0.036$. Combining the richness measurements with Sunyaev-Zel'dovich observations of MaDCoWS clusters, we also present a preliminary mass-richness relation that can be used to infer the approximate mass distribution of the full sample. The estimated median mass for the WISE--Pan-STARRS catalog is $M_{500}=1.6^{+0.7}_{-0.8}\times10^{14} \mathrm{M}_\odot$, with the Sunyaev-Zel'dovich data confirming that we detect clusters with masses up to $M_{500}\sim5\times10^{14} \mathrm{M}_\odot$ $(M_{200}\sim10^{15} \mathrm{M}_\odot)$.

We present the results from a pilot study with the Karl G. Jansky Very Large Array (JVLA) to determine the radio morphologies of extended radio sources and the properties of their host-galaxies in 10 massive galaxy clusters at z~1, an epoch in which clusters are assembling rapidly. These clusters are drawn from a parent sample of WISE-selected galaxy clusters that were cross-correlated with the VLA Faint Images of the Radio Sky at Twenty-Centimeters survey (FIRST) to identify extended radio sources within 1$^{\prime}$ of the cluster centers. Out of the ten targeted sources, six are FR II sources, one is an FR I source, and three sources have undetermined morphologies. Eight radio sources have associated Spitzer data, 75% presenting infrared counterparts. A majority of these counterparts are consistent with being massive galaxies. The angular extent of the FR sources exhibits a strong correlation with the cluster-centric radius, which warrants further investigation with a larger sample.

We have observed the environments of a population of 33 heavily dust obscured, ultra-luminous, high-redshift galaxies, selected using WISE and NVSS at $z>$1.3 with the Infra-Red Array Camera on the $Spitzer$ Space Telescope over $\rm5.12\,'\times5.12\,'$ fields. Colour selections are used to quantify any potential overdensities of companion galaxies in these fields. We find no significant excess of galaxies with the standard colour selection for IRAC colours of $\rm[3.6]-[4.5]>-0.1$ consistent with galaxies at $z>$1.3 across the whole fields with respect to wide-area $Spitzer$ comparison fields, but there is a $\rm>2\sigma$ statistical excess within $\rm0.25\,'$ of the central radio-WISE galaxy. Using a colour selection of $\rm[3.6]-[4.5]>0.4$, 0.5 magnitudes redder than the standard method of selecting galaxies at $z>$1.3, we find a significant overdensity, in which $\rm76\%$ ($\rm33\%$) of the 33 fields have a surface density greater than the $\rm3\sigma$ ($\rm5\sigma$) level. There is a statistical excess of these redder galaxies within $\rm0.5\,'$, rising to a central peak $\rm\sim2$--4 times the average density. This implies that these galaxies are statistically linked to the radio-WISE selected galaxy, indicating similar structures to those traced by red galaxies around radio-loud AGN.

Extended narrow-line regions (NLRs) around active galactic nuclei (AGN) are shaped by the distribution of gas in the host galaxy and by the geometry of the circumnuclear obscuration, and thus they can be used to test the AGN unification model. In this work, we quantify the morphologies of the narrow-line regions in 308 nearby AGNs ($z=0-0.14$, \lbol $\sim 10^{42.4-44.1}$ \erg{}) from the MaNGA survey. Based on the narrow-line region maps, we find that a large fraction (81\%) of these AGN have bi-conical NLR morphology. The distribution of their measured opening angles suggests that the intrinsic opening angles of the ionization cones has a mean value of 85--98$^\circ$ with a finite spread of 39-44$^\circ$ (1-$\sigma$). Our inferred opening angle distribution implies a number ratio of type I to type II AGN of 1:1.6--2.3, consistent with other measurements of the type I / type II ratio at low AGN luminosities. Combining these measurements with the WISE photometry data, we find that redder mid-IR color (lower effective temperature of dust) corresponds to stronger and narrower photo-ionized bicones. This relation is in agreement with the unification model that suggests that the bi-conical narrow-line regions are shaped by a toroidal dusty structure within a few pc from the AGN. Furthermore, we find a significant alignment between the minor axis of host galaxy disks and AGN ionization cones. Together, these findings suggest that obscuration on both circumnuclear ($\sim $pc) and galactic ($\sim$ kpc) scales are important in shaping and orienting the AGN narrow-line regions.

We present an analysis of the radial distribution of Active Galactic Nuclei (AGN) in $2300$ galaxy clusters from the Massive and Distant Clusters of {\it WISE} Survey (MaDCoWS). MaDCoWS provides the largest coverage of the extragalactic sky for a cluster sample at $z\sim1$. We use literature catalogs of AGN selected via optical, mid-infrared (MIR), and radio data, and by optical-to-MIR (OIR) color. Stacking the radial distribution of AGN within the $6\arcmin$ of the centers of MaDCoWS galaxy clusters, we find a distinct overdensity of AGN within $1\arcmin$ of the galaxy cluster center for AGN of all selection methods. The fraction of red galaxies that host AGN as a function of clustercentric distance is, however, dependent on the AGN selection. The fraction of red galaxies in cluster environments that host AGN selected by optical signatures or blue OIR color is at a deficit compared to the field, while MIR-selected and red OIR color AGN are enhanced in the centers of clusters when compared to field levels. The radio-selected AGN fraction is more than $2.5$ times that of the field, implying that the centers of clusters are conducive to the triggering of radio emission in AGN. We do not find a statistically significant change in the AGN fraction as a function of cluster richness. We also investigate the correlation of central radio activity with other AGN in galaxy clusters. Clusters with radio activity have more central AGN than radio-inactive clusters, implying that central cluster radio activity and AGN triggering may be linked.

Negative feedback from accretion onto super-massive black holes (SMBHs), that is to remove gas and suppress star formation in galaxies, has been widely suggested. However, for Seyfert galaxies which harbor less active, moderately accreting SMBHs in the local universe, the feedback capability of their black hole activity is elusive. We present spatially-resolved H$\alpha$ measurements to trace ongoing star formation in Seyfert galaxies and compare their specific star formation rate with a sample of star-forming galaxies whose global galaxy properties are controlled to be the same as the Seyferts. From the comparison we find that the star formation rates within central kpc of Seyfert galaxies are mildly suppressed as compared to the matched normal star forming galaxies. This suggests that the feedback of moderate SMBH accretion could, to some extent, regulate the ongoing star formation in these intermediate to late type galaxies under secular evolution.

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.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

We report spectroscopic results from our 40-orbit $Hubble~Space~Telescope$ slitless grism spectroscopy program observing the 20 densest CARLA candidate galaxy clusters at $1.4 < z < 2.8$. These candidate rich structures, among the richest and most distant known, were identified on the basis of $[3.6]-[4.5]$ color from a $408~$hr multi-cycle $Spitzer$ program targeting $420$ distant radio-loud AGN. We report the spectroscopic confirmation of $16$ distant structures at $1.4 < z < 2.8$ associated with the targeted powerful high-redshift radio-loud AGN. We also report the serendipitous discovery and spectroscopic confirmation of seven additional structures at $0.87 < z < 2.12$ not associated with the targeted radio-loud AGN. We find that $10^{10} - 10^{11}\,M_{\odot}$ member galaxies of our confirmed CARLA structures form significantly fewer stars than their field counterparts at all redshifts within $1.4\leq z\leq 2$. We also observe higher star-forming activity in the structure cores up to $z = 2$, finding similar trends as cluster surveys at slightly lower redshifts ($1.0 < z < 1.5$). By design, our efficient strategy of obtaining just two grism orbits per field only obtains spectroscopic confirmation of emission-line galaxies. Deeper spectroscopy will be required to study the population of evolved, massive galaxies in these (forming) clusters. Lacking multi-band coverage of the fields, we adopt a very conservative approach of calling all confirmations "structures", although we note that a number of features are consistent with some of them being bona fide galaxy clusters. Together, this survey represents a unique and large homogenous sample of spectroscopically confirmed structures at high redshifts, potentially more than doubling the census of confirmed, massive clusters at $z > 1.4$.

Different mechanisms can drive outflows in active galactic nuclei (AGN), but it is often unclear which mechanism dominates, if any. To quantify the impact of AGN feedback on galaxy evolution, the driving mechanism of outflows must be better understood.

We present new evidence for AGN feedback in a subset of 69 quenched low-mass galaxies ($M_{\star} \lesssim 5\times10^{9}$ M$_{\odot}$, $M_{\rm{r}} > -19$) selected from the first two years of the SDSS-IV MaNGA survey. The majority (85 per cent) of these quenched galaxies appear to reside in a group environment. We find 6 galaxies in our sample that appear to have an active AGN that is preventing on-going star-formation; this is the first time such a feedback mechanism has been observed in this mass range. Interestingly, five of these six galaxies have an ionised gas component that is kinematically offset from their stellar component, suggesting the gas is either recently accreted or outflowing. We hypothesise these six galaxies are low-mass equivalents to the "red geysers" observed in more massive galaxies. Of the other 63 galaxies in the sample, we find 8 do appear for have some low-level, residual star formation, or emission from hot, evolved stars. The remaining galaxies in our sample have no detectable ionised gas emission throughout their structures, consistent with them being quenched. This work shows the potential for understanding the detailed physical properties of dwarf galaxies through spatially resolved spectroscopy.

Quasars may have played a key role in limiting the stellar mass of massive galaxies. Identifying those quasars in the process of removing star formation fuel from their hosts is an exciting ongoing challenge in extragalactic astronomy. In this paper we present X-ray observations of eleven extremely red quasars (ERQs) with $L_{\rm bol}\sim 10^{47}$ erg s$^{-1}$ at $z=1.5-3.2$ with evidence for high-velocity ($v > 1000$ km s$^{-1}$) [OIII]$\lambda$5007\AA\ outflows. X-rays allow us to directly probe circumnuclear obscuration and to measure the instantaneous accretion luminosity. We detect ten out of eleven extremely red quasars available in targeted and archival data. Using a combination of X-ray spectral fitting and hardness ratios, we find that all of the ERQs show signs of absorption in the X-rays with inferred column densities of $N_{\rm H}\approx 10^{23}$ cm$^{-2}$, including four Compton-thick candidates ($N_{\rm H} > 10^{24}$ cm$^{-2}$). We stack the X-ray emission of the seven weakly detected sources, measuring an average column density of $N_{\rm H}\sim 8\times 10^{23}$ cm$^{-2}$. The absorption-corrected (intrinsic) $2-10$ keV X-ray luminosity of the stack is $2.7\times 10^{45}$ erg s$^{-1}$, consistent with X-ray luminosities of type 1 quasars of the same infrared luminosity. Thus, we find that ERQs are a highly obscured, borderline Compton-thick population, and based on optical and infrared data we suggest that these objects are partially hidden by their own equatorial outflows. However, unlike some quasars with known outflows, ERQs do not appear to be intrinsically underluminous in X-rays for their bolometric luminosity. Our observations indicate that low X-rays are not necessary to enable some types of radiatively driven winds.

In this paper, we investigate 2727 galaxies observed by MaNGA as of June 2016 to develop spatially resolved techniques for identifying signatures of active galactic nuclei (AGN). We identify 303 AGN candidates. The additional spatial dimension imposes challenges in identifying AGN due to contamination from diffuse ionized gas, extra-planar gas and photoionization by hot stars. We show that the combination of spatially-resolved line diagnostic diagrams and additional cuts on H$\alpha$ surface brighness and H$\alpha$ equivalent width can distinguish between AGN-like signatures and high-metallicity galaxies with LINER-like spectra. Low mass galaxies with high specific star formation rates are particularly difficult to diagnose and routinely show diagnostic line ratios outside of the standard star-formation locus. We develop a new diagnostic -- the distance from the standard diagnostic line in the line-ratios space -- to evaluate the significance of the deviation from the star-formation locus. We find 173 galaxies that would not have been selected as AGN candidates based on single-fibre spectral measurements but exhibit photoionization signatures suggestive of AGN activity in the MaNGA resolved observations, underscoring the power of large integral field unit (IFU) surveys. A complete census of these new AGN candidates is necessary to understand their nature and probe the complex co-evolution of supermassive black holes and their hosts.

We report the characterization of the first $62$ MaNGA Active Galactic Nuclei (AGN) hosts in the Fifth Product Launch (MPL-5) and the definition of a control sample of non-active galaxies. This control sample - comprising two galaxies for each AGN - was selected in order to match the AGN hosts in terms of stellar mass, redshift, visual morphology and inclination. The stellar masses are in the range $9.4<\log(M/M_\odot)<11.5$, and most objects have redshifts $\leq 0.08$. The AGN sample is mostly comprised of low-luminosity AGN, with only 17 nuclei with $L([OIII]\lambda 5007)\geq 3.8\times 10^{40}$ erg s$^{-1}$ (that we call "strong AGN"). The stellar population of the control sample galaxies within the inner $1$-$3$ kpc is dominated by the old ($\sim$ $4$ - $13$ Gyr) age component, with a small contribution of intermediate age ($\sim 640$-$940$ Myr) and young stars ($\leq 40$ Myr) to the total light at $5700\AA$. While the weaker AGN show a similar age distribution to that of the control galaxies, the strong AGN show an increased contribution of younger stars and a decreased contribution of older stars. Examining the relationship between the AGN stellar population properties and $L([OIII])$, we find that with increasing $L([OIII])$, the AGN exhibit a decreasing contribution from the oldest ($>4$ Gyr) stellar population relative to control galaxies, but have an increasing contribution from the younger components with ages $\sim 40$ Myr. We also find a correlation of the mean age differences (AGN - control) with $L([OIII])$, in the sense that more luminous AGN are younger than the control objects, while the low-luminosity AGN seem to be older. These results support a connection between the growth of the galaxy bulge via formation of new stars and the growth of the Supermassive Black Hole via matter accretion in the AGN phase.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

It is not yet known if the properties of molecular gas in distant protocluster galaxies are significantly affected by their environment as galaxies are in local clusters. Through a deep, 64 hours of effective on-source integration with the Australian Telescope Compact Array (ATCA), we discovered a massive, M_mol=2.0+-0.2x10^11 M_sun, extended, ~40 kpc, CO(1-0)-emitting disk in the protocluster surrounding the radio galaxy, MRC1138-262. The galaxy, at z_CO=2.1478, is a clumpy, massive disk galaxy, M_star~5x10^11 M_sun, which lies 250 kpc in projection from MRC1138-262 and is a known H-alpha emitter, HAE229. HAE229 has a molecular gas fraction of ~30%. The CO emission has a kinematic gradient along its major axis, centered on the highest surface brightness rest-frame optical emission, consistent with HAE229 being a rotating disk. Surprisingly, a significant fraction of the CO emission lies outside of the UV/optical emission. In spite of this, HAE229 follows the same relation between star-formation rate and molecular gas mass as normal field galaxies. HAE229 is the first CO(1-0) detection of an ordinary, star-forming galaxy in a protocluster. We compare a sample of cluster members at z>0.4 that are detected in low-order CO transitions with a similar sample of sources drawn from the field. We confirm findings that the CO-luminosity and FWHM are correlated in starbursts and show that this relation is valid for normal high-z galaxies as well as those in overdensities. We do not find a clear dichotomy in the integrated Schmidt-Kennicutt relation for protocluster and field galaxies. Not finding any environmental dependence in the "star-formation efficiency" or the molecular gas content, especially for such an extended CO disk, suggests that environmentally-specific processes such as ram pressure stripping are not operating efficiently in (proto)clusters. (abridged)

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

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.

Ionised gas outflows driven by active galactic nuclei (AGN) are ubiquitous in high luminosity AGN with outflow speeds apparently correlated with the total bolometric luminosity of the AGN. This empirical relation and theoretical work suggest that in the range L_bol ~ 10^43-45 erg/s there must exist a threshold luminosity above which the AGN becomes powerful enough to launch winds that will be able to escape the galaxy potential. In this paper, we present pilot observations of two AGN in this transitional range that were taken with the Gemini North Multi-Object Spectrograph Integral Field Unit (IFU). Both sources have also previously been observed within the Sloan Digital Sky Survey-IV (SDSS) Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. While the MaNGA IFU maps probe the gas fields on galaxy-wide scales and show that some regions are dominated by AGN ionization, the new Gemini IFU data zoom into the centre with four times better spatial resolution. In the object with the lower L_bol we find evidence of a young or stalled biconical AGN-driven outflow where none was obvious at the MaNGA resolution. In the object with the higher L_bol we trace the large-scale biconical outflow into the nuclear region and connect the outflow from small to large scales. These observations suggest that AGN luminosity and galaxy potential are crucial in shaping wind launching and propagation in low-luminosity AGN. The transition from small and young outflows to galaxy-wide feedback can only be understood by combining large-scale IFU data that trace the galaxy velocity field with higher resolution, small scale IFU maps.

Using HST slitless grism data, we report the spectroscopic confirmation of two distant structures at $z \sim 2$ associated with powerful high-redshift radio-loud AGN. These rich structures, likely (forming) clusters, are among the most distant currently known and were identified on the basis of Spitzer/IRAC [3.6] - [4.5] color. We spectroscopically confirm 9 members in the field of MRC 2036-254, comprising eight star-forming galaxies and the targeted radio galaxy. The median redshift is $z = 2.000$. We spectroscopically confirm 10 members in the field of B3 0756+406, comprising eight star-forming galaxies and two AGN, including the targeted radio-loud quasar. The median redshift is $z = 1.986$. All confirmed members are within 500 kpc (1 arcmin) of the targeted AGN. We derive median (mean) star-formation rates of $\sim 35~M_{\odot}\rm ~ yr^{-1}$ ($\sim 50~M_{\odot}\rm ~ yr^{-1}$) for the confirmed star-forming members of both structures based on their [OIII]$\lambda5007$ luminosities, and estimate average galaxy stellar masses $\lesssim 1 \times 10^{11} ~M_{\odot}$ based on mid-infrared fluxes and SED modeling. Most of our confirmed members are located above the star-forming main-sequence towards starburst galaxies, consistent with clusters at these early epochs being the sites of significant levels of star formation. The structure around MRC 2036-254 shows an overdensity of IRAC-selected candidate galaxy cluster members consistent with being quiescent galaxies, while the structure around B3 0756+406 shows field values, albeit with many lower limits to colors that could allow an overdensity of faint red quiescent galaxies. The structure around MRC 2036-254 shows a red sequence of passive galaxy candidates.

Feedback from active galactic nuclei (AGN) is widely considered to be the main driver in regulating the growth of massive galaxies through heating or driving gas out of the galaxy, preventing further increase in stellar mass. Observational proof for this scenario has, however, been scarce. We have assembled a sample of 132 radio-quiet type-2 and red AGN at 0.1<z<1. We measure the kinematics of the AGN-ionized gas, the host galaxies' stellar masses and star formation rates and investigate the relationships between AGN luminosities, specific star formation rates (sSFR) and outflow strengths W90 -- the 90\% velocity width of the [OIII]5007 line power and a proxy for the AGN-driven outflow speed. Outflow strength is independent of sSFR for AGN selected on their mid-IR luminosity, in agreement with previous work demonstrating that star formation is not sufficient to produce the observed ionized gas outflows which have to be powered by AGN activity. More importantly, we find a negative correlation between W90 and sSFR in the AGN hosts with the highest SFRs, i.e., with the highest gas content, where presumably the coupling of the AGN-driven wind to the gas is strongest. This implies that AGN with strong outflow signatures are hosted in galaxies that are more `quenched' than galaxies with weaker outflow signatures. Despite the galaxies' high SFRs, we demonstrate that the outflows are not star-formation driven but indeed due to AGN-powering. This observation is consistent with the AGN having a net suppression, `negative' impact, through feedback on the galaxies' star formation history.

Black hole feedback is now a standard component of galaxy formation models. These models predict that the impact of black hole activity on its host galaxy likely peaked at z=2-3, the epoch of strongest star formation activity and black hole accretion activity in the Universe. We used XShooter on the Very Large Telescope to measure rest-frame optical spectra of four z~2.5 extremely red quasars with infrared luminosities ~10^47 erg/sec. We present the discovery of very broad (full width at half max= 2600-5000 km/sec), strongly blue-shifted (by up to 1500 km/sec) [OIII]5007A emission lines in these objects. In a large sample of obscured and red quasars, [OIII] kinematics are positively correlated with infrared luminosity, and the four objects in our sample are on the extreme end both in [OIII] kinematics and infrared luminosity. We estimate that ~3% of the bolometric luminosity in these objects is being converted into the kinetic power of the observed wind. Photo-ionization estimates suggest that the [OIII] emission might be extended on a few kpc scales, which would suggest that the extreme outflow is affecting the entire host galaxy of the quasar. These sources may be the signposts of the most extreme form of quasar feedback at the peak epoch of galaxy formation, and may represent an active "blow-out" phase of quasar evolution.

Luminous type-2 quasars in which the glow from the central black hole is obscured by dust are ideal targets for studying their host galaxies and the quasars' effect on galaxy evolution. Such feedback appears ubiquitous in luminous obscured quasars where high velocity ionized nebulae have been found. We present rest-frame yellow-band (~5000 Angstroms) observations using the Hubble Space Telescope for a sample of 20 luminous quasar host galaxies at 0.2 < z < 0.6 selected from the Sloan Digital Sky Survey. For the first time, we combine host galaxy observations with geometric measurements of quasar illumination using blue-band HST observations and [OIII] integral field unit observations probing the quasar winds. The HST images reveal bright merger signatures in about half the galaxies; a significantly higher fraction than in comparison inactive ellipticals. We show that the host galaxies are primarily bulge-dominated, with masses close to M*, but belong to < 30% of elliptical galaxies that are highly star-forming at z ~ 0.5. Ionized gas signatures are uncorrelated with faint stellar disks (if present), confirming that the ionized gas is not concentrated in a disk. Scattering cones and [OIII] ionized gas velocity field are aligned with the forward scattering cones being co-spatial with the blue-shifted side of the velocity field, suggesting the high velocity gas is indeed photo-ionized by the quasar. Based on the host galaxies' high star-formation rates and bright merger signatures, we suggest that this low-redshift outbreak of luminous quasar activity is triggered by recent minor mergers. Combining these novel observations, we present new quasar unification tests, which are in agreement with expectations of the orientation-based unification model for quasars.

We analyze Hubble Space Telescope observations of scattering regions in 20 luminous obscured quasars at $0.24<z<0.65$ (11 new observations and 9 archival ones) observed at rest-frame $\sim 3000$\AA. We find spectacular $5-10$ kpc-scale scattering regions in almost all cases. The median scattering efficiency at this wavelength (the ratio of observed to estimated intrinsic flux) is 2.3\%, and 73\% of the observed flux at this wavelength is due to scattered light, which if unaccounted for may strongly bias estimates of quasar hosts' star formation rates. Modeling these regions as illuminated dusty cones, we estimate the radial density distributions of the interstellar medium as well as the geometric properties of circumnuclear quasar obscuration -- inclinations and covering factors. Small derived opening angles (median half-angle and standard deviation 27\dg$\pm$9\dg) are inconsistent with a 1:1 type 1 / type 2 ratio. We suggest that quasar obscuration is patchy and that the observer has a $\sim 40\%$ chance of seeing a type 1 source even through the obscuration. We estimate median density profile of the scattering medium to be $n_{\rm H}=0.04-0.5$ $(1{\rm kpc}/r)^2$ cm$^{-3}$, depending on the method. Quasars in our sample likely exhibit galaxy-wide winds, but if these consist of optically thick clouds then only a small fraction of the wind mass ($\la 10\%$) contributes to scattering.

We report on the discovery of a z=1.58 mature cluster around the high-redshift radio galaxy 7C1753+6311, first identified in the Clusters Around Radio-Loud AGN survey. Two-thirds of the excess galaxies within the central 1Mpc lie on a red sequence with a colour that is consistent with an average formation redshift of zf~3. We show that 80+/-6% of the red sequence galaxies in the cluster core are quiescent, while the remaining 20% are red due to dusty star formation. We demonstrate that the cluster has an enhanced quiescent galaxy fraction that is three times that of the control field. We also show that this enhancement is mass dependent: 91+/-9% of the M* >10^{10.5}Msun cluster galaxies are quiescent, compared to only 36+/-2% of field galaxies, whereas the fraction of quiescent galaxies with lower masses is the same in the cluster and field environments. The presence of a dense core and a well-formed, quiescent red sequence suggest that this is a mature cluster. This means that distant radio galaxies do not solely reside in young, uncollapsed protoclusters, rather they can be found in clusters in a wide range of evolutionary states.

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.

We present confirmation of the cluster MOO J1142+1527, a massive galaxy cluster discovered as part of the Massive and Distant Clusters of WISE Survey. The cluster is confirmed to lie at $z=1.19$, and using the Combined Array for Research in Millimeter-wave Astronomy we robustly detect the Sunyaev-Zel'dovich (SZ) decrement at 13.2$\sigma$. The SZ data imply a mass of $\mathrm{M}_{200m}=(1.1\pm0.2)\times10^{15}$ $\mathrm{M}_\odot$, making MOO J1142+1527 the most massive galaxy cluster known at $z>1.15$ and the second most massive cluster known at $z>1$. For a standard $\Lambda$CDM cosmology it is further expected to be one of the $\sim 5$ most massive clusters expected to exist at $z\ge1.19$ over the entire sky. Our ongoing Spitzer program targeting $\sim1750$ additional candidate clusters will identify comparably rich galaxy clusters over the full extragalactic sky.

We use a sample of 37 of the densest clusters and protoclusters across $1.3 \le z \le 3.2$ from the Clusters Around Radio-Loud AGN (CARLA) survey to study the formation of massive cluster galaxies. We use optical $i'$-band and infrared 3.6$\mu$m and 4.5$\mu$m images to statistically select sources within these protoclusters and measure their median observed colours; $\langle i'-[3.6] \rangle$. We find the abundance of massive galaxies within the protoclusters increases with decreasing redshift, suggesting these objects may form an evolutionary sequence, with the lower redshift clusters in the sample having similar properties to the descendants of the high redshift protoclusters. We find that the protocluster galaxies have an approximately unevolving observed-frame $i'-[3.6]$ colour across the examined redshift range. We compare the evolution of the $\langle i'-[3.6] \rangle$ colour of massive cluster galaxies with simplistic galaxy formation models. Taking the full cluster population into account, we show that the formation of stars within the majority of massive cluster galaxies occurs over at least 2Gyr, and peaks at $z \sim 2$-3. From the median $i'-[3.6]$ colours we cannot determine the star formation histories of individual galaxies, but their star formation must have been rapidly terminated to produce the observed red colours. Finally, we show that massive galaxies at $z>2$ must have assembled within 0.5Gyr of them forming a significant fraction of their stars. This means that few massive galaxies in $z>2$ protoclusters could have formed via dry mergers.

We report the detection of extended warm ionized gas in two powerful high-redshift radio galaxies, NVSS J210626-314003 at z=2.10 and TXS 2353-003 at z=1.49, that does not appear to be associated with the radio jets. This is contrary to what would be expected from the alignment effect, a characteristic feature of distant, powerful radio galaxies at z> 0.6. The gas also has smaller velocity gradients and line widths than most other high-z radio galaxies with similar data. Both galaxies are part of a systematic study of 50 high-redshift radio galaxies with SINFONI, and are the only two that are characterized by the presence of high surface-brightness gas not associated with the jet axis and by the absence of such gas aligned with the jet. Both galaxies are spatially resolved with ISAAC broadband imaging covering the rest-frame R band, and have extended wings that cannot be attributed to line contamination. We argue that the gas and stellar properties of these galaxies are more akin to gas-rich brightest cluster galaxies in cool-core clusters than the general population of high-redshift radio galaxies at z>2. In support of this interpretation, one of our sources, TXS 2353-003, for which we have H\alpha\ narrowband imaging, is associated with an overdensity of candidate H\alpha\ emitters by a factor of 8 relative to the field at z=1.5. We discuss possible scenarios of the evolutionary state of these galaxies and the nature of their emission line gas within the context of cyclical AGN feedback.

We present CARMA 30 GHz Sunyaev-Zel'dovich (SZ) observations of five high-redshift ($z \gtrsim 1$), infrared-selected galaxy clusters discovered as part of the all-sky Massive and Distant Clusters of WISE Survey (MaDCoWS). The SZ decrements measured toward these clusters demonstrate that the MaDCoWS selection is discovering evolved, massive galaxy clusters with hot intracluster gas. Using the SZ scaling relation calibrated with South Pole Telescope clusters at similar masses and redshifts, we find these MaDCoWS clusters have masses in the range $M_{200} \approx 2-6 \times 10^{14}$ $M_\odot$. Three of these are among the most massive clusters found to date at $z\gtrsim 1$, demonstrating that MaDCoWS is sensitive to the most massive clusters to at least $z = 1.3$. The added depth of the AllWISE data release will allow all-sky infrared cluster detection to $z \approx 1.5$ and beyond.

The WISE mission has unveiled a rare population of high-redshift ($z=1-4.6$), dusty, hyper-luminous galaxies, with infrared luminosities $L_{\rm IR} > 10^{13}~L_{\odot}$, and sometimes exceeding $10^{14}~L_{\odot}$. Previous work has shown that their dust temperatures and overall far-IR SEDs are significantly hotter than expected for star-formation. We present here an analysis of the rest-frame optical through mid-IR SEDs for a large sample of these so-called "Hot, Dust-Obscured Galaxies" (Hot DOGs). We find that the SEDs of Hot DOGs are generally well modeled by the combination of a luminous, yet obscured AGN that dominates the rest-frame emission at $\lambda > 1\mu\rm m$ and the bolometric luminosity output, and a less luminous host galaxy that is responsible for the bulk of the rest optical/UV emission. Even though the stellar mass of the host galaxies may be as large as $10^{11}-10^{12}~M_{\odot}$, the AGN emission, with luminosities comparable to those of the most luminous QSOs known, require that either Hot DOGs have black hole masses significantly in excess of the local relations, or that they radiate significantly above the Eddington limit. We show that, while rare, the number density of Hot DOGs is comparable to that of equally luminous but unobscured (i.e., Type 1) QSOs. This is inconsistent with the trend of a diminishing fraction of obscured objects with increasing luminosity found for less luminous QSOs, possibly indicating a reversal in this relation at high luminosity, and that Hot DOGs are not the torus-obscured counterparts of the known optically selected, largely unobscured Hyper-Luminous QSOs. Hot DOGs may represent a different type of galaxy and thus a new component of the galaxy evolution paradigm. Finally, we discuss the environments of Hot DOGs and show that these objects are in regions as dense as those of known high-redshift proto-clusters.(Abridged)

We present APEX LABOCA 870 micron observations of the field around the high-redshift radio galaxy MRC1138-262 at z=2.16. We detect 16 submillimeter galaxies in this ~140 square arcmin bolometer map with flux densities in the range 3-11 mJy. The raw number counts indicate a density of submillimeter galaxies (SMGs) that is up to four times that of blank field surveys. Based on an exquisite multiwavelength database, including VLA 1.4 GHz radio and infrared observations, we investigate whether these sources are members of the protocluster structure at z=2.2. Using Herschel PACS+SPIRE and Spitzer MIPS photometry, we derive reliable far-infrared photometric redshifts for all sources. Follow-up VLT ISAAC and SINFONI near-infrared spectra confirm that four of these SMGs have redshifts of z=2.2. We also present evidence that another SMG in this field, detected earlier at 850 micron, has a counterpart that exhibits Halpha and CO(1-0) emission at z=2.15. Including the radio galaxy and two SMGs with far-IR photometric redshifts at z=2.2, we conclude that at least eight submm sources are part of the protocluster at z=2.16 associated with the radio galaxy MRC1138-262. We measure a star formation rate density SFRD ~1500 Msun yr^-1 Mpc^-3, four magnitudes higher than the global SFRD of blank fields at this redshift. Strikingly, these eight sources are concentrated within a region of 2 Mpc (the typical size of clusters in the local universe) and are distributed within the filaments traced by the Halpha emitters at z=2.2. This concentration of massive, dusty starbursts is not centered on the submillimeter-bright radio galaxy which could support the infalling of these sources into the cluster center. Approximately half (6/11) of the SMGs that are covered by the Halpha imaging data are associated with Halpha emitters, demonstrating the potential of tracing SMG counterparts with this population (abridged).

Distant powerful radio-loud active galactic nuclei (RLAGN) tend to reside in dense environments and are commonly found in proto-clusters at z > 1.3. We examine whether this occurs because RLAGN are hosted by massive galaxies, which preferentially reside in rich environments. We compare the environments of powerful RLAGN at 1.3 < z < 3.2 from the CARLA survey to a sample of radio-quiet galaxies matched in mass and redshift. We find the environments of RLAGN are significantly denser than those of radio-quiet galaxies, implying that not more than 50% of massive galaxies in this epoch can host powerful radio-loud jets. This is not an observational selection effect as we find no evidence to suggest it is easier to observe the radio emission when the galaxy resides in a dense environment. We therefore suggest that the dense Mpc-scale environment fosters the formation of a radio-jet from an AGN. We show that the number density of potential RLAGN host galaxies is consistent with every > 10^14 solar mass cluster having experienced powerful radio-loud feedback of duration ~60 Myr during 1.3 < z < 3.2. This feedback could heat the intracluster medium to the extent of 0.5-1 keV per gas particle, which could limit the amount of gas available for further star formation in the proto-cluster galaxies.

We present optical and infrared imaging and optical spectroscopy of galaxy clusters which were identified as part of an all-sky search for high-redshift galaxy clusters, the Massive and Distant Clusters of WISE Survey (MaDCoWS). The initial phase of MaDCoWS combined infrared data from the all-sky data release of the Wide-field Infrared Survey Explorer (WISE) with optical data from the Sloan Digital Sky Survey (SDSS) to select probable z ~ 1 clusters of galaxies over an area of 10,000 deg^2. Our spectroscopy confirms 19 new clusters at 0.7 < z < 1.3, half of which are at z > 1, demonstrating the viability of using WISE to identify high-redshift galaxy clusters. The next phase of MaDCoWS will use the greater depth of the AllWISE data release to identify even higher redshift cluster candidates.

Herschel Space Observatory photometry and extensive multiwavelength followup have revealed that the powerful radio galaxy 3C 220.3 at z=0.685 acts as a gravitational lens for a background submillimeter galaxy (SMG) at z=2.221. At an observed wavelength of 1mm, the SMG is lensed into three distinct images. In the observed near infrared, these images are connected by an arc of 1.8" radius forming an Einstein half-ring centered near the radio galaxy. In visible light, only the arc is apparent. 3C 220.3 is the only known instance of strong galaxy-scale lensing by a powerful radio galaxy not located in a galaxy cluster and therefore it offers the potential to probe the dark matter content of the radio galaxy host. Lens modeling rejects a single lens, but two lenses centered on the radio galaxy host A and a companion B, separated by 1.5", provide a fit consistent with all data and reveal faint candidates for the predicted fourth and fifth images. The model does not require an extended common dark matter halo, consistent with the absence of extended bright X-ray emission on our Chandra image. The projected dark matter fractions within the Einstein radii of A (1.02") and B (0.61") are about 0.4 +/- 0.3 and 0.55 +/- 0.3. The mass to i-band light ratios of A and B, M/L ~ 8 +/- 4 Msun/Lsun, appear comparable to those of radio-quiet lensing galaxies at the same redshift in the CASTLES, LSD, and SL2S samples. The lensed SMG is extremely bright with observed f(250um) = 440mJy owing to a magnification factor mu~10. The SMG spectrum shows luminous, narrow CIV 154.9nm emission, revealing that the SMG houses a hidden quasar in addition to a violent starburst. Multicolor image reconstruction of the SMG indicates a bipolar morphology of the emitted ultraviolet (UV) light suggestive of cones through which UV light escapes a dust-enshrouded nucleus.

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 4.5 {\mu}m luminosity functions for galaxies identified in 178 candidate galaxy clusters at 1.3 < z < 3.2. The clusters were identified as Spitzer/IRAC color-selected overdensities in the Clusters Around Radio-Loud AGN (CARLA) project, which imaged 421 powerful radio-loud AGN at z > 1.3. The luminosity functions are derived for different redshift and richness bins, and the IRAC imaging reaches depths of m*+2, allowing us to measure the faint end slopes of the luminosity functions. We find that {\alpha} = -1 describes the luminosity function very well in all redshifts bins and does not evolve significantly. This provides evidence that the rate at which the low mass galaxy population grows through star formation, gets quenched and is replenished by in-falling field galaxies does not have a major net effect on the shape of the luminosity function. Our measurements for m* are consistent with passive evolution models and high formation redshifts z_f ~ 3. We find a slight trend towards fainter m* for the richest clusters, implying that the most massive clusters in our sample could contain older stellar populations, yet another example of cosmic downsizing. Modelling shows that a contribution of a star-forming population of up to 40% cannot be ruled out. This value, found from our targeted survey, is significantly lower than the values found for slightly lower redshift, z ~ 1, clusters found in wide-field surveys. The results are consistent with cosmic downsizing, as the clusters studied here were all found in the vicinity of radio-loud AGNs -- which have proven to be preferentially located in massive dark matter halos in the richest environments at high redshift -- they may therefore be older and more evolved systems than the general protocluster population.

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 the spectroscopic confirmation of a structure of galaxies surrounding the radio galaxy MRC0156-252 at z = 2.02. The structure was initially discovered as an overdensity of both near-infrared selected z > 1.6 and mid-infrared selected z > 1.2 galaxy candidates. We used the VLT/FORS2 multi-object spectrograph to target ~80 high-redshift galaxy candidates, and obtain robust spectroscopic redshifts for more than half the targets. The majority of the confirmed sources are star-forming galaxies at z > 1.5. In addition to the radio galaxy, two of its close-by companions (< 6'') also show AGN signatures. Ten sources, including the radio galaxy, lie within |z - 2.020 | < 0.015 (i.e., velocity offsets < 1500 km/s) and within projected 2 Mpc comoving of the radio galaxy. Additional evidence suggests not only that the galaxy structure associated with MRC0156-252 is a forming galaxy cluster but also that this structure is most probably embedded in a larger scale structure.

We report the first results from the Clusters Around Radio-Loud AGN (CARLA) program, a Cycle 7 and 8 Spitzer Space Telescope snapshot program to investigate the environments of a large sample of obscured and unobscured luminous radio-loud AGN at 1.2 < z < 3.2. These data, obtained for 387 fields, reach 3.6 and 4.5 um depths of [3.6] (AB) = 22.6 and [4.5] (AB) = 22.9 at the 95% completeness level, which is two to three times fainter than L* in this redshift range. By using the color cut [3.6]-[4.5] > -0.1 (AB), which efficiently selects high-redshift (z > 1.3) galaxies of all types, we identify galaxy cluster member candidates in the fields of the radio-loud AGN. The local density of these IRAC-selected sources is compared to the density of similarly selected sources in blank fields. We find that 92% of the radio-loud AGN reside in environments richer than average. The majority (55%) of the radio-loud AGN fields are found to be overdense at a > 2 {\sigma} level; 10% are overdense at a > 5 {\sigma} level. A clear rise in surface density of IRAC-selected sources towards the position of the radio-loud AGN strongly supports an association of the majority of the IRAC-selected sources with the radio-loud AGN. Our results provide solid statistical evidence that radio-loud AGN are likely beacons for finding high-redshift galaxy (proto-)clusters. We investigate how environment depends on AGN type (unobscured radio-loud quasars vs. obscured radio galaxies), radio luminosity and redshift, finding no correlation with either AGN type or radio luminosity. We find a decrease in density with redshift, consistent with galaxy evolution for this uniform, flux-limited survey. These results are consistent with expectations from the orientation-driven AGN unification model, at least for the high radio luminosity regimes considered in this sample.

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

We report new, sensitive observations of two z ~ 3-3.5 FIR-luminous radio galaxies, 6C1909+72 and B3J2330+3927, in 12CO J=1-0 with the Karl Jansky VLA and at 100-500um using Herschel, alongside new and archival 12CO J=4-3 observations from IRAM PdBI. We introduce a new colour-colour diagnostic plot to constrain the redshifts of several distant, dusty galaxies in our target fields. A bright SMG near 6C1909+72 likely shares the same node or filament as the signpost AGN, but it is not detected in CO despite ~20,000 km/s of velocity coverage. Also in the 6C1909+72 field, a large, red dust feature spanning ~500 kpc is aligned with the radio jet. We suggest several processes by which metal-rich material may have been transported, favouring a collimated outflow reminiscent of the jet-oriented metal enrichment seen in local cluster environments. Our interferometric imaging reveals a gas-rich companion to B3J2330+3927; indeed, all bar one of the eight z >~ 2 radio galaxies (or companions) detected in CO provide some evidence that starburst activity in radio-loud AGN at high redshift is driven by the interaction of two or more gas-rich systems in which a significant mass of stars has already formed, rather than via steady accretion of cold gas from the cosmic web. We find that the CO Tb ratios in radio-loud AGN host galaxies are significantly higher than those seen in similarly intense starbursts where AGN activity is less pronounced. Our most extreme example, where L'(CO4-3)/L'(CO1-0) > 2.7, provides evidence that significant energy is being deposited rapidly into the molecular gas via X-rays and/or mechanical (`quasar-mode') feedback from the AGN, leading to a high degree of turbulence globally and a low optical depth in 12CO - feedback that may lead to the cessation of star formation on a timescale commensurate with that of the jet activity, <~10 Myr.