VLA Low-Frequency Sky Survey Research Articles

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.

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.

BLAZAR SPECTRAL PROPERTIES AT 74 MHz

Blazars are the most extreme class of active galactic nuclei (AGNs). Despite a previous investigation at 102 MHz for a small sample of BL Lacs and our recent analysis of blazars detected in the Westerbork Northern Sky Survey (WENSS), a systematic study of the blazar spectral properties at frequencies below 100 MHz has been never carried out. In this paper, we present the first analysis of the radio spectral behavior of blazars based on the recent Very Large Array Low-Frequency Sky Survey (VLSS) at 74 MHz. We search for blazar counterparts in the VLSS catalog confirming that they are detected at 74 MHz. We then show that blazars present radio flat spectra {(i.e., radio spectral indices ~0.5)} when evaluated also about an order of magnitude in frequency lower than previous analyses. Finally, we discuss the implications of our findings in the context of the blazars - radio galaxies connection since the low frequency radio data provide a new diagnostic tool to verify the expectations of the unification scenario for radio-loud active galaxies.

Climatology of midlatitude ionospheric disturbances from the Very Large Array Low‐frequency Sky Survey

The results of a climatological study of ionospheric disturbances derived from observations of cosmic sources from the Very Large Array (VLA) Low‐frequency Sky Survey (VLSS) are presented. We have used the ionospheric corrections applied to the 74 MHz interferometric data within the VLSS imaging process to obtain fluctuation spectra for the total electron content (TEC) gradient on spatial scales from a few to hundreds of kilometers and temporal scales from less than one minute to nearly an hour. The observations sample nearly all times of day and all seasons. They also span latitudes and longitudes from 28°N to 40°N and 95°W to 114°W, respectively. We have binned and averaged the fluctuation spectra according to time of day, season, and geomagnetic (Kp index) and solar (F10.7) activity. These spectra provide a detailed, multiscale account of seasonal and intraday variations in ionospheric activity with wavelike structures detected at wavelengths between about 35 and 250 km. In some cases, trends between spectral power and Kp index and/or F10.7 are also apparent. In addition, the VLSS observations allow for measurements of the turbulent power spectrum down to periods of 40 s (scales of ∼0.4 km at the height of the E‐region). While the level of turbulent activity does not appear to have a strong dependence on eitherKp index or F10.7, it does appear to be more pronounced during the winter daytime, summer nighttime, and near dusk during the spring.

VLSS redux: Software improvements applied to the Very Large Array Low‐Frequency Sky Survey

We present details of improvements to data processing and analysis which were recently used for a re‐reduction of the Very Large Array (VLA) Low‐Frequency Sky Survey (VLSS) data. Algorithms described are implemented in the data‐reduction package Obit, and include smart‐windowing to reduce clean bias, improved automatic radio frequency interference removal, improved bright‐source peeling, and higher‐order Zernike fits to model the ionospheric phase contributions. An additional, but less technical improvement was using the original VLSS catalog as a same‐frequency/same‐resolution reference for calculating ionospheric corrections, allowing more accuracy and a higher percentage of data for which solutions are found. We also discuss new algorithms for extracting a source catalog and analyzing ionospheric fluctuations present in the data. The improved reduction techniques led to substantial improvements including images of six previously unpublished fields (1% of the survey area) and reducing the clean bias by 50%. The largest angular size imaged has been roughly doubled, and the number of cataloged sources is increased by 35% to 95,000.

Radio galaxies in the Sloan Digital Sky Survey: spectral index-environment correlations

We analyze optical and radio properties of radiogalaxies detected in the Sloan Digital Sky Survey (SDSS). The sample of radio sources are selected from the catalogue of Kimball & Ivezi\'c (2008) with flux densities at 325, 1400 and 4850 MHz, using WENSS, NVSS and GB6 radio surveys and from flux measurements at 74 MHz taken from VLA Low-frequency Sky Survey \citep{cohen}. We study radiogalaxy spectral properties using radio colour-colour diagrams and find that our sample follows a single power law from 74 to 4850 MHz. The spectral index vs. spectroscopic redshift relation ($\alpha-z$) is not significant for our sample of radio sources. We analyze a subsample of radio sources associated with clusters of galaxies identified from the maxBCG catalogue and find that about 40% of radio sources with ultra steep spectra (USS, $\alpha<-1$, where $S_\nu \propto \nu^{\alpha}$) are associated with galaxy clusters or groups of galaxies. We construct a Hubble diagram of USS radio sources in the optical $r$ band up to $z\sim0$.8 and compare our results with those for normal galaxies selected from different optical surveys and find that USS radio sources are around as luminous as the central galaxies in the maxBCG cluster sample and typically more than 4 magnitudes brighter than normal galaxies at $z\sim0$.3. We study correlations between spectral index, richness and luminosity of clusters associated with radio sources. We find that USS at low redshift are rare, most of them reside in regions of unusually high ambient density, such of those found in rich cluster of galaxies. Our results also suggest that clusters of galaxies associated with steeper than the average spectra have higher richness counts and are populated by luminous galaxies in comparison with those environments associated to radio sources with flatter than the average spectra. A plausible explanation for our results is that radio emission is more pressure confined in higher gas density environments such as those found in rich clusters of galaxies and as a consequence radio lobes in rich galaxy clusters will expand adiabatically and lose energy via synchrotron and inverse Compton losses, resulting in a steeper radio spectra.

A BLIND SEARCH FOR MAGNETOSPHERIC EMISSIONS FROM PLANETARY COMPANIONS TO NEARBY SOLAR-TYPE STARS

This paper reports a blind search for magnetospheric emissions from planets around nearby stars. Young stars are likely to have much stronger stellar winds than the Sun, and because planetary magnetospheric emissions are powered by stellar winds, stronger stellar winds may enhance the radio luminosity of any orbiting planets. Using various stellar catalogs, we selected nearby stars (<~ 30 pc) with relatively young age estimates (< 3 Gyr). We constructed different samples from the stellar catalogs, finding between 100 and several hundred stars. We stacked images from the 74-MHz (4-m wavelength) VLA Low-frequency Sky Survey (VLSS), obtaining 3\sigma limits on planetary emission in the stacked images of between 10 and 33 mJy. These flux density limits correspond to average planetary luminosities less than 5--10 x 10^{23} erg/s. Using recent models for the scaling of stellar wind velocity, density, and magnetic field with stellar age, we estimate scaling factors for the strength of stellar winds, relative to the Sun, in our samples. The typical kinetic energy carried by the stellar winds in our samples is 15--50 times larger than that of the Sun, and the typical magnetic energy is 5--10 times larger. If we assume that every star is orbited by a Jupiter-like planet with a luminosity larger than that of the Jovian decametric radiation by the above factors, our limits on planetary luminosities from the stacking analysis are likely to be a factor of 10--100 above what would be required to detect the planets in a statistical sense. Similar statistical analyses with observations by future instruments, such as the Low Frequency Array (LOFAR) and the Long Wavelength Array (LWA), offer the promise of improvements by factors of 10--100.

PROBING FINE-SCALE IONOSPHERIC STRUCTURE WITH THE VERY LARGE ARRAY RADIO TELESCOPE

High resolution (~1 arcminute) astronomical imaging at low frequency (below 150 MHz) has only recently become practical with the development of new calibration algorithms for removing ionospheric distortions. In addition to opening a new window in observational astronomy, the process of calibrating the ionospheric distortions also probes ionospheric structure in an unprecedented way. Here we explore one aspect of this new type of ionospheric measurement, the differential refraction of celestial source pairs as a function of their angular separation. This measurement probes variations in the spatial gradient of the line-of-sight total electron content (TEC) to 0.001 TECU/km accuracy over spatial scales of under 10 km to over 100 km. We use data from the VLA Low-frequency Sky Survey (VLSS; Cohen et al. 2007, AJ 134, 1245), a nearly complete 74 MHz survey of the entire sky visible to the Very Large Array (VLA) telescope in Socorro, New Mexico. These data comprise over 500 hours of observations, all calibrated in a standard way. While ionospheric spatial structure varies greatly from one observation to the next, when analyzed over hundreds of hours, statistical patterns become apparent. We present a detailed characterization of how the median differential refraction depends on source pair separation, elevation and time of day. We find that elevation effects are large, but geometrically predictable and can be "removed" analytically using a "thin-shell" model of the ionosphere. We find significantly greater ionospheric spatial variations during the day than at night. These diurnal variations appear to affect the larger angular scales to a greater degree indicating that they come from disturbances on relatively larger spatial scales (100s of km, rather than 10s of km).

A deep 150 MHz Giant Metrewave Radio Telescope survey in Eridanus

We present results of a 150 MHz survey of a field centred on e Eridani, undertaken with the Giant Metrewave Radio Telescope (GMRT). The survey covers an area with a diameter of ∼2°, has a spatial resolution of ∼30 arcsec and a noise level of σ= 3.1 mJy at the pointing centre. These observations provide a deeper and higher resolution view of the 150 MHz radio sky than the 7C survey (although the 7C survey covers a much larger area). A total of 113 sources were detected, most are point-like, but 20 are extended. We present an analysis of these sources, in conjunction with the National Radio Astronomy Observatory (NRAO) Very Large Array (VLA) telescope Sky Survey (NVSS) (at 1.4 GHz) and VLA Low-frequency Sky Survey (VLSS) (at 74 MHz). This process allowed us to identify five Ultra-Steep Spectrum (USS) radio sources that are candidate high-redshift radio galaxies (HzRGs). In addition, we have derived the dN/dS distribution for these observations and compare our results with other low-frequency radio surveys.

Obit: A Development Environment for Astronomical Algorithms

ABSTRACTObit is a development environment for astronomical algorithms. It implements multiple “native” data formats including FITS and AIPS. The shared data model with AIPS means almost complete interoperability, allowing algorithmic development in a production environment. The Obit package supports both prepackaged, compiled tasks and a python interface to the major class functionality to allow rapid prototyping using python scripts. Developed features are predominantly radio astronomical, but support for astronomical images and tables is generic and should allow a wider range of applications. The design is intended to be readily extensible. An example is given that explores the effect of “CLEAN bias” in the VLA Low-Frequency Sky Survey and compares its method of CLEAN to an automated CLEAN window selection scheme that should reduce the bias. This technique shows much promise for reducing this systematic effect.

The VLA Low-Frequency Sky Survey

The VLA Low-frequency Sky Survey (VLSS) has imaged 95% of the 3π sr of sky north of δ = -30° at a frequency of 74 MHz (4 m wavelength). The resolution is 80'' (FWHM) throughout, and the typical rms noise level is ⟨σ⟩ ≈ 0.1 Jy beam-1. The typical point-source detection limit is 0.7 Jy beam-1, and so far nearly 70,000 sources have been cataloged. This survey used the 74 MHz system added to the VLA in 1998. It required new imaging algorithms to remove the large ionospheric distortions at this very low frequency throughout the entire ∼11.9° field of view. This paper describes the observation and data reduction methods used for the VLSS and presents the survey images and source catalog. All of the calibrated images and the source catalog are available from the VLSS Web site for use by the astronomical community.

Radio Frequency Spectra of 388 Bright 74 MHz Sources

As a service to the community, we have compiled radio frequency spectra from the literature for all sources within the VLA Low Frequency Sky Survey (VLSS) that are brighter than 15 Jy at 74 MHz. Over 160 references were used to maximize the amount of spectral data used in the compilation of the spectra, while also taking care to determine the corrections needed to put the flux densities from all reference on the same absolute flux density scale. With the new VLSS data, we are able to vastly improve upon previous efforts to compile spectra of bright radio sources to frequencies below 100 MHz because (1) the VLSS flux densities are more reliable than those from some previous low frequency surveys and (2) the VLSS covers a much larger area of the sky (declination >-30 deg.) than many other low frequency surveys (e.g., the 8C survey). In this paper, we discuss how the spectra were constructed and how parameters quantifying the shapes of the spectra were derived. Both the spectra and the shape parameters are made available here to assist in the calibration of observations made with current and future low frequency radio facilities.