Search

Aims. The orientation-based unification scheme of radio-loud active galactic nuclei (AGNs) asserts that radio galaxies and quasars are essentially the same type of object, but viewed from different angles. To test this unification model, we compared the environments of radio galaxies and quasars, which would reveal similar properties when an accurate model is utilized. Methods. Using the second data release of the LOFAR Two-metre Sky Survey (LoTSS DR2), we constructed a sample of 26,577 radio galaxies and 2028 quasars at 0.08 < z < 0.4. For radio galaxies with optical spectra, we further classified them as 3631 low-excitation radio galaxies (LERGs) and 1143 high-excitation radio galaxies (HERGs). We crossmatched these samples with two galaxy cluster catalogs from the Sloan Digital Sky Survey (SDSS). Results. We find that $17.1 \pm 0.2%$ of the radio galaxies and $4.1 \pm 0.4%$ of the quasars are associated with galaxy clusters. Luminous quasars are very rare in clusters, while $18.7 \pm 0.7%$ LERGs and $15.2 \pm 1.1%$ HERGs reside in clusters. We also note that in radio galaxies, both HERGs and LERGs tend to reside in the centers of clusters, while quasars do not show a strong preference for their positions in clusters. Conclusions. This study shows that local quasars and radio galaxies exist in different environments, challenging the orientation-based unification model. This means that factors other than orientation may play an important role in distinguishing radio galaxies from quasars. The future WEAVE-LOFAR survey will offer high-quality spectroscopic data for a large number of radio sources and allow for a more comprehensive exploration of the environments of radio galaxies and quasars., Comment: 9 pages, 6 figures. Accepted for publication by A&A

Diffuse radio emission in galaxy clusters is a tracer of ultra-relativistic particles and $\mu$G-level magnetic fields, and is thought to be triggered by cluster merger events. In the distant Universe (i.e. $z>0.6$), such sources have been observed only in a handful of systems, and their study is important to understand the evolution of large-scale magnetic fields over the cosmic time. Previous studies of nine {\it Planck} clusters up to $z\sim0.9$ suggest a fast amplification of cluster-scale magnetic fields, at least up to half of the current Universe's age, and steep spectrum cluster scale emission, in line with particle re-acceleration due to turbulence. In this paper, we investigate the presence of diffuse radio emission in a larger sample of galaxy clusters reaching even higher redshifts (i.e. $z\gtrsim1$). We selected clusters from the Massive and Distant Clusters of {\it WISE} Survey (MaDCoWS) with richness $\lambda_{15}>40$ covering the area of the second data release of the LOFAR Two-Meter Sky Survey (LoTSS-DR2) at 144 MHz. These selected clusters are in the redshift range $0.78-1.53$ (with a median value of 1.05). We detect the possible presence of diffuse radio emission, with the largest linear sizes of $350-500$ kpc, in 5 out of the 56 clusters in our sample. If this diffuse radio emission is due to a radio halo, these radio sources lie on or above the scatter of the $P_\nu-M_{500}$ radio halo correlations (at 150 MHz and 1.4 GHz) found at $z<0.6$, depending on the mass assumed. We also find that these radio sources are at the limit of the detection by LoTSS, and therefore deeper observations will be important for future studies., Comment: 20 pages, 9 figures; accepted for publication in Astronomy & Astrophysics

The thermal and non-thermal components in galaxy clusters have properties that, although shaped from different physical phenomena, can share some similarities, mainly driven by their halo mass and the accretion processes. Scaling relations have been proven to exist for both components and studied in X-ray (thermal) and radio (non-thermal) bands. At the radio wavelength, such investigations are so far limited to the integrated quantities (e.g. total power and mass). We aimed to investigate the scaling relations between the mass of a galaxy cluster and its radio emission at low frequencies, treating both the integrated and the spatially resolved quantities for a sample of well-selected targets. We crossmatched LoTSS DR2 and CHEX-MATE datasets in order to get the deepest and most homogeneous radio data of a representative sample of objects. We analytically derived the expected relation between the radio power ($P_{\nu}$) and radio surface brightness profile, and performed a comparison with observational results. We obtained that properly accounting for the mass and redshift dependence in the radio profile can reduce the overall scatter by a factor of $\sim 4$, with an evident residual dependence on the cluster dynamical status. We showed that assuming no relation between the halo size ($R_{H}$) and the cluster mass ($M$) allowed us to reconcile the observed radio profile mass scaling and the one predicted starting from the $P_{\nu}-M$ relation. We discuss the implications of a lack of $R_H-M$ relation, assessing possible systematics and biases in the analyses, and interpreting it as a natural consequence of the structure formation process. Finally, we also considered the role of the magnetic field in the $P_{\nu}-M$ relation, putting constraints on its dependence upon the cluster mass and finding consistent results with expectations from our radio power mass scaling., Comment: 15 pages, 10 figures, Accepted by A&A

Nest200047 is a clear example of multiple radio bubbles from an Active Galactic Nucleus (AGN) in a galaxy group, featuring non-thermal filaments likely shaped by buoyancy, gas motions, and stabilized by magnetic fields. This study presents high-quality data obtained from uGMRT, MeerKAT, and VLA, alongside existing LOFAR data, to analyze the system's morphology and spectrum over a broad frequency range (53-1518 MHz). Our findings reveal new filamentary emission in the inner 60 kpc, surrounding and extending from the inner bubbles and jets, suggesting complex dynamical evolution of the non-thermal plasma in the group core. The filaments have widths of a few kpc and lengths from tens to hundreds of kpc, with a steep and curved radio spectrum ($\rm \alpha=1\sim2$). They exhibit a constant spectral index profile along their length, implying particles are either (re-)accelerated together or move at super-Alfvenic speeds. Spectral aging analysis yields jet active times between 50 and 100 Myr with short inactive phases, suggesting continuous energy injection typical of AGN feedback in galaxy groups. This study highlights the potential of combining high-quality radio data to understand recurrent jet activity and feedback, with implications for future research with the SKA observatory., Comment: 25 pages, 17 figures, accepted by A&A

Diffuse radio sources, known as mini-halos and halos, are detected at the centres of galaxy clusters. These centralized diffuse sources are typically observed individually, with both appearing together only in rare cases. The origin of the diffuse radio sources in such systems remains unclear. We investigate the formation of large-scale radio emission in the most X-ray luminous, massive galaxy cluster RXJ~1347.5-1145 which is known to host a mini-halo at its centre and possibly additional more extended emission. We conduct deep multi-frequency observations of the galaxy cluster using the MeerKAT at 1.28 GHz and the uGMRT (upgraded Giant Metrewave Radio Telescope) at 1.26 GHz and 700 MHz. We characterize the brightness and spectral properties of the central diffuse sources and combine our radio observations with \textit{Chandra} X-ray data to explore the correlation between the cluster's non-thermal and thermal emissions. We confirm the presence of the diffuse emission and find that it extends up to 1~Mpc in size. Our multi-wavelength data reveal that the central diffuse emission consists of two distinct components: a mini-halo located in the cluster core and a larger radio halo extending around it. The correlation between radio and X-ray surface brightness in both sources indicates a strong connection between the non-thermal and thermal properties of the ICM. The differing slopes in the $I_R-I_X$ and $\alpha-I_X$ relations suggest that distinct mechanisms are responsible for the formation of the mini-halo and halo. The properties of the halo align with the turbulent model, while both turbulent and hadronic processes may contribute to the formation of the mini-halo.

Head-tail radio galaxies are characterized by a head, corresponding to an elliptical galaxy, and two radio jets sweeping back from the head, forming an extended structure behind the host galaxy that is moving through the intracluster medium (ICM). This morphology arises from the interaction between the diffuse radio-emitting plasma and the surrounding environment. Sometimes revived fossil plasma is found in galaxy clusters, tracing old active galactic nucleus ejecta with a very steep spectrum re-energized through processes in the ICM, unrelated to the progenitor galaxy. We aim to study the central region of Abell 1775, a galaxy cluster in an unclear dynamical state at z = 0.072. It hosts two giant radio-loud elliptical galaxies, the head-tail radio galaxy that "breaks" at the position of a cold front detected in the X-rays, filamentary revived fossil plasma, and central diffuse emission. This study aims to investigate and constrain the spectral properties and trends along the head-tail, as well as the revived fossil plasma, to better understand the formation process of the non-thermal phenomena in A1775. We make use of LOFAR (144 MHz), and new deep uGMRT observations (400 and 650 MHz). We observe an overall steepening along the tail of the head-tail radio galaxy. In the radio colour-colour diagram, ageing models reproduce the emission of the head-tail. An unexpected brightness increase at the head of the tail suggests a complex bending of the jets. We derived the equipartition magnetic field and minimum pressure along the tail. We recovered the structure of the revived fossil plasma, which appears as thin filaments with ultra-steep spectra. We show that high-sensitivity, high-resolution observations at low frequencies are essential for detecting the full extent of the tail, enabling a deeper spectral analysis and resolving the structure and spectral properties of revived fossil plasma., Comment: 14 pages, 14 figures (including Appendices). Accepted for publication in A&A on 30 Jan 2025. Abstract abridged for arXiv submission

Producing ultra-deep high-angular-resolution images with current and next-generation radio interferometers introduces significant computational challenges. In particular, the imaging is so demanding that processing large datasets, accumulated over hundreds of hours on the same pointing, is likely infeasible in the current data reduction schemes. In this paper, we revisit a solution to this problem that was considered in the past but is not being used in modern software: sidereal visibility averaging (SVA). This technique combines individual observations taken at different sidereal days into one much smaller dataset by averaging visibilities at similar baseline coordinates. We present our method and validated it using four separate 8-hour observations of the ELAIS-N1 deep field, taken with the International LOw Frequency ARray (LOFAR) Telescope (ILT) at 140~MHz. Additionally, we assessed the accuracy constraints imposed by Earth's orbital motion relative to the observed pointing when combining multiple datasets. We find, with four observations, data volume reductions of a factor of 1.8 and computational time improvements of a factor of 1.6 compared to standard imaging. These factors will increase when more observations are combined with SVA. For instance, with 3000~hours of LOFAR data aimed at achieving sensitivities of the order of {\mu}Jy/beam at sub-arcsecond resolutions, we estimate data volume reductions of up to a factor of 169 and a 14-fold decrease in computing time using our current algorithm. This advancement for imaging large deep interferometric datasets will benefit current generation instruments, such as LOFAR, and upcoming instruments such as the Square Kilometre Array (SKA), provided the calibrated visibility data of the individual observations are retained.

We present the final 6'' resolution data release of the ELAIS-N1 field from the LOw-Frequency ARray (LOFAR) Two-metre Sky Survey Deep Fields project (LoTSS Deep). The 144MHz images are the most sensitive achieved to date at this frequency and were created from 290 TB of data obtained from 505 hrs on-source observations taken over 7.5 years. The data were processed following the strategies developed for previous LoTSS and LoTSS Deep data releases. The resulting images span 24.53 square degrees and, using a refined source detection approach, we identified 154,952 radio sources formed from 182,184 Gaussian components within this area. The maps reach a noise level of 10.7 $\mu$Jy/beam at 6'' resolution where approximately half of the noise is due to source confusion. In about 7.4% of the image our limited dynamic range around bright sources results in a further > 5% increase in the noise. The images have a flux density scale accuracy of about 9% and the standard deviation of offsets between our source positions and those from Pan-STARRS is 0.2'' in RA and Dec for high significance detections. We searched individual epoch images for variable sources, identifying 39 objects with considerable variation. We also searched for circularly polarised sources achieving three detections of previously known emitters (two stars and one pulsar) whilst constraining the typical polarisation fraction plus leakage to be less than 0.045%., Comment: Accepted for publication in A&A. 16 figures, 1 table and 20 pages. The catalogues and images associated with this data release are publicly available via https://lofar-surveys.org/

Some galaxy clusters contain non-thermal synchrotron emitting plasma permeating the intracluster medium (ICM). The spectral properties of this radio emission are not well characterized at decameter wavelengths ({\nu} < 30 MHz), primarily due to the severe corrupting effects of the ionosphere. Using a recently developed calibration strategy, we present LOFAR images below 30 MHz of the low mass galaxy cluster Abell 655, which was serendipitously detected in an observation of the bright calibrator 3C 196. We combine this observation with LOFAR data at 144 MHz, and new Band 4 Giant Metrewave Radio Telescope observations centered at 650 MHz. In the 15-30 MHz LOFAR image, diffuse emission is seen with a physical extent of about 700 kpc. We argue that the diffuse emission detected in this galaxy cluster likely has multiple origins. At higher frequencies (650 MHz), the diffuse emission resembles a radio halo, while at lower frequencies the emission seems to consist of several components and bar-like structures. It suggests that most low-frequency emission in this cluster comes from re-energized fossil plasma from old AGN outbursts, coexisting with the radio halo component. By counting the number of cluster radio detections in the decameter band, we estimate that around a quarter of the Planck clusters host re-energised fossil plasma that is detectable in the decameter band with LOFAR., Comment: Accepted for publication in A&A. 11 pages, 3 tables, 8 figures

Galaxy clusters and groups are the last link in the chain of hierarchical structure formation. Their environments can be significantly affected by outbursts from AGN, especially in groups where the medium density is lower and the gravitational potential shallower. The interaction between AGN and group weather can therefore greatly impact their evolution. We investigate the non-thermal radio emission in Abell 1213, a galaxy group which is part of a larger sample of ~50 systems (X-GAP) recently granted XMM-Newton observations. We exploit proprietary LOFAR 54 MHz and uGMRT 380 MHz observations, complementing them with 144 MHz LOFAR survey and XMM-Newton archival data. A1213 hosts a bright AGN associated with one of the central members, 4C 29.41, which was previously optically identified as a dumb-bell galaxy. Observations at 144 MHz at a resolution of 0.3'' allow us to resolve the central radio galaxy. From this source, a ~500 kpc-long tail extends North-East. Our analysis suggests that the tail likely originated from a past outburst of 4C 29.41, and its current state might be the result of the interaction with the surrounding environment. The plateau of the spectral index distribution in the Easternmost part of the tail suggests mild particle re-acceleration, that could have re-energised seed electrons from the past activity of the AGN. While we observe a spatial and physical correlation of the extended, central emission with the thermal plasma, which might hint at a mini-halo, current evidence cannot conclusively prove this. A1213 is only the first group, among the X-GAP sample, that we are able to investigate through low-frequency radio observations. Its complex environment once again demonstrates the significant impact that the interplay between thermal and non-thermal processes can have on galaxy groups., Comment: 15 pages, 10 figures

Both star formation (SF) and Active Galactic Nuclei (AGN) play an important role in galaxy evolution. Statistically quantifying their relative importance can be done using radio luminosity functions. Until now these relied on galaxy classifications, where sources with a mixture of radio emission from SF and AGN are labelled as either a star-forming galaxy or an AGN. This can cause the misestimation of the relevance of AGN. Brightness temperature measurements at 144 MHz with the International LOFAR telescope can separate radio emission from AGN and SF. We use the combination of sub-arcsec and arcsec resolution imaging of 7,497 sources in the Lockman Hole and ELAIS-N1 fields to identify AGN components in the sub-arcsec resolution images and subtract them from the total flux density, leaving flux density from SF only. We construct, for the first time, radio luminosity functions by physical process, either SF or AGN activity, revealing a hidden AGN population at $L_{\textrm{144MHz}}$$<10^{24}$ W$\,$Hz$^{-1}$ . This population is 1.56$\pm$0.06 more than expected for $0.5

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