We present the source associations, cross-identifications, and multi-wavelength properties of the faint radio source population detected in the deep tier of the LOFAR Two Metre Sky Survey (LoTSS): ...the LoTSS Deep Fields. The first LoTSS Deep Fields data release consists of deep radio imaging at 150 MHz of the ELAIS-N1, Lockman Hole, and Boötes fields, down to RMS sensitives of around 20, 22, and 32
μ
Jy beam
−1
, respectively. These fields are some of the best studied extra-galactic fields in the northern sky, with existing deep, wide-area panchromatic photometry from X-ray to infrared wavelengths, covering a total of ≈26 deg
2
. We first generated improved multi-wavelength catalogues in ELAIS-N1 and Lockman Hole; combined with the existing catalogue for Boötes, we present forced, matched aperture photometry for over 7.2 million sources across the three fields. We identified multi-wavelength counterparts to the radio detected sources, using a combination of the Likelihood Ratio method and visual classification, which greatly enhances the scientific potential of radio surveys and allows for the characterisation of the photometric redshifts and the physical properties of the host galaxies. The final radio-optical cross-matched catalogue consists of 81 951 radio-detected sources, with counterparts identified and multi-wavelength properties presented for 79 820 (>97%) sources. We also examine the properties of the host galaxies, and through stacking analysis find that the radio population with no identified counterpart is likely dominated by active galactic nuclei (AGN) at
z
~ 3−4. This dataset contains one of the largest samples of radio-selected star-forming galaxies and AGN at these depths, making it ideal for studying the history of star-formation, and the evolution of galaxies and AGN across cosmic time.
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
≥ 5 ...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 ∼83 μ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 four) between 4.9 ≤
z
≤ 6.6. Out of the 24 quasars, 21 satisfy the traditional radio-loudness criterion of
R
=
f
5 GHz
/
f
4400 Å
> 10, with the full sample spanning
R
∼ 6–1000, thereby more than doubling the sample of known radio-loud quasars at
z
≥ 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
α
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.
The cosmic radio dipole is an anisotropy in the number counts of radio sources and is analogous to the dipole seen in the cosmic microwave background (CMB). Measurements of source counts of large ...radio surveys have shown that, although the radio dipole is generally consistent in direction with the CMB dipole, the amplitudes are in tension. These observations present an intriguing puzzle, namely the cause of this discrepancy, with a true anisotropy breaking with the assumptions of the cosmological principle, invalidating the most common cosmological models that are built on these assumptions. We present a novel set of Bayesian estimators to determine the cosmic radio dipole and compare the results with those of commonly used methods applied to the Rapid ASKAP Continuum Survey (RACS) and the NRAO VLA Sky Survey (NVSS) radio surveys. In addition, we adapt the Bayesian estimators to take into account systematic effects known to influence large radio surveys of this kind, folding information such as the local noise floor or array configuration directly into the parameter estimation. The enhancement of these estimators allows us to greatly increase the number of sources used in the parameter estimation, yielding tighter constraints on the cosmic radio dipole estimation than previously achieved with NVSS and RACS. We extend the estimators further to work on multiple catalogues simultaneously, leading to a combined parameter estimation using both NVSS and RACS. The result is a dipole estimate that perfectly aligns with the CMB dipole in terms of direction but with an amplitude that is three times as large, and a significance of 4.8
σ
. This new dipole measurement is made to an unprecedented level of precision for radio sources, which is only matched by recent results using infrared quasars.
We present a novel method for identifying candidate high-redshift quasars (HzQs;
z
≳ 5.5) –which are unique probes of supermassive black hole growth in the early Universe– from large-area optical ...and infrared photometric surveys. Using Gaussian mixture models to construct likelihoods and incorporating informed priors based on population statistics, our method uses a Bayesian framework to assign posterior probabilities that differentiate between HzQs and contaminating sources. We additionally include deep radio data to obtain informed priors. Using existing HzQ data in the literature, we set a posterior threshold that accepts ∼90% of known HzQs while rejecting > 99% of contaminants such as dwarf stars or lower redshift galaxies. Running the probability selection on test samples of simulated HzQs and contaminants, we find that the efficacy of the probability method is higher than traditional colour cuts, decreasing the fraction of accepted contaminants by 86% while retaining a similar fraction of HzQs. As a test, we apply our method to the Pan-STARRS Data Release 1 (PS1) source catalogue within the HETDEX Spring field area on the sky, covering 400 sq. deg. and coinciding with deep radio data from the LOFAR Two-metre Sky Survey Data Release 1. From an initial sample of ∼5 × 10
5
sources in PS1, our selection shortlists 251 candidate HzQs, which are further reduced to 63 after visual inspection. Shallow spectroscopic follow-up of 13 high-probability HzQs resulted in the confirmation of a previously undiscovered quasar at
z
= 5.66 with photometric colours
i
−
z
= 1.4, lying outside the typically probed regions when selecting HzQs based on colours. This discovery demonstrates the efficacy of our probabilistic HzQ selection method in selecting more complete HzQ samples, which holds promise when employed on large existing and upcoming photometric data sets.
ABSTRACT
We present MeerKAT Absorption Line Survey (MALS) observations of the H i gas in the Klemola 31 galaxy group (z = 0.029), located along the line of sight to the radio-loud quasar PKS 2020−370 ...(z = 1.048). Four galaxies of the group are detected in H i emission, and H i absorption is also detected in front of PKS 2020−370 in Klemola 31A. The emission and absorption are somewhat compensating on the line of sight of the quasar, and the derived column density of the absorption appears underestimated, with respect to the neighbouring emission. A symmetric tilted-ring model of Klemola 31A, assuming the absorbing gas in regular rotation in the plane, yields a rather high spin temperature of 530 K. An alternative interpretation is that the absorbing gas is extra-planar, which will also account for its non-circular motion. The Na i/Ca ii ratio also suggests that the absorbing gas is unrelated to cold H i disc. Two of the galaxies in the Klemola group are interacting with a small companion, and reveal typical tidal tails, and velocity perturbations. Only one of the galaxies, ESO 400−13, reveals a strong H i deficiency, and a characteristic ram-pressure stripping, with a total asymmetry in the distribution of its gas. Since a small galaxy group as Klemola 31 is not expected to host a dense intra-group gas, this galaxy must be crossing the group at a very high velocity, mostly in the sky plane.
The BL Lac object PKS 1413+135 was observed by the Large Survey Project MeerKAT Absorption Line Survey (MALS) in the
L
-band, at 1139 MHz and 1293–1379 MHz, targeting the H
I
and OH lines in ...absorption at
z
= 0.24671. The radio continuum might come from the nucleus of the absorbing galaxy or from a background object at redshift lower than 0.5, as suggested by the absence of gravitational images. The H
I
absorption line is detected at a high signal-to-noise ratio, with a narrow central component, and with a red wing, confirming previous results. The OH 1720 MHz line is clearly detected in (maser) emission, peaking at a velocity shifted by −10 to −15 km s
−1
with respect to the H
I
peak. The 1612 MHz line is lost due to radio frequency interference. The OH 1667 MHz main line is tentatively detected in absorption, but not the 1665 MHz line. Over 30 years a high variability is observed in optical depths, due to the rapid changes of the line of sight caused by the superluminal motions of the radio knots. The H
I
line has varied by 20% in depth, while the OH-1720 MHz depth has varied by a factor of ∼3. The position of the central velocity and the widths also varied. The absorbing galaxy is an early-type spiral (maybe S0) seen edge-on, with a prominent dust lane, covering the whole disk. Given the measured mass concentration and the radio continuum size at centimeter wavelengths (100 mas corresponding to 400 pc at
z
= 0.25), the width of the absorption lines from the nuclear regions are expected up to 250 km s
−1
. The narrowness of the observed lines (< 15 km s
−1
) suggests that the absorption comes from an outer gas ring, as frequently observed in S0 galaxies. The millimetric lines are even narrower (< 1 km s
−1
), which corresponds to the continuum size restricted to the core. The radio core is covered by individual 1 pc molecular clouds, whose column density is a few 10
22
cm
−2
, which is compatible with the gas screen detected in X-rays.
The number counts of homogeneous samples of radio sources are a tried and true method of probing the large-scale structure of the Universe, as most radio sources outside the Galactic plane are at ...cosmological distances. As such, they are expected to trace the cosmic radio dipole, an anisotropy analogous to the dipole seen in the cosmic microwave background (CMB). Results have shown that although the cosmic radio dipole matches the direction of the CMB dipole, it has a significantly larger amplitude. This unexplained result challenges our assumption of the Universe being isotropic, which can have large repercussions for the current cosmological paradigm. Though significant measurements have been made, sensitivity to the radio dipole is generally hampered by systematic effects that can cause large biases in the measurement. Here we assess these systematics with data from the MeerKAT Absorption Line Survey (MALS), a blind search for absorption lines with pointings centred on bright radio sources. With the sensitivity and field of view of MeerKAT, thousands of sources are observed in each pointing, allowing for the possibility of measuring the cosmic radio dipole given enough pointings. We present the analysis of ten MALS pointings, focusing on systematic effects that could lead to an inhomogeneous catalogue. We describe the calibration and creation of full band continuum images and catalogues, producing a combined catalogue containing 16 307 sources and covering 37.5 square degrees of sky down to a sensitivity of 10 μJy beam
−1
. We measure the completeness, purity, and flux recovery statistics for these catalogues using simulated data. We investigate different source populations in the catalogues by looking at flux densities and spectral indices and how they might influence source counts. Using the noise characteristics of the pointings, we find global measures that can be used to correct for the incompleteness of the catalogue, producing corrected number counts down to 100–200 μJy. We show that we can homogenise the catalogues and properly account for systematic effects. We determine that we can measure the dipole to 3
σ
significance with 100 MALS pointings.
We present the serendipitous discovery of eight distant (> 50 pc) late M dwarfs with plausible associated radio emission at 144 MHz. The M-dwarf nature of our sources has been confirmed with optical ...spectroscopy performed using HET/LRS2 and Subaru/FOCAS, and their radio flux densities are within the range of 0.5−1.0 mJy at 144 MHz. Considering the radio-optical source separation and source densities of the parent catalogues, we suggest that it is statistically probable that the M dwarfs are associated with the radio emission. However, it remains plausible that for some of the sources the radio emission originates from an optically faint and red galaxy hiding behind the M dwarf. The isotropic radio luminosities (∼10
17 − 18
erg s
−1
Hz
−1
) of the M dwarfs suggest that if the association is real, the radio emission is likely driven by a coherent emission process produced via plasma or electron-cyclotron maser instability processes, which is potentially caused by binary interaction. Long-term monitoring in the radio and high-resolution radio follow-up observations are necessary to search for any variability and pinpoint the radio emission to determine whether our tentative conclusion that these ultracool dwarfs are radio emitting is correct. If the low-frequency radio emission is conclusively associated with the M dwarfs, this would reveal a new population of optically faint and distant (> 50 pc) radio-emitting M dwarfs.
Abstract
The MeerKAT Absorption Line Survey (MALS) has observed 391 telescope pointings at the
L
band (900–1670 MHz) at
δ
≲ +20°. We present radio continuum images and a catalog of 495,325 (240,321) ...radio sources detected at a signal-to-noise ratio (S/N) > 5 over an area of 2289 deg
2
(1132 deg
2
) at 1006 MHz (1381 MHz). Every MALS pointing contains a central bright radio source (
S
1 GHz
≳ 0.2 Jy). The median spatial resolution is 12″ (8″). The median rms noise away from the pointing center is 25
μ
Jy beam
−1
(22
μ
Jy beam
−1
) and is within ∼15% of the achievable theoretical sensitivity. The flux density scale ratio and astrometric accuracy deduced from multiply observed sources in MALS are <1% (8% scatter) and 1″, respectively. Through comparisons with NVSS and FIRST at 1.4 GHz, we establish the catalog’s accuracy in the flux density scale and astrometry to be better than 6% (15% scatter) and 0.″8, respectively. The median flux density offset is higher (9%) for an alternate beam model based on holographic measurements. The MALS radio source counts at 1.4 GHz are in agreement with literature. We estimate spectral indices (
α
) of a subset of 125,621 sources (S/N > 8), confirm the flattening of spectral indices with decreasing flux density, and identify 140 ultra-steep-spectrum (
α
< −1.3) sources as prospective high-
z
radio galaxies (
z
> 2). We have identified 1308 variable and 122 transient radio sources comprising primarily active galactic nuclei that demonstrate long-term (26 yr) variability in their observed flux densities. The MALS catalogs and images are publicly available at
https://mals.iucaa.in
.
The cosmic radio dipole is an anisotropy in the number counts of radio sources, analogous to the dipole seen in the cosmic microwave background (CMB). Measurements of source counts of large radio ...surveys have shown that though the radio dipole is generally consistent in direction with the CMB dipole, the amplitudes are in tension. These observations present an intriguing puzzle as to the cause of this discrepancy, with a true anisotropy breaking with the assumptions of the cosmological principle, invalidating the most common cosmological models that are built on these assumptions. We present a novel set of Bayesian estimators to determine the cosmic radio dipole and compare the results with commonly used methods on the Rapid ASKAP Continuum Survey (RACS) and the NRAO VLA Sky Survey (NVSS) radio surveys. In addition, we adapt the Bayesian estimators to take into account systematic effects known to affect such large radio surveys, folding information such as the local noise floor or array configuration directly into the parameter estimation. The enhancement of these estimators allows us to greatly increase the amount of sources used in the parameter estimation, yielding tighter constraints on the cosmic radio dipole estimation than previously achieved with NVSS and RACS. We extend the estimators further to work on multiple catalogues simultaneously, leading to a combined parameter estimation using both NVSS and RACS. The result is a dipole estimate that perfectly aligns with the CMB dipole in terms of direction but with an amplitude that is three times as large, and a significance of 4.8\(\sigma\). This new dipole measurement is made to an unprecedented level of precision for radio sources, which is only matched by recent results using infrared quasars.