Low-frequency Sky Research Articles

Impact of extragalactic point sources on the low-frequency sky spectrum and cosmic dawn global 21-cm measurements

Abstract Contribution of resolved and unresolved extragalactic point sources to the low-frequency sky spectrum is a potentially non-negligible part of the astrophysical foregrounds for cosmic dawn 21-cm experiments. The clustering of such point sources on the sky, combined with the frequency-dependence of the antenna beam, can also make this contribution chromatic. By combining low-frequency measurements of the luminosity function and the angular correlation function of extragalactic point sources, we develop a model for the contribution of these sources to the low-frequency sky spectrum. Using this model, we find that the contribution of sources with flux density >10−6 Jy to the sky-averaged spectrum is smooth and of the order of a few kelvins at 50–200 MHz. We combine this model with measurements of the galactic foreground spectrum and weigh the resultant sky by the beam directivity of the conical log-spiral antenna planned as part of the Radio Experiment for the Analysis of Cosmic Hydrogen (REACH) project. We find that the contribution of point sources to the resultant spectrum is $\sim 0.4\%$ of the total foregrounds, but still larger by at least an order of magnitude than the standard predictions for the cosmological 21-cm signal. As a result, not accounting for the point-source contribution leads to a systematic bias in 21-cm signal recovery. We show, however, that in the REACH case, this reconstruction bias can be removed by modelling the point-source contribution as a power law with a running spectral index. We make our code publicly available as a Python package labelled epspy.

Lower Ionospheric Remote Sensing Using Low-Frequency Sky Waves

Lower Ionospheric Remote Sensing Using Low-Frequency Sky Waves

A Model Estimator for Noisy Compact Emission Recovery in Radio Synthesis Imaging

Reconstruction of a noisy compact emission must consider not only the point-spread function but also the effect of noise. However, the traditional threshold method in widely-used CLEAN-based algorithms finds it difficult to effectively prevent noise in the model image during noisy compact-emission reconstruction. This significantly limits the performance in noisy compact-emission reconstruction, such as deep field imaging. There are two major difficulties in the accurate reconstruction of a Stokes-I image of compact emission: first, the threshold method that has been used in practice is difficult to use to separate compact emission and noise; and second, over-subtraction makes it difficult for the reconstructed Stokes-I model image to remain positive. Therefore, a filter-based denoizing mechanism is introduced in the search phase of the model components to separate signal and noise so that the signal can be effectively extracted. The relatively larger loop gain for positive components means that the reconstructed model is in line with astrophysics. This will reduce the errors between the true sky image and the model image. The new model estimator is tested on a simulated JVLA observation with realistic source distributions from the VLA Low-Frequency Sky Survey project and the SKADS/SCubed simulation. The experiments show that it is very effective when used to separate signal and noise to lower the noise in the model image. This work explores the use of existing common software CASA to achieve high dynamic range imaging, which is an important step toward square kilometer array data processing.

Constraining a Model of the Radio Sky below 6 MHz Using the Parker Solar Probe/FIELDS Instrument in Preparation for Upcoming Lunar-based Experiments

We present a Bayesian analysis of data from the FIELDS instrument on board the Parker Solar Probe (PSP) spacecraft with the aim of constraining low-frequency (≲6 MHz) sky in preparation for several upcoming lunar-based experiments. We utilize data recorded during PSP’s coning roll maneuvers, in which the axis of the spacecraft is pointed 45° off of the Sun. The spacecraft then rotates about a line between the Sun and the spacecraft with a period of 24 minutes. We reduce the data into two formats: roll-averaged, in which the spectra are averaged over the roll, and phase-binned, in which the spectra are binned according to the phase of the roll. We construct a forward model of the FIELDS observations that includes numerical simulations of the antenna beam, an analytic emissivity function of the galaxy, and estimates of the absorption due to free electrons. Fitting 5 parameters, we find that the roll-averaged data can be fit well by this model, and we obtain posterior parameter constraints that are in general agreement with previous estimates. The model is not, however, able to fit the phase-binned data well, likely due to limitations such as the lack of nonsmooth emission structure at both small and large scales, enforced symmetry between the northern and southern galactic hemispheres, and large uncertainties in the free electron density. This suggests that significant improvement in the low-frequency sky model is needed in order to fully and accurately represent the sky at frequencies below 6 MHz.

The Southern-sky MWA Rapid Two-metre (SMART) pulsar survey—I. Survey design and processing pipeline

Abstract We present an overview of the Southern-sky MWA Rapid Two-metre (SMART) pulsar survey that exploits the Murchison Widefield Array’s large field of view and voltage-capture system to survey the sky south of 30 $^{\circ}$ in declination for pulsars and fast transients in the 140–170 MHz band. The survey is enabled by the advent of the Phase II MWA’s compact configuration, which offers an enormous efficiency in beam-forming and processing costs, thereby making an all-sky survey of this magnitude tractable with the MWA. Even with the long dwell times employed for the survey (4800 s), data collection can be completed in $<$ 100 h of telescope time, while still retaining the ability to reach a limiting sensitivity of $\sim$ 2–3 mJy (at 150 MHz, near zenith), which is effectively 3–5 times deeper than the previous-generation low-frequency southern-sky pulsar survey, completed in the 1990s. Each observation is processed to generate $\sim$ 5000–8000 tied-array beams that tessellate the full $\sim 610\, {\textrm{deg}^{2}}$ field of view (at 155 MHz), which are then processed to search for pulsars. The voltage-capture recording of the survey also allows a multitude of post hoc processing options including the reprocessing of data for higher time resolution and even exploring image-based techniques for pulsar candidate identification. Due to the substantial computational cost in pulsar searches at low frequencies, the survey data processing is undertaken in multiple passes: in the first pass, a shallow survey is performed, where 10 min of each observation is processed, reaching about one-third of the full-search sensitivity. Here we present the system overview including details of ongoing processing and initial results. Further details including first pulsar discoveries and a census of low-frequency detections are presented in a companion paper. Future plans include deeper searches to reach the full sensitivity and acceleration searches to target binary and millisecond pulsars. Our simulation analysis forecasts $\sim$ 300 new pulsars upon the completion of full processing. The SMART survey will also generate a complete digital record of the low-frequency sky, which will serve as a valuable reference for future pulsar searches planned with the low-frequency Square Kilometre Array.

Deep sub-arcsecond wide-field imaging of the Lockman Hole field at 144 MHz

High quality low-frequency radio surveys have the promise of advancing our understanding of many important topics in astrophysics, including the life cycle of active galactic nuclei (AGN), particle acceleration processes in jets, the history of star formation, and exoplanet magnetospheres. Currently leading low-frequency surveys reach an angular resolution of a few arcseconds. However, this resolution is not yet sufficient to study the more compact and distant sources in detail. Sub-arcsecond resolution is therefore the next milestone in advancing these fields. The biggest challenge at low radio frequencies is the ionosphere. If not adequately corrected for, ionospheric seeing blurs the images to arcsecond or even arcminute scales. Additionally, the required image size to map the degree-scale field of view of low-frequency radio telescopes at this resolution is far greater than what typical soft- and hardware is currently capable of handling. Here we present for the first time (to the best of our knowledge) widefield sub-arcsecond imaging at low radio frequencies. We derive ionospheric corrections in a few dozen individual directions and apply those during imaging efficiently using a recently developed imaging algorithm (arXiv:1407.1943, arXiv:1909.07226). We demonstrate our method by applying it to an eight hour observation of the International LOw Frequency ARray (LOFAR) Telescope (ILT) (arXiv:1305.3550). Doing so we have made a sensitive $7.4\ \mathrm{deg}^2$ $144\ \mathrm{MHz}$ map at a resolution of $0.3''$ reaching $25\ \mu\mathrm{Jy\ beam}^{-1}$ near the phase centre. The estimated $250,000$ core hours used to produce this image, fit comfortably in the budget of available computing facilities. This result will enable future mapping of the entire northern low-frequency sky at sub-arcsecond resolution.

Imaging the southern sky at 159 MHz using spherical harmonics with the engineering development array 2

Abstract One of the major priorities of international radio astronomy is to study the early universe through the detection of the 21 cm HI line from the epoch of reionisation (EoR). Due to the weak nature of the 21 cm signal, an important part in the detection of the EoR is removing contaminating foregrounds from our observations as they are multiple orders of magnitude brighter. In order to achieve this, sky maps spanning a wide range of frequencies and angular scales are required for calibration and foreground subtraction. Complementing the existing low-frequency sky maps, we have constructed a Southern Sky map through spherical harmonic transit interferometry utilising the Engineering Development Array 2 (EDA2), a Square Kilometre Array (SKA) low-frequency array prototype system. We use the m-mode formalism to create an all-sky map at 159 MHz with an angular resolution of 3 degrees, with data from the EDA2 providing information over +60 degrees to –90 degrees in declination. We also introduce a new method for visualising and quantifying how the baseline distribution of an interferometer maps to the spherical harmonics and discuss how prior information can be used to constrain spherical harmonic components that the interferometer is not sensitive to.

Pilot Study and Early Results of the Cosmic Filaments and Magnetism Survey with Nenufar: The Coma Cluster Field

NenuFAR, the New Extension in Nancay Upgrading LOFAR, is currently in its early science phase. It is in this context that the Cosmic Filaments and Magnetism Pilot Survey is observing sources with the array as it is still under construction—with 57 (56 core, 1 distant) out of a total planned 102 (96 core, 6 distant) mini-arrays online at the time of observation—to get a first look at the low-frequency sky with NenuFAR. One of its targets is the Coma galaxy cluster: a well-known object, host of the prototype radio halo. It also hosts other features of scientific import, including a radio relic, along with a bridge of emission connecting it with the halo. It is thus a well-studied object.In this paper, we show the first confirmed NenuFAR detection of the radio halo and radio relic of the Coma cluster at 34.4 MHz, with associated intrinsic flux density estimates: we find an integrated flux value of 106.3 ± 3.5 Jy for the radio halo, and 102.0 ± 7.4 Jy for the radio relic. These are upper bound values, as they do not include point-source subtraction. We also give an explanation of the technical difficulties encountered in reducing the data, along with steps taken to resolve them. This will be helpful for other scientific projects which will aim to make use of standalone NenuFAR imaging observations in the future.

Thermal radio absorption as a tracer of the interaction of SNRs with their environments

We present new images and continuum spectral analysis for 14 resolved Galactic supernova remnants (SNRs) selected from the 74 MHz Very Large Array Low-Frequency Sky Survey Redux (VLSSr). We combine new integrated measurements from the VLSSr with, when available, flux densities extracted from the Galactic and Extragalactic All-Sky Murchison Widefield Array Survey and measurements from the literature to generate improved integrated continuum spectra sampled from ~15 MHz to ~217 GHz. We present the VLSSr images. When possible we combine them with publicly available images at 1.4 GHz, to analyse the resolved morphology and spectral index distribution across each SNR. We interpret the results and look for evidence of thermal absorption caused by ionised gas either proximate to the SNR itself, or along its line of sight. Three of the SNRs, G4.5+6.8 (Kepler), G28.6−0.1, and G120.1+1.4 (Tycho), have integrated spectra which can be adequately fit with simple power laws. The resolved spectral index map for Tycho confirms internal absorption which was previously detected by the Low Frequency Array, but it is insufficient to affect the fit to the integrated spectrum. Two of the SNRs are pulsar wind nebulae, G21.5−0.9 and G130.7+3.1 (3C 58). For those we identify high-frequency spectral breaks at 38 and 12 GHz, respectively. For the integrated spectra of the remaining nine SNRs, a low frequency spectral turnover is necessary to adequately fit the data. In all cases we are able to explain the turnover by extrinsic thermal absorption. For G18.8+0.3 (Kes 67), G21.8−0.6 (Kes 69), G29.7−0.3 (Kes 75), and G41.1−0.3 (3C 397), we attribute the absorption to ionised gas along the line of sight, possibly from extended H II region envelopes. For G23.3−0.3 (W41) the absorption can be attributed to H II regions located in its immediate proximity. Thermal absorption from interactions at the ionised interface between SNR forward shocks and the surrounding medium were previously identified as responsible for the low frequency turnover in SNR G31.9+0.0 (3C 391); our integrated spectrum is consistent with the previous results. We present evidence for the same phenomenon in three additional SNRs G27.4+0.0 (Kes 73), G39.2–0.3 (3C 396), and G43.3–0.2 (W49B), and derive constraints on the physical properties of the interaction. This result indicates that interactions between SNRs and their environs should be readily detectable through thermal absorption by future low frequency observations of SNRs with improved sensitivity and resolution.

The LOFAR Two-meter Sky Survey: Deep Fields Data Release 1

The Low Frequency Array (LOFAR) is an ideal instrument to conduct deep extragalactic surveys. It has a large field of view and is sensitive to large-scale and compact emission. It is, however, very challenging to synthesize thermal noise limited maps at full resolution, mainly because of the complexity of the low-frequency sky and the direction dependent effects (phased array beams and ionosphere). In this first paper of a series, we present a new calibration and imaging pipeline that aims at producing high fidelity, high dynamic range images with LOFAR High Band Antenna data, while being computationally efficient and robust against the absorption of unmodeled radio emission. We apply this calibration and imaging strategy to synthesize deep images of the Boötes and Lockman Hole fields at ~150 MHz, totaling ~80 and ~100 h of integration, respectively, and reaching unprecedented noise levels at these low frequencies of ≲30 and ≲23μJy beam−1in the inner ~3 deg2. This approach is also being used to reduce the LOTSS-wide data for the second data release.

Cassiopeia A, Cygnus A, Taurus A, and Virgo A at ultra-low radio frequencies

Context. The four persistent radio sources in the northern sky with the highest flux density at metre wavelengths are Cassiopeia A, Cygnus A, Taurus A, and Virgo A; collectively they are called the A-team. Their flux densities at ultra-low frequencies (< 100 MHz) can reach several thousands of janskys, and they often contaminate observations of the low-frequency sky by interfering with image processing. Furthermore, these sources are foreground objects for all-sky observations hampering the study of faint signals, such as the cosmological 21 cm line from the epoch of reionisation. Aims. We aim to produce robust models for the surface brightness emission as a function of frequency for the A-team sources at ultra-low frequencies. These models are needed for the calibration and imaging of wide-area surveys of the sky with low-frequency interferometers. This requires obtaining images at an angular resolution better than 15″ with a high dynamic range and good image fidelity. Methods. We observed the A-team with the Low Frequency Array (LOFAR) at frequencies between 30 MHz and 77 MHz using the Low Band Antenna system. We reduced the datasets and obtained an image for each A-team source. Results. The paper presents the best models to date for the sources Cassiopeia A, Cygnus A, Taurus A, and Virgo A between 30 MHz and 77 MHz. We were able to obtain the aimed resolution and dynamic range in all cases. Owing to its compactness and complexity, observations with the long baselines of the International LOFAR Telescope will be required to improve the source model for Cygnus A further.

Ultimate frontier of 21-cm cosmology

We present the most detailed computation to date of the 21-cm global signal and fluctuations at $z\gtrsim 500$. Our calculations include a highly precise estimate of the Wouthuysen-Field (WF) effect and the first explicit calculation of the impact of free-free processes, the two dominant components of the signal at $z\gtrsim 800$. We implement a new high-resolution Ly$\alpha$ radiative transfer calculation, coupled to a state-of-the-art primordial recombination code. Using these tools, we find a global signal from 21-cm processes alone of roughly 0.01mK at $z\sim1000$, slightly larger than it would be without the WF effect, but much weaker than previous estimates including this effect. We also find that this signal is swamped by a smooth $1-2$ mK signal due to free-free absorption at high redshift by the partially ionized gas along the line of sight. In addition, we estimate the amplitude of 21-cm fluctuations, of order $\sim 10^{-7}$ mK at $z\sim1000$. Unfortunately, we find that due to the brightness of the low-frequency sky, these fluctuations will not be observable beyond $z\sim$ a few hundred by even extremely futuristic observations. The 21 cm fluctuations are exponentially suppressed at higher redshifts by the large free-free optical depth, making this the ultimate upper redshift limit for 21-cm surveys.

G51.04+0.07 and its environment: Identification of a new Galactic supernova remnant at low radio frequencies

We have identified a new supernova remnant (SNR), G51.04+0.07, using observations at 74 MHz from the Very Large Array Low-Frequency Sky Survey Redux (VLSSr). Earlier, higher frequency radio continuum, recombination line, and infrared data had correctly inferred the presence of nonthermal radio emission within a larger, complex environment including ionised nebulae and active star formation. However, our observations have allowed us to redefine at least one SNR as a relatively small source (7.′5 × 3′in size) located at the southern periphery of the originally defined SNR candidate G51.21+0.11. The integrated flux density of G51.04+0.07 at 74 MHz is 6.1 ± 0.8 Jy, while its radio continuum spectrum has a slope α = −0.52 ± 0.05 (S v ∝ vα), typical of a shell-type remnant. We also measured spatial variations in the spectral index between 74 and 1400 MHz across the source, ranging from a steeper spectrum (α = −0.50 ± 0.04) coincident with the brightest emission to a flatter component (α = −0.30 ± 0.07) in the surrounding fainter region. To probe the interstellar medium into which the redefined SNR is likely evolving, we have analysed the surrounding atomic and molecular gas using the 21 cm neutral hydrogen (HI) and 13CO(J = 1 − 0) emissions. We found that G51.04+0.07 is confined within an elongated HI cavity and that its radio emission is consistent with the remains of a stellar explosion that occurred ~6300 yr ago at a distance of 7.7 ± 2.3 kpc. Kinematic data suggest that the newly discovered SNR lies in front of HII regions in the complex, consistent with the lack of a turnover in the low frequency continuum spectrum. The CO observations revealed molecular material that traces the central and northern parts of G51.04+0.07. The interaction between the cloud and the radio source is not conclusive and motivates further study. The relatively low flux density (~1.5 Jy at 1400 MHz) of G51.04+0.07 is consistent with this and many similar SNRs lying hidden along complex lines of sight towards inner Galactic emission complexes. It would also not be surprising if the larger complex studied here hosted additional SNRs.

The effect of the ionosphere on ultra-low-frequency radio-interferometric observations

Context. The ionosphere is the main driver of a series of systematic effects that limit our ability to explore the low-frequency (<1 GHz) sky with radio interferometers. Its effects become increasingly important towards lower frequencies and are particularly hard to calibrate in the low signal-to-noise ratio (S/N) regime in which low-frequency telescopes operate. Aims. In this paper we characterise and quantify the effect of ionospheric-induced systematic errors on astronomical interferometric radio observations at ultra-low frequencies (<100 MHz). We also provide guidelines for observations and data reduction at these frequencies with the LOw Frequency ARray (LOFAR) and future instruments such as the Square Kilometre Array (SKA). Methods. We derive the expected systematic error induced by the ionosphere. We compare our predictions with data from the Low Band Antenna (LBA) system of LOFAR. Results. We show that we can isolate the ionospheric effect in LOFAR LBA data and that our results are compatible with satellite measurements, providing an independent way to measure the ionospheric total electron content (TEC). We show how the ionosphere also corrupts the correlated amplitudes through scintillations. We report values of the ionospheric structure function in line with the literature. Conclusions. The systematic errors on the phases of LOFAR LBA data can be accurately modelled as a sum of four effects (clock, ionosphere first, second, and third order). This greatly reduces the number of required calibration parameters, and therefore enables new efficient calibration strategies.

A Molecular Line Survey around Orion at Low Frequencies with the MWA

Abstract The low-frequency sky may reveal some of the secrets yet to be discovered. Until recently, molecules had never been detected within interstellar clouds at frequencies below 700 MHz. Following the pilot survey toward the Galactic center at 103–133 MHz with the Murchison Widefield Array, we surveyed 400 deg2 centered on the Orion KL nebula from 99 to 170 MHz. Orion is a nearby region of active star formation and known to be a chemically rich environment. In this paper, we present tentative detections of nitric oxide and its isotopologues, singularly deuterated formic acid, molecular oxygen, and several unidentified transitions. The three identified molecules are particularly interesting, as laboratory experiments have suggested that these molecules are precursors to the formation of amines.

Redundant Calibration: breaking the constraints of limited sky information

AbstractThe latest generation of low frequency radio interferometers, e.g. LOFAR, MWA, PAPER, has been pushing down the detection limits on the hydrogen signal from the Epoch of Reionisation. However, due to the challenges posed by foregrounds and instrumental systematics the signal has eluded detection thus far. To overcome these challenges we require a detailed understanding of the calibration of these relatively new telescopes. This led to a renewed interest in redundant calibration. Classical calibration schemes depend on sky models based on limited knowledge of the low frequency sky. Redundant calibration, however, allows us to escape our ignorance as it is sky model independent. We will review the field of redundant calibration, and present work we have undertaken to understand the limitations of this calibration method.

The LWA1 Low Frequency Sky Survey

We present a survey of the radio sky accessible from the first station of the Long Wavelength Array (LWA1). Images are presented at nine frequencies between 35 and 80 MHz with spatial resolutions ranging from $4.7^\circ$ to $2.0^\circ$, respectively. The maps cover the sky north of a declination of $-40^\circ$ and represent the most modern systematic survey of the diffuse Galactic emission within this frequency range. We also combine our survey with other low frequency sky maps to create an updated model of the low frequency sky. Due to the low frequencies probed by our survey, the updated model better accounts for the effects of free-free absorption from Galactic ionized Hydrogen. A longer term motivation behind this survey is to understand the foreground emission that obscures the redshifted 21 cm transition of neutral hydrogen from the cosmic dark ages ($z>10$) and, at higher frequencies, the epoch of reionization ($z>6$).

INTEGRATED RADIO CONTINUUM SPECTRA OF GALAXIES

We investigate the spectral shape of the total continuum radiation, between 74 MHz and 5 GHz (400 to 6 cm in wavelength), for a large sample of bright galaxies. We take advantage of the overlapping survey coverage of the VLA Low-Frequency Sky Survey, the Westerbork Northern Sky Survey, the NRAO VLA Sky Survey, and the Green Bank 6 cm survey to achieve significantly better resolution, sensitivity, and sample size compared to prior efforts of this nature. For our sample of 250 bright galaxies we measure a mean spectral index, ${\alpha}$, of -0.69 between 1.4 and 4.85 GHz, -0.55 between 325 MHz and 1.4 GHz, and -0.45 between 74 and 325 MHz, which amounts to a detection of curvature in the mean spectrum. The magnitude of this curvature is approximately ${\Delta \alpha}$ = -0.2 per logarithmic frequency decade when fit with a generalized function having constant curvature. No trend in low frequency spectral flattening versus galaxy inclination is evident in our data, suggesting that free-free absorption is not a satisfying explanation for the observed curvature. The ratio of thermal to non-thermal emission is estimated by two independent methods, (1) using the IRAS far-IR fluxes, and (2) with the value of the total spectral index. Method (1) results in a distribution of 1.4 GHz thermal fractions of 9% ${\pm}$ 3%, which is consistent with previous studies, while method (2) produces a mean 1.4 GHz thermal fraction of 51% with dispersion 26%. The highly implausible values produced by method (2) indicate that the sum of typical power-law thermal and non-thermal components is not a viable model for the total spectral index between 325 and 1.4 GHz.

The Very Large Array Low-frequency Sky Survey Redux (VLSSr)

We present the results of a recent re-reduction of the data from the Very Large Array (VLA) Low-frequency Sky Survey (VLSS). We used the VLSS catalog as a sky model to correct the ionospheric distortions in the data and create a new set of sky maps and corresponding catalog at 73.8 MHz. The VLSS Redux (VLSSr) has a resolution of 75 arcsec, and an average map RMS noise level of $\sigma\sim0.1$ Jy beam$^{-1}$. The clean bias is $0.66\times\sigma$, and the theoretical largest angular size is 36 arcmin. Six previously un-imaged fields are included in the VLSSr, which has an unbroken sky coverage over 9.3 sr above an irregular southern boundary. The final catalog includes 92,964 sources. The VLSSr improves upon the original VLSS in a number of areas including imaging of large sources, image sensitivity, and clean bias; however the most critical improvement is the replacement of an inaccurate primary beam correction which caused source flux errors which vary as a function of radius to nearest pointing center in the VLSS.

Unveiling the nature of the unidentified gamma-ray sources: blazar counterparts at low radio frequencies

AbstractAbout one third of the gamma-ray sources detected by Fermi have still no firmly established counterpart at lower energies. Here we propose a new approach to find candidate counterparts for the unidentified gamma-ray sources (UGSs) based on the 325 MHz radio survey performed with Westerbork Synthesis Radio Telescope (WSRT) in the northern hemisphere. First we investigate the low-frequency radio properties of blazars, the largest known population of gamma-ray sources; then we search for sources with similar radio properties combining the information derived from the Westerbork Northern Sky Survey (WENSS) with those of the NRAO VLA Sky survey (NVSS). We present a list of candidate counterparts for 32 UGSs with at least one counterpart in the WENSS. We also performed an extensive research in literature to look for infrared and optical counterparts of the gamma-ray blazar candidates selected with the low-frequency radio observations to confirm their nature. On the basis of our multifrequency research we identify 23 new gamma-ray blazar candidates out of 32 UGSs investigated. I will also present the first analysis of very low frequency radio emission of blazars based on the recent Very Large Array Low-Frequency Sky Survey (VLSS) at 74 MHz. I show that blazars present radio flat spectra when evaluated at 74 MHz, about an order of magnitude in frequency lower than previous analyses. The implications of these findings in the contest of the blazars – radio galaxies connection will be discussed.