Search

In this paper, we investigate the polarization properties of the double radio relics in PSZ2 G096.88+24.18 using the rotation measure synthesis, and try to constrain the characteristics of the magnetic field that reproduce the observed beam depolarization. Our aim is to understand the nature of the low polarization fraction that characterizes the southern relic with respect to the northern relic. Using new 1-2 GHz VLA observations, we derive the rotation measure and polarization of the two relics by applying the RM synthesis technique, thus solving for bandwidth depolarization in the wide observing bandwidth. To study the effect of beam depolarization, we degraded the image resolution and studied the decreasing trend of polarization fraction with increasing beam size. Finally, we performed 3D magnetic field simulations using multiple models for the magnetic field power spectrum over a wide range of scales, in order to constrain the characteristics of the cluster magnetic field that can reproduce the observed beam depolarization trend. Using RM synthesis, we obtained a polarization fraction of ($18.6 \pm 0.3$)% for the norther relic and ($14.6 \pm 0.1$)% for the southern one. Having corrected for bandwidth depolarization, we infer that the nature of the depolarization for the southern relic is external, and possibly related to the turbulent gas distribution within the cluster, or to the complex spatial structure of the relic. The best-fit magnetic field power spectrum, that reproduces the observed depolarization trend for the southern relic, is obtained for a turbulent magnetic field model, described by a power spectrum derived from cosmological simulations, and defined within the scales of $\Lambda_{\rm{min}}=35~\rm{kpc}$ and $\Lambda_{\rm{max}}=400~\rm{kpc}$. This yields an average magnetic field of the cluster within 1$~\rm{Mpc}^3$ volume of $\sim 2~\rm{\mu G}$., Comment: 17 pages, 10 figures, accepted for publication in Astronomy & Astrophysics on 13th August 2024

A fraction of galaxy clusters harbor diffuse radio sources known as radio halos. The currently adopted scenario for their formation is based on second-order Fermi re-acceleration of seed electrons that is driven by merger-driven turbulence in the intra-cluster medium. This mechanism is expected to be inefficient, which implies that a significant fraction of halos should have very steep ($\alpha < -1.5$) energy spectra. We start investigating the potential and current limitations of the combination of the two surveys conducted by LOFAR, LoTSS (144 MHz) and LoLSS (54 MHz), to probe the origin of radio halos. We follow up the 20 radio halos detected in the DR1 of LoTSS, which covers the HETDEX field, with the LoLSS survey, and we study their spectral properties between 54 and 144 MHz. After the removal of compact sources, 9 halos were excluded due to unreliable halo flux density measurements at 54 MHz. Our main finding is that 7 out of 11 ($\sim$ 64%) exhibit an ultra-steep spectrum ($\alpha < -1.5$), which is a key prediction of turbulent re-acceleration models. We also note a tentative trend for more massive systems to host flatter halos, although the currently poor statistics does not allow for a deeper analysis. Our sample suffers from low angular resolution at 54 MHz, which limits the accuracy of the compact-sources subtraction. Nevertheless, this study is the first step towards providing compelling evidence for the existence of a large fraction of radio halos with very steep spectrum, which is a fundamental prediction of turbulent re-acceleration models. In this regard, the forthcoming second data release of LoLSS, along with the integration of LOFAR international stations and the instrumental upgrade to LOFAR2.0, will improve both the statistics and the low-frequency angular resolution, allowing to conclusively determine the origin of radio halos in galaxy clusters., Comment: A&A Letter; 11 pages, 12 figures

We present the discovery of a 100 kpc low-frequency radio tail behind the nearby group galaxy, NGC 2276. The extent of this tail is a factor of ten larger than previously reported from higher-frequency radio and X-ray imaging. The radio morphology of the galaxy disc and tail suggest that the tail was produced via ram-pressure stripping, cementing NGC 2276 as the clearest known example of ram-pressure stripping in a low-mass group. With multi-frequency imaging, we extract radio continuum spectra between ~50 MHz and 1.2 GHz as a function of projected distance along the tail. All of the spectra are well fit by a simple model of spectral ageing due to synchrotron and inverse-Compton losses. From these fits we estimate a velocity of 870 km/s for the stripped plasma across the plane of the sky, and a three-dimensional orbital velocity of 970 km/s for NGC 2276. The orbital speed that we derive is in excellent agreement with the previous estimate from Rasmussen et al., despite it being derived with a completely independent methodology., Comment: 12 pages, 9 figures, 1 appendix, accepted for publication in A&A

Radio continuum emission from galaxies at gigahertz frequencies can be used as an extinction-free tracer of star formation. However, at frequencies of a few hundred megahertz, there is evidence for low-frequency spectral flattening. We wish to better understand the origin of this low-frequency flattening and, to this end, perform a spatially resolved study of the nearby spiral galaxy M 51. We explore the different effects that can cause flattening of the spectrum towards lower frequencies, such as free-free absorption and cosmic-ray ionisation losses. We used radio continuum intensity maps between 54 and 8350 MHz at eight different frequencies, with observations at 240 MHz from the Giant Metrewave Radio Telescope presented for the first time. We corrected for contribution from thermal free-free emission using an H$\alpha$ map that has been extinction-corrected with 24 $\mu$m data. We fitted free-free absorption models to the radio spectra to determine the emission measure (EM) as well as polynomial functions to measure the non-thermal spectral curvature. The non-thermal low-frequency radio continuum spectrum between 54 and 144 MHz is very flat and even partially inverted, particularly in the spiral arms; contrary, the spectrum at higher frequencies shows the typical non-thermal radio continuum spectrum. However, we do not find any correlation between the EMs calculated from radio and from H$\alpha$ observations; instead, the non-thermal spectral curvature weakly correlates with the HI gas mass surface density. This suggests that cosmic-ray ionisation losses play an important role in the low-frequency spectral flattening. The observed spectral flattening towards low frequencies in M 51 is caused by a combination of ionisation losses and free-free absorption. The reasons for this flattening need to be understood in order to use sub-GHz frequencies as a star-formation tracer.

The evolution of Active Galactic Nuclei (AGN) is closely connected to their host galaxies and surroundings. Via feedback processes, AGN can counteract the cooling of the intracluster medium (ICM) and suppress star formation in their host galaxies. Radio observations at low frequencies provide a glimpse into the history of AGN activity. The Virgo cluster is a substantial reservoir of nearby galaxies and provides an ideal laboratory for the study of AGN as well as their feedback mechanisms. The aim of our work is to characterise the AGN population within the Virgo cluster down to low radio luminosities, constrain the AGN duty cycle and investigate environmental feedback in cluster member galaxies. We analyse 144 MHz and 1.3 GHz radio observations of early-type galaxies from the ACS Virgo Cluster Survey (ACSVCS) taken with LOFAR and MeerKAT. We detect 12 of these galaxies at 144 MHz, 5 of which show clearly extended radio emission. The radio luminosity shows a strong dependence on the stellar mass of the host galaxy, in agreement with previous results. As a notable outlier, the massive elliptical galaxy NGC 4365 ($M_* = 2.2 \times 10^{11} M_\odot$) is not detected as a compact source in the LOFAR observations. Instead, it is surrounded by diffuse, low-surface brightness emission, which hints towards a past phase of stronger nuclear activity. Furthermore, we find a cavity in NGC 4472 (= M 49) inflated by the wide-angle tail only visible in the LOFAR data, which implies that the cavity was created by a past outburst. The corresponding cavity power is of the same order of magnitude as the jet power in the present duty cycle of the AGN., Comment: 17 pages, 10 figures

Merging galaxy clusters often host spectacular diffuse radio synchrotron sources. These sources can be explained by a non-thermal pool of relativistic electrons accelerated by shocks and turbulence in the intracluster medium. The origin of the pool and details of the cosmic ray transport and acceleration mechanisms in clusters are still open questions. Due to the often extremely steep spectral indices of diffuse radio emission, it is best studied at low frequencies. However, the lowest frequency window available to ground-based telescopes (10-30 MHz) has remained largely unexplored, as radio frequency interference and calibration problems related to the ionosphere become severe. Here, we present LOFAR observations from 16 to 168 MHz targeting the famous cluster Abell 2256. In the deepest-ever images at decametre wavelengths, we detect and resolve the radio halo, radio shock and various steep spectrum sources. We measure standard single power-law behaviour for the radio halo and radio shock spectra and find significant spectral index and curvature fluctuations across the radio halo, indicating an inhomogeneous emitting volume. In contrast to the straight power-law spectra of the large-scale diffuse sources, the various AGN-related sources often show extreme steepening towards higher frequencies and flattening towards low frequencies. We also discover a new fossil plasma source with a steep spectrum between 23 and 144 MHz, with $\alpha=-1.9\pm 0.1$. Finally, by comparing radio and gamma-ray observations, we rule out purely hadronic models for the radio halo origin in Abell 2256, unless the magnetic field strength in the cluster is exceptionally high, which is unsupportable by energetic arguments and inconsistent with the knowledge of other cluster magnetic fields., Comment: Accepted for publication in Astronomy & Astrophysics

The largely unexplored decameter radio band (10-30 MHz) provides a unique window for studying a range of astronomical topics, such as auroral emission from exoplanets, inefficient cosmic ray acceleration mechanisms, fossil radio plasma, and free-free absorption. The scarcity of low-frequency studies is mainly due to the severe perturbing effects of the ionosphere. Here we present a calibration strategy that can correct for the ionosphere in the decameter band. We apply this to an observation from the Low Frequency Array (LOFAR) between 16 to 30 MHz . The resulting image covers 330 square degrees of sky at a resolution of 45", reaching a sensitivity of 12 mJy/beam. Residual ionospheric effects cause additional blurring ranging between 60 to 100". This represents an order of magnitude improvement in terms of sensitivity and resolution compared to previous decameter band observations. In the region we surveyed, we have identified four fossil plasma sources. These rare sources are believed to contain old, possibly re-energised, radio plasma originating from previous outbursts of active galactic nuclei. At least three of them are situated near the center of low-mass galaxy clusters. Notably, two of these sources display the steepest radio spectral index among all the sources detected at 23 MHz. This indicates that fossil plasma sources constitute the primary population of steep-spectrum sources at these frequencies, emphasising the large discovery potential of ground-based decameter observations., Comment: 19 pages, 9 figures, 1 table. Submitted version, full version is published by Nature Astronomy

Recent low-frequency radio observations at 140 MHz discovered a 3 Mpc-long bridge of diffuse emission connecting the galaxy clusters Abell 0399 and Abell 0401. We present follow-up observations at 60 MHz to constrain the spectral index of the bridge, which so far has only been detected at 140 and 144 MHz. We analysed deep (~18 hours) LOw Frequency ARray (LOFAR) Low Band Antenna (LBA) data at 60 MHz to detect the bridge at very low frequencies. We then conducted a multi-frequency study with LOFAR HBA data at 144 MHz and uGMRT data at 400 MHz. Assuming second-order Fermi mechanisms for the re-acceleration of relativistic electrons driven by turbulence in the radio bridge regions, we compare the observed radio spectrum with theoretical synchrotron models. The bridge is detected in the 75'' resolution LOFAR image at 60 MHz and its emission fully connects the region between the two galaxy clusters. Between 60 MHz and 144 MHz we found an integrated spectral index value of -1.44 +\- 0.16 for the bridge emission. For the first time, we produced spectral index and related uncertainties maps for a radio bridge. We produce a radio spectrum, which show significant steepening between 144 and 400 MHz. This detection at low frequencies provides important information on the models of particle acceleration and magnetic field structure on very extended scales. The spectral index gives important clues to the origin of inter-cluster diffuse emission. The steepening of the spectrum above 144 MHz can be explained in a turbulent re-acceleration framework, assuming that the acceleration timescales are longer than ~200 Myr., Comment: Letter to the Editor: 7 pages, 4 figures, 2 tables; accepted for publication in Astronomy & Astrophysics (A&A)

Radio observations of a few cool-core galaxy clusters have revealed the presence of diffuse emission on cluster scales, similar to what was found in merging clusters in the form of radio halos. These sources might suggest that a minor merger, while not sufficiently energetic to disrupt the cool core, could still trigger particle acceleration in the intracluster medium on scales of hundreds of kpc. We observed with LOFAR at 144 MHz a sample of twelve cool-core galaxy clusters presenting some level of dynamical disturbances, according to X-ray data. We also performed a systematic search of cold fronts in these clusters, re-analysing archival Chandra data. The clusters PSZ1G139.61+24, A1068 (new detection), MS 1455.0+2232, and RX J1720.1+2638 present diffuse radio emission on a cluster scale. This emission is characterised by a double component: a central mini-halo confined by cold fronts and diffuse emission on larger scales, whose radio power at 144 MHz is comparable to that of radio halos detected in merging systems. The cold fronts in A1068 are a new detection. We also found a candidate plasma depletion layer in this cluster. No sloshing features are found in the other eight clusters. Two of them present a mini-halo, with diffuse radio emission confined to the cluster core. We also found a new candidate mini-halo. Whereas, for the remaining five clusters, we did not detect halo-like emission. For clusters without cluster-scale halos, we derived upper limits to the radio halo power. We found that cluster-scale diffuse radio emission is not present in all cool-core clusters when observed at a low frequency, but it is correlated to the presence of cold fronts. This morphology requires a specific configuration of the merger and so it puts some constraints on the turbulence, which deserves to be investigated in the future with theoretical works., Comment: 34 pages, 30 figures. Accepted for publication in A&A

Radio jets are present in a diverse sample of AGN. However, the mechanisms of jet powering are not fully understood, and it is yet unclear to what extent they obey mass-invariant scaling relations, similar to those found for the triggering and fuelling of X-ray selected AGN. We study the incidence of eROSITA/eFEDS X-ray and LOFAR radio AGN as a function of several stellar mass normalised AGN power indicators. A new sample of radio AGN from the LOFAR-eFEDS survey is defined and we publicly release this catalogue, including host galaxy counterparts from the Legacy Survey DR9, LOFAR radio morphologies and host galaxy properties from the complete, spectroscopic (z<0.4) GAMA09 survey. The fraction of GAMA09 galaxies hosting radio, X-ray and both radio and X-ray AGN are calculated as a function of the specific black hole kinetic ($\lambda_{\rm Jet}$) and radiative ($\lambda_{\rm Edd}$) power. The incidence of eFEDS X-ray AGN as a function of $\lambda_{\rm Edd}$ shows the same mass-invariance as found in past studies. Meanwhile, radio AGN, regardless of their morphology, are more likely to be hosted in more massive galaxies, at all $\lambda_{\rm Jet}$. Across the stellar mass range, the compact radio AGN incidence follows the same power-law distribution, showing that it is not only high mass galaxies that host high power radio AGN and vice versa. On the other hand, the incidence of compact and complex radio AGN is boosted at the highest jet powers, diverging from a simple power-law. Interestingly, this increased incidence cannot be explained by more powerful radio AGN lying in more dense environments which could naturally boost their radio luminosity. Overall, we show that statistical incidence studies are a powerful method to probe disk-jet coupling for different AGN accretion modes, although future work on a more reliable determination of jet power for diverse samples of radio AGN is needed., Comment: 30 pages, 32 figures, accepted for publication in A&A

We present deep total intensity and polarization observations of the Coma cluster at 1.4 and 6.6 GHz performed with the Sardinia Radio Telescope. By combining the single-dish 1.4 GHz data with archival Very Large Array observations we obtain new images of the central radio halo and of the peripheral radio relic where we properly recover the brightness from the large scale structures. At 6.6 GHz we detect both the relic and the central part of the halo in total intensity and polarization. These are the highest frequency images available to date for these radio sources in this galaxy cluster. In the halo, we find a localized spot of polarized signal, with fractional polarization of about 45%. The polarized emission possibly extends along the north-east side of the diffuse emission. The relic is highly polarized, up to 55%, as usually found for these sources. We confirm the halo spectrum is curved, in agreement with previous single-dish results. The spectral index is alpha=1.48 +/- 0.07 at a reference frequency of 1 GHz and varies from alpha ~1.1, at 0.1 GHz, up to alpha ~ 1.8, at 10 GHz. We compare the Coma radio halo surface brightness profile at 1.4 GHz (central brightness and e-folding radius) with the same properties of the other halos, and we find that it has one of the lowest emissivities observed so far. Reanalyzing the relic's spectrum in the light of the new data, we obtain a refined radio Mach number of M=2.9 +/- 0.1., Comment: 27 pages, 18 figures, this article has been accepted for publication in MNRAS Published by Oxford University Press on behalf of the Royal Astronomical Society

Low-frequency radio observations show an increasing number of radio galaxies located in galaxy clusters that display peculiar morphologies and spectral profiles. This is the result of the dynamical interaction of the galaxy with the surrounding medium. Studying this phenomenon is key to understanding the evolution of low-energy relativistic particles in the intracluster medium. We present a multi-frequency study of the three head-tail (HT) radio galaxies and the radio halo in the galaxy cluster ZwCl0634.1+4747. We make use of observations at four frequencies performed with LOFAR LBA (53 MHz), HBA (144 MHz), GMRT (323 MHz) and VLA (1518 MHz) data. The use of extremely low radio frequency observations, such as LOFAR at 53 and 144 MHz, allowed us to detect the extension of the tails up to a distance of ~ 1 Mpc. We extracted spectral profiles along the tails in order to identify possible departures from a pure ageing model, such as the Jaffe-Perola (JP) model, which only involves synchrotron and inverse-Compton losses. We found clear evidence of departures from this simple ageing model, such as surface brightness enhancement and spectral flattening along all of the tails. This can be interpreted as the consequence of particle re-acceleration along the tails. Possible explanations for this behaviour include the interaction between a shock and the radio tails or a turbulence-driven re-acceleration mechanism. We show that the latter scenario is able to reproduce the characteristic features that we observed in our profiles.

We present the discovery of a very extended (550 kpc) and low-surface-brightness ($ 3.3 \mu \mathrm{Jy} \, arcsec^{-2} $ at 144 MHz) radio emission region in Abell 1318. These properties are consistent with its characterisation as an active galactic nucleus (AGN) remnant radio plasma, based on its morphology and radio spectral properties. We performed a broad-band (54 - 1400 MHz) radio spectral index and curvature analysis using LOFAR, uGMRT, and WSRT-APERTIF data. We also derived the radiative age of the detected emission, estimating a maximum age of 250 Myr. The morphology of the source is remarkably intriguing, with two larger, oval-shaped components and a thinner, elongated, and filamentary structure in between, plausibly reminiscent of two aged lobes and a jet. Based on archival {\it Swift} as well as SDSS data we performed an X-ray and optical characterisation of the system, whose virial mass was estimated to be $ \sim 7.4 \times 10^{13} \, \mathrm{M} _{\odot}$. This places A1318 in the galaxy group regime. Interestingly, the radio source does not have a clear optical counterpart embedded in it, thus, we propose that it is most likely an unusual AGN remnant of previous episode(s) of activity of the AGN hosted by the brightest group galaxy ($ \sim 2.6 \times 10^{12} \, \mathrm{M} _{\odot}$), which is located at a projected distance of $\sim$170 kpc in the current epoch. This relatively high offset may be a result of IGrM sloshing sourced by a minor merger. The filamentary morphology of the source may suggest that the remnant plasma has been perturbed by the system dynamics, however, only future deeper X-ray observations will be able to address this question., Comment: 16 pages, 9 figures, accepted by A&A; edited some affiliations

Environmental effects such as ram-pressure stripping (RPS) shape the evolution of galaxies in dense regions. We use the nearby Virgo cluster as a laboratory to study environmental effects on the non-thermal components of star-forming galaxies. We constructed a sample of 17 RPS galaxies in the Virgo cluster and a statistical control sample of 119 nearby galaxies from the Herschel Reference Survey. All objects in these samples are detected in LOFAR 144 MHz observations and come with H$\alpha$ and/or far-UV star formation rate (SFR) estimates. We derived the radio-SFR relations, confirming a clearly super-linear slope of $\approx1.4$. We found that Virgo cluster RPS galaxies have radio luminosities that are a factor of 2-3 larger than galaxies in our control sample. We also investigated the total mass-spectral index relation, where we found a relation for the Virgo cluster RPS galaxies that is shifted to steeper spectral index values by $0.17\pm0.06$. Analyzing the spatially resolved ratio between the observed and the expected radio emission based on the hybrid near-UV + 100$\,\mu$m SFR surface density, we generally observe excess radio emission all across the disk with the exception of a few leading-edge radio-deficient regions. The radio excess and the spectral steepening for the RPS sample could be explained by an increased magnetic field strength if the disk-wide radio enhancement is due to projection effects. For the galaxies that show the strongest radio excesses (NGC 4330, NGC 4396, NGC 4522), a rapid decline of the SFR ($t_\mathrm{quench} \leq 100$ Myr) could be an alternative explanation. We disfavor shock acceleration of electrons as cause for the radio excess since it cannot easily explain the spectral steepening and radio morphology., Comment: 14 pages, 10 figures

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

Turbulence introduced into the intra-cluster medium (ICM) through cluster merger events transfers energy to non-thermal components, and can trigger the formation of diffuse synchrotron radio sources. Typical diffuse sources in the forms of giant radio halos and mini-halos are found in merging and relaxed cool core galaxy clusters, respectively. On the other hand, recent observations have revealed an increasing complexity of the non-thermal phenomenology. Abell 2142 (A2142) is a mildly disturbed cluster that exhibits uncommon thermal and non-thermal properties. It is known to host a hybrid halo consisting of two components (H1 and H2), namely a mini-halo-like and an enigmatic elongated radio halo-like structure. We aim to investigate the properties, origin, and connections of each component. We present deep LOFAR observations of A2142 in the frequency ranges $30-78$ MHz and $120-168$ MHz. With complementary multi-frequency radio and X-ray data, we analyse the radio spectral properties of the halo and assess the connection between the non-thermal and thermal components of the ICM. We detected a third radio component (H3), which extends over the cluster volume on scales $\sim 2$ Mpc, embeds H1 and H2, and has a morphology that roughly follows the thermal ICM distribution. The radio spectral index is moderately steep in H1 ($\alpha=1.09\pm 0.02$) and H2 ($\alpha=1.15\pm 0.02$), but is steeper ($\alpha=1.57\pm 0.20$) in H3. The analysis of the thermal and non-thermal properties allowed us to discuss possible formation scenarios for each radio component. Turbulence from sloshing motions of low-entropy gas on different scales may be responsible for the origin of H1 and H2. We classified H3 as a giant ultra-steep spectrum radio halo, which could trace the residual activity from an old energetic merger and/or inefficient turbulent re-acceleration induced by ongoing minor mergers., Comment: 24 pages (including Appendix), 15 Figures (plus 4 Figures in Appendix), Accepted for publication by A&A

We present the first results of a 21 cm HI line pilot observation carried out with MeerKAT in preparation for the ViCTORIA project, an untargeted survey of the Virgo galaxy cluster. The extraordinary quality of the data in terms of sensitivity and angular resolution (rms~0.65 mJy beam^-1 at ~27"x39" and 11 km/s resolution) allowed us to detect an extended (~10 kpc projected length) low column density (N(HI) < 2.5x10^20 cm^-2) HI gas tail associated with the dwarf irregular galaxy NGC4523 at the northern edge of the cluster. The morphology of the tail and of the stellar disc suggest that the galaxy is suffering a hydrodynamic interaction with the surrounding hot intracluster medium (ICM; ram pressure stripping). The orientation of the trailing tail, the gradient in the HI gas column density at the interface between the cold ISM and the hot ICM, the velocity of the galaxy with respect to that of the cluster, and its position indicate that NGC4523 is infalling for the first time into Virgo from the NNW background of the cluster. Using a grid of hydrodynamic simulations we derive the impact parameters with the surrounding ICM, and estimate that the galaxy will be at pericentre (D~500-600 kpc) in ~1 Gyr, where ram pressure stripping will be able to remove most, if not all, of its gas. The galaxy is located on the star formation main sequence when its star formation rate is derived using Halpha images obtained during the VESTIGE survey, suggesting that NGC4523 is only at the beginning of its interaction with the surrounding environment. A few HII regions are detected in the Halpha images within the HI gas tail outside the stellar disc. Their ages, derived by comparing their Halpha, FUV, NUV, and optical colours with the predictions of SED fitting models, are <30 Myr, and suggest that these HII regions have formed within the stripped gas., Comment: Accepted for publication in A&A

The Virgo cluster is the nearest massive galaxy cluster and thus a prime target to study astrophysical processes in dense large-scale environments. In the radio band, we can probe the non-thermal components of the inter-stellar medium (ISM), intracluster medium (ICM) and of active galactic nuclei (AGN). With the ViCTORIA (Virgo Cluster multi-Telescope Observations in Radio of Interacting galaxies and AGN) project, we are carrying out multiple wide-field surveys of the Virgo cluster at different frequencies. We aim to investigate the impact of the environment on the evolution of galaxies and the contribution of AGN to the ICM-heating, from the inner cluster regions out to beyond the virial radius. We present a survey of the cluster at 120-168 MHz using LOFAR. We image a 132 deg$^2$ region of the cluster, reaching an order of magnitude greater sensitivity than existing wide-field radio surveys of this field at three times higher spatial resolution compared to other low-frequency observations. We developed a tailored data processing strategy to subtract the bright central radio galaxy M87 from the data. This allowed us to correct for the systematic effects due to ionospheric variation as a function of time and direction. In the final mosaic with a resolution of 9"x5", we reach a median noise level of 140 ${\mu}$Jy/beam inside the virial radius and 280 ${\mu}$Jy/beam for the full area. We detect 112 Virgo member galaxies and 114 background galaxies. In at least 18 cases, the radio morphology of the cluster member galaxies shows clear signs of ram-pressure stripping. This includes three previously unreported candidates. In addition, we reveal for the first time 150 kpc long tails from a previous epoch of AGN activity for NGC 4472 (M 49). While no cluster-scale diffuse radio sources are discovered, we find the presence of an extended radio signature of the W$'$-group., Comment: 18 pages, 12 figures, accepted for publication in A&A

Many galaxy clusters show diffuse cluster-scale emission in the form of radio halos, showing that magnetic fields and relativistic electrons are mixed in with the intra-cluster medium (ICM). There is general agreement that the origin of radio halos is connected to turbulence, generated during cluster mergers. Statistical studies of large samples of galaxy clusters in the radio band have the potential to unveil the connection between the properties of radio halos and the mass and dynamics of the host clusters. Previous studies have been limited to massive clusters and based on a small number of radio halos. The aim of this paper is to investigate the scaling relation between the radio power of radio halos and the mass of the host clusters at low frequencies and down to lower cluster masses. We analysed the clusters from the second catalogue of Planck Sunyaev Zel'dovich sources that lie within the 5634 sq deg covered by the second Data Release of the LOFAR Two-meter Sky Survey. We derived the correlation between the radio power and the mass of the host clusters and we investigated the distribution of clusters without radio halos with respect to the correlation. We use X-ray observations to classify the dynamical state of clusters and investigate its role on the power of radio halos. We found a correlation between the power of radio halos at 150 MHz and the mass of the host clusters down to 3e14 Msun. This correlation has a large scatter, part of which can be attributed to the different dynamical states of host clusters. We used two statistical test to show that the distribution of clusters with and without (upper limits) radio halos in the mass-radio power diagram is not compatible with a single correlation and that it is also not compatible with clusters being uniformly distributed below an upper envelope constituted by the correlation., Comment: 12 pages, 7 figures, submitted to A&A

Weakly interacting massive particles (WIMPs) can self-annihilate, thus providing us with a way to indirectly detect dark matter (DM). Dwarf spheroidal (dSph) galaxies are excellent places to search for annihilation signals because they are rich in DM and background emission is low. If magnetic fields in dSph galaxies exist, the particles produced in DM annihilation emit synchrotron radiation. We used the non-detection of 150 MHz radio continuum emission from dSph galaxies with the LOw Frequency ARray (LOFAR) to derive constraints on the annihilation cross section of WIMPs in electron-positron pairs. Our main underlying assumption is that the transport of the cosmic rays can be described by the diffusion approximation, which necessitates the existence of magnetic fields. We used observations of six dSph galaxies in the LOFAR Two-metre Sky Survey (LoTSS). The data were reimaged and radial profiles were generated. We also used stacking to increase the sensitivity. In order to derive upper limits on the WIMP cross section, we injected fake Gaussian sources into the data, which were then detected with 2$\sigma$ significance in the radial profile. These sources represent the lowest emission we would have been able to detect. We present limits from the observations of individual galaxies as well as from stacking. We explored the uncertainty due to the choice of diffusion and magnetic field parameters by constructing three different model scenarios: optimistic (OPT), intermediate (INT), and pessimistic (PES). Assuming monochromatic annihilation into electron-positron pairs, the limits from the INT scenario exclude thermal WIMPs below 20 GeV, and the limits from the OPT scenario even exclude thermal WIMPs below 70 GeV. The INT limits can compete with limits set by Fermi-LAT using $\gamma$-ray observations of multiple dwarf galaxies, and they are especially strong for low WIMP masses., Comment: 15 pages, 10 figures

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

Diffuse cluster-scale synchrotron radio emission is discovered in an increasing number of galaxy clusters in the form of radio halos (RHs), probing the presence of relativistic electrons and magnetic fields in the intra-cluster medium. The favoured scenario to explain their origin is that they trace turbulent regions generated during cluster mergers where particles are re-accelerated. In this framework, RHs are expected to probe cluster dynamics and are predicted to be more frequent in massive systems. Statistical studies are important to study the connection of RHs with cluster dynamics and to constrain theoretical models. Furthermore, low-frequency surveys can shed light on the existence of RHs with very steep radio-spectra, a key prediction of turbulent models. We study the properties of RHs from clusters of the second catalog of Planck Sunyaev Zel'dovich detected sources that lie within the 5634 deg^2 covered by the second Data Release (DR2) of the LOFAR Two-meter Sky Survey. We find that the number of observed RHs, their radio flux density and redshift distributions are in line with what is expected in the framework of the re-acceleration scenario. In addition, the fraction of clusters with RHs increases with the cluster mass, confirming the leading role of the gravitational process of cluster formation in the generation of RHs. These models predict a large fraction of RHs with very steep spectrum in the DR2 Planck sample, this will be tested in future studies, yet a comparison of the occurrence of halos in GMRT and LOFAR samples indeed shows a larger occurrence of RHs at lower frequencies suggesting the presence of a number of very steep spectrum RH that is preferentially detected by LOFAR. Using morphological information we confirm that RHs are preferentially located in merging systems and that the fraction of newly LOFAR discovered RHs is larger in less disturbed systems., Comment: 12 pages, 10 figures, accepted for publication in A&A

Context. It is well-established that shock waves in the intracluster medium launched by galaxy cluster mergers can produce synchrotron emission, which is visible to us at radio frequencies as radio relics. However, the particle acceleration mechanism producing these relics is still not fully understood. It is also unclear how relics relate to radio halos, which trace merger-induced turbulence in the intracluster medium. Aims. We aim to perform the first statistical analysis of radio relics in a mass-selected sample of galaxy clusters, using homogeneous observations. Methods. We analysed all relics observed by the Low Frequency Array Two Metre Sky Survey Data Release 2 (LoTSS DR2) at 144 MHz, hosted by galaxy clusters in the second Planck catalogue of SZ sources (PSZ2). We measured and compared the relic properties in a uniform, unbiased way. In particular, we developed a method to describe the characteristic downstream width in a statistical manner. Additionally, we searched for differences between radio relic-hosting clusters with and without radio halos. Results. We find that, in our sample, $\sim$ 10% of galaxy clusters host at least one radio relic. We confirm previous findings, at higher frequencies, of a correlation between the relic-cluster centre distance and the longest linear size, as well as the radio relic power and cluster mass. However, our findings suggest that we are still missing a population of low-power relics. We also find that relics are wider than theoretically expected, even with optimistic downstream conditions. Finally, we do not find evidence of a single property that separates relic-hosting clusters with and without radio halos., Comment: 28 pages, 17 figures. Accepted for publication in A&A

Context. The details of cosmic-ray transport have a strong impact on galaxy evolution. The peak of the cosmic-ray energy distribution is observable in the radio continuum using the electrons as proxy. Aims. We measure the length that the cosmic-ray electrons (CRE) are transported during their lifetime in the nearby galaxy M 51 across one order of magnitude in cosmic-ray energy (approximately 1-10 GeV). To this end we use new ultra-low frequency observations from the LOw Frequency ARay (LOFAR) at 54 MHz and ancillary data between 144 and 8350 MHz. Methods. As the the CRE originate from supernova remnants, the radio maps are smoothed in comparison to the distribution of the star formation. By convolving the map of the star-formation rate (SFR) surface density with a Gaussian kernel, we can linearise the radio-SFR relation. The best-fitting convolution kernel is then our estimate of the CRE transport length. Results. We find that the CRE transport length increases at low frequencies, as expected since the CRE have longer lifetimes. The CRE transport length is $l_{\rm CRE} = \sqrt{4Dt_{\rm syn}}$, where $D$ is the isotropic diffusion coefficient and $t_{\rm syn}$ is the CRE lifetime as given by synchrotron and inverse Compton losses. We find that the data can be well fitted by diffusion, where $D=(2.14\pm 0.13) \times 10^{28}~\rm cm^2\,s^{-1}$. With $D\propto E^{0.001\pm 0.185}$, the diffusion coefficient is independent of the CRE energy $E$ in the range considered. Conclusions. Our results suggest that the transport of GeV-cosmic ray electrons in the star-forming discs of galaxies is governed by energy-independent diffusion., Comment: Accepted to Astronomy and Astrophysics. 11 pages, 5 figures, 2 tables

The hot plasma within merging galaxy clusters is predicted to be filled with shocks and turbulence that may convert part of their kinetic energy into relativistic electrons and magnetic fields generating synchrotron radiation. Analyzing Low Frequency Array (LOFAR) observations of the galaxy cluster Abell 2255, we show evidence of radio synchrotron emission distributed over very large scales of at least 5 megaparsec. The pervasive radio emission witnesses that shocks and turbulence efficiently transfer kinetic energy into relativistic particles and magnetic fields in a region that extends up to the cluster outskirts. The strength of the emission requires a magnetic field energy density at least 100 times higher than expected from a simple compression of primordial fields, presumably implying that dynamo operates efficiently also in the cluster periphery. It also suggests that nonthermal components may contribute substantially to the pressure of the intracluster medium in the cluster periphery., Comment: Published in Science Advances on Nov 2nd, 2022. 8 pages, 5 images

We present a sample of 373 peaked-spectrum (PS) sources with spectral peaks around 150MHz, selected using a subset of two LOFAR all-sky surveys, the LOFAR Two Meter Sky Survey and the LOFAR LBA Sky Survey. These surveys are the most sensitive low-frequency widefield surveys to date, allowing us to select low-luminosity PS sources. Our sample increases the number of known PS sources in our survey area by a factor 50. The 5GHz luminosity distribution of our PS sample shows we sample the lowest luminosity PS sources to-date by nearly an order of magnitude. Since high-frequency PS sources and compact steep-spectrum sources are hypothesised to be the precursors to large radio galaxies, we investigate whether this is also the case for our sample of low-frequency PS sources. Using optical line emission criteria, we find that our PS sources are predominately high-excitation radio galaxies instead of low-excitation radio galaxies, corresponding to a quickly evolving population. We compute the radio source counts of our PS sample, and find they are scaled down by a factor of $\sim$40 compared to a general sample of radio-loud active galactic nuclei (AGN). This implies that the lifetimes of PS sources are 40 times shorter than large scale radio galaxies, if their luminosity functions are identical. To investigate this, we compute the first radio luminosity function for a homogeneously-selected PS sample. We find that for 144MHz luminosities $\gtrsim 10^{25}$W Hz$^{-1}$, the PS luminosity function has the same shape as an unresolved radio-loud AGN population but shifted down by a factor of $\sim$10. We interpret this as strong evidence that these high-luminosity PS sources evolve into large-scale radio-loud AGN. For local, low-luminosity PS sources, there is a surplus of PS sources, which we hypothesise to be the addition of frustrated PS sources that do not evolve into large-scale AGN., Comment: 21 pages, accepted in Astronomy & Astrophysics

The footprint of LoTSS-DR2 covers 309 PSZ2 galaxy clusters, 83 of which host a radio halo and 26 host a radio relic(s). It provides us an excellent opportunity to statistically study the properties of extended cluster radio sources, especially their connection with merging activities. We aim to quantify cluster dynamic states to investigate their relation with the occurrence of extended radio sources. We also search for connections between intracluster medium (ICM) turbulence and nonthermal characteristics of radio halos in the LoTSS-DR2. We analyzed XMM-Newton and Chandra archival X-ray data and computed concentration parameters and centroid shifts that indicate the dynamic states of the clusters. We also performed a power spectral analysis of the X-ray surface brightness (SB) fluctuations to investigate large-scale density perturbations and estimate the turbulent velocity dispersion. The power spectral analysis results in a large scatter density fluctuation amplitude. We therefore only found a marginal anticorrelation between density fluctuations and cluster relaxation state, and we did not find a correlation between density fluctuations and radio halo power. Nevertheless, the injected power for particle acceleration calculated from turbulent dissipation is correlated with the radio halo power, where the best-fit unity slope supports the turbulent (re)acceleration scenario. Two different acceleration models, transit-time damping and adiabatic stochastic acceleration, cannot be distinguished due to the large scatter of the estimated turbulent Mach number. We introduced a new quantity $[kT\cdot Y_X]_{r_\mathrm{RH}}$, which is proportional to the turbulent acceleration power assuming a constant Mach number. This quantity is strongly correlated with radio halo power, where the slope is also unity., Comment: 20 pages, 13 figures, 5 tables includig appendix. Accepted for publication in A&A. The abstract on arXiv has been shorten to meet the word limit

The central regions of galaxy clusters are permeated by magnetic fields and filled with relativistic electrons. When clusters merge, the magnetic fields are amplified and relativistic electrons are re-accelerated by turbulence in the intra cluster medium. These electrons reach energies of 1 -- 10 GeV and, in the presence of magnetic fields, produce diffuse radio halos that typically cover an area of ~1 square Mpc. Here we report observations of four clusters whose radio halos are embedded in much more extended, diffuse radio emission, filling a volume 30 times larger than that of radio halos. The emissivity in these larger features is about 20 times lower than the emissivity in radio halos. We conclude that relativistic electrons and magnetic fields extend far beyond radio halos, and that the physical conditions in the outer regions of the clusters are quite different from those in the radio halos., Comment: 21 pages, 8 figures, Nature in press

Abell 1033 is a merging galaxy cluster of moderate mass ($M_{500}=3.2\times10^{14} \mathrm{M}_\odot$). It hosts a broad variety of diffuse radio sources linked to different astrophysical phenomena. The most peculiar one is an elongated feature with an ultra-steep spectrum that is the prototype of the category of gently reenergized tails (GReET). Furthermore, the cluster hosts sources that were previously classified as a radio phoenix and a radio halo. We aim to improve the understanding of the cosmic-ray acceleration mechanisms in galaxy clusters in a frequency and mass range that has been poorly explored so far. To investigate the ultra-steep synchrotron emission in the cluster, we performed a full calibration of a LOFAR observation centered at 54 MHz. We analyzed this observation together with re-calibrated data of LoTSS at 144 MHz and an archival GMRT observation at 323 MHz. We performed a spectral study of the radio galaxy tail connected to the GReET to test if the current interpretation of the source is in agreement with observational evidence below 100 MHz. Additionally, we studied the radio halo at different frequencies. We report an extreme spectral curvature for the GReET, the spectral index flattens from $\alpha_{144}^{323} = -4$ to $\alpha_{54}^{144} = -2$. This indicates the presence of a cut-off in the electron energy spectrum. At the cluster center, we detect the radio halo at 54, 144 and at lower significance at 323 MHz. We categorize it as an ultra-steep spectrum radio halo with a spectral index $\alpha = -1.65 \pm 0.17$. Additionally, it is found to be significantly above the radio power-to-cluster mass correlations reported in the literature. Furthermore, the synchrotron spectrum of the halo is found to further steepen between 144 and 323 MHz, in agreement with the presence of a break in the electron spectrum, which is a prediction of homogeneous re-acceleration models., Comment: 15 pages, 13 figures, accepted for publication in A&A

Radio relics are Mpc-size synchrotron sources located in the outskirts of some merging galaxy clusters. Binary-merging systems with favorable orientation may host two almost symmetric relics, named double radio relics. Double radio relics are seen preferentially edge-on and, thus, constitute a privileged sample for statistical studies. Their polarization and Faraday rotation properties give direct access to the relics origin and magnetic fields. In this paper, we present a polarization and Rotation Measure (RM) synthesis study of four clusters hosting double radio relics, namely 8C 0212+703, Abell 3365, PLCK G287.0+32.9, previously missing polarization studies, and ZwCl 2341+0000, for which conflicting results have been reported. We used 1-2 GHz Karl G. Jansky Very Large Array observations. We also provide an updated compilation of known double radio relics with important observed quantities. We studied their polarization and Faraday rotation properties at 1.4 GHz and we searched for correlations between fractional polarization and physical resolution, distance from the cluster center, and shock Mach number. The weak correlations found between these quantities are well reproduced by state-of-the-art magneto-hydrodynamical simulations of radio relics, confirming that merger shock waves propagate in a turbulent medium with tangled magnetic fields. Both external and internal Faraday depolarization should play a fundamental role in determining the polarization properties of radio relics at 1.4 GHz. Although the number of double radio relics with RM information is still low, their Faraday rotation properties (i.e., rest-frame RM and RM dispersion below 40 rad m$^{-2}$ and non-Gaussian RM distribution) can be explained in the scenario in which shock waves with Mach numbers larger than 2.5 propagate along the plane of the sky and compress the turbulent intra-cluster medium., Comment: 23 pages, 5 tables, 14 Figures, 1 appendix, accepted for publication by A&A

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.

Galaxy mergers have been observed to trigger nuclear activity by feeding gas to the central supermassive black hole. One such class of objects are Ultra Luminous InfraRed Galaxies (ULIRGs), which are mostly late stage major mergers of gas-rich galaxies. Recently, large-scale ($\sim$100 kpc) radio continuum emission has been detected in a select number of ULIRGs, all of which also harbour powerful Active Galactic Nuclei (AGN). This hints at the presence of large-scale radio emission being evidence for nuclear activity. Exploring the origin of this radio emission and its link to nuclear activity requires high sensitivity multi-frequency data. We present such an analysis of the ULIRG Mrk 273. Using the International LOFAR telescope (ILT), we detected spectacular large-scale arcs in this system. This detection includes, for the first time, a giant $\sim$190 kpc arc in the north. We propose these arcs are fuelled by a low power radio AGN triggered by the merger. We also identified a bright $\sim$45 kpc radio ridge, which is likely related to the ionised gas nebula in that region. We combined this with high sensitivity data from APERture Tile In Focus (Apertif) and archival data from the Very Large Array (VLA) to explore the spectral properties. The ILT simultaneously allowed us to probe the nucleus at a resolution of $\sim$0.3 arcsec, where we detected three components, and, for the first time, diffuse emission around these components. Combining this with archival high frequency VLA images of the nucleus allowed us to detect absorption in one component, and a steep spectrum radio AGN in another. We then extrapolate from this case study to the importance of investigating the presence of radio emission in more ULIRGs and what it can tell us about the link between mergers and the presence of radio activity., Comment: Accepted for publication in A&A

We study diffuse radio emission in the galaxy cluster A1550, with the aim of constraining particle re-acceleration in the intra-cluster medium. We exploit observations at four different frequencies: 54, 144, 400 and 1400 MHz. To complement our analysis, we make use of archival Chandra X-ray data. At all frequencies we detect an ultra-steep spectrum radio halo ($S_\nu \propto \nu^{-1.6}$) with an extent of 1.2 Mpc at 54 MHz. Its morphology follows the distribution of the thermal intra-cluster medium inferred from the Chandra observation. West of the centrally located head-tail radio galaxy, we detect a radio relic with projected extent of 500 kpc. From the relic, a 600 kpc long bridge departs and connect it to the halo. Between the relic and the radio galaxy, we observe what is most likely a radio phoenix, given its curved spectrum. The phoenix is connected to the tail of the radio galaxy through two arms, which show a nearly constant spectral index for 300 kpc. The halo could be produced by turbulence induced by a major merger, with its axis lying in the NE-SW direction. This is supported by the position of the relic, whose origin could be attributed to a shock propagating along the merger axis. It is possible that the same shock has also produced the phoenix through adiabatic compression, while the bridge could be generated by electrons which were pre-accelerated by the shock, and then re-accelerated by turbulence. Finally, we detect hints of gentle re-energisation in the two arms which depart from the tail of the radio galaxy., Comment: 16 pages, 10 figures

Context. Cosmic rays and magnetic fields are key ingredients in galaxy evolution, regulating both stellar feedback and star formation. Their properties can be studied with low-frequency radio continuum observations, free from thermal contamination. Aims. We define a sample of 76 nearby (< 30 Mpc) galaxies, with rich ancillary data in the radio continuum and infrared from the CHANG-ES and KINGFISH surveys, which will be observed with the LOFAR Two-metre Sky Survey (LoTSS) at 144 MHz. Methods. We present maps for 45 of them as part of the LoTSS data release 2 (LoTSS-DR2), where we measure integrated flux densities and study integrated and spatially resolved radio spectral indices. We investigate the radio-SFR relation, using star-formation rates (SFR) from total infrared and H $\alpha$ + 24-$\mu$m emission. Results. The radio-SFR relation at 144 MHz is clearly super-linear with $L_{144} \propto SFR^{1.4-1.5}$. The mean integrated radio spectral index between 144 and $\approx$1400 MHz is $\langle \alpha\rangle = -0.56 \pm 0.14$, in agreement with the injection spectral index for cosmic ray electrons (CRE). However, the radio spectral index maps show a variation of spectral indices with flatter spectra associated with star-forming regions and steeper spectra in galaxy outskirts and, in particular, in extra-planar regions. We found that galaxies with high star-formation rates (SFR) have steeper radio spectra; we find similar correlations with galaxy size, mass, and rotation speed. Conclusions. Galaxies that are larger and more massive are better electron calorimeters, meaning that the CRE lose a higher fraction of their energy within the galaxies. This explains the super-linear radio-SFR relation, with more massive, star-forming galaxies being radio bright. We propose a semi-calorimetric radio-SFR relation, which employs the galaxy mass as a proxy for the calorimetric efficiency., Comment: Accepted to Astronomy & Astrophysics. 16 pages, 8 figures, 3 tables + 45-page image atlas

We present LOw Frequency ARray observations of the Coma cluster field at 144\,MHz. The cluster hosts one of the most famous radio halos, a relic, and a low surface-brightness bridge. We detect new features that allow us to make a step forward in the understanding of particle acceleration in clusters. The radio halo extends for more than 2 Mpc, which is the largest extent ever reported. To the North-East of the cluster, beyond the Coma virial radius, we discover an arc-like radio source that could trace particles accelerated by an accretion shock. To the West of the halo, coincident with a shock detected in the X-rays, we confirm the presence of a radio front, with different spectral properties with respect to the rest of the halo. We detect a radial steepening of the radio halo spectral index between 144 MHz and 342 MHz, at $\sim 30^{\prime}$ from the cluster centre, that may indicate a non constant re-acceleration time throughout the volume. We also detect a mild steepening of the spectral index towards the cluster centre. For the first time, a radial change in the slope of the radio-X-ray correlation is found, and we show that such a change could indicate an increasing fraction of cosmic ray versus thermal energy density in the cluster outskirts. Finally, we investigate the origin of the emission between the relic and the source NGC 4789, and we argue that NGC4789 could have crossed the shock originating the radio emission visible between its tail and the relic., Comment: submitted to ApJ

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

Radio observations represent a powerful probe of the physics occurring in the intracluster medium (ICM) because they trace the relativistic cosmic rays in the cluster magnetic fields, or within galaxies themselves. By probing the low-energy cosmic rays, low-frequency radio observations are especially interesting because they unveil emission powered by low-efficiency particle acceleration processes, which are believed to play a crucial role in the origin of diffuse radio emission. We investigate the origin of the radio mini-halo at the center of the galaxy cluster 2A0335+096 and its connection to the central galaxy and the sloshing cool core. We also study the properties of the head-tail galaxy GB6 B0335+096 hosted in the cluster to explore the lifecycle of the relativistic electrons in its radio tails. We use new LOw Frequency ARray (LOFAR) observations from the LOFAR Two-meter Sky Survey at 144 MHz to map the low-frequency emission with a high level of detail. The new data were combined with archival Giant Metrewave Radio Telescope and Chandra observations to carry out a multi-wavelength study. We have made the first measurement of the spectral index of the mini-halo ($\alpha=-1.2\pm0.1$ between 144 MHz and 1.4 GHz) and the lobes of the central source ($\alpha\simeq-1.5\pm0.1$ between 144 and 610 MHz). Based on the low-frequency radio emission morphology with respect to the thermal ICM, we propose that the origin of the diffuse radio emission is linked to the sloshing of the cool core. The new data revealed the presence of a Mpc-long radio tail associated with GB6 B0335+096. The observed projected length is a factor 3 longer than the expected cooling length, with evidence of flattening in the spectral index trend along the tail. Therefore, we suggest that the electrons toward the end of the tail are kept alive by the ICM gentle re-acceleration., Comment: 11 pages, 7 figures. Accepted for publication on A&A on 29/11/2021

During their lifetime, galaxy clusters grow through the accretion of matter from the filaments of the large scale structure and from mergers with other clusters. These mergers release a large amount of energy into the intracluster medium (ICM) through merger shocks and turbulence. These phenomena are associated with the formation of radio sources known as radio relics and radio halos, respectively. Radio relics and halos are unique proxies to study the complex properties of these dynamically active regions of clusters and in general the micro physics of the ICM. Abell 3667 is a spectacular examples of a merging system hosting a large pair of radio relics. Due to its proximity (z=0.0553) and large mass, the system enables the study of these sources to a uniquely high level of detail. We observed Abell 3667 with MeerKAT as part of the MeerKAT Galaxy Cluster Legacy Survey. We used these data to study the large scale emission of the cluster, including its polarisation and spectral properties. We present the most detailed view of the radio relic system in Abell 3667 to date, with a resolution reaching 3 kpc. The relics are filled with a network of filaments with different spectral and polarisation properties that are likely associated with multiple regions of particle acceleration and local enhancements of the magnetic field. Conversely, the magnetic field in the space between filaments has strengths close to that expected in unperturbed regions at the same cluster-centric distance. Comparisons with MHD simulations supports the idea of filaments as multiple acceleration sites. Our observations also confirm the presence of an elongated radio halo, developed in the wake of the bullet-like sub-cluster that merged from the South-East. Finally, we associate the process of magnetic draping to a thin polarised radio source surrounding the remnant of the bullet's cool core., Comment: 28 pages, 30 figures, accepted A&A

MeerKAT's large number of antennas, spanning 8 km with a densely packed 1 km core, create a powerful instrument for wide-area surveys, with high sensitivity over a wide range of angular scales. The MeerKAT Galaxy Cluster Legacy Survey (MGCLS) is a programme of long-track MeerKAT L-band (900-1670 MHz) observations of 115 galaxy clusters, observed for $\sim$6-10 hours each in full polarisation. The first legacy product data release (DR1), made available with this paper, includes the MeerKAT visibilities, basic image cubes at $\sim$8" resolution, and enhanced spectral and polarisation image cubes at $\sim$8" and 15" resolutions. Typical sensitivities for the full-resolution MGCLS image products are $\sim$3-5 {\mu}Jy/beam. The basic cubes are full-field and span 4 deg^2. The enhanced products consist of the inner 1.44 deg^2 field of view, corrected for the primary beam. The survey is fully sensitive to structures up to $\sim$10' scales and the wide bandwidth allows spectral and Faraday rotation mapping. HI mapping at 209 kHz resolution can be done at $0

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

Active Galactic Nuclei (AGN) inject large amounts of energy into their host galaxies and surrounding environment, shaping their properties and evolution. In particular, AGN jets inflate cosmic-ray lobes, which can rise buoyantly as light `bubbles' in the surrounding medium, displacing and heating the encountered thermal gas and thus halting its spontaneous cooling. These bubbles have been identified in a wide range of systems. However, due to the short synchrotron lifetime of electrons, the most advanced phases of their evolution have remained observationally unconstrained, preventing us to fully understand their coupling with the external medium, and thus AGN feedback. Simple subsonic hydrodynamic models predict that the pressure gradients, naturally present around the buoyantly rising bubbles, transform them into toroidal structures, resembling mushroom clouds in a stratified atmosphere. The way and timescales on which these tori will eventually disrupt depend on various factors including magnetic fields and plasma viscosity. Here we report LOFAR observations below 200 MHz, sensitive to the oldest radio-emitting particles, showing the late evolution of multiple generations of cosmic-ray AGN bubbles in a galaxy group with unprecedented level of detail. The bubbles' buoyancy power can efficiently offset the radiative cooling of the intragroup medium. However, the bubbles have still not thoroughly mixed with the thermal gas, after hundreds of million years, likely under the action of magnetic fields., Comment: 34 pages, 9 figures, 5 tables, published in Nature Astronomy

Coherent low-frequency ($\lesssim 200$ MHz) radio emission from stars encodes the conditions of the outer corona, mass-ejection events, and space weather. Previous low-frequency searches for radio emitting stellar systems have lacked the sensitivity to detect the general population, instead largely focusing on targeted studies of anomalously active stars. Here we present 19 detections of coherent radio emission associated with known M~dwarfs from a blind flux-limited low-frequency survey. Our detections show that coherent radio emission is ubiquitous across the M~dwarf main sequence, and that the radio luminosity is independent of known coronal and chromospheric activity indicators. While plasma emission can generate the low-frequency emission from the most chromospherically active stars of our sample, the origin of the radio emission from the most quiescent sources is yet to be ascertained. Large-scale analogues of the magnetospheric processes seen in gas-giant planets likely drive the radio emission associated with these quiescent stars. The slowest-rotating stars of this sample are candidate systems to search for star-planet interaction signatures., Comment: Accepted to Nature Astronomy; 31 pages, 9 figures, and 2 tables

Diffuse radio emission at the centre of galaxy clusters has been observed both in merging clusters on scales of Mpc, called giant radio haloes, and in relaxed systems with a cool-core on smaller scales, named mini haloes. Giant radio haloes and mini haloes are thought to be distinct classes of sources. However, recent observations have revealed the presence of diffuse radio emission on Mpc scales in clusters that do not show strong dynamical activity. RX J1720.1+2638 is a cool-core cluster, presenting both a bright central mini halo and a fainter diffuse, steep-spectrum emission extending beyond the cluster core that resembles giant radio halo emission. In this paper, we present new observations performed with the LOFAR Low Band Antennas (LBA) at 54 MHz. These observations, combined with data at higher frequencies, allow us to constrain the spectral properties of the radio emission. The large-scale emission presents an ultra-steep spectrum with $\alpha_{54}^{144}\sim3.2$. The radio emission inside and outside the cluster core have strictly different properties, as there is a net change in spectral index and they follow different radio-X-ray surface brightness correlations. We argue that the large-scale diffuse emission is generated by particles re-acceleration after a minor merger. While for the central mini halo we suggest that it could be generated by secondary electrons and positrons from hadronic interactions of relativistic nuclei with the dense cool-core gas, as an alternative to re-acceleration models., Comment: 13 pages, 10 figures; accepted for publication in MNRAS

We present the first sub-mJy ($\approx0.7$ mJy beam$^{-1}$) survey to be completed below 100 MHz, which is over an order of magnitude deeper than previously achieved for widefield imaging of any field at these low frequencies. The high resolution ($15 \times 15$ arcsec) image of the Bo\"otes field at 34-75 MHz is made from 56 hours of observation with the LOw Frequency ARray (LOFAR) Low Band Antenna (LBA) system. The observations and data reduction, including direction-dependent calibration, are described here. We present a radio source catalogue containing 1,948 sources detected over an area of $23.6$ deg$^2$, with a peak flux density threshold of $5\sigma$. Using existing datasets, we characterise the astrometric and flux density uncertainties, finding a positional uncertainty of $\sim1.2$ arcsec and a flux density scale uncertainty of about 5 per cent. Using the available deep 144-MHz data, we identified 144-MHz counterparts to all the 54-MHz sources, and produced a matched catalogue within the deep optical coverage area containing 829 sources. We calculate the Euclidean-normalised differential source counts and investigate the low-frequency radio source spectral indices between 54 and 144 MHz, both of which show a general flattening in the radio spectral indices for lower flux density sources, from $\sim-0.75$ at 144-MHz flux densities between 100-1000 mJy to $\sim-0.5$ at 144-MHz flux densities between 5-10 mJy, due to a growing population of star forming galaxies and compact core-dominated AGN., Comment: 15 pages, 17 figures. Accepted for publication by A&A. Data available online at lofar-surveys.org

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

Relatively little information is available about the Universe at ultra-low radio frequencies, i.e. below 50 MHz (ULF), although the ULF spectral window contains a wealth of unique diagnostics for studying galactic and extragalactic phenomena. Sub-arcsecond resolution imaging at these frequencies is extremely difficult, due to the long baselines (>1000 km) required and large ionospheric perturbations. We have conducted a pilot project to investigate the ULF performance and potential of the International LOFAR Telescope (ILT), a trans-European interferometric array with baselines up to ~2000 km and observing frequencies down to 10 MHz. We have successfully produced images with sub-arcsecond resolution for 6 radio sources at frequencies down to 30 MHz. This is more than an order of magnitude better resolution than pre-ILT observations at similar frequencies. The six targets that we have imaged (3C 196, 3C 225, 3C 273, 3C 295, 3C 298 and 3C 380) are bright radio sources with compact structures. By comparing our data of 3C 196 and 3C 273 with observations at higher frequencies, we investigate their spatially resolved radio spectral properties. Our success shows that at frequencies down to 30 MHz, sub-arcsecond imaging with the ILT is possible. Further analysis is needed to determine the feasibility of observations of fainter sources or sources with less compact emission., Comment: Accepted to a special issue of A&A on sub-arcsecond imaging with LOFAR. 11 pages, 8 figures, 3 tables

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