In this data release from the ongoing LOw-Frequency ARray (LOFAR) Two-metre Sky Survey (LoTSS) we present 120-168 MHz images covering 27% of the northern sky. Our coverage is split into two regions centred at approximately 12h45m +44 • 30 and 1h00m +28 • 00 and spanning 4178 and 1457 square degrees respectively. The images were derived from 3,451 hrs (7.6 PB) of LOFAR High Band Antenna data which were corrected for the direction-independent instrumental properties as well as direction-dependent ionospheric distortions during extensive, but fully automated, data processing. A catalogue of 4,396,228 radio sources is derived from our total intensity (Stokes I) maps, where the majority of these have never been detected at radio wavelengths before. At 6 resolution, our full bandwidth Stokes I continuum maps with a central frequency of 144 MHz have: a median rms sensitivity of 83 µJy/beam; a flux density scale accuracy of approximately 10%; an astrometric accuracy of 0.2 ; and we estimate the point-source completeness to be 90% at a peak brightness of 0.8 mJy/beam. By creating three 16 MHz bandwidth images across the band we are able to measure the in-band spectral index of many sources, albeit with an error on the derived spectral index of > ±0.2 which is a consequence of our flux-density scale accuracy and small fractional bandwidth. Our circular polarisation (Stokes V) 20 resolution 120-168 MHz continuum images have a median rms sensitivity of 95 µJy/beam, and we estimate a Stokes I to Stokes V leakage of 0.056%. Our linear polarisation (Stokes Q and Stokes U) image cubes consist of 480 × 97.6 kHz wide planes and have a median rms sensitivity per plane of 10.8 mJy/beam at 4 and 2.2 mJy/beam at 20 ; we estimate the Stokes I to Stokes Q/U leakage to be approximately 0.2%. Here we characterise and publicly release our Stokes I, Q, U and V images in addition to the calibrated uv-data to facilitate the thorough scientific exploitation of this unique dataset.

We present a study of the stellar host galaxy, CO (1–0) molecular gas distribution and AGN emission on 50–500 pc-scales of the gravitationally lensed dust-obscured AGN MG J0751+2716 and JVAS B1938+666 at redshifts 3.200 and 2.059, respectively. By correcting for the lensing distortion using a grid-based lens modelling technique, we spatially locate the different emitting regions in the source plane for the first time. Both AGN host galaxies have 300–500 pc-scale size and surface brightness consistent with a bulge/pseudo-bulge, and 2 kpc-scale AGN radio jets that are embedded in extended molecular gas reservoirs that are 5–20 kpc in size. The CO (1–0) velocity fields show structures possibly associated with discs (elongated velocity gradients) and interacting objects (off-axis velocity components). There is evidence for a decrement in the CO (1–0) surface brightness at the location of the host galaxy, which may indicate radiative feedback from the AGN, or offset star formation. We find CO–H2 conversion factors of around αCO = 1.5 ± 0.5 (K km s−1 pc2)−1, molecular gas masses of >3 × 1010 M⊙, dynamical masses of ∼1011 M⊙, and gas fractions of around 60 per cent. The intrinsic CO line luminosities are comparable to those of unobscured AGN and dusty star-forming galaxies at similar redshifts, but the infrared luminosities are lower, suggesting that the targets are less efficient at forming stars. Therefore, they may belong to the AGN feedback phase predicted by galaxy formation models, because they are not efficiently forming stars considering their large amount of molecular gas.

No abstract

The International LOFAR Telescope is an interferometer with stations spread across Europe. With baselines of up to ∼ 2,000 km, LOFAR has the unique capability of achieving sub-arcsecond resolution at frequencies below 200 MHz. However, it is technically and logistically challenging to process LOFAR data at this resolution. To date only a handful of publications have exploited this capability. Here we present a calibration strategy that builds on previous high-resolution work with LOFAR. It is implemented in a pipeline using mostly dedicated LOFAR software tools and the same processing framework as the LOFAR Two-metre Sky Survey (LoTSS). We give an overview of the calibration strategy and discuss the special challenges inherent to enacting high-resolution imaging with LOFAR, and describe the pipeline, which is publicly available, in detail. We demonstrate the calibration strategy by using the pipeline on P205+55, a typical LoTSS pointing with an 8 hour observation and 13 international stations. We perform in-field delay calibration, solution referencing to other calibrators in the field, self-calibration of these calibrators, and imaging of example directions of interest in the field. We find that for this specific field and these ionospheric conditions, dispersive delay solutions can be transferred between calibrators up to ∼1.5 degrees away, while phase solution transferral works well over ∼1 degree. We also demonstrate a check of the astrometry and flux density scale with the in-field delay calibrator source. Imaging in 17 directions, we find the restoring beam is typically ∼0.3 ×0.2 although this varies slightly over the entire 5 square degree field of view. We find we can achieve ∼80 to 300 µJy bm −1 image rms noise, which is dependent on the distance from the phase centre; typical values are ∼ 90 µJy bm −1 for the 8 hour observation with 48 MHz of bandwidth. Seventy percent of processed sources are detected, and from this we estimate that we should be able to image roughly 900 sources per LoTSS pointing. This equates to ∼3 million sources in the northern sky, which LoTSS will entirely cover in the next several years. Future optimisation of the calibration strategy for efficient post-processing of LoTSS at high resolution (LoTSS-HR) makes this estimate a lower limit.

Active galactic nuclei (AGNs) show episodic activity, which can be evident in galaxies that exhibit restarted radio jets. These restarted jets can interact with their environment, leaving signatures on the radio spectral energy distribution. Tracing these signatures is a powerful way to explore the life of radio galaxies. This requires resolved spectral index measurements over a broad frequency range including low frequencies. We present such a study for the radio galaxy 3C 293, which has long been thought to be a restarted galaxy on the basis of its radio morphology. Using the International LOFAR telescope (ILT) we probed spatial scales as fine as ∼0.2 at 144 MHz, and to constrain the spectrum we combined these data with Multi-Element Radio Linked Interferometer Network (MERLIN) and Very Large Array (VLA) archival data at frequencies up to 8.4 GHz that have a comparable resolution. In the inner lobes (∼2 kpc), we detect the presence of a spectral turnover that peaks at ∼225 MHz and is most likely caused by free-free absorption from the rich surrounding medium. We confirm that these inner lobes are part of a jet-dominated young radio source (spectral age 0.17 Myr), which is strongly interacting with the rich interstellar medium (ISM) of the host galaxy. The diffuse emission surrounding these lobes on scales of up to ∼4.5 kpc shows steeper spectral indices (∆α ∼ 0.2-0.5, S ∝ ν −α ) and a spectral age of 0.27 Myr. The outer lobes (extending up to ∼100 kpc) have a spectral index of α ∼ 0.6-0.8 from 144-4850 MHz with a remarkably uniform spatial distribution and only mild spectral curvature (∆α 0.2). We propose that intermittent fuelling and jet flow disruptions are powering the mechanisms that keep the spectral index in the outer lobes from steepening and maintain the spatial uniformity of the spectral index. Overall, it appears that 3C 293 has gone through multiple (two to three) epochs of activity. This study adds 3C 293 to the new sub-group of restarted galaxies with short interruption time periods. This is the first time a spatially resolved study has been performed that simultaneously studies a young source as well as the older outer lobes at such low frequencies. This illustrates the potential of the International LOFAR telescope to expand such studies to a larger sample of radio galaxies.

The prominent radio source Hercules A features complex structures in its radio lobes. Although it is one of the most comprehensively studied sources in the radio sky, the origin of the ring structures in the Hercules A radio lobes remains an open question. We present the first sub-arcsecond angular resolution images at low frequencies (<300 MHz) of Hercules A, made with the International LOFAR Telescope. With the addition of data from the Karl G. Jansky Very Large Array, we mapped the structure of the lobes from 144 MHz to 7 GHz. We explore the origin of the rings within the lobes of Hercules A, and test whether their properties are best described by a shock model, where shock waves are produced by the jet propagating in the radio lobe, or by an inner-lobe model, where the rings are formed by decelerated jetted plasma. From spectral index mapping our large frequency coverage reveals that the curvature of the different ring spectra increases with distance away from the central active galactic nucleus. We demonstrate that the spectral shape of the rings is consistent with synchrotron aging, which speaks in favor of an inner-lobe model where the rings are formed from the deposition of material from past periods of intermittent core activity.

The Low-Frequency Array (LOFAR) Long-Baseline Calibrator Survey (LBCS) was conducted between 2014 and 2019 in order to obtain a set of suitable calibrators for the LOFAR array. In this paper we present the complete survey, building on the preliminary analysis published in 2016 which covered approximately half the survey area. The final catalogue consists of 30006 observations of 24713 sources in the northern sky, selected for a combination of high low-frequency radio flux density and flat spectral index using existing surveys (WENSS, NVSS, VLSS, and MSSS). Approximately one calibrator per square degree, suitable for calibration of ≥ 200 km baselines is identified by the detection of compact flux density, for declinations north of 30 • and away from the Galactic plane, with a considerably lower density south of this point due to relative difficulty in selecting flat-spectrum candidate sources in this area of the sky. The catalogue contains indicators of degree of correlated flux on baselines between the Dutch core and each of the international stations, involving a maximum baseline length of nearly 2000 km, for all of the observations. Use of the VLBA calibrator list, together with statistical arguments by comparison with flux densities from lower-resolution catalogues, allow us to establish a rough flux density scale for the LBCS observations, so that LBCS statistics can be used to estimate compact flux densities on scales between 300 mas and 2 , for sources observed in the survey. The survey is used to estimate the phase coherence time of the ionosphere for the LOFAR international baselines, with median phase coherence times of about 2 minutes varying by a few tens of percent between the shortest and longest baselines. The LBCS can be used to assess the structures of point sources in lower-resolution surveys, with significant reductions in the degree of coherence in these sources on scales between 2 and 300 mas. The LBCS survey sources show a greater incidence of compact flux density in quasars than in radio galaxies, consistent with unified schemes of radio sources. Comparison with samples of sources from interplanetary scintillation (IPS) studies with the Murchison Widefield Array (MWA) shows consistent patterns of detection of compact structure in sources observed both interferometrically with LOFAR and using IPS.

It has long been known that radio sources with the steepest spectra are preferentially associated with the most distant galaxies, the α − z relation, but the reason for this relation is an open question. The spatial distribution of spectra in high-z radio sources can be used to study this relation, and low-frequency observations are particularly important in understanding the particle acceleration and injection mechanisms. However, the small angular sizes of high-z sources together with the inherently low resolution of lowfrequency radio telescopes until now has prevented high angular resolution low-frequency observations of distant objects. Here we present subarcsecond observations of a z = 2.4 radio galaxy at frequencies between 121 MHz and 166 MHz. We measure the spatial distribution of spectra, and discuss the implications for models of the α − z relation. We targeted 4C 43.15 with the High Band Antennas (HBAs) of the International LOFAR Telescope (ILT) with a range of baselines up to 1300 km. At the central frequency of 143 MHz we achieve an angular resolution of ∼ 0.3 . By complementing our data with archival Very Large Array (VLA) data we study the spectral index distribution across 4C 43.15 between 55 MHz and 8.4 GHz at resolutions of 0.4 and 0.9 . With a magnetic field strength of B = 5.2 nT and fitted injection indices of α north inj = −0.8 and α south inj = −0.6, fitting a Tribble spectral ageing model results in a spectral age of τ spec = 1.1 ± 0.1 Myr. We conclude that our data on 4C 43.15 indicates that inverse Compton losses could become comparable to or exceed synchrotron losses at higher redshifts and that inverse Compton losses could be a viable explanation for the α − z relation. Statistical studies of these objects will become possible in the future when wide-area subarcsecond surveys start.