Search

We present a detailed analysis of the thermal X-ray emission from the intracluster medium (ICM) in the cool-core galaxy cluster Abell 2667 ($z=0.23$). Our goal is to detect low-temperature ($<2$ keV) X-ray emitting gas, potentially associated to a cooling flow that connects the hot ICM reservoir to the cold gas phase responsible for star formation and supermassive black hole feeding. We use new deep XMM-Newton EPIC and RGS data, combined with archival Chandra data, to perform a spectral analysis for the core region. We find 1$\sigma$ upper limits to the cooling gas fraction of $\sim$40 $\rm M_{\odot}yr^{-1}$ and $\sim$50-60 $\rm M_{\odot}yr^{-1}$ in the temperature ranges 0.5-1 keV and 1-2 keV, respectively. The lack of OVII, FeXXI-FeXXII, and FeXVII emission lines in the RGS spectra suggest that the fraction of gas cooling below 1 keV is limited to a few tens of $\rm M_{\odot}yr^{-1}$ at most. However, we detect several lines (e.g. SiXIV, MgXII, FeXXIII/FeXXIV, NeX, OVIII$\alpha$) that allow us to estimate a 1$\sigma$ upper limit for turbulent broadening of $\sim$320 km $\rm s^{-1}$, higher that other cool-core clusters such as Abell 1835, implying mechanisms that boost turbulence in Abell 2667's atmosphere. Imaging analysis of Chandra data suggests the presence of a cold front, possibly lined to sloshing or ICM cavities. However, current data do not clearly identify the physical mechanism driving turbulence. These finding indicate that Abell 2667 is similar to other low-redshift cool-core clusters, though the large upper limit on turbulence hints at significant ICM heating, which may suppress cooling for extended periods and contribute to future condensation events.

The Perseus cluster is the brightest X-ray cluster in the sky and is known as a cool-core galaxy cluster. Being a very nearby cluster, it has been extensively studied. This has provided a comprehensive view of the physical processes that operate in the intracluster medium (ICM), including feedback from the AGN 3C84 and measurements of ICM turbulence. Additionally, the Perseus cluster contains a central radio mini-halo. This diffuse radio source traces cosmic ray electrons (re-)accelerated in-situ in the ICM. Here we report on LOFAR high-band antenna 120-168 MHz observations of the Perseus cluster that probe a range of four orders of magnitude in angular scales. In our 0.3 arcsec resolution image, we find that the northern extension of the 3C84 lobe consists of several narrow 1.5-3 kpc parallel strands of emission. In addition, we detect steep-spectrum filaments associated with a previous outburst of the central AGN radio emission filling two known X-ray ghost cavities. At 7 arcsec resolution, our images show a complex structured radio mini-halo, with several edges and filaments. At resolutions of 26 arcsec and 80 arcsec, we discover diffuse radio emission with a 1.1 Mpc extent. We classify this emission as a giant radio halo and its properties are distinct from the inner mini-halo. We also detect two diffuse sources at projected cluster centric radii of 0.7 and 1.0 Mpc. Finally, we observe a 0.9 Mpc long trail of radio emission from the cluster member galaxy IC310, connecting it with the giant radio halo. Together with other recent studies of relaxed clusters, our LOFAR observations indicate that cluster-wide radio emission could be (more) common in cool-core clusters. In the case of the Perseus cluster, a past off-axis merger event that preserved the cool core might have generated enough turbulence to produce an extended radio halo observable at low frequencies., Comment: accepted for publication in A&A, 18 pages, 19 figures

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

Context. Tailed radio galaxies are shaped by ram pressure owing to the high-velocity motion of their host through the intracluster medium (ICM). Recent works have reported on the increasing complexity of the phenomenology of tailed galaxies, with departures from theoretical ageing models and evidence of re-energising mechanisms, which are yet unclear. Aims. The nearby (z = 0.0894) galaxy cluster Abell 2142 hosts two tailed galaxies, namely T1 and T2, which exhibit peculiar morphological features. We aim to investigate the properties of T1 and T2 and constrain their spectral evolution, dynamics, and interactions with the ICM. Methods. We combined LOw Frequency Array (LOFAR), upgraded Giant Metrewave Radio Telescope (uGMRT), Very Large Array (VLA), and MeerKAT data (from 30 MHz to 6.5 GHz) to carry out a detailed spectral analysis of T1 and T2. We analysed surface brightness profiles, measured integrated and spatially-resolved spectral indices, and performed a comparison with single injection ageing models. Chandra X-ray data were used to search for discontinuities in the ICM properties in the direction of the targets. Results. The spectral properties of T1 at low frequencies are predicted by ageing models, and provide constraints on the 3D dynamics of the host by assuming a constant velocity. However, sharp transitions along sub-regions of the tail, local surface brightness enhancements, and a spectral shape at high frequencies that is not predicted by models suggest a more complex scenario, possibly involving hydrodynamical instabilities and particle mixing. T2 exhibits unusual morphological and surface brightness features, and its spectral behaviour is not predicted by standard models. Two AGN outburst events during the infall of T2 towards the cluster centre could explain its properties., Comment: 23 pages (including Appendices). 15 figures. Accepted for publication in A&A

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

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

Abell 754 is a rich galaxy cluster at $z=0.0543$ and is considered the prototype of a major cluster merger. Like many dynamically unrelaxed systems, it hosts diffuse radio emission on Mpc-scales. Extended synchrotron sources in the intra-cluster medium (ICM) are commonly interpreted as evidence that a fraction of the gravitational energy released during cluster mergers is dissipated into nonthermal components. Here, we use new MeerKAT UHF- and L-band observations to study nonthermal phenomena in Abell 754. These data are complemented with archival XMM-Newton observations to investigate the resolved spectral properties of both the radio and X-ray cluster emission.For the first time, we employed the pipeline originally developed to calibrate LOFAR data to MeerKAT observations. This allowed us to perform a direction-dependent calibration and obtain highly sensitive radio images in UHF- and L-bands which capture the extended emission with unprecedented detail. By using a large XMM-Newton mosaic, we produced thermodynamic maps of the ICM. Our analysis reveals that the radio halo in the cluster center is bounded by the well-known shock in the eastern direction. Furthermore, in the southwest periphery, we discover an extended radio source that we classify as a radio relic which is possibly tracing a shock driven by the squeezed gas compressed by the merger, outflowing in perpendicular directions. The low-luminosity of this relic appears compatible with direct acceleration of thermal pool electrons. We interpret the observed radio and X-ray features in the context of a major cluster merger with a nonzero impact parameter. Abell 754 is a remarkable galaxy cluster showcasing exceptional features associated with the ongoing merger event. The high quality of the new MeerKAT data motivates further work on this system., Comment: 18 pages, 14 figures, 3 tables (excluding Appendix); Updated to match the accepted version in A&A

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

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-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 scaling relation between $Y_{\text{SZ}}$, the integrated Sunyaev-Zeldovich signal and $M_{500}^{Y_{\text{X}}}$, the cluster mass derived from X-ray observations, 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 obtained from XMM-Newton observations of a subsample of the ESZ. We calibrated a mass bias on a subset of the Planck cosmological cluster sample using published weak-lensing data from CCCP and MENeaCS, 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 Chandra-calibrated scaling relation, and $(1-b)=0.76\pm0.04$ for the XMM-Newton-calibrated scaling relation. We apply the scaling relations we derived to the full Planck cosmological cluster sample and obtain identical cosmological constraints regardless of the X-ray sample used, with $\sigma_8 =0.77\pm0.02$, $\Omega_m=0.31\pm0.02$, and $S_8= \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 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 relaxing the assumption of redshift evolution and fitting it to the data, we find no significant tension with results from either late-time probes or the CMB., Comment: 21 pages, 20 figures, accepted by 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.

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

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

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

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

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

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

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

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

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

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

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

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.

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