Search

The clusters Abell 2061 and Abell 2067 in the Corona Borealis supercluster have been studied at different radio frequencies and are both known to host diffuse radio emission. The aim of this work is to investigate the radio emission in between them, suggested by low resolution observations. We analyse deep LOFAR HBA observations at 144 MHz to follow up on the possible intercluster filament suggested by previous 1.4 GHz observations. We investigate the radial profiles and the point-to-point surface brightness correlation of the emission in A2061 with radio and Xray observations, to describe the nature of the diffuse emission. We report the detection of diffuse radio emission on 800 kpc scale, more extended than previously known, reaching beyond the radio halo in A2061 towards A2067 and over the separation outside the two clusters R500 radii. We confirm the presence of a radio halo in A2061, while do not find evidence of diffuse emission in A2067. The surface brightness profile from the centre of A2061 shows an excess of emission with respect to the azimuthally averaged radio halo profile and X-ray background. We explore three different dynamical scenario to explain the nature of the diffuse emission. We analyse a trail of emission of 760 kpc between the radio halo and radio relic in A2061. This pre merger system closely resembles the two other cluster pairs where radio bridges connecting the radio halos on Mpc scales have been detected. The diffuse emission extends beyond each cluster R500 radius but in this unique case, the absence of the radio halo in A2067 is likely the reason for the observed 'gap' between the two systems. However, the point-to-point correlation results are challenging to explain. The classification of the emission remains unclear, and detailed spectral analysis and further Xray observations are required to understand the origin of the diffuse emission., Comment: Accepted for publication on Astronomy & Astrophysics 12 pages, 6 figures, 1 appendix

High-precision pulsar timing is highly dependent on precise and accurate modeling of any effects that impact the data. It was shown that commonly used Solar Wind models do not accurately account for variability in the amplitude of the Solar wind on both short and long time scales. In this study, we test and validate a new, cutting-edge Solar wind modeling method included in the \texttt{enterprise} software suite through extended simulations, and we apply it to investigate temporal variability in LOFAR data. Our model testing scheme in itself provides an invaluable asset for pulsar timing array (PTA) experiments. As improperly accounting for the solar wind signature in pulsar data can induce false-positive signals, it is of fundamental importance to include in any such investigations. We employ a Bayesian approach utilizing a continuously varying Gaussian process to model the solar wind referred to as Solar Wind Gaussian Process (SWGP). We conduct noise analysis on eight pulsars from the LOFAR dataset with most pulsars having a timespan of $\sim 11$ years encompassing one full solar activity cycle. Our analysis reveals a strong correlation between the electron density at 1 AU and the ecliptic latitude (ELAT) of the pulsar. Pulsars with $|ELAT|< 3^{\circ}$ exhibit significantly higher average electron densities. We observe distinct temporal patterns in electron densities in different pulsars. In particular, pulsars within $|ELAT|< 3^{\circ}$ exhibit similar temporal variations, while the electron densities of those outside this range correlate with the solar activity cycle. The continuous variability in electron density offered in this model represents a substantial improvement over previous models, which assume a single value for piece-wise bins of time. This advancement holds promise for solar wind modeling in future International Pulsar Timing Array data combinations., Comment: Accepted in Section 9. Sun and the Heliosphere, Astronomy and Astrophysics

In the development of any type of magnetic field sensor based on magnetic films, special consideration must be given to the magnetic layer component. The presented work investigates the use of flux closing magnetostrictive multilayers for inverse magnetoelectric sensors. In such a type of magnetic field sensor, highly sensitive AC and DC field detection relies on strong excitation of the incorporated magnetic layers by piezoelectrically driven cantilever oscillation at mechanical resonances. The provoked periodic flux change is influenced by the magnetic field to be measured and is picked up by a coil, which generates the measured output. The effect of the magnetic multilayer on linearity, noise behavior, and detection limit of DC and AC signals is investigated. This study demonstrates the next step for inverse magnetoelectric thin film sensors, which achieve one order of magnitude improved detection limits with less than $8 pT/Hz^{1/2}$ at $10 Hz$ and $18 pT/Hz^{1/2}$ at $DC$ using exchange bias stabilized magnetic multilayers for obtaining flux closure., Comment: 12 pages, 5 figures

The Next Generation Very Large Array (ngVLA) is a planned radio interferometer providing unprecedented sensitivity at wavelengths between 21 cm and 3 mm. Its 263 antenna element array will be spatially distributed across North America to enable both superb low surface brightness recovery and sub-milliarcsecond angular resolution imaging. The project was developed by the international astronomy community under the lead of the National Radio Astronomy Observatory (NRAO), and is anticipated to be built between 2027 and 2037. Two workshops have been held in 2022 and 2023 with the goal to discuss and consolidate the scientific interests in the ngVLA within the German astronomical community. This community paper constitutes a collection of 48 science ideas which the German community aims to pursue with the ngVLA in the 2030s. This is not a complete list and the ideas are not developed at the level of a "Science Book", such that the present document is mainly meant provide a basis for further discussion within the community. As such, additional contributions are welcome, and will be considered for inclusion in future revisions., Comment: Version 2.0 (status June 18, 2024): 169 pages, comments and future contributions welcome [v2.0: 7 new science cases added, some minor revisions to other chapters]

In this work, we investigated the interplay between the X-ray and radio emission of the cluster PSZ2G113.91-37.01 (z = 0.371) using the high-quality XMM-Newton observations of the CHEX-MATE project, and the images of the LoTSS-DR2. The cluster is undergoing a merger along the north-south axis, and shows a central radio halo and two radio relics, one in the southern and one in the northern regions. The analysis of the intracluster medium distribution revealed the presence of a northern surface brightness jump associated to the merger event. By extracting spectra across this discontinuity, we classified the edge as a cold front. Furthermore, we made use of upgraded Giant Metrewave Radio Telescope observations that allowed us to perform a spectral analysis of the G113 radio emission. We found evidence of re-acceleration of particles in the northern relic, and we measured an associated Mach number of M = 1.95 $\pm$ 0.01, as inferred from radio observations. We then performed a point-to-point analysis of the X-ray and radio emission both in the halo and in the northern relic regions. We found a strong correlation for the halo and an anti-correlation for the relic. The former behaviour is in agreement with previous studies. The relic anti-correlation is likely related to the reverse radial distribution of the X-ray (increasing towards the cluster centre) and radio (decreasing towards the cluster centre) emissions. Finally, we performed a point-to-point analysis of the radio emission and the residuals obtained by subtracting a double beta model to the X-ray emission. We found a strong correlation between the two quantities. This behaviour suggests the presence of a connection between the process responsible for the radio emission and the one that leaves fluctuations in the X-ray observations., Comment: 13 pages, 8 figures, 1 table. Accepted for publication in A&A

The second data release of the LOFAR Two-Metre Sky Survey (LoTSS) covers 27% of the northern sky, with a total area of $\sim 5,700$ deg$^2$. The high angular resolution of LOFAR with Dutch baselines (6 arcsec) allows us to carry out optical identifications of a large fraction of the detected radio sources without further radio followup; however, the process is made more challenging by the many extended radio sources found in LOFAR images as a result of its excellent sensitivity to extended structure. In this paper we present source associations and identifications for sources in the second data release based on optical and near-infrared data, using a combination of a likelihood-ratio cross-match method developed for our first data release, our citizen science project Radio Galaxy Zoo: LOFAR, and new approaches to algorithmic optical identification, together with extensive visual inspection by astronomers. We also present spectroscopic or photometric redshifts for a large fraction of the optical identifications. In total 4,116,934 radio sources lie in the area with good optical data, of which 85% have an optical or infrared identification and 58% have a good redshift estimate. We demonstrate the quality of the dataset by comparing it with earlier optically identified radio surveys. This is by far the largest ever optically identified radio catalogue, and will permit robust statistical studies of star-forming and radio-loud active galaxies., Comment: 29 pages. Accepted by A&A; data products available at https://lofar-surveys.org/dr2_release.html

We present the peculiar case of PSZ2G091.83+26.11 at z=0.822. This cluster hosts a Mpc-scale radio halo and an elongated radio source, whose location with the respect to the intracluster medium (ICM) distribution and to the cluster centre is not consistent with a simple merger scenario. We use VLA data at 1-4 GHz to investigate the spectral and polarisation properties of the diffuse radio emission. We combine them with previously published data from LOFAR n the 120-168 MHz band, and from the uGMRT at 250-500 and 550-900 MHz. We also complement the radio data with Chandra X-ray observations to compare the thermal and non-thermal emission of the cluster. The elongated radio emission is visible up to 3.0 GHz and has an integrated spectral index of $-1.24\pm0.03$, with a steepening from $-0.89\pm0.03$ to $-1.39\pm0.03$. These values correspond to Mach numbers $\mathcal{M}_{\rm radio,int}=3.0\pm0.19$ and $\mathcal{M}_{\rm radio,inj}=2.48\pm0.15$. Chandra data reveals a surface brightness discontinuity at the location of the radio source, with a compression factor of $\mathcal{C}=2.22^{+0.39}_{-0.30}$ (i.e. $\mathcal{M}_{\rm Xray}=1.93^{+0.42}_{-0.32}$). We also find that the source is polarised at GHz frequencies. We estimate an intrinsic polarisation fraction of $\sim0.2$, a Rotation Measure of $\sim50~{\rm rad~m^{-2}}$ (including the Galactic contribution) and an external depolarisation of $\sim60~{\rm rad~m^{-2}}$. The $B$-vectors are aligned with the major axis of the source, suggesting magnetic field compression. Hence, we classify this source as a radio relic. We also find a linear/super-linear correlation between the non-thermal and thermal emission. We propose an off-axis merger and/or multiple merger events to explain the position and orientation of the relic. Given the properties of the radio relic, we speculate that PSZ2G091.83+26.11 is in a fairly young merger state., Comment: 21 pages, 17 figures, 9 tables; accepted for publication in Astronomy & Astrophysics

We present the detection of 68 sources from the most sensitive radio survey in circular polarisation conducted to date. We use the second data release of the 144 MHz LOFAR Two-metre Sky Survey to produce circularly-polarised maps with median 140 $\mu$Jy beam$^{-1}$ noise and resolution of 20$''$ for $\approx$27% of the northern sky (5634 deg$^{2}$). The leakage of total intensity into circular polarisation is measured to be $\approx$0.06%, and our survey is complete at flux densities $\geq1$ mJy. A detection is considered reliable when the circularly-polarised fraction exceeds 1%. We find the population of circularly-polarised sources is composed of four distinct classes: stellar systems, pulsars, active galactic nuclei, and sources unidentified in the literature. The stellar systems can be further separated into chromospherically-active stars, M dwarfs, and brown dwarfs. Based on the circularly-polarised fraction and lack of an optical counterpart, we show it is possible to infer whether the unidentified sources are likely unknown pulsars or brown dwarfs. By the completion of this survey of the northern sky, we expect to detect 300$\pm$100 circularly-polarised sources., Comment: Accepted for publication in A&A, 15 pages, 8 figures, 1 table. The catalogue will be publicly available at http://lofar-surveys.org/ and via Vizier shortly

We present the first detailed analysis of the radio halo in the merging galaxy cluster Abell 2256 using the LOFAR, uGMRT, and VLA. These observations combined with archival X-ray data allowed us to study the halo emission with unprecedented detail. The integrated radio emission from the entire halo is characterized by an ultra-steep spectrum, which can be described by a power law with $\alpha_{144 \rm MHz}^{1.5 \rm GHz}=-1.63\pm0.03$, and a radial steepening in the outer regions. The halo is significantly underluminous according to the scaling relations between radio power and mass at 1.4 GHz but not at 150 MHz; ultra-steep spectrum halos are predicted to be statistically underluminous. Despite the complex structure of this system, the radio halo morphology is remarkably similar to that of the X-ray emission. The radio surface brightness distribution across the halo is strongly correlated with the X-ray brightness of the intracluster medium. The derived correlations show sublinear slopes and there are distinct structures: the core is $\rm I_{R}\propto I_{X}^{1.51}$, the outermost region $\rm I_{R}\propto I_{X}^{0.41}$, and we find radio morphological connections with X-ray discontinuities. We also find a strong anti-correlation between the radio spectral index and the X-ray surface brightness, implying radial steepening. We suggests that the halo core is either related to old plasma from previous AGN activity, being advected, compressed and re-accelerated by mechanisms activated by the cold front or less turbulent with strong magnetic field in the core. The change in the radio vs X-ray correlation slopes in the outer regions of the halo could be due to a radial decline of magnetic field, increase in the number density of seed particles or increasing turbulence. Our findings suggest that that the emitting volume is not homogenous according to turbulence re-acceleration models., Comment: 19 pages, 15 figures, accepted for publication in A&A

Abell 3266 is a massive and complex merging galaxy cluster that exhibits significant substructure. We present new, highly sensitive radio continuum observations of Abell 3266 performed with the Australian Square Kilometre Array Pathfinder (0.8$-$1.1 GHz) and the Australia Telescope Compact Array (1.1$-$3.1 GHz). These deep observations provide new insights into recently-reported diffuse non-thermal phenomena associated with the intracluster medium, including a 'wrong-way' relic, a fossil plasma source, and an as-yet unclassified central diffuse ridge, which we reveal comprises the brightest part of a large-scale radio halo detected here for the first time. The 'wrong-way' relic is highly atypical of its kind: it exhibits many classical signatures of a shock-related radio relic, while at the same time exhibiting strong spectral steepening. While radio relics are generally consistent with a quasi-stationary shock scenario, the 'wrong-way' relic is not. We study the spectral properties of the fossil plasma source; it exhibits an ultra-steep and highly curved radio spectrum, indicating an extremely aged electron population. The larger-scale radio halo fills much of the cluster centre, and presents a strong connection between the thermal and non-thermal components of the intracluster medium, along with evidence of substructure. Whether the central diffuse ridge is simply a brighter component of the halo, or a mini-halo, remains an open question. Finally, we study the morphological and spectral properties of the multiple complex radio galaxies in this cluster in unprecedented detail, tracing their evolutionary history., Comment: Main article contains 26 pages, 15 figures, 6 tables. Supplementary material contains figures illustrating the improvement gain by applying third-generation calibration and imaging techniques to ATCA and ASKAP data. Manuscript published Open Access in MNRAS on August 1st

The presence of large-scale magnetic fields and ultra-relativistic electrons in the intra-cluster medium (ICM) is confirmed through the detection of diffuse radio synchrotron sources, so-called radio halos and relics. Due to their steep-spectrum nature, these sources are rarely detected at frequencies above a few GHz, especially in low-mass systems. The aim of this study is to discover and characterise diffuse radio sources in low-mass galaxy clusters in order to understand their origin and their scaling with host cluster properties. We searched for cluster-scale radio emission from low-mass galaxy clusters in the Low Frequency Array (LOFAR) Two-metre Sky Survey - Data Release 2 (LoTSS-DR2) fields. We made use of existing optical (Abell, DESI, WHL) and X-ray (comPRASS, MCXC) catalogues. The LoTSS-DR2 data were processed further to improve the quality of the images that are used to detect and characterize diffuse sources. We have detected diffuse radio emission in 28 galaxy clusters. The number of confirmed (candidates) halos and relics are six (seven) and 10 (three), respectively. Among these, 11 halos and 10 relics, including candidates, are newly discovered by LOFAR. Beside these, five diffuse sources are detected in tailed radio galaxies and are probably associated with mergers during the formation of the host clusters. We are unable to classify other 13 diffuse sources. We compare our newly detected, diffuse sources to known sources by placing them on the scaling relation between the radio power and the mass of the host clusters.

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/

Relativistic electrons and magnetic fields permeate the intra-cluster medium (ICM) and manifest themselves as diffuse sources of synchrotron emission observable at radio wavelengths, namely radio halos and radio relics. Although there is broad consensus that the formation of these sources is connected to turbulence and shocks in the ICM, the details of the required particle acceleration, the strength and morphology of the magnetic field in the cluster volume, and the influence of other sources of high-energy particles are poorly known. Sufficiently large samples of radio halos and relics, which would allow us to examine the variation among the source population and pinpoint their commonalities and differences, are still missing. At present, large numbers of these sources are easiest to detect at low radio frequencies, where they shine brightly. We examined the low-frequency radio emission from all 309 clusters in the second catalog of Planck Sunyaev Zel'dovich detected sources that lie within the 5634 deg$^2$ covered by the Second Data Release of the LOFAR Two-meter Sky Survey (LoTSS-DR2). We produced LOFAR images at different resolutions, with and without discrete sources subtracted, and created overlays with optical and X-ray images before classifying the diffuse sources in the ICM, guided by a decision tree. Overall, we found 83 clusters that host a radio halo and 26 that host one or more radio relics (including candidates). About half of them are new discoveries. The detection rate of clusters hosting a radio halo and one or more relics in our sample is $30\pm11$% and $10\pm6$%, respectively. Extrapolating these numbers, we anticipate that once LoTSS covers the entire northern sky it will provide the detection of $251\pm92$ clusters hosting a halo and $83\pm50$ clusters hosting at least one relic from Planck clusters alone., Comment: Matched to the A&A published version. Abstract abridged to meet arXiv requirements. All FITS images and tables produced in this work are publicly available on the project website https://lofar-surveys.org/planck_dr2.html

We present deep and high fidelity images of the merging galaxy cluster Abell 2256 at low frequencies, using the upgraded Giant Metrewave Radio Telescope (uGMRT) and LOw-Frequency ARray (LOFAR). This cluster hosts one of the most prominent known relics, with a remarkably spectacular network of filamentary substructures. The new uGMRT (300-850 MHz) and LOFAR (120-169 MHz) observations, combined with the archival Karl G. Jansky Very Large Array (VLA; 1-4 GHz) data, allowed us to carry out the first spatially resolved spectral analysis of the exceptional relic emission down to 6 arcsec resolution over a broad range of frequencies. Our new sensitive radio images confirm the presence of complex filaments of magnetized relativistic plasma also at low frequencies. We find that the integrated spectrum of the relic is consistent with a single power law, without any sign of spectral steepening, at least below 3 GHz. Unlike previous claims, the relic shows an integrated spectral index of $-1.07\pm0.02$ between 144 MHz and 3 GHz, which is consistent with the (quasi)stationary shock approximation. The spatially resolved spectral analysis suggests that the relic surface very likely traces the complex shock front, with a broad distribution of Mach numbers propagating through a turbulent and dynamically active intracluster medium. Our results show that the northern part of the relic is seen edge-on and the southern part close to face-on. We suggest that the complex filaments are regions where higher Mach numbers dominate the (re-)acceleration of electrons that are responsible for the observed radio emission., Comment: 27 pages, 14 Figures, accepted for publication in ApJ

We present wideband polarimetric observations, obtained with the Karl G. Jansky Very Large Array (VLA), of the merging galaxy cluster MACS J0717.5+3745, which hosts one of the most complex known radio relic and halo systems. We use both Rotation Measure Synthesis and QU-fitting, and find a reasonable agreement of the results obtained with these methods, in particular, when the Faraday distribution is simple and the depolarization is mild. The relic is highly polarized over its entire length reaching a fractional polarization ${>}30\%$ in some regions. We also observe a strong wavelength-dependent depolarization for some regions of the relic. The northern part of the relic shows a complex Faraday distribution suggesting that this region is located in or behind the intracluster medium (ICM). Conversely, the southern part of the relic shows a Rotation Measure very close to the Galactic foreground, with a rather low Faraday dispersion, indicating very little magnetoionic material intervening the line-of-sight. From spatially resolved polarization analysis, we find that the scatter of Faraday depths correlates with the depolarization, indicating that the tangled magnetic field in the ICM causes the depolarization. At the position of a well known narrow-angle-tailed galaxy (NAT), we find evidence of two components clearly separated in Faraday space. The high Faraday dispersion component seems to be associated with the NAT, suggesting the NAT is embedded in the ICM while the southern part of the relic lies in front of it. The magnetic field orientation follows the relic structure indicating a well-ordered magnetic field. We also detect polarized emission in the halo region; however the absence of significant Faraday rotation and a low value of Faraday dispersion suggests the polarized emission, previously considered as the part of the halo, has a shock(s) origin., Comment: 22 pages, 19 figures, accepted for publication in A&A

Context. Radio halos are megaparsec-scale diffuse radio sources{ mostly} located at the centres of merging galaxy clusters. The common mechanism invoked to explain their origin is the re-acceleration of relativistic particles caused by large-scale turbulence. Aims. Current re-acceleration models predict that a significant number of halos at high redshift should be characterised by very steep spectra ($\alpha<-1.5$) because of increasing inverse Compton energy losses. In this paper, we investigate the spectral index properties of a sample of nine clusters selected from the second Planck Sunyaev-Zel'dovich catalogue showing diffuse radio emission with the Low Frequency Array (LOFAR) in the 120-168 MHz band. This is the first time that radio halos discovered at low frequencies are followed up at higher frequencies. Methods. We analysed upgraded Giant Metrewave Radio Telescope (uGMRT) observations in Bands 3 and 4, that is, 250-500 and 550-900 MHz respectively. These observations were combined with existing LOFAR data to obtain information on the spectral properties of the diffuse radio emission. Results. We find diffuse radio emission in the uGMRT observations for five of the nine high-$z$ radio halos previously discovered with LOFAR. For those, we measure spectral indices in the range of $-1$ to $-1.4$. For the uGMRT non-detections, we estimated that the halos should have a spectral index steeper than $-1.5$. We also confirm the presence of one candidate relic. Conclusions. Despite the small number of clusters, we find evidence that about half of the massive and merging clusters at high redshift host radio halos with a very steep spectrum. This is in line with theoretical predictions, although larger statistical samples are necessary to test models., Comment: 23 pages, 22 figures, 4 tables; accepted for publication in A&A

[abridged] 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, although this is technically and logistically challenging. Here we present a calibration strategy that builds on previous high-resolution work with LOFAR. 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 LOFAR Two-metre Sky Survey (LoTSS) pointing. We perform in-field delay calibration, solution referencing to other calibrators, self-calibration, and imaging of example directions of interest in the field. 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 demonstrate a check of the astrometry and flux density scale. Imaging in 17 directions, the restoring beam is typically 0.3" x 0.2" although this varies slightly over the entire 5 square degree field of view. We achieve ~80 to 300 $\mu$Jy/bm image rms noise, which is dependent on the distance from the phase centre; typical values are ~90 $\mu$Jy/bm 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 ~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., Comment: Accepted to a special issue of A&A on sub-arcsecond imaging with LOFAR. 24 pages, 16 figures

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 $\geq$ 200 km baselines is identified by the detection of compact flux density, for declinations north of 30 degrees 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. 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 arcsec, for sources observed in the survey. 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 arcsec 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., Comment: Accepted to a special issue of A&A on sub-arcsecond imaging with LOFAR

Radio relics are diffuse radio sources in galaxy clusters that are associated with merger shock waves. Detailed observations of radio relics in total intensity and in polarisation show complex structures on kiloparsec scales. The relation between the observed features and the underlying morphology of the magnetic field is not clear. Using three-dimensional magneto-hydrodynamical simulations, we study the polarised emission produced by a shock wave that propagates through a turbulent medium that resembles the intracluster medium. We model the polarised synchrotron emission on the basis of diffusive shock acceleration of cosmic-ray electrons. We find that the synchrotron emission produced in a shocked turbulent medium can reproduce some of the observed features in radio relics. Shock compression can give rise to a high polarisation fraction at the shock front and a partial alignment of the polarisation $E$-vectors with the shock normal. Our work confirms that radio relics can also be formed in an environment with a tangled magnetic field. We also discuss the effect of Faraday Rotation intrinsic to the source, and how our results depend on the angular resolution of observations., Comment: Accepted for publication in MNRAS, 22 pages, 23 figures

Mega-parsec scale radio sources in the form of halos and relics are often detected in dynamically disturbed galaxy clusters. Although they are associated with merger-induced turbulence and shocks, respectively, their formation is not fully understood. We aim to identify the mechanisms responsible for particle acceleration and magnetic field amplification in the halo and relics of the galaxy cluster ClG 0217+70. We observed ClG 0217+70 with LOFAR at 141 MHz and VLA at 1.5 GHz, and combine these observations with VLA 1.4 GHz archival data to study the morphological and spectral properties of the diffuse sources. We add Chandra archival data to examine the thermal and non-thermal properties of the halo. Our LOFAR and VLA data confirm the presence of a giant radio halo in the cluster centre and multiple relics in the outskirts. The radio and X-ray emission from the halo are correlated, implying a tight relation between the thermal and non-thermal components. The diffuse radio structure in the south east with a projected size of 3.5 Mpc is the most extended radio relic detected to date. The spectral index across the relic width steepens towards the cluster centre, suggesting the electron ageing in the post-shock regions. The shock Mach numbers for the relics derived from the spectral index map range between 2.0 and 3.2. However, the integrated spectral indices lead to increasingly high Mach numbers for the relics farther from the cluster centre. This discrepancy could be because the relation between injection and integrated spectra does not hold for distant shocks, suggesting that the cooling time for the radio-emitting electrons is longer than the crossing time of the shocks. The variations in the surface brightness of the relics and the low Mach numbers imply that the radio-emitting electrons are re-accelerated from fossil plasma that could originate in active galactic nuclei.

The steep spectrum radio source VLSSJ2217.5+5943 shows a complex, filamentary morphology and a curved spectrum. Therefore, the source has previously been classified as a radio phoenix. However, no galaxy cluster associated with this radio source has been confidently detected so far because the source is located in the direction of the innermost zone of the Galactic Plane at b = +2.4 degr (innermost Zone of Avoidance, ZoA). We analysed archival observations in the near infrared (UKIDSS) and mid infrared (Spitzer) to select the galaxies in the immediate neighbourhood of the radio source. A sample of 23 galaxies was selected as candidate cluster members. Furthermore, we carried out deep integral field spectroscopy covering 6450 to 10500 AA with the red unit of the Hobby-Eberly Telescope second generation low resolution spectrograph (LRS2-R). We also reanalysed archival GMRT observations at 325 and 610 MHz. We selected 23 galaxies within a radius of 2.5 arcmin, centered on RA=22:17.5, DEC=+59:43 (J2000). Spectra were obtained for three of the brightest galaxies. For two galaxies we derived redshifts of z = 0.165 and z = 0.161, based on NaD absorption and TiO band heads. Their spectra correspond to E-type galaxies. Both galaxies are spatially associated with VLSSJ2217.5+5943. The spectrum of the third galaxy, which is slightly more distant from the radio source, indicates a LINER at z = 0.042. It is apparently a foreground galaxy with respect to the cluster we identified. VLSSJ2217.5+5943 is associated with a massive galaxy cluster at redshift z = 0.163 +- .003, supporting its classification as radio phoenix., Comment: 12 pages, 11 figures, Astronomy & Astrophysics in press

We present the first deep low frequency radio observations of the massive and highly disturbed galaxy cluster Abell 2744 using the upgraded Giant Metrewave Radio Telescope (uGMRT). The cluster is experiencing a very complex multiple merger and hosts a giant halo and four radio relics. The uGMRT observations, together with existing VLA and Chandra observations, allow us to study the complexity of the physical mechanisms active in this system. Our new images reveal that the central halo emission is more extended toward low frequencies. We find that the integrated spectrum of the halo follows a power-law between 150 MHz and 3 GHz, while its subregions show significantly different spectra, also featuring high frequency spectral steepening. The halo also shows local regions in which the spectral index is significantly different from the average value. Our results highlight that an overall power-law spectrum, as observed in many radio halos, may also arise from the superposition of different subcomponents. The comparison of the radio surface brightness and spectral index with the X-ray brightness and temperature reveals for the first time different trends, indicating that the halo consists of two main components. All four relics in this system follow a power-law radio spectrum, compatible with shocks with Mach numbers in the range $3.0-4.5$. All relics are also highly polarized from 1-4 GHz and show low Faraday dispersion measures, suggesting that they are located in the outermost regions of the cluster. The complexity in the distribution and properties of nonthermal components in Abell 2744 supports a multiple merger scenario, as also highlighted by previous X-ray and lensing studies. Our unique results demonstrate the importance of sensitive and high-resolution, multi-frequency radio observations for understanding the interplay between the thermal and non-thermal components of the ICM., Comment: 24 pages, 19 figures, accepted for publication in A&A

The low-frequency radio spectra of the hotspots within powerful radio galaxies can provide valuable information about the physical processes operating at the site of the jet termination. These processes are responsible for the dissipation of jet kinetic energy, particle acceleration, and magnetic-field generation. Here we report new observations of the powerful radio galaxy Cygnus A using the Low Frequency Array (LOFAR) between 109 and 183 MHz, at an angular resolution of ~3.5 arcsec. The radio emission of the lobes is found to have a complex spectral index distribution, with a spectral steepening found towards the centre of the source. For the first time, a turnover in the radio spectrum of the two main hotspots of Cygnus A has been directly observed. By combining our LOFAR imaging with data from the Very Large Array at higher frequencies, we show that the very rapid turnover in the hotspot spectra cannot be explained by a low-energy cut-off in the electron energy distribution, as has been previously suggested. Thermal (free-free) absorption or synchrotron self absorption models are able to describe the low-frequency spectral shape of the hotspots, however, as with previous studies, we find that the implied model parameters are unlikely, and interpreting the spectra of the hotspots remains problematic., Comment: 9 pages, 3 figures, accepted for publication in MNRAS

Thermal gas in the center of galaxy clusters can show substantial motions that generate surface-brightness and temperature discontinuities known as cold fronts. The motions may be triggered by minor or off-axis mergers that preserve the cool core of the system. The dynamics of the thermal gas can also generate radio emission from the intra-cluster medium (ICM) and impact the evolution of clusters' radio sources. We aim to study the central region of Abell 1775, a system in an ambiguous dynamical state at $z=0.072$ which is known to host an extended head-tail radio galaxy, with the goal of investigating the connection between thermal and nonthermal components in its center. We made use of a deep (100 ks) Chandra observation accompanied by LOFAR 144 MHz, GMRT 235 MHz and 610 MHz, and VLA 1.4 GHz radio data. We find a spiral-like pattern in the X-ray surface brightness that is mirrored in the temperature and pseudo-entropy maps. Additionally, we characterize an arc-shaped cold front in the ICM. We interpret these features in the context of a slingshot gas tail scenario. The structure of the head-tail radio galaxy "breaks" at the position of the cold front, showing an extension that is detected only at low frequencies, likely due to its steep and curved spectrum. We speculate that particle reacceleration is occurring in the outer region of this tail, which in total covers a projected size of $\sim800$ kpc. We also report the discovery of revived fossil plasma with ultra-steep spectrum radio emission in the cluster core together with a central diffuse radio source that is bounded by the arc-shaped cold front. The results reported in this work demonstrate the interplay between thermal and nonthermal components in the cluster center and the presence of ongoing particle reacceleration in the ICM on different scales., Comment: 17 pages, 17 figures, 2 tables (including appendix); Updated to match the accepted version in A&A

Radio relics are elongated sources related to shocks driven by galaxy cluster merger events. Although these objects are highly polarized at GHz frequencies ($\gtrsim 20\%$), high-resolution studies of their polarization properties are still lacking. We present the first high-resolution and high-sensitivity polarimetry study of the merging galaxy cluster CIZA J2242.8+5301 in the 1-4 GHz frequency band. We use the $QU$-fitting approach to model the Stokes $I$, $Q$ and $U$ emission, obtaining best-fit intrinsic polarization fraction ($p_0$), intrinsic polarization angle ($\chi_0$), Rotation Measure (RM) and wavelength-dependent depolarization ($\sigma_{\rm RM}$) maps of the cluster. Our analysis focuses on the northern relic (RN). For the first time in a radio relic, we observe a decreasing polarization fraction in the downstream region. Our findings are possibly explained by geometrical projections and/or by decreasing of the magnetic field anisotropy towards the cluster center. From the amount of depolarization of the only detected background radio galaxy, we estimate a turbulent magnetic field strength of $B_{\rm turb}\sim5.6~\mu$Gauss in the relic. Finally, we observe Rotation Measure fluctuations of about 30 rad m$^{-2}$ around at the median value of 140.8 rad m$^{-2}$ at the relic position., Comment: 26 pages, 23 Figures, 5 Tables. Accepted for publication in ApJ

We present new LOFAR observations of the massive merging galaxy cluster MACS J0717.5+3745. The cluster hosts the most powerful radio halo known to date. These new observations, in combination with published uGMRT (300$-$850 MHz) and VLA (1$-$6.5 GHz) data, reveal that the halo is more extended than previously thought, with a largest linear size of $\sim2.2 \rm Mpc$. The halo shows a steep spectrum ($\alpha_{144\,\text{MHz}}^{1.5\,\text{GHz}}\sim-1.4$) and a steepening ($\alpha_{1.5 \text{GHz}}^{5.5 \text{GHz}}\sim-1.9$) above 1.5 GHz. We find a strong scattering in spectral index maps on scales of 50$-$100 kpc. We suggest that such a strong scattering may be a consequence of the regime where inverse Compton dominate the energy losses of electrons. The spectral index becomes steeper and shows an increased curvature in the outermost regions of the halo. We combined the radio data with \textit{Chandra} observations to investigate the connection between the thermal and non-thermal components of the intracluster medium (ICM). Despite a significant substructure in the halo emission, the radio brightness correlates strongly with the X-ray brightness at all observed frequencies. The radio-versus-X-ray brightness correlation slope steepens at a higher radio frequency (from $b_{144 \text{MHz}}=0.67\pm0.05$ to $b_{3.0 \text{GHz}}=0.98\pm0.09$) and the spectral index shows a significant anti correlation with the X-ray brightness. Both pieces of evidence further support a spectral steepening in the external regions. The compelling evidence for a steep spectral index, the existence of a spectral break above 1.5 GHz, and the dependence of radio and X-ray surface brightness correlation on frequency are interpreted in the context of turbulent reacceleration models. Under this scenario, our results allowed us to constrain that the turbulent kinetic pressure of the ICM is up to 10%., Comment: 16 pages, 12 figures, accepted for publication in A&A

High-precision pulsar timing requires accurate corrections for dispersive delays of radio waves, parametrized by the dispersion measure (DM), particularly if these delays are variable in time. In a previous paper we studied the Solar-wind (SW) models used in pulsar timing to mitigate the excess of DM annually induced by the SW, and found these to be insufficient for high-precision pulsar timing. Here we analyze additional pulsar datasets to further investigate which aspects of the SW models currently used in pulsar timing can be readily improved, and at what levels of timing precision SW mitigation is possible. Our goals are to verify: a) whether the data are better described by a spherical model of the SW with a time-variable amplitude rather than a time-invariant one as suggested in literature, b) whether a temporal trend of such a model's amplitudes can be detected. We use the pulsar-timing technique on low-frequency pulsar observations to estimate the DM and quantify how this value changes as the Earth moves around the Sun. Specifically, we monitor the DM in weekly to monthly observations of 14 pulsars taken with LOFAR across time spans of up to 6 years. We develop an informed algorithm to separate the interstellar variations in DM from those caused by the SW and demonstrate the functionality of this algorithm with extensive simulations. Assuming a spherically symmetric model for the SW density, we derive the amplitude of this model for each year of observations. We show that a spherical model with time-variable amplitude models the observations better than a spherical model with constant amplitude, but that both approaches leave significant SW induced delays uncorrected in a number of pulsars in the sample. The amplitude of the spherical model is found to be variable in time, as opposed to what has been previously suggested., Comment: 14 pages, 9 figures. Accepted for publication in Astronomy and Astrophysics

Radio relics are diffuse, extended synchrotron sources that originate from shock fronts generated during cluster mergers. The massive merging galaxy cluster MACS J0717.5+3745 hosts one of the more complex relics known to date. We present upgraded Giant Metrewave Radio Telescope band 3 (300-500 MHz) and band 4 (550-850 MHz) observations. These new observations, combined with published VLA and the new LOFAR HBA data, allow us to carry out a detailed, high spatial resolution spectral analysis of the relic over a broad range of frequencies. The integrated spectrum of the relic closely follows a power-law between 144 MHz and 5.5 GHz with a mean spectral slope $\alpha=-1.16\pm0.03$. Despite its complex morphology, the subregions of the relic and the other isolated filaments also follow power-law behaviors, and show similar spectral slopes. Assuming Diffusive Shock Acceleration, we estimate a dominant Mach number of $\sim 3.7$ for the shocks that make up the relic. Comparison with recent numerical simulations suggests that in the case of radio relics, the slopes of the integrated radio spectra are determined by the Mach number of the accelerating shock, with $\alpha$ nearly constant, namely between $-1.13$ and $-1.17$, for Mach numbers $3.5 - 4.0$. The spectral shapes inferred from spatially resolved regions show curvature, we speculate that the relic is inclined along the line-of-sight. The locus of points in the simulated color-color plots changes significantly with the relic viewing angle. We conclude that projection effects and inhomogeneities in the shock Mach number dominate the observed spectral properties of the relic in this complex system. Based on the new observations we raise the possibility that the relic and a narrow-angle-tailed radio galaxy are two different structures projected along the same line-of-sight., Comment: 14 pages, 12 figures

Radio pulses from pulsars are affected by plasma dispersion, which results in a frequency-dependent propagation delay. Variations in the magnitude of this effect lead to an additional source of red noise in pulsar timing experiments, including pulsar timing arrays that aim to detect nanohertz gravitational waves. We aim to quantify the time-variable dispersion with much improved precision and characterise the spectrum of these variations. We use the pulsar timing technique to obtain highly precise dispersion measure (DM) time series. Our dataset consists of observations of 36 millisecond pulsars, which were observed for up to 7.1 years with the LOFAR telescope at a centre frequency of ~150 MHz. Seventeen of these sources were observed with a weekly cadence, while the rest were observed at monthly cadence. We achieve a median DM precision of the order of 10^-5 cm^-3 pc for a significant fraction of our sources. We detect significant variations of the DM in all pulsars with a median DM uncertainty of less than 2x10^-4 cm^-3 pc. The noise contribution to pulsar timing experiments at higher frequencies is calculated to be at a level of 0.1-10 us at 1.4 GHz over a timespan of a few years, which is in many cases larger than the typical timing precision of 1 us or better that PTAs aim for. We found no evidence for a dependence of DM on radio frequency for any of the sources in our sample. The DM time series we obtained using LOFAR could in principle be used to correct higher-frequency data for the variations of the dispersive delay. However, there is currently the practical restriction that pulsars tend to provide either highly precise times of arrival (ToAs) at 1.4 GHz or a high DM precision at low frequencies, but not both, due to spectral properties. Combining the higher-frequency ToAs with those from LOFAR to measure the infinite-frequency ToA and DM would improve the result., Comment: 16 pages, 8 figures, 2 tables

Diffuse cluster radio sources, in the form of radio halos and relics, reveal the presence of cosmic rays and magnetic fields in the intracluster medium (ICM). These cosmic rays are thought to be (re-)accelerated through ICM turbulence and shock waves generated by cluster merger events. Here we characterize the presence of diffuse radio emission in known galaxy clusters in the HETDEX Spring Field, covering 424 deg$^2$. For this, we developed a method to extract individual targets from LOFAR observations processed with the LoTSS DDF-pipeline. This procedure enables improved calibration and joint imaging and deconvolution of multiple pointings of selected targets. The calibration strategy can also be used for LOFAR Low-Band Antenna (LBA) and international-baseline observations. The fraction of Planck PSZ2 clusters with any diffuse radio emission apparently associated with the ICM is $73\pm17\%$. We detect a total of 10 radio halos and 12 candidate halos in the HETDEX Spring Field. Five clusters host radio relics. The fraction of radio halos in Planck PSZ2 clusters is $31\pm11\%$, and $62\pm15\%$ when including the candidate radio halos. Based on these numbers, we expect that there will be at least $183 \pm 65$ radio halos found in the LoTSS survey in PSZ2 clusters, in agreement with predictions. The integrated flux densities for the radio halos were computed by fitting exponential models to the radio images. From these flux densities, we determine the cluster mass (M$_{500}$) and Compton Y parameter (Y$_{500}$) 150 MHz radio power (P$_{\rm{150 MHz}}$) scaling relations for Planck PSZ2-detected radio halos. We find that the slopes of these relations are steeper than those determined from the 1.4 GHz radio powers. However, considering the uncertainties this is not a statistically significant result., Comment: submitted, 35 pages and 45 figures

Diffuse radio emission has been found in many galaxy clusters, predominantly in massive systems which are in the state of merging. The radio emission can usually be classified as relic or halo emission, which are believed to be related to merger shocks or volume-filling turbulence, respectively. Recent observations have revealed radio bridges for some pairs of very closeby galaxy clusters. The mechanisms that may allow to explain the high specific density of relativistic electrons, necessary to explain the radio luminosity of these bridge regions, are poorly explored. We analyse the galaxy cluster Abell 1430 with LoTSS data in detail and complement it with recent JVLA L-band observations, XMM-Newton, Chandra, and SDSS data. Moreover, we compare our results to clusters extracted from the "The Three Hundred Project" cosmological simulation. We find that Abell 1430 consists of two components, namely A1430-A and A1430-B. We speculate that the two components undergo an off-axis merger. The more massive component shows diffuse radio emission which can be classified as radio halo showing a low radio power given the mass of the cluster. Most interestingly, there is extended diffuse radio emission, dubbed as the `Pillow', which is apparently related to A1430-B and thus related to low density intracluster or intergalactic medium. To date, a only few examples for emission originating from such regions are known. These discoveries are crucial to constrain possible acceleration mechanisms, which may allow to explain the presence of relativistic electrons in these regions. Our results indicate a spectral index of $\alpha_{144\,\text{MHz}}^{1.5\,\text{GHz}}=-1.4\pm0.5$ for the Pillow. If future observations confirm a slope as flat as the central value of -1.4 or even flatter, this would pose a severe challenge for the electron acceleration scenarios., Comment: 14 pages, 12 figures, accepted for publication in A&A

Radio relics trace shock fronts generated in the intracluster medium (ICM) during cluster mergers. The particle acceleration mechanism at the shock fronts is not yet completely understood. We observed the Toothbrush relic with the Effelsberg and Sardinia Radio Telescope at 14.25 GHz and 18.6 GHz, respectively. Unlike previously claimed, the integrated spectrum of the relic closely follows a power law over almost three orders of magnitude in frequency, with a spectral index of $\alpha_{\rm 58\,MHz}^{\rm 18.6\,GHz}=-1.16\pm0.03$. Our finding is consistent with a power-law injection spectrum, as predicted by diffusive shock acceleration theory. The result suggests that there is only little magnetic field strength evolution downstream to the shock. From the lack of spectral steepening, we find that either the Sunyaev-Zeldovich decrement produced by the pressure jump is less extended than $\sim$ 600\,kpc along the line of sight or, conversely, that the relic is located far behind in the cluster. For the first time, we detect linearly polarized emission from the "brush" at 18.6 GHz. Compared to 8.3 GHz, the degree of polarization across the brush increases at 18.6 GHz, suggesting a strong Faraday depolarization towards lower frequencies. The observed depolarization is consistent with an intervening magnetized screen that arise from the dense ICM containing turbulent magnetic fields. The depolarization, corresponding to a standard deviation of the Rotation Measures as high as $\sigma_{\rm RM}= 212\pm23\rm \,rad\,m^{-2}$, suggests that the brush is located in or behind the ICM. Our findings indicate that the Toothbrush can be consistently explained by the standard scenario for relic formation, Comment: Accepted for publication in A&A Letters

We investigate the galaxy population in a field of the Perseus cluster that roughly covers the virial radius of the cluster. The galaxies were selected on Schmidt CCD images in B and H alpha in combination with SDSS images. We present a catalogue of 1294 galaxies. Morphological information was obtained for 90% of the galaxies from the `eyeball' inspection, partly supported by the surface brightness profile analysis. Redshifts were taken from SDSS, literature sources, and own spectroscopic observations and are available for 24% of the catalogues galaxies. The galaxy catalogue is used to derive cluster properties, such as radial profiles, indications of sub-structure, virial mass, and viral radius and to study the cluster galaxy population with regard to morphological types and peculiarities, star formation rates and active galactic nuclei. In addition to the statistical approach, we present brief individual descriptions of 18 cluster galaxies with conspicuous morphological peculiarities. (Abstract modified to match the arXiv format.), Comment: Accepted for publication in A&A

We present LOFAR observations of one of the most spectacular objects in the radio sky: Abell 2255. This is a nearby ($z = 0.0806$) merging galaxy cluster hosting one of the first radio halos ever detected in the intra-cluster medium (ICM). The deep LOFAR images at 144 MHz of the central $\sim10$ Mpc$^2$ region show a plethora of emission on different scales, from tens of kpc to above Mpc sizes. In this work, we focus on the innermost region of the cluster. Among the numerous interesting features observed, we discover remarkable bright and filamentary structures embedded in the radio halo. We incorporate archival WSRT 1.2 GHz data to study the spectral properties of the diffuse synchrotron emission and find a very complex spectral index distribution in the halo spanning a wide range of values. We combine the radio data with Chandra observations to investigate the connection between the thermal and non-thermal components by quantitatively comparing the radio and X-ray surface brightness and the spectral index of the radio emission with the thermodynamical quantities of the ICM. Despite the multitude of structures observed in the radio halo, we find that the X-ray and radio emission are overall well correlated. The fact that the steepest spectrum emission is located in the cluster center and traces regions with high entropy possibly suggests the presence of seed particles injected by radio galaxies that are spread in the ICM by turbulence generating the extended radio halo., Comment: 17 pages, 13 figures, 4 tables (including appendix). Accepted for publication in ApJ

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., Comment: 9 pages, 2 figures, accepted for publication in MNRAS

We report on the multi-frequency timing observations of 21 pulsars discovered in the LOFAR Tied-Array All-Sky Survey (LOTAAS). The timing data were taken at central frequencies of 149 MHz (LOFAR) as well as 334 and 1532 MHz (Lovell Telecope). The sample of pulsars includes 20 isolated pulsars and the first binary pulsar discovered by the survey, PSR J1658$+$3630. We modelled the timing properties of the pulsars, which showed that they have, on average, larger characteristic ages. We present the pulse profiles of the pulsars across the three observing bands, where PSR J1643$+$1338 showed profile evolution that appears not to be well-described by the radius-to-frequency-mapping model. Furthermore, we modelled the spectra of the pulsars across the same observing bands, using a simple power law, and found an average spectral index of $-1.9 \pm 0.5$. Amongst the pulsars studied here, PSR J1657$+$3304 showed large flux density variations of a factor of 10 over 300 days, as well as mode changing and nulling on timescales of a few minutes. We modelled the rotational and orbital properties of PSR J1658$+$3630, which has a spin period of 33 ms in a binary orbit of 3.0 days with a companion of minimum mass of 0.87$M_{\odot}$, likely a Carbon-Oxygen or Oxygen-Neon-Magnesium type white dwarf. PSR J1658$+$3630 has a dispersion measure of 3.0 pc cm$^{-3}$, making it possibly one of the closest binary pulsars known., Comment: 21 pages, 12 figures, Accepted for publication in the Monthly Notices of the Royal Astronomical Society

Radio relics are diffuse extended synchrotron sources that originate from shock fronts induced by galaxy cluster mergers. The actual particle acceleration mechanism at the shock fronts is still under debate. The galaxy cluster 1RXS J0603.3+4214 hosts one of the most intriguing examples of radio relics, known as the Toothbrush. We present new wideband radio continuum observations made with uGMRT and VLA. Our new observations, in combination with previously published data, allowed us to carry out a detailed high spatial resolution spectral and curvature analysis of the known diffuse radio emission sources, over a broad range of frequencies. The integrated spectrum of the Toothbrush follows closely a power law over close to 2 decades in frequency, with a spectral index of $-1.16\pm0.02$. We do not find any evidence of spectral steepening below 8 GHz. The subregions of the Toothbrush show an identical spectral slopes, suggesting that observed spectral index is rather set by the distribution of Mach numbers which may have a similar shape at different parts of the shock front. Indeed, numerical simulations show an intriguing similar spectral index, indicating that the radio spectrum is dominated by the average over the inhomogeneities within the shock, with most of the emission coming from the tail of the Mach number distribution. In contrast to the Toothbrush, the spectrum of the fainter relics show a high frequency steepening. The integrated spectrum of the halo follows a power law from 150 MHz to 3 GHz with a spectral index of $-1.16\pm0.04$. We do not find any evidence for spectral curvature, not even in subareas of the halo. This suggest a homogeneous acceleration throughout the cluster volume. Between the brush region of the Toothbrush and the halo, the color-color analysis revealed emission that was consistent with an overlap between the two different spectral regions., Comment: Accepted for publication in A&A

The new generation of high-resolution broad-band radio telescopes, like the Low Frequency Array (LOFAR), produces, depending on the level of compression, between 1 to 10 TB of data per hour after correlation. Such a large amount of scientific data demand powerful computing resources and efficient data handling strategies to be mastered. The LOFAR Two-metre Sky Survey (LoTSS) is a Key Science Project (KSP) of the LOFAR telescope. It aims to map the entire northern hemisphere at unprecedented sensitivity and resolution. The survey consist of 3 168 pointings, requiring about 30 PBytes of storage space. As a member of the German Long Wavelength Consortioum (GLOW) the Forschungszentrum J\"ulich (FSZ) stores in the Long Term Archive (LTA) about 50% of all LoTSS observations conducted to date. In collaboration with SURFsara in Amsterdam we developed service tools that enable the KSP to process LOFAR data stored in the J\"ulich LTA at the supercomputer JUWELS in an automated and robust fashion. Through our system more than 500 out of 800 existing LoTSS observations have already been processed with the prefactor pipeline. This pipeline calibrates the direction-independent instrumental and ionospheric effects and furthermore reduces the data size significantly. For continuum imaging, this processing pipeline is the standard pipeline that is executed before more advanced processing and image reconstruction methods are applied., Comment: 10 pages, 6 figures, NIC Symposium 2020, accepted

Radio relics are sites of electron (re)acceleration in merging galaxy clusters but the mechanism of acceleration and the topology of the magnetic field in and near relics are yet to be understood. We are carrying out an observational campaign on double relic galaxy clusters starting with RXC J1314.4-2515. With $Jansky Very Large Array$ multi-configuration observations in the frequency range 1-4 GHz, we perform both spectral and polarization analyses, using the Rotation Measure synthesis technique. We use archival $XMM-Newton$ observations to constrain the properties of the shocked region. We discover a possible connection between the activity of a radio galaxy and the emission of the eastern radio relic. In the northern elongated arc of the western radio relic, we detect polarized emission with an average polarization fraction of $31 \ \%$ at 3 GHz and we derive the Mach number of the underlying X-ray shock. Our observations reveal low levels of fractional polarization and Faraday-complex structures in the southern region of the relic, which point to the presence of thermal gas and filamentary magnetic field morphology inside the radio emitting volume. We measured largely different Rotation Measure dispersion from the two relics. Finally, we use cosmological magneto-hydrodynamical simulations to constrain the magnetic field, viewing angle, and to derive the acceleration efficiency of the shock. We find that the polarization properties of RXC J1314.4-2515 are consistent with a radio relic observed at $70^{\circ}$ with respect to the line of sight and that efficient re-acceleration of fossil electrons has taken place., Comment: 24 pages, 19 figures, 5 tables, accepted for publication in MNRAS

Galaxy clusters are the most massive gravitationally bound structures in the Universe. They grow by accreting smaller structures in a merging process that produces shocks and turbulence in the intra-cluster gas. We observed a ridge of radio emission connecting the merging galaxy clusters Abell 0399 and Abell 0401 with the Low Frequency Array (LOFAR) at 140 MHz. This emission requires a population of relativistic electrons and a magnetic field located in a filament between the two galaxy clusters. We performed simulations to show that a volume-filling distribution of weak shocks may re-accelerate a pre-existing population of relativistic particles, producing emission at radio wavelengths that illuminates the magnetic ridge., Comment: Published in Science

We present LOFAR data from 110--180~MHz of the merging galaxy cluster Abell 1682, alongside archival optical, radio and X-ray data. Our 6 arc-second resolution images at low frequencies reveal new structures associated with numerous radio galaxies in the cluster. At 20 arc-second resolution we see diffuse emission throughout the cluster over hundreds of kpc, indicating particle acceleration mechanisms are in play as a result of the cluster merger event and powerful active galactic nuclei. We show that a significant part of the cluster emission is from an old radio galaxy with very steep spectrum emission (having a spectral index of $\alpha < -2.5$). Furthermore we identify a new region of diffuse steep spectrum emission ($\alpha < -1.1$) as a candidate for a radio halo which is co-spatial with the centre of the cluster merger. We suggest its origin as a population of old and mildly relativistic electrons left over from radio galaxies throughout the cluster which have been re-accelerated to higher energies by shocks and turbulence induced by the cluster merger event. We also note the discovery of six new giant radio galaxies in the vicinity of Abell 1682., Comment: Accepted in A&A; see journal version for uncompressed figures

Dispersive delays due to the Solar wind introduce excess noise in high-precision pulsar timing experiments, and must be removed in order to achieve the accuracy needed to detect, e.g., low-frequency gravitational waves. In current pulsar timing experiments, this delay is usually removed by approximating the electron density distribution in the Solar wind either as spherically symmetric, or with a two-phase model that describes the contributions from both high- and low-speed phases of the Solar wind. However, no dataset has previously been available to test the performance and limitations of these models over extended timescales and with sufficient sensitivity. Here we present the results of such a test with an optimal dataset of observations of pulsar J0034-0534, taken with the German stations of LOFAR. We conclude that the spherical approximation performs systematically better than the two-phase model at almost all angular distances, with a residual root-mean-square (rms) given by the two-phase model being up to 28% larger than the result obtained with the spherical approximation. Nevertheless, the spherical approximation remains insufficiently accurate in modelling the Solar-wind delay (especially within 20 degrees of angular distance from the Sun), as it leaves timing residuals with rms values that reach the equivalent of 0.3 microseconds at 1400 MHz. This is because a spherical model ignores the large daily variations in electron density observed in the Solar wind. In the short term, broadband observations or simultaneous observations at low frequencies are the most promising way forward to correct for Solar-wind induced delay variations., Comment: Accepted for publication in MNRAS. 17 pages, 13 figures

Context. High-precision pulsar-timing experiments are affected by temporal variations of the Dispersion Measure (DM), which are related to spatial variations in the interstellar electron content. Correcting for DM variations relies on the cold-plasma dispersion law which states that the dispersive delay varies with the squared inverse of the observing frequency. This may however give incorrect measurements if the probed electron content (and therefore the DM) varies with observing frequency, as is predicted theoretically. Aims. We study small-scale density variations in the ionised interstellar medium. These structures may lead to frequency-dependent DMs in pulsar signals and could inhibit the use of lower-frequency pulsar observations to correct time-variable interstellar dispersion in higher-frequency pulsar-timing data. Methods. We used high-cadence, low-frequency observations with three stations from the German LOng-Wavelength (GLOW) consortium, which are part of the LOw Frequency ARray (LOFAR). Specifically, 3.5 years of weekly observations of PSR J2219+4754 are presented. Results. We present the first detection of frequency-dependent DMs towards any interstellar object and a precise multi-year time-series of the time- and frequency-dependence of the measured DMs. The observed DM variability is significant and may be caused by extreme scattering events. Potential causes for frequency-dependent DMs are quantified and evaluated. Conclusions. We conclude that frequency-dependence of DMs has been reliably detected and is caused by small-scale (up to 10s of AUs) but steep density variations in the interstellar electron content. We find that long-term trends in DM variability equally affect DMs measured at both ends of our frequency band and hence the negative impact on long-term high-precision timing projects is expected to be limited., Comment: to be published in A&A (accepted 2019-02-06), 11 pages, 7 figures, update: A&A language editing, typos

Faraday rotation of polarized emission from pulsars measured at radio frequencies provides a powerful tool to investigate the interstellar and interplanetary magnetic fields. However, besides being sensitive to the astrophysical media, pulsar observations in radio are affected by the highly time-variable ionosphere. In this article, the amount of ionospheric Faraday rotation has been computed by assuming a thin layer model. For this aim, ionospheric maps of the free electron density (based on Global Positioning System data) and semi-empirical geomagnetic models are needed. Through the data of five highly polarized pulsars observed with the individual German LOw-Frequency ARray stations, we investigate the performances of the ionospheric modelling. In addition, we estimate the parameters of the systematics and the correlated noise generated by the residual unmodelled ionospheric effects, and show the comparison of the different free-electron density maps. For the best ionospheric maps, we have found that the rotation measure corrections on one-year timescales after subtraction of diurnal periodicity are accurate to $\sim$ 0.06--0.07 rad m$^{-2}$., Comment: 18 pages, 13 figures, 4 tables

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., Comment: 16 figures, 1 table and 22 pages. This paper is part of the LOFAR surveys data release 1 and has been accepted for publication in a special edition of A&A that will appear in Feb 2019, volume 622. The catalogues and images from the data release will be publicly available on lofar-surveys.org upon publication of the journal

The LOFAR Two-metre Sky Survey (LoTSS) is a sensitive, high-resolution 120-168 MHz survey of the Northern sky. The LoTSS First Data Release (DR1) presents 424 square degrees of radio continuum observations over the HETDEX Spring Field (10h45m00s $<$ right ascension $<$ 15h30m00s and 45$^\circ$00$'$00$'$ $<$ declination $<$ 57$^\circ$00$'$00$''$) with a median sensitivity of 71$\mu$Jy/beam and a resolution of 6$''$. In this paper we present photometric redshifts (photo-$z$) for 94.4% of optical sources over this region that are detected in the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) 3$\pi$ steradian survey. Combining the Pan-STARRS optical data with mid-infrared photometry from the Wide-field Infrared Survey Explorer, we estimate photo-$z$s using a novel hybrid photometric redshift methodology optimised to produce the best possible performance for the diverse sample of radio continuum selected sources. For the radio-continuum detected population, we find an overall scatter in the photo-$z$ of 3.9% and an outlier fraction ($\left | z_{\rm{phot}} - z_{\rm{spec}} \right | / (1+z_{\rm{spec}}) > 0.15$) of 7.9%. We also find that, at a given redshift, there is no strong trend in photo-$z$ quality as a function of radio luminosity. However there are strong trends as a function of redshift for a given radio luminosity, a result of selection effects in the spectroscopic sample and/or intrinsic evolution within the radio source population. Additionally, for the sample of sources in the LoTSS First Data Release with optical counterparts, we present rest-frame optical and mid-infrared magnitudes based on template fits to the consensus photometric (or spectroscopic when available) redshift., Comment: 16 figures, 1 table and 18 pages. This paper is part of the LOFAR surveys Data Release 1 and has been accepted for publication in a special edition of A&A that will appear in Feb 2019, Volume 622. The catalogues and images from the data release will be publicly available at http://www.lofar-surveys.org upon publication of the volume

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., Comment: 8 figures, 5 tables and 23 pages. This paper is part of the LOFAR surveys data release 1 and has been accepted for publication in a special edition of A&A that will appear in Feb 2019, volume 622. The catalogues and images from the data release will be publicly available at https://lofar-surveys.org/ upon publication in the journal

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

Despite progress in understanding radio relics, there are still open questions regarding the underlying particle acceleration mechanisms. In this paper we present deep 1--4 GHz VLA observations of CIZA\,J2242.8+5301 ($z=0.1921$), a double radio relic cluster characterized by small projection on the plane of the sky. Our VLA observations reveal, for the first time, the complex morphology of the diffuse sources and the filamentary structure of the northern relic. We discover new faint diffuse radio emission extending north of the main northern relic. Our Mach number estimates for the northern and southern relics, based on the radio spectral index map obtained using the VLA observations and existing LOFAR and GMRT data, are consistent with previous radio and X-ray studies ($\mathcal{M}_{\rm RN}=2.58\pm0.17$ and $\mathcal{M}_{\rm RS}=2.10\pm0.08$). However, color-color diagrams and modelings suggest a flatter injection spectral index than the one obtained from the spectral index map, indicating that projection effects might be not entirely negligible. The southern relic consists of five "arms". Embedded in it, we find a tailed radio galaxy which seems to be connected to the relic. A spectral index flattening, where the radio tail connects to the relic, is also measured. We propose that the southern relic may trace AGN fossil electrons that are re-accelerated at a shock, with an estimated strength of $\mathcal{M}=2.4$. High-resolution mapping of other tailed radio galaxies also supports a scenario where AGN fossil electrons are revived by the merger event and could be related to the formation of some diffuse cluster radio emission., Comment: 29 pages, 19 figures, 5 tables, accepted for publication in the Astrophysical Journal

We present LOFAR $120-168$ MHz images of the merging galaxy cluster Abell 1240 that hosts double radio relics. In combination with the GMRT $595-629$ MHz and VLA $2-4$ GHz data, we characterised the spectral and polarimetric properties of the radio emission. The spectral indices for the relics steepen from their outer edges towards the cluster centre and the electric field vectors are approximately perpendicular to the major axes of the relics. The results are consistent with the picture that these relics trace large-scale shocks propagating outwards during the merger. Assuming diffusive shock acceleration (DSA), we obtain shock Mach numbers of $\mathcal{M}=2.4$ and $2.3$ for the northern and southern shocks, respectively. For $\mathcal{M}\lesssim3$ shocks, a pre-existing population of mildly relativistic electrons is required to explain the brightness of the relics due to the high ($>10$ per cent) particle acceleration efficiency required. However, for $\mathcal{M}\gtrsim4$ shocks the required efficiency is $\gtrsim1\%$ and $\gtrsim0.5\%$, respectively, which is low enough for shock acceleration directly from the thermal pool. We used the fractional polarization to constrain the viewing angle to $\geqslant(53\pm3)^\circ$ and $\geqslant(39\pm5)^\circ$ for the northern and southern shocks, respectively. We found no evidence for diffuse emission in the cluster central region. If the halo spans the entire region between the relics ($\sim1.8\,\text{Mpc}$) our upper limit on the power is $P_\text{1.4 GHz}=(1.4\pm0.6)\times10^{23}\,\text{W}\text{Hz}^{-1}$ which is approximately equal to the anticipated flux from a cluster of this mass. However, if the halo is smaller than this, our constraints on the power imply that the halo is underluminous.