We present Atacama Large Millimeter/submillimeter Array observations of the brightest cluster galaxy Hydra-A, a nearby ($z=0.054$) giant elliptical galaxy with powerful and extended radio jets. The observations reveal CO(1-0), CO(2-1), $^{13}$CO(2-1), CN(2-1), SiO(5-4), HCO$^{+}$(1-0), HCO$^{+}$(2-1), HCN(1-0), HCN(2-1), HNC(1-0) and H$_{2}$CO(3-2) absorption lines against the galaxy's bright and compact active galactic nucleus. These absorption features are due to at least 12 individual molecular clouds which lie close to the centre of the galaxy and have velocities of approximately $-50$ to $+10$ km/s relative to its recession velocity, where positive values correspond to inward motion. The absorption profiles are evidence of a clumpy interstellar medium within brightest cluster galaxies composed of clouds with similar column densities, velocity dispersions and excitation temperatures to those found at radii of several kpc in the Milky Way. We also show potential variation in a $\sim 10$ km/s wide section of the absorption profile over a two year timescale, most likely caused by relativistic motions in the hot spots of the continuum source which change the background illumination of the absorbing clouds.

SURF is the collaborative ICT organisation for Dutch education & research and coordinates the national e-infrastructure (surf.nl). To accelerate scientific discovery, SURF invests in, operates and explores high-end IT solutions for and with researchers in the Netherlands. In this ADASS 2019 contribution we present our latest developments in high performance and high throughput cloud computing. These developments are particularly relevant for Astronomy, as this science domain: (i) generates large (Petabyte sized) data collections, (ii) uses rapid release and deployment schemes for their software, (iii) requires flexible and interactive test and staging environments and (iv) needs to execute complex workflows on diverse data structures. We highlight our new OpenStack-based, cloud infrastructure layer and focus on Spider, our new addition to SURF's high throughput data processing platforms.

We present low-radio-frequency follow-up observations of AT 2017gfo, the electromagnetic counterpart of GW170817, which was the first binary neutron star merger to be detected by Advanced LIGO-Virgo. These data, with a central frequency of 144 MHz, were obtained with LOFAR, the Low-Frequency Array. The maximum elevation of the target is just 13.7 degrees when observed with LOFAR, making our observations particularly challenging to calibrate and significantly limiting the achievable sensitivity. On time-scales of 130-138 and 371-374 days after the merger event, we obtain 3$\sigma$ upper limits for the afterglow component of 6.6 and 19.5 mJy beam$^{-1}$, respectively. Using our best upper limit and previously published, contemporaneous higher-frequency radio data, we place a limit on any potential steepening of the radio spectrum between 610 and 144 MHz: the two-point spectral index $\alpha^{610}_{144} \gtrsim -2.5$. We also show that LOFAR can detect the afterglows of future binary neutron star merger events occurring at more favourable elevations.

A radio interferometer uses time delays to maximize its response to radiation coming from a particular direction. These time delays compensate for differences in the time of arrival of the wavefront at the different elements of the interferometer, and for delays in the instrument's signal chain. If the radio interferometer is operated as a phased array (tied array), the time delays cannot be accounted for after an observation, so they must be determined in advance. Our aim is to characterize the time delays between the stations in the core of the LOw Frequency ARray (LOFAR). We used radio holography to determine the time delays for the core stations of LOFAR (innermost 3.5 km). Using the multibeaming capability of LOFAR we map the voltage beam faster than with a raster scan, while simultaneously calibrating the observed beam continuously. For short radio holographic observations (60 s and 600 s) of 3C196, 3C147, and 3C48 we are able to derive time delays with errors of less than one nanosecond. After applying the derived time delays to the beamformer, the beam shows residuals of less than $20\%$ with respect to the theoretical beam shape. Tied-array holography could be a way towards semi-real-time beam calibration for the Square Kilometer Array.

The new generation of high-resolution broad-band radio telescopes, like the Low Frequency Array (LOFAR), produces, depending on the level of compression, between 1 to 10 TB of data per hour after correlation. Such a large amount of scientific data demand powerful computing resources and efficient data handling strategies to be mastered. The LOFAR Two-metre Sky Survey (LoTSS) is a Key Science Project (KSP) of the LOFAR telescope. It aims to map the entire northern hemisphere at unprecedented sensitivity and resolution. The survey consist of 3 168 pointings, requiring about 30 PBytes of storage space. As a member of the German Long Wavelength Consortioum (GLOW) the Forschungszentrum J\"ulich (FSZ) stores in the Long Term Archive (LTA) about 50% of all LoTSS observations conducted to date. In collaboration with SURFsara in Amsterdam we developed service tools that enable the KSP to process LOFAR data stored in the J\"ulich LTA at the supercomputer JUWELS in an automated and robust fashion. Through our system more than 500 out of 800 existing LoTSS observations have already been processed with the prefactor pipeline. This pipeline calibrates the direction-independent instrumental and ionospheric effects and furthermore reduces the data size significantly. For continuum imaging, this processing pipeline is the standard pipeline that is executed before more advanced processing and image reconstruction methods are applied.

We present CO(1-0) and CO(3-2) ALMA observations of the molecular gas in PKS 1549-79, as well as mm and VLBI 2.3-GHz continuum observations of its radio jet. PKS 1549-79 is one of the closest young, radio-loud quasars caught in an on-going merger in which the AGN is in the first phases of its evolution. We detect three structures tracing the accretion and the outflow of molecular gas: kpc-scale tails of gas accreting onto PKS 1549-79, a circumnuclear disc (CND) in the inner few hundred parsec, and a very broad (>2300 \kms) component detected in CO(1-0) at the position of the AGN. Thus, in PKS 1549-79 we see the co-existence of accretion and the ejection of gas. The line ratio CO(1-0)/CO(3-2) suggests that the gas in the CND has both high densities and high kinetic temperatures. We estimate a mass outflow rate of at least 650 msun/yr. This massive outflow is confined to r < 120 pc, which suggests that the AGN drives the outflow. Considering the amount of molecular gas available in CND and the observed outflow rate, we estimate a time scale of ~10^5 yr over which the AGN would be able to destroy the CND, although gas from the merger may come in from larger radii, rebuilding this disc at the same time. The AGN appears to self-regulate gas accretion onto the super-massive black hole. From a comparison with HST data, we find that the ionised gas outflow is more extended. Nevertheless, the warm outflow is about two orders of magnitude less massive than the molecular outflow. PKS 1549-79 does not seem to follow the scaling relation between bolometric luminosity and the relative importance of warm ionised and molecular outflows claimed to exist for other AGN. We argue that, although PKS 1549-79 hosts a powerful quasar nucleus and an ultra-fast outflow, the radio jet plays a significant role in producing the outflow.

(abridged) Radio recombination lines (RRLs) at frequencies $\nu$ < 250 MHz trace the cold, diffuse phase of the ISM. Next generation low frequency interferometers, such as LOFAR, MWA and the future SKA, with unprecedented sensitivity, resolution, and large fractional bandwidths, are enabling the exploration of the extragalactic RRL universe. We observed the radio quasar 3C 190 (z~1.2) with the LOFAR HBA. In reducing this data for spectroscopic analysis, we have placed special emphasis on bandpass calibration. We devised cross-correlation techniques to significantly identify the presence of RRLs in a low frequency spectrum. We demonstrate the utility of this method by applying it to existing low-frequency spectra of Cassiopeia A and M 82, and to the new observations of 3C 190. RRLs have been detected in the foreground of 3C 190 at z = 1.12355 (assuming a carbon origin), owing to the first detection of RRLs outside of the local universe (first reported in Emig et al. 2019). Towards the Galactic supernova remnant Cas A, we uncover three new detections: (1) C$\epsilon$-transitions ($\Delta$n = 5) for the first time at low radio frequencies, (2) H$\alpha$-transitions at 64 MHz with a FWHM of 3.1 km/s, the most narrow and one of the lowest frequency detections of hydrogen to date, and (3) C$\alpha$ at v$_{LSR}$ = 0 km/s in the frequency range 55-78 MHz for the first time. Additionally we recover C$\alpha$, C$\beta$, C$\gamma$, and C$\delta$ from the -47 km/s and -38 km/s components. In the nearby starburst galaxy, M 82, we do not find a significant feature. Our current searches for RRLs in LOFAR observations are limited to narrow (< 100 km/s) features, owing to the relatively small number of channels available for continuum estimation. Future strategies making use of larger contiguous frequency coverage would aid calibration to deeper sensitivities and broader features.

To advance our understanding of the fuelling and feedback processes which power the Universe's most massive black holes, we require a significant increase in our knowledge of the molecular gas which exists in their immediate surroundings. However, the behaviour of this gas is poorly understood due to the difficulties associated with observing it directly. We report on a survey of 18 brightest cluster galaxies lying in cool cores, from which we detect molecular gas in the core regions of eight via carbon monoxide (CO), cyanide (CN) and silicon monoxide (SiO) absorption lines. These absorption lines are produced by cold molecular gas clouds which lie along the line of sight to the bright continuum sources at the galaxy centres. As such, they can be used to determine many properties of the molecular gas which may go on to fuel supermassive black hole accretion and AGN feedback mechanisms. The absorption regions detected have velocities ranging from -45 to 283 km s$^{-1}$ relative to the systemic velocity of the galaxy, and have a bias for motion towards the host supermassive black hole. We find that the CN N = 0 - 1 absorption lines are typically 10 times stronger than those of CO J = 0 - 1. This is due to the higher electric dipole moment of the CN molecule, which enhances its absorption strength. In terms of molecular number density CO remains the more prevalent molecule with a ratio of CO/CN $\sim 10$, similar to that of nearby galaxies. Comparison of CO, CN and HI observations for these systems shows many different combinations of these absorption lines being detected.

The Engineering Development Array (EDA) is a single test station for Square Kilometre Array (SKA) precursor technology. We have used the EDA to detect low-frequency radio recombination lines (RRLs) from the Galactic Centre region. Low-frequency RRLs are an area of interest for future low-frequency SKA work as these lines provide important information on the physical properties of the cold neutral medium. In this project we investigate the EDA, its bandpass and the radio frequency interference environment for low-frequency spectroscopy. We present line spectra from 30 to 325 MHz for the Galactic Centre region. The decrease in sensitivity for the EDA at the low end of the receiver prevents carbon and hydrogen RRLs to be detected below 40 and 60 MHz respectively. RFI strongly affects frequencies in the range 276-292, 234-270, 131-138, 95-102 and below 33 MHz. Cn$\alpha$ RRLs were detected in absorption for quantum levels n = 378 to 550 (39-121 MHz) and in emission for n = 272 to 306 (228-325 MHz). Cn$\beta$ lines were detected in absorption for n = 387 to 696 (39-225 MHz). Hn$\alpha$ RRLs were detected in emission for n = 272 to 480 (59-325 MHz). Hn$\beta$ lines were detected for n = 387 to 453 (141-225 MHz). The stacked Hn$\alpha$ detection at 63 MHz is the lowest frequency detection made for hydrogen RRLs and shows that a cold (partially) ionized medium exists along the line of sight to the Galactic Centre region. The size and velocity of this cold H$^{+}$ gas indicates that it is likely associated with the nearby Riegel-Crutcher cloud.

We present analytical expressions for direct evaluation of $\ell$-mixing rate coefficients in proton-excited hydrogen atom collisions and describe a software package for efficient numerical evaluation of the collisional rate coefficients. Comparisons between rate coefficients calculated with various levels of approximation are discussed, highlighting their range of validity. These rate coefficients are benchmarked via radio recombination lines for hydrogen, evaluating the corresponding departure coefficients from local thermal equilibrium.

Context. The combined use of carbon radio recombination lines (CRRLs) and the 158 $\mu$m-[CII] line is a powerful tool for the study of the energetics and physical conditions (e.g., temperature and density) of photodissociation regions (PDRs). However, there are few observational studies that exploit this synergy. Aims. Here we explore the relation between CRRLs and the 158 $\mu$m-[CII] line in light of new observations and models. Methods. We present new and existing observations of CRRLs in the frequency range 0.15--230 GHz with ALMA, VLA, the GBT, Effelsberg 100m, and LOFAR towards Orion~A (M42). We complement these observations with SOFIA observations of the 158 $\mu$m-[CII] line. We studied two PDRs: the foreground atomic gas, known as the Veil, and the dense PDR between the HII region and the background molecular cloud. Results. In the Veil we are able to determine the gas temperature and electron density, which we use to measure the ionization parameter and the photoelectric heating efficiency. In the dense PDR, we are able to identify a layered PDR structure at the surface of the molecular cloud to the south of the Trapezium cluster. There we find that the radio lines trace the colder portion of the ionized carbon layer, the C$^{+}$/C/CO interface. By modeling the emission of the $158$~$\mu$m-[CII] line and CRRLs as arising from a PDR we derive a thermal pressure $>5\times10^{7}$ K cm$^{-3}$ and a radiation field $G_{0}\approx10^{5}$ close to the Trapezium. Conclusions. This work provides additional observational support for the use of CRRLs and the 158 $\mu$m-[CII] line as complementary tools to study dense and diffuse PDRs, and highlights the usefulness of CRRLs as probes of the C$^{+}$/C/CO interface.

We present the first LOFAR observations of the radio jet in the quasar 4C+19.44 (a.k.a. PKS 1354+19) obtained with the long baselines. The achieved resolution is very well matched to that of archival Jansky Very Large Array (JVLA) observations at higher radio frequencies as well as the archival X-ray images obtained with {\it Chandra}. We found that, for several knots along the jet, the radio flux densities measured at hundreds of MHz lie well below the values estimated by extrapolating the GHz spectra. This clearly indicates the presence of spectral curvature. Radio spectral curvature has been already observed in different source classes and/or extended radio structures and it has been often interpreted as due to intrinsic processes, as a curved particle energy distribution, rather than absorption mechanisms ({ Razin-Tsytovich} effect, free-free or synchrotron self absorption to name a few). Here we discuss our results according to the scenario where particles undergo stochastic acceleration mechanisms also in quasar jet knots.

Active galactic nuclei play a crucial role in the accretion and ejection of gas in galaxies. Although their outflows are well studied, finding direct evidence of accretion has proved very difficult and has so far been done for very few sources. A promising way to study the significance of cold accretion is by observing the absorption of an active galactic nucleus's extremely bright radio emission by the cold gas lying along the line-of-sight. As such, we present ALMA CO(1-0) and CO(2-1) observations of the Hydra-A brightest cluster galaxy (z=0.054) which reveal the existence of cold, molecular gas clouds along the line-of-sight to the galaxy's extremely bright and compact mm-continuum source. They have apparent motions relative to the central supermassive black hole of between -43 and -4 km s$^{-1}$ and are most likely moving along stable, low ellipticity orbits. The identified clouds form part of a $\sim$$10^{9}$ $\text{M}_{\odot}$, approximately edge-on disc of cold molecular gas. With peak CO(2-1) optical depths of $\tau$=0.88 $^{+0.06}_{-0.06}$, they include the narrowest and by far the deepest absorption of this type which has been observed to date in a brightest cluster galaxy. By comparing the relative strengths of the lines for the most strongly absorbing region, we are able to estimate a gas temperature of $42^{+25}_{-11}$ K and line-of-sight column densities of $N_{CO}=2^{+3}_{-1}\times 10 ^{17} cm^{-2}$ and $N_{ H_{2} }=7^{+10}_{-4}\times 10 ^{20} cm^{-2}$.

The relative importance of the physical processes shaping the thermodynamics of the hot gas permeating rotating, massive early-type galaxies is expected to be different from that in non-rotating systems. Here, we report the results of the analysis of XMM-Newton data for the massive, lenticular galaxy NGC 7049. The galaxy harbours a dusty disc of cool gas and is surrounded by an extended hot X-ray emitting gaseous atmosphere with unusually high central entropy. The hot gas in the plane of rotation of the cool dusty disc has a multi-temperature structure, consistent with ongoing cooling. We conclude that the rotational support of the hot gas is likely capable of altering the multiphase condensation regardless of the $t_{\rm cool}/t_{\rm ff}$ ratio, which is here relatively high, $\sim 40$. However, the measured ratio of cooling time and eddy turnover time around unity ($C$-ratio $\approx 1$) implies significant condensation, and at the same time, the constrained ratio of rotational velocity and the velocity dispersion (turbulent Taylor number) ${\rm Ta_t} > 1$ indicates that the condensing gas should follow non-radial orbits forming a disc instead of filaments. This is in agreement with hydrodynamical simulations of massive rotating galaxies predicting a similarly extended multiphase disc.

Recombination lines involving high principal quantum numbers populate the radio spectrum in large numbers. Low-frequency (<1 GHz) observations of radio recombination lines (RRLs) primarily from carbon and hydrogen offer a new, if not unique, way to probe cold, largely atomic gas and warm, ionised gas in other galaxies. Furthermore, RRLs can be used to determine the physical state of the emitting regions, such as temperature and density. These properties make RRLs, potentially, a powerful tool of extragalactic ISM physics. At low radio frequencies, its conceivable to detect RRLs out to cosmological distances when illuminated by a strong radio continuum. However, they are extremely faint (tau ~ 1e-3 -- 1e-4) and have so far eluded detection outside of the local universe. With LOFAR observations of the radio quasar 3C 190 (z=1.1946), we aim to demonstrate that the ISM can be explored out to great distances through low-frequency RRLs. We report the detection of RRLs in the frequency range 112--163 MHz in the spectrum of 3C 190. Stacking 13 a-transitions with principal quantum numbers n=266-301, a peak 6sigma feature of optical depth, tau(peak) = (1.0 +\- 0.2) x 1e-3 and FWHM = 31.2 +/- 8.3 km/s was found at z=1.124. This corresponds to a velocity offset of -9965 km/s with respect to the systemic redshift of 3C 190. We consider three interpretations of the origin of the RRL emission: an intervening dwarf-like galaxy, an AGN-driven outflow, and the inter-galactic medium. We argue that the RRLs most likely originate in a dwarf-like galaxy (M ~ 1e9 Msun) along the line of sight, although we cannot rule out an AGN-driven outflow. We do find the RRLs to be inconsistent with an inter-galactic medium origin. With this detection, we have opened up a new way to study the physical properties of cool, diffuse gas out to cosmological distances.

The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168MHz survey of the entire northern sky for which observations are now 20% complete. We present our first full-quality public data release. For this data release 424 square degrees, or 2% of the eventual coverage, in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45$^\circ$00$'$00$''$ to 57$^\circ$00$'$00$''$) were mapped using a fully automated direction-dependent calibration and imaging pipeline that we developed. A total of 325,694 sources are detected with a signal of at least five times the noise, and the source density is a factor of $\sim 10$ higher than the most sensitive existing very wide-area radio-continuum surveys. The median sensitivity is S$_{\rm 144 MHz} = 71\,\mu$Jy beam$^{-1}$ and the point-source completeness is 90% at an integrated flux density of 0.45mJy. The resolution of the images is 6$''$ and the positional accuracy is within 0.2$''$. This data release consists of a catalogue containing location, flux, and shape estimates together with 58 mosaic images that cover the catalogued area. In this paper we provide an overview of the data release with a focus on the processing of the LOFAR data and the characteristics of the resulting images. In two accompanying papers we provide the radio source associations and deblending and, where possible, the optical identifications of the radio sources together with the photometric redshifts and properties of the host galaxies. These data release papers are published together with a further $\sim$20 articles that highlight the scientific potential of LoTSS.

Mixed-morphology supernova remnants (MM SNRs) are a mysterious class of objects that display thermal X-ray emission within their radio shell. They are an older class of SNRs, and as such are profoundly affected by the environment into which they evolve. VRO 42.05.01 is a MM SNR of puzzling morphology in the direction of the Galactic anticentre. Low-frequency radio observations of supernova remnants are sensitive to synchrotron electrons accelerated in the shock front. We aim to compare the low-frequency emission to higher frequency observations to understand the environmental and shock acceleration conditions that have given rise to the observed properties of this source. We present a LOFAR High Band Antenna map centred at 143 MHz of the region of the Galactic plane centred at $l = 166^\rm{o}, b = 3.5^\rm{o}$ at 143 MHz, with a resolution of 148'' and an rms noise of 4.4 mJy bm$^{-1}$ . Our map is sensitive to scales as large as $6^\rm{o}$. We compared the LOw Frequency ARay (LOFAR) observations to archival higher frequency radio, infrared, and optical data to study the emission properties of the source in different spectral regimes. We did this both for the SNR and for OA 184, an H II region within our field of view. We find that the radio spectral index of VRO 42.05.01 increases at low radio frequencies; i.e. the LOFAR flux is higher than expected from the measured spectral index value at higher radio frequencies. This observed curvature in the low-frequency end of the radio spectrum occurs primarily in the brightest regions of the source, while the fainter regions present a roughly constant power-law behaviour between 143 MHz and 2695 MHz. We favour an explanation for this steepening whereby radiative shocks have high compression ratios and electrons of different energies probe different length scales across the shocks, therefore sampling regions of different compression ratios.

We present new, deep (245 ks) Chandra observations of the galaxy cluster Abell 1664 ($z = 0.1283$). These images reveal rich structure, including elongation and accompanying compressions of the X-ray isophotes in the NE-SW direction, suggesting that the hot gas is sloshing in the gravitational potential. This sloshing has resulted in cold fronts, at distances of 55, 115 and 320 kpc from the cluster center. Our results indicate that the core of A1664 is highly disturbed, as the global metallicity and cooling time flatten at small radii, implying mixing on large scales. The central AGN appears to have recently undergone a mechanical outburst, as evidenced by our detection of cavities. These cavities are the X-ray manifestations of radio bubbles inflated by the AGN, and may explain the motion of cold molecular CO clouds previously observed with ALMA. The estimated mechanical power of the AGN, using the minimum energy required to inflate the cavities as a proxy, is $P_{\rm cav} = (1.1 \pm 1.0) \times 10^{44} $ erg s$^{-1}$, which may be enough to drive the molecular gas flows, and offset the cooling luminosity of the ICM, at $L_{\rm cool} = (1.90 \pm0.01)\times 10^{44}$ erg s$^{-1}$. This mechanical power is orders of magnitude higher than the measured upper limit on the X-ray luminosity of the central AGN, suggesting that its black hole may be extremely massive and/or radiatively inefficient. We map temperature variations on the same spatial scale as the molecular gas, and find that the most rapidly cooling gas is mostly coincident with the molecular gas reservoir centered on the BCG's systemic velocity observed with ALMA and may be fueling cold accretion onto the central black hole.

The LOFAR radio telescope creates Petabytes of data per year. This data is important for many scientific projects. The data needs to be efficiently processed within the timespan of these projects in order to maximize the scientific impact. We present a workflow orchestration system that integrates LOFAR processing with a distributed computing platform. The system is named Automated Grid-enabled LOFAR Workflows (AGLOW). AGLOW makes it fast and easy to develop, test and deploy complex LOFAR workflows, and to accelerate them on a distributed cluster architecture. AGLOW provides a significant reduction in time for setting up complex workflows: typically, from months to days. We lay out two case studies that process the data from the LOFAR Surveys Key Science Project. We have implemented these into the AGLOW environment. We also describe the capabilities of AGLOW, paving the way for use by other LOFAR science cases. In the future, AGLOW will automatically produce multiple science products from a single dataset, serving several of the LOFAR Key Science Projects.

We present ALMA and MUSE observations of the Brightest Cluster Galaxy in Abell 2597, a nearby (z=0.0821) cool core cluster of galaxies. The data map the kinematics of a three billion solar mass filamentary nebula that spans the innermost 30 kpc of the galaxy's core. Its warm ionized and cold molecular components are both cospatial and comoving, consistent with the hypothesis that the optical nebula traces the warm envelopes of many cold molecular clouds that drift in the velocity field of the hot X-ray atmosphere. The clouds are not in dynamical equilibrium, and instead show evidence for inflow toward the central supermassive black hole, outflow along the jets it launches, and uplift by the buoyant hot bubbles those jets inflate. The entire scenario is therefore consistent with a galaxy-spanning "fountain", wherein cold gas clouds drain into the black hole accretion reservoir, powering jets and bubbles that uplift a cooling plume of low-entropy multiphase gas, which may stimulate additional cooling and accretion as part of a self-regulating feedback loop. All velocities are below the escape speed from the galaxy, and so these clouds should rain back toward the galaxy center from which they came, keeping the fountain long-lived. The data are consistent with major predictions of chaotic cold accretion, precipitation, and stimulated feedback models, and may trace processes fundamental to galaxy evolution at effectively all mass scales.

Our view of the central regions of AGN has been enriched by the discovery of fast and massive outflows of HI and molecular gas. Here we present a brief summary of results obtained for young (and restarted) radio AGN. We find that HI outflows tend to be particularly common in this group of objects. This supports the idea that the jet, expanding in a clumpy medium, plays a major role in driving these outflows. The clumpiness of the medium is confirmed by VLBI and ALMA observations. The HI observations reveal that, at least part of the gas, is distributed in clouds with sizes up to a few tens of pc and mass ~10^4 Msun. A change of the conditions in the outflow, with an increasing fraction of diffuse components, as the radio jets grow, is suggested by the high resolution HI observations. The molecular gas completes the picture, showing that the radio plasma jet can couple well with the ISM, strongly affecting the kinematics, but also the physical conditions of the molecular gas. This is confirmed by numerical simulations reproducing, to first order, the kinematics of the gas.

Modern astronomical data processing requires complex software pipelines to process ever growing datasets. For radio astronomy, these pipelines have become so large that they need to be distributed across a computational cluster. This makes it difficult to monitor the performance of each pipeline step. To gain insight into the performance of each step, a performance monitoring utility needs to be integrated with the pipeline execution. In this work we have developed such a utility and integrated it with the calibration pipeline of the Low Frequency Array, LOFAR, a leading radio telescope. We tested the tool by running the pipeline on several different compute platforms and collected the performance data. Based on this data, we make well informed recommendations on future hardware and software upgrades. The aim of these upgrades is to accelerate the slowest processing steps for this LOFAR pipeline. The pipeline collector suite is open source and will be incorporated in future LOFAR pipelines to create a performance database for all LOFAR processing.

The Low Frequency Array (LOFAR) radio telescope is an international aperture synthesis radio telescope used to study the Universe at low frequencies. One of the goals of the LOFAR telescope is to conduct deep wide-field surveys. Here we will discuss a framework for the processing of the LOFAR Two Meter Sky Survey (LoTSS). This survey will produce close to 50 PB of data within five years. These data rates require processing at locations with high-speed access to the archived data. To complete the LoTSS project, the processing software needs to be made portable and moved to clusters with a high bandwidth connection to the data archive. This work presents a framework that makes the LOFAR software portable, and is used to scale out LOFAR data reduction. Previous work was successful in preprocessing LOFAR data on a cluster of isolated nodes. This framework builds upon it and and is currently operational. It is designed to be portable, scalable, automated and general. This paper describes its design and high level operation and the initial results processing LoTSS data.

Cassiopeia A is one of the best-studied supernova remnants. Its shocked ejecta emits brightly in radio and X-rays. Its unshocked ejecta can be studied through infrared emission, the radio-active decay of $^{44}$Ti, and low frequency free-free absorption due to cold gas internal to the shell. Free-free absorption is affected by the mass, geometry, temperature, and ionisation conditions in the absorbing gas. Observations at the lowest radio frequencies constrain a combination of these properties. We use LOFAR LBA observations at 30-77 MHz and L-band VLA observations to compare $u-v$-matched images with a common resolution of 17". We simultaneously fit, per pixel, for the emission measure and the ratio of the emission from the unabsorbed front of the shell versus the absorbed back of the shell. We explore the effects that low temperatures and a high degree of clumping can have on the derived physical properties, such as mass and density. We also compile published radio flux measurements, fit for the absorption processes that occur in the radio band, and consider how they affect the secular decline of the source. We find a mass in the unshocked ejecta of $M = 2.95 \pm {0.48} \,M_{\odot}$ for an assumed gas temperature of $T=100$ K. This estimate is reduced for colder gas temperatures and if the ejecta are clumped. We measure the reverse shock to have a radius of $114$" $\pm $6". We also find that a decrease in the amount of mass in the unshocked ejecta (as more and more material meets the reverse shock and heats up) cannot account for the observed low frequency behaviour of the secular decline rate. To reconcile our low frequency absorption measurements with models that predict little mass in the unshocked ejecta we need the ejecta to be very clumped, or the temperature in the cold gas to be low ($\sim10$ K). Both conditions can jointly contribute to the high absorption.

Quantitative understanding of the interstellar medium requires knowledge of its physical conditions. Low frequency carbon radio recombination lines (CRRLs) trace cold interstellar gas, and can be used to determine its physical conditions (e.g., electron temperature and density). In this work we present spatially resolved observations of the low frequency ($\leq390$ MHz) CRRLs centered around C$268\alpha$, C$357\alpha$, C$494\alpha$ and C$539\alpha$ towards Cassiopeia A on scales of $\leq1.2$ pc. We compare the spatial distribution of CRRLs with other ISM tracers. This comparison reveals a spatial offset between the peak of the CRRLs and other tracers, which is very characteristic for photodissociation regions and that we take as evidence for CRRLs being preferentially detected from the surfaces of molecular clouds. Using the CRRLs we constrain the gas electron temperature and density. These constraints on the gas conditions suggest variations of less than a factor of two in pressure over $\sim1$ pc scales, and an average hydrogen density of $200$-$470$ cm$^{-3}$. From the electron temperature and density maps we also constrain the ionized carbon emission measure, column density and path length. Based on these, the hydrogen column density is larger than $10^{22}$ cm$^{-2}$, with a peak of $\sim4\times10^{22}$ cm$^{-2}$ towards the South of Cassiopeia A. Towards the southern peak the line of sight length is $\sim40$ pc over a $\sim2$ pc wide structure, which implies that the gas is a thin surface layer on a large (molecular) cloud that is only partially intersected by Cassiopeia A. These observations highlight the utility of CRRLs as tracers of low density extended HI and CO-dark gas halo's around molecular clouds.

We present ALMA observations of the $^{12}$CO (2--1) line of the newly born ($t_\mathrm{radio}\sim10^2$ years) active galactic nucleus (AGN), PKS B1718-649. These observations reveal that the carbon monoxide in the innermost 15 kpc of the galaxy is distributed in a complex warped disk. In the outer parts of this disk, the CO gas follows the rotation of the dust lane and of the stellar body of the galaxy hosting the radio source. In the innermost kiloparsec, the gas abruptly changes orientation and forms a circumnuclear disk ($r\lesssim700$ pc) with its major axis perpendicular to that of the outer disk. Against the compact radio emission of PKS B1718-649 ($r\sim 2$ pc), we detect an absorption line at red-shifted velocities with respect to the systemic velocity ($\Delta v = +365\pm22$\kms). This absorbing CO gas could trace molecular clouds falling onto the central super-massive black hole. A comparison with the near-infra red H$_{\,2}$ 1-0 S(1) observations shows that the clouds must be close to the black hole ($r\lesssim 75$ pc). The physical conditions of these clouds are different from the gas at larger radii, and are in good agreement with the predictions for the conditions of the gas when cold chaotic accretion triggers an active galactic nucleus. These observations on the centre of PKS B1718-649 provide one of the best indications that a population of cold clouds is falling towards a radio AGN, likely fuelling its activity.

We present the results from ALMA observations centred $40^{\prime\prime}$ (3 kpc in projection) southeast of the nucleus of M87. We report the detection of extended CO (2-1) line emission with a total flux of $(5.5 \pm 0.6) \times 10^{-18}$ erg s$^{-1}$ cm$^{-2}$ and corresponding molecular gas mass $M_{H_2}=(4.7 \pm 0.4) \times 10^5 M_\odot$, assuming a Galactic CO to H$_2$ conversion factor. ALMA data indicate a line-of-sight velocity of $-129\pm3$ km s$^{-1}$, in good agreement with measurements based on the [CII] and H$\alpha$+[NII] lines, and a velocity dispersion of $\sigma=27\pm3$ km s$^{-1}$. The CO(2-1) emission originates only outside the radio lobe of the AGN seen in the 6~cm VLA image, while the filament prolongs further inwards at other wavelengths. The molecular gas in M87 appears to be destroyed or excited by AGN activity, either by direct interaction with the radio plasma, or by the shock driven by the lobe into the X-ray emitting atmosphere. This is an important piece of the puzzle in understanding the impact of the central AGN on the amount of the coldest gas from which star formation can proceed.

We present a detailed study of the molecular gas in the fast AGN-driven outflow in the nearby radio-loud Seyfert galaxy IC 5063. Using ALMA observations of a number of tracers (12CO(1-0), 12CO(2-1), 12CO(3-2), 13CO(2-1) and HCO+(4-3)), we map the differences in excitation, density and temperature of the gas. The results show that in the immediate vicinity of the radio jet, a fast outflow, with velocities up to 800 km/s, is occurring of which the gas has high excitation temperatures in the range 30-55 K, demonstrating the direct impact of the jet on the ISM. The relative brightness of the CO lines show that the outflow is optically thin. We estimate the mass of the molecular outflow to be 1.2 x 10^6 Msol and likely to be a factor 2-3 larger. This is similar to that of the outflow of atomic gas, but much larger than that of the ionised outflow, showing that the outflow is dominated by cold gas. The total mass outflow rate we estimate to be ~12 Msol/yr. The mass of the outflow is much smaller than the total gas mass of the ISM of IC 5063. Therefore, although the influence of the radio jet is very significant in the inner regions, globally speaking the impact will be very modest. We use RADEX modelling to explore the physical conditions of the molecular gas in the outflow. Models with the outflowing gas being quite clumpy give the most consistent results and our preferred solutions have kinetic temperatures in the range 20-100 K and densities between 10^5 and 10^6 cm^-3. The resulting pressures are 10^6-10^7.5 K cm^-3, about two orders of magnitude higher than in the outer quiescent disk. The results strongly suggest that the outflow is driven by the radio jet expanding into a clumpy medium, creating a cocoon of gas which is pushed away from the jet axis resulting in a lateral outflow, very similar to what is predicted by numerical simulations.

We present a study of carbon radio recombination lines towards Cassiopeia A using LOFAR observations in the frequency range 10-33 MHz. Individual carbon $\alpha$ lines are detected in absorption against the continuum at frequencies as low as 16 MHz. Stacking several C$\alpha$ lines we obtain detections in the 11-16 MHz range. These are the highest signal-to-noise measurements at these frequencies. The peak optical depth of the C$\alpha$ lines changes considerably over the 11-33 MHz range with the peak optical depth decreasing from 4$\times10^{-3}$ at 33 MHz to 2$\times10^{-3}$ at 11 MHz, while the line width increases from 20 km s$^{-1}$ to 150 km s$^{-1}$. The combined change in peak optical depth and line width results in a roughly constant integrated optical depth. We interpret this as carbon atoms close to local thermodynamic equilibrium. In this work we focus on how the 11-33 MHz carbon radio recombination lines can be used to determine the gas physical conditions. We find that the ratio of the carbon radio recombination lines to that of the 158 $\mu$m [CII] fine-structure line is a good thermometer, while the ratio between low frequency carbon radio recombination lines provides a good barometer. By combining the temperature and pressure constraints with those derived from the line width we are able to constrain the gas properties (electron temperature and density) and radiation field intensity. Given the 1$\sigma$ uncertainties in our measurements these are; $T_{e}\approx68$-$98$ K, $n_{e}\approx0.02$-$0.035$ cm$^{-3}$ and $T_{r,100}\approx1500$-$1650$ K. Despite challenging RFI and ionospheric conditions, our work demonstrates that observations of carbon radio recombination lines in the 10-33 MHz range can provide insight into the gas conditions.

According to optical stellar kinematics observations, an over-massive black hole candidate has been reported by van den Bosch et al. (2012) in the normal early-type galaxy NGC 1277. This galaxy is located in the central region of the Perseus cluster. Westerbork Synthesis Radio Telescope (WSRT) observations have shown that NGC 1277 and other early-type galaxies in the neighbourhood have radio counterparts. These nuclear radio sources have stable flux densities on time scale of years. In order to investigate the origin of the radio emission from these normal galaxies, we selected five sources (NGC 1270, NGC 1272, NGC 1277, NGC 1278 and VZw 339) residing in the central 10 arcminute region of the Perseus cluster and requested to re-correlate the data of an existing very long baseline interferometry (VLBI) experiment at these new positions. With the re-correlation data provided by the European VLBI Network (EVN), we imaged the five sources with a resolution of about eight milliarcseconds and detected all of them with a confidence level above 5{\sigma} at 1.4 GHz. They show compact structure and brightness temperatures above $10^7$ K, which implies that the radio emission is non-thermal. We rule out ongoing nuclear star formation and conclude that these VLBI-detected radio sources are parsec-scale jet activity associated with the supermassive black holes in low luminosity active galactic nuclei (LLAGNs), although there are no clear signs of nuclear activity observed in the optical and infrared bands. Using the fundamental plane relation in black holes, we find no significant evidence for or against an extremely massive black hole hiding in NGC 1277.

The LOFAR Two-metre Sky Survey (LoTSS) is a deep 120-168 MHz imaging survey that will eventually cover the entire Northern sky. Each of the 3170 pointings will be observed for 8 hrs, which, at most declinations, is sufficient to produce ~5arcsec resolution images with a sensitivity of ~0.1mJy/beam and accomplish the main scientific aims of the survey which are to explore the formation and evolution of massive black holes, galaxies, clusters of galaxies and large-scale structure. Due to the compact core and long baselines of LOFAR, the images provide excellent sensitivity to both highly extended and compact emission. For legacy value, the data are archived at high spectral and time resolution to facilitate subarcsecond imaging and spectral line studies. In this paper we provide an overview of the LoTSS. We outline the survey strategy, the observational status, the current calibration techniques, a preliminary data release, and the anticipated scientific impact. The preliminary images that we have released were created using a fully-automated but direction-independent calibration strategy and are significantly more sensitive than those produced by any existing large-area low-frequency survey. In excess of 44,000 sources are detected in the images that have a resolution of 25arcsec, typical noise levels of less than 0.5 mJy/beam, and cover an area of over 350 square degrees in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45d00m00s to 57d00m00s).

We use the Low Frequency Array to perform a systematic high spectral resolution investigation of the low-frequency 33-78 MHz spectrum along the line of sight to Cassiopeia A. We complement this with a 304-386 MHz Westerbork Synthesis Radio telescope observation. In this first paper we focus on the carbon radio recombination lines. We detect Cn$\alpha$ lines at -47 and -38 km s$^{-1}$ in absorption for quantum numbers n=438-584 and in emission for n=257-278 with high signal to noise. These lines are associated with cold clouds in the Perseus spiral arm component. Hn$\alpha$ lines are detected in emission for n=257-278. In addition, we also detect Cn$\alpha$ lines at 0 km s$^{-1}$ associated with the Orion arm. We analyze the optical depth of these transitions and their line width. Our models show that the carbon line components in the Perseus arm are best fit with an electron temperature 85 K and an electron density 0.04 cm$^{-3}$ and can be constrained to within 15\%. The electron pressure is constrained to within 20\%. We argue that much of these carbon radio recombination lines arise in the CO-dark surface layers of molecular clouds where most of the carbon is ionized but hydrogen has made the transition from atomic to molecular. The hydrogen lines are clearly associated with the carbon line emitting clouds, but the low-frequency upperlimits indicate that they likely do not trace the same gas. Combining the hydrogen and carbon results we arrive at a firm lower limit to the cosmic ray ionization rate of 2.5$\times$10$^{-18}$ s$^{-1}$, but the actual value is likely much larger.

Using CO (4-3) and (2-1) Atacama Large Millimeter Array (ALMA) data, we prove that the molecular gas in the jet-driven winds of the galaxy IC5063 is more highly excited than the rest of the molecular gas in the disk of the same galaxy. On average, the CO(4-3)/CO(2-1) flux ratio is 1 for the disk and 5 for the jet accelerated or impacted gas. Spatially-resolved maps reveal that in regions associated with winds, the CO(4-3)/CO(2-1) flux ratio significantly exceeds the upper limit of 4 for optically thick gas. It frequently takes values between 5 and 11, and it occasionally further approaches the upper limit of 16 for optically thin gas. Excitation temperatures of 30-100 K are common for the molecules in these regions. If all of the outflowing molecular gas is optically thin, at 30-50 K, then its mass is 2*10^6 M_sun. This lower mass limit is an order of magnitude below the mass derived from the CO(2-1) flux in the case of optically thick gas. Molecular winds can thus be less massive, but more easily detectable at high z than they were previously thought to be.

We present ALMA CO(1-0) and CO(3-2) observations of the brightest cluster galaxy (BCG) in the 2A 0335+096 galaxy cluster (z = 0.0346). The total molecular gas mass of (1.13+/-0.15) x 10^9 M_sun is divided into two components: a nuclear region and a 7 kpc long dusty filament. The central molecular gas component accounts for (3.2+/-0.4) x 10^8 M_sun of the total supply of cold gas. Instead of forming a rotationally-supported ring or disk, it is composed of two distinct, blueshifted clumps south of the nucleus and a series of low-significance redshifted clumps extending toward a nearby companion galaxy. The velocity of the redshifted clouds increases with radius to a value consistent with the companion galaxy, suggesting that an interaction between these galaxies <20 Myr ago disrupted a pre-existing molecular gas reservoir within the BCG. Most of the molecular gas, (7.8+/-0.9) x 10^8 M_sun, is located in the filament. The CO emission is co-spatial with a 10^4 K emission-line nebula and soft X-rays from 0.5 keV gas, indicating that the molecular gas has cooled out of the intracluster medium over a period of 25-100 Myr. The filament trails an X-ray cavity, suggesting that the gas has cooled from low entropy gas that has been lifted out of the cluster core and become thermally unstable. We are unable to distinguish between inflow and outflow along the filament with the present data. Cloud velocities along the filament are consistent with gravitational free-fall near the plane of the sky, although their increasing blueshifts with radius are consistent with outflow.

In the second paper of the series, we have modeled low frequency carbon radio recombination lines (CRRL) from the interstellar medium. Anticipating the LOw Frequency ARray (LOFAR) survey of Galactic CRRLs, we focus our study on the physical conditions of the diffuse cold neutral medium (CNM). We have used the improved departure coefficients computed in the first paper of the series to calculate line-to-continuum ratios. The results show that the line width and integrated optical depths of CRRL are sensitive probes of the electron density, gas temperature, and the emission measure of the cloud. Furthermore, the ratio of CRRL to the [CII] at 158 $\mu$m line is a strong function of the temperature and density of diffuse clouds. Guided by our calculations, we analyze CRRL observations and illustrate their use with data from the literature.

In the first paper of this series, we study the level population problem of recombining carbon ions. We focus our study on high quantum numbers anticipating observations of Carbon Radio Recombination Lines to be carried out by the LOw Frequency ARray (LOFAR). We solve the level population equation including angular momentum levels with updated collision rates up to high principal quantum numbers. We derive departure coefficients by solving the level population equation in the hydrogenic approximation and including low temperature dielectronic recombination effects. Our results in the hydrogenic approximation agree well with those of previous works. When comparing our results including dielectronic recombination we find differences which we ascribe to updates in the atomic physics (e.g., collision rates) and to the approximate solution method of the statistical equilibrium equations adopted in previous studies. A comparison with observations is discussed in an accompanying article, as radiative transfer effects need to be considered.

We report the detection, using observations of the CO(2-1) line performed with the Atacama Pathfinder EXperiment (APEX), of molecular gas in the region of the outer filament of Centaurus A, a complex region known to show various signatures of an interaction between the radio jet, an HI cloud, and ionised gas filaments. We detect CO(2-1) at all observed locations, which were selected to represent regions with very different physical conditions. The H_2 masses of the detections range between 0.2 x 10^6 and 1.1 x 10^6 \msun, for conservative choices of the CO to H_2 conversion factor. Surprisingly, the stronger detections are not coincident with the HI cloud, but instead are in the region of the ionised filaments. We also find variations in the widths of the CO(2-1) lines throughout the region, with broader lines in the region of the ionised gas, i.e. where the jet--cloud interaction is strongest, and with narrow profiles in the HI cloud. This may indicate that the molecular gas in the region of the ionised gas has the momentum of the jet-cloud interaction encoded in it, in the same way as the ionised gas does. These molecular clouds may therefore be the result of very efficient cooling of the down-stream gas photo- or shock-ionised by the interaction. On the other hand, the molecular clouds with narrower profiles, which are closer to or inside the HI cloud, could be pre-existing cold H_2 cores which manage to survive the effects of the passing jet.

The correlation between radio spectral index and redshift has been exploited to discover high redshift radio galaxies, but its underlying cause is unclear. It is crucial to characterise the particle acceleration and loss mechanisms in high redshift radio galaxies to understand why their radio spectral indices are steeper than their local counterparts. Low frequency information on scales of $\sim$1 arcsec are necessary to determine the internal spectral index variation. In this paper we present the first spatially resolved studies at frequencies below 100 MHz of the $z = 2.4$ radio galaxy 4C 43.15 which was selected based on its ultra-steep spectral index ($\alpha < -1$; $S_{\nu} \sim \nu^{\alpha}$ ) between 365 MHz and 1.4 GHz. Using the International Low Frequency Array (LOFAR) Low Band Antenna we achieve sub-arcsecond imaging resolution at 55 MHz with VLBI techniques. Our study reveals low-frequency radio emission extended along the jet axis, which connects the two lobes. The integrated spectral index for frequencies $<$ 500 MHz is -0.83. The lobes have integrated spectral indices of -1.31$\pm$0.03 and -1.75$\pm$0.01 for frequencies $\geq$1.4 GHz, implying a break frequency between 500 MHz and 1.4 GHz. These spectral properties are similar to those of local radio galaxies. We conclude that the initially measured ultra-steep spectral index is due to a combination of the steepening spectrum at high frequencies with a break at intermediate frequencies.

Supermassive black holes in galaxy centres can grow by the accretion of gas, liberating energy that might regulate star formation on galaxy-wide scales. The nature of the gaseous fuel reservoirs that power black hole growth is nevertheless largely unconstrained by observations, and is instead routinely simplified as a smooth, spherical inflow of very hot gas. Recent theory and simulations instead predict that accretion can be dominated by a stochastic, clumpy distribution of very cold molecular clouds - a departure from the "hot mode" accretion model - although unambiguous observational support for this prediction remains elusive. Here we report observations that reveal a cold, clumpy accretion flow towards a supermassive black hole fuel reservoir in the nucleus of the Abell 2597 Brightest Cluster Galaxy (BCG), a nearby (redshift z=0.0821) giant elliptical galaxy surrounded by a dense halo of hot plasma. Under the right conditions, thermal instabilities can precipitate from this hot gas, producing a rain of cold clouds that fall toward the galaxy's centre, sustaining star formation amid a kiloparsec-scale molecular nebula that inhabits its core. The observations show that these cold clouds also fuel black hole accretion, revealing "shadows" cast by the molecular clouds as they move inward at about 300 kilometres per second towards the active supermassive black hole in the galaxy centre, which serves as a bright backlight. Corroborating evidence from prior observations of warmer atomic gas at extremely high spatial resolution, along with simple arguments based on geometry and probability, indicate that these clouds are within the innermost hundred parsecs of the black hole, and falling closer towards it.

We present a detailed study of the complex ionization structure in a small (~250 pc) extended narrow line region (ENLR) cloud near Centaurus A using the Multi Unit Spectroscopic Explorer. This cloud is located in the so-called outer filament of ionized gas (about 15 kpc from the nucleus) where jet-induced star formation has been suggested to occur by different studies. We find that, despite the small size, a mixture of ionization mechanisms is operating, resulting in considerable complexity in the spatial ionization structure. The area includes two H II regions where star formation is occurring and another location where star formation must have ceased very recently. Interestingly, the extreme Balmer decrement of one of the star forming regions (H_alpha/H_beta~6) indicates that it is still heavily embedded in its natal cocoon of gas and dust. At all three locations a continuum counterpart is found with spectra matching those of O/B stars local to Centaurus A. The H II regions are embedded in a larger gas complex which is photoionized by the radiation of the central active galactic nucleus (AGN), but the O/B stars affect the spatial ionization pattern in the ENLR cloud very locally. In particular, in the surroundings of the youngest star forming region, we can isolate a tight mixing sequence in the diagnostic diagram going from gas with ionization due to a pure stellar continuum to gas only photoionized by the AGN. These results emphasize the complexity and the mixture of processes occurring in star forming regions under the influence of an AGN radiation. This is relevant for our understanding of AGN-induced star formation suggested to occur in a number of objects, including this region of Centaurus A. They also illustrate that these young stars influence the gas over only a limited region.

We present ALMA observations of the CO(1-0) and CO(3-2) line emission tracing filaments of cold molecular gas in the central galaxy of the cluster PKS0745-191. The total molecular gas mass of 4.6 +/- 0.3 x 10^9 solar masses, assuming a Galactic X_{CO} factor, is divided roughly equally between three filaments each extending radially 3-5 kpc from the galaxy centre. The emission peak is located in the SE filament roughly 1 arcsec (2 kpc) from the nucleus. The velocities of the molecular clouds in the filaments are low, lying within +/-100 km/s of the galaxy's systemic velocity. Their FWHMs are less than 150 km/s, which is significantly below the stellar velocity dispersion. Although the molecular mass of each filament is comparable to a rich spiral galaxy, such low velocities show that the filaments are transient and the clouds would disperse on <10^7 yr timescales unless supported, likely by the indirect effect of magnetic fields. The velocity structure is inconsistent with a merger origin or gravitational free-fall of cooling gas in this massive central galaxy. If the molecular clouds originated in gas cooling even a few kpc from their current locations their velocities would exceed those observed. Instead, the projection of the N and SE filaments underneath X-ray cavities suggests they formed in the updraft behind bubbles buoyantly rising through the cluster atmosphere. Direct uplift of the dense gas by the radio bubbles appears to require an implausibly high coupling efficiency. The filaments are coincident with low temperature X-ray gas, bright optical line emission and dust lanes indicating that the molecular gas could have formed from lifted warmer gas that cooled in situ.

We present new SINFONI VLT observations of molecular hydrogen (H2) in the central regions (< 2.5 kpc) of the youngest and closest radio source PKS B1718-649. We study the distribution of the H2 traced by the 1-0 S(1) ro-vibrational line, revealing a double disk structure with the kinematics of both disks characterised by rotation. An outer disk (r > 650 pc) is aligned with other components of the galaxy (atomic hydrogen, stars, dust), while the inner disk (r< 600 pc) is perpendicular to it and is polar with respect to the stellar distribution. However, in the innermost 75 pc, the data show the presence of H2 gas redshifted with respect to the rotating inner disk ($\Delta v\,$+150 km/s) which may trace gas falling into the super-massive black hole associated with the central radio source. Along the same line of sight, earlier observations had shown the presence in the central regions of PKS B1718-649 of clouds of atomic hydrogen with similar unsettled kinematics. The range of velocities and mass of these unsettled clouds of HI and H2 suggest they may be actively contributing in fuelling the central newly-born radio source.

The interaction between the radio plasma ejected by the active nucleus of a galaxy and the surrounding medium is a key process that can have a strong impact on the interstellar medium of the galaxy and hence on galaxy evolution. The closest laboratory where we can observe and investigate this phenomenon is the radio galaxy Centaurus A. About 15 kpc north-east of this galaxy, a particularly complex region is found: the so-called Outer Filament where jet-cloud interactions have been proposed to occur. We investigate the presence of signatures of jet-ISM interaction by a detailed study of the kinematics of the ionized gas, expanding on previous results obtained from the HI. We observed two regions of the outer filament with VLT/VIMOS in the IFU observing mode. Emission from Hbeta and [OIII]4959,5007\AA\ is detected in both pointings. We found two distinct kinematical components of ionized gas that well match the kinematics of the nearby HI cloud. One component follows the regular kinematics of the rotating gas while the second shows similar velocities to those of the nearby HI component thought to be disturbed by an interaction with the radio jet. We suggest that the ionized and atomic gas are part of the same dynamical gas structure originating as result of the merger that shaped Centaurus A and which is regularly rotating around Centaurus A as proposed by other authors. The gas (ionized and HI) with anomalous velocities is tracing the interaction of the Large-Scale radio Jet with the ISM, suggesting that, although poorly collimated as structure, the jet is still active. However, we can exclude that a strong shock is driving the ionization of the gas. It is likely that a combination of jet entrainment and photoionization by the UV continuum from the central engine is needed in order to explain both the ionization and the kinematics of the gas in the Outer Filament.

We present neutral hydrogen (HI) and warm molecular hydrogen (H2) observations of the young (10^2 years) radio galaxy PKS B1718-649. We study the morphology and the kinematics of both gas components, focusing, in particular, on their properties in relation to the triggering of the radio activity. The regular kinematics of the large scale HI disk, seen in emission, suggests that an interaction event occurred too long ago to be responsible for the recent triggering of the radio activity. In absorption, we detect two absorption lines along the narrow line of sight of the compact (r<2 pc) radio source. The lines trace two clouds with opposite radial motions. These may represent a population of clouds in the very inner regions of the galaxy, which may be involved in triggering the radio activity. The warm molecular hydrogen (H2 1-0 S(1) ro-vibrational line) in the innermost kilo-parsec of the galaxy appears to be distributed in a circum-nuclear disk following the regular kinematics of the HI and of the stellar component. An exception to this behaviour arises only in the very centre, where a highly dispersed component is detected. These particular HI and H2 features suggest that a strong interplay between the radio source and the surrounding ISM is on-going. The physical properties of the cold gas in the proximity of the radio source may regulate the accretion recently triggered in this AGN.

Newly born and young radio sources are in a delicate phase of their life. Their jets are fighting their way through the surrounding gaseous medium, strongly experiencing this interaction while, at the same time, impacting and affecting the interstellar medium (ISM). Here we present the results from two studies of HI (in absorption) and molecular gas illustrating what can be learned from these phases of the gas. We first describe a statistical study with the WSRT. The study shows that the young radio sources not only have an higher detection rate of HI, but also systematically broader and more asymmetric HI profiles, most of them blueshifted. This supports the idea that we are looking at young radio jets making their way through the surrounding ISM, which also appears to be, on average, richer in gas than in evolved radio sources. Signatures of the impact of the jet are seen in the kinematics of the gas. However, even among the young sources, we identify a population that remains undetected in HI even after stacking their profiles. Orientation effects can only partly explain the result. These objects either are genuinely gas-poor or have different conditions of the medium, e.g. higher spin temperature. We further present the ALMA study of molecular gas in IC5063 to trace in detail the jet impacting the ISM. The kinematics of the cold, molecular gas co-spatial with the radio plasma shows this process in action. The ALMA data reveal a fast outflow of molecular gas extending along the entire radio jet (~1 kpc), with the highest outflow velocities at the location of the brighter hot-spot. We propose a scenario where the radio jet is expanding into a clumpy medium, interacting directly with the clouds and inflating a cocoon that drives a lateral outflow into the ISM.

Fast outflows of gas, driven by the interaction between the radio-jets and ISM of the host galaxy, are being observed in an increasing number of galaxies. One such example is the nearby radio galaxy 3C293. In this paper we present Integral Field Unit (IFU) observations taken with OASIS on the William Herschel Telescope (WHT), enabling us to map the spatial extent of the ionised gas outflows across the central regions of the galaxy. The jet-driven outflow in 3C293 is detected along the inner radio lobes with a mass outflow rate ranging from $\sim 0.05-0.17$ solar masses/yr (in ionised gas) and corresponding kinetic power of $\sim 0.5-3.5\times 10^{40}$ erg/s. Investigating the kinematics of the gas surrounding the radio jets (i.e. not directly associated with the outflow), we find line-widths broader than $300$ km/s up to 5 kpc in the radial direction from the nucleus (corresponding to 3.5 kpc in the direction perpendicular to the radio axis at maximum extent). Along the axis of the radio jet line-widths $>400$ km/s are detected out to 7 kpc from the nucleus and line-widths of $>500$ km/s at a distance of 12 kpc from the nucleus, indicating that the disturbed kinematics clearly extend well beyond the high surface brightness radio structures of the jets. This is suggestive of the cocoon structure seen in simulations of jet-ISM interaction and implies that the radio jets are capable of disturbing the gas throughout the central regions of the host galaxy in all directions.

We use high-resolution (0.5 arcsec) CO(2-1) observations performed with ALMA to trace the kinematics of the molecular gas in the Seyfert 2 galaxy IC5063. A fast outflow of molecular gas extends along the entire radio jet, with the highest outflow velocities about 0.5kpc from the nucleus, at the location of the brighter hot-spot in the W lobe. The data show that a massive, fast outflow with velocities up to 650 km/s of cold molecular gas is present, in addition to one detected earlier in warm H2, HI and ionised gas. Both the central AGN and the radio jet could energetically drive the outflow. However, the characteristics of the outflowing gas point to the radio jet being the main driver. This is important, because IC5063, although one of the most powerful Seyfert galaxies, is a relatively weak radio source (P = 3x10^23 W/Hz). All the observed characteristics can be described by a scenario of a radio plasma jet expanding into a clumpy medium, interacting directly with the clouds and inflating a cocoon that drives a lateral outflow into the interstellar medium. This model is consistent with results obtained by recent simulations such as those of Wagner et al.. A stronger, direct interaction between the jet and a gas cloud is present at the location of the brighter W lobe. Even assuming the most conservative values for the conversion factor CO-to-H2, the mass of the outflowing gas is between 1.9 and 4.8x10^7 Msun. These amounts are much larger than those of the outflow of warm gas (molecular and ionized) and somewhat larger than of the HI outflow. This suggests that most of the observed cold molecular outflow is due to fast cooling after being shocked. This gas is the end product of the cooling process. Our CO observations demonstrate that fast outflows of molecular gas can be driven by relativistic jets.

High-resolution (0.5 arcsec) CO(2-1) observations performed with the Atacama Large Millimetre/submillimetre Array have been used to trace the kinematics of the molecular gas in the Seyfert 2 galaxy{IC~5063}. Although one of the most radio-loud Seyfert galaxy, IC~5063 is a relatively weak radio source (P_1.4GHz = 3 x 10^23 W Hz^-1). The data reveal that the kinematics of the gas is very complex. A fast outflow of molecular gas extends along the entire radio jet (~ 1 kpc), with the highest outflow velocities about 0.5 kpc from the nucleus, at the location of the brighter hot-spot in the W lobe. All the observed characteristics can be described by a scenario of a radio plasma jet expanding into a clumpy medium, interacting directly with the clouds and inflating a cocoon that drives a lateral outflow into the interstellar medium. This suggests that most of the observed cold molecular outflow is due to fast cooling of the gas after the passage of a shock and that it is the end product of the cooling process.

We investigate signatures of a jet-interstellar medium (ISM) interaction using optical integral-field observations of the so-called outer filament near Centaurus A, expanding on previous results obtained on a more limited area. Using the Multi Unit Spectroscopic Explorer (MUSE) on the VLT during science verification, we observed a significant fraction of the brighter emitting gas across the outer filament. The ionized gas shows complex morphology with compact blobs, arc-like structures and diffuse emission. Based on the kinematics, we identified three main components. The more collimated component is oriented along the direction of the radio jet. The other two components exhibit diffuse morphology together with arc-like structures also oriented along the radio jet direction. Furthermore, the ionization level of the gas is found to decrease from the more collimated component to the more diffuse components. The morphology and velocities of the more collimated component confirm our earlier results that the outer filament and the nearby HI cloud are likely partially shaped by the lateral expansion of the jet. The arc-like structures embedded within the two remaining components are the clearest evidence of a smooth jet-ISM interaction along the jet direction. This suggests that, although poorly collimated, the radio jet is still active and has an impact on the surrounding gas. This result indicates that the effect on the ISM of even low-power radio jets should be considered when studying the influence Active Galactic Nuclei can have on their host galaxy.

The Square Kilometre Array (SKA) will transform our understanding of the role of the cold, atomic gas in galaxy evolution. The interstellar medium (ISM) is the repository of stellar ejecta and the birthsite of new stars and, hence, a key factor in the evolution of galaxies over cosmic time. Cold, diffuse, atomic clouds are a key component of the ISM, but so far this phase has been difficult to study, because its main tracer, the HI 21 cm line, does not constrain the basic physical information of the gas (e.g., temperature, density) well. The SKA opens up the opportunity to study this component of the ISM through a complementary tracer in the form of low-frequency (<350 MHz) carbon radio recombination lines (CRRL). These CRRLs provide a sensitive probe of the physical conditions in cold, diffuse clouds. The superb sensitivity, large field of view, frequency resolution and coverage of the SKA allows for efficient surveys of the sky, that will revolutionize the field of low-frequency recombination line studies. By observing these lines with the SKA we will be able determine the thermal balance, chemical enrichment, and ionization rate of the cold, atomic medium from degree-scales down to scales corresponding to individual clouds and filaments in our Galaxy, the Magellanic Clouds and beyond. Furthermore, being sensitive only to the cold, atomic gas, observations of low-frequency CRRLs with the SKA will aid in disentangling the warm and cold constituents of the HI 21 cm emission.

While the cooling of the hot intra-cluster medium (ICM) in the cores of galaxy clusters is mostly counteracted by heating from the central active galactic nucleus (AGN), the balance is not perfect. This can lead to residual cooling flows and low-level star formation, the physics of which is not well understood. Here we present a detailed study of the residual cooling flow in the center of the low mass galaxy cluster Sersic 159-03 (A S1101) using far-ultraviolet imaging from the Hubble Space Telescope and far-infrared (FIR) spectroscopy and photometry from the Herschel space observatory, along with a wealth of archival data. We detect extended emission at UV, FIR, and [CII], indicating a star formation rate of ~1-3 Msun/yr, depending on the indicator and assumptions made. The most recently formed stars appear spatially coincident with the lowest entropy ICM. We speculate that this low-entropy gas has been displaced by the central AGN ~7.5 kpc north of the cD galaxy. These data demonstrate that the displacement of the cooling core from the direct vicinity of the central AGN can temporarily break the feedback cycle and lead to cooling and star formation that is offset from the center of the galaxy. We find an abundance (~10^7 Msun) of cold (20K) dust in the center of the cluster and a second FIR peak ~30kpc to the north of the central galaxy. If confirmed to be associated with the cooling filaments, this would be the most extended complex of dust yet found in a cool core cluster.

Carbon radio recombination lines (RRLs) at low frequencies (<=500 MHz) trace the cold, diffuse phase of the interstellar medium, which is otherwise difficult to observe. We present the detection of carbon RRLs in absorption in M82 with LOFAR in the frequency range of 48-64 MHz. This is the first extragalactic detection of RRLs from a species other than hydrogen, and below 1 GHz. Since the carbon RRLs are not detected individually, we cross-correlated the observed spectrum with a template spectrum of carbon RRLs to determine a radial velocity of 219 +- 9 km/s . Using this radial velocity, we stack 22 carbon-{\alpha} transitions from quantum levels n = 468-508 to achieve an 8.5 sigma detection. The absorption line profile exhibits a narrow feature with peak optical depth of 0.003 and FWHM of 31 km/s. Closer inspection suggests that the narrow feature is superimposed on a broad, shallow component. The total line profile appears to be correlated with the 21 cm H I line profile reconstructed from H I absorption in the direction of supernova remnants in the nucleus. The narrow width and centroid velocity of the feature suggests that it is associated with the nuclear starburst region. It is therefore likely that the carbon RRLs are associated with cold atomic gas in the direction of the nucleus of M82.

We present LOFAR Low Band observations of the Bootes and 3C295 fields. Our images made at 34, 46, and 62 MHz reach noise levels of 12, 8, and 5 mJy beam$^{-1}$, making them the deepest images ever obtained in this frequency range. In total, we detect between 300 and 400 sources in each of these images, covering an area of 17 to 52 deg$^{2}$. From the observations we derive Euclidean-normalized differential source counts. The 62 MHz source counts agree with previous GMRT 153 MHz and VLA 74 MHz differential source counts, scaling with a spectral index of $-0.7$. We find that a spectral index scaling of $-0.5$ is required to match up the LOFAR 34 MHz source counts. This result is also in agreement with source counts from the 38 MHz 8C survey, indicating that the average spectral index of radio sources flattens towards lower frequencies. We also find evidence for spectral flattening using the individual flux measurements of sources between 34 and 1400 MHz and by calculating the spectral index averaged over the source population. To select ultra-steep spectrum ($\alpha < -1.1$) radio sources, that could be associated with massive high redshift radio galaxies, we compute spectral indices between 62 MHz, 153 MHz and 1.4 GHz for sources in the Bo\"otes field. We cross-correlate these radio sources with optical and infrared catalogues and fit the spectral energy distribution to obtain photometric redshifts. We find that most of these ultra-steep spectrum sources are located in the $ 0.7 \lesssim z \lesssim 2.5$ range.

Comparison of observations of radio recombination lines in the interstellar medium with theoretical models can be used to constrain electron temperature and density of the gas. An important component of the models is spontaneous transition rates between bound levels. Calculating these rates relies on accurate bound-bound oscillator strengths, which can be cast in terms of the Gaunt factor. The Gaunt factor contains terminating hypergeometric functions that cannot be calculated with sufficient accuracy for high quantum levels ($n \gtrsim 50$) by standard machine-precision methods. Methods to overcome the accuracy problem have been developed, which include asymptotic expansions and recursion relations. These methods, used in astrophysical models to calculate oscillator strengths, can introduce errors, sometimes up to as much as $\sim 8$ per cent. Detections of radio recombination lines with the new Low Frequency Array (LOFAR) has prompted an examination of theoretical models of the interstellar medium. We revisit the calculation of the Gaunt factor, employing modern arbitrary-precision computational methods to tabulate the Gaunt factor for transitions up to quantum level $n=2000$, sufficient to model low frequency Carbon radio recombination lines. The calculations provide a relative error of $\sim3\times 10^{-4}$ when compared to more detailed calculations including relativistic corrections. Our values for the Gaunt factor are provided for download in a tabular format to be used for a wide range of applications.

Massive outflows driven by active galactic nuclei (AGN) are widely recognised to play a key role in the evolution of galaxies, heating the ambient gas, expelling it from the nuclear regions, and thereby affecting the star formation histories of the galaxy bulges. It has been proposed that the powerful jets of relativistic particles launched by some AGN can both accelerate and heat the molecular gas, which often dominates the mass budgets of the outflows. However, clear evidence for this mechanism in the form of detailed associations between the molecular gas kinematics and features in the radio-emitting jets has been lacking. Here we show that the warm molecular hydrogen gas in the western radio lobe of the Seyfert galaxy IC5063 is moving at high velocities - up to 600 km/s - relative to the galaxy disk. This suggests that the molecules have been accelerated by fast shocks driven into the interstellar medium (ISM) by the expanding radio jets. These results demonstrate the general feasibility of accelerating molecular outflows in fast shocks driven by AGN.

We report ALMA Early Science observations of the Abell 1835 brightest cluster galaxy (BCG) in the CO (3-2) and CO (1-0) emission lines. We detect $5\times 10^{10}~\rm M_\odot$ of molecular gas within 10 kpc of the BCG. Its ensemble velocity profile width of $\sim 130 ~\rm km~s^{-1}$ FWHM is too narrow for the molecular cloud sto be supported in the galaxy by dynamic pressure. The gas may instead be supported in a rotating, turbulent disk oriented nearly face-on. Roughly $10^{10}~\rm M_\odot$ of molecular gas is projected $3-10 ~\rm kpc$ to the north-west and to the east of the nucleus with line of sight velocities lying between $-250 ~\rm km~s^{-1}$ to $+480 ~\rm km~s^{-1}$ with respect to the systemic velocity. The high velocity gas may be either inflowing or outflowing. However, the absence of high velocity gas toward the nucleus that would be expected in a steady inflow, and its bipolar distribution on either side of the nucleus, are more naturally explained as outflow. Star formation and radiation from the AGN are both incapable of driving an outflow of this magnitude. If so, the molecular outflow may be associated a hot outflow on larger scales reported by Kirkpatrick and colleagues. The molecular gas flow rate of approximately $200~\rm M_\odot ~yr^{-1}$ is comparable to the star formation rate of $100-180~\rm M_\odot ~yr^{-1}$ in the central disk. How radio bubbles would lift dense molecular gas in their updrafts, how much gas will be lost to the BCG, and how much will return to fuel future star formation and AGN activity are poorly understood. Our results imply that radio-mechanical (radio mode) feedback not only heats hot atmospheres surrounding elliptical galaxies and BCGs, it is able to sweep higher density molecular gas away from their centers.

We present the first detection of carbon radio recombination line absorption along the line of sight to Cygnus A. The observations were carried out with the LOw Frequency ARray in the 33 to 57 MHz range. These low frequency radio observations provide us with a new line of sight to study the diffuse, neutral gas in our Galaxy. To our knowledge this is the first time that foreground Milky Way recombination line absorption has been observed against a bright extragalactic background source. By stacking 48 carbon $\alpha$ lines in the observed frequency range we detect carbon absorption with a signal-to-noise ratio of about 5. The average carbon absorption has a peak optical depth of 2$\times$10$^{-4}$, a line width of 10 km s$^{-1}$ and a velocity of +4 km s$^{-1}$ with respect to the local standard of rest. The associated gas is found to have an electron temperature $T_{e}\sim$ 110 K and density $n_{e}\sim$ 0.06 cm$^{-3}$. These properties imply that the observed carbon $\alpha$ absorption likely arises in the cold neutral medium of the Orion arm of the Milky Way. Hydrogen and helium lines were not detected to a 3$\sigma$ peak optical depth limit of 1.5$\times$10$^{-4}$ for a 4 km s$^{-1}$ channel width. Radio recombination lines associated with Cygnus A itself were also searched for, but are not detected. We set a 3$\sigma$ upper limit of 1.5$\times$10$^{-4}$ for the peak optical depth of these lines for a 4 km s$^{-1}$ channel width.

We present a multi-wavelength study of the interstellar medium in eight nearby, X-ray and optically bright, giant elliptical galaxies. Using Herschel PACS, we map the cold gas in the lines of [CII], [OI], and [OIb]. Additionally, we present Ha+[NII] imaging of warm ionized gas with the SOAR telescope, and a study of the hot X-ray emitting plasma with Chandra. All systems with extended Ha emission in our sample (6/8 galaxies) display significant [CII] line emission indicating the presence of cold gas. This emission is co-spatial with the Ha+[NII] emitting nebulae and the lowest entropy X-ray emitting plasma. The entropy profiles of the hot galactic atmospheres show a clear dichotomy, with the systems displaying extended emission line nebulae having lower entropies beyond r~1 kpc than the cold-gas-poor systems. We show that while the hot atmospheres of the cold-gas-poor galaxies are thermally stable outside of their innermost cores, the atmospheres of the cold-gas-rich systems are prone to cooling instabilities. This result indicates that the cold gas is produced chiefly by thermally unstable cooling from the hot phase. We show that cooling instabilities may develop more easily in rotating systems and discuss an alternative condition for thermal instability for this case. The hot atmospheres of cold-gas-rich galaxies display disturbed morphologies indicating that the accretion of clumpy multiphase gas in these systems may result in variable power output of the AGN jets, potentially triggering sporadic, larger outbursts. In the two cold-gas-poor, X-ray morphologically relaxed galaxies of our sample, NGC 1399 and NGC 4472, powerful AGN outbursts may have destroyed or removed most of the cold gas from the cores, allowing the jets to propagate and deposit most of their energy further out, increasing the entropy of the hot galactic atmospheres and leaving their cores relatively undisturbed.

We present multi-frequency observations of the radio galaxy Hydra-A (3C218) located in the core of a massive, X-ray luminous galaxy cluster. IFU spectroscopy is used to trace the kinematics of the ionised and warm molecular hydrogen which are consistent with a ~ 5 kpc rotating disc. Broad, double-peaked lines of CO(2-1), [CII]157 $\mu$m and [OI]63 $\mu$m are detected. We estimate the mass of the cold gas within the disc to be M$_{gas}$ = 2.3 $\pm$ 0.3 x 10$^9$ M$_{\odot}$. These observations demonstrate that the complex line profiles found in the cold atomic and molecular gas are related to the rotating disc or ring of gas. Finally, an HST image of the galaxy shows that this gas disc contains a substantial mass of dust. The large gas mass, SFR and kinematics are consistent with the levels of gas cooling from the ICM. We conclude that the cold gas originates from the continual quiescent accumulation of cooled ICM gas. The rotation is in a plane perpendicular to the projected orientation of the radio jets and ICM cavities hinting at a possible connection between the kpc-scale cooling gas and the accretion of material onto the black hole. We discuss the implications of these observations for models of cold accretion, AGN feedback and cooling flows.

We report ALMA Early Science CO(1-0) and CO(3-2) observations of the brightest cluster galaxy (BCG) in Abell 1664. The BCG contains 1.1x10^{10} solar masses of molecular gas divided roughly equally between two distinct velocity systems: one from -250 to +250 km/s centred on the BCG's systemic velocity and a high velocity system blueshifted by 570 km/s with respect to the systemic velocity. The BCG's systemic component shows a smooth velocity gradient across the BCG center with velocity proportional to radius suggestive of solid body rotation about the nucleus. However, the mass and velocity structure are highly asymmetric and there is little star formation coincident with a putative disk. It may be an inflow of gas that will settle into a disk over several 10^8 yr. The high velocity system consists of two gas clumps, each ~2 kpc across, located to the north and southeast of the nucleus. Each has a line of sight velocity spread of 250-300 km/s. The velocity of the gas in the high velocity system tends to increase towards the BCG center and could signify a massive high velocity flow onto the nucleus. However, the velocity gradient is not smooth and these structures are also coincident with low optical-UV surface brightness regions, which could indicate dust extinction associated with each clump. If so, the high velocity gas would be projected in front of the BCG and moving toward us along the line of sight in a massive outflow most likely driven by the AGN. A merger origin is unlikely but cannot be ruled out.

We report ALMA Early Science observations of the Abell 1835 brightest cluster galaxy (BCG) in the CO (3-2) and CO (1-0) emission lines. We detect 5E10 solar masses of molecular gas within 10 kpc of the BCG. Its velocity width of ~130 km/s FWHM is too narrow to be supported by dynamical pressure. The gas may instead be supported in a rotating, turbulent disk oriented nearly face-on. The disk is forming stars at a rate of 100-180 solar masses per year. Roughly 1E10 solar masses of molecular gas is projected 3-10 kpc to the north-west and to the east of the nucleus with line of sight velocities lying between -250 km/s to +480 km/s with respect to the systemic velocity. Although inflow cannot be ruled out, the rising velocity gradient with radius is consistent with a broad, bipolar outflow driven by radio jets or buoyantly rising X-ray cavities. The molecular outflow may be associated with an outflow of hot gas in Abell 1835 seen on larger scales. Molecular gas is flowing out of the BCG at a rate of approximately 200 solar masses per year, which is comparable to its star formation rate. How radio bubbles lift dense molecular gas in their updrafts, how much gas will be lost to the BCG, and how much will return to fuel future star formation and AGN activity are poorly understood. Our results imply that radio-mechanical (radio mode) feedback not only heats hot atmospheres surrounding elliptical galaxies and BCGs, it is able to sweep higher density molecular gas away from their centers.

We present a multi-wavelength study of the emission-line nebulae located southeast of the nucleus of M87, the central dominant galaxy of the Virgo Cluster. We report the detection of far-infrared (FIR) [CII] line emission from the nebulae using observations made with Herschel PACS. The infrared line emission is extended and cospatial with optical H{\alpha}+[NII], far-ultraviolet CIV lines, and soft X-ray emission. The filamentary nebulae evidently contain multi-phase material spanning a temperature range of at least 5 orders of magnitude, from ~100 K to ~10^7 K. This material has most likely been uplifted by the AGN from the center of M87. The thermal pressure of the 10^4 K phase appears to be significantly lower than that of the surrounding hot intra-cluster medium (ICM) indicating the presence of additional turbulent and magnetic pressure in the filaments. If the turbulence in the filaments is subsonic then the magnetic field strength required to balance the pressure of the surrounding ICM is B~30-70 {\mu}G. The spectral properties of the soft X-ray emission from the filaments indicate that it is due to thermal plasma with kT~0.5-1 keV, which is cooling by mixing with the cold gas and/or radiatively. Charge exchange can be ruled out as a significant source of soft X-rays. Both cooling and mixing scenarios predict gas with a range of temperatures. This is at first glance inconsistent with the apparent lack of X-ray emitting gas with kT<0.5 keV. However, we show that the missing very soft X-ray emission could be absorbed by the cold gas in the filaments with an integrated hydrogen column density of ~1.6x10^21 cm^-2, providing a natural explanation for the apparent temperature floor to the X-ray emission at kT~0.5 keV. The FIR through ultra-violet line emission is most likely primarily powered by the ICM particles penetrating the cold gas following a shearing induced mixing process.

Cassiopeia A was observed using the Low-Band Antennas of the LOw Frequency ARray (LOFAR) with high spectral resolution. This allowed a search for radio recombination lines (RRLs) along the line-of-sight to this source. Five carbon-alpha RRLs were detected in absorption between 40 and 50 MHz with a signal-to-noise ratio of > 5 from two independent LOFAR datasets. The derived line velocities (v_LSR ~ -50 km/s) and integrated optical depths (~ 13 s^-1) of the RRLs in our spectra, extracted over the whole supernova remnant, are consistent within each LOFAR dataset and with those previously reported. For the first time, we are able to extract spectra against the brightest hotspot of the remnant at frequencies below 330 MHz. These spectra show significantly higher (15-80 %) integrated optical depths, indicating that there is small-scale angular structure on the order of ~1 pc in the absorbing gas distribution over the face of the remnant. We also place an upper limit of 3 x 10^-4 on the peak optical depths of hydrogen and helium RRLs. These results demonstrate that LOFAR has the desired spectral stability and sensitivity to study faint recombination lines in the decameter band.

Some merging galaxy clusters host diffuse extended radio emission, so-called radio halos and relics. Here we present observations between 147 MHz and 4.9 GHz of a new radio-selected galaxy cluster 1RXS J0603.3+4214 (z=0.225). The cluster is also detected as an extended X-ray source in the RASS. It hosts a large bright 1.9 Mpc radio relic, an elongated ~2 Mpc radio halo, and two smaller radio relics. The large radio relic has a peculiar linear morphology. For this relic we observe a clear spectral index gradient, in the direction towards the cluster center. We performed Rotation Measure (RM) Synthesis between 1.2 and 1.7 GHz. The results suggest that for the west part of the large relic some of the Faraday rotation is caused by ICM and is not only due to galactic foregrounds. We also carried out a detailed spectral analysis of this radio relic and created radio color-color diagrams. We find (i) an injection spectral index of -0.6 to -0.7, (ii) steepening spectral index and increasing spectral curvature in the post-shock region, and (iii) an overall power-law spectrum between 74 MHz and 4.9 GHz with \alpha=-1.10 \pm 0.02. Mixing of emission in the beam from regions with different spectral ages is probably the dominant factor that determines the shape of the radio spectra. Changes in the magnetic field, total electron content, or adiabatic gains/losses do not play a major role. A model in which particles are (re)accelerated in a first order Fermi process at the front of the relic provides the best match to the observed spectra. We speculate that in the post-shock region particles are re-accelerated by merger induced turbulence to form the radio halo as the relic and halo are connected. The 1RXS J0603.3+4214 merger is probably more complex than the "simple'" binary merger events that are thought to give rise to symmetric double radio relics.

We present Herschel observations of the core of the Perseus cluster of galaxies. The brightest cluster galaxy, NGC 1275, is surrounded by a network of filaments previously imaged extensively in H{\alpha} and CO. In this work, we report detections of FIR lines with Herschel. All but one of the lines are spatially extended, with the [CII] line emission extending up to 25 kpc from the core. There is spatial and kinematical correlation among [CII], H{\alpha} and CO, which gives us confidence to model the different components of the gas with a common heating model. With the help of FIR continuum Herschel measurements, together with a suite of coeval radio, submm and infrared data, we performed a SED fitting of NGC 1275 using a model that contains contributions from dust emission as well as synchrotron AGN emission. The data indicate a low dust emissivity index, beta ~ 1, a total dust mass close to 10^7 solar mass, a cold dust component with temperature 38 \pm 2 K and a warm dust component with temperature of 116 \pm 9 K. The FIR-derived star formation rate (SFR) is 24 \pm 1 solar mass per yr, in close agreement with the FUV-derived SFR. We investigated in detail the source of the Herschel FIR and H{\alpha} emissions emerging from a core region 4 kpc in radius. Based on simulations conducted using the radiative transfer code, CLOUDY, a heating model comprising old and young stellar populations is sufficient to explain these observations. We have also detected [CII] in three well-studied regions of the filaments. We find a [OI]/[CII] ratio about 1 dex smaller than predicted by the otherwise functional Ferland (2009) model. The line ratio suggests that the lines are optically thick, as is typical of galactic PDRs, and implies that there is a large reservoir of cold atomic gas. [abridged]

IRAS 09104+4109 is a rare example of a dust enshrouded type 2 QSO in the centre of a cool-core galaxy cluster. Previous observations of this z=0.44 system showed that as well as powering the hyper-luminous infrared emission of the cluster-central galaxy, the QSO is associated with a double-lobed radio source. However, the steep radio spectral index and misalignment between the jets and ionised optical emission suggested that the orientation of the QSO had recently changed. We use a combination of new, multi-band Giant Metrewave Radio Telescope observations and archival radio data to confirm that the jets are no longer powered by the QSO, and estimate their age to be 120-160 Myr. This is in agreement with the ~70-200 Myr age previously estimated for star-formation in the galaxy. Previously unpublished Very Long Baseline Array data reveal a 200 pc scale double radio source in the galaxy core which is more closely aligned with the current QSO axis and may represent a more recent period of jet activity. These results suggest that the realignment of the QSO, the cessation of jet activity, and the onset of rapid star-formation may have been caused by a gas-rich galaxy merger. A Chandra X-ray observation confirms the presence of cavities associated with the radio jets, and we estimate the energy required to inflate them to be ~7.7x10^60 erg. The mechanical power of the jets is sufficient to balance radiative cooling in the cluster, provided they are efficiently coupled to the intra-cluster medium (ICM). We find no evidence of direct radiative heating and conclude that the QSO either lacks the radiative luminosity to heat the ICM, or that it requires longer than 100-200 Myr to significantly impact its environment. [Abridged]

We present new Chandra X-ray observations of the brightest cluster galaxy (BCG) in the cool core cluster Abell 2597. The data reveal an extensive kpc-scale X-ray cavity network as well as a 15 kpc filament of soft-excess gas exhibiting strong spatial correlation with archival VLA radio data. In addition to several possible scenarios, multiwavelength evidence may suggest that the filament is associated with multiphase (10^3 - 10^7 K) gas that has been entrained and dredged-up by the propagating radio source. Stemming from a full spectral analysis, we also present profiles and 2D spectral maps of modeled X-ray temperature, entropy, pressure, and metal abundance. The maps reveal an arc of hot gas which in projection borders the inner edge of a large X-ray cavity. Although limited by strong caveats, we suggest that the hot arc may be (a) due to a compressed rim of cold gas pushed outward by the radio bubble or (b) morphologically and energetically consistent with cavity-driven active galactic nucleus (AGN) heating models invoked to quench cooling flows, in which the enthalpy of a buoyant X-ray cavity is locally thermalized as ambient gas rushes to refill its wake. If confirmed, this would be the first observational evidence for this model.

New Chandra X-ray and Herschel FIR observations enable a multiwavelength study of active galactic nucleus (AGN) heating and intracluster medium (ICM) cooling in the brightest cluster galaxy of Abell 2597. The new Chandra observations reveal the central < 30 kiloparsec X-ray cavity network to be more extensive than previously thought, and associated with enough enthalpy to theoretically inhibit the inferred classical cooling flow. Nevertheless, we present new evidence, consistent with previous results, that a moderately strong residual cooling flow is persisting at 4%-8% of the classically predicted rates in a spatially structured manner amid the feedback-driven excavation of the X-ray cavity network. New Herschel observations are used to estimate warm and cold dust masses, a lower-limit gas-to-dust ratio, and a star formation rate consistent with previous measurements. The cooling time profile of the ambient X-ray atmosphere is used to map the locations of the observational star formation entropy threshold as well as the theoretical thermal instability threshold. Both lie just outside the < 30 kpc central region permeated by X-ray cavities, and star formation as well as ionized and molecular gas lie interior to both. The young stars are distributed in an elongated region that is aligned with the radio lobes, and their estimated ages are both younger and older than the X-ray cavity network, suggesting both jet-triggered as well as persistent star formation over the current AGN feedback episode. Bright X-ray knots that are coincident with extended Ly-alpha and FUV continuum filaments motivate a discussion of structured cooling from the ambient hot atmosphere along a projected axis that is perpendicular to X-ray cavity and radio axis. We conclude that the cooling ICM is the dominant contributor of the cold gas reservoir fueling star formation and AGN activity in the Abell 2597 BCG.

We use the Hubble Space Telescope ACS/SBC and Very Large Telescope FORS cameras to observe the Brightest Cluster Galaxies in Abell 2597 and Abell 2204 in the far-ultraviolet (FUV) F150LP and optical U, B, V, R, I Bessel filters. The FUV and U band emission is enhanced in bright, filamentary structures surrounding the BCG nuclei. These filaments can be traced out to 20 kpc from the nuclei in the FUV. Excess FUV and U band light is determined by removing emission due to the underlying old stellar population and mapped with 1 arcsec spatial resolution over the central 20 kpc regions of both galaxies. We find the FUV and U excess emission to be spatially coincident and a stellar interpretation requires the existence of a significant amount of 10000-50000 K stars. Correcting for nebular continuum emission and dust intrinsic to the BCG further increases the FUV to U band emission ratio and implies that stars alone may not suffice to explain the observations. However, lack of detailed information on the gas and dust distribution and extinction law in these systems prevents us from ruling out a purely stellar origin. Non-stellar processes, such as the central AGN, Scattering, Synchrotron and Bremsstrahlung emission are investigated and found to not be able to explain the FUV and U band measurements in A2597. Contributions from non-thermal processes not treated here should be investigated. Comparing the FUV emission to the optical H-alpha line emitting nebula shows good agreement on kpc-scales in both A2597 and A2204. In concordance with an earlier investigation by O'Dea et al. (2004) we find that O-stars can account for the ionising photons necessary to explain the observed H-alpha line emission.

In order to reconstruct the global SEDs of the Magellanic Clouds over eight decades in spectral range, we combined literature flux densities representing the entire LMC and SMC respectively, and complemented these with maps extracted from the WMAP and COBE databases covering the missing the 23--90 GHz (13--3.2 mm) and the poorly sampled 1.25--250 THz (240--1.25 micron). We have discovered a pronounced excess of emission from both Magellanic Clouds, but especially the SMC, at millimeter and sub-millimeter wavelengths. We also determined accurate thermal radio fluxes and very low global extinctions for both LMC and SMC. Possible explanations are briefly considered but as long as the nature of the excess emission is unknown, the total dust masses and gas-to-dust ratios of the Magellanic Clouds cannot reliably be determined.

The question of how much gas cools in the cores of clusters of galaxies has been the focus of many, multiwavelength studies in the past 30 years. In this letter we present the first detections of the strongest atomic cooling lines, [C II], [O I] and [N I] in two strong cooling flow clusters, A1068 and A2597, using Herschel PACS. These spectra indicate that the substantial mass of cold molecular gas (>10^9 Mo) known to be present in these systems is being irradiated by intense UV radiation, most probably from young stars. The line widths of these FIR lines indicate that they share dynamics similar but not identical to other ionised and molecular gas traced by optical, near-infrared and CO lines. The relative brightness of the FIR lines compared to CO and FIR luminosity is consistent with other star-forming galaxies indicating that the properties of the molecular gas clouds in cluster cores and the stars they form are not unusual. These results provide additional evidence for a reservoir of cold gas that is fed by the cooling of gas in the cores of the most compact clusters and provide important diagnostics of the temperature and density of the dense clouds this gas resides in.

The dust destruction timescales in the cores of clusters of galaxies are relatively short given their high central gas densities. However, substantial mid-infrared and sub-mm emission has been detected in many brightest cluster galaxies. In this letter we present Herschel PACS and SPIRE photometry of the brightest cluster galaxy in three strong cooling flow clusters, A1068, A2597 and Zw3146. This photometry indicates that a substantial mass of cold dust is present (>3 x 10^7 Mo) at temperatures significantly lower (20-28K) than previously thought based on limited MIR and/or sub-mm results. The mass and temperature of the dust appear to match those of the cold gas traced by CO with a gas-to-dust ratio of 80-120.

We have used the SINFONI integral field spectrograph to map the near-infrared K-band emission lines of molecular and ionised hydrogen in the central regions of two cool core galaxy clusters, Abell 2597 and Sersic 159-03. Gas is detected out to 20 kpc from the nuclei of the brightest cluster galaxies and found to be distributed in clumps and filaments around it. The ionised and molecular gas phases trace each other closely in extent and dynamical state. Both gas phases show signs of interaction with the active nucleus. Within the nuclear regions the kinetic luminosity of this gas is found to be somewhat smaller than the current radio luminosity. Outside the nuclear region the gas has a low velocity dispersion and shows smooth velocity gradients. There is no strong correlation between the intensity of the molecular and ionised gas emission and either the radio or X-ray emission. The molecular gas in Abell 2597 and Sersic 159-03 is well described by a gas in local thermal equilibrium (LTE) with a single excitation temperature T_exc ~ 2300 K. The emission line ratios do not vary strongly as function of position, with the exception of the nuclear regions where the ionised to molecular gas ratio is found decrease. These constant line ratios imply a single source of heating and excitation for both gas phases.

Context: Hierarchical galaxy formation models indicate that during their lifetime galaxy clusters undergo several mergers. Here we report on the discovery of three diffuse radio sources in the periphery of Abell 2256, using the GMRT. Aims: The aim of the observations was to search for diffuse ultra-steep spectrum radio sources within the galaxy cluster Abell 2256. Methods: We have carried out GMRT 325 MHz radio continuum observations of Abell 2256. V, R and I band images of the cluster were taken with the 4.2m WHT. Results: We have discovered three diffuse elongated radio sources located about 1 Mpc from the cluster center. Two are located to the west of the cluster center, and one to the southeast. The sources have a measured physical extent of 170, 140 and 240 kpc, respectively. The two western sources are also visible in deep low-resolution 115-165 MHz WSRT images, although they are blended into a single source. For the combined emission of the blended source we find an extreme spectral index of -2.05\pm 0.14 between 140 and 351 MHz. The extremely steep spectral index suggests these two sources are most likely the result of adiabatic compression of fossil radio plasma due to merger shocks. Conclusions: The discovery of the steep spectrum sources implies the existence of a population of faint diffuse radio sources in (merging) clusters with such steep spectra that they have gone unnoticed in higher frequency (\gtrsim 1 GHz) observations. An exciting possibility therefore is that such sources will determine the general appearance of clusters in low-frequency high resolution radio maps as will be produced by for example LOFAR or LWA.