Search

The turbulent ionosphere causes phase shifts to incoming radio waves on a broad range of temporal and spatial scales. When an interferometer is not sufficiently calibrated for the direction-dependent ionospheric effects, the time-varying phase shifts can cause the signal to decorrelate. The ionosphere's influence over various spatiotemporal scales introduces a baseline-dependent effect on the interferometric array. We study the impact of baseline-dependent decorrelation on high-redshift observations with the Low Frequency Array (LOFAR). Datasets with a range of ionospheric corruptions are simulated using a thin-screen ionosphere model, and calibrated using the state-of-the-art LOFAR Epoch of Reionisation pipeline. For the first time ever, we show the ionospheric impact on various stages of the calibration process including an analysis of the transfer of gain errors from longer to shorter baselines using realistic end-to-end simulations. We find that direction-dependent calibration for source subtraction leaves excess power of up to two orders of magnitude above the thermal noise at the largest spectral scales in the cylindrically averaged auto-power spectrum under normal ionospheric conditions. However, we demonstrate that this excess power can be removed through Gaussian process regression, leaving no excess power above the ten per cent level for a $5~$km diffractive scale. We conclude that ionospheric errors, in the absence of interactions with other aggravating effects, do not constitute a dominant component in the excess power observed in LOFAR Epoch of Reionisation observations of the North Celestial Pole. Future work should therefore focus on less spectrally smooth effects, such as beam modelling errors., Comment: 22 pages, 20 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS)

We present the deepest wide-field 115-166 MHz image at sub-arcsecond resolution spanning an area of 2.5 by 2.5 degrees centred at the ELAIS-N1 deep field. To achieve this, we improved the calibration for the International LOFAR Telescope. This enhancement enabled us to efficiently process 32 hrs of data from four different 8-hr observations using the high-band antennas (HBAs) of all 52 stations, covering baselines up to approximately 2,000 km across Europe. The DI calibration was improved by using an accurate sky model and refining the series of calibration steps on the in-field calibrator, while the DD calibration was improved by adopting a more automated approach for selecting the DD calibrators and inspecting the self-calibration on these sources. We also added an additional round of self-calibration for the Dutch core and remote stations in order to refine the solutions for shorter baselines. To complement our highest resolution at 0.3", we also made intermediate resolution wide-field images at 0.6" and 1.2". Our resulting wide-field images achieve a central noise level of 14 muJy/beam at 0.3", doubling the depth and uncovering four times more objects than the Lockman Hole deep field image at comparable resolution but with only 8 hrs of data. Compared to LOFAR imaging without the international stations, we note that due to the increased collecting area and the absence of confusion noise, we reached a point-source sensitivity comparable to a 500-hr ELAIS-N1 6" image with 16 times less observing time. Importantly, we have found that the computing costs for the same amount of data are almost halved (to about 139,000 CPU hrs per 8 hrs of data) compared to previous efforts, though they remain high. Our work underscores the value and feasibility of exploiting all Dutch and international LOFAR stations to make deep wide-field images at sub-arcsecond resolution., Comment: Images and catalogues will be available on https://lofar-surveys.org/hd-en1.html

Current radio interferometers output multi-petabyte-scale volumes of data per year making the storage, transfer, and processing of this data a sizeable challenge. This challenge is expected to grow with the next-generation telescopes such as the Square Kilometre Array. Lossy compression of interferometric data post-correlation can abate this challenge. However, since high-redshift 21-cm studies impose strict precision requirements, the impact of such lossy data compression on the 21-cm signal power spectrum statistic should be understood. We apply Dysco visibility compression, a technique to normalize and quantize specifically designed for radio interferometric data. We establish the level of the compression noise in the power spectrum in comparison to the thermal noise as well as its coherency behavior. Finally, for optimal compression results, we compare the compression noise obtained from different compression settings to a nominal 21-cm signal power. From a single night of observation, we find that the noise introduced due to the compression is more than five orders of magnitude lower than the thermal noise level in the power spectrum. The noise does not affect calibration. The compression noise shows no correlation with the sky signal and has no measurable coherent component. The level of compression error in the power spectrum ultimately depends on the compression settings. Dysco visibility compression is found to be of insignificant concern for 21-cm power spectrum studies. Hence, data volumes can be safely reduced by factors of $\sim 4$ and with insignificant bias to the final power spectrum. Data from SKA-low will likely be compressible by the same factor as LOFAR, owing to the similarities of the two instruments. The same technique can be used to compress data from other telescopes, but a small adjustment of the compression parameters might be required., Comment: 12 pages, 9 figures, and 3 tables; submitted to Astronomy and Astrophysics (A&A)

One of the main obstacles preventing the detection of the redshifted 21-cm signal from neutral hydrogen in the early Universe is the astrophysical foreground emission, which is several orders of magnitude brighter than the signal. The foregrounds, due to their smooth spectra, are expected to predominantly occupy a region in the cylindrical power spectrum known as the foreground wedge. However, the conventional equations describing the extent of the foreground wedge are derived under a flat-sky approximation. This assumption breaks down for tracking wide-field instruments, thus rendering these equations inapplicable in these situations. In this paper, we derive equations for the full sky foreground wedge and show that the foregrounds can potentially extend far beyond what the conventional equations suggest. We also derive the equations that describe a specific bright source in the cylindrical power spectrum space. The validity of both sets of equations is tested against numerical simulations. Many current and upcoming interferometers (e.g., LOFAR, NenuFAR, MWA, SKA) are wide-field phase-tracking instruments. These equations give us new insights into the nature of foreground contamination in the cylindrical power spectra estimated using wide-field instruments. Additionally, they allow us to accurately associate features in the power spectrum to foregrounds or instrumental effects. The equations are also important for correctly selecting the "EoR window" for foreground avoidance analyses, and for planning 21-cm observations. In future analyses, it is recommended to use these updated horizon lines to indicate the foreground wedge in the cylindrical power spectrum accurately. The new equations for generating the updated wedge lines are made available in a Python library, pslines., Comment: 20 pages, 13 figures, and 2 tables; accepted for publication in Astronomy and Astrophysics (A&A)

The origin of the radio synchrotron background (RSB) is currently unknown. Its understanding might have profound implications in fundamental physics or might reveal a new class of radio emitters. In this work, we consider the scenario in which the RSB is due to extragalactic radio sources and measure the angular cross-correlation of LOFAR images of the diffuse radio sky with matter tracers at different redshifts, provided by galaxy catalogs and CMB lensing. We compare these measured cross-correlations to those expected for models of RSB sources. We find that low-redshift populations of discrete sources are excluded by the data, while higher redshift explanations are compatible with available observations. We also conclude that at least 20\% of the RSB surface brightness level must originate from populations tracing the large-scale distribution of matter in the universe, indicating that at least this fraction of the RSB is of extragalactic origin. Future measurements of the correlation between the RSB and tracers of high-redshift sources will be crucial to constraining the source population of the RSB., Comment: 12 pages, 8 figures. Version accepted for publication in MNRAS

The redshifted 21 cm signal from neutral hydrogen is a direct probe of the physics of the early universe and has been an important science driver of many present and upcoming radio interferometers. In this study we use a single night of observations with the New Extension in Nan\c{c}ay Upgrading LOFAR (NenuFAR) to place upper limits on the 21 cm power spectrum from cosmic dawn at a redshift of $z$ = 20.3. NenuFAR is a new low-frequency radio interferometer, operating in the 10-85 MHz frequency range, currently under construction at the Nan\c{c}ay Radio Observatory in France. It is a phased array instrument with a very dense uv coverage at short baselines, making it one of the most sensitive instruments for 21 cm cosmology analyses at these frequencies. Our analysis adopts the foreground subtraction approach, in which sky sources are modeled and subtracted through calibration and residual foregrounds are subsequently removed using Gaussian process regression. The final power spectra are constructed from the gridded residual data cubes in the uv plane. Signal injection tests are performed at each step of the analysis pipeline, the relevant pipeline settings are optimized to ensure minimal signal loss, and any signal suppression is accounted for through a bias correction on our final upper limits. We obtain a best 2$\sigma$ upper limit of $2.4\times 10^7$ $\text{mK}^{2}$ at $z$ = 20.3 and $k$ = 0.041 $h\,\text{cMpc}^{-1}$. We see a strong excess power in the data, making our upper limits two orders of magnitude higher than the thermal noise limit. We investigate the origin and nature of this excess power and discuss further improvements to the analysis pipeline that can potentially mitigate it and consequently allow us to reach thermal noise sensitivity when multiple nights of observations are processed in the future., Comment: 27 pages, 21 figures, and 6 tables; accepted for publication in Astronomy and Astrophysics (A&A); language edits implemented; typos corrected

Measurement of the redshifted 21-cm signal of neutral hydrogen from the Cosmic Dawn (CD) and Epoch of Reionisation (EoR) promises to unveil a wealth of information about the astrophysical processes during the first billion years of evolution of the universe. The AARTFAAC Cosmic Explorer (ACE) utilises the AARTFAAC wide-field imager of LOFAR to measure the power spectrum of the intensity fluctuations of the redshifted 21-cm signal from the CD at z~18. The RFI from various sources contaminates the observed data and it is crucial to exclude the RFI-affected data in the analysis for reliable detection. In this work, we investigate the impact of non-ground-based transient RFI using cross-power spectra and cross-coherence metrics to assess the correlation of RFI over time and investigate the level of impact of transient RFI on the ACE 21-cm power spectrum estimation. We detected moving sky-based transient RFI sources that cross the field of view within a few minutes and appear to be mainly from aeroplane communication beacons at the location of the LOFAR core in the 72-75 MHz band, by inspecting filtered images. This transient RFI is mostly uncorrelated over time and is only expected to dominate over the thermal noise for an extremely deep integration time of 3000 hours or more with a hypothetical instrument that is sky temperature dominated at 75 MHz. We find no visible correlation over different k-modes in Fourier space in the presence of noise for realistic thermal noise scenarios. We conclude that the sky-based transient RFI from aeroplanes, satellites and meteorites at present does not pose a significant concern for the ACE analyses at the current level of sensitivity and after integrating over the available 500 hours of observed data. However, it is crucial to mitigate or filter such transient RFI for more sensitive experiments aiming for significantly deeper integration., Comment: 13 pages, 9 figures, and 3 tables; accepted for publication in Astronomy and Astrophysics (A&A)

We present a new method, called "forced-spectrum fitting", for physically-based spectral modelling of radio sources during deconvolution. This improves upon current common deconvolution fitting methods, which often produce inaccurate spectra. Our method uses any pre-existing spectral index map to assign spectral indices to each model component cleaned during the multi-frequency deconvolution of WSClean, where the pre-determined spectrum is fitted. The component magnitude is evaluated by performing a modified weighted linear least-squares fit. We test this method on a simulated LOFAR-HBA observation of the 3C196 QSO and a real LOFAR-HBA observation of the 4C+55.16 FRI galaxy. We compare the results from the forced-spectrum fitting with traditional joined-channel deconvolution using polynomial fitting. Because no prior spectral information was available for 4C+55.16, we demonstrate a method for extracting spectral indices in the observed frequency band using "clustering". The models generated by the forced-spectrum fitting are used to improve the calibration of the datasets. The final residuals are comparable to existing multi-frequency deconvolution methods, but the output model agrees with the provided spectral index map, embedding correct spectral information. While forced-spectrum fitting does not solve the determination of the spectral information itself, it enables the construction of accurate multi-frequency models that can be used for wide-band calibration and subtraction., Comment: 17 pages, 9 figures, 5 tables. Accepted for publication in MNRAS

We present a new image of a 26.5 square degree region in the Bo\"otes constellation obtained at 1.4 GHz using the Aperture Tile in Focus (Apertif) system on the Westerbork Synthesis Radio Telescope. We use a newly developed processing pipeline which includes direction-dependent self-calibration which provides a significant improvement of the quality of the images compared to those released as part of the Apertif first data release. For the Bo\"otes region, we mosaic 187 Apertif images and extract a source catalog. The mosaic image has an angular resolution of 27${\times}$11.5 arcseconds and a median background noise of 40 ${\mu}$Jy/beam. The catalog has 8994 sources and is complete down to the 0.3 mJy level. We combine the Apertif image with LOFAR images of the Bo\"otes field at 54 and 150 MHz to study spectral properties of the sources. We find a spectral flattening towards low flux density sources. Using the spectral index limits from Apertif non-detections we derive that up to 9 percent of the sources have ultra-steep spectra with a slope steeper than -1.2. Steepening of the spectral index with increasing redshift is also seen in the data showing a different dependency for the low-frequency spectral index and the high frequency one. This can be explained by a population of sources having concave radio spectra with a turnover frequency around the LOFAR band. Additionally, we discuss cases of individual extended sources with an interesting resolved spectral structure. With the improved pipeline, we aim to continue processing data from the Apertif wide-area surveys and release the improved 1.4 GHz images of several famous fields., Comment: 13 pages, 9 figures; to be published in A&A

We present the largest low frequency (120~MHz) arcminute resolution image of the radio synchrotron background (RSB) to date, and its corresponding angular power spectrum of anisotropies (APS) with angular scales ranging from $3^\circ$ to $0.3^\prime$. We show that the RSB around the North Celestial Pole has a significant excess anisotropy power at all scales over a model of unclustered point sources based on source counts of known source classes. This anisotropy excess, which does not seem attributable to the diffuse Galactic emission, could be linked to the surface brightness excess of the RSB. To better understand the information contained within the measured APS, we model the RSB varying the brightness distribution, size, and angular clustering of potential sources. We show that the observed APS could be produced by a population of faint clustered point sources only if the clustering is extreme and the size of the Gaussian clusters is $\lesssim 1'$. We also show that the observed APS could be produced by a population of faint diffuse sources with sizes $\lesssim 1'$, and this is supported by features present in our image. Both of these cases would also cause an associated surface brightness excess. These classes of sources are in a parameter space not well probed by even the deepest radio surveys to date., Comment: 13 pages, 14 figures. Accepted for publication in MNRAS

Radio spectroscopy provides a unique inspection perspective for solar and space weather research, which can reveal the plasma and energetic electron information in the solar corona and inner heliosphere. However, Radio-Frequency Interference (RFI) from human activities affects sensitive radio telescopes, and significantly affects the quality of observation. Thus, RFI detection and mitigation for the observations is necessary to obtain high quality, science-ready data. The flagging of RFI is particularly challenging for the solar and space weather observations at low frequency, because the solar radio bursts can be brighter than the RFI, and may show similar temporal behavior. In this work, we investigate RFI flagging methods for solar and space weather observations, including a strategy for AOFlagger, and a novel method that makes use of a morphology convolution. These algorithms can effectively flag RFI while preserving solar radio bursts.

The Low Frequency Array (LOFAR) is the only existing radio interferometer able to observe at ultra-low frequencies (<100 MHz) with high resolution (<15") and high sensitivity (<1 mJy/beam). To exploit these capabilities, the LOFAR Surveys Key Science Project is using the LOFAR Low Band Antenna (LBA) to carry out a sensitive wide-area survey at 41-66 MHz named the LOFAR LBA Sky Survey (LoLSS). LoLSS is covering the whole northern sky above declination 24 deg with a resolution of 15" and a sensitivity of 1-2 mJy/beam (1 sigma) depending on declination, field properties, and observing conditions. Here we present the first data release. An automated pipeline was used to reduce the 95 fields included in this data release. The data reduction procedures developed for this project have general application and are currently being used to process LOFAR LBA interferometric observations. Compared to the preliminary release, direction-dependent errors have been corrected for during the calibration process. This results in a typical sensitivity of 1.55 mJy/beam at the target resolution of 15". The first data release of the LOFAR LBA Sky Survey covers 650 sqdeg in the HETDEX spring field. The resultant data products released to the community include mosaic images (I and V Stokes) of the region, and a catalogue of 42463 detected sources and related Gaussian components used to describe sources' morphologies. Separate catalogues for 6 in-band frequencies are also released. The first data release of LoLSS shows that, despite the influences of the ionosphere, LOFAR can conduct large-scale surveys in the frequency window 42-66 MHz with unprecedentedly high sensitivity and resolution. The data can be used to derive unique information on the low-frequency spectral properties of many thousands of sources with a wide range of applications in extragalactic and galactic astronomy., Comment: 20 pages, 22 figures, images and catalogues available at https://www.lofar-surveys.org/lolss.html

Context. Apertif is a multi-beam receiver system for the Westerbork Synthesis Radio Telescope that operates at 1.1-1.5 GHz, which overlaps with various radio services, resulting in contamination of astronomical signals with radio-frequency interference (RFI). Aims. We analyze approaches to mitigate Apertif interference and design an automated detection procedure for its imaging mode. Using this approach, we present long-term RFI detection results of over 300 Apertif observations. Methods. Our approach is based on the AOFlagger detection approach. We introduce several new features, including ways to deal with ranges of invalid data (e.g. caused by shadowing) in both the SumThreshold and scale-invariant rank operator steps; pre-calibration bandpass calibration; auto-correlation flagging; and HI flagging avoidance. These methods are implemented in a new framework that uses the Lua language for scripting, which is new in AOFlagger version 3. Results. Our approach removes RFI fully automatically, and is robust and effective enough for further calibration and (continuum) imaging of these data. Analysis of 304 observations show an average of 11.1% of lost data due to RFI with a large spread. We observe 14.6% RFI in auto-correlations. Computationally, AOFlagger achieves a throughput of 370 MB/s on a single computing node. Compared to published machine learning results, the method is one to two orders of magnitude faster., Comment: 15 pages, 10 figures. Accepted for publication in A&A

Context. The occultation of a radio source by the plasma tail of a comet can be used to probe structure and dynamics in the tail. Such occultations are rare, and the occurrence of scintillation, due to small-scale density variations in the tail, remains somewhat controversial. Aims. A detailed observation taken with the Low-Frequency Array (LOFAR) of a serendipitous occultation of the compact radio source 3C196 by the plasma tail of comet C/2020 F3 (Neowise) is presented. 3C196 tracked almost perpendicularly behind the tail, providing a unique profile cut only a short distance downstream from the cometary nucleus itself. Methods. Interplanetary scintillation (IPS) is observed as the rapid variation of the intensity received of a compact radio source due to density variations in the solar wind. IPS in the signal received from 3C196 was observed for five hours, covering the full transit behind the plasma tail of comet C/2020 F3 (Neowise) on 16 July 2020, and allowing an assessment of the solar wind in which the comet and its tail are embedded. Results. The results reveal a sudden and strong enhancement in scintillation which is unequivocally attributable to the plasma tail. The strongest scintillation is associated with the tail boundaries, weaker scintillation is seen within the tail, and previously-unreported periodic variations in scintillation are noted, possibly associated with individual filaments of plasma. Furthermore, contributions from the solar wind and comet tail are separated to measure a sharp decrease in the velocity of material within the tail, suggesting a steep velocity shear resulting in strong turbulence along the tail boundary, Comment: Accepted for publication in Astronomy and Astrophysics, 8 pages, 9 figures

Detecting the 21-cm signal from the Epoch of Reionisation (EoR) is challenging due to the strong astrophysical foregrounds, ionospheric effects, radio frequency interference and instrumental effects. Understanding and calibrating these effects are crucial for the detection. In this work, we introduce a newly developed direction-dependent (DD) calibration algorithm DDECAL and compare its performance with an existing algorithm, SAGECAL, in the context of the LOFAR-EoR 21-cm power spectrum experiment. In our data set, the North Celestial Pole (NCP) and its flanking fields were observed simultaneously. We analyse the NCP and one of its flanking fields. The NCP field is calibrated by the standard pipeline, using SAGECAL with an extensive sky model and 122 directions, and the flanking field is calibrated by DDECAL and SAGECAL with a simpler sky model and 22 directions. Additionally, two strategies are used for subtracting Cassiopeia A and Cygnus A. The results show that DDECAL performs better at subtracting sources in the primary beam region due to the application of a beam model, while SAGECAL performs better at subtracting Cassiopeia A and Cygnus A. This indicates that including a beam model during DD calibration significantly improves the performance. The benefit is obvious in the primary beam region. We also compare the 21-cm power spectra on two different fields. The results show that the flanking field produces better upper limits compared to the NCP in this particular observation. Despite the minor differences between DDECAL and SAGECAL due to the beam application, we find that the two algorithms yield comparable 21-cm power spectra on the LOFAR-EoR data after foreground removal. Hence, the current LOFAR-EoR 21-cm power spectrum limits are not likely to depend on the DD calibration method., Comment: 28 pages, 14 figures, accepted for publication in A&A; a few equations in section 3 were corrected in this version

(Abridged) Apertif is a phased-array feed system for WSRT, providing forty instantaneous beams over 300 MHz of bandwidth. A dedicated survey program started on 1 July 2019, with the last observations taken on 28 February 2022. We describe the release of data products from the first year of survey operations, through 30 June 2020. We focus on defining quality control metrics for the processed data products. The Apertif imaging pipeline, Apercal, automatically produces non-primary beam corrected continuum images, polarization images and cubes, and uncleaned spectral line and dirty beam cubes for each beam of an Apertif imaging observation. For this release, processed data products are considered on a beam-by-beam basis within an observation. We validate the continuum images by using metrics that identify deviations from Gaussian noise in the residual images. If the continuum image passes validation, we release all processed data products for a given beam. We apply further validation to the polarization and line data products. We release all raw observational data from the first year of survey observations, for a total of 221 observations of 160 independent target fields, covering approximately one thousand square degrees of sky. Images and cubes are released on a per beam basis, and 3374 beams are released. The median noise in the continuum images is 41.4 uJy/bm, with a slightly lower median noise of 36.9 uJy/bm in the Stokes V polarization image. The median angular resolution is 11.6"/sin(Dec). The median noise for all line cubes, with a spectral resolution of 36.6 kHz, is 1.6 mJy/bm, corresponding to a 3-sigma HI column density sensitivity of 1.8 x 10^20 atoms cm^-2 over 20 km/s (for a median angular resolution of 24" x 15"). We also provide primary beam images for each individual Apertif compound beam. The data are made accessible using a Virtual Observatory interface., Comment: Accepted for publication in A&A, updated Figure 1

The first data release of Apertif survey contains 3074 radio continuum images covering a thousand square degrees of the sky. The observations were performed during August 2019 to July 2020. The continuum images were produced at a central frequency 1355 MHz with the bandwidth of $\sim$150 MHz and angular resolution reaching 10". In this work we introduce and apply a new method to obtain a primary beam model using a machine learning approach, Gaussian process regression. The primary beam models obtained with this method are published along with the data products for the first Apertif data release. We apply the method to the continuum images, mosaic them and extract the source catalog. The catalog contains 249672 radio sources many of which are detected for the first time at these frequencies. We cross-match the coordinates with the NVSS, LOFAR/DR1/value-added and LOFAR/DR2 catalogs resulting in 44523, 22825 and 152824 common sources respectively. The first sample provides a unique opportunity to detect long term transient sources which have significantly changed their flux density for the last 25 years. The second and the third ones combined together provide information about spectral properties of the sources as well as the redshift estimates., Comment: 12 pages, 9 figures; accepted for publication in A&A

The detection of the 21 cm signal of neutral hydrogen from the Epoch of Reionization (EoR) is challenging due to bright foreground sources, radio frequency interference (RFI), the ionosphere, and instrumental effects. Even after correcting for these effects in the calibration step and applying foreground removal techniques, the remaining residuals in the observed 21 cm power spectra are still above the thermal noise, which is referred to as the "excess variance." We study potential causes of this excess variance based on 13 nights of data obtained with the Low-Frequency Array (LOFAR). We focused on the impact of gain errors, the sky model, and ionospheric effects on the excess variance by correlating the relevant parameters such as the gain variance over time or frequency, local sidereal time (LST), diffractive scale, and phase structure-function slope with the level of excess variance. Our analysis shows that excess variance has an LST dependence, which is related to the power from the sky. And the simulated Stokes I power spectra from bright sources and the excess variance show a similar progression over LST with the minimum power appearing at LST bin 6h to 9h. This LST dependence is also present in sky images of the residual Stokes I of the observations. In very-wide sky images, we demonstrate that the extra power comes exactly from the direction of bright and distant sources Cassiopeia A and Cygnus A with the array beam patterns. These results suggest that the level of excess variance in the 21 cm signal power spectra is related to sky effects and, hence, it depends on LST. In particular, very bright and distant sources such as Cassiopeia A and Cygnus A can dominate the effect. This is in line with earlier studies and offers a path forward toward a solution since the correlation between the sky-related effects and the excess variance is non-negligible., Comment: 29 pages, 14 figures, accepted for publication in A&A

In this data release from the LOFAR Two-metre Sky Survey (LoTSS) we present 120-168MHz images covering 27% of the northern sky. Our coverage is split into two regions centred at approximately 12h45m +44$^\circ$30' and 1h00m +28$^\circ$00' and spanning 4178 and 1457 square degrees respectively. The images were derived from 3,451hrs (7.6PB) 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 144MHz have: a median rms sensitivity of 83$\mu$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.8mJy/beam. By creating three 16MHz 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-168MHz continuum images have a median rms sensitivity of 95$\mu$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 x 97.6 kHz wide planes and have a median rms sensitivity per plane of 10.8mJy/beam at 4' and 2.2mJy/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., Comment: 23 figures, 1 table and 29 pages. The catalogues, images and uv-data associated with this data release are publicly available via https://lofar-surveys.org/

High quality low-frequency radio surveys have the promise of advancing our understanding of many important topics in astrophysics, including the life cycle of active galactic nuclei (AGN), particle acceleration processes in jets, the history of star formation, and exoplanet magnetospheres. Currently leading low-frequency surveys reach an angular resolution of a few arcseconds. However, this resolution is not yet sufficient to study the more compact and distant sources in detail. Sub-arcsecond resolution is therefore the next milestone in advancing these fields. The biggest challenge at low radio frequencies is the ionosphere. If not adequately corrected for, ionospheric seeing blurs the images to arcsecond or even arcminute scales. Additionally, the required image size to map the degree-scale field of view of low-frequency radio telescopes at this resolution is far greater than what typical soft- and hardware is currently capable of handling. Here we present for the first time (to the best of our knowledge) widefield sub-arcsecond imaging at low radio frequencies. We derive ionospheric corrections in a few dozen individual directions and apply those during imaging efficiently using a recently developed imaging algorithm (arXiv:1407.1943, arXiv:1909.07226). We demonstrate our method by applying it to an eight hour observation of the International LOw Frequency ARray (LOFAR) Telescope (ILT) (arXiv:1305.3550). Doing so we have made a sensitive $7.4\ \mathrm{deg}^2$ $144\ \mathrm{MHz}$ map at a resolution of $0.3''$ reaching $25\ \mu\mathrm{Jy\ beam}^{-1}$ near the phase centre. The estimated $250,000$ core hours used to produce this image, fit comfortably in the budget of available computing facilities. This result will enable future mapping of the entire northern low-frequency sky at sub-arcsecond resolution., Comment: 19 pages, 4 figures; This preprint has not undergone peer review or any post-submission improvements or corrections. The Version of Record of this article is published in Nature Astronomy, and is available online at https://doi.org/10.1038/s41550-021-01573-z

Apertif (APERture Tile In Focus) is one of the Square Kilometre Array (SKA) pathfinder facilities. The Apertif project is an upgrade to the 50-year-old Westerbork Synthesis Radio Telescope (WSRT) using phased-array feed technology. The new receivers create 40 individual beams on the sky, achieving an instantaneous sky coverage of 6.5 square degrees. The primary goal of the Apertif Imaging Survey is to perform a wide survey of 3500 square degrees (AWES) and a medium deep survey of 350 square degrees (AMES) of neutral atomic hydrogen (up to a redshift of 0.26), radio continuum emission and polarisation. Each survey pointing yields 4.6 TB of correlated data. The goal of Apercal is to process this data and fully automatically generate science ready data products for the astronomical community while keeping up with the survey observations. We make use of common astronomical software packages in combination with Python based routines and parallelisation. We use an object oriented module-based approach to ensure easy adaptation of the pipeline. A Jupyter notebook based framework allows user interaction and execution of individual modules as well as a full automatic processing of a complete survey observation. If nothing interrupts processing, we are able to reduce a single pointing survey observation on our five node cluster with 24 physical cores and 256 GB of memory each within 24h keeping up with the speed of the surveys. The quality of the generated images is sufficient for scientific usage for 44 % of the recorded data products with single images reaching dynamic ranges of several thousands. Future improvements will increase this percentage to over 80 %. Our design allowed development of the pipeline in parallel to the commissioning of the Apertif system.

Aims: Contamination from bright diffuse Galactic thermal and non-thermal radio emission poses crucial challenges in experiments aiming to measure the 21-cm signal of neutral hydrogen from the Cosmic Dawn and Epoch of Reionization. If not included in calibration, this diffuse emission can severely impact the analysis and signal extraction in 21-cm experiments. We examine large-scale diffuse Galactic emission at 122~MHz, around the North Celestial Pole, using the AARTFAAC-HBA system. Methods: In this pilot project, we present the first-ever wide-field image produced with a single sub-band of the data recorded with the AARTFAAC-HBA. We demonstrate two methods: multiscale CLEAN and shapelet decomposition, to model the diffuse emission revealed in the image. We use angular power spectrum metrics to quantify different components of the emission and compare the performance of the two diffuse structure modelling approaches. Results: We observe that the point sources dominate the angular power spectrum ($\ell(\ell+1)C_{\ell}/2\pi\equiv \Delta^2(\ell)$) of the emission in the field on scales $\ell\gtrsim 60$ ($\lesssim 3$~degree). The angular power spectrum after subtraction of compact sources is flat within $20\lesssim \ell \lesssim200$ range, suggesting that the residual power is dominated by the diffuse emission on scales $\ell\lesssim200$. The residual diffuse emission has a brightness temperature variance of $\Delta^2_{\ell=180} = (145.64 \pm 13.61)~{\rm K}^2$ at 122~MHz on angular scales of 1~degree, and is consistent with a power-law following $C_{\ell}\propto \ell^{-2.0}$ in $20\lesssim \ell \lesssim200$ range. We also find that, in the current setup, the multiscale CLEAN is suitable to model the compact and diffuse structures on a wide range of angular scales, whereas the shapelet decomposition method better models the large scales, which are of the order of a few degrees and wider., Comment: 11 pages, 9 figures, accepted (A&A)

Supermassive black holes and supernovae explosions at the centres of active galaxies power cycles of outflowing and inflowing gas that affect galactic evolution and the overall structure of the Universe. While simulations and observations show that this must be the case, the range of physical scales (over ten orders of magnitude) and paucity of available tracers, make both the simulation and observation of these effects difficult. By serendipity, there lies an active galaxy, Centaurus A (NGC 5128), at such a close proximity as to allow its observation over this entire range of scales and across the entire electromagnetic spectrum. In the radio band, however, details on scales of 10-100 kpc from the supermassive black hole have so far been obscured by instrumental limitations. Here we report low-frequency radio observations that overcome these limitations and show evidence for a broad, bipolar outflow with velocity 1100 km per s and mass outflow rate of 2.9 solar masses per year on these scales. We combine our data with the plethora of multi-scale, multi-wavelength historical observations of Centaurus A to probe a unified view of feeding and feedback, which we show to be consistent with the Chaotic Cold Accretion self-regulation scenario., Comment: 15 pages, 5 figures, 2 tables

We investigate systematic effects in direction dependent gain calibration in the context of the Low-Frequency Array (LOFAR) 21-cm Epoch of Reionization (EoR) experiment. The LOFAR EoR Key Science Project aims to detect the 21-cm signal of neutral hydrogen on interferometric baselines of $50-250 \lambda$. We show that suppression of faint signals can effectively be avoided by calibrating these short baselines using only the longer baselines. However, this approach causes an excess variance on the short baselines due to small gain errors induced by overfitting during calibration. We apply a regularised expectation-maximisation algorithm with consensus optimisation (sagecal-co) to real data with simulated signals to show that overfitting can be largely mitigated by penalising spectrally non-smooth gain solutions during calibration. This reduces the excess power with about a factor 4 in the simulations. Our results agree with earlier theoretical analysis of this bias-variance trade off and support the gain-calibration approach to the LOFAR 21-cm signal data., Comment: This is a pre-copyedited, author-produced PDF of an article accepted for publication in MNRAS following peer review

We present the first targeted measurement of the power spectrum of anisotropies of the radio synchrotron background, at 140 MHz where it is the overwhelmingly dominant photon background. This measurement is important for understanding the background level of radio sky brightness, which is dominated by steep-spectrum synchrotron radiation at frequencies below 0.5 GHz and has been measured to be significantly higher than that which can be produced by known classes of extragalactic sources and most models of Galactic halo emission. We determine the anisotropy power spectrum on scales ranging from 2 degrees to 0.2 arcminutes with LOFAR observations of two 18 square degree fields -- one centered on the Northern hemisphere coldest patch of radio sky where the Galactic contribution is smallest and one offset from that location by 15 degrees. We find that the anisotropy power is higher than that attributable to the distribution of point sources above 100 micro-Jy in flux. This level of radio anisotropy power indicates that if it results from point sources, those sources are likely at low fluxes and incredibly numerous, and likely clustered in a specific manner., Comment: 8 pages, 5 figures, published in MNRAS, updated to published version and typos removed

We describe the APERture Tile In Focus (Apertif) system, a phased array feed (PAF) upgrade of the Westerbork Synthesis Radio Telescope which has transformed this telescope into a high-sensitivity, wide field-of-view L-band imaging and transient survey instrument. Using novel PAF technology, up to 40 partially overlapping beams can be formed on the sky simultaneously, significantly increasing the survey speed of the telescope. With this upgraded instrument, an imaging survey covering an area of 2300 deg2 is being performed which will deliver both continuum and spectral line data sets, of which the first data has been publicly released. In addition, a time domain transient and pulsar survey covering 15,000 deg2 is in progress. An overview of the Apertif science drivers, hardware and software of the upgraded telescope is presented, along with its key performance characteristics., Comment: 29 pages, 42 figures, accepted for publication by A&A