We present LOFAR observations at 120-168 MHz of 42 systems with possible X-ray cavities in their hot atmosphere, of which 17 are groups or ellipticals, 19 are nearby clusters (z<0.3), and six are higher-redshift clusters (z>0.3). The X-ray cavities, formed by the radio lobes of the central active galactic nucleus (AGN), are evidence of radio-mode AGN feedback. In the groups and ellipticals sample, more than half of the systems have X-ray cavities for which no associated lobe emission was detected. Conversely, we report the discovery of large radio lobes in NGC 6338, extending far beyond the cavities reported previously in the literature. In the case of the nearby clusters, our observations show that there is little low-frequency radio emission that extends beyond the cavities (e.g., MS 0735.6+7421 and A2052). For the first time, we report secure cavity-radio associations in 2A 0335+096, ZwCl 2701, and ZwCl 8276 that strengthens their interpretation as AGN-created cavities. However, in some known cavity systems (e.g., A1795 and ZwCl 3146) we report the lack of detectable low-frequency radio emission associated with the cavities. Our sample of higher-redshift systems is small, and unfortunately the present LOFAR observations are not able to resolve the lobes in many of them. Nevertheless, our sample represents one of the best available for investigating the connection between radio and jet power in radio-mode AGN feedback.

We have examined the giant radio galaxy 3C~236 using LOFAR at 143 MHz down to an angular resolution of 7", in combination with observations at higher frequencies. We have used the low frequency data to derive spectral index maps with the highest resolution yet at these low frequencies. We confirm a previous detection of an inner hotspot in the north-west lobe and for the first time observe that the south-east lobe hotspot is in fact a triple hotspot, which may point to an intermittent source activity. Also, the spectral index map of 3C 236 shows that the spectral steepening at the inner region of the northern lobe is prominent at low frequencies. The outer regions of both lobes show spectral flattening, in contrast with previous high frequency studies. We derive spectral age estimates for the lobes, as well as particle densities of the IGM at various locations. We propose that the morphological differences between the lobes are driven by variations in the ambient medium density as well as the source activity history.

The largest galaxy clusters are observed still to be forming through major cluster-cluster mergers, often showing observational signatures such as radio relics and giant radio haloes. Using LOFAR Two-meter Sky Survey data, we present new detections of both a radio halo (with a spectral index of $\alpha_{143}^{1400}=1.48^{+0.06}_{-0.23}$) and a likely radio relic in Abell 959, a massive cluster at a redshift of z=0.288. Using a sample of clusters with giant radio haloes from the literature (80 in total), we show that the radio halo in A959 lies reasonably well on the scaling relations between the thermal and non-thermal power of the system. Additionally, we find evidence that steep-spectrum haloes tend to reside in clusters with high X-ray luminosities relative to those expected from cluster LM scaling relations, indicating that such systems may preferentially lie at an earlier stage of the merger, consistent with the theory that some steep-spectrum haloes result from low-turbulence mergers. Lastly, we find that halo systems containing radio relics tend to lie at lower X-ray luminosities, relative to those expected from cluster LM scaling relations, for a given halo radio power than those without relics, suggesting that the presence of relics indicates a later stage of the merger, in line with simulations.

Centrally located diffuse radio emission has been observed in both merging and non-merging galaxy clusters. Depending on their morphology and size, we distinguish between giant radio haloes, which occur predominantly in merging clusters, and mini haloes, which are found in non-merging, cool-core clusters. Low-frequency sensitive observations are required to assess whether the emission discovered in these few cases is common in galaxy clusters or not. With this aim, we carried out a campaign of observations with the LOw Frequency ARay (LOFAR) in the frequency range 120 - 168 MHz of nine massive clusters selected from the \textit{Planck} SZ catalogue, which had no sign of major mergers. In this paper, we discuss the results of the observations that have led to the largest cluster sample studied within the LOFAR Two-metre Sky Survey, and we present Chandra X-ray data used to investigate the dynamical state of the clusters, verifying that the clusters are currently not undergoing major mergers, and to search for traces of minor or off-axis mergers. We discover large-scale steep-spectrum emission around mini haloes in the cool-core clusters PSZ1G139 and RXJ1720, which is not observed around the mini halo in the non-cool-core cluster A1413. We also discover a new 570 kpc-halo in the non-cool-core cluster RXCJ0142. We derived new upper limits to the radio power for clusters in which no diffuse radio emission was found, and we discuss the implication of our results to constrain the cosmic-ray energy budget in the ICM. We conclude that radio emission in non-merging massive clusters is not common at the sensitivity level reached by our observations and that no clear connection with the cluster dynamical state is observed. Our results might indicate that the sloshing of a dense cool core could trigger particle acceleration on larger scales and generate steep-spectrum radio emission.

A number of radio observations have revealed the presence of large synchrotron-emitting sources associated with the intra-cluster medium. There is strong observational evidence that the emitting particles have been (re-)accelerated by shocks and turbulence generated during merger events. The particles that are accelerated are thought to have higher initial energies than those in the thermal pool but the origin of such mildly relativistic particles remains uncertain and needs to be further investigated. The galaxy cluster Abell 1914 is a massive galaxy cluster in which X-ray observations show clear evidence of merging activity. We carried out radio observations of this cluster with the LOw Frequency ARay (LOFAR) at 150 MHz and the Giant Metrewave Radio Telescope (GMRT) at 610 MHz. We also analysed Very Large Array (VLA) 1.4 GHz data, archival GMRT 325 MHz data, CFHT weak lensing data and Chandra observations. Our analysis shows that the ultra-steep spectrum source (4C38.39; $\alpha \lesssim -2$), previously thought to be part of a radio halo, is a distinct source with properties that are consistent with revived fossil plasma sources. Finally, we detect some diffuse emission to the west of the source 4C38.39 that could belong to a radio halo.

A number of merging galaxy clusters show the presence of large-scale radio emission associated with the intra-cluster medium (ICM). These synchrotron sources are generally classified as radio haloes and radio relics. Whilst it is commonly accepted that mergers play a crucial role in the formation of radio haloes and relics, not all the merging clusters show the presence of giant diffuse radio sources and this provides important information concerning current models. The Abell 781 complex is a spectacular system composed of an apparent chain of clusters on the sky. Its main component is undergoing a merger and hosts peripheral emission that is classified as a candidate radio relic and a disputed radio halo. We used new LOw Frequency ARay (LOFAR) observations at 143 MHz and archival Giant Metrewave Radio Telescope (GMRT) observations at 325 and 610 MHz to study radio emission from non-thermal components in the ICM of Abell 781. Complementary information came from XMM-Newton data, which allowed us to investigate the connection with the thermal emission and its complex morphology. The origin of the peripheral emission is still uncertain. We speculate that it is related to the interaction between a head tail radio galaxy and shock. However, the current data allow us only to set an upper limit of $\mathcal{M} < 1.4$ on the Mach number of this putative shock. Instead, we successfully characterise the surface brightness and temperature jumps of a shock and two cold fronts in the main cluster component of Abell 781. Their positions suggest that the merger is involving three substructures. We do not find any evidence for a radio halo either at the centre of this system or in the other clusters of the chain. We place an upper limit to the diffuse radio emission in the main cluster of Abell 781 that is a factor of 2 below the current radio power-mass relation for giant radio haloes.

Massive, merging galaxy clusters often host giant, diffuse radio sources that arise from shocks and turbulence; hence, radio observations can be useful for determining the merger state of a cluster. In preparation for a larger study, we selected three clusters -- Abell 1319, Abell 1314, and RXC J1501.3+4220 (Z7215) -- making use of the new LOFAR Two-Metre Sky Survey (LoTSS) at 120-168 MHz, and together with archival data, show that these clusters appear to be in pre-merging, merging, and post-merging states, respectively. We argue that Abell 1319 is likely in its pre-merging phase, where three separate cluster components are about to merge. There are no radio halos nor radio relics detected in this system. Abell 1314 is a highly-disturbed, low-mass cluster which is likely in the process of merging. This low-mass system does not show a radio halo, however, we argue that the merger activates mechanisms that cause electron re-acceleration in the large 800 kpc radio tail associated with IC~711. In the cluster Z7215 we discover diffuse radio emission at the cluster center, and we classify this emission as a radio halo, although it is dimmer and smaller than expected by the radio halo power versus cluster mass correlation. We suggest that the disturbed cluster Z7215 is in its post-merging phase. Systematic studies of this kind over a larger sample of clusters observed with LoTSS will help constrain the time scales involved in turbulent re-acceleration and the subsequent energy losses of the underlying electrons.

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.

We examine the possible acceleration mechanisms of the relativistic particles responsible for the extended radio emission in Abell 520. We used new LOFAR 145 MHz, archival GMRT 323 MHz and VLA 1.5 GHz data to study the morphological and spectral properties of extended cluster emission. The observational properties are discussed in the framework of particle acceleration models associated with cluster merger turbulence and shocks. In Abell 520, we confirm the presence of extended synchrotron radio emission that has been classified as a radio halo. The comparison between the radio and X-ray brightness suggests that the halo might originate in a cocoon rather than from the central X-ray bright regions of the cluster. The halo spectrum is roughly uniform on the scale of 66 kpc. There is a hint of spectral steepening from the SW edge towards the cluster centre. Assuming DSA, the radio data are suggestive of a shock of $\mathcal{M}_{SW}=2.6_{-0.2}^{+0.3}$ that is consistent with the X-ray derived estimates. This is in line with the scenario in which relativistic electrons in the SW radio edge gain their energies at the shock front via acceleration of either thermal or fossil electrons. We do not detect extended radio emission ahead of the SW shock that is predicted if the emission is the result of adiabatic compression. An X-ray surface brightness discontinuity is detected towards the NE region that may be a counter shock of $\mathcal{M}_{NE}^{X}=1.52\pm0.05$. This is lower than the value predicted from the radio emission ($\mathcal{M}_{NE}=2.1\pm0.2$). Our observations indicate that the SW radio emission in Abell 520 is likely effected by the prominent X-ray detected shock in which radio emitting particles are (re-)accelerated through the Fermi-I mechanism. The NE X-ray discontinuity that is approximately collocated with an edge in the radio emission hints at the presence of a counter shock.

Radio continuum (RC) emission in galaxies allows us to measure star formation rates (SFRs) unaffected by extinction due to dust, of which the low-frequency part is uncontaminated from thermal (free-free) emission. We calibrate the conversion from the spatially resolved 140 MHz RC emission to the SFR surface density ($\Sigma_{\rm SFR}$) at 1 kpc scale. We used recent observations of three galaxies (NGC 3184, 4736, and 5055) from the LOFAR Two-metre Sky Survey (LoTSS), and archival LOw-Frequency ARray (LOFAR) data of NGC 5194. Maps were created with the facet calibration technique and converted to radio $\Sigma_{\rm SFR}$ maps using the Condon relation. We compared these maps with hybrid $\Sigma_{\rm SFR}$ maps from a combination of GALEX far-ultraviolet and Spitzer 24 $\mu\rm m$ data using plots tracing the relation at $1.2\times 1.2$-kpc$^2$ resolution. The RC emission is smoothed with respect to the hybrid $\Sigma_{\rm SFR}$ owing to the transport of cosmic-ray electrons (CREs). This results in a sublinear relation $(\Sigma_{\rm SFR})_{\rm RC} \propto [(\Sigma_{\rm SFR})_{\rm hyb}]^{a}$, where $a=0.59\pm 0.13$ (140 MHz) and $a=0.75\pm 0.10$ (1365 MHz). Both relations have a scatter of $\sigma = 0.3~\rm dex$. If we restrict ourselves to areas of young CREs ($\alpha > -0.65$; $I_\nu \propto \nu^\alpha$), the relation becomes almost linear at both frequencies with $a\approx 0.9$ and a reduced scatter of $\sigma = 0.2~\rm dex$. We then simulate the effect of CRE transport by convolving the hybrid $\Sigma_{\rm SFR}$ maps with a Gaussian kernel until the RC-SFR relation is linearised; CRE transport lengths are $l=1$-5 kpc. Solving the CRE diffusion equation, we find diffusion coefficients of $D=(0.13$-$1.5) \times 10^{28} \rm cm^2\,s^{-1}$ at 1 GeV. A RC-SFR relation at $1.4$ GHz can be exploited to measure SFRs at redshift $z \approx 10$ using $140$ MHz observations.

Context: New generation low-frequency telescopes are exploring a new parameter space in terms of depth and resolution. The data taken with these interferometers, for example with the LOw Frequency ARray (LOFAR), are often calibrated in a low signal-to-noise ratio regime and the removal of critical systematic effects is challenging. The process requires an understanding of their origin and properties. Aim: In this paper we describe the major systematic effects inherent to next generation low-frequency telescopes, such as LOFAR. With this knowledge, we introduce a data processing pipeline that is able to isolate and correct these systematic effects. The pipeline will be used to calibrate calibrator observations as the first step of a full data reduction process. Methods: We processed two LOFAR observations of the calibrator 3C196: the first using the Low Band Antenna (LBA) system at 42-66 MHz and the second using the High Band Antenna (HBA) system at 115-189 MHz. Results: We were able to isolate and correct for the effects of clock drift, polarisation misalignment, ionospheric delay, Faraday rotation, ionospheric scintillation, beam shape, and bandpass. The designed calibration strategy produced the deepest image to date at 54 MHz. The image has been used to confirm that the spectral energy distribution of the average radio source population tends to flatten at low frequencies. Conclusions: We prove that LOFAR systematic effects can be described by a relatively small number of parameters. Furthermore, the identification of these parameters is fundamental to reducing the degrees of freedom when the calibration is carried out on fields that are not dominated by a strong calibrator.

Diffuse synchrotron radio emission from cosmic-ray electrons is observed at the center of a number of galaxy clusters. These sources can be classified either as giant radio halos, which occur in merging clusters, or as mini halos, which are found only in cool-core clusters. In this paper, we present the first discovery of a cool-core cluster with an associated mini halo that also shows ultra-steep-spectrum emission extending well beyond the core that resembles radio halo emission. The large-scale component is discovered thanks to LOFAR observations at 144 MHz. We also analyse GMRT observations at 610 MHz to characterise the spectrum of the radio emission. An X-ray analysis reveals that the cluster is slightly disturbed, and we suggest that the steep-spectrum radio emission outside the core could be produced by a minor merger that powers electron re-acceleration without disrupting the cool core. This discovery suggests that, under particular circumstances, both a mini and giant halo could co-exist in a single cluster, opening new perspectives for particle acceleration mechanisms in galaxy clusters.

We present results from LOFAR and GMRT observations of the galaxy cluster MACS$\,$J0717.5$+$3745. The cluster is undergoing a violent merger involving at least four sub-clusters, and it is known to host a radio halo. LOFAR observations reveal new sources of radio emission in the Intra-Cluster Medium: (i) a radio bridge that connects the cluster to a head-tail radio galaxy located along a filament of galaxies falling into the main cluster, (ii) a 1.9 Mpc radio arc, that is located North West of the main mass component, (iii) radio emission along the X-ray bar, that traces the gas in the X-rays South West of the cluster centre. We use deep GMRT observations at 608 MHz to constrain the spectral indices of these new radio sources, and of the emission that was already studied in the literature at higher frequency. We find that the spectrum of the radio halo and of the relic at LOFAR frequency follows the same power law as observed at higher frequencies. The radio bridge, the radio arc, and the radio bar all have steep spectra, which can be used to constrain the particle acceleration mechanisms. We argue that the radio bridge could be caused by the re-acceleration of electrons by shock waves that are injected along the filament during the cluster mass assembly. Despite the sensitivity reached by our observations, the emission from the radio halo does not trace the emission of the gas revealed by X-ray observations. We argue that this could be due to the difference in the ratio of kinetic over thermal energy of the intra-cluster gas, suggested by X-ray observations.

The ionosphere is the main driver of a series of systematic effects that limit our ability to explore the low frequency (<1 GHz) sky with radio interferometers. Its effects become increasingly important towards lower frequencies and are particularly hard to calibrate in the low signal-to-noise ratio regime in which low-frequency telescopes operate. In this paper we characterize and quantify the effect of ionospheric-induced systematic errors on astronomical interferometric radio observations at ultra-low frequencies (<100 MHz). We also provide guidelines for observations and data reduction at these frequencies with the Low Frequency Array (LOFAR) and future instruments such as the Square Kilometre Array (SKA). We derive the expected systematic error induced by the ionosphere. We compare our predictions with data from the Low Band Antenna (LBA) system of LOFAR. We show that we can isolate the ionospheric effect in LOFAR LBA data and that our results are compatible with satellite measurements, providing an independent way to measure the ionospheric total electron content (TEC). We show how the ionosphere also corrupts the correlated amplitudes through scintillations. We report values of the ionospheric structure function in line with the literature. The systematic errors on the phases of LOFAR LBA data can be accurately modelled as a sum of four effects (clock, ionosphere 1st, 2nd, and 3rd order). This greatly reduces the number of required calibration parameters, and therefore enables new efficient calibration strategies.

Extended, steep-spectrum radio synchrotron sources are pre-dominantly found in massive galaxy clusters as opposed to groups. LOFAR Two-Metre Sky Survey images have revealed a diffuse, ultra-steep spectrum radio source in the low-mass cluster Abell 1931. The source has a fairly irregular morphology with a largest linear size of about 550 kpc. The source is only seen in LOFAR observations at 143 MHz and GMRT observations at 325 MHz. The spectral index of the total source between 143 MHz and 325 MHz is $\alpha_{143}^{325} = -2.86 \pm 0.36$. The source remains invisible in Very Large Array (1-2 GHz) observations as expected given the spectral index. Chandra X-ray observations of the cluster revealed a bolometric luminosity of $L_X = (1.65 \pm 0.39) \times 10^{43}$ erg s$^{-1}$ and a temperature of $2.92_{-0.87}^{+1.89}$ keV which implies a mass of around $\sim 10^{14} M_{\odot}$. We conclude that the source is a remnant radio galaxy that has shut off around 200 Myr ago. The brightest cluster galaxy, a radio-loud elliptical galaxy, could be the source for this extinct source. Unlike remnant sources studied in the literature, our source has a steep spectrum at low radio frequencies. Studying such remnant radio galaxies at low radio frequencies is important for understanding the scarcity of such sources and their role in feedback processes.

We use 2.0 Msec of Chandra observations to investigate the cocoon shocks of Cygnus A and some implications for its lobes and jet. Measured shock Mach numbers vary in the range 1.18-1.66 around the cocoon. We estimate a total outburst energy of $\simeq 4.7\times10^{60}\rm\ erg$, with an age of $\simeq 2 \times 10^{7}\rm\ yr$. The average postshock pressure is found to be $8.6 \pm 0.3 \times 10^{-10}\rm\ erg\ cm^{-3}$, which agrees with the average pressure of the thin rim of compressed gas between the radio lobes and shocks, as determined from X-ray spectra. However, average rim pressures are found to be lower in the western lobe than in the eastern lobe by $\simeq 20\%$. Pressure estimates for hotspots A and D from synchrotron self-Compton models imply that each jet exerts a ram pressure $\gtrsim$ 3 times its static pressure, consistent with the positions of the hotspots moving about on the cocoon shock over time. A steady, one-dimensional flow model is used to estimate jet properties, finding mildly relativistic flow speeds within the allowed parameter range. Models in which the jet carries a negligible flux of rest mass are consistent with with the observed properties of the jets and hotspots. This favors the jets being light, implying that the kinetic power and momentum flux are carried primarily by the internal energy of the jet plasma rather than by its rest mass.

We use new and archival Chandra observations of Cygnus A, totalling $\sim$1.9 Ms, to investigate the distribution and temperature structure of gas lying within the projected extent of the cocoon shock and exhibiting a rib-like structure. We confirm that the X-rays are dominated by thermal emission with an average temperature of around 4 keV, and have discovered an asymmetry in the temperature gradient, with the southwestern part of the gas cooler than the rest by up to 2 keV. Pressure estimates suggest that the gas is a coherent structure of single origin located inside the cocoon, with a mass of roughly $2\times10^{10} M_{\odot}$. We conclude that the gas is debris resulting from disintegration of the cool core of the Cygnus A cluster after the passage of the jet during the early stages of the current epoch of activity. The 4 keV gas now lies on the central inside surface of the hotter cocoon rim. The temperature gradient could result from an offset between the centre of the cluster core and the Cygnus A host galaxy at the switch-on of current radio activity.

Merging galaxy clusters produce low Mach number shocks in the intracluster medium. These shocks can accelerate electrons to relativistic energies that are detectable at radio frequencies. MACS J0744.9+3927 is a massive ($M_{500} = (11.8 \pm 2.8) \times 10^{14} M_{\odot}$), high-redshift ($z=0.6976$) cluster where a Bullet-type merger is presumed to have taken place. Sunyaev-Zel'dovich maps from MUSTANG indicate that a shock, with Mach number $\mathcal{M} = 1.0-2.9$ and an extension of $\sim 200$ kpc, sits near the centre of the cluster. The shock is also detected as a brightness and temperature discontinuity in X-ray observations. To search for diffuse radio emission associated with the merger, we have imaged the cluster with the LOw Frequency ARray (LOFAR) at 120-165 MHz. Our LOFAR radio images reveal previously undetected AGN emission, but do not show clear cluster-scale diffuse emission in the form of a radio relic nor a radio halo. The region of the shock is on the western edge of AGN lobe emission from the brightest cluster galaxy. Correlating the flux of known shock-induced radio relics versus their size, we find that the radio emission overlapping the shocked region in MACS J0744.9+3927 is likely of AGN origin. We argue against the presence of a relic caused by diffusive shock acceleration and suggest that the shock is too weak to accelerate electrons from the intracluster medium.

Low-mass galaxies are subject to strong galactic outflows, in which cosmic rays may play an important role, they can be best traced with low-frequency radio continuum observations, which are less affected by spectral ageing. We present a study of the nearby star burst dwarf irregular galaxy IC 10 using observations at 140 MHz with the LOw-Frequency ARray (LOFAR), at 1580 MHz with the Very Large Array (VLA) and at 6200 MHz with the VLA and the 100-m Effelsberg telescope. We find that IC 10 has a low-frequency radio halo, which manifests itself as a second component (thick disc) in the minor axis profiles of the non-thermal radio continuum emission at 140 and 1580 MHz. These profiles are then fitted with 1D cosmic-ray transport models for pure diffusion and advection. We find that a diffusion model fits best, with a diffusion coefficient of $D=(0.4$-$0.8) \times 10^{26}(E/{\rm GeV})^{0.5}~{\rm cm^2\,s^{-1}}$, which is at least an order of magnitude smaller than estimates both from anisotropic diffusion and the diffusion length. In contrast, advection models, which cannot be ruled out due to the mild inclination, while providing poorer fits, result in advection speeds close to the escape velocity of $\approx$$50~\rm km\,s^{-1}$, as expected for a cosmic-ray driven wind. Our favoured model with an accelerating wind provides a self-consistent solution, where the magnetic field is in energy equipartition with both the warm neutral and warm ionized medium with an important contribution from cosmic rays. Consequently, cosmic rays can play a vital role for the launching of galactic winds in the disc--halo interface.

Feedback by radio-loud active galactic nuclei (AGN) in galaxy groups is not fully understood. Open questions include the duty cycle of the AGN, the spatial extent of the radio lobes, the effect they have on the intragroup medium, and the fate of the cosmic rays. We present the discovery of a 650 kpc-radio galaxy embedded in steep diffuse emission at $z = 0.18793 \pm 5 \times 10^{-5}$ located at the center of the galaxy group MaxBCG J199.31832+51.72503 using an observation from the LOFAR Two-meter Sky Survey (LoTSS) at the central frequency of 144 MHz. Subsequently, we performed a GMRT observation at the central frequency of 607 MHz to study the spectral properties of the source. The observations reveal a radio galaxy with a total radio power $P_{\rm tot, 1.4} \sim 2.1 \times 10^{24}$ W Hz$^{-1}$, exhibiting two asymmetrical jets and lobes. The derived spectral index map shows a steepening toward the inner regions and a steep-spectrum core region. We model the integrated radio spectrum, providing two possible interpretations: the radio source is evolved but still active or it is just at the end of its active phase. Finally, in the same field of view we have discovered Mpc-sized emission surrounding a close pair of AGN located at a redshift $z = 0.0587 \pm 2 \times 10^{-4}$ (SDSS J131544.56+521213.2 and SDSS J131543.99+521055.7) which could be a radio remnant source.

LOFAR observations at 144 MHz have revealed large-scale radio sources in the unrelaxed galaxy cluster Abell 1132. The cluster hosts diffuse radio emission on scales of $\sim$650 kpc near the cluster center and a head-tail (HT) radio galaxy, extending up to 1 Mpc, South of the cluster center. The central diffuse radio emission is not seen in NVSS, FIRST, WENSS, nor in C & D array VLA observations at 1.4 GHz, but is detected in our follow-up GMRT observations at 325 MHz. Using LOFAR and GMRT data, we determine the spectral index of the central diffuse emission to be $\alpha=-1.75\pm0.19$ ($S\propto\nu^{\alpha}$). We classify this emission as an ultra-steep spectrum radio halo and discuss the possible implications for the physical origin of radio halos. The HT radio galaxy shows narrow, collimated emission extending up to 1 Mpc and another 300 kpc of more diffuse, disturbed emission, giving a full projected linear size of 1.3 Mpc - classifying it as a giant radio galaxy (GRG) and making it the longest HT found to date. The head of the GRG coincides with an elliptical galaxy (SDSS J105851.01$+$564308.5) belonging to Abell 1132. In our LOFAR image, there appears to be a connection between the radio halo and the GRG. The turbulence that may have produced the halo may have also affected the tail of the GRG. In turn, the GRG may have provided seed electrons for the radio halo.

We present a study of AGN feedback at higher redshifts ($0.3<z<1.2$) using Sunyaev-Zel'dovich (SZ) selected samples of clusters from the South-Pole Telescope and Atacama Cosmology Telescope surveys. In contrast to studies of nearby systems, we do not find a separation between cooling flow clusters and non-cooling flow clusters based on the radio luminosity of the central radio source. This lack may be due to the increased incidence of galaxy-galaxy mergers at higher redshift that triggers AGN activity. In support of this scenario, we find evidence for evolution in the radio luminosity function of the central radio source: while the lower-luminosity sources do not evolve much, the higher-luminosity sources show a strong increase in the frequency of their occurrence at higher redshifts. We interpret this evolution as an increase in high-excitation radio galaxies (HERGs) in massive clusters at $z>0.6$, implying a transition from HERG-mode accretion to lower-power low-excitation radio galaxy (LERG)-mode accretion at intermediate redshifts. Additionally, we use local radio-to-jet power scaling relations to estimate feedback power and find that half of the cooling flow systems in our sample probably have enough heating to balance cooling. However, we postulate that the local relations are likely not well suited to predict feedback power in high-luminosity HERGs, as they are derived from samples composed mainly of lower-luminosity LERGs.

We present a new analysis of the widely used relation between cavity power and radio luminosity in clusters of galaxies with evidence for strong AGN feedback. We study the correlation at low radio frequencies using two new surveys - the First Alternative Data Release of the TIFR GMRT Sky Survey (TGSS ADR1) at 148 MHz and LOFAR's first all-sky survey, the Multifrequency Snapshot Sky Survey (MSSS) at 140 MHz. We find a scaling relation $P_{\rm cav} \propto L_{148}^{\beta}$, with a logarithmic slope of $\beta = 0.51 \pm 0.14$, which is in good agreement with previous results based on data at 327 MHz. The large scatter present in this correlation confirms the conclusion reached at higher frequencies that the total radio luminosity at a single frequency is a poor predictor of the total jet power. We show that including measurements at 148 MHz alone is insufficient to reliably compute the bolometric radio luminosity and reduce the scatter in the correlation. For a subset of four well-resolved sources, we examine the detected extended structures at low frequencies and compare with the morphology known from higher frequency images and Chandra X-ray maps. In Perseus we discuss details in the structures of the radio mini-halo, while in the 2A 0335+096 cluster we observe new diffuse emission associated with multiple X-ray cavities and likely originating from past activity. For A2199 and MS 0735.6+7421, we confirm that the observed low-frequency radio lobes are confined to the extents known from higher frequencies. This new low-frequency analysis highlights the fact that existing cavity power to radio luminosity relations are based on a relatively narrow range of AGN outburst ages. We discuss how the correlation could be extended using low frequency data from the LOFAR Two-metre Sky Survey (LoTSS) in combination with future, complementary deeper X-ray observations.

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 present LOFAR High-Band Array (HBA) observations of the Herschel-ATLAS North Galactic Pole survey area. The survey we have carried out, consisting of four pointings covering around 142 square degrees of sky in the frequency range 126--173 MHz, does not provide uniform noise coverage but otherwise is representative of the quality of data to be expected in the planned LOFAR wide-area surveys, and has been reduced using recently developed `facet calibration' methods at a resolution approaching the full resolution of the datasets ($\sim 10 \times 6$ arcsec) and an rms off-source noise that ranges from 100 $\mu$Jy beam$^{-1}$ in the centre of the best fields to around 2 mJy beam$^{-1}$ at the furthest extent of our imaging. We describe the imaging, cataloguing and source identification processes, and present some initial science results based on a 5-$\sigma$ source catalogue. These include (i) an initial look at the radio/far-infrared correlation at 150 MHz, showing that many Herschel sources are not yet detected by LOFAR; (ii) number counts at 150 MHz, including, for the first time, observational constraints on the numbers of star-forming galaxies; (iii) the 150-MHz luminosity functions for active and star-forming galaxies, which agree well with determinations at higher frequencies at low redshift, and show strong redshift evolution of the star-forming population; and (iv) some discussion of the implications of our observations for studies of radio galaxy life cycles.

LOFAR, the Low-Frequency Array, is a powerful new radio telescope operating between 10 and 240 MHz. LOFAR allows detailed sensitive high-resolution studies of the low-frequency radio sky. At the same time LOFAR also provides excellent short baseline coverage to map diffuse extended emission. However, producing high-quality deep images is challenging due to the presence of direction dependent calibration errors, caused by imperfect knowledge of the station beam shapes and the ionosphere. Furthermore, the large data volume and presence of station clock errors present additional difficulties. In this paper we present a new calibration scheme, which we name facet calibration, to obtain deep high-resolution LOFAR High Band Antenna images using the Dutch part of the array. This scheme solves and corrects the direction dependent errors in a number of facets that cover the observed field of view. Facet calibration provides close to thermal noise limited images for a typical 8 hr observing run at $\sim$ 5arcsec resolution, meeting the specifications of the LOFAR Tier-1 northern survey.

We present deep LOFAR observations between 120-181 MHz of the "Toothbrush" (RX J0603.3+4214), a cluster that contains one of the brightest radio relic sources known. Our LOFAR observations exploit a new and novel calibration scheme to probe 10 times deeper than any previous study in this relatively unexplored part of the spectrum. The LOFAR observations, when combined with VLA, GMRT, and Chandra X-ray data, provide new information about the nature of cluster merger shocks and their role in re-accelerating relativistic particles. We derive a spectral index of $\alpha = -0.8 \pm 0.1$ at the northern edge of the main radio relic, steepening towards the south to $\alpha \approx - 2$. The spectral index of the radio halo is remarkably uniform ($\alpha = -1.16$, with an intrinsic scatter of $\leq 0.04$). The observed radio relic spectral index gives a Mach number of $\mathcal{M} = 2.8^{+0.5}_{-0.3}$, assuming diffusive shock acceleration (DSA). However, the gas density jump at the northern edge of the large radio relic implies a much weaker shock ($\mathcal{M} \approx 1.2$, with an upper limit of $\mathcal{M} \approx 1.5$). The discrepancy between the Mach numbers calculated from the radio and X-rays can be explained if either (i) the relic traces a complex shock surface along the line of sight, or (ii) if the radio relic emission is produced by a re-accelerated population of fossil particles from a radio galaxy. Our results highlight the need for additional theoretical work and numerical simulations of particle acceleration and re-acceleration at cluster merger shocks.

Steep spectrum radio sources associated with active galactic nuclei (AGN) may contain remnants of past AGN activity episodes. Novel instruments like the LOw Frequency ARray (LOFAR) are enabling studies of these fascinating structures to be made at tens to hundreds of MHz with sufficient resolution to analyse their complex morphology. Our goal is to characterize the integrated and resolved spectral properties of VLSS J1431+1331 and estimate source ages based on synchrotron radio emission models, thus putting constraints on the AGN duty cycle. Using a broad spectral coverage, we have derived spectral and curvature maps, and used synchrotron ageing models to determine the time elapsed from the last time the source plasma was energized. We used LOFAR, Giant Metrewave Radio Telescope (GMRT) and Jansky Very Large Array (VLA) data. Based on our ageing analysis, we infer that the AGN that created this source currently has very low levels of activity or that it is switched off. The derived ages for the larger source component range from around 60 to 130 Myr, hinting that the AGN activity decreased or stopped around 60 Myr ago. Our analysis suggests that VLSS J1431.8+1331 is an intriguing, two-component source. The larger component seems to host a faint radio core, suggesting that the source may be an AGN radio relic. The spectral index we observe from the smaller component is distinctly flatter at lower frequencies than the spectral index of the larger component, suggesting the possibility that the smaller component may be a shocked plasma bubble. From the integrated source spectrum, we deduce that its shape and slope can be used as tracers of the activity history of this type of steep spectrum radio source.

The Evolutionary Map of the Universe (EMU) is a proposed radio continuum survey of the Southern Hemisphere up to declination +30 deg., with the Australian Square Kilometre Array Pathfinder (ASKAP). EMU will use an automated source identification and measurement approach that is demonstrably optimal, to maximise the reliability, utility and robustness of the resulting radio source catalogues. As part of the process of achieving this aim, a "Data Challenge" has been conducted, providing international teams the opportunity to test a variety of source finders on a set of simulated images. The aim is to quantify the accuracy of existing automated source finding and measurement approaches, and to identify potential limitations. The Challenge attracted nine independent teams, who tested eleven different source finding tools. In addition, the Challenge initiators also tested the current ASKAPsoft source-finding tool to establish how it could benefit from incorporating successful features of the other tools. Here we present the results of the Data Challenge, identifying the successes and limitations for this broad variety of the current generation of radio source finding tools. As expected, most finders demonstrate completeness levels close to 100% at 10sigma dropping to levels around 10% by 5sigma. The reliability is typically close to 100% at 10sigma, with performance to lower sensitivities varying greatly between finders. All finders demonstrate the usual trade-off between completeness and reliability, whereby maintaining a high completeness at low signal-to-noise comes at the expense of reduced reliability, and vice-versa. We conclude with a series of recommendations for improving the performance of the ASKAPsoft source-finding tool.

We present the Multifrequency Snapshot Sky Survey (MSSS), the first northern-sky LOFAR imaging survey. In this introductory paper, we first describe in detail the motivation and design of the survey. Compared to previous radio surveys, MSSS is exceptional due to its intrinsic multifrequency nature providing information about the spectral properties of the detected sources over more than two octaves (from 30 to 160 MHz). The broadband frequency coverage, together with the fast survey speed generated by LOFAR's multibeaming capabilities, make MSSS the first survey of the sort anticipated to be carried out with the forthcoming Square Kilometre Array (SKA). Two of the sixteen frequency bands included in the survey were chosen to exactly overlap the frequency coverage of large-area Very Large Array (VLA) and Giant Metrewave Radio Telescope (GMRT) surveys at 74 MHz and 151 MHz respectively. The survey performance is illustrated within the "MSSS Verification Field" (MVF), a region of 100 square degrees centered at J2000 (RA,Dec)=(15h,69deg). The MSSS results from the MVF are compared with previous radio survey catalogs. We assess the flux and astrometric uncertainties in the catalog, as well as the completeness and reliability considering our source finding strategy. We determine the 90% completeness levels within the MVF to be 100 mJy at 135 MHz with 108" resolution, and 550 mJy at 50 MHz with 166" resolution. Images and catalogs for the full survey, expected to contain 150,000-200,000 sources, will be released to a public web server. We outline the plans for the ongoing production of the final survey products, and the ultimate public release of images and source catalogs.

PSR B0823+26, a 0.53-s radio pulsar, displays a host of emission phenomena over timescales of seconds to (at least) hours, including nulling, subpulse drifting, and mode-changing. Studying pulsars like PSR B0823+26 provides further insight into the relationship between these various emission phenomena and what they might teach us about pulsar magnetospheres. Here we report on the LOFAR discovery that PSR B0823+26 has a weak and sporadically emitting 'quiet' (Q) emission mode that is over 100 times weaker (on average) and has a nulling fraction forty-times greater than that of the more regularly-emitting 'bright' (B) mode. Previously, the pulsar has been undetected in the Q-mode, and was assumed to be nulling continuously. PSR B0823+26 shows a further decrease in average flux just before the transition into the B-mode, and perhaps truly turns off completely at these times. Furthermore, simultaneous observations taken with the LOFAR, Westerbork, Lovell, and Effelsberg telescopes between 110 MHz and 2.7 GHz demonstrate that the transition between the Q-mode and B-mode occurs within one single rotation of the neutron star, and that it is concurrent across the range of frequencies observed.

Using observations obtained with the LOw Fequency ARray (LOFAR), the Westerbork Synthesis Radio Telescope (WSRT) and archival Very Large Array (VLA) data, we have traced the radio emission to large scales in the complex source 4C 35.06 located in the core of the galaxy cluster Abell 407. At higher spatial resolution (~4"), the source was known to have two inner radio lobes spanning 31 kpc and a diffuse, low-brightness extension running parallel to them, offset by about 11 kpc (in projection). At 62 MHz, we detect the radio emission of this structure extending out to 210 kpc. At 1.4 GHz and intermediate spatial resolution (~30"), the structure appears to have a helical morphology. We have derived the characteristics of the radio spectral index across the source. We show that the source morphology is most likely the result of at least two episodes of AGN activity separated by a dormant period of around 35 Myr. The AGN is hosted by one of the galaxies located in the cluster core of Abell 407. We propose that it is intermittently active as it moves in the dense environment in the cluster core. Using LOFAR, we can trace the relic plasma from that episode of activity out to greater distances from the core than ever before. Using the the WSRT, we detect HI in absorption against the center of the radio source. The absorption profile is relatively broad (FWHM of 288 km/s), similar to what is found in other clusters. Understanding the duty cycle of the radio emission as well as the triggering mechanism for starting (or restarting) the radio-loud activity can provide important constraints to quantify the impact of AGN feedback on galaxy evolution. The study of these mechanisms at low frequencies using morphological and spectral information promises to bring new important insights in this field.

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.

Faint undetected sources of radio-frequency interference (RFI) might become visible in long radio observations when they are consistently present over time. Thereby, they might obstruct the detection of the weak astronomical signals of interest. This issue is especially important for Epoch of Reionisation (EoR) projects that try to detect the faint redshifted HI signals from the time of the earliest structures in the Universe. We explore the RFI situation at 30-163 MHz by studying brightness histograms of visibility data observed with LOFAR, similar to radio-source-count analyses that are used in cosmology. An empirical RFI distribution model is derived that allows the simulation of RFI in radio observations. The brightness histograms show an RFI distribution that follows a power-law distribution with an estimated exponent around -1.5. With several assumptions, this can be explained with a uniform distribution of terrestrial radio sources whose radiation follows existing propagation models. Extrapolation of the power law implies that the current LOFAR EoR observations should be severely RFI limited if the strength of RFI sources remains strong after time integration. This is in contrast with actual observations, which almost reach the thermal noise and are thought not to be limited by RFI. Therefore, we conclude that it is unlikely that there are undetected RFI sources that will become visible in long observations. Consequently, there is no indication that RFI will prevent an EoR detection with LOFAR.

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

Mechanical feedback via Active Galactic Nuclei (AGN) jets in the centres of galaxy groups and clusters is a crucial ingredient in current models of galaxy formation and cluster evolution. Jet feedback is believed to regulate gas cooling and thus star formation in the most massive galaxies, but a robust physical understanding of this feedback mode is currently lacking. The large collecting area, excellent spectral resolution and high spatial resolution of Athena+ will provide the breakthrough diagnostic ability necessary to develop this understanding, via: (1) the first kinematic measurements on relevant spatial scales of the hot gas in galaxy, group and cluster haloes as it absorbs the impact of AGN jets, and (2) vastly improved ability to map thermodynamic conditions on scales well-matched to the jets, lobes and gas disturbances produced by them. Athena+ will therefore determine for the first time how jet energy is dissipated and distributed in group and cluster gas, and how a feedback loop operates in group/cluster cores to regulate gas cooling and AGN fuelling. Athena+ will also establish firmly the cumulative impact of powerful radio galaxies on the evolution of baryons from the epoch of group/cluster formation to the present day.

LOFAR, the LOw-Frequency ARray, is a new-generation radio interferometer constructed in the north of the Netherlands and across europe. Utilizing a novel phased-array design, LOFAR covers the largely unexplored low-frequency range from 10-240 MHz and provides a number of unique observing capabilities. Spreading out from a core located near the village of Exloo in the northeast of the Netherlands, a total of 40 LOFAR stations are nearing completion. A further five stations have been deployed throughout Germany, and one station has been built in each of France, Sweden, and the UK. Digital beam-forming techniques make the LOFAR system agile and allow for rapid repointing of the telescope as well as the potential for multiple simultaneous observations. With its dense core array and long interferometric baselines, LOFAR achieves unparalleled sensitivity and angular resolution in the low-frequency radio regime. The LOFAR facilities are jointly operated by the International LOFAR Telescope (ILT) foundation, as an observatory open to the global astronomical community. LOFAR is one of the first radio observatories to feature automated processing pipelines to deliver fully calibrated science products to its user community. LOFAR's new capabilities, techniques and modus operandi make it an important pathfinder for the Square Kilometre Array (SKA). We give an overview of the LOFAR instrument, its major hardware and software components, and the core science objectives that have driven its design. In addition, we present a selection of new results from the commissioning phase of this new radio observatory.

The Chandra X-ray Observatory has revealed X-ray bubbles in the intracluster medium (ICM) of many nearby cooling flow clusters. The bubbles trace feedback that is thought to couple the central active galactic nucleus (AGN) to the ICM, helping to stabilize cooling flows and govern the evolution of massive galaxies. However, the prevalence and duty cycle of such AGN outbursts is not well understood. To this end, we study how cooling is balanced by bubble heating for complete samples of clusters (the Brightest 55 clusters of galaxies, hereafter B55, and the HIghest X-ray FLUx Galaxy Cluster Sample, HIFLUGCS). We find that the radio luminosity of the central galaxy only exceeds 2.5 x 10^30 erg s^-1 Hz^-1 in cooling flow clusters. This result implies a connection between the central radio source and the ICM, as expected if AGN feedback is operating. Additionally, we find a duty cycle for radio mode feedback, the fraction of time that a system possesses bubbles inflated by its central radio source, of > 69 per cent for B55 and > 63 per cent for HIFLUGCS. These duty cycles are lower limits since some bubbles are likely missed in existing images. We used simulations to constrain the bubble power that might be present and remain undetected in the cooling flow systems without detected bubbles. Among theses systems, almost all could have significant bubble power. Therefore, our results imply that the duty cycle of AGN outbursts with the potential to heat the gas significantly in cooling flow clusters is at least 60 per cent and could approach 100 per cent.

We report on an analysis of new Chandra data of the galaxy group HCG 62, well known for possessing cavities in its intragroup medium (IGM) that were inflated by the radio lobes of its central active galactic nucleus (AGN). With the new data, a factor of three deeper than previous Chandra data, we re-examine the energetics of the cavities and determine new constraints on their contents. We confirm that the ratio of radiative to mechanical power of the AGN outburst that created the cavities is less than 10^-4, among the lowest of any known cavity system, implying that the relativistic electrons in the lobes can supply only a tiny fraction of the pressure required to support the cavities. This finding implies additional pressure support in the lobes from heavy particles (e.g., protons) or thermal gas. Using spectral fits to emission in the cavities, we constrain any such volume-filling thermal gas to have a temperature kT > 4.3 keV. For the first time, we detect X-ray emission from the central AGN, with a luminosity of L(2-10 keV) = (1.1 +/- 0.4) x 10^39 erg s^-1 and properties typical of a low-luminosity AGN. Lastly, we report evidence for a recent merger from the surface brightness, temperature, and metallicity structure of the IGM.

M87 is a giant elliptical galaxy located in the centre of the Virgo cluster, which harbours a supermassive black hole of mass 6.4x10^9 M_sun, whose activity is responsible for the extended (80 kpc) radio lobes that surround the galaxy. The energy generated by matter falling onto the central black hole is ejected and transferred to the intra-cluster medium via a relativistic jet and morphologically complex systems of buoyant bubbles, which rise towards the edges of the extended halo. Here we present the first observations made with the new Low-Frequency Array (LOFAR) of M87 at frequencies down to 20 MHz. Images of M87 were produced at low radio frequencies never explored before at these high spatial resolution and dynamic range. To disentangle different synchrotron models and place constraints on source magnetic field, age and energetics, we also performed a detailed spectral analysis of M87 extended radio-halo using these observations together with archival data. We do not find any sign of new extended emissions; on the contrary the source appears well confined by the high pressure of the intra-cluster medium. A continuous injection of relativistic electrons is the model that best fits our data, and provides a scenario in which the lobes are still supplied by fresh relativistic particles from the active galactic nuclei. We suggest that the discrepancy between the low-frequency radio-spectral slope in the core and in the halo implies a strong adiabatic expansion of the plasma as soon as it leaves the core area. The extended halo has an equipartition magnetic field strength of ~10 uG, which increases to ~13 uG in the zones where the particle flows are more active. The continuous injection model for synchrotron ageing provides an age for the halo of ~40 Myr, which in turn provides a jet kinetic power of 6-10x10^44 erg/s.

Deep radio observations of galaxy clusters have revealed the existence of diffuse radio sources related to the presence of relativistic electrons and weak magnetic fields in the intracluster volume. The role played by this non-thermal intracluster component on the thermodynamical evolution of galaxy clusters is debated, with important implications for cosmological and astrophysical studies of the largest gravitationally bound structures of the Universe. The low surface brightness and steep spectra of diffuse cluster radio sources make them more easily detectable at low-frequencies. LOFAR is the first instrument able to detect diffuse radio emission in hundreds of massive galaxy clusters up to their formation epoch. We present the first observations of clusters imaged by LOFAR and the huge perspectives opened by this instrument for non-thermal cluster studies.

Aims: This paper discusses the spectral occupancy for performing radio astronomy with the Low-Frequency Array (LOFAR), with a focus on imaging observations. Methods: We have analysed the radio-frequency interference (RFI) situation in two 24-h surveys with Dutch LOFAR stations, covering 30-78 MHz with low-band antennas and 115-163 MHz with high-band antennas. This is a subset of the full frequency range of LOFAR. The surveys have been observed with a 0.76 kHz / 1 s resolution. Results: We measured the RFI occupancy in the low and high frequency sets to be 1.8% and 3.2% respectively. These values are found to be representative values for the LOFAR radio environment. Between day and night, there is no significant difference in the radio environment. We find that lowering the current observational time and frequency resolutions of LOFAR results in a slight loss of flagging accuracy. At LOFAR's nominal resolution of 0.76 kHz and 1 s, the false-positives rate is about 0.5%. This rate increases approximately linearly when decreasing the data frequency resolution. Conclusions: Currently, by using an automated RFI detection strategy, the LOFAR radio environment poses no perceivable problems for sensitive observing. It remains to be seen if this is still true for very deep observations that integrate over tens of nights, but the situation looks promising. Reasons for the low impact of RFI are the high spectral and time resolution of LOFAR; accurate detection methods; strong filters and high receiver linearity; and the proximity of the antennas to the ground. We discuss some strategies that can be used once low-level RFI starts to become apparent. It is important that the frequency range of LOFAR remains free of broadband interference, such as DAB stations and windmills.

Abell 2256 is one of the best known examples of a galaxy cluster hosting large-scale diffuse radio emission that is unrelated to individual galaxies. It contains both a giant radio halo and a relic, as well as a number of head-tail sources and smaller diffuse steep-spectrum radio sources. The origin of radio halos and relics is still being debated, but over the last years it has become clear that the presence of these radio sources is closely related to galaxy cluster merger events. Here we present the results from the first LOFAR Low band antenna (LBA) observations of Abell 2256 between 18 and 67 MHz. To our knowledge, the image presented in this paper at 63 MHz is the deepest ever obtained at frequencies below 100 MHz in general. Both the radio halo and the giant relic are detected in the image at 63 MHz, and the diffuse radio emission remains visible at frequencies as low as 20 MHz. The observations confirm the presence of a previously claimed ultra-steep spectrum source to the west of the cluster center with a spectral index of -2.3 \pm 0.4 between 63 and 153 MHz. The steep spectrum suggests that this source is an old part of a head-tail radio source in the cluster. For the radio relic we find an integrated spectral index of -0.81 \pm 0.03, after removing the flux contribution from the other sources. This is relatively flat which could indicate that the efficiency of particle acceleration at the shock substantially changed in the last \sim 0.1 Gyr due to an increase of the shock Mach number. In an alternative scenario, particles are re-accelerated by some mechanism in the downstream region of the shock, resulting in the relatively flat integrated radio spectrum. In the radio halo region we find indications of low-frequency spectral steepening which may suggest that relativistic particles are accelerated in a rather inhomogeneous turbulent region.

We present an analysis of deep multiwavelength data for z ~ 0.3-3 starburst galaxies selected by their 70 um emission in the Extended-Chandra Deep Field-South and Extended Groth Strip. We identify active galactic nuclei (AGNs) in these infrared sources through their X-ray emission and quantify the fraction that host an AGN. We find that the fraction depends strongly on both the mid-infrared color and rest-frame mid-infrared luminosity of the source, rising to ~ 50-70% at the warmest colors and highest mid-infrared luminosities (corresponding to ultraluminous infrared galaxies), similar to the trends found locally. Additionally, we find that the AGN fraction depends strongly on the star formation rate of the host galaxy (inferred from the observed 70 um luminosity after subtracting the estimated AGN contribution), particularly for more luminous AGNs (L_X > 10^43 erg/s). At the highest star formation rates (~ 1000 M_sun/yr), the fraction of galaxies with an X-ray detected AGN rises to ~ 30%, roughly consistent with that found in high-redshift submillimeter galaxies. Assuming that the AGN fraction is driven by the star formation rate (rather than stellar mass or redshift, for which our sample is largely degenerate), this result implies that the duty cycle of luminous AGN activity increases with the star formation rate of the host galaxy: specifically, we find that luminous X-ray detected AGNs are at least ~ 5-10 times more common in systems with high star formation rates (> 300 M_sun/yr) than in systems with lower star formation rates (< 30 M_sun/yr). Lastly, we investigate the ratio between the supermassive black hole accretion rate and the growth rate of the host galaxy and find that, for sources with detected AGNs and star formation, this ratio in distant starbursts agrees well with that expected from the local scaling relation assuming the black holes and bulges grew at the same epoch.

(abridged) We identify a numerically significant population of heavily obscured AGNs at z~0.5-1 in the Chandra Deep Field-South (CDF-S) and Extended Chandra Deep Field-South by selecting 242 X-ray undetected objects with infrared-based star formation rates (SFRs) substantially higher (a factor of 3.2 or more) than their SFRs determined from the UV after correcting for dust extinction. An X-ray stacking analysis of 23 candidates in the central CDF-S region using the 4 Ms Chandra data reveals a hard X-ray signal with an effective power-law photon index of Gamma=0.6_{-0.4}^{+0.3}, indicating a significant contribution from obscured AGNs. Based on Monte Carlo simulations, we conclude that 74+-25% of the selected galaxies host obscured AGNs, within which ~95% are heavily obscured and ~80% are Compton-thick (CT; NH>1.5x10^{24} cm^{-2}). The heavily obscured objects in our sample are of moderate intrinsic X-ray luminosity [ ~ (0.9-4)x10^{42} erg/s in the 2-10 keV band]. The space density of the CT AGNs is (1.6+-0.5)x10^{-4} Mpc^{-3}. The z~0.5-1 CT objects studied here are expected to contribute ~1% of the total XRB flux in the 10-30 keV band, and they account for ~5-15% of the emission in this energy band expected from all CT AGNs according to population-synthesis models. In the 6--8 keV band, the stacked signal of the 23 heavily obscured candidates accounts for <5% of the unresolved XRB flux, while the unresolved ~25% of the XRB in this band can probably be explained by a stacking analysis of the X-ray undetected optical galaxies in the CDF-S (a 2.5 sigma stacked signal). We discuss prospects to identify such heavily obscured objects using future hard X-ray observatories.

At very low frequencies, the new pan-European radio telescope LOFAR is opening the last unexplored window of the electromagnetic spectrum for astrophysical studies. The revolutionary APERTIF phased arrays that are about to be installed on the Westerbork radio telescope (WSRT) will dramatically increase the survey speed for the WSRT. Combined surveys with these two facilities will deeply chart the northern sky over almost two decades in radio frequency from \sim 15 up to 1400 MHz. Here we briefly describe some of the capabilities of these new facilities and what radio surveys are planned to study fundamental issues related the formation and evolution of galaxies and clusters of galaxies. In the second part we briefly review some recent observational results directly showing that diffuse radio emission in clusters traces shocks due to cluster mergers. As these diffuse radio sources are relatively bright at low frequencies, LOFAR should be able to detect thousands of such sources up to the epoch of cluster formation. This will allow addressing many question about the origin and evolution of shocks and magnetic fields in clusters. At the end we briefly review some of the first and very preliminary LOFAR results on clusters.

The Low Frequency Array (LOFAR) is under construction in the Netherlands and in several surrounding European countries. In this contribution, we describe the layout and design of the telescope, with a particular emphasis on the imaging characteristics of the array when used in its "standard imaging" mode. After briefly reviewing the calibration and imaging software used for LOFAR image processing, we show some recent results from the ongoing imaging commissioning efforts. We conclude by summarizing future prospects for the use of LOFAR in observing the little-explored low frequency Universe.

We use the 4Ms CDF-S survey to place direct X-ray constraints on the ubiquity of z~2 heavily obscured AGNs in K<22 BzK galaxies. Forty seven of the 222 BzK galaxies in the central region of the CDF-S are detected at X-ray energies, 11 of which have hard X-ray spectral slopes (Gamma<1) indicating the presence of heavily obscured AGN activity. The other 36 X-ray detected BzK galaxies appear to be relatively unobscured AGNs and starburst galaxies; we use X-ray variability analyses over a rest-frame baseline of ~3 years to further confirm the presence of AGN activity in many of these systems. The majority (7 out of 11) of the heavily obscured AGNs have excess IR emission over that expected from star formation (termed "IR-excess galaxies"). However, we find that X-ray detected heavily obscured AGNs only comprise ~25% of the IR-excess galaxy population, which is otherwise composed of relatively unobscured AGNs and starburst galaxies. We find that the typical X-ray spectrum of the heavily obscured AGNs is better characterized by a pure reflection model than an absorbed power-law model, suggesting extreme Compton-thick absorption in some systems. We verify this result by producing a composite rest-frame 2-20 keV spectrum, which has a similar shape as a reflection-dominated X-ray spectrum and reveals an emission feature at rest-frame energy ~6.4 keV, likely to be due to Fe K. These heavily obscured AGNs are likely to be the distant analogs of the reflection-dominated AGNs recently identified at z~0 with >10 keV observatories. On the basis of these analyses we estimate the space density for typical (intrinsic X-ray luminosities of L_X>1E43 erg/s) heavily obscured and Compton-thick AGNs at z~2. Our space-density constraints are conservative lower limits but they are already consistent with the range of predictions from X-ray background models.

The Low Frequency Array (LOFAR) will operate between 10 and 250 MHz, and will observe the low frequency Universe to an unprecedented sensitivity and angular resolution. The construction and commissioning of LOFAR is well underway, with over 27 of the Dutch stations and five International stations routinely performing both single-station and interferometric observations over the frequency range that LOFAR is anticipated to operate at. Here, we summarize the capabilities of LOFAR and report on some of the early commissioning imaging of Cygnus A.

[abridged] We present point-source catalogs for the 4Ms Chandra Deep Field-South (CDF-S), which is the deepest Chandra survey to date and covers an area of 464.5 arcmin^2. We provide a main source catalog, which contains 740 X-ray point sources that are detected with wavdetect at a false-positive probability threshold of 1E-5 and also satisfy a binomial-probability source-selection criterion of P<0.004; this approach is designed to maximize the number of reliable sources detected. A total of 300 main-catalog sources are new compared to the previous 2Ms CDF-S main-catalog sources. We also provide a supplementary catalog, which consists of 36 sources that are detected with wavdetect at 1E-5, satisfy 0.004< P<0.1, and have an optical counterpart with R<24. Multiwavelength identifications, basic optical/infrared/radio photometry, and spectroscopic/photometric redshifts are provided for the X-ray sources. Basic analyses of the X-ray and multiwavelength properties of the sources indicate that >75% of the main-catalog sources are AGNs; of the 300 new main-catalog sources, about 35% are likely normal and starburst galaxies, reflecting the rise of normal and starburst galaxies at the very faint flux levels uniquely accessible to the 4Ms CDF-S. Near the center of the 4Ms CDF-S, the observed AGN and galaxy source densities have reached ~9800 and 6900 per square degree, respectively. The 4 Ms CDF-S reaches on-axis flux limits of ~9.1E-18 and 5.5E-17 erg/cm^2/s for the soft and hard bands, respectively. An increase in the CDF-S exposure by a factor of ~2-2.5 would provide further significant gains and probe key unexplored discovery space.

Utilizing $\sim 50$ ks of Chandra X-ray Observatory imaging, we present an analysis of the intracluster medium (ICM) and cavity system in the galaxy cluster RBS 797. In addition to the two previously known cavities in the cluster core, the new and deeper X-ray image has revealed additional structure associated with the active galactic nucleus (AGN). The surface brightness decrements of the two cavities are unusually large, and are consistent with elongated cavities lying close to our line-of-sight. We estimate a total AGN outburst energy and mean jet power of $\approx 3 - 6 \times 10^{60}$ erg and $\approx 3 - 6 \times 10^{45}$ erg s$^{-1}$, respectively, depending on the assumed geometrical configuration of the cavities. Thus, RBS 797 is apparently among the the most powerful AGN outbursts known in a cluster. The average mass accretion rate needed to power the AGN by accretion alone is $\sim 1 M_{\odot}$ yr$^{-1}$. We show that accretion of cold gas onto the AGN at this level is plausible, but that Bondi accretion of the hot atmosphere is probably not. The BCG harbors an unresolved, non-thermal nuclear X-ray source with a bolometric luminosity of $\approx 2 \times 10^{44}$ erg s$^{-1}$. The nuclear emission is probably associated with a rapidly-accreting, radiatively inefficient accretion flow. We present tentative evidence that star formation in the BCG is being triggered by the radio jets and suggest that the cavities may be driving weak shocks ($M \sim 1.5$) into the ICM, similar to the process in the galaxy cluster MS 0735.6+7421.

We present the results of a program to acquire high-quality optical spectra of X-ray sources detected in the E-CDF-S and its central area. New spectroscopic redshifts are measured for 283 counterparts to Chandra sources with deep exposures (t~2-9 hr per pointing) using multi-slit facilities on both the VLT and Keck thus bringing the total number of spectroscopically-identified X-ray sources to over 500 in this survey field. We provide a comprehensive catalog of X-ray sources detected in the E-CDF-S including the optical and near-infrared counterparts, and redshifts (both spectroscopic and photometric) that incorporate published spectroscopic catalogs thus resulting in a final sample with a high fraction (80%) of X-ray sources having secure identifications. We demonstrate the remarkable coverage of the Lx-z plane now accessible from our data while emphasizing the detection of AGNs that contribute to the faint end of the luminosity function at 1.5<z<3. Our redshift catalog includes 17 type 2 QSOs that significantly increases such samples (2x). With our deepest VIMOS observation, we identify "elusive" optically-faint galaxies (R~25) at z~2-3 based upon the detection of interstellar absorption lines; we highlight one such case, an absorption-line galaxy at z=3.208 having no obvious signs of an AGN in its optical spectrum. In addition, we determine distances to eight galaxy groups with extended X-ray emission. Finally, we measure the physical extent of known large-scale structures (z~0.7) evident in the CDF-S. While a thick sheet (radial size of 67.7 Mpc) at z~0.67 extends over the full field, the z~0.73 structure is thin (18.8 Mpc) and filamentary as traced by both AGNs and galaxy groups. In the appendix, we provide spectroscopic redshifts for 49 counterparts to fainter X-ray sources detected only in the 1 and 2 Ms catalogs, and 48 VLA radio sources not detected in X-rays.

One of the science drivers of the new Low Frequency Array (LOFAR) is large-area surveys of the low-frequency radio sky. Realizing this goal requires automated processing of the interferometric data, such that fully calibrated images are produced by the system during survey operations. The LOFAR Imaging Pipeline is the tool intended for this purpose, and is now undergoing significant commissioning work. The pipeline is now functional as an automated processing chain. Here we present several recent LOFAR images that have been produced during the still ongoing commissioning period. These early LOFAR images are representative of some of the science goals of the commissioning team members.

[Abridged] We extend color-magnitude relations for moderate-luminosity X-ray AGN hosts and non-AGN galaxies through the galaxy formation epoch in the Chandra Deep Fields. We utilized analyses of color-magnitude diagrams (CMDs) to assess the role of moderate-luminosity AGNs in galaxy evolution. First, we confirm some previous results and extend them to higher redshifts, e.g., there is no apparent color bimodality for AGN hosts from z~0-2, but non-AGN galaxy color bimodality exists up to z~3; most AGNs reside in massive hosts and the AGN fraction rises strongly toward higher stellar mass up to z~2-3; and the colors of both AGN hosts and non-AGN galaxies become redder as the stellar mass increases up to z~2-3. Second, we point out that it is critical to use mass-matched samples to examine color-magnitude relations of AGN hosts and non-AGN galaxies. We show that for mass-matched samples up to z~2-3, AGN hosts lie in the same CMD region as non-AGN galaxies; i.e., there is no specific clustering of AGN hosts around the red sequence, the top of the blue cloud, or the green valley in between. The AGN fraction (~10%) is mostly independent of host-galaxy color, providing an indication of the duty cycle of SMBH growth in typical massive galaxies. These results are in contrast to those obtained with non-mass-matched samples where there is apparent AGN clustering in the CMD and the AGN fraction generally increases as the color becomes redder. We also find, for mass-matched samples, that the SFRs of AGN hosts are typically a factor of ~2-3 larger than those of non-AGN galaxies at z~0-1, whereas this difference diminishes at z~1-3. Our results can be reasonably explained by two main ingredients, color-mass correlation and passive or secular evolution of galaxies; and thus tightly constrain any effects from moderate-luminosity AGN feedback upon color-magnitude properties over the ~80% of cosmic time.

[Abridged] We present reliable multiwavelength identifications and high-quality photometric redshifts for the 462 X-ray sources in the ~2 Ms Chandra Deep Field-South. Source identifications are carried out using deep optical-to-radio multiwavelength catalogs, and are then combined to create lists of primary and secondary counterparts for the X-ray sources. We identified reliable counterparts for 446 (96.5%) of the X-ray sources, with an expected false-match probability of ~6.2%. A likelihood-ratio method is used for source matching, which effectively reduces the false-match probability at faint magnitudes compared to a simple error-circle matching method. We construct a master photometric catalog for the identified X-ray sources including up to 42 bands of UV-to-infrared data, and then calculate their photometric redshifts (photo-z's). The reliability of the photo-z's is evaluated using the subsample of 220 sources with secure spectroscopic redshifts. We achieve an accuracy of ~1% [|Delta z|/(1+z)] and an outlier fraction of ~1.4% for sources with spectroscopic redshifts. We performed blind tests to derive a more realistic estimate of the photo-z quality for sources without spectroscopic redshifts. We expect there are ~9% outliers for the relatively brighter sources (R<26), and the outlier fraction will increase to ~15-25% for the fainter sources (R>26). The typical photo-z accuracy is ~6-7%. The outlier fraction and photo-z accuracy do not appear to have a redshift dependence (for z~0-4). These photo-z's appear to be the best obtained so far for faint X-ray sources, and they have been significantly (>50%) improved compared to previous estimates of the photo-z's for the X-ray sources in the ~2 Ms Chandra Deep Field-North and ~1 Ms Chandra Deep Field-South.

We study the co-existence of star formation and AGN activity in X-ray selected AGN by analyzing stacked 870um submm emission from a deep and wide map of the Extended Chandra Deep Field South, obtained with LABOCA at the APEX telescope. The total X-ray sample of 895 sources with median redshift z~1 is detected at a mean submm flux of 0.49+-0.04mJy, corresponding to a typical star formation rate around 30Msun/yr for a T=35K, beta=1.5 greybody far-infrared SED. The good S/N permits stacking analyses for subgroups. We observe a trend of star formation rate increasing with redshift. An increase of star formation rate with AGN luminosity is indicated at the highest L_2-10>~1E44erg/s luminosities only. Increasing trends with X-ray obscuration as expected in some AGN evolutionary scenarios are not observed for the bulk of the X-ray AGN sample but may be present for the highest intrinsic luminosity objects. This suggests a transition between two modes in the coexistence of AGN activity and star formation. For the bulk of the sample, the X-ray luminosity and obscuration of the AGN are not intimately linked to the global star formation rate of their hosts. The hosts are likely massive and forming stars secularly, at rates similar to the pervasive star formation seen in massive galaxies without an AGN at similar redshifts. The change indicated towards more intense star formation, and a more pronounced increase in star formation rates between unobscured and obscured AGN at highest luminosities suggests that luminous AGN follow an evolutionary path on which obscured AGN activity and intense star formation are linked, possibly via merging. Comparison to local hard X-ray selected AGN supports this interpretation. [Abridged]

We investigate the correlation between far-infrared (FIR) and radio luminosities in distant galaxies, a lynchpin of modern astronomy. We use data from the Balloon-borne Large Aperture Submillimetre Telescope (BLAST), Spitzer, the Large Apex BOlometer CamerA (LABOCA), the Very Large Array (VLA) and the Giant Metre-wave Radio Telescope (GMRT) in the Extended Chandra Deep Field South (ECDFS). For a catalogue of BLAST 250-micron-selected galaxies, we re-measure the 70--870-micron flux densities at the positions of their most likely 24-micron counterparts, which have a median [interquartile] redshift of 0.74 [0.25, 1.57]. From these, we determine the monochromatic flux density ratio, q_250 = log_10 (S_250micron / S_1400MHz), and the bolometric equivalent, q_IR. At z ~= 0.6, where our 250-micron filter probes rest-frame 160-micron emission, we find no evolution relative to q_160 for local galaxies. We also stack the FIR and submm images at the positions of 24-micron- and radio-selected galaxies. The difference between q_IR seen for 250-micron- and radio-selected galaxies suggests star formation provides most of the IR luminosity in ~< 100-uJy radio galaxies, but rather less for those in the mJy regime. For the 24-micron sample, the radio spectral index is constant across 0 < z < 3, but q_IR exhibits tentative evidence of a steady decline such that q_IR is proportional to (1+z)^(-0.15 +/- 0.03) - significant evolution, spanning the epoch of galaxy formation, with major implications for techniques that rely on the FIR/radio correlation. We compare with model predictions and speculate that we may be seeing the increase in radio activity that gives rise to the radio background.

Chandra X-ray Observatory has revealed X-ray cavities in many nearby cooling flow clusters. The cavities trace feedback from the central active galactic nulceus (AGN) on the intracluster medium (ICM), an important ingredient in stabilizing cooling flows and in the process of galaxy formation and evolution. But, the prevalence and duty cycle of such AGN outbursts is not well understood. To this end, we study how the cooling is balanced by the cavity heating for a complete sample of clusters (the Brightest 55 clusters of galaxies, hereafter B55). In the B55, we found 33 cooling flow clusters, 20 of which have detected X-ray bubbles in their ICM. Among the remaining 13, all except Ophiuchus could have significant cavity power yet remain undetected in existing images. This implies that the duty cycle of AGN outbursts with significant heating potential in cooling flow clusters is at least 60 % and could approach 100 %, but deeper data is required to constrain this further.

Powering the 10^62 erg nuclear outburst in the MS0735.6+7421 cluster central galaxy by accretion implies that its supermassive black hole (SMBH) grew by ~6x10^8 solar masses over the past 100 Myr. We place upper limits on the amount of cold gas and star formation near the nucleus of <10^9 solar masses and <2 solar masses per year, respectively. These limits imply that an implausibly large fraction of the preexisting cold gas in the bulge must have been consumed by its SMBH at the rate of ~3-5 solar masses per year while leaving no trace of star formation. Such a high accretion rate would be difficult to maintain by stellar accretion or the Bondi mechanism, unless the black hole mass approaches 10^11 solar masses. Its feeble nuclear luminosities in the UV, I, and X-ray bands compared to its enormous mechanical power are inconsistent with rapid accretion onto a ~5x10^9 solar mass black hole. We suggest instead that the AGN outburst is powered by a rapidly-spinning black hole. A maximally-spinning, 10^9 solar mass black hole contains enough rotational energy, ~10^62 erg, to quench a cooling flow over its lifetime and to contribute significantly to the excess entropy found in the hot atmospheres of groups and clusters. Two modes of AGN feedback may be quenching star formation in elliptical galaxies centered in cooling halos at late times. An accretion mode that operates in gas-rich systems, and a spin mode operating at modest accretion rates. The spin conjecture may be avoided in MS0735 by appealing to Bondi accretion onto a central black hole whose mass greatly exceeds 10^10 solar mass. The host galaxy's unusually large, 3.8 kpc stellar core radius (light deficit) may witness the presence of an ultramassive black hole.

The brightest cluster galaxy (BCG) in the Abell 1664 cluster is unusually blue and is forming stars at a rate of ~ 23 M_{\sun} yr^{-1}. The BCG is located within 5 kpc of the X-ray peak, where the cooling time of 3.5x10^8 yr and entropy of 10.4 keV cm^2 are consistent with other star-forming BCGs in cooling flow clusters. The center of A1664 has an elongated, "bar-like" X-ray structure whose mass is comparable to the mass of molecular hydrogen, ~ 10^{10} M_{\sun} in the BCG. We show that this gas is unlikely to have been stripped from interloping galaxies. The cooling rate in this region is roughly consistent with the star formation rate, suggesting that the hot gas is condensing onto the BCG. We use the scaling relations of Birzan et al. 2008 to show that the AGN is underpowered compared to the central X-ray cooling luminosity by roughly a factor of three. We suggest that A1664 is experiencing rapid cooling and star formation during a low-state of an AGN feedback cycle that regulates the rates of cooling and star formation. Modeling the emission as a single temperature plasma, we find that the metallicity peaks 100 kpc from the X-ray center, resulting in a central metallicity dip. However, a multi-temperature cooling flow model improves the fit to the X-ray emission and is able to recover the expected, centrally-peaked metallicity profile.

We present results from an extensive survey of 64 cavities in the X-ray halos of clusters, groups and normal elliptical galaxies. We show that the evolution of the size of the cavities as they rise in the X-ray atmosphere is inconsistent with the standard model of adiabatic expansion of purely hydrodynamic models. We also note that the majority of the observed bubbles should have already been shredded apart by Rayleigh-Taylor and Richtmyer-Meshkov instabilities if they were of purely hydrodynamic nature. Instead we find that the data agrees much better with a model where the cavities are magnetically dominated and inflated by a current-dominated magneto-hydrodynamic jet model, recently developed by Li et al. (2006) and Nakamura et al. (2006). We conduct complex Monte-Carlo simulations of the cavity detection process including incompleteness effects to reproduce the cavity sample's characteristics. We find that the current-dominated model agrees within 1sigma, whereas the other models can be excluded at >5sigma confidence. To bring hydrodynamic models into better agreement, cavities would have to be continuously inflated. However, these assessments are dependent on our correct understanding of the detectability of cavities in X-ray atmospheres, and will await confirmation when automated cavity detection tools become available in the future. Our results have considerable impact on the energy budget associated with active galactic nucleus feedback.

Using broadband optical imaging and Chandra X-ray data for a sample of 46 cluster central dominant galaxies (CDGs), we investigate the connection between star formation, the intracluster medium (ICM), and the central active galactic nucleus (AGN). We report the discovery of a remarkably sharp threshold for the onset of star formation that occurs when the central cooling time of the hot atmosphere falls below ~ 5x10^8 yr, or equivalently when the central entropy falls below ~ 30 keV cm^2. In addition to this criterion, star formation in cooling flows also appears to require that the X-ray and galaxy centroids lie within ~ 20 kpc of each other, and that the jet (cavity) power is smaller than the X-ray cooling luminosity. These three criteria, together with the high ratio of cooling time to AGN outburst (cavity) age across our sample, directly link the presence of star formation and AGN activity in CDGs to cooling instabilities in the intracluster plasma. Our results provide compelling evidence that AGN feedback into the hot ICM is largely responsible for regulating cooling and star formation in the cores of clusters, leading to the significant growth of supermassive black holes in CDGs at late times.

Current models of galaxy evolution suggest that feedback from active galactic nuclei is needed to explain the high-luminosity cutoff in the galaxy luminosity function. Exactly how an AGN outflow couples with the ambient medium and suppresses star formation remains poorly understood. However, we have recently uncovered an important clue to how that coupling might work. Observations of H-alpha emission and blue light from the universe's most luminous galaxies, which occupy the centers of galaxy clusters, show that star formation happens only if the minimum specific entropy of the intracluster gas is less than 30 keV cm^2. Here we suggest that this threshold for star formation is set by the physics of electron thermal conduction, implying that conduction is critical for channeling AGN energy input toward incipient star-forming regions and limiting the progress of star formation.

We present results from magneto-hydrodynamic jet simulations in the magnetically dominated regime and find that they inflate cavities into the intracluster medium. Mock Chandra X-ray observations of our simulations show a strong resemblance to real cavities observed in X-ray images of clusters, both in terms of their morphology and thermodynamic structure. An analysis of the evolution of bubble sizes in the multi-cavity system Hydra A, as well as in a large sample of 64 cavities in 32 clusters shows that bubbles tend to expand much faster than expected in the purely adiabatic regime. Instead, we find that the bubbles follow more closely a trend predicted by our current-carrying jet models.

X-ray observations of cavities and shock fronts produced by jets streaming through hot halos have significantly advanced our understanding of the energetics and dynamics of extragalactic radio sources. Radio sources at the centers of clusters have dynamical ages between ten and several hundred million years. They liberate between 1E58-1E62 erg per outburst, which is enough energy to regulate cooling of hot halos from galaxies to the richest clusters. Jet power scales approximately with the radio synchrotron luminosity to the one half power. However, the synchrotron efficiency varies widely from nearly unity to one part in 10,000, such that relatively feeble radio source can have quasar-like mechanical power. The synchrotron ages of cluster radio sources are decoupled from their dynamical ages, which tend to be factors of several to orders of magnitude older. Magnetic fields and particles in the lobes tend to be out of equipartition. The lobes may be maintained by heavy particles (e.g., protons), low energy electrons, a hot, diffuse thermal gas, or possibly magnetic (Poynting) stresses. Sensitive X-ray images of shock fronts and cavities can be used to study the dynamics of extragalactic radio sources.

Central cluster galaxies (cDs) in cooling flows are growing rapidly through gas accretion and star formation. At the same time, AGN outbursts fueled by accretion onto supermassive black holes are generating X-ray cavity systems and driving outflows that exceed those in powerful quasars. We show that the resulting bulge and black hole growth follows a trend that is roughly consistent with the slope of the local (Magorrian) relation between bulge and black hole mass for nearby quiescent ellipticals. However, a large scatter suggests that cD bulges and black holes do not always grow in lock-step. New measurements made with XMM, Chandra, and FUSE of the condensation rates in cooling flows are now approaching or are comparable to the star formation rates, alleviating the need for an invisible sink of cold matter. We show that the remaining radiation losses can be offset by AGN outbursts in more than half of the systems in our sample, indicating that the level of cooling and star formation is regulated by AGN feedback.

Three comments are made on AGN heating of cooling flows. A simple physical argument is used to show that the enthalpy of a buoyant radio lobe is converted to heat in its wake. Thus, a significant part of ``cavity'' enthalpy is likely to end up as heat. Second, the properties of the repeated weak shocks in M87 are used to argue that they can plausibly prevent gas close to the AGN from cooling. As the most significant heating mechanism at work closest to the AGN, shock heating probably plays a critical role in the feedback mechanism. Third, results are presented from a survey of AGN heating rates in nearby giant elliptical galaxies. With inactive systems included, the overall AGN heating rate is reasonably well matched to the total cooling rate for the sample. Thus, intermittent AGN outbursts are energetically capable of preventing the hot atmospheres of these galaxies from cooling and forming stars.

We present an analysis of the growth of black holes through accretion and bulges through star formation in 33 galaxies at the centers of cooling flows. Most of these systems show evidence of cavities in the intracluster medium (ICM) inflated by radio jets emanating from their active galactic nuclei (AGN). We present a new and extensive analysis of X-ray cavities in these systems. We find that AGN are energetically able to balance radiative losses (cooling) from the ICM in more than half of our sample. Using a subsample of 17 systems, we examine the relationship between cooling and star formation. We find that the star formation rates are approaching or are comparable to X-ray and far UV limits on the rates of gas condensation onto the central galaxy. The remaining radiative losses could be offset by AGN feedback. The vast gulf between radiative losses and the sink of cooling material, which has been the primary objection to cooling flows, has narrowed and, in some cases, is no longer a serious issue. Using the cavity (jet) powers, we place strong lower limits on the rate of growth of supermassive black holes in central galaxies, and we find that they are growing at an average rate of ~ 0.1 solar masses per year, with some systems growing as quickly as ~ 1 solar mass per year. We find a trend between bulge growth (star formation) and black hole growth that is approximately in accordance with the slope of the local (Magorrian) relation between black hole and bulge mass. However, the large scatter in the trend suggests that bulges and black holes do not always grow in lock step. (Abridged)

We present an optical, X-ray, and radio analysis of the starburst in the Abell 1835 cluster's central cD galaxy. The dense gas surrounding the galaxy is radiating X-rays with a luminosity of ~1E45 erg/s consistent with a cooling rate of ~1000-2000 solar masses per year. However, Chandra and XMM-Newton observations find less than 200 solar masses per year of gas cooling below ~2 keV, a level that is consistent with the cD's current star formation rate of 100-180 solar masses per year. One or more heating agents (feedback) must then be replenishing the remaining radiative losses. The heat fluxes from supernova explosions and thermal conduction alone are unable to do so. However, a pair of X-ray cavities from an AGN outburst has deposited ~1.7E60 erg into the surrounding gas over the past 40 Myr. The corresponding jet power 1.4E45 erg/sec is enough to offset most of the radiative losses from the cooling gas. The jet power exceeds the radio synchrotron power by ~4000 times, making this one of the most radiatively inefficient radio sources known. The large jet power implies that the cD's supermassive black hole accreted at a mean rate of ~0.3 solar masses per year over the last 40 Myr or so, which is a small fraction of the Eddington accretion rate for a 10E9 solar mass black hole. The ratio of the bulge growth rate through star formation and the black hole growth rate through accretion is consistent with the slope of the (Magorrian) relationship between bulge and central black hole mass in nearby quiescent galaxies. The consistency between net cooling, heating (feedback), and the cooling sink (star formation) in this system resolves the primary objection to traditional cooling flow models. (abridged)

The radio source Hercules A resides at the center of a cooling flow cluster of galaxies at redshift z = 0.154. A Chandra X-ray image reveals a shock front in the intracluster medium (ICM) surrounding the radio source, about 160 kpc from the active galactic nucleus (AGN) that hosts it. The shock has a Mach number of 1.65, making it the strongest of the cluster-scale shocks driven by an AGN outburst found so far. The age of the outburst ~5.9e7 y, its energy about 3e61 erg and its mean power ~1.6e46 erg/s. As for the other large AGN outbursts in cooling flow clusters, this outburst overwhelms radiative losses from the ICM of the Hercules A cluster by a factor of ~100. It adds to the case that AGN outbursts are a significant source of preheating for the ICM. Unless the mechanical efficiency of the AGN in Hercules A exceeds 10%, the central black hole must have grown by more than 1.7e8 Msun to power this one outburst.

We report the discovery of a pair of giant cavities, each nearly 200 kpc in diameter, embedded in the gaseous halo of the redshift z=0.216 galaxy cluster MS0735.6+7421. The cavities appear as surface brightness depressions in an image of the cluster taken with the Chandra X-ray Observatory. The X-ray cavities are filled with radio emission and are bounded by an elliptical shock structure resembling a ``radio cocoon.'' The shock's energy and age, ~6E61 erg and ~100 Myr, respectively, give a time averaged power of ~1.7E46 erg/s, making it the most powerful radio outburst known. The outburst is apparently heating the cluster by 1/3 keV per particle, which is sufficient to quench a large cooling flow for several Gyr and to supply a substantial fraction of the energy required to "preheat" the cluster. If the outburst was powered by accretion onto a supermassive black hole, more than 3E8 solar masses of material would be required to do so. This implies that the putative supermassive black hole lying at the center of the cD galaxy has accreted a substantial fraction of its own mass in less than 100 Myr.

We present an analysis of sixteen galaxy clusters, one group and one galaxy drawn from the Chandra X-ray Observatory's data archive. These systems possess prominent X-ray surface brightness depressions associated with cavities or bubbles that were created by interactions between powerful radio sources and the surrounding hot gas. The minimum energy associated with the cavities ranges between pV~10^55 erg in galaxies, groups, and poor clusters to pV~10^60 erg in rich clusters. We evaluate the hypothesis that cooling in the hot gas can be quenched by energy injected into the surrounding gas by the rising bubbles. Nearly half of the systems in this study may have instantaneous mechanical luminosities large enough to balance cooling, at least for a short period of time, if the cavities are filled with a relativistic gas. We find a trend or upper envelope in the distribution of central X-ray luminosity versus instantaneous mechanical luminosity with the sense that the most powerful cavities are found in the most X-ray luminous systems. Such a trend would be expected if many of these systems produce bubbles at a rate that scales in proportion to the cooling rate of the surrounding gas. Finally, we use the X-ray cavities to measure the mechanical power of radio sources over six decades of radio luminosity, independently of the radio properties themselves. We find that the ratio of the instantaneous mechanical (kinetic) luminosity to the 1.4 GHz synchrotron luminosity ranges from a few to roughly a thousand. This wide range implies that the 1.4 GHz synchrotron luminosity is an unreliable gauge of the mechanical power of radio sources.