Search

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-mode feedback associated with the active galactic nuclei (AGN) at the cores of galaxy clusters injects large amount of energy into the intracluster medium (ICM), offsetting radiative losses through X-ray emission. This mechanism prevents the ICM from rapidly cooling down and fueling extreme starburst activity as it accretes onto the central galaxies, and is therefore a key ingredient in the evolution of galaxy clusters. However, the influence and mode of feedback at high redshifts (z~1) remains largely unknown. Low-frequency sub-arcsecond resolution radio observations taken with the International LOFAR Telescope have demonstrated their ability to assist X-ray observations with constraining the energy output from the AGNs (or "cavity power") in galaxy clusters, thereby enabling research at higher redshifts than before. In this pilot project, we test this hybrid method on a high redshift (0.6

The connection between the thermal and non-thermal properties in galaxy clusters hosting radio halos seems fairly well established. However, a comprehensive analysis of such a connection has been made only for integrated quantities (e.g. $L_X - P_{radio}$ relation). In recent years new-generation radio telescopes have enabled the unprecedented possibility to study the non-thermal properties of galaxy clusters on a spatially resolved basis. Here, we perform a pilot study to investigate the mentioned properties on five targets, by combining X-ray data from the CHEX-MATE project with the second data release from the LOFAR Two meter Sky survey. We find a strong correlation ($r_s \sim 0.7$) with a slope less than unity between the radio and X-ray surface brightness. We also report differences in the spatially resolved properties of the radio emission in clusters which show different levels of dynamical disturbance. In particular, less perturbed clusters (according to X-ray parameters) show peaked radio profiles in the centre, with a flattening in the outer regions, while the three dynamically disturbed clusters have steeper profiles in the outer regions. We fit a model to the radio emission in the context of turbulent re-acceleration with a constant ratio between thermal and non-thermal particles energy density and a magnetic field profile linked to the thermal gas density as $B(r) \propto n_{th}^{0.5}$. We found that this simple model cannot reproduce the behaviour of the observed radio emission., Comment: 23 pages, 10 figures, accepted for publication in A&A

We provide a new $Y_{SZ} - M_{500}^{Y_X}$ scaling relation using a sample of clusters from the Planck Early Sunyaev-Zeldovich (ESZ) catalogue observed in X-rays by Chandra, and compare it to the results of the Planck collaboration from XMM-Newton observations of a subsample of the ESZ. We calibrate a mass bias for a subset of the Planck cosmological cluster sample using published weak-lensing data from the Canadian Cluster Cosmology Project and Multi Epoch Nearby Cluster Survey, for the new scaling relation as well as that from the Planck collaboration. We propose a novel method to account for selection effects and find a mass bias of $(1-b)=0.89\pm0.04$ for the new scaling relation, and $(1-b)=0.76\pm0.04$ for the Planck scaling relation. With these mass biases, we obtain $Y_{SZ} - M_{500}$ scaling relations that we apply to the full Planck cosmological cluster sample, to obtain new constraints on the cosmological parameters. We find fully consistent constraints regardless of the X-ray sample used, with $\sigma_8 = 0.77\pm0.02$, $\Omega_m = 0.31\pm0.02$ and $S_8 \equiv \sigma_8 \sqrt{\Omega_m / 0.3}=0.78\pm0.02$. We also provide constraints with a redshift evolution of the scaling relation fitted from the data instead of fixing it to the self-similar value. We find a redshift evolution significantly deviating from the self-similar value, leading to a higher value of $S_8=0.81\pm0.02$. We compare our results to those from recent analyses based on various cosmological probes, and find that our $S_8$ constraints are competitive with the tightest constraints from the literature. When assuming a self-similar redshift evolution, our constraints are in agreement with most late time probes and in tension with constraints from the CMB primary anisotropies. When fitting the redshift evolution from the data, we find no significant tension with results from either late time probes or the CMB., Comment: 16 pages, 17 figures, submitted to A&A

The JWST North Ecliptic Pole (NEP) Time Domain Field (TDF) is a $>$14 arcmin diameter field optimized for multi-wavelength time-domain science with JWST. It has been observed across the electromagnetic spectrum both from the ground and from space, including with the Hubble Space Telescope (HST). As part of HST observations over 3 cycles (the "TREASUREHUNT" program), deep images were obtained with ACS/WFC in F435W and F606W that cover almost the entire JWST NEP TDF. Many of the individual pointings of these programs partially overlap, allowing an initial assessment of the potential of this field for time-domain science with HST and JWST. The cumulative area of overlapping pointings is ~88 arcmin$^2$, with time intervals between individual epochs that range between 1 day and 4$+$ years. To a depth of $m_{AB}$ $\simeq$ 29.5 mag (F606W), we present the discovery of 12 transients and 190 variable candidates. For the variable candidates, we demonstrate that Gaussian statistics are applicable, and estimate that ~80 are false positives. The majority of the transients will be supernovae, although at least two are likely quasars. Most variable candidates are AGN, where we find 0.42% of the general $z$ $<$ 6 field galaxy population to vary at the $~3\sigma$ level. Based on a 5-year timeframe, this translates into a random supernova areal density of up to ~0.07 transients per arcmin$^2$ (~245 deg$^{-2}$) per epoch, and a variable AGN areal density of ~1.25 variables per arcmin$^2$ (~4500 deg$^{-2}$) to these depths., Comment: 32 pages, 11 figures, 5 tables, 1 Appendix

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.

WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, saw first light in late 2022. WEAVE comprises a new 2-deg field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable ‘mini’ integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrograph covering the wavelength range 366–959 nm at R ∼ 5000, or two shorter ranges at R ∼ 20 000. After summarizing the design and implementation of WEAVE and its data systems, we present the organization, science drivers, and design of a five- to seven-year programme of eight individual surveys to: (i) study our Galaxy’s origins by completing Gaia’s phase-space information, providing metallicities to its limiting magnitude for ∼3 million stars and detailed abundances for ∼1.5 million brighter field and open-cluster stars; (ii) survey ∼0.4 million Galactic-plane OBA stars, young stellar objects, and nearby gas to understand the evolution of young stars and their environments; (iii) perform an extensive spectral survey of white dwarfs; (iv) survey ∼400 neutral-hydrogen-selected galaxies with the IFUs; (v) study properties and kinematics of stellar populations and ionized gas in z < 0.5 cluster galaxies; (vi) survey stellar populations and kinematics in ∼25 000 field galaxies at 0.3 ≲ z ≲ 0.7; (vii) study the cosmic evolution of accretion and star formation using >1 million spectra of LOFAR-selected radio sources; and (viii) trace structures using intergalactic/circumgalactic gas at z > 2. Finally, we describe the WEAVE Operational Rehearsals using the WEAVE Simulator.

In this work, we analysed new LOw Frequency ARray observations of the mini halo in the cluster RBS797, together with archival Very Large Array observations and the recent Chandra results. This cluster is known to host a powerful active galactic nucleus (AGN) at its centre, with two pairs of jets propagating in orthogonal directions. Recent X-ray observations have detected three pairs of shock fronts, connected with the activity of the central AGN. Our aim is to investigate the connection between the mini halo emission and the activity of the central source. We find that the diffuse radio emission is elongated in different directions at 144 MHz (east-west) with respect to 1.4 GHz (north-south), tracing the orientation of the two pairs of jets. The mini halo emission is characterised by an average spectral index $\alpha=-1.02\pm 0.05$. The spectral index profile of the mini halo shows a gradual flattening from the centre to the periphery. Such a trend is unique among the mini halos studied to date, and resembles the spectral index trend typical of particles re-accelerated by shocks. However, the estimated contribution to the radio brightness profile coming from shock re-acceleration is found to be insufficient to account for the radial brightness profile of the mini halo. We propose three scenarios that could explain the observed trend: (i) the AGN-driven shocks are propagating onto an already existing mini halo, re-energising the electrons. We estimate that the polarisation induced by the shocks could be detected at 6 GHz and above; (ii) we could be witnessing turbulent re-acceleration in a high magnetic field cluster; and (iii) the mini halo could have a hadronic origin, in which the particles are injected by Future observations in polarisation would be fundamental to understand the role of shocks and the magnetic field., Comment: A&A, accepted

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

Galaxy cluster mergers are natural consequences of the structure formation in the Universe. Such events involve a large amount of energy ($\sim 10^{63}$ erg) dissipated during the process. Part of this energy can be channelled in particle acceleration and magnetic field amplification, enhancing non-thermal emission of the intra- and inter-cluster environment. Recently, low-frequency observations have detected a bridge of diffuse synchrotron emission connecting two merging galaxy clusters, Abell 399 and Abell 401. Such a result provides clear observational evidence of relativistic particles and magnetic fields in-between clusters. In this work, we have used LOw Frequency ARray (LOFAR) observations at 144 MHz to study for the first time the polarized emission in the A399-A401 bridge region. No polarized emission was detected from the bridge region. Assuming a model where polarization is generated by multiple shocks, depolarization can be due to Faraday dispersion in the foreground medium with respect to the shocks. We constrained its Faraday dispersion to be greater than 0.10 rad m$^{-2}$ at 95% confidence level, which corresponds to an average magnetic field of the bridge region greater than 0.46 nG (or 0.41 nG if we include regions of the Faraday spectrum that are contaminated by Galactic emission). This result is largely consistent with the predictions from numerical simulations for Mpc regions where the gas density is $\sim 300$ times larger than the mean gas density., Comment: 12 pages, 13 figures; accepted in A&A

We present deep XMM-Newton, Karl Jansky Very Large Array, and upgraded Giant Metrewave Radio Telescope observations of Abell 746, a cluster that hosts a plethora of diffuse emission sources that provide evidence for the acceleration of relativistic particles. Our new XMM-Newton images reveal a complex morphology of the thermal gas with several substructures. We observe an asymmetric temperature distribution across the cluster: the southern regions exhibit higher temperatures, reaching ~9 keV, while the northern regions have lower temperatures (below 4 keV), likely due to a complex merger. We find evidence of four surface brightness edges, of which three are merger-driven shock fronts. Combining our new data with the published LOw-Frequency ARray observations has unveiled the nature of diffuse sources in this system. The bright northwest relic shows thin filaments and high degree of polarization with aligned magnetic field vectors. We detect a density jump, aligned with the fainter relic to the north. To the south, we detect high-temperature regions, consistent with shock-heated regions and density jump coincident with the northern tip of the southern radio structure. Its integrated spectrum shows a high-frequency steepening. Lastly, we find that the cluster hosts large-scale radio halo emission. The comparison of the thermal and nonthermal emission reveals an anticorrelation between the bright radio and X-ray features at the center. Our findings suggest that Abell 746 is a complex system that involves multiple mergers., Comment: 22 pages, 13 figures, accepted for publication in ApJ

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

Many relaxed cool-core clusters host diffuse radio emission on scales of hundreds of kiloparsecs: mini-haloes. However, the mechanism responsible for generating them, as well as their connection with central active galactic nuclei, is elusive and many questions related to their physical properties and origins remain unanswered. This paper presents new radio observations of the galaxy cluster Abell 1413 performed with MeerKAT (L-band; 872 to 1712 MHz) and LOFAR HBA (120 to 168 MHz) as part of a statistical and homogeneous census of mini-haloes. Abell 1413 is unique among mini-halo clusters as it is a moderately-disturbed non-cool-core cluster. Our study reveals an asymmetric mini-halo up to 584 kpc in size at 1283 MHz, twice as large as first reported at similar frequencies. The spectral index is flatter than previously reported, with an integrated value of $\alpha = -1.01 \pm 0.06$, shows significant spatial variation, and a tentative radial steepening. We studied the point-to-point X-ray/radio surface brightness correlation to investigate the thermal/non-thermal connection: our results show a strong connection between these components, with a super-linear slope of $b = 1.63 \pm 0.10$ at 1283 MHz and $b = 1.20 \pm 0.12$ at 145 MHz. We also explore the X-ray surface brightness/radio spectral index correlation, finding a slope of $b = 0.59 \pm 0.11$. Both investigations support the evidence of spectral steepening. Finally, in the context of understanding the particle acceleration mechanism, we present a simple theoretical model which demonstrates that hybrid scenarios - secondary electrons (re-)accelerated by turbulence - reproduce a super-linear correlation slope., Comment: 19 pages, 11 figures, 4 tables. Accepted for publication in MNRAS

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

It has long been known that luminous, ultra-steep spectrum radio sources are preferentially associated with massive galaxies at high redshifts. Here we describe a pilot project directed at such objects, to demonstrate the feasibility and importance of using LOFAR to study the most distant forming massive galaxies and protoclusters. We have successfully imaged four high-redshift ($z>2$) high-luminosity radio galaxies with sub-arcsecond resolution, at 144 MHz, using the International LOFAR Telescope (ILT). Our targets were 4C 41.17 ($z=3.8$), the "Anthill", B2 0902+34 ($z=3.4$), 4C 34.34 ($z=2.4$) and 4C 43.15 ($z=2.5$). Their low-frequency morphologies and the spatial distributions of their low-frequency spectral indices have been mapped, and compared with available optical, infrared, and X-ray images. Both for the Anthill at $z = 3.8$ and B2 0902+34 at $z=3.4$, the location of the steepest radio emission coincides with the Ly$\alpha$ emitting ionized gas halo. Our pilot project demonstrates that, because of its outstanding sensitivity and high angular resolution at low frequencies, the ILT is a unique facility for studying the co-evolution and interaction of massive galaxies, galaxy clusters, and supermassive black holes in the early Universe., Comment: 10 pages, 10 figures, 2 tables

Cosmological simulations predict the presence of warm hot thermal gas in the cosmic filaments that connect galaxy clusters. This gas is thought to constitute an important part of the missing baryons in the Universe. In addition to the thermal gas, cosmic filaments could contain a population of relativistic particles and magnetic fields. A detection of magnetic fields in filaments can constrain early magnetogenesis in the cosmos. So far, the resulting diffuse synchrotron emission has only been indirectly detected. We present our search for thermal and non-thermal diffuse emission from inter-cluster regions of 106 paired galaxy clusters by stacking the $0.6-2.3$~keV X-ray and 144~MHz radio data obtained with the eROSITA telescope on board the Spectrum-Roentgen-Gamma (SRG) observatory and LOw Frequency ARray (LOFAR), respectively. The stacked data do not show the presence of X-ray and radio diffuse emission in the inter-cluster regions. This could be due to the sensitivity of the data sets and/or the limited number of cluster pairs used in this study. Assuming a constant radio emissivity in the filaments, we find that the mean radio emissivity is not higher than $1.2\times10^{-44}\,{\rm erg \, s^{-1} \, cm^{-3} \, Hz^{-1}}$. Under equipartition conditions, our upper limit on the mean emissivity translates to an upper limit of $\sim75\,{\rm nG}$ for the mean magnetic field strength in the filaments, depending on the spectral index and the minimum energy cutoff. We discuss the constraint for the magnetic field strength in the context of the models for the formation of magnetic fields in cosmic filaments., Comment: Accepted for publication on MNRAS on June 5, 2023

Radio galaxies are a subclass of active galactic nuclei in which accretion onto the supermassive black hole releases energy via relativistic jets. The jets are not constantly active throughout the life of the host galaxy and alternate between active and quiescent phases. Remnant radio galaxies are detected during a quiescent phase and define a class of unique sources to constrain the AGN duty cycle. We present, a spatially resolved radio analysis of the radio galaxy associated with NGC 6086 and constraints on the spectral age of the diffuse emission to investigate the duty cycle and evolution of the source. We use three new low-frequency, high-sensitivity observations, performed with the Low Frequency Array at 144 MHz and with the upgraded Giant Metrewave Radio Telescope at 400 MHz and 675 MHz. To these, we add two Very Large Array archival observations at 1400 and 4700 MHz. In the new observations, we detect a second pair of larger lobes and three regions with a filamentary morphology. We analyse the spectral index trend in the inner remnant lobes and see systematic steeper values at the lower frequencies compared to the GHz ones. Steeper spectral indices are found in the newly detected outer lobes (up to 2.1), as expected if they trace a previous phase of activity of the AGN. However, the differences between the spectra suggest different dynamical evolution within the intragroup medium during their expansion and/or different magnetic field values. We place constraints on the age of the inner and outer lobes and derive the duty cycle of the source. This results in a total active time of $\sim$39%. The filamentary structures have a steep spectral index ($\sim$1) without any spectral index trend and only one of them shows a steepening in the spectrum. Their origin is not yet clear, but they may have formed due to the compression of the plasma or due to magnetic field substructures., Comment: 19 pages, 11(+5) figures, 6 tables. Astronomy & Astrophysics: Received: 30 May 2022 - Accepted: 21 June 2022

Wide-field radio continuum observations of galaxy clusters are revealing an increasing number of spiral galaxies hosting tens of kpc-long radio tails produced by the nonthermal interstellar medium being displaced by the ram pressure. We present a semi-empirical model for the multi-frequency radio continuum emission from ram pressure stripped tails based on the pure synchrotron cooling of a radio plasma moving along the stripping direction with a uniform velocity. We combine LOFAR and uGMRT observations at 144 and 400 MHz to study the flux density and spectral index profiles of the radio tails of 7 galaxies in Abell 2255, and use the model to reproduce the flux density and spectral index profiles, and infer the stripped radio plasma velocity. For 5 out of 7 galaxies we observe monotonic decrease in both flux density and spectral index up to $~30$ kpc from their stellar disk. Our model reproduces the observed trends with a radio plasma bulk projected velocity between 160 and 430 km s$^{-1}$. This result represents the first indirect measure of the stripped, nonthermal interstellar medium velocity. The observed spectral index trends indicate that the synchrotron cooling is faster than the adiabatic expansion losses, thus suggesting that the stripped radio plasma can survive for a few tens of Myr outside of the stellar disk. This provides a lower limit for the lifetime of the stripped ISM outside of the disk. As a proof of concept, we use the best-fit velocities to constrain the galaxies' 3D velocity in the cluster to be in the 300-1300 km s$^{-1}$. We estimate the ram pressure affecting these galaxies to be between 0.1 and 2.9 $\times10^{-11}$ erg cm$^{-3}$, and measure the inclination between their stellar disk and the ram pressure wind., Comment: 15 pages, 9 figures, 5 tables. Accepted for publication on A&A on May 31st 2023

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

Source classifications, stellar masses and star formation rates are presented for 80,000 radio sources from the first data release of the Low Frequency Array Two-metre Sky Survey (LoTSS) Deep Fields, which represents the widest deep radio survey ever undertaken. Using deep multi-wavelength data spanning from the ultraviolet to the far-infrared, spectral energy distribution (SED) fitting is carried out for all of the LoTSS-Deep host galaxies using four different SED codes, two of which include modelling of the contributions from an active galactic nucleus (AGN). Comparing the results of the four codes, galaxies that host a radiative AGN are identified, and an optimised consensus estimate of the stellar mass and star-formation rate for each galaxy is derived. Those galaxies with an excess of radio emission over that expected from star formation are then identified, and the LoTSS-Deep sources are divided into four classes: star-forming galaxies, radio-quiet AGN, and radio-loud high-excitation and low-excitation AGN. Ninety-five per cent of the sources can be reliably classified, of which more than two-thirds are star-forming galaxies, ranging from normal galaxies in the nearby Universe to highly-starbursting systems at z>4. Star-forming galaxies become the dominant population below 150-MHz flux densities of about 1 mJy, accounting for 90 per cent of sources at a 150-MHz flux density of 100 microJy. Radio-quiet AGN comprise around 10 per cent of the overall population. Results are compared against the predictions of the SKADS and T-RECS radio sky simulations, and improvements to the simulations are suggested., Comment: Accepted for publication in MNRAS. Catalogues available at www.lofar-surveys.org/deepfields.html. 27 pages

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

Models of radio galaxy physics have been primarily based on high frequency ($\geqslant$1 GHz) observations of their jets, hotspots, and lobes. Without highly resolved low frequency observations, which provide information on older plasma, our understanding of the dynamics of radio galaxies and their interaction with their environment is limited. Here, we present the first sub-arcsecond (0.3") resolution images at 144 MHz of two powerful radio galaxies situated in rich cluster environments, namely 3C 34 and 3C 320, using the International Low Frequency Array Telescope. We detect for the first time at low frequencies a plethora of structures in these objects, including strikingly large filaments across the base of the lobes in both sources, which are spatially associated with dense regions in the ambient medium. For 3C 34 we report a spectral flattening in the region of the central filament, suggesting that the origin of the filaments is related to the presence of large-scale ordered magnetic fields. We also report periodic total intensity and spectral index banding of diffuse emission in the eastern lobe, seen for the first time in radio galaxy lobes. The hotspot complexes are resolved into multiple fragments of varying structure and spectral index; we discuss the implications for particle acceleration and jet termination models. We find at most smooth gradients in the spectral behaviour of the hotspot structure suggesting that particle acceleration, if present, may be occurring throughout the complex, in contrast to simple models, but different jet termination models may apply to both sources., Comment: 16 pages, 7 figures, 2 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society, 2023 February 01

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

Extended radio sources in the sky require a dense sampling of short baselines to be properly imaged by interferometers. This problem arises in many areas of radio astronomy, such as in the study of galaxy clusters, which may host Mpc-scale diffuse synchrotron sources in the form of radio halos. In clusters where no radio halos are detected, owing to intrinsic absence of emission or extrinsic (instrumental and/or observational) effects, it is possible to determine upper limits. We consider a sample of Planck galaxy clusters from the Second Data Release of the LOFAR Two Meter Sky Survey (LoTSS-DR2) where no radio halos are detected. We use this sample to test the capabilities of LOFAR to recover diffuse extended emission and derive upper limits. Through the injection technique, we simulate radio halos with various surface brightness profiles. We then predict the corresponding visibilities and image them along with the real visibilities. This method allows us to test the fraction of flux density losses owing to inadequate uv-coverage and obtain thresholds at which the mock emission becomes undetectable by visual inspection. The dense uv-coverage of LOFAR at short spacings allows to recover $\gtrsim90\%$ of the flux density of targets with sizes up to $\sim 15'$. We find a relation that provides upper limits based on the image noise and extent (in terms of number of beams) of the mock halo. This relation can be safely adopted to obtain upper limits without injecting when artifacts introduced by the subtraction of the discrete sources are negligible in the central region of the cluster. Otherwise, the injection process and visual inspection of the images are necessary to determine more reliable limits. Through these methods, we obtain upper limits for 75 clusters to be exploited in ongoing statistical studies., Comment: 20 pages (including Appendices). 21 Figures. 4 Tables. Accepted for publication by A&A. For the project Webpage see https://lofar-surveys.org/planck_dr2.html

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

A Faraday rotation measure (RM) catalogue, or RM Grid, is a valuable resource for the study of cosmic magnetism. Using the second data release (DR2) from the LOFAR Two-metre Sky Survey (LoTSS), we have produced a catalogue of 2461 extragalactic high-precision RM values across 5720 deg$^{2}$ of sky (corresponding to a polarized source areal number density of $\sim$0.43 deg$^{-2}$). The linear polarization and RM properties were derived using RM synthesis from the Stokes $Q$ and $U$ channel images at an angular resolution of 20'' across a frequency range of 120 to 168 MHz with a channel bandwidth of 97.6 kHz. The fraction of total intensity sources ($>1$ mJy beam$^{-1}$) found to be polarized was $\sim$0.2%. The median detection threshold was 0.6 mJy beam$^{-1}$ ($8\sigma_{QU}$), with a median RM uncertainty of 0.06 rad m$^{-2}$ (although a systematic uncertainty of up to 0.3 rad m$^{-2}$ is possible, after the ionosphere RM correction). The median degree of polarization of the detected sources is 1.8%, with a range of 0.05% to 31%. Comparisons with cm-wavelength RMs indicate minimal amounts of Faraday complexity in the LoTSS detections, making them ideal sources for RM Grid studies. Host galaxy identifications were obtained for 88% of the sources, along with redshifts for 79% (both photometric and spectroscopic), with the median redshift being 0.6. The focus of the current catalogue was on reliability rather than completeness, and we expect future versions of the LoTSS RM Grid to have a higher areal number density. In addition, 25 pulsars were identified, mainly through their high degrees of linear polarization., Comment: 20 pages, 23 figures, 3 tables. Accepted for publication in MNRAS. Go to https://lofar-mksp.org/data/ to download the catalogue and associated data products

We present the first deep (72 hours of observations) radio image of the Euclid Deep Field North (EDFN) obtained with the LOw-Frequency ARray (LOFAR) High Band Antenna (HBA) at 144 MHz. The EDFN is the latest addition to the LOFAR Two-Metre Sky Survey (LoTSS) Deep Fields and these observations represent the first data release for this field. The observations produced a 6" resolution image with a central r.m.s. noise of $32\,\mu$Jy\,beam$^{-1}$. A catalogue of $\sim 23,000$ radio sources above a signal-to-noise ratio (SNR) threshold of 5 is extracted from the inner circular 10 deg$^2$ region. We discuss the data analysis and we provide a detailed description of how we derived the catalogue of radio sources and on the issues related to direction-dependent calibration and their effects on the final products. Finally, we derive the radio source counts at 144 MHz in the EDFN using catalogues of mock radio sources to derive the completeness correction factors. The source counts in the EDFN are consistent with those obtained from the first data release of the other LoTSS Deep Fields (ELAIS-N1, Lockman Hole and Bootes), despite the different method adopted to construct the final catalogue and to assess its completeness., Comment: accepted for publication in A&A

Radio-loud active galactic nuclei (RLAGN) play an important role in the evolution of galaxies through the effects on their environment. The two major morphological classes are core-bright (FRI) and edge-bright (FRII) sources. With the LOw-Frequency ARray (LOFAR) we compare the FRI and FRII evolution down to lower flux densities and with larger samples than before with the aim to examine the cosmic space density evolution for FRIs and FRIIs by analyzing their space density evolution between L_150~10^24.5 W/Hz and L_150~10^28.5 W/Hz and up to z=2.5. We construct radio luminosity functions (RLFs) from FRI and FRII catalogues based on recent data from LOFAR at 150MHz to study the space densities as a function of radio luminosity and redshift. To partly correct for selection biases and completeness, we simulate how sources appear at a range of redshifts. We report a space density enhancement from low to high redshift for FRI and FRII sources brighter than L_150~10^27 W/Hz. This is possibly related to the higher gas availability in the earlier denser universe. The constant FRI/FRII space density ratio evolution as a function of radio luminosity and redshift in our results suggests that the jet-disruption of FRIs might be primarily caused by events occurring on scales within the host galaxy, rather than being driven by changes in the overall large-scale environment. Remaining selection biases in our results also highlight the need to resolve more sources at angular scales below 40 arcsec and therefore strengthens the motivation for the further development and automation of the calibration and imaging pipeline of LOFAR data to produce images at sub-arcsecond resolution.

$Aims:$ We used the nearby Coma Cluster as a laboratory in order to probe the impact of ram pressure on star formation as well as to constrain the characteristic timescales and velocities for the stripping of the non-thermal ISM. $Methods:$ We used high-resolution ($6.5'' \approx 3\,\mathrm{kpc}$), multi-frequency ($144\,\mathrm{MHz} - 1.5\,\mathrm{GHz}$) radio continuum imaging of the Coma Cluster to resolve the low-frequency radio spectrum across the discs and tails of 25 ram pressure stripped galaxies. With resolved spectral index maps across these galaxy discs, we constrained the impact of ram pressure perturbations on galaxy star formation. We measured multi-frequency flux-density profiles along each of the ram pressure stripped tails in our sample. We then fit the resulting radio continuum spectra with a simple synchrotron aging model. $Results:$ We showed that ram pressure stripped tails in Coma have steep ($-2 \lesssim \alpha \lesssim -1$) spectral indices. The discs of galaxies undergoing ram pressure stripping have integrated spectral indices within the expected range for shock acceleration from supernovae ($-0.8 \lesssim \alpha \lesssim -0.5$), though there is a tail towards flatter values. In a resolved sense, there are gradients in spectral index across the discs of ram pressure stripped galaxies in Coma. These gradients are aligned with the direction of the observed radio tails, with the flattest spectral indices being found on the `leading half'. From best-fit break frequencies we estimated the projected plasma velocities along the tail to be on the order of hundreds of kilometers per second, with the precise magnitude depending on the assumed magnetic field strength., Comment: 18 pages, 10 figures, 2 appendices, accepted for publication in A&A

We report on the detection of diffuse radio emission with peculiar morphology in the central region of the galaxy cluster Abell 2657. The most striking feature identified in our 144 MHz LOFAR image is a bifurcated radio arc that extends for a projected size of 150-200 kpc. From the analysis of XMM-Newton data, we find clear evidence of gas sloshing in the cluster and a possible dip in X-ray surface brightness between the two radio arcs which deserves confirmation. Interestingly, the synchrotron emission of the bifurcated radio arc is stretched along the sloshing spiral. We compare our observational results with numerical simulations of non-thermal components interacting with gas motions. We suggest that the detected emission may trace a radio bubble shredded by gas sloshing, where relativistic electrons and magnetic fields are expected to be stretched and stirred as a consequence of tangential flows induced by the spiralling gas motion. Lastly, we report on the presence of two thin (6-7 kpc in width) and parallel strands of radio emission embedded in the outer arc that are morphologically similar to the emerging population of non-thermal filaments observed in galaxy clusters, radio galaxies, and the Galactic centre. While this work further demonstrates the complex interplay between thermal and non-thermal components in the intracluster medium, follow-up observations in radio and X-rays are required to firmly determine the origin of the features observed in Abell 2657., Comment: 12 pages, 9 figures, 1 table. Accepted for publication in MNRAS

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

MS 0735.6+7421 ($z = 0.216$) is a massive cool core galaxy cluster hosting one of the most powerful active galactic nuclei (AGN) outbursts known. The radio jets of the AGN have carved out an unusually large pair of X-ray cavities, each reaching a diameter of $200$ kpc. This makes MS 0735.6+7421 a unique case to investigate active galactic nuclei feedback processes, as well as other cluster astrophysics at radio wavelengths. We present new low-radio-frequency observations of MS 0735.6+7421 taken with the Karl G. Jansky Very Large Array (VLA): 5 hours of P-band ($224-480$ MHz) and 5 hours of L-band ($1-2$ GHz) observations, both in C configuration. Our VLA P-band ($224-480$ MHz) observations reveal the presence of a new diffuse radio component reaching a scale of $\sim$ $900$ kpc in the direction of the jets and of $\sim$ $500$ kpc in the direction perpendicular to the jets. This component is centered on the cluster core and has a radio power scaled at $1.4$ GHz of $P_{1.4\text{ GHz}} = (4\pm2)\times 10^{24}$ WHz$^{-1}$. Its properties are consistent with those expected from a radio mini-halo as seen in other massive cool core clusters, although it may also be associated with radio plasma that has diffused out of the X-ray cavities. Observations at higher spatial resolution are needed to fully characterize the properties and nature of this component. We also suggest that if radio mini-halos originate from jetted activity, we may be witnessing the early stages of this process., Comment: 11 pages, 7 figures, submitted to MNRAS

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

Jets from active galactic nuclei (AGN) are known to recurrently enrich their surrounding medium with mildly-relativistic particles and magnetic fields. Here, we present a detailed multi-frequency analysis of the nearby (z=0.01646) galaxy group NGC 507. In particular, we present new high-sensitivity and high spatial resolution radio images in the frequency range 144-675 MHz obtained using LOFAR and uGMRT observations. These reveal the presence of previously undetected diffuse radio emission with complex, filamentary morphology, likely related to a previous outburst of the central galaxy. Based on spectral ageing considerations, we derived that the plasma was first injected by the AGN 240-380 Myr ago and is now cooling. Our analysis of deep archival XMM-Newton data confirms that the system is dynamically disturbed, as previously suggested. We detect two discontinuities in the X-ray surface brightness distribution (in East and South direction) tracing a spiral pattern, which we interpret as cold fronts produced by sloshing motions. The remarkable spatial coincidence observed between the newly-detected arc-like radio filament and the southern concave X-ray discontinuity strongly suggests that the remnant plasma has been displaced by the sloshing motions on large scales. Overall, NGC 507 represents one of the clearest examples known to date in which a direct interaction between old AGN remnant plasma and the external medium is observed in a galaxy group. Our results are consistent with simulations, which suggest that filamentary emission can be created by the cluster/group weather disrupting AGN lobes and spreading their relativistic content into the surrounding medium., Comment: 21 pages, 15 figures, 3 tables, A&A accepted

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 (∼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 ∼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 ∼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 ∼0.3″, 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. © P. Kukreti et al. 2022.

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

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

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

We examine the particle acceleration mechanism in the Mpc-scale bridge between Abell 399 and Abell 401 and assess in particular if the synchrotron emission originates from first-order or second-order Fermi re-acceleration. We use deep (~40 hours) LOw-Frequency ARray (LOFAR) observations from Abell 399 and Abell 401 and apply improved direction-dependent calibration to produce deep radio images at three different resolutions at 144 MHz. With a point-to-point analysis we find in the bridge trends between the radio emission from our new maps and X-ray emission from an XMM Newton observation. By analyzing our observations and results, we argue that second-order Fermi re-acceleration is currently the most favoured process to explain the emission from the radio bridge, where past AGN activity may be responsible for the supply of fossil plasma needed for in-situ re-acceleration. The radio halos from Abell 401 and Abell 399 are also consistent with a second-order Fermi re-acceleration model.

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

Identifying active galactic nuclei (AGN) and isolating their contribution to a galaxy's energy budget is crucial for studying the co-evolution of AGN and their host galaxies. Brightness temperature ($T_b$) measurements from high-resolution radio observations at GHz frequencies are widely used to identify AGN. Here we investigate using new sub-arcsecond imaging at 144 MHz with the International LOFAR Telescope to identify AGN using $T_b$ in the Lockman Hole field. We use ancillary data to validate the 940 AGN identifications, finding 83 percent of sources have AGN classifications from SED fitting and/or photometric identifications, yielding 160 new AGN identifications. Considering the multi-wavelength classifications, brightness temperature criteria select over half of radio-excess sources, 32 percent of sources classified as radio-quiet AGN, and 20 percent of sources classified as star-forming galaxies. Infrared colour-colour plots and comparison with what we would expect to detect based on peak brightness in 6 arcsec LOFAR maps, imply that the star-forming galaxies and sources at low flux densities have a mixture of star-formation and AGN activity. We separate the radio emission from star-formation and AGN in unresolved, $T_b$-identified AGN with no significant radio excess and find the AGN comprises $0.49\pm 0.16$ of the radio luminosity. Overall the non-radio excess AGN show evidence for having a variety of different radio emission mechanisms, which can provide different pathways for AGN and galaxy co-evolution. This validation of AGN identification using brightness temperature at low frequencies opens the possibility for securely selecting AGN samples where ancillary data is inadequate., Comment: Accepted for publication in MNRAS. 19 pages, 15 figures, 2 tables

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

It has been well established that galaxy clusters have magnetic fields. The exact properties and origin of these magnetic fields are still uncertain even though these fields play a key role in many astrophysical processes. Various attempts have been made to derive the magnetic field strength and structure of nearby galaxy clusters using Faraday rotation of extended cluster radio sources. This approach needs to make various assumptions that could be circumvented when using background radio sources. However, because the number of polarised radio sources behind clusters is low, at the moment such a study can only be done statistically. In this paper, we investigate the depolarisation of radio sources inside and behind clusters in a sample of 124 massive clusters at $z<0.35$ observed with the Karl G. Jansky Very Large Array. We detect a clear depolarisation trend with the cluster impact parameter, with sources at smaller projected distances to the cluster centre showing more depolarisation. By combining the radio observations with ancillary X-ray data from Chandra, we compare the observed depolarisation with expectations from cluster magnetic field models using individual cluster density profiles. The best-fitting models have a central magnetic field strength of $5-10\,\mu$G with power-law indices between $n=1$ and $n=4$. We find no strong difference in the depolarisation trend between sources embedded in clusters and background sources located at similar projected radii, although the central region of clusters is still poorly probed by background sources. We also examine the depolarisation trend as a function of cluster properties such as the dynamical state, mass, and redshift. Our findings show that the statistical depolarisation of radio sources is a good probe of cluster magnetic field parameters. [abridged], Comment: Replaced with Published version (A&A). 32 pages, 34 figures

We report the discovery of a ring of low surface brightness radio emission around the Calvera pulsar, a high Galactic latitude, isolated neutron star, in the LOFAR Two-metre Sky Survey (LoTSS). It is centered at $\alpha=14\mathrm{h}11\mathrm{m}12.6\mathrm{s}$, $\delta=+79^\mathrm{o}23'15"$, has inner and outer radii of $14.2'$ and $28.4'$, and an integrated flux density at 144 MHz of $1.08\pm0.15$ Jy. The ring center is offset by $4.9'$ from the location of the Calvera pulsar. H$\alpha$ observations with the Isaac Newton Telescope show no coincident optical emission, but do show a small ($\sim20"$) optical structure internal to the ring. We consider three possible interpretations for the ring: that it is an H~II region, a supernova remnant (SNR), or an Odd Radio Circle (ORC). The positional coincidence of the ring, the pulsar, and an X-ray-emitting non-equilibrium ionisation plasma previously detected, lead us to prefer the SNR interpretation. If the source is indeed a SNR and its association with the Calvera pulsar is confirmed, then Calvera's SNR, or G118.4+37.0, will be one of few SNRs in the Galactic halo., Comment: Accepted in A&A

We present high-resolution, high-sensitivity upgraded Giant Metrewave Radio Telescope observations of the Coma cluster (A1656) at 250-500 MHz and 550-850 MHz. At 250-500 MHz, 135 sources have extensions $>$ 0.45 arcmin (with peak-to-local-noise ratio $> 4$). Of these, 24 sources are associated with Coma-member galaxies. In addition, we supplement this sample of 24 galaxies with 20 ram pressure stripped galaxies from Chen et al. (2020, eight are included in the original extended radio source sample) and an additional five are detected and extended. We present radio morphologies, radio spectra, spectral index maps, and equipartition properties for these two samples. In general, we find the equipartition properties lie within a narrow range (e.g., $P_{\rm min}$ = 1-3 dynes cm$^{-2}$). Only NGC 4874, one of the two brightest central Coma cluster galaxies, has a central energy density and pressure about five times higher and a radio source age about 50 % lower than that of the other Coma galaxies. We find a diffuse tail of radio emission trailing the dominant galaxy of the merging NGC 4839 group that coincides with the "slingshot" tail, seen in X-rays. The southwestern radio relic, B1253$+$275, has a large extent $\approx$ 32$^\prime$ $\times$ 10$^\prime$ ($\simeq$ 1.08 $\times$ 0.34 Mpc$^2$). For NGC 4789, whose long radio tails merge into the relic and may be a source of its relativistic seed electrons, and we find a transverse radio spectral gradient, a steepening from southwest to northeast across the width of the radio source. Finally, radio morphologies of the extended and RPS samples suggest that these galaxies are on their first infall into Coma on (predominantly) radial orbits., Comment: 35 pages (includes Appendices), 15 figures (includes 2 figures in Appendix, some Figures are of lower quality), 6 tables (includes 1 table in Appendix), accepted for publication in the Astrophysical Journal