We report on the detection of extreme giant pulses (GPs) from one of the oldest-known pulsars, the highly variable PSR B0950+08, with the Amsterdam-ASTRON Radio Transient Facility And Analysis Centre (AARTFAAC), a parallel transient detection instrument operating as a subsystem of the LOw Frequency ARray (LOFAR). During processing of our Northern Hemisphere survey for low-frequency radio transients, a sample of 275 pulses with fluences ranging from 42 to 177 kJy ms were detected in one-second snapshot images. The brightest pulses are an order of magnitude brighter than those previously reported at 42 and 74 MHz, on par with the levels observed in a previous long-term study at 103 MHz. Both their rate and fluence distribution differ between and within the various studies done to date. The GP rate is highly variable, from 0 to 30 per hour, with only two three-hour observations accounting for nearly half of the pulses detected in the 96 h surveyed. It does not vary significantly within a few-hour observation, but can vary strongly one from day to the next. The spectra appear strongly and variably structured, with emission sometimes confined to a single 195.3 kHz subband, and the pulse spectra changing on a timescale of order 10 min.

Restarted radio galaxies represent a unique tool to investigate the duty cycle of the jet activity in Active Galactic Nuclei. Due to a sharp discontinuity of the GHz spectral index distribution within its lobes, the radio galaxy 3C388 has for long being claimed to be a peculiar example of AGN with multi-epoch activity. In this work we have performed, for the first time, a spatially resolved study of the radio spectrum of this source down to MHz-frequencies, aimed at investigating its radiative age and duty cycle. We have used new low frequency observations at 144 MHz performed with LOFAR and at 350 MHz performed with the VLA combined with archival data at higher frequencies. We find that the spectral indices in the lower frequency range 144-614 MHz have flatter values ($\sim$0.55-1.14) than those observed in the higher frequency range 1400-4850 MHz ($\sim$0.75-1.57) but follow the same distribution across the lobes, with a systematic steepening towards the edges. However, the spectral shape throughout the source is not uniform and often deviates from standard models. This suggests that mixing of different particle populations is occurring, although it remains difficult to understand whether this is caused by observational limitations (insufficient spatial resolution and/or projection effects) or by the intrinsic presence of multiple particle populations, possibly related to the two different outbursts. By using single-injection radiative models we compute that the total source age is $\lesssim$80 Myr and that the duty cycle is about $\rm t_{on}/t_{tot}\sim$ 60%, which is enough to prevent the intracluster medium from cooling according to X-ray estimates. While to date the radio spectral distribution of 3C388 remains a rare case among radio galaxies, upcoming multi-frequency surveys will soon allow us to investigate whether more sources with the same characteristics do actually exist.

We present LOFAR observations at 150 MHz of the borderline FRI/FRII giant radio galaxy NGC 6251. This paper presents the most sensitive and highest-resolution images of NGC 6251 at these frequencies to date, revealing for the first time a low-surface-brightness extension to the northern lobe, and a possible backflow associated with the southern lobe. The integrated spectra of components of NGC 6251 are consistent with previous measurements at higher frequencies, similar to results from other LOFAR studies of nearby radio galaxies. We find the outer structures of NGC 6251 to be either at equipartition or slightly electron dominated, similar to those of FRII sources rather than FRIs; but this conclusion remains tentative because of uncertainties associated with the geometry and the extrapolation of X-ray measurements to determine the external pressure distribution on the scale of the outer lobes. We place lower limits on the ages of the extension of the northern lobe and the backflow of the southern lobe of $t \gtrsim 250$ Myr and $t \gtrsim 210$ Myr respectively. We present the first detection of polarisation at 150 MHz in NGC 6251. Taking advantage of the high Faraday resolution of LOFAR, we place an upper limit on the magnetic field in the group of $B < 0.2 (\Lambda_B / 10 {\rm kpc})^{-0.5} \mu$G for a coherence scale of $\Lambda_B < 60 {\rm kpc}$ and $B < 13 \mu$G for $\Lambda_B = 240$ kpc.

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

We report the experimental setup and overall results of the AARTFAAC wide-field radio survey, which consists of observing the sky within 50 deg of Zenith with a bandwidth of 3.2 MHz at a cadence of 1 s, for 545 h. This yielded nearly 4 million snapshots, two per second, of on average 4800 square degrees and a sensitivity of around 60 Jy. We find two populations of transient events, one originating from PSR B0950+08 and one from strong ionospheric lensing events, as well as a possible single extragalactic transient. We also set a strong upper limit of 1.1 all-sky per day to the rate of any other populations of fast, bright transients. Lastly, we constrain some previously detected types of transient source by comparing our detections and limits with other low-frequency radio transient surveys.

We report the discovery of bright, fast, radio flares lasting tens of seconds with the AARTFAAC high-cadence all-sky survey at 60 MHz. The vast majority of these coincide with known, bright radio sources that brighten by factors of up to 100 during such an event. We attribute them to magnification events induced by plasma near the Earth, most likely in the densest parts of the ionosphere. They can occur both in relative isolation, during otherwise quiescent ionospheric conditions, and in large clusters during more turbulent ionospheric conditions. We also find one event that has many properties in common with these events, but has a very well-determined dispersion delay as a function of frequency consistent with a very distant Galactic or extragalactic origin. If astrophysical, it represents extreme emission with a minimum brightness temperature of $10^{28-34}$ K.

This paper presents the results from one of the first observations of ionospheric scintillation taken using the Low-Frequency Array (LOFAR). The observation was of the strong natural radio source Cas A, taken overnight on 18-19 August 2013, and exhibited moderately strong scattering effects in dynamic spectra of intensity received across an observing bandwidth of 10-80MHz. Delay-Doppler spectra (the 2-D FFT of the dynamic spectrum) from the first hour of observation showed two discrete parabolic arcs, one with a steep curvature and the other shallow, which can be used to provide estimates of the distance to, and velocity of, the scattering plasma. A cross-correlation analysis of data received by the dense array of stations in the LOFAR "core" reveals two different velocities in the scintillation pattern: a primary velocity of ~30m/s with a north-west to south-east direction, associated with the steep parabolic arc and a scattering altitude in the F-region or higher, and a secondary velocity of ~110m/s with a north-east to south-west direction, associated with the shallow arc and a scattering altitude in the D-region. Geomagnetic activity was low in the mid-latitudes at the time, but a weak sub-storm at high latitudes reached its peak at the start of the observation. An analysis of Global Navigation Satellite Systems (GNSS) and ionosonde data from the time reveals a larger-scale travelling ionospheric disturbance (TID), possibly the result of the high-latitude activity, travelling in the north-west to south-east direction, and, simultaneously, a smaller--scale TID travelling in a north-east to south-west direction, which could be associated with atmospheric gravity wave activity. The LOFAR observation shows scattering from both TIDs, at different altitudes and propagating in different directions. To the best of our knowledge this is the first time that such a phenomenon has been reported.

The four persistent radio sources in the northern sky with the highest flux density at metre wavelengths are Cassiopeia A, Cygnus A, Taurus A, and Virgo A; collectively they are called the A-team. Their flux densities at ultra-low frequencies (<100 MHz) can reach several thousands of janskys, and they often contaminate observations of the low-frequency sky by interfering with image processing. Furthermore, these sources are foreground objects for all-sky observations hampering the study of faint signals, such as the cosmological 21 cm line from the epoch of reionisation. We aim to produce robust models for the surface brightness emission as a function of frequency for the A-team sources at ultra-low frequencies. These models are needed for the calibration and imaging of wide-area surveys of the sky with low-frequency interferometers. This requires obtaining images at an angular resolution better than 15 arcsec with a high dynamic range and good image fidelity. We observed the A-team with the Low Frequency Array (LOFAR) at frequencies between 30 MHz and 77 MHz using the Low Band Antenna (LBA) system. We reduced the datasets and obtained an image for each A-team source. The paper presents the best models to date for the sources Cassiopeia A, Cygnus A, Taurus A, and Virgo A between 30 MHz and 77 MHz. We were able to obtain the aimed resolution and dynamic range in all cases. Owing to its compactness and complexity, observations with the long baselines of the International LOFAR Telescope will be required to improve the source model for Cygnus A further.

Giant radio galaxies (GRGs) are a subclass of radio galaxies which have grown to megaparsec scales. GRGs are much rarer than normal sized radio galaxies (< 0.7 Mpc) and the reason for their gigantic sizes is still debated. Here, we report the biggest sample of GRGs identified to date. These objects were found in the LOFAR Two-metre Sky Survey (LoTSS) first data release images, which cover a 424 square degrees region. Of the 239 GRGs found, 225 are new discoveries. The GRGs in our sample have sizes ranging from 0.7 to 3.5 Mpc and have redshifts (z) between 0.1 and 2.3. Seven GRGs have sizes above 2 Mpc and one has a size of ~ 3.5 Mpc. The sample contains 40 GRGs hosted by spectroscopically confirmed quasars. Here, we present the search techniques employed and the resulting catalogue of the newly discovered large sample of GRGs along with their radio properties. We, here also show for the first time that the spectral index of GRGs is similar to that of normal sized radio galaxies, indicating that most of the GRG population is not dead or is not like remnant type radio galaxy. We find 20/239 GRGs in our sample are located at the centres of clusters and present our analysis on their cluster environment and radio morphology.

The nature of the central engines of gamma-ray bursts (GRBs) and the composition of their relativistic jets are still under debate. If the jets are Poynting flux dominated rather than baryon dominated, a coherent radio flare from magnetic re-connection events might be expected with the prompt gamma-ray emission. There are two competing models for the central engines of GRBs; a black hole or a newly formed milli-second magnetar. If the central engine is a magnetar it is predicted to produce coherent radio emission as persistent or flaring activity. In this paper, we present the deepest limits to date for this emission following LOFAR rapid response observations of GRB 180706A. No emission is detected to a 3$\sigma$ limit of 1.7 mJy beam$^{-1}$ at 144 MHz in a two-hour LOFAR observation starting 4.5 minutes after the gamma-ray trigger. A forced source extraction at the position of GRB 180706A provides a marginally positive (1 sigma) peak flux density of $1.1 \pm 0.9$ mJy. The data were time-sliced into different sets of snapshot durations to search for FRB like emission. No short duration emission was detected at the location of the GRB. We compare these results to theoretical models and discuss the implications of a non-detection.

Low-luminosity radio-loud active galactic nuclei (AGN) are of importance in studies concerning feedback from radio AGN since a dominant fraction of AGN belong to this class. We report high-resolution Very Large Array (VLA) and European VLBI Network (EVN) observations of HI-21cm absorption from a young, compact steep-spectrum radio source, B2 0258+35, nested in the early-type galaxy NGC 1167, which contains a 160 kpc HI disc. Our VLA and EVN HI absorption observations, modelling, and comparison with molecular gas data suggest that the cold gas in the centre of NGC 1167 is very turbulent (with a velocity dispersion of ~ 90 km/s) and that this turbulence is induced by the interaction of the jets with the interstellar medium (ISM). Furthermore, the ionised gas in the galaxy shows evidence of shock heating at a few kpc from the radio source. These findings support the results from numerical simulations of radio jets expanding into a clumpy gas disc, which predict that the radio jets in this case percolate through the gas disc and drive shocks into the ISM at distances much larger than their physical extent. These results expand the number of low-luminosity radio sources found to impact the surrounding medium, thereby highlighting the possible relevance of these AGN for 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.

We present an overview of the LOFAR Tied-Array All-Sky Survey (LOTAAS) for radio pulsars and fast transients. The survey uses the high-band antennas of the LOFAR Superterp, the dense inner part of the LOFAR core, to survey the northern sky (dec > 0 deg) at a central observing frequency of 135 MHz. A total of 219 tied-array beams (coherent summation of station signals, covering 12 square degrees), as well as three incoherent beams (covering 67 square degrees) are formed in each survey pointing. For each ofthe 222 beams, total intensity is recorded at 491.52 us time resolution. Each observation integrates for 1 hr and covers 2592 channels from 119 to 151 MHz. This instrumental setup allows LOTAAS to reach a detection threshold of 1 to 5 mJy for periodic emission. Thus far, the LOTAAS survey has resulted in the discovery of 73 radio pulsars. Among these are two mildly recycled binary millisecond pulsars (P = 13 and 33 ms), as well as the slowest-spinning radio pulsar currently known (P = 23.5 s). The survey has thus far detected 311 known pulsars, with spin periods ranging from 4 ms to 5.0 s and dispersion measures from 3.0 to 217 pc/cc. Known pulsars are detected at flux densities consistent with literature values. We find that the LOTAAS pulsar discoveries have, on average, longer spin periods than the known pulsar population. This may reflect different selection biases between LOTAAS and previous surveys, though it is also possible that slower-spinning pulsars preferentially have steeper radio spectra. LOTAAS is the deepest all-sky pulsar survey using a digital aperture array; we discuss some of the lessons learned that can inform the approach for similar surveys using future radio telescopes such as the Square Kilometre Array.

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.

Double-double radio galaxies (DDRGs) represent a short but unique phase in the life-cycle of some of the most powerful radio-loud active galactic nuclei (RLAGN). These galaxies display large-scale remnant radio plasma in the intergalactic medium left behind by a past episode of active galactic nuclei (AGN) activity, and meanwhile, the radio jets have restarted in a new episode. The knowledge of what causes the jets to switch off and restart is crucial to our understanding of galaxy evolution, while it is important to know if DDRGs form a host galaxy dichotomy relative to RLAGN. We utilised the LOFAR Two-Metre Sky Survey DR1, using a visual identification method to compile a sample of morphologically selected candidate DDRGs, showing two pairs of radio lobes. To confirm the restarted nature in each of the candidate sources, we obtained follow-up observations with the VLA at higher resolution to observe the inner lobes or restarted jets, the confirmation of which created a robust sample of 33 DDRGs. We created a comparison sample of 777 RLAGN from the DR1 catalogue, and compared the optical and infrared magnitudes and colours of their host galaxies. We find that there is no statistically significant difference in the brightness of the host galaxies between double-doubles and single-cycle RLAGN. The DDRG and RLAGN samples also have similar distributions in WISE mid-infrared colours, indicating similar ages of stellar populations and dust levels in the hosts of DDRGs. We conclude that DDRGs and 'normal' RLAGN are hosted by galaxies of the same type, and that DDRG activity is simply a normal part of the life cycle of RLAGN. Restarted jets, particularly for the class of low-excitation radio galaxies, rather than being a product of a particular event in the life of a host galaxy, must instead be caused by smaller scale changes, such as in the accretion system surrounding the black hole.

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.

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.

Historically, the blazar population has been poorly understood at low frequencies because survey sensitivity and angular resolution limitations have made it difficult to identify megahertz counterparts. We used the LOFAR Two-Metre Sky Survey (LoTSS) first data release value-added catalogue (LDR1) to study blazars in the low-frequency regime with unprecedented sensitivity and resolution. We identified radio counterparts to all $98$ known sources from the Third \textit{Fermi}-LAT Point Source Catalogue (3FGL) or Roma-BZCAT Multi-frequency Catalogue of Blazars ($5^{\mathrm{th}}$ edition) that fall within the LDR1 footprint. Only the 3FGL unidentified $\gamma$-ray sources (UGS) could not be firmly associated with an LDR1 source; this was due to source confusion. We examined the redshift and radio luminosity distributions of our sample, finding flat-spectrum radio quasars (FSRQs) to be more distant and more luminous than BL Lacertae objects (BL Lacs) on average. Blazars are known to have flat spectra in the gigahertz regime but we found this to extend down to $144$ MHz, where the radio spectral index, $\alpha$, of our sample is $-0.17 \pm 0.14$. For BL Lacs, $\alpha = -0.13 \pm 0.16$ and for FSRQs, $\alpha = -0.15 \pm 0.17$. We also investigated the radio-to-$\gamma$-ray connection for the $30$ $\gamma$-ray-detected sources in our sample. We find Pearson's correlation coefficient is $0.45$ ($p = 0.069$). This tentative correlation and the flatness of the spectral index suggest that the beamed core emission contributes to the low-frequency flux density. We compare our sample distribution with that of the full LDR1 on colour-colour diagrams, and we use this information to identify possible radio counterparts to two of the four UGS within the LDR1 field. We will refine our results as LoTSS continues.

The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168 MHz survey of the Northern sky with diverse and ambitious science goals. Many of the scientific objectives of LoTSS rely upon, or are enhanced by, the association or separation of the sometimes incorrectly catalogued radio components into distinct radio sources, and the identification and characterisation of the optical counterparts to these sources. Here we present the source associations and optical and/or IR identifications for sources in the first data release, which are made using a combination of statistical techniques and visual association and identification. We document in detail the colour- and magnitude-dependent likelihood ratio method used for statistical identification as well as the Zooniverse project, called LOFAR Galaxy Zoo, used for the visual classification. We describe the process used to select which of these two different methods is most appropriate for each LoTSS source. The final LoTSS-DR1-IDs value-added catalogue presented contains 318,520 radio sources, of which 231,716 (73%) have optical and/or IR identifications in Pan-STARRS and WISE. The value-added catalogue is available online at https://lofar-surveys.org/, as part of this data release.

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.

We constructed a sample of 23,344 radio-loud active galactic nuclei (RLAGN) from the catalogue derived from the LOFAR Two-Metre Sky Survey (LoTSS) survey of the HETDEX Spring field. Although separating AGN from star-forming galaxies remains challenging, the combination of spectroscopic and photometric techniques we used gives us one of the largest available samples of candidate RLAGN. We used the sample, combined with recently developed analytical models, to investigate the lifetime distribution of RLAGN. We show that large or giant powerful RLAGN are probably the old tail of the general RLAGN population, but that the low-luminosity RLAGN candidates in our sample, many of which have sizes $<100$ kpc, either require a very different lifetime distribution or have different jet physics from the more powerful objects. We then used analytical models to develop a method of estimating jet kinetic powers for our candidate objects and constructed a jet kinetic luminosity function based on these estimates. These values can be compared to observational quantities, such as the integrated radiative luminosity of groups and clusters, and to the predictions from models of RLAGN feedback in galaxy formation and evolution. In particular, we show that RLAGN in the local Universe are able to supply all the energy required per comoving unit volume to counterbalance X-ray radiative losses from groups and clusters and thus prevent the hot gas from cooling. Our computation of the kinetic luminosity density of local RLAGN is in good agreement with other recent observational estimates and with models of galaxy formation.

The shape of low-frequency radio continuum spectra of normal galaxies is not well understood, the key question being the role of physical processes such as thermal absorption in shaping them. In this work we take advantage of the LOFAR Multifrequency Snapshot Sky Survey (MSSS) to investigate such spectra for a large sample of nearby star-forming galaxies. Using the measured 150MHz flux densities from the LOFAR MSSS survey and literature flux densities at various frequencies we have obtained integrated radio spectra for 106 galaxies. The spectra are explained through the use of a three-dimensional model of galaxy radio emission, and radiation transfer dependent on the galaxy viewing angle and absorption processes. Spectra of our galaxies are generally flatter at lower compared to higher frequencies but as there is no tendency for the highly inclined galaxies to have more flattened low-frequency spectra, we argue that the observed flattening is not due to thermal absorption, contradicting the suggestion of Israel & Mahoney (1990). According to our modelled radio maps for M51-like galaxies, the free-free absorption effects can be seen only below 30MHz and in the global spectra just below 20MHz, while in the spectra of starburst galaxies, like M82, the flattening due to absorption is instead visible up to higher frequencies of about 150MHz. Locally, within galactic disks, the absorption effects are distinctly visible in M51-like galaxies as spectral flattening around 100-200MHz in the face-on objects, and as turnovers in the edge-on ones, while in M82-like galaxies there are strong turnovers at frequencies above 700MHz, regardless of viewing angle. Our modelling suggests that the weak spectral flattening observed in the nearby galaxies studied here results principally from synchrotron spectral curvature due to cosmic ray energy losses and propagation effects.

Radio loud Active Galactic Nuclei are episodic in nature, cycling through periods of activity and quiescence. In this work we investigate the duty cycle of the radio galaxy B2~0258+35, which was previously suggested to be a restarted radio galaxy based on its morphology. The radio source consists of a pair of kpc-scale jets embedded in two large-scale lobes (~240 kpc) with relaxed shape and very low surface brightness, which resemble remnants of a past AGN activity. We have combined new LOFAR data at 145 MHz and new SRT data at 6600 MHz with available WSRT data at 1400 MHz to investigate the spectral properties of the outer lobes and derive their age. Interestingly, the spectrum of both the outer Northern and Southern lobes is not ultra-steep as expected for an old ageing plasma with spectral index values equal to $\rm \alpha_{1400}^{145}=0.48\pm0.11$ and $\rm \alpha_{6600}^{1400}=0.69\pm0.20$ in the outer Northern lobe, and $\rm \alpha_{1400}^{145}=0.73\pm0.07$ in the outer Southern lobe. Moreover, despite the wide frequency coverage available for the outer Northern lobe (145-6600~MHz), we do not identify a significant spectral curvature (SPC$\simeq$0.2$\pm0.2$). While mechanisms such as in-situ particle reacceleration, mixing or compression can temporarily play a role in preventing the spectrum from steepening, in no case seem the outer lobes to be compatible with being very old remnants of past activity as previously suggested (with age $\gtrsim$ 80 Myr). We conclude that either the large-scale lobes are still fuelled by the nuclear engine or the jets have switched off no more than a few tens of Myr ago. Our study shows the importance of combining morphological and spectral properties to reliably classify the evolutionary stage of low surface brightness, diffuse emission that low frequency observations are revealing around a growing number of radio sources.

Type II radio bursts are evidence of shocks in the solar atmosphere and inner heliosphere that emit radio waves ranging from sub-meter to kilometer lengths. These shocks may be associated with CMEs and reach speeds higher than the local magnetosonic speed. Radio imaging of decameter wavelengths (20-90 MHz) is now possible with LOFAR, opening a new radio window in which to study coronal shocks that leave the inner solar corona and enter the interplanetary medium and to understand their association with CMEs. To this end, we study a coronal shock associated with a CME and type II radio burst to determine the locations at which the radio emission is generated, and we investigate the origin of the band-splitting phenomenon.

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.

We present a deep, low-frequency radio continuum study of the nearby Fanaroff--Riley class I (FR I) radio galaxy 3C 31 using a combination of LOw Frequency ARray (LOFAR; 30--85 and 115--178 MHz), Very Large Array (VLA; 290--420 MHz), Westerbork Synthesis Radio Telescope (WSRT; 609 MHz) and Giant Metre Radio Telescope (GMRT; 615 MHz) observations. Our new LOFAR 145-MHz map shows that 3C 31 has a largest physical size of $1.1$ Mpc in projection, which means 3C 31 now falls in the class of giant radio galaxies. We model the radio continuum intensities with advective cosmic-ray transport, evolving the cosmic-ray electron population and magnetic field strength in the tails as functions of distance to the nucleus. We find that if there is no in-situ particle acceleration in the tails, then decelerating flows are required that depend on radius $r$ as $v\propto r^{\beta}$ ($\beta\approx -1$). This then compensates for the strong adiabatic losses due to the lateral expansion of the tails. We are able to find self-consistent solutions in agreement with the entrainment model of Croston & Hardcastle, where the magnetic field provides $\approx$$1/3$ of the pressure needed for equilibrium with the surrounding intra-cluster medium (ICM). We obtain an advective time-scale of $\approx$$190$ Myr, which, if equated to the source age, would require an average expansion Mach number ${\cal M} \approx 5$ over the source lifetime. Dynamical arguments suggest that instead, either the outer tail material does not represent the oldest jet plasma or else the particle ages are underestimated due to the effects of particle acceleration on large scales.

The phase of radio galaxy evolution after the jets have switched off, often referred to as the remnant phase, is poorly understood and very few sources in this phase are known. In this work we present an extensive search for remnant radio galaxies in the Lockman Hole, a well-studied extragalactic field. We create mock catalogues of low-power radio galaxies based on Monte Carlo simulations to derive first-order predictions of the fraction of remnants in radio flux limited samples for comparison with our Lockman-Hole sample. We have combined LOFAR observations at 150 MHz with public surveys at higher frequencies to perform a complete selection and have used, for the first time, a combination of spectral criteria (e.g. the classical ultra-steep spectral index and high spectral curvature) as well as morphological criteria (e.g. low radio core prominence and relaxed shapes). Mock catalogues of radio galaxies are created based on existing spectral and dynamical evolution models combined with observed source properties. We have identified 23 candidate remnant radio galaxies which cover a variety of morphologies and spectral characteristics. We suggest that these different properties are related to different stages of the remnant evolution. We find that ultra-steep spectrum remnants represent only a fraction of our remnant sample suggesting a very rapid luminosity evolution of the radio plasma. Results from mock catalogues demonstrate the importance of dynamical evolution in the remnant phase of low-power radio galaxies to obtain fractions of remnant sources consistent with our observations. Moreover, these results confirm that ultra-steep spectrum remnants represent only a subset of the entire population ($\sim$50%) when frequencies higher than 1400 MHz are not included in the selection process, and that they are biased towards old ages.

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.

Low radio frequency surveys are important for testing unified models of radio-loud quasars and radio galaxies. Intrinsically similar sources that are randomly oriented on the sky will have different projected linear sizes. Measuring the projected linear sizes of these sources provides an indication of their orientation. Steep-spectrum isotropic radio emission allows for orientation-free sample selection at low radio frequencies. We use a new radio survey of the Bo\"otes field at 150 MHz made with the Low Frequency Array (LOFAR) to select a sample of radio sources. We identify 44 radio galaxies and 16 quasars with powers $P>10^{25.5}$ W Hz$^{-1}$ at 150 MHz using cross-matched multi-wavelength information from the AGN and Galaxy Evolution Survey (AGES), which provides spectroscopic redshifts. We find that LOFAR-detected radio sources with steep spectra have projected linear sizes that are on average 4.4$\pm$1.4 larger than those with flat spectra. The projected linear sizes of radio galaxies are on average 3.1$\pm$1.0 larger than those of quasars (2.0$\pm$0.3 after correcting for redshift evolution). Combining these results with three previous surveys, we find that the projected linear sizes of radio galaxies and quasars depend on redshift but not on power. The projected linear size ratio does not correlate with either parameter. The LOFAR data is consistent within the uncertainties with theoretical predictions of the correlation between the quasar fraction and linear size ratio, based on an orientation-based unification scheme.

In this paper, the second in a series investigating FR II radio galaxies at low frequencies, we use LOFAR and VLA observations between 117 and 456 MHz in addition to archival data to determine the dynamics and energetics of two radio galaxies, 3C452 and 3C223, through fitting of spectral ageing models on small spatial scales. We provide improved measurements for the physical extent of the two sources, including a previously unknown low surface brightness extension to the northern lobe of 3C223, and revised energetics based on these values. We find spectral ages of $77.05^{+9.22}_{-8.74}$ and $84.96^{+15.02}_{-13.83}$ Myr for 3C452 and 3C223 respectively suggesting a characteristic advance speed for the lobes of around one per cent the speed of light. For 3C452 we show that, even for a magnetic field strength not assumed to be in equipartition, a disparity of factor of approximately 2 exists between the spectral age and that determined from a dynamical standpoint. We confirm that the injection index of both sources (as derived from the lobe emission) remains steeper than classically assumed values even when considered on well resolved scales at low frequencies, but find an unexpected sharp discontinuity between the spectrum of the hotspots and the surrounding lobe emission. We suggest that this discrepancy is due to the absorption of hotspot emission and/or non-homogeneous and additional acceleration mechanisms and, as such, hotspots should not be used in the determination of the underlying initial electron energy distribution.

We report on the discovery in the LOFAR Multifrequency Snapshot Sky Survey (MSSS) of a giant radio galaxy (GRG) with a projected size of $2.56 \pm 0.07$ Mpc projected on the sky. It is associated with the galaxy triplet UGC 9555, within which one is identified as a broad-line galaxy in the Sloan Digital Sky Survey (SDSS) at a redshift of $0.05453 \pm 1 \times 10^{-5} $, and with a velocity dispersion of $215.86 \pm 6.34$ km/s. From archival radio observations we see that this galaxy hosts a compact flat-spectrum radio source, and we conclude that it is the active galactic nucleus (AGN) responsible for generating the radio lobes. The radio luminosity distribution of the jets, and the broad-line classification of the host AGN, indicate this GRG is orientated well out of the plane of the sky, making its physical size one of the largest known for any GRG. Analysis of the infrared data suggests that the host is a lenticular type galaxy with a large stellar mass ($\log~\mathrm{M}/\mathrm{M}_\odot = 11.56 \pm 0.12$), and a moderate star formation rate ($1.2 \pm 0.3~\mathrm{M}_\odot/\mathrm{year}$). Spatially smoothing the SDSS images shows the system around UGC 9555 to be significantly disturbed, with a prominent extension to the south-east. Overall, the evidence suggests this host galaxy has undergone one or more recent moderate merger events and is also experiencing tidal interactions with surrounding galaxies, which have caused the star formation and provided the supply of gas to trigger and fuel the Mpc-scale radio lobes.

We have observed the steep spectrum radio source B2 0924+30 using the LOw Frequency ARray (LOFAR) telescope. Hosted by a z=0.026 elliptical galaxy, it has a relatively large angular size of 12' (corresponding to 360 kpc projected linear size) and a morphology reminiscent of a remnant Fanaroff-Riley type II (FRII) radio galaxy. We combine LOFAR and archival images to study the spectral properties at a spatial resolution of 1'. We derive low frequency spectral index maps and use synchrotron ageing models to infer ages for different regions of the source. Using radiative ageing model fitting we show that the AGN activity ceased around 50 Myr ago. We note that the outer regions of the lobes are younger than the inner regions which is interpreted as a sign that those regions are remnant hotspots. We demonstrate the usefulness of maps of AGN radio remnants taken at low frequencies and suggest caution over the interpretation of spectral ages derived from integrated flux density measurements versus age mapping. The spectral index properties as well as the derived ages of B2 0924+30 are consistent with it being an FRII AGN radio remnant. LOFAR data are proving to be instrumental in extending our studies to the lowest radio frequencies and enabling analysis of the oldest source regions.

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

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

(abridged). We outline LBCS (the LOFAR Long-Baseline Calibrator Survey), whose aim is to identify sources suitable for calibrating the highest-resolution observations made with the International LOFAR Telescope, which include baselines >1000 km. Suitable sources must contain significant correlated flux density (50-100mJy) at frequencies around 110--190~MHz on scales of a few hundred mas. At least for the 200--300-km international baselines, we find around 1 suitable calibrator source per square degree over a large part of the northern sky, in agreement with previous work. This should allow a randomly selected target to be successfully phase calibrated on the international baselines in over 50% of cases. Products of the survey include calibrator source lists and fringe-rate and delay maps of wide areas -- typically a few degrees -- around each source. The density of sources with significant correlated flux declines noticeably with baseline length over the range 200--600~km, with good calibrators on the longest baselines appearing only at the rate of 0.5 per square degree. Coherence times decrease from 1--3 minutes on 200-km baselines to about 1 minute on 600-km baselines, suggesting that ionospheric phase variations contain components with scales of a few hundred kilometres. The longest median coherence time, at just over 3 minutes, is seen on the DE609 baseline, which at 227km is close to being the shortest. We see median coherence times of between 80 and 110 seconds on the four longest baselines (580--600~km), and about 2 minutes for the other baselines. The success of phase transfer from calibrator to target is shown to be influenced by distance, in a manner that suggests a coherence patch at 150-MHz of the order of 1 degree.

Due to their steep spectra, low-frequency observations of FR II radio galaxies potentially provide key insights in to the morphology, energetics and spectrum of these powerful radio sources. However, limitations imposed by the previous generation of radio interferometers at metre wavelengths has meant that this region of parameter space remains largely unexplored. In this paper, the first in a series examining FR IIs at low frequencies, we use LOFAR observations between 50 and 160 MHz, along with complementary archival radio and X-ray data, to explore the properties of two FR II sources, 3C452 and 3C223. We find that the morphology of 3C452 is that of a standard FR II rather than of a double-double radio galaxy as had previously been suggested, with no remnant emission being observed beyond the active lobes. We find that the low-frequency integrated spectra of both sources are much steeper than expected based on traditional assumptions and, using synchrotron/inverse-Compton model fitting, show that the total energy content of the lobes is greater than previous estimates by a factor of around 5 for 3C452 and 2 for 3C223. We go on to discuss possible causes of these steeper than expected spectra and provide revised estimates of the internal pressures and magnetic field strengths for the intrinsically steep case. We find that the ratio between the equipartition magnetic field strengths and those derived through synchrotron/inverse-Compton model fitting remains consistent with previous findings and show that the observed departure from equipartition may in some cases provide a solution to the spectral versus dynamical age disparity.

We present the results of a four-month campaign searching for low-frequency radio transients near the North Celestial Pole with the Low-Frequency Array (LOFAR), as part of the Multifrequency Snapshot Sky Survey (MSSS). The data were recorded between 2011 December and 2012 April and comprised 2149 11-minute snapshots, each covering 175 deg^2. We have found one convincing candidate astrophysical transient, with a duration of a few minutes and a flux density at 60 MHz of 15-25 Jy. The transient does not repeat and has no obvious optical or high-energy counterpart, as a result of which its nature is unclear. The detection of this event implies a transient rate at 60 MHz of 3.9 (+14.7, -3.7) x 10^-4 day^-1 deg^-2, and a transient surface density of 1.5 x 10^-5 deg^-2, at a 7.9-Jy limiting flux density and ~10-minute time-scale. The campaign data were also searched for transients at a range of other time-scales, from 0.5 to 297 min, which allowed us to place a range of limits on transient rates at 60 MHz as a function of observation duration.

Compact Steep Spectrum, Gigahertz Peaked Spectrum and High Frequency Peak (CSS, GPS, HFP) sources are considered to be young radio sources but the details of their duty cycle are not well understood. In some cases they are thought to develop in large radio galaxies, while in other cases their jets may experience intermittent activity or die prematurely and remain confined within the host galaxy. By studying in a systematic way the presence and the properties of any extended emission surrounding these compact sources we can provide firmer constraints on their evolutionary history and on the timescales of activity of the radio source. Remnant emission from previous outbursts is supposed to have very low surface brightness and to be brighter at low frequency. Taking advantage of the unprecedented sensitivity and resolution provided by the Low Frequency Array (LOFAR) we have started a systematic search of new CSS, GPS and HFP sources with extended emission, as well as a more detailed study of some well-known of these sources. Here we present the key points of our search in the LOFAR fields and a more in-depth analysis on the source B2 0258+35, a CSS source surrounded by a pair of large, diffuse radio lobes.

Remnant radio galaxies represent the final dying phase of radio galaxy evolution, in which the jets are no longer active. Due to their rarity in flux limited samples and the difficulty of identification, this dying phase remains poorly understood and the luminosity evolution largely unconstrained. Here we present the discovery and detailed analysis of a large (700 kpc) remnant radio galaxy with a low surface brightness that has been identified in LOFAR images at 150 MHz. By combining LOFAR data with new follow-up Westerbork observations and archival data at higher frequencies, we investigated the source morphology and spectral properties from 116 to 4850 MHz. By modelling the radio spectrum we probed characteristic timescales of the radio activity. The source has a relatively smooth, diffuse, amorphous appearance together with a very weak central compact core which is associated with the host galaxy located at z=0.051. From our ageing and morphological analysis it is clear that the nuclear engine is currently switched off or, at most, active at a very low power state. The host galaxy is currently interacting with another galaxy located at a projected separation of 15 kpc and a radial velocity offset of 300 km/s. This interaction may have played a role in the triggering and/or shut down of the radio jets. The spectral shape of this remnant radio galaxy differs from the majority of the previously identified remnant sources, which show steep or curved spectra at low to intermediate frequencies. In light of this finding and in preparation for new-generation deep low-frequency surveys, we discuss the selection criteria to be used to select representative samples of these sources.

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.

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.

We present radio observations of the Moon between $35$ and $80$ MHz to demonstrate a novel technique of interferometrically measuring large-scale diffuse emission extending far beyond the primary beam (global signal) for the first time. In particular, we show that (i) the Moon appears as a negative-flux source at frequencies $35<\nu<80$ MHz since it is `colder' than the diffuse Galactic background it occults, (ii) using the (negative) flux of the lunar disc, we can reconstruct the spectrum of the diffuse Galactic emission with the lunar thermal emission as a reference, and (iii) that reflected RFI (radio-frequency interference) is concentrated at the center of the lunar disc due to specular nature of reflection, and can be independently measured. Our RFI measurements show that (i) Moon-based Cosmic Dawn experiments must design for an Earth-isolation of better than $80$ dB to achieve an RFI temperature $<1$ mK, (ii) Moon-reflected RFI contributes to a dipole temperature less than $20$ mK for Earth-based Cosmic Dawn experiments, (iii) man-made satellite-reflected RFI temperature exceeds $20$ mK if the aggregate cross section of visible satellites exceeds $80$ m$^2$ at $800$ km height, or $5$ m$^2$ at $400$ km height. Currently, our diffuse background spectrum is limited by sidelobe confusion on short baselines (10-15% level). Further refinement of our technique may yield constraints on the redshifted global $21$-cm signal from Cosmic Dawn ($40>z>12$) and the Epoch of Reionization ($12>z>5$).

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

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.

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.

Outlined is the discovery of a very faint, diffuse, low surface-brightness (0.5 \mJybeam, 1.4 \mJyarcminsq on average) structure around the radio source B2 0258+35 hosted by an HI-rich early-type galaxy (NGC 1167). Since B2 0258+35 is a young Compact Steep Spectrum (CSS) source, the newly discovered structure could represent a remnant from an earlier stage of AGN activity. We go on by explaining in detail all the possibilities for triggering the radio activity in B2 0258+35 regarding gas accretion in a recurrent AGN activity framework. NGC 1167 hosts a very regular, extended and massive \HI\ disc that has been studied in great detail. Previous studies of the \HI\ closer to the core seem to go against the assumption of a circum-nuclear disc of \HI\ as the source of the accreting gas. We consider the cooling of gas from the hot, X-ray halo as a possible alternative option for the fueling of the AGN, as suggested in the case of other sources of similar radio power as B2 0258+35. Estimates are given for the age of the faint diffuse emission as well as for the current accretion rate, which are in good agreement with literature values. If our assumptions about the accretion mechanism are correct, similar large-scale, relic-like structures should be more commonly found around early-type galaxies and this will be hopefully confirmed by the next generation of sensitive, low-frequency radio surveys.

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.

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.

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.

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.