ABSTRACT

Kardashev Type III civilisations have by definition energy requirements that are likely to generate strong excess emission in the mid-infrared (MIR) that is associated with the waste heat they generate. For those civilisations that capture starlight via Dyson sphere like constructions, the Optical to MIR flux ratio of the host galaxies is also expected to be unusually low. Since a wide range of galaxy types adhere to the infrared-radio correlation (IRC), galaxies hosting Type III civilisations should also strongly deviate from this relation. Radio data can therefore play a crucial role in breaking the degeneracy between the effects of dust obscuration in a galaxy and the signature of an energy-intensive civilisation. We have used the newly released LoTSS-DR1 value-added catalogue to create a sample of 16 367 z < 0.2 sources with radio and MIR flux densities, optical photometry and (photometric) redshifts. We calculated the IRC parameter |$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| = log(S150 MHz/|$S_{22~{\rm \mu m}}$|⁠) and identified 21 sources with exceptionally high values of |$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$|⁠, an indication of MIR emission enhanced by a factor of ∼10. Out of the 21 high q-value sources, 4 sources have very red MIR colours, and appear to have relatively low optical/MIR ratios. Two of the four sources are not well known in the literature, they are considered as potential hosts of Type III civilisations. These sources deserve further study and investigation. Finally, we note that extending our analysis to the full LoTSS survey area can place very strong constraints on the incidence of Type III civilisations in the Universe.

1 INTRODUCTION

The well-established far-infrared-radio correlation applies to a wide variety of galaxies in the local Universe (see Solarz et al. 2019 and references therein). Similar trends are also seen in the mid-infrared (MIR)-radio domain (Donley et al. 2005; Park et al. 2008). As demonstrated by Garrett (2002) and Elbaz et al. (2002), the correlation continues to apply at cosmological distances and over 7 orders of magnitude in luminosity. The tightness of the Infrared Radio Correlation (hereafter IRC) is usually interpreted as the non-thermal radio emission and the thermal IR emission, both being related to processes associated with massive star formation (see Magnelli et al. 2015 for a review). However, the IRC also applies to non-star-forming galaxies, including many radio-quiet AGN. Radio-loud AGN (RL AGN; L1.4 GHz ≳ 1024 W Hz−1 in the local Universe) deviate from the IRC, presenting very low values of q, the latter being defined as the logarithmic ratio between the IR and radio flux densities. The value of q has often been used to disentangle AGN activity from star-forming processes in extragalactic objects, e.g. Condon, Cotton & Broderick (2002), Sargent et al. (2012), and Delhaize et al. (2017). On the other hand, radio-quiet AGN (RQ AGN), or very young star-forming galaxies (SFGs) with comparably low radio emission may cause elevated values of q (Roussel et al. 2003; Costagliola et al. 2016). So far, a satisfactory physical model that is consistent with all observed constraints and that fully explains the IRC and its outliers is still to be fully developed (Read et al. 2018; Wang et al. 2019). A better understanding of the astrophysical nature of sources with very high values of q is useful in progressing studies in this field.

In addition to natural astrophysical objects like RQ AGNs and young SFGs possessing high values of q, it has been proposed that artificial sources, such as Kardashev Type III civilisations (Kardashev 1964) can also present the same unusual signature Garrett (2015). A Kardashev Type III civilisation is defined as being able to draw on a large fraction of the stellar energy available from its host galaxy (∼4 × 1037 W; Kardashev 1964). As suggested by Wright et al. (2014a,b), Type III civilisations may be identified via their waste heat signature which is expected to peak at MIR wavelengths, assuming peak temperatures of ∼100−600 K (Carrigan 2009). A significant fraction of the optical emission in the galaxy is also expected to be absorbed by artificial structures such as Dyson Spheres or Swarms (Dyson 1960). This leads to the identification of three characteristics that a galaxy hosting a Type III civilisation might be expected to have: (1) extremely red MIR colours, (2) high values of q, and (3) unusually low optical/IR luminosity ratios.

To characterize the first effect in the MIR, Wright et al. (2014a) developed a formalism to define the fractional starlight re-radiated as alien waste heat (γ) by fitting the measured WISE (Wide-field Infrared Survey Explorer) MIR data to a spectral energy distribution (SED) model for normal spiral galaxies. This formalism reveals the mechanism of how the waste heat from a Type III civilisation will potentially generate excessive MIR emission with very red colours.

Since the radio emission of a galaxy is largely expected to remain unaffected by the presence of a Type III civilisation, such systems should present large values of q, well in excess of the average value (Garrett 2015). As suggested by Wright et al. (2014a), with an approximation to a blackbody with temperature T* and the Rayleigh–Jeans limit (λ ≫ hc/kT*), the expected SED of a home galaxy of Type III civilisation can be written as
(1)
where Twaste is the characteristic temperature of the waste heat; and Fλ(γ = 0) is the flux density of a normal galaxy.

For example, taking T* ∼ 6000 K, and Twaste ∼ 100–600 K, and assuming values of γ ∼ 0.01−0.1, the MIR flux will be enhanced by a factor of 10−100. If the typical value of q is ∼0.5, then values in excess of ∼1.5–2.5 might be anticipated.

These estimates are likely to hold for quiescent galaxies with little to no star formation (Ilbert et al. 2013). Quiescent galaxies account for about half of the stellar mass in the Universe since z ∼ 1, so they may be important as hosts of advanced civilisations. However, the Rayleigh-Jeans limit would fail in typical star-forming galaxies, in which the MIR radiation is mostly contributed by cold dust. We quantify the MIR enhancement for such sources in Section 3.1.

Griffith et al. (2015) constructed a sample of 93 sources (the Ĝ sample) with γ > 0.25 out of ∼100 000 WISE sources with good photometry. Many of these sources has little or no presence in the scientific literature and were proposed as potential advanced Type III civilisations. However, Garrett (2015) cross-matched the Ĝ sample with the NRAO/VLA Sky Survey (NVSS) catalogue and found that the vast majority of the sources have q values consistent with the IRC. This suggests that the MIR excess and prominent red colours of the sources was likely to be associated with natural astrophysical processes. It turns out that samples solely based on MIR data produce many false-positives – star-forming galaxies naturally present a MIR excess as a consequence of the reprocessing of optical/UV emission by dust. Indeed, searches for anomalies in astronomical data that target deviations from galaxy scaling laws in general (e.g. outliers in the Tully–Fisher relation; Annis 1999) have similar shortcomings. Radio data can therefore play a crucial role in breaking this degeneracy between the effects of dust obscuration in a galaxy and the signature of an energy-intensive civilisation.

The newly released LOw Frequency ARray (LOFAR; van Haarlem et al. 2013) Two-metre Sky Survey first data release (LoTSS-DR1; Shimwell et al. 2017) value-added catalogue (Duncan et al. 2019; Shimwell et al. 2019; Williams et al. 2019) contains ∼300 000 radio sources cross-matched with the WISE catalogue. This is an excellent resource for the investigation of the nature of sources with high values of q and the identification of potential Kardashev Type III civilisations. Optical measurements are also available for these sources, in particular Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) i-band photometry data. In this work, we perform an IRC and colour distribution analysis based on the LoTSS-DR1 value-added catalogue to build a sample of sources which display high values of q, extreme red MIR colours, and unusually low optical/MIR flux ratios. It is notable however, our radio-selected sample might be biased to starbursts and AGN with strong radio emission. To calculate q-values, the measured LoTSS radio flux is required, therefore host galaxies of Type III civilisations with radio flux below the detection threshold (0.6 mJy) of LoTSS-DR1 would be missed in this methodology. On the other hand, very hot Dyson spheres like the “Matrioshka brain” (Bradbury, Cirkovic & Dvorsky 2011) could reach a waste heat of ∼1000 K, however, since Type III civilisations with Twaste ≳ 1000 K would peak at a wavelength below 3.4 |$\mu$|m and therefore be excluded in the colour distribution analysis, we only look for Type III civilisations with Twaste ∼ 100–600 K.

Our approach is complimentary to the earlier work of Wright et al. (2014a,b) and Griffith et al. (2015) but takes advantage of the new large-scale radio surveys now being conducted by telescopes such as LOFAR and the other SKA pathfinders/precursors. We introduce the LoTSS-DR1 value-added catalogue in Section 2, followed by a description of our high-q source selection in Section 3. In Section 4, we present MIR colour and optical magnitude analysis to find potential host galaxies of Type III civilisations. A summary and suggestion for future steps are given in Section 5.

2 THE LOFAR DR1 VALUE-ADDED CATALOGUE

LoTSS is an ambitious low-frequency radio survey with a wide range of science goals that aims to observe the entire northern sky at frequencies of 120–168 MHz with a resolution of ∼6 arcsec. Using the High Band Antenna (HBA) system of LOFAR, the survey reaches a sensitivity of less than 0.1 mJy beam−1. The first LoTSS data release covers 424 deg2 over a field centred on the Hobby–Eberly Telescope Dark Energy Experiment (HETDEX; Hill et al. 2008). In a series of papers, Shimwell et al. (2019), Williams et al. (2019), and Duncan et al. (2019) built up a LoTSS-DR1 “value-added” catalogue1 including 318 520 discrete radio sources identified by the Python Blob Detector and Source Finder (PyBDSF; Mohan & Rafferty 2015). The value-added catalogue was cross-matched with the Pan-STARRS 3π Steradian survey (Chambers et al. 2016), and the All-sky WISE (All-WISE) catalogue (Cutri et al. 2021). The source associations and optical and/or IR identifications were done by using a combination of statistical techniques and visual inspections (Williams et al. 2019). Photometric redshift information is given for 94.4 per cent of the LoTSS-DR1 sources (Duncan et al. 2019).

The All-WISE catalogue includes photometric data from four MIR channels at 3.4, 4.6, 12, and 22 |$\mu$|m (W1, W2, W3, and W4) for more than 747 million sources over the full sky. The WISE images have a typical point spread function of 6−6.5 arcsec for W1, W2 and 12 arcsec for W3 and W4.

The Pan-STARRS 3π Steradian Survey was observed in five bands (grizy) covering the entire northern sky of ð > 30°. The mean 5σ point source limiting sensitivities for the five bands are 23.3, 23.2, 23.1, 22.3, and 21.4 mag, respectively. There is a systematic uncertainty of 7–12 millimag in the photometric calibration depending on the bandpass.

There are 231 716 (73 per cent) sources in the final LoTSS-DR1 value-added catalogue having a WISE and/or Pan-STARRS counterpart, amongst them, 55 694 were detected in WISE W4 and have profile fitted W4 flux density data, which can be used to measure the IRC parameter |$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| (see Section 3). While the LoTSS value-added catalogue is dominated by extragalactic sources, we cannot rule out the possibility that our sample may be contaminated by a few Galactic objects. However, considering that the HETDEX field is ∼60° off of the Galactic plane, we do not expect this to be a major problem for our analysis.

In this paper, the prefix “ILT” before the sourcenames refers to the LoTSS catalogue.

3 SAMPLE OF HIGH Q-VALUE SOURCES

3.1 Candidate selection

In Section 1, we use the SED model of a Type III civilisation with Rayleigh–Jeans limit in Wright et al. (2014b) to estimate the enhancement of their MIR emission. It should be noted that these estimates are likely to hold for quiescent galaxies but may not be applicable to all galaxy types. Many quiescent galaxies are early-type galaxies (ETGs). As shown by Nyland et al. (2017), ∼19−50 per cent of ETGs have higher IR-radio ratios than what is typical for SFGs. Possible explanations for this include the presence of nascent SF in these systems or specific properties of the galaxy e.g. a bottom-heavy stellar IMF, enhanced cosmic ray escape, weak magnetic fields, or other environmental effects. Looking for sources with q-value +1 than normal would filter out most of the high-q ETG outliers, but there may still be false-positives. For star-forming galaxies, however, the majority of the 22 |$\mu$|m emission is not from stars but from cold dust (∼10−20 K) and small grains, the Rayleigh–Jeans limit would not be applicable here. According to fig. 15 of Draine & Li (2007), the value of νLν (ν =  24 ≈ 22 |$\mu$|m) in such galaxies is about 3|$-6{{\ \rm per\ cent}}$| of the total flux. As shown in Fig. 1, the value of νLν (ν = 22 |$\mu$|m)/Ltotal peaks at ∼0.7 at 160 K. If an advanced civilisation captured all starlight with Dyson spheres before it can be absorbed by dust, νLν (ν = 22 |$\mu$|m)/Ltotal raises to ∼10−70 per cent for Dyson spheres with a blackbody temperature of 100−600 K. In other words, these galaxies would see their MIR luminosity enhanced by a factor of 1.7−23.3, resulting in |$q_{\rm 22~\mu m}$|-values increased by up to ∼+1.4 dex in excess of normal values. For star-forming systems with an initial νLν (ν = 22 |$\mu$|m)/Ltotal = 0.03 and a Dyson sphere with Twaste ≲ 100 or ≳ 400 K, the enhanced νLν (ν = 22 |$\mu$|m)/Ltotal value would be ≲ 0.3, which is less than 10 times to the original value, these sources may be largely contaminated by AGN.

Fractional 22 $\mu$m radiation to the total flux of a blackbody (Dyson sphere), the horizontal lines denote where Lν (ν  =  22 $\mu$m)/Ltotal = 0.3 and 0.6, the lowest values required for star-forming galaxies with νLν (ν = 22 $\mu$m)/Ltotal = 0.03 and 0.06 to be enlarged by at least 10 times.
Figure 1.

Fractional 22 |$\mu$|m radiation to the total flux of a blackbody (Dyson sphere), the horizontal lines denote where Lν (ν  =  22 |$\mu$|m)/Ltotal = 0.3 and 0.6, the lowest values required for star-forming galaxies with νLν (ν = 22 |$\mu$|m)/Ltotal = 0.03 and 0.06 to be enlarged by at least 10 times.

In this paper, we only look for Type III civilisations with their MIR radiation at least 10 times enlarged, this would allow us to find Type III civilisations in quiescent galaxies with γ ∼ 0.01 and Twaste ∼ 100–600 K or star-forming galaxies with Twaste ∼ 100–400 K. Galaxies in which the MIR emission is being enhanced by a factor of 10, should present q values raised linearly by +1.

We use the WISE 22 |$\mu$|m and the 120−160 MHz LOFAR measurements to calculate the IRC parameter: |$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$|⁠, which is defined as log(⁠|$S_{22~{\rm \mu m}}$|/S150 MHz). We restricted ourselves to galaxies with redshifts <0.2 in order to have good auxiliary data for sources that warranted more detailed investigation. Another consequence of this limited redshift range is that a k-correction need not be applied to the data. To check this assumption, we divided the sources into two sub-samples: z < 0.1 and 0.1 < z < 0.2. Out of the ∼300 000 LoTSS-DR1 sources, there are 5790 and 10 577 sources with WISE W4 photometric flux data, and redshifts of z < 0.1 and 0.1 < z < 0.2, respectively. The z < 0.1 sources have an average |$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| value and a standard deviation of ∼0.49 ± 0.37, while the sources with 0.1 < z < 0.2 have a mean |$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| of ∼0.48 ± 0.30. Fig. 2 shows the distribution of |$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| value versus calculated WISE 22 |$\mu$|m luminosity. As shown in Fig. 2, |$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| tends to increase along the 22 |$\mu$|m luminosity, this trend probably arises from the higher number of starburst galaxies in the more luminous sample, for which a more intense radiation field and hot dust are expected. In addition, starburst galaxies tend to have slightly flatter radio spectra and more high-frequency emission at >5 GHz due to free–free from hot gas. As illustrated by left-hand panel of Fig. 2, in the low 22 |$\mu$|m luminosity and low redshift regime, the fit seems much stronger than any actual trend. To filter out false-positive in low-z and low-|$L_{\rm 22~\mu m}$| regime, we plot the median q22 value of sources in Yun, Reddy & Condon (2001) calculated by Nyland et al. (2017), and a ± log105 constraints for normal galaxies, corresponding to five times the IR-excess, which applies to most of galaxies in Yun et al. (2001). The median q22 value of Yun et al. (2001) sample of 0.99 was defined as the logarithmic ratio between |$L_{\rm 22~\mu m}$| and L1.4 GHz (Nyland et al. 2017), we convert their median q22-value to |$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| assuming a spectral index between L1.4 GHz and L150 MHz of −0.755, which is the averaged value for AGN and SFGs measured in Calistro Rivera et al. (2017). As a result, in the scenario of Yun et al. (2001), most galaxies would have a |$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| value within ∼−0.45 to 0.95, and a median value of ∼0.25. Fig. 3 is a histogram of distance of |$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| value to the linear fit over |$L_{\rm 22~\mu m}$|⁠. For clarity, the source count scale is logarithmic. As demonstrated by Figs 2 and 3, the majority (∼99.9 per cent) of the sources fell below the linear fits +1 and 0.25 + log105 (= 0.95).

$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$ vs WISE W4 luminosity for sources with (a) z < 0.1 and (b) 0.1 < z < 0.2. The solid and dashed lines denote the least square linear fitting over the sources, the dashed line restrain the region between ± 1 to the fitting line, the solid and dashed red lines define the median $q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$ value of ∼0.25 for sources in Yun et al. (2001), and the five-time IR-excess/deficit.
Figure 2.

|$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| vs WISE W4 luminosity for sources with (a) z < 0.1 and (b) 0.1 < z < 0.2. The solid and dashed lines denote the least square linear fitting over the sources, the dashed line restrain the region between ± 1 to the fitting line, the solid and dashed red lines define the median |$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| value of ∼0.25 for sources in Yun et al. (2001), and the five-time IR-excess/deficit.

The histogram of distance of $q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$ value to the linear fit over $L_{\rm 22~\mu m}$ in a logarithmic scale, the red bins represent sources with $q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$ values larger than the linear fitting +1 and higher than 0.25 + log105.
Figure 3.

The histogram of distance of |$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| value to the linear fit over |$L_{\rm 22~\mu m}$| in a logarithmic scale, the red bins represent sources with |$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| values larger than the linear fitting +1 and higher than 0.25 + log105.

We classify sources with |$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| larger than the linear fits +1 and higher than 0.25 + log105 as high q-value sources. There are 20 sources selected in the z < 0.1 sample and 1 source in the 0.1 < z < 0.2 sample. We therefore identify these 21 (20 + 1) sources as parent sample of possible Kardashev Type III candidates.

We investigate the physical properties of these high-|$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| value sources in the following sections.

3.2 Source classification

We used SIMBAD2 to search within 10 arcsec around the positions of the 21 high-|$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| value sources looking for object identifiers in the literature. An overview of the source classification is shown in Table 1. The majority (19/21) of the sources are previously identified in the literature. Of the 19 sources with identifiers, nine of them are identified as AGN, four are classified as emission-line galaxies, as well as one composite galaxy showing both star-forming and AGN features, the other six sources are confirmed to be galaxies but with no specific type designated. We have calculated the radio luminosity of the sources – they are all radio-faint with L150 MHz ≲ 1023 W Hz−1. The source properties, including radio and MIR flux densities, radio luminosity, |$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| values, and redshift information, as well as the names of the identifiers in the literature for the high-|$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| value sources are listed in Table 2.

Table 1.

Numbers of each type of sources with high value of |$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$|⁠.

Typea (total counts)AGNb (9)EmGc (4)Comp (1)Other (5)No identifier (2)Total sources
ClassificationdSy1Sy2LINAGNQSOH2GEmGCompGGiP
Counts4221132141221
Typea (total counts)AGNb (9)EmGc (4)Comp (1)Other (5)No identifier (2)Total sources
ClassificationdSy1Sy2LINAGNQSOH2GEmGCompGGiP
Counts4221132141221
Notes. a

Counts include all selected sources with redshift up to 0.2.

b

This sub-sample contains a multi-associated sources: ILT J112637.84+513422.9 (SIMBAD identifier: MCG+09-19-078) that it classified both as a Seyfert 1 Galaxy in Dong et al. (2012) (source ID: J112637.74 + 513423.0) and Seyfert 2 Galaxy in Toba et al. (2014) (source ID: MCG +09-19-078).

c

This sub-sample contains a multi-associated source, ILT J140956.92+545649.5 (SIMBAD identifier: LEDA 2816082), that is classified as an emission-line galaxy in (Shim & Chary 2013; source ID: SDSS 587735696987389965) and an H ii galaxy in Shirazi & Brinchmann (2012) (source ID: 1325-52762-412.)

d

Sy1: Seyfert 1 Galaxy; Sy2: Seyfert 2 Galaxy; LIN: LINER-type AGN; QSO: quasar object; EmG: Emission-line galaxy; H2G: H ii galaxy; Comp: Composite galaxy with both star-forming and AGN activity; G: Galaxy; GiP: Galaxy in pair of galaxies.

Table 1.

Numbers of each type of sources with high value of |$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$|⁠.

Typea (total counts)AGNb (9)EmGc (4)Comp (1)Other (5)No identifier (2)Total sources
ClassificationdSy1Sy2LINAGNQSOH2GEmGCompGGiP
Counts4221132141221
Typea (total counts)AGNb (9)EmGc (4)Comp (1)Other (5)No identifier (2)Total sources
ClassificationdSy1Sy2LINAGNQSOH2GEmGCompGGiP
Counts4221132141221
Notes. a

Counts include all selected sources with redshift up to 0.2.

b

This sub-sample contains a multi-associated sources: ILT J112637.84+513422.9 (SIMBAD identifier: MCG+09-19-078) that it classified both as a Seyfert 1 Galaxy in Dong et al. (2012) (source ID: J112637.74 + 513423.0) and Seyfert 2 Galaxy in Toba et al. (2014) (source ID: MCG +09-19-078).

c

This sub-sample contains a multi-associated source, ILT J140956.92+545649.5 (SIMBAD identifier: LEDA 2816082), that is classified as an emission-line galaxy in (Shim & Chary 2013; source ID: SDSS 587735696987389965) and an H ii galaxy in Shirazi & Brinchmann (2012) (source ID: 1325-52762-412.)

d

Sy1: Seyfert 1 Galaxy; Sy2: Seyfert 2 Galaxy; LIN: LINER-type AGN; QSO: quasar object; EmG: Emission-line galaxy; H2G: H ii galaxy; Comp: Composite galaxy with both star-forming and AGN activity; G: Galaxy; GiP: Galaxy in pair of galaxies.

Table 2.

Information for the sources with high values of |$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$|⁠.

LOFAR ID (J2000)S150 MHzL150 MHz|$S_{22~{\rm \mu m}}$||$L_{22~{\rm \mu m}}$|zmagW1magW2magW3magW4magi|$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$|IdentifierSource TypeReference
[RA, DEC][mJy][W Hz−1][μ Jy][W Hz−1]
(1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)(12)(13)(14)(15)
J104725.50+535338.50.5165.399 × 102125298.9092.648 × 10230.06417.27717.31915.07212.89217.5631.691LEDA 139444LIN AGN1
J112637.84+513422.90.6791.156 × 102123979.2504.080 × 10220.02615.22815.46314.45312.95016.0601.548MCG+09-19-078Sy1/21, 2
J112804.87+555630.81.6451.412 × 102160454.8325.188 × 10220.01915.84416.27413.92011.94616.0021.5652MASX J11280486+5556303G1
J115039.20+552112.70.6401.509 × 101913153.9043.104 × 10200.00315.10915.76414.44513.60215.7141.313NGC 3913G3
J120522.51+502110.86.4639.774 × 1019351854.5925.321 × 10210.00312.97913.37411.18110.03415.0061.736NGC 4085GiP4
J120658.79+524126.51.3182.532 × 1022112931.8142.170 × 10240.08517.38616.23312.88711.26818.2211.933SDSS J120658.79+524126.4Sy11
J122615.87+482937.71.4833.472 × 101853937.7281.263 × 10200.00116.99616.97914.45612.07017.4191.561SDSS J122615.72+482938.6H2G5
J123803.78+461819.90.4511.187 × 102247285.6351.246 × 10240.09918.20216.85414.35012.21318.3312.021SDSS J123803.75+461820.1EmG6
J124505.63+561430.20.6447.668 × 10166134.2637.300 × 10170.000217.98217.12815.27314.43118.4940.979J124505.6+561430QSO7
J130708.14+461834.22.2033.254 × 1021105655.4861.560 × 10230.02515.37615.74613.27711.34015.8841.6812MASX J13070814+4618346AGN8
J131048.89+544111.50.3353.493 × 102116538.9281.722 × 10230.06415.93715.82514.44113.35416.6021.6932MASX J13104883+5441104Sy21
J131459.85+511621.40.2684.042 × 10207364.1171.110 × 10220.02515.77416.37215.48214.23215.9641.4392MASX J13145978+5116213LIN AGN1
J131506.48+461409.70.4531.180 × 102225132.7246.547 × 10230.09815.78915.44413.99712.89917.1421.7442MASS J13150651+4614101Sy19
J131856.77+474609.61.2595.009 × 102155490.3032.208 × 10230.04015.86316.17713.82412.03916.3761.6442MASX J13185676+4746092G1
J134649.72+542621.70.4693.741 × 102014262.2381.138 × 10220.01815.88915.45014.53613.51518.5041.483
J140411.25+542519.61.1662.358 × 101821906.9504.432 × 10190.00117.37717.47115.07313.04918.2891.274SDSS J140411.24+542518.6H2G5
J140956.92+545649.50.7291.142 × 102237215.4665.832 × 10230.07717.63417.09914.38012.47317.4671.708LEDA 2816082EmG/H2G5, 6
J141424.94+465348.40.3462.227 × 102227003.1921.740 × 10240.15015.67915.37114.23012.82117.1471.8932MASSJ14142489+4653485Sy110
J143551.22+501118.20.4987.584 × 102119785.4773.012 × 10230.07617.13317.49515.30713.15917.2691.599SDSS J143551.27+501118.3G1
J144632.38+490038.20.5811.987 × 102192241.2233.153 × 10230.03716.53616.62713.54611.48816.8202.201SDSS J144632.39+490038.0Comp11
J145757.90+565323.80.8064.262 × 102151843.2542.740 × 10230.04617.27617.19014.29612.11317.6181.808
LOFAR ID (J2000)S150 MHzL150 MHz|$S_{22~{\rm \mu m}}$||$L_{22~{\rm \mu m}}$|zmagW1magW2magW3magW4magi|$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$|IdentifierSource TypeReference
[RA, DEC][mJy][W Hz−1][μ Jy][W Hz−1]
(1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)(12)(13)(14)(15)
J104725.50+535338.50.5165.399 × 102125298.9092.648 × 10230.06417.27717.31915.07212.89217.5631.691LEDA 139444LIN AGN1
J112637.84+513422.90.6791.156 × 102123979.2504.080 × 10220.02615.22815.46314.45312.95016.0601.548MCG+09-19-078Sy1/21, 2
J112804.87+555630.81.6451.412 × 102160454.8325.188 × 10220.01915.84416.27413.92011.94616.0021.5652MASX J11280486+5556303G1
J115039.20+552112.70.6401.509 × 101913153.9043.104 × 10200.00315.10915.76414.44513.60215.7141.313NGC 3913G3
J120522.51+502110.86.4639.774 × 1019351854.5925.321 × 10210.00312.97913.37411.18110.03415.0061.736NGC 4085GiP4
J120658.79+524126.51.3182.532 × 1022112931.8142.170 × 10240.08517.38616.23312.88711.26818.2211.933SDSS J120658.79+524126.4Sy11
J122615.87+482937.71.4833.472 × 101853937.7281.263 × 10200.00116.99616.97914.45612.07017.4191.561SDSS J122615.72+482938.6H2G5
J123803.78+461819.90.4511.187 × 102247285.6351.246 × 10240.09918.20216.85414.35012.21318.3312.021SDSS J123803.75+461820.1EmG6
J124505.63+561430.20.6447.668 × 10166134.2637.300 × 10170.000217.98217.12815.27314.43118.4940.979J124505.6+561430QSO7
J130708.14+461834.22.2033.254 × 1021105655.4861.560 × 10230.02515.37615.74613.27711.34015.8841.6812MASX J13070814+4618346AGN8
J131048.89+544111.50.3353.493 × 102116538.9281.722 × 10230.06415.93715.82514.44113.35416.6021.6932MASX J13104883+5441104Sy21
J131459.85+511621.40.2684.042 × 10207364.1171.110 × 10220.02515.77416.37215.48214.23215.9641.4392MASX J13145978+5116213LIN AGN1
J131506.48+461409.70.4531.180 × 102225132.7246.547 × 10230.09815.78915.44413.99712.89917.1421.7442MASS J13150651+4614101Sy19
J131856.77+474609.61.2595.009 × 102155490.3032.208 × 10230.04015.86316.17713.82412.03916.3761.6442MASX J13185676+4746092G1
J134649.72+542621.70.4693.741 × 102014262.2381.138 × 10220.01815.88915.45014.53613.51518.5041.483
J140411.25+542519.61.1662.358 × 101821906.9504.432 × 10190.00117.37717.47115.07313.04918.2891.274SDSS J140411.24+542518.6H2G5
J140956.92+545649.50.7291.142 × 102237215.4665.832 × 10230.07717.63417.09914.38012.47317.4671.708LEDA 2816082EmG/H2G5, 6
J141424.94+465348.40.3462.227 × 102227003.1921.740 × 10240.15015.67915.37114.23012.82117.1471.8932MASSJ14142489+4653485Sy110
J143551.22+501118.20.4987.584 × 102119785.4773.012 × 10230.07617.13317.49515.30713.15917.2691.599SDSS J143551.27+501118.3G1
J144632.38+490038.20.5811.987 × 102192241.2233.153 × 10230.03716.53616.62713.54611.48816.8202.201SDSS J144632.39+490038.0Comp11
J145757.90+565323.80.8064.262 × 102151843.2542.740 × 10230.04617.27617.19014.29612.11317.6181.808

Note. (1) Source ID in the LoTSS value-added catalogue, (2) the total, integrated Stokes I 150 MHz radio flux density measured by the LOFAR HBA, (3) radio luminosity derived from LoTSS 150 MHz radio flux density, (4) WISE W4 profile-fit photometry flux, (5) calculated WISE 22μ luminosity, (6) Best available redshift estimate, (7)−(10) WISE W1−W4 profile-fit photometry magnitude, (11) Pan-STARSS forced aperture AB magnitude, (12) IRC parameter between the WISE 22 |$\mu$|m and LoTSS 150 MHz flux density, (13) Identifier archived in the SIMBAD catalogue, (14) Source type, G: Galaxy; Sy2: Seyfert 2 Galaxy; Sy1: Seyfert 1 Galaxy; Comp: Composite galaxy that likely have both star-forming and AGN activity; GiC: Galaxy in cluster of galaxies; LIN AGN: LINER-type AGN; QSO: Quasar object; EmG: Emission-line galaxy; X: X-ray source. (15) References; 1: Toba et al. (2014); 2: Dong et al. (2012); 3: Lavaux & Hudson (2011); 4: Crook et al. (2007); 5: Shirazi & Brinchmann (2012); 6: Shim & Chary (2013); 7: Allen et al. (2011); 8: Liu et al. (2011); 9: Oh et al. (2015); 10: Zhou et al. (2006); 11: Brinchmann, Kunth & Durret (2008).

Table 2.

Information for the sources with high values of |$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$|⁠.

LOFAR ID (J2000)S150 MHzL150 MHz|$S_{22~{\rm \mu m}}$||$L_{22~{\rm \mu m}}$|zmagW1magW2magW3magW4magi|$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$|IdentifierSource TypeReference
[RA, DEC][mJy][W Hz−1][μ Jy][W Hz−1]
(1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)(12)(13)(14)(15)
J104725.50+535338.50.5165.399 × 102125298.9092.648 × 10230.06417.27717.31915.07212.89217.5631.691LEDA 139444LIN AGN1
J112637.84+513422.90.6791.156 × 102123979.2504.080 × 10220.02615.22815.46314.45312.95016.0601.548MCG+09-19-078Sy1/21, 2
J112804.87+555630.81.6451.412 × 102160454.8325.188 × 10220.01915.84416.27413.92011.94616.0021.5652MASX J11280486+5556303G1
J115039.20+552112.70.6401.509 × 101913153.9043.104 × 10200.00315.10915.76414.44513.60215.7141.313NGC 3913G3
J120522.51+502110.86.4639.774 × 1019351854.5925.321 × 10210.00312.97913.37411.18110.03415.0061.736NGC 4085GiP4
J120658.79+524126.51.3182.532 × 1022112931.8142.170 × 10240.08517.38616.23312.88711.26818.2211.933SDSS J120658.79+524126.4Sy11
J122615.87+482937.71.4833.472 × 101853937.7281.263 × 10200.00116.99616.97914.45612.07017.4191.561SDSS J122615.72+482938.6H2G5
J123803.78+461819.90.4511.187 × 102247285.6351.246 × 10240.09918.20216.85414.35012.21318.3312.021SDSS J123803.75+461820.1EmG6
J124505.63+561430.20.6447.668 × 10166134.2637.300 × 10170.000217.98217.12815.27314.43118.4940.979J124505.6+561430QSO7
J130708.14+461834.22.2033.254 × 1021105655.4861.560 × 10230.02515.37615.74613.27711.34015.8841.6812MASX J13070814+4618346AGN8
J131048.89+544111.50.3353.493 × 102116538.9281.722 × 10230.06415.93715.82514.44113.35416.6021.6932MASX J13104883+5441104Sy21
J131459.85+511621.40.2684.042 × 10207364.1171.110 × 10220.02515.77416.37215.48214.23215.9641.4392MASX J13145978+5116213LIN AGN1
J131506.48+461409.70.4531.180 × 102225132.7246.547 × 10230.09815.78915.44413.99712.89917.1421.7442MASS J13150651+4614101Sy19
J131856.77+474609.61.2595.009 × 102155490.3032.208 × 10230.04015.86316.17713.82412.03916.3761.6442MASX J13185676+4746092G1
J134649.72+542621.70.4693.741 × 102014262.2381.138 × 10220.01815.88915.45014.53613.51518.5041.483
J140411.25+542519.61.1662.358 × 101821906.9504.432 × 10190.00117.37717.47115.07313.04918.2891.274SDSS J140411.24+542518.6H2G5
J140956.92+545649.50.7291.142 × 102237215.4665.832 × 10230.07717.63417.09914.38012.47317.4671.708LEDA 2816082EmG/H2G5, 6
J141424.94+465348.40.3462.227 × 102227003.1921.740 × 10240.15015.67915.37114.23012.82117.1471.8932MASSJ14142489+4653485Sy110
J143551.22+501118.20.4987.584 × 102119785.4773.012 × 10230.07617.13317.49515.30713.15917.2691.599SDSS J143551.27+501118.3G1
J144632.38+490038.20.5811.987 × 102192241.2233.153 × 10230.03716.53616.62713.54611.48816.8202.201SDSS J144632.39+490038.0Comp11
J145757.90+565323.80.8064.262 × 102151843.2542.740 × 10230.04617.27617.19014.29612.11317.6181.808
LOFAR ID (J2000)S150 MHzL150 MHz|$S_{22~{\rm \mu m}}$||$L_{22~{\rm \mu m}}$|zmagW1magW2magW3magW4magi|$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$|IdentifierSource TypeReference
[RA, DEC][mJy][W Hz−1][μ Jy][W Hz−1]
(1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)(12)(13)(14)(15)
J104725.50+535338.50.5165.399 × 102125298.9092.648 × 10230.06417.27717.31915.07212.89217.5631.691LEDA 139444LIN AGN1
J112637.84+513422.90.6791.156 × 102123979.2504.080 × 10220.02615.22815.46314.45312.95016.0601.548MCG+09-19-078Sy1/21, 2
J112804.87+555630.81.6451.412 × 102160454.8325.188 × 10220.01915.84416.27413.92011.94616.0021.5652MASX J11280486+5556303G1
J115039.20+552112.70.6401.509 × 101913153.9043.104 × 10200.00315.10915.76414.44513.60215.7141.313NGC 3913G3
J120522.51+502110.86.4639.774 × 1019351854.5925.321 × 10210.00312.97913.37411.18110.03415.0061.736NGC 4085GiP4
J120658.79+524126.51.3182.532 × 1022112931.8142.170 × 10240.08517.38616.23312.88711.26818.2211.933SDSS J120658.79+524126.4Sy11
J122615.87+482937.71.4833.472 × 101853937.7281.263 × 10200.00116.99616.97914.45612.07017.4191.561SDSS J122615.72+482938.6H2G5
J123803.78+461819.90.4511.187 × 102247285.6351.246 × 10240.09918.20216.85414.35012.21318.3312.021SDSS J123803.75+461820.1EmG6
J124505.63+561430.20.6447.668 × 10166134.2637.300 × 10170.000217.98217.12815.27314.43118.4940.979J124505.6+561430QSO7
J130708.14+461834.22.2033.254 × 1021105655.4861.560 × 10230.02515.37615.74613.27711.34015.8841.6812MASX J13070814+4618346AGN8
J131048.89+544111.50.3353.493 × 102116538.9281.722 × 10230.06415.93715.82514.44113.35416.6021.6932MASX J13104883+5441104Sy21
J131459.85+511621.40.2684.042 × 10207364.1171.110 × 10220.02515.77416.37215.48214.23215.9641.4392MASX J13145978+5116213LIN AGN1
J131506.48+461409.70.4531.180 × 102225132.7246.547 × 10230.09815.78915.44413.99712.89917.1421.7442MASS J13150651+4614101Sy19
J131856.77+474609.61.2595.009 × 102155490.3032.208 × 10230.04015.86316.17713.82412.03916.3761.6442MASX J13185676+4746092G1
J134649.72+542621.70.4693.741 × 102014262.2381.138 × 10220.01815.88915.45014.53613.51518.5041.483
J140411.25+542519.61.1662.358 × 101821906.9504.432 × 10190.00117.37717.47115.07313.04918.2891.274SDSS J140411.24+542518.6H2G5
J140956.92+545649.50.7291.142 × 102237215.4665.832 × 10230.07717.63417.09914.38012.47317.4671.708LEDA 2816082EmG/H2G5, 6
J141424.94+465348.40.3462.227 × 102227003.1921.740 × 10240.15015.67915.37114.23012.82117.1471.8932MASSJ14142489+4653485Sy110
J143551.22+501118.20.4987.584 × 102119785.4773.012 × 10230.07617.13317.49515.30713.15917.2691.599SDSS J143551.27+501118.3G1
J144632.38+490038.20.5811.987 × 102192241.2233.153 × 10230.03716.53616.62713.54611.48816.8202.201SDSS J144632.39+490038.0Comp11
J145757.90+565323.80.8064.262 × 102151843.2542.740 × 10230.04617.27617.19014.29612.11317.6181.808

Note. (1) Source ID in the LoTSS value-added catalogue, (2) the total, integrated Stokes I 150 MHz radio flux density measured by the LOFAR HBA, (3) radio luminosity derived from LoTSS 150 MHz radio flux density, (4) WISE W4 profile-fit photometry flux, (5) calculated WISE 22μ luminosity, (6) Best available redshift estimate, (7)−(10) WISE W1−W4 profile-fit photometry magnitude, (11) Pan-STARSS forced aperture AB magnitude, (12) IRC parameter between the WISE 22 |$\mu$|m and LoTSS 150 MHz flux density, (13) Identifier archived in the SIMBAD catalogue, (14) Source type, G: Galaxy; Sy2: Seyfert 2 Galaxy; Sy1: Seyfert 1 Galaxy; Comp: Composite galaxy that likely have both star-forming and AGN activity; GiC: Galaxy in cluster of galaxies; LIN AGN: LINER-type AGN; QSO: Quasar object; EmG: Emission-line galaxy; X: X-ray source. (15) References; 1: Toba et al. (2014); 2: Dong et al. (2012); 3: Lavaux & Hudson (2011); 4: Crook et al. (2007); 5: Shirazi & Brinchmann (2012); 6: Shim & Chary (2013); 7: Allen et al. (2011); 8: Liu et al. (2011); 9: Oh et al. (2015); 10: Zhou et al. (2006); 11: Brinchmann, Kunth & Durret (2008).

AGN are the most common type of sources in our high-|$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| value sample. The vast majority of AGN are RQ AGN with radio powers that are at least 3 orders of magnitude smaller than the much brighter RL AGN. The radio emission in these systems is potentially generated by several different processes, in addition to non-thermal emission from accretion and the ejection of low-powered jets, accretion disc coronal activity, free–free emission from photo-ionized gas, and star-formation on sub-kpc and kpc scales are also expected to be major sources of radio emission in this class of object (Padovani et al. 2015; Panessa et al. 2019). The majority of RQ AGN will typically follow the IRC but some galaxies such as those in our sub-samples can have high values of q when their IR emission is enhanced by the presence of dust, usually located in the central regions surrounding the AGN. In addition to AGN, there are also four high-|$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| value sources that are identified as emission-line galaxies which are possibly dominated by star-forming processes. Additionally, there is one source that is classified as a composite galaxies with both AGN and star-forming activity. As suggested by Roussel et al. (2003), Costagliola et al. (2016), and Galvin et al. (2016), radio continuum emission is a delayed tracer of star formation due to the long time-scales involved (∼107 yr) for electrons to diffuse within galaxies on kpc scales. For galaxies with a very recent onset of star formation, this can lead to higher values of q being measured. There are five sources identified as galaxies of one type or another in our sub-samples but two sources (ILT J134649.72+542621.7 and ILT J145757.90+56532) have very little presence in the literature and have not been well-studied to date.

3.3 Additional q-value distributions

In addition to |$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| that is derived from the WISE W4 22 |$\mu$|m and 150 MHz LOFAR flux densities, we can also calculate additional q-values for the high-q sub-samples using the W1, W2, and W3 WISE MIR bands. Within the 5790 z < 0.1 and 10577 0.1 < z < 0.2 sources that have WISE W4 photometric flux, there are 5775 and 10 559 of them that also have WISE channel 1, 2, and 3 data available. As suggested by Wright et al. (2014a,b), the waste heat of a Type III civilisation is expected to peak at 100−600 K, corresponding to an excess MIR emission that peaks between 5−30 |$\mu$|m. Since the radio flux densities used to measure the values of q are constant for each source, the q-value distribution curves for the sources reflect the shape of the source SEDs in the MIR. Fig. 4 shows the q-values for the high-|$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| value sources versus wavelengths of the four WISE channels. As shown in Fig. 4, our sources generally have an SED that rises towards longer MIR wavelengths, which is similar to the model of a typical spiral galaxy built by Silva et al. (1998) ∓ no explicit features of Type III civilisations can be seen in either sub-sample.

The SED reflected by the various q-values for the (a) 20 (z < 0.1) and (b) 1 (0.1 < z < 0.2) high-$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$ value sources.
Figure 4.

The SED reflected by the various q-values for the (a) 20 (z < 0.1) and (b) 1 (0.1 < z < 0.2) high-|$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| value sources.

4 CANDIDATE TYPEII/III CIVILISATIONS

In this section, we look at the MIR colours and optical properties of sources exhibiting high |$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| value (see Section 3.1) to look for potential hosts of Type III civilisations.

4.1 MIR colour distribution

In order to find sources with extreme MIR colours among our high-|$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| sample, we plot the WISE w1 ∓ w2 colour distribution over the two sub-samples in Fig. 5. The criteria of W1 ∓ W2 > mean + 2σ is applied on the sources, however, we note that these sources could be the upper tail of the colour distribution, rather than a new population. As shown in Fig. 5(a) and (b), there are nine high-|$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| sources that appear to be outliers in the top right-hand of each panel – eight with (z < 0.1) and one with (0.1 < z < 0.2). We investigate the optical properties of these sources in the next section.

The WISE MIR colour distribution of W1 ∓ W2 versus 22 $\mu$m luminosity. The red points represent the (a) 20 (z < 0.1) and (b) 1 (0.1 < z < 0.2) high-$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$ value sources. The blue lines in both panels denote the averaged values and ± 2σ away from it.
Figure 5.

The WISE MIR colour distribution of W1 ∓ W2 versus 22 |$\mu$|m luminosity. The red points represent the (a) 20 (z < 0.1) and (b) 1 (0.1 < z < 0.2) high-|$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| value sources. The blue lines in both panels denote the averaged values and ± 2σ away from it.

4.2 Optical versus MIR source characteristics

In Fig. 6, we plot the Pan-STARRS i-band monochromatic AB magnitude (see Tonry et al. 2012 for details) versus WISE W4 luminosity for the two sub-samples, including the nine extremely red high-|$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| sources. Within the 5790 z < 0.1 and 10577 0.1 < z < 0.2 sources that have WISE W4 photometric fluxes, there are 5708 and 10 541 of them that also have PanSTARRs i-band forced aperture photometry magnitude data available. The Pan-STARRS i-band is the most sensitive band of the survey, with a bandpass of 690−819 nm, and a capture cross-section of ∼8.5 m2e photon−1 at 1.2 airmass (Tonry et al. 2012).

The distributions of Pan-STARRS i-band forced aperture photometry magnitude. The red points represent the high-$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$ value sources. The yellow triangles represent extremely red sources within the high $q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$ sample. The blue lines in both panels denote the least square linear fit over the whole sample and values of ± 5σ away from it.
Figure 6.

The distributions of Pan-STARRS i-band forced aperture photometry magnitude. The red points represent the high-|$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| value sources. The yellow triangles represent extremely red sources within the high |$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| sample. The blue lines in both panels denote the least square linear fit over the whole sample and values of ± 5σ away from it.

As shown in Fig. 6, the optical magnitude of the sources correlates with increasing 22 |$\mu$|m flux density, the solid orange lines show the linear fits to 22 |$\mu$|m over the samples. We select sources with their Mi 5σ under the linear fit as the optical-faint sample. We note that these sources could be the optically faint tail of the sources rather than a new population. Of the nine sources that have relatively high values of q and are extremely red in the MIR, four of these are underluminous, all four sources have a redshift of <0.1.

We present the radio, MIR, and optical images of these four sources in Fig. 7, the field of view (FoV) of each image is 0.5 arcsec. The LoTSS 150 MHz images retrieved from the LoTSS-DR1 image cutout service3 have a pixel size of 1.5 arcsec. The WISE multicolour images were retrieved from the NASA/IPAC Infrared Science Archive4 and have a pixel size of 1.25 arcsec. The latter are combined images, summing by the 3.4 (blue), 4.6 (green), and 12 (red) |$\mu$|m images together. The Pan-STARRS optical images have a pixel size of 0.25 arcsec, and are combined images using the Pan-STARRS y (red),i (green), g (blue) filters – these images are retrieved from the Pan-STARRS Data Release 1 Image Cutout Server.5 All the sub-figures in Fig. 7 are centered on the target coordinates drawn from the LoTSS DR1 catalogue. As shown in this figure, the sources appear to be point-like and very red in the MIR. These four sources are therefore candidate Type III civilisations.

Multiband images for the four very red, high $q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$ value, and optically faint sources. The LoTSS HBA 150 MHz radio images (left-hand column) have a pixel size of 1.5 arcsec, the WISE (middle column) MIR multicolour images are combined by 3.4 (blue), 4.6(green), and 12 (red) $\mu$m images, and have a pixel size of 1.5 arcsec, the Pan-STARRS optical optical images (right-hand column) are combined by data from the Pan-STARRS y (red), i (green), g (blue) filters and have a pixel size of 0.25 arcsec. All the sub-figures have the same FoV of 50 arcsec and are centered on the target coordinates drawn from the LoTSS DR1 catalogue.
Figure 7.

Multiband images for the four very red, high |$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| value, and optically faint sources. The LoTSS HBA 150 MHz radio images (left-hand column) have a pixel size of 1.5 arcsec, the WISE (middle column) MIR multicolour images are combined by 3.4 (blue), 4.6(green), and 12 (red) |$\mu$|m images, and have a pixel size of 1.5 arcsec, the Pan-STARRS optical optical images (right-hand column) are combined by data from the Pan-STARRS y (red), i (green), g (blue) filters and have a pixel size of 0.25 arcsec. All the sub-figures have the same FoV of 50 arcsec and are centered on the target coordinates drawn from the LoTSS DR1 catalogue.

Among the four candidate sources, two sources were classified as natural sources in the literature. ILT J120658.79+524126.5 (SIMBAD identifier: SDSS J120658.79+524126.4), was identified as a Seyfert 1 galaxy in (Toba et al. 2014; source ID: SDSS J120658.79+524126.4) based on its location on the BPT diagram (Baldwin, Phillips & Terlevich 1981). Another source−−ILT J123803.78+461819.9 (SIMBAD identifier: SDSS J123803.75+461820.1) was classified as an emission-line galaxy in (Shim & Chary 2013; source ID: SDSS 588017109686943854) with Hα line width measurement of ≳500 Å. In addition, there are two sources (ILT J134649.72+542621.7 and ILT J145757.90+565323.8) that have little presence in the scientific literature. The faint optical emission of these two sources may be explained by the presence of dust that generates excessive MIR emission and causes extinction in the optical band but they deserve further study in the context of being Type II/III civilisations.

5 SUMMARY

We have created a sample of 16 367 sources selected from the LoTSS-DR1 value-added catalogue that includes radio and MIR flux density measurements, optical photomertry, and photometric redshifts. We calculated |$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| values for 16 367 sources with z < 0.2, in an attempt to identify galaxies in the local Universe displaying the characteristics of energy intensive Kardasehv Type III civilisations. In particular, we looked for outliers in the IRC, searching for sources with |$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| > the linear fits +1 and |$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| > 0.25 + log105, one of the signatures expected of Type III civilisations generating significant amount of waste heat emission. Twenty one sources with such large q-values were found in the sample, corresponding to candidate Type III civilisations with their MIR radiation ∼10 times enhanced by waste heat are present in the sample.

We searched for identifiers for the 21 high-|$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| value sources from the SIMBAD catalogue – 19 out of 21 of the sources are clearly identified as natural sources of cosmic emission, 2 sources have little presence in the scientific literature. The majority of the previously identified sources with high q-values are likely to be obscured AGN or star-forming ETG. All these high-|$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| value sources are radio-faint (L150 MHz ≲ 1023 W Hz−1) and have an SED similar to that of a typical spiral galaxy.

With the expectation of MIR excess caused by waste heat in this subsample of 21 sources, we looked for sources with SEDs peaking in the MIR. By doing a colour distribution analysis, within the 21 sources with high-|$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| value, we identified 9 sources with the most extreme MIR colours that is 2σ above the average. In addition to MIR-excess characteristics, we also looked for optical-faint features within the nine extremely red high-q sources. Amongst them, we found four outliers that showed significantly low optical/MIR luminosity ratio in the optical magnitude distribution, these sources have optical magnitudes that are 5σ fainter than the linear fitting of the whole sample. Two of these four appear to be reasonably well-studied, with characteristics that can be explained by natural astrophysical processes – they are unlikely to be Type III civilisations. The other two sources (ILT J134649.72+542621.7 and ILT J145757.90+565323.8) warrant further observations. In addition, some of the sources in the two sub-samples with less extreme MIR colours are still poorly observed and understood.

The construction of a sample of high-|$q^{{\rm 22~\mu m}}_{{\rm 150~MHz}}$| values and very red, optically faint sources not only provides a good starting point for surveying advanced civilisations but also serves as an important input to studies of IRC modelling, AGN dust obscuration, and galaxy evolution. We note that this very direct method of using extreme q-values to identify Kardashev Type III candidates can be extended to the full LoTSS survey area and also to sources located at z > 0.2. This approach would generate a sample with an estimated 1 million sources (about 100 orders of magnitude greater than the sample presented here) and can place very strong constraints on the incidence of Type III civilisations in the Universe.

ACKNOWLEDGEMENTS

This work is supported by National Natural Science Foundation of China grant no. 11988101. This research made use of astropy, a community-developed core python package for Astronomy (Astropy Collaboration et al. 2013, 2018).

DATA AVAILABILITY

The data underlying this article are available in the article.

Footnotes

1

The value-added catalogue is available online at https://lofar-surveys.org/releases.html

REFERENCES

Allen
J. T.
,
Hewett
P. C.
,
Maddox
N.
,
Richards
G. T.
,
Belokurov
V.
,
2011
,
MNRAS
,
410
,
860

Annis
J.
,
1999
,
J. Br. Interplanet. Soc.
,
52
,
33

Astropy Collaboration
et al. .,
2013
,
A&A
,
558
,
A33

Astropy Collaboration
et al. .,
2018
,
AJ
,
156
,
123

Baldwin
J. A.
,
Phillips
M. M.
,
Terlevich
R.
,
1981
,
PASP
,
93
,
5

Bradbury
R. J.
,
Cirkovic
M. M.
,
Dvorsky
G.
,
2011
,
J. Br. Interplanet. Soc.
,
64
,
156

Brinchmann
J.
,
Kunth
D.
,
Durret
F.
,
2008
,
A&A
,
485
,
657

Calistro Rivera
G.
et al. ,
2017
,
MNRAS
,
469
,
3468

Carrigan,

Richard
A. J.
,
2009
,
ApJ
,
698
,
2075

Chambers
K. C.
et al. ,
2016
,
preprint (arXiv:1612.05560)

Condon
J. J.
,
Cotton
W. D.
,
Broderick
J. J.
,
2002
,
AJ
,
124
,
675

Costagliola
F.
,
Herrero-Illana
R.
,
Lohfink
A.
,
Pérez-Torres
M.
,
Aalto
S.
,
Muller
S.
,
Alberdi
A.
,
2016
,
A&A
,
594
,
A114

Crook
A. C.
,
Huchra
J. P.
,
Martimbeau
N.
,
Masters
K. L.
,
Jarrett
T.
,
Macri
L. M.
,
2007
,
ApJ
,
655
,
790

Cutri
R. M.
et al. ,
2021
,
VizieR Online Data Catalog
, p.
II/328

Delhaize
J.
et al. ,
2017
,
A&A
,
602
,
A4

Dong
X.-B.
,
Ho
L. C.
,
Yuan
W.
,
Wang
T.-G.
,
Fan
X.
,
Zhou
H.
,
Jiang
N.
,
2012
,
ApJ
,
755
,
167

Donley
J. L.
,
Rieke
G. H.
,
Rigby
J. R.
,
Pérez-González
P. G.
,
2005
,
ApJ
,
634
,
169

Draine
B. T.
,
Li
A.
,
2007
,
ApJ
,
657
,
810

Duncan
K. J.
et al. ,
2019
,
A&A
,
622
,
A3

Dyson
F. J.
,
1960
,
Science
,
131
,
1667

Elbaz
D.
,
Cesarsky
C. J.
,
Chanial
P.
,
Aussel
H.
,
Franceschini
A.
,
Fadda
D.
,
Chary
R. R.
,
2002
,
A&A
,
384
,
848

Galvin
T. J.
,
Seymour
N.
,
Filipović
M. D.
,
Tothill
N. F. H.
,
Marvil
J.
,
Drouart
G.
,
Symeonidis
M.
,
Huynh
M. T.
,
2016
,
MNRAS
,
461
,
825

Garrett
M. A.
,
2002
,
A&A
,
384
,
L19

Garrett
M. A.
,
2015
,
A&A
,
581
,
L5

Griffith
R. L.
,
Wright
J. T.
,
Maldonado
J.
,
Povich
M. S.
,
Sigurđsson
S.
,
Mullan
B.
,
2015
,
ApJS
,
217
,
25

Hill
G. J.
et al. ,
2008
,
ASP Conf. Ser. Vol. 4, The Hobby-Eberly Telescope Dark Energy Experiment (HETDEX): Description and Early Pilot Survey Results
.
Astron. Soc. Pac
,
San Francisco
, p.
115

Ilbert
O.
et al. ,
2013
,
A&A
,
556
,
A55

Kardashev
N. S.
,
1964
,
Soviet Ast.
,
8
,
217

Lavaux
G.
,
Hudson
M. J.
,
2011
,
MNRAS
,
416
,
2840

Liu
X.
,
Shen
Y.
,
Strauss
M. A.
,
Hao
L.
,
2011
,
ApJ
,
737
,
101

Magnelli
B.
et al. ,
2015
,
A&A
,
573
,
A45

Mohan
N.
,
Rafferty
D.
,
2015
,
Astrophysics Source Code Library, record ascl:1502.007

Nyland
K.
et al. ,
2017
,
MNRAS
,
464
,
1029

Oh
K.
,
Yi
S. K.
,
Schawinski
K.
,
Koss
M.
,
Trakhtenbrot
B.
,
Soto
K.
,
2015
,
ApJS
,
219
,
1

Padovani
P.
,
Bonzini
M.
,
Kellermann
K. I.
,
Miller
N.
,
Mainieri
V.
,
Tozzi
P.
,
2015
,
MNRAS
,
452
,
1263

Panessa
F.
,
Baldi
R. D.
,
Laor
A.
,
Padovani
P.
,
Behar
E.
,
McHardy
I.
,
2019
,
Nature Astron.
,
3
,
387

Park
S. Q.
et al. ,
2008
,
ApJ
,
678
,
744

Read
S. C.
et al. ,
2018
,
MNRAS
,
480
,
5625

Roussel
H.
,
Helou
G.
,
Beck
R.
,
Condon
J. J.
,
Bosma
A.
,
Matthews
K.
,
Jarrett
T. H.
,
2003
,
ApJ
,
593
,
733

Sargent
M. T.
,
Béthermin
M.
,
Daddi
E.
,
Elbaz
D.
,
2012
,
ApJ
,
747
,
L31

Shim
H.
,
Chary
R.-R.
,
2013
,
ApJ
,
765
,
26

Shimwell
T. W.
et al. ,
2017
,
A&A
,
598
,
A104

Shimwell
T. W.
et al. ,
2019
,
A&A
,
622
,
A1

Shirazi
M.
,
Brinchmann
J.
,
2012
,
MNRAS
,
421
,
1043

Silva
L.
,
Granato
G. L.
,
Bressan
A.
,
Danese
L.
,
1998
,
ApJ
,
509
,
103

Solarz
A.
,
Pollo
A.
,
Bilicki
M.
,
Pȩpiak
A.
,
Takeuchi
T. T.
,
Pia̧tek
P.
,
2019
,
PASJ
,
71
,
28

Toba
Y.
et al. ,
2014
,
ApJ
,
788
,
45

Tonry
J. L.
et al. ,
2012
,
ApJ
,
750
,
99

van Haarlem
M. P.
et al. ,
2013
,
A&A
,
556
,
A2

Wang
L.
et al. ,
2019
,
A&A
,
631
,
A109

Williams
W. L.
et al. ,
2019
,
A&A
,
622
,
A2

Wright
J. T.
,
Mullan
B.
,
Sigurdsson
S.
,
Povich
M. S.
,
2014a
,
ApJ
,
792
,
26

Wright
J. T.
,
Griffith
R. L.
,
Sigurdsson
S.
,
Povich
M. S.
,
Mullan
B.
,
2014b
,
ApJ
,
792
,
27

Yun
M. S.
,
Reddy
N. A.
,
Condon
J. J.
,
2001
,
ApJ
,
554
,
803

Zhou
H.
,
Wang
T.
,
Yuan
W.
,
Lu
H.
,
Dong
X.
,
Wang
J.
,
Lu
Y.
,
2006
,
ApJS
,
166
,
128

This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://academic.oup.com/journals/pages/open_access/funder_policies/chorus/standard_publication_model)