Search

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

We present deep, wideband multifrequency radio observations (144 MHz$-$8 GHz) of the remarkable galaxy group NGC 741, which yield crucial insights into the interaction between the infalling head-tail radio galaxy (NGC 742) and the main group. Our new data provide an unprecedentedly detailed view of the NGC 741-742 system, including the shock cone, disrupted jets from NGC 742, the long ($\sim$ 255 kpc) braided southern radio tail, and eastern lobe-like structure, and reveal, for the first time, complex radio filaments throughout the tail and lobe, and a likely vortex ring behind the shock cone. The cone traces the bow shock caused by the supersonic ($\mathcal{M}\sim2$) interaction between the head-tail radio galaxy NGC 742 and the intragroup medium (IGrM) while the ring may have been formed by interaction between the NGC 742 shock and a previously existing lobe associated with NGC 741. This interaction plausibly compressed and re-accelerated the radio plasma. We estimate that shock-heating by NGC 742 has likely injected $\sim$2-5$\times$10$^{57}$ erg of thermal energy into the central 10 kpc cooling region of the IGrM, potentially affecting the cooling and feedback cycle of NGC 741. A comparison with Chandra X-ray images shows that some of the previously detected thermal filaments align with radio edges, suggesting compression of the IGrM as the relativistic plasma of the NGC 742 tail interacts with the surrounding medium. Our results highlight that multi-frequency observations are key to disentangling the complex, intertwined origins of the variety of radio features seen in the galaxy group NGC 741, and the need for simulations to reproduce all the detected features., Comment: 26 pages, 15 figures, submitted to ApJ

In spite of their key role in galaxy evolution and several decades of observational efforts, the census of supernova remnants (SNRs) in our Galaxy remains incomplete. Theoretical predictions based on the local supernova rate estimate the expected number of SNRs in the Galaxy to be $\gtrsim$ 1000. By contrast, the number of detected SNRs amounts to about 300. High-resolution, wide-area radio surveys at low frequencies are ideal tools with which to find missing SNRs, given the prominence of these sources at low radio frequencies. We aim to find missing SNRs using proprietary data from the LOFAR Two-Metre Sky Survey (LoTSS) at 144~MHz. We used LoTSS total intensity maps of two Galactic regions, one with $39^\mathrm{o} < l < 66^\mathrm{o}$ and $|b|< 2.5^\mathrm{o}$, and the other with $145^\mathrm{o} < l < 150^\mathrm{o}$ and $|b| < 3^\mathrm{o}$, in addition to mid-infrared (MIR) data from the Wide-field Infrared Survey Explorer (WISE) all-sky survey to search for SNR candidates. We report the discovery of 14 new SNR candidates selected on the basis of their morphology at 144 MHz and a lack of MIR emission. We also follow up on 24 previously reported SNR candidates, inferring their spectral index between the LoTSS frequency (144 MHz) and the frequency at which they were reported. The high resolution and sensitivity of LoTSS observations has resulted in the detection of 14 new SNR candidates. In order to unambiguously confirm the SNR nature of these candidates, follow-up X-ray observations are required with facilities such as eROSITA., Comment: Accepted in 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

We report on the serendipitous discovery of an unprecedented interaction between the radio lobe of a radio galaxy and a spiral galaxy. The discovery was made thanks to LOFAR observations at 144 MHz of the galaxy cluster Abell 160 ($z=0.04317$) provided by the LOFAR Two-metre Sky Survey. The new low-frequency observations revealed that one of the radio plumes of the central galaxy GIN 049 overlaps with the spiral galaxy JO36. Previous studies carried out with MUSE revealed that the warm ionized gas in the disk of JO36, traced by the H$\alpha$ emission, is severely truncated with respect to the stellar disk. We further explore this unique system by including new uGMRT observations at 675 MHz to map the spectral index. The emerging scenario is that JO36 has interacted with the radio plume in the past 200-500 Myr. The encounter resulted in a positive feedback event for JO36 in the form of a star formation rate burst of $\sim14$ $M_\odot$ yr$^{-1}$. In turn, the galaxy passage left a trace in the radio-old plasma by re-shaping the old relativistic plasma via magnetic draping., Comment: 20 pages, 11 figures. Accepted for publication on ApJ on September 4th, 2023

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

The programming language Prolog makes declarative programming possible, at least to a substantial extent. Programs may be written and reasoned about in terms of their declarative semantics. All the advantages of declarative programming are however lost when it comes to program debugging. This is because the Prolog debugger is based solely on the operational semantics. Declarative methods of diagnosis (i.e. locating errors in programs) exist, but are neglected. This paper discusses their possibly main weaknesses and shows how to overcome them. We argue that useful ways of declarative diagnosis of logic programs exist, and should be usable in actual programming., Comment: In Proceedings ICLP 2023, arXiv:2308.14898

We discuss how to implement backjumping (or intelligent backtracking) in Prolog by using the built-ins throw/1 and catch/3. We show that it is impossible in a general case, contrary to a claim that ``backjumping is exception handling". We provide two solutions. One works for binary programs; in a general case it imposes a restriction on where backjumping may originate. The other restricts the class of backjump targets. We also discuss implementing backjumping by using backtracking and the Prolog database. Additionally, we explain the semantics of Prolog exception handling in the presence of coroutining., Comment: 14 pages, 2 figures. Accepted for ICLP 2023. A former version of this paper is arxiv:2008.09610

We present a novel natural language processing (NLP) approach to deriving plain English descriptors for science cases otherwise restricted by obfuscating technical terminology. We address the limitations of common radio galaxy morphology classifications by applying this approach. We experimentally derive a set of semantic tags for the Radio Galaxy Zoo EMU (Evolutionary Map of the Universe) project and the wider astronomical community. We collect 8,486 plain English annotations of radio galaxy morphology, from which we derive a taxonomy of tags. The tags are plain English. The result is an extensible framework which is more flexible, more easily communicated, and more sensitive to rare feature combinations which are indescribable using the current framework of radio astronomy classifications., Comment: 17 pages, 11 Figures, Accepted at MNRAS

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

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

Radio galaxies are luminous structures created by the jets of supermassive black holes, and consist of atomic nuclei, relativistic electrons, and magnetic fields. In exceptional cases, radio galaxies attain cosmological, megaparsec extents - and thus turn into giants. Giants embody the most extreme known mechanism through which galaxies can impact the Cosmic Web around them. The triggers of giant growth remain a mystery. Excitingly, new sensitive low-frequency sky surveys hold promise to change this situation. In this work, we perform a precision measurement of the distribution of giant growth's central dynamical quantity: total length. We first construct a statistical geometric framework for radio galaxies that is both rigorous and practical. We then search the LOFAR Two-metre Sky Survey for giants, discovering 2060 previously unknown specimina: more than have been found in all preceding literature combined. Spectacular discoveries include the longest giant hosted by an elliptical galaxy, the longest giant hosted by a spiral galaxy, and 13 giants with an angular length larger than that of the full Moon. By combining theory and observations - and carefully forward modelling selection effects - we infer that giant radio galaxy lengths are well described by a Pareto distribution with tail index $-3.5 \pm 0.5$. This finding is a new observational constraint for models and simulations of radio galaxy growth. In addition, for the first time, we determine the comoving number density of giants, $5 \pm 2\ (100\ \mathrm{Mpc})^{-3}$, and the volume-filling fraction of giant radio galaxy lobes in clusters and filaments, $5\substack{+8\\-2}\cdot 10^{-6}$. We conclude that giants are truly rare - not only from an observational perspective, but also from a cosmological one. At any moment in time, most clusters and filaments - the building blocks of the modern Cosmic Web - do not harbour giants., Comment: 41 pages, 28 figures, 5 tables, published in Astronomy & Astrophysics

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

High redshift quasars ($z>5$) that also shine brightly at radio wavelengths are unique signposts of supermassive black hole activity in the early universe. However, bright radio sources at $z\ge5$ are extremely rare and therefore we have started a campaign to search for new high-$z$ quasars by combining an optical dropout selection driven by the $g$, $r$, and $z$ bands from the Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Surveys with low-frequency radio observations from the LOFAR Two-metre Sky Survey (LoTSS). Currently, LoTSS covers a large fraction of the northern sky (5720 deg$^2$) to such a depth (median noise level of 83 $\mu$Jy beam$^{-1}$) that about 30% of the general quasar population is detected $-$ which is a factor of 5-10 more than previous large sky radio surveys such as NVSS and FIRST, respectively. In this paper, we present the discovery of 20 new quasars (and the independent confirmation of 4) between $4.9\leq z\leq 6.6$. Out of the 24 quasars, 21 satisfy the traditional radio-loudness criterion of $R=f_{5\text{GHz}}/f_{4400A} > 10$, with the full sample spanning $R\sim$6-1000, thereby more than doubling the sample of known radio-loud quasars at $z \ge 5$. Our radio detection requirement strongly decreases the contamination of stellar sources and allows one to select these quasars in a broad redshift range. Despite selecting our quasar candidates using fewer and less conservative colour restrictions, both the optical and near-infrared colours, Ly$\alpha$ emission line properties, and dust reddening, $E(B-V)$, measurements of our quasar sample do not deviate from the known radio-quiet quasar population, suggesting similar optical quasar properties of the radio-loud and radio-quiet quasar population at high-$z$. Our campaign demonstrates the potential for discovering new high-$z$ quasar populations through next generation radio continuum surveys., Comment: 19 pages, 11 figures, Tab A1 corrected, accepted for publication in A&A

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

Context. Magnetic fields are key to understand galaxy evolution, regulating stellar feedback and star formation in galaxies. Aims. We probe the origin of magnetic fields in late-type galaxies, measuring magnetic field strengths, exploring whether magnetic fields are only passive constituents of the interstellar medium, or whether they are active constituents being part of the local energy equilibrium. Methods. We measure equipartition magnetic field strengths in 39 galaxies from LoTSS-DR2 using LOFAR observations at 144 MHz with 6 arcsec angular resolution which (0.1-0.7 kpc). For a subset of 9 galaxies, we obtain atomic and molecular mass surface densities using HI and CO(2-1) data, from the THINGS and HERACLES surveys, respectively. These data are at 13 arcsec angular resolution, which corresponds to 0.3-1.2 kpc at the distances of our galaxies. We measure kinetic energy densities using HI and CO velocity dispersions. Results. We found a mean magnetic field strength of 3.6-12.5 $\mu$G with a mean of $7.9 \pm 2.0$ $\mu$G across the full sample. The magnetic field strength has the tightest and steepest relation with the total gas surface density with $B\propto \Sigma_{\rm HI+H2}^{0.309\pm0.006}$. The relation with the star-formation rate surface density and molecular gas surface density has significantly flatter slopes. After accounting for the influence of cosmic-ray transport, we found an even steeper relation of $B\propto \Sigma_{\rm HI+H2}^{0.393\pm0.009}$. Conclusions. These results suggest that the magnetic field is regulated by a $B$-$\rho$ relation, which has its origin in the saturation of the small-scale dynamo. This is borne out by an agreement of kinetic and magnetic energy densities although local deviations do exist in particular in areas of high kinetic energy densities where the magnetic field is sub-dominant., Comment: Accepted to Astronomy and Astrophysics. 13 pages plus appendix including atlas of magnetic fields

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

Broad Absorption Line Quasars (BALQSOs) show strong signatures of powerful outflows, with the potential to alter the cosmic history of their host galaxies. These signatures are only seen in ~10% of optically selected quasars, although the fraction significantly increases in IR and radio selected samples. A proven physical explanation for this observed fraction has yet to be found, along with a determination of why this fraction increases at radio wavelengths. We present the largest sample of radio matched BALQSOs using the LOFAR Two-metre Sky Survey Data Release 2 and employ it to investigate radio properties of BALQSOs. Within the DR2 footprint, there are 3537 BALQSOs from Sloan Digital Sky Survey DR12 with continuum signal to noise >5. We find radio-detections for 1108 BALQSOs, with an important sub-population of 120 LoBALs, an unprecedented sample size for radio matched BALQSOs given the LoTSS sky coverage to date. BALQSOs are a radio-quiet population that show an increase of $\times 1.50$ radio-detection fraction compared to non-BALQSOs. LoBALs show an increase of $\times 2.22$ that of non-BALQSO quasars. We show that this detection fraction correlates with wind-strength, reddening and C_{IV} emission properties of BALQSOs and that these features may be connected, although no single property can fully explain the enhanced radio detection fraction. We create composite spectra for sub-classes of BALQSOs based on wind strength and colour, finding differences in the absorption profiles of radio-detected and radio-undetected sources, particularly for LoBALs. Overall, we favour a wind-ISM interaction explanation for the increased radio-detection fraction of BALQSOs., Comment: 16 pages, 14 figures

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

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

The effect of Active Galactic Nuclei (AGN) on their host galaxies -- in particular their levels of star formation -- remains one of the key outstanding questions of galaxy evolution. Successful cosmological models of galaxy evolution require a fraction of energy released by an AGN to be redistributed into the interstellar medium to reproduce the observed stellar mass and luminosity function and to prevent the formation of over-massive galaxies. Observations have confirmed that the radio-AGN population is energetically capable of heating and redistributing gas at all phases, however, direct evidence of AGN enhancing or quenching star formation remains rare. With modern, deep radio surveys and large integral field spectroscopy (IFS) surveys, we can detect fainter synchrotron emission from AGN jets and accurately probe the star-forming properties of galaxies, respectively. In this paper, we combine data from the LOw Frequency ARray Two-meter Sky Survey with data from one of the largest optical IFS surveys, Mapping Nearby Galaxies at Apache Point Observatory to probe the star-forming properties of 307 local (z $<$ 0.15) galaxies that host radio-detected AGN (RDAGN). We compare our results to a robust control sample of non-active galaxies that each match the stellar mass, redshift, visual morphology, and inclination of a RDAGN host. We find that RDAGN and control galaxies have broad SFR distributions, typically lie below the star-forming main-sequence, and have negative stellar light-weighted age gradients. These results indicate that AGN selected based on their current activity are not responsible for suppressing their host galaxies' star formation. Rather, our results support the maintenance mode role that radio AGN are expected to have in the local Universe., Comment: Accepted for publication in A&A. 24 pages, 14 figures

Studies of the most distant AGNs allow us to test our current understanding of the physics present in radio-jetted AGNs across a range of environments. The decrease in apparent luminosity with distance is the primary difficulty to overcome in the study of these distant AGNs, which requires highly sensitive instruments. Our goal is to employ new long wavelength radio data to better parametrise the broad-band SED of GB 1508+5714, a high-redshift (z=4.30) AGN. Its high redshift, high intrinsic luminosity, and classification as a blazar allow us to test emission models that consider the efficient cooling of jet electrons via inverse Compton losses in interactions with the dense CMB photon field at high redshifts. A significant detection of this effect in GB 1508+5714 may partly explain the apparent sparsity of high-redshift radio galaxies in wide-field surveys; detections of this kind are only becoming possible with the current generation of SKA precursors. We used international LOFAR telescope to image the long wavelength radio emission around the high-redshift blazar GB 1508+5714 on arcsecond scales at frequencies between 128 MHz and 160 MHz. This allowed us to compare the spatially resolved structure with higher frequency observations, and to construct spectral index maps. The LOFAR image shows a compact unresolved core and two resolved emission regions around 2 arcsec to the east and to the west of the radio core. We find structure consistent with previous VLA observations, as well as a previously unreported emission region to the east. We interpret the arcsecond-scale radio structure of GB 1508+5714 as a FR II-like radio galaxy at a small viewing angle. Our SED modelling shows that a scenario featuring significant quenching effects caused by interaction with the CMB provides a good description of the data, and notably explains the suppressed radio emission., Comment: 11 pages, 10 figures, 2 tables

We weaken the notion of "not subject to occur-check" (NSTO), on which most known results on avoiding the occur-check in logic programming are based. NSTO means that unification is performed only on such pairs of atoms for which the occur-check never succeeds in any run of a nondeterministic unification algorithm. Here we show that "any run" can be weakened to "some run". We present some related sufficient conditions under which the occur-check may be safely omitted. We show examples for which the proposed approach provides more general results than the approaches based on well-moded and nicely moded programs (this includes cases to which the latter approaches are inapplicable). We additionally present a sufficient condition based on NSTO, working for arbitrary selection rules., Comment: 21 pages. Extended version of arXiv:2109.07914 in "Proceedings 37th ICLP (Technical Communications)". This version: corrections, extensions, simplifications

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

Observations of synchrotron emission at low radio frequencies reveal a labyrinth of polarised Galactic structures. However, the explanation for the wealth of structures remains uncertain due to the complex interactions between the interstellar medium and the magnetic field. A multi-tracer approach to the analysis of large sky areas is needed. This paper aims to use polarimetric images from the LOFAR Two metre Sky Survey (LoTSS) to produce the biggest mosaic of polarised emission in the northern sky at low radio frequencies (150 MHz) to date. The large area this mosaic covers allows for detailed morphological and statistical studies of polarised structures in the high-latitude outer Galaxy, including the well-known Loop III region. We produced a 3100 square degree Faraday tomographic cube using a rotation measure synthesis tool. We calculated the statistical moments of Faraday spectra and compared them with data sets at higher frequencies (1.4 GHz) and with a map of a rotation measure derived from extragalactic sources. The mosaic is dominated by polarised emission connected to Loop III. Additionally, the mosaic reveals an abundance of other morphological structures, mainly {narrow and extended} depolarisation canals, which are found to be ubiquitous. We find a correlation between the map of an extragalactic rotation measure and the LoTSS first Faraday moment image. The ratio of the two deviates from a simple model of a Burn slab (Burn 1966) along the line of sight, which highlights the high level of complexity in the magnetoionic medium that can be studied at these frequencies., Comment: 20 pages, 25 figures, accepted for publication in A&A

In this data release from the LOFAR Two-metre Sky Survey (LoTSS) we present 120-168MHz images covering 27% of the northern sky. Our coverage is split into two regions centred at approximately 12h45m +44$^\circ$30' and 1h00m +28$^\circ$00' and spanning 4178 and 1457 square degrees respectively. The images were derived from 3,451hrs (7.6PB) of LOFAR High Band Antenna data which were corrected for the direction-independent instrumental properties as well as direction-dependent ionospheric distortions during extensive, but fully automated, data processing. A catalogue of 4,396,228 radio sources is derived from our total intensity (Stokes I) maps, where the majority of these have never been detected at radio wavelengths before. At 6" resolution, our full bandwidth Stokes I continuum maps with a central frequency of 144MHz have: a median rms sensitivity of 83$\mu$Jy/beam; a flux density scale accuracy of approximately 10%; an astrometric accuracy of 0.2"; and we estimate the point-source completeness to be 90% at a peak brightness of 0.8mJy/beam. By creating three 16MHz bandwidth images across the band we are able to measure the in-band spectral index of many sources, albeit with an error on the derived spectral index of +/-0.2 which is a consequence of our flux-density scale accuracy and small fractional bandwidth. Our circular polarisation (Stokes V) 20" resolution 120-168MHz continuum images have a median rms sensitivity of 95$\mu$Jy/beam, and we estimate a Stokes I to Stokes V leakage of 0.056%. Our linear polarisation (Stokes Q and Stokes U) image cubes consist of 480 x 97.6 kHz wide planes and have a median rms sensitivity per plane of 10.8mJy/beam at 4' and 2.2mJy/beam at 20"; we estimate the Stokes I to Stokes Q/U leakage to be approximately 0.2%. Here we characterise and publicly release our Stokes I, Q, U and V images in addition to the calibrated uv-data., Comment: 23 figures, 1 table and 29 pages. The catalogues, images and uv-data associated with this data release are publicly available via https://lofar-surveys.org/

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

We study the properties and the origin of the radio emission in the most luminous early-type galaxies (ETGs) in the nearby Universe (MK<-25, recession velocity < 7,500 km/s) as seen by the 150 MHz Low-Frequency ARray (LOFAR) observations. LOFAR images are available for 188 of these giant ETGs (gETGs) and 146 (78%) of them are detected above a typical luminosity of ~10E21 W/Hz. They show a large spread in power, reaching up to ~10E26 W/Hz. We confirm a positive link between the stellar luminosity of gETGs and their median radio power, the detection rate, and the fraction of extended sources. About two-thirds (91) of the detected gETGs are unresolved, with sizes <4 kpc, confirming the prevalence of compact radio sources in local sources. Forty-six gETGs show extended emission on scales ranging from 4 to 340 kpc, at least 80% of which have a FRI class morphology. Based on the morphology and spectral index of the extended sources, ~30% of them might be remnant or restarted sources but further studies are needed to confirm this. Optical spectroscopy (available for 44 gETGs) indicates that for seven of them the nuclear gas is ionized by young stars suggesting a contribution to their radio emission from star forming regions. Their radio luminosities correspond to a star formation rate (SFR) in the range 0.1-8 Msun/yr and a median specific SFR of 0.8x10E-12 yr-1. The gas flowing toward the center of gETGs can accrete onto the supermassive black hole but also stall at larger radii and form new stars, an indication that feedback does not completely quench star formation. The most luminous gETGs (25 galaxies with MK < -25.8) are all detected at 150 MHz however they are not all currently turned on: at least four of them are remnant sources and at least one is likely powered by star formation., Comment: Pre-proofs version - Accepted for publication in A&A

The empirical relationship between the non-thermal 5GHz radio luminosity and the soft X-ray luminosity of active stellar coronae, canonically called the G\"udel-Benz relationship (G\"udel & Benz 1993), has been a cornerstone of stellar radio astronomy as it explicitly ties the radio emission to the coronal heating mechanisms. The relationship extends from microflares on the Sun to the coronae of the most active stars suggesting that active coronae are heated by a flare-like process (Benz & G\"udel 1994). The relationship is thought to originate from a consistent partition of the available flare energy into relativistic charges, that emit in the radio-band via the incoherent gyrosynchrotron mechanism, and heating of the bulk coronal plasma, that emits in the X-ray band via the Bremsstrahlung mechanism. Consequently, coherent emission from stellar and sub-stellar objects is not expected to adhere to this empirical relationship, as is observed in ultracool dwarf stars and brown dwarfs. Here we report a population of radio-detected chromospherically active stars that surprisingly follows the G\"udel-Benz relationship despite their radio emission being classified as coherent emission by virtue of its high circularly polarised fraction and high brightness temperature. Our results prompt a re-examination of the physics behind the G\"udel-Benz relationship, its implication for the mechanism of coronal heating and particle acceleration in active stars and the phenomenological connection between solar and stellar flares., Comment: Under review, ApJL

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

In a recent work, a large sample of double-peak (DP) emission line galaxies have been identified from the SDSS. While the two peaks could represent two kinematic components, they may be linked to the large bulges which their host galaxies tend to have. Star-forming DP galaxies display a central star formation enhancement and have been discussed to be compatible with a sequence of recent minor mergers. In order to probe merger induced star formation mechanisms, we conducted observations of the molecular gas content of 35 star-forming DP galaxies in the upper part of the main sequence (MS) of star formation with the IRAM 30m telescope. Including similar galaxies 0.3 dex above MS and with existing molecular gas observations from the literature, we finally obtain a sample of 52 such galaxies. We succeed in fitting the same kinematic parameters to the optical ionised and molecular gas emission lines for 10 (19%) galaxies. We find a central star formation enhancement resulting most likely from a galaxy merger or galaxy interaction, which is indicated by an excess of gas extinction found in the centre. This star formation is traced by radio continuum emissions at 150MHz, 1.4 and 3GHz, which are all three linearly correlated in log with the CO luminosity with the same slope. The 52 DP galaxies are found to have a significantly larger amount of molecular gas and larger depletion times, hence a lower star formation efficiency, than the expected values at their distance of the MS. This is consistent with a scenario of minor mergers increasing the mass of bulges and driving gas to the centre. We also exclude the inwards directed gas migration and central star formation enhancement to be the origin of a bar morphology. Hence, these 52 DP galaxies could be the results of recent minor mergers that funnelled molecular gas towards their centre, triggering star formation, but with a moderate efficiency., Comment: 33 pages, 19 figures, accepted by A&A

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

Optically luminous quasars at $z > 5$ are important probes of super-massive black hole (SMBH) formation. With new and future radio facilities, the discovery of the brightest low-frequency radio sources in this epoch would be an important new probe of cosmic reionization through 21-cm absorption experiments. In this work, we systematically study the low-frequency radio properties of a sample of 115 known spectroscopically confirmed $z>5$ quasars using the second data release of the Low Frequency Array (LOFAR) Two Metre Sky survey (LoTSS-DR2), reaching noise levels of $\sim$80 $\mu$Jy beam$^{-1}$ (at 144 MHz) over an area of $\sim5720$ deg$^2$. We find that 41 sources (36%) are detected in LoTSS-DR2 at $>2 \sigma$ significance and we explore the evolution of their radio properties (power, spectral index, and radio loudness) as a function of redshift and rest-frame ultra-violet properties. We obtain a median spectral index of $-0.29^{+0.10}_{-0.09}$ by stacking 93 quasars using LoTSS-DR2 and Faint Images of the Radio Sky at Twenty Centimetres (FIRST) data at 1.4 GHz, in line with observations of quasars at $z<3$. We compare the radio loudness of the high-$z$ quasar sample to a lower-$z$ quasar sample at $z\sim2$ and find that the two radio loudness distributions are consistent with no evolution, although the low number of high-z quasars means that we cannot rule out weak evolution. Furthermore, we make a first order empirical estimate of the $z=6$ quasar radio luminosity function, which is used to derive the expected number of high-$z$ sources that will be detected in the completed LoTSS survey. This work highlights the fact that new deep radio observations can be a valuable tool in selecting high-$z$ quasar candidates for follow-up spectroscopic observations by decreasing contamination of stellar dwarfs and reducing possible selection biases introduced by strict colour cuts., Comment: 11 pages, 10 figures, accepted for publication in A&A

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

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

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

Most known results on avoiding the occur-check are based on the notion of "not subject to occur-check" (NSTO). It means that unification is performed only on such pairs of atoms for which the occur-check never succeeds in any run of a nondeterministic unification algorithm. Here we show that this requirement is too strong. We show how to weaken it, and present some related sufficient conditions under which the occur-check may be safely omitted. We show examples for which the proposed approach provides more general results than the approaches based on well-moded and nicely moded programs (this includes cases to which the latter approaches are inapplicable)., Comment: In Proceedings ICLP 2021, arXiv:2109.07914

We present 150 MHz, 1.4 GHz, and 3 GHz radio imaging (LoTSS, FIRST and VLASS) and spatially resolved ionized gas characteristics (SDSS IV-MaNGA) for 140 local ($z<0.1$) early-type "red geyser" galaxies. These galaxies have low star formation activity (SFR $\sim \rm 0.01\ M_{\odot} yr^{-1}$), but show unique extended patterns in spatially-resolved emission line maps that have been interpreted as large-scale ionized winds driven by active galactic nuclei (AGN). In this work we confirm that red geysers host low-luminosity radio sources ($\rm L_{1.4GHz} \sim 10^{22} W Hz^{-1}$). Out of 42 radio-detected red geysers, 32 are spatially resolved in LoTSS and FIRST, with radio sizes varying between $\sim 5-25$ kpc. Three sources have radio sizes exceeding 40 kpc. A majority display a compact radio morphology and are consistent with either low-power compact radio sources ("FR0" galaxies) or "radio-quiet quasars". They may be powered by small-scale AGN-driven jets which remain unresolved at the current $5"$ resolution of radio data. The extended radio sources, not belonging to the "compact" morphological class, exhibit steeper spectra with a median spectral index of $-0.67$ indicating the dominance of lobed components. The red geysers hosting extended radio sources also have the lowest specific star formation rates, suggesting they either have a greater impact on the surrounding interstellar medium or are found in more massive halos on average. The degree of alignment of the ionized wind cone and the extended radio features are either 0$^{\circ}$ or 90$^{\circ}$, indicating possible interaction between the interstellar medium and the central radio AGN., Comment: 23 pages, 14 figures. Accepted for publication in The Astrophysical Journal

Thom Fr\"uhwirth presented a short, elegant and efficient Prolog program for the n queens problem. However the program may be seen as rather tricky and one may not be convinced about its correctness. This paper explains the program in a declarative way, and provides proofs of its correctness and completeness. The specification and the proofs are declarative, i.e. they abstract from any operational semantics. The specification is approximate, it is unnecessary to describe the program's semantics exactly. Despite the program works on non-ground terms, this work employs the standard semantics, based on logical consequence and Herbrand interpretations. Another purpose of the paper is to present an example of precise declarative reasoning about the semantics of a logic program. Under consideration in Theory and Practice of Logic Programming (TPLP)., Comment: Under consideration in Theory and Practice of Logic Programming (TPLP). 13 pages, 1 figure. A re-worked version of arXiv:1909.07479

[abridged] The International LOFAR Telescope is an interferometer with stations spread across Europe. With baselines of up to ~2,000 km, LOFAR has the unique capability of achieving sub-arcsecond resolution at frequencies below 200 MHz, although this is technically and logistically challenging. Here we present a calibration strategy that builds on previous high-resolution work with LOFAR. We give an overview of the calibration strategy and discuss the special challenges inherent to enacting high-resolution imaging with LOFAR, and describe the pipeline, which is publicly available, in detail. We demonstrate the calibration strategy by using the pipeline on P205+55, a typical LOFAR Two-metre Sky Survey (LoTSS) pointing. We perform in-field delay calibration, solution referencing to other calibrators, self-calibration, and imaging of example directions of interest in the field. For this specific field and these ionospheric conditions, dispersive delay solutions can be transferred between calibrators up to ~1.5 degrees away, while phase solution transferral works well over 1 degree. We demonstrate a check of the astrometry and flux density scale. Imaging in 17 directions, the restoring beam is typically 0.3" x 0.2" although this varies slightly over the entire 5 square degree field of view. We achieve ~80 to 300 $\mu$Jy/bm image rms noise, which is dependent on the distance from the phase centre; typical values are ~90 $\mu$Jy/bm for the 8 hour observation with 48 MHz of bandwidth. Seventy percent of processed sources are detected, and from this we estimate that we should be able to image ~900 sources per LoTSS pointing. This equates to ~3 million sources in the northern sky, which LoTSS will entirely cover in the next several years. Future optimisation of the calibration strategy for efficient post-processing of LoTSS at high resolution (LoTSS-HR) makes this estimate a lower limit., Comment: Accepted to a special issue of A&A on sub-arcsecond imaging with LOFAR. 24 pages, 16 figures

The Low-Frequency Array (LOFAR) Long-Baseline Calibrator Survey (LBCS) was conducted between 2014 and 2019 in order to obtain a set of suitable calibrators for the LOFAR array. In this paper we present the complete survey, building on the preliminary analysis published in 2016 which covered approximately half the survey area. The final catalogue consists of 30006 observations of 24713 sources in the northern sky, selected for a combination of high low-frequency radio flux density and flat spectral index using existing surveys (WENSS, NVSS, VLSS, and MSSS). Approximately one calibrator per square degree, suitable for calibration of $\geq$ 200 km baselines is identified by the detection of compact flux density, for declinations north of 30 degrees and away from the Galactic plane, with a considerably lower density south of this point due to relative difficulty in selecting flat-spectrum candidate sources in this area of the sky. Use of the VLBA calibrator list, together with statistical arguments by comparison with flux densities from lower-resolution catalogues, allow us to establish a rough flux density scale for the LBCS observations, so that LBCS statistics can be used to estimate compact flux densities on scales between 300 mas and 2 arcsec, for sources observed in the survey. The LBCS can be used to assess the structures of point sources in lower-resolution surveys, with significant reductions in the degree of coherence in these sources on scales between 2 arcsec and 300 mas. The LBCS survey sources show a greater incidence of compact flux density in quasars than in radio galaxies, consistent with unified schemes of radio sources. Comparison with samples of sources from interplanetary scintillation (IPS) studies with the Murchison Widefield Array (MWA) shows consistent patterns of detection of compact structure in sources observed both interferometrically with LOFAR and using IPS., Comment: Accepted to a special issue of A&A on sub-arcsecond imaging with LOFAR

Coherent radio emission from stars can be used to constrain fundamental coronal plasma parameters, such as plasma density and magnetic field strength. It is also a probe of chromospheric and magnetospheric acceleration mechanisms. Close stellar binaries, such as RS Canum Venaticorum (RS CVn) systems, are particularly interesting as their heightened level of chromospheric activity and possible direct magnetic interaction make them a unique laboratory to study coronal and magnetospheric acceleration mechanisms. RS CVn binaries are known to be radio-bright but coherent radio emission has only conclusively been detected previously in one system. Here, we present a population of 14 coherent radio emitting RS CVn systems. We identified the population in the ongoing LOFAR Two Metre Sky Survey as circularly polarised sources at 144MHz that are astrometrically associated with known RS CVn binaries. We show that the observed emission is powered by the electron cyclotron maser instability. We use numerical calculations of the maser's beaming geometry to argue that the commonly invoked 'loss-cone' maser cannot generate the necessary brightness temperature in some of our detections and that a more efficient instability, such as the shell or horseshoe maser, must be invoked. Such maser configurations are known to operate in planetary magnetospheres. We also outline two acceleration mechanisms that could produce coherent radio emission, one where the acceleration occurs in the chromosphere and one where the acceleration is due to an electrodynamic interaction between the stars. We propose radio and optical monitoring observations that can differentiate between these two mechanisms., Comment: This article consists of 14 pages, 10 normal figures and 4 gifs. The gifs themselves each consist of 32 figures. This article is to be published in Astronomy and Astrophysics, section 8. 'Stellar Atmospheres'

We present observations of planetary nebulae with the LOw Frequency ARray (LOFAR) between 120 and 168 MHz. The images show thermal free-free emission from the nebular shells. We have determined the electron temperatures for spatially resolved, optically thick nebulae. These temperatures are 20 to 60% lower than those estimated from collisionally excited optical emission lines. This strongly supports the existence of a cold plasma component, which co-exists with hot plasma in planetary nebulae. This cold plasma does not contribute to the collisionally excited lines, but does contribute to recombination lines and radio flux. Neither of the plasma components are spatially resolved in our images, although we infer that the cold plasma extends to the outer radii of planetary nebulae. However, more cold plasma appears to exist at smaller radii. The presence of cold plasma should be taken into account in modeling of radio emission of planetary nebulae. Modelling of radio emission usually uses electron temperatures calculated from collisionally excited optical and/or infrared lines. This may lead to an underestimate of the ionized mass and an overestimate of the extinction correction from planetary nebulae when derived from the radio flux alone. The correction improves the consistency of extinction derived from the radio fluxes when compared to estimates from the Balmer decrement flux ratios., Comment: 20 pages, 10 figures, ApJ, accepted

We present a search for transient radio sources on timescales of 2 to 9 yr at 150 MHz. This search is conducted by comparing the first Alternative Data Release of the TIFR GMRT Sky Survey (TGSS ADR1) and the second data release of the LOFAR Two-metre Sky Survey (LoTSS DR2). The overlapping survey area covers 5570 $\rm{deg}^2$ on the sky, or 14 per cent of the total sky. We introduce a method to compare the source catalogues that involves a pair match of sources, a flux density cutoff to meet the survey completeness limit and a newly developed compactness criterion. This method is used to identify both transient candidates in the TGSS source catalogue that have no counterpart in the LoTSS catalogue and transient candidates in LoTSS without a counterpart in TGSS. We find that imaging artefacts and uncertainties and variations in the flux density scales complicate the transient search. Our method to search for transients by comparing two different surveys, while taking into account imaging artefacts around bright sources and misaligned flux scales between surveys, is universally applicable to future radio transient searches. No transient sources were identified, but we are able to place an upper limit on the transient surface density of $<5.4 \cdot 10^{-4}\ \text{deg}^{-2}$ at 150 MHz for compact sources with an integrated flux density over 100 mJy. Here we define a transient as a compact source with flux density greater than 100 mJy that appears in the catalogue of one survey without a counterpart in the other survey., Comment: 14 pages, 11 figures