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.

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.

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.

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.

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.

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

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.