Search

The MeerKAT Absorption Line Survey (MALS) has observed 391 telescope pointings at L-band (900 - 1670 MHz) at $\delta\lesssim$ $+20\deg$. We present radio continuum images and a catalog of 495,325 (240,321) radio sources detected at a signal-to-noise ratio (SNR) $>$5 over an area of 2289 deg$^2$ (1132 deg$^2$) at 1006 MHz (1381 MHz). Every MALS pointing contains a central bright radio source ($S_{1\,\mathrm{GHz}} \gtrsim 0.2$ Jy). The median spatial resolution is $12^{\prime\prime}$ ($8^{\prime\prime}$). The median rms noise away from the pointing center is 25 $\mu$Jy beam$^{-1}$ (22 $\mu$Jy beam$^{-1}$) and is within $\sim$ 15% of the achievable theoretical sensitivity. The flux density scale ratio and astrometric accuracy deduced from multiply observed sources in MALS are less than 1% (8% scatter) and $1^{\prime\prime}$, respectively. Through comparisons with NVSS and FIRST at 1.4 GHz, we establish the catalog's accuracy in the flux density scale and astrometry to be better than 6% (15% scatter) and $0.8^{\prime\prime}$, respectively. The median flux density offset is higher (9%) for an alternate beam model based on holographic measurements. The MALS radio source counts at 1.4 GHz are in agreement with literature. We estimate spectral indices ($\alpha$) of a subset of 125,621 sources (SNR$>$8), confirm the flattening of spectral indices with decreasing flux density and identify 140 ultra steep-spectrum ($\alpha<-1.3$) sources as prospective high-$z$ radio galaxies ($z>2$). We have identified 1308 variable and 122 transient radio sources comprising primarily of AGN that demonstrate long-term (26 years) variability in their observed flux densities. The MALS catalogs and images are publicly available at https://mals.iucaa.in., Comment: 64 pages, 25 figures, accepted for publication in the ApJS (full version of the paper with complete tables is available at DR1 release notes)

We present X-ray and radio observations of the recently-discovered bow shock pulsar wind nebula associated with PSR J0002+6216, characterizing the PWN morphology, which was unresolved in previous studies. The multi-frequency, multi-epoch Very Large Array radio observations reveal a cometary tail trailing the pulsar and extending up to 5.3', with multiple kinks along the emission. The presented radio continuum images from multi-configuration broadband VLA observations are one of the first results from the application of multi-term multi-frequency synthesis deconvolution in combination with the awproject gridder implemented in the Common Astronomy Software Applications package (CASA). The X-ray emission observed with Chandra extends to only 21'', fades quickly, and has some hot spots present along the extended radio emission. These kinks could indicate the presence of density variation in the local ISM or turbulence. The bow shock standoff distance estimates a small bow shock region with a size 0.003-0.009 pc, consistent with the pulsar spin-down power of Edot=1.51x10^35 ergs/s estimated from timing. The high-resolution radio image reveals the presence of an asymmetry in the bow shock region which is also present in the X-ray image. The broadband radio image shows an unusually steep spectrum along with a flat-spectrum sheath, which could indicate varying opacity or energy injection into the region. Spatially-resolved X-ray spectra provide marginal evidence of synchrotron cooling along the extended tail. Our analysis of the X-ray data also shows that this pulsar has a low spin-down power and one of the lowest X-ray efficiencies observed in these objects., Comment: 15 pages, 7 figures Accepted for publication in ApJ

We present details of the Automated Radio Telescope Imaging Pipeline (ARTIP) and results of a sensitive blind search for HI and OH absorbers at $z<0.4$ and $z<0.7$, respectively. ARTIP is written in Python 3.6, extensively uses the Common Astronomy Software Application (CASA) tools and tasks, and is designed to enable the geographically-distributed MeerKAT Absorption Line Survey (MALS) team to collaboratively process large volumes of radio interferometric data. We apply it to the first MALS dataset obtained using the 64-dish MeerKAT radio telescope and 32K channel mode of the correlator. With merely 40 minutes on target, we present the most sensitive spectrum of PKS1830-211 ever obtained and characterize the known HI ($z=0.19$) and OH ($z=0.89$) absorbers. We further demonstrate ARTIP's capabilities to handle realistic observing scenarios by applying it to a sample of 72 bright radio sources observed with the upgraded Giant Metrewave Radio Telescope (uGMRT) to blindly search for HI and OH absorbers. We estimate the numbers of HI and OH absorbers per unit redshift to be $n_{21}(z\sim0.18)<$0.14 and $n_{\rm OH}(z\sim0.40)<$0.12, respectively, and constrain the cold gas covering factor of galaxies at large impact parameters (50 kpc $<\rho<$ 150 kpc) to be less than 0.022. Due to the small redshift path, $\Delta z\sim$13 for HI with column density$>5.4\times10^{19}$ cm$^{-2}$, the survey has probed only the outskirts of star-forming galaxies at $\rho>30$ kpc. MALS with the expected $\Delta z\sim10^{3-4}$ will overcome this limitation and provide stringent constraints on the cold gas fraction of galaxies in diverse environments over $0

CASA, the Common Astronomy Software Applications package, is the primary data processing software for the Atacama Large Millimeter/submillimeter Array (ALMA) and NSF's Karl G. Jansky Very Large Array (VLA), and is frequently used also for other radio telescopes. The CASA software can process data from both single-dish and aperture-synthesis telescopes, and one of its core functionalities is to support the data reduction and imaging pipelines for ALMA, VLA and the VLA Sky Survey (VLASS). CASA has recently undergone several exciting new developments, including an increased flexibility in Python (CASA 6), support of Very Long Baseline Interferometry (VLBI), performance gains through parallel imaging, data visualization with the new Cube Analysis Rendering Tool for Astronomy (CARTA), enhanced reliability and testing, and modernized documentation. These proceedings of the 2019 Astronomical Data Analysis Software & Systems (ADASS) conference give an update of the CASA project, and detail how these new developments will enhance user experience of CASA., Comment: 4 pages, to appear in proceedings of ADASS XXIX, ASP Conf. Series

Frequency dependant complex magnetic permeability is used to understand RF-microwave behaviour of magnetic nanoparticles in the frequency range 250 MHz to 3 GHz. The stable dispersions of Fe3O4 nanoparticles with mean size varying between 11 to 16 nm are prepared for this purpose. The effect of mean particle size and external static magnetic field over microwave absorption properties of magnetic fluid is studied. It is observed that frequency of ferrimagnetic resonance, frequency of maximum absorption, loss tangent and reflection loss (RL) can be controlled by modifying mean particle size and strength of applied external static magnetic field. This kind of study can be useful for radio-microwave devices like tunable attenuator, EM sheilder, and other applications like Hyperthermia., Comment: 13 pages, 10 figures, 2 tables

We present the DSA-2000: a world-leading radio survey telescope and multi-messenger discovery engine for the next decade. The array will be the first true radio camera, outputting science-ready image data over the 0.7 - 2 GHz frequency range with a spatial resolution of 3.5 arcsec. With 2000 x 5 m dishes, the DSA-2000 will have an equivalent point-source sensitivity to SKA1-mid, but with ten times the survey speed. The DSA-2000 is envisaged as an all-sky survey instrument complementary to the ngVLA, and as a counterpart to the LSST (optical), SPHEREx (near-infrared) and SRG/eROSITA (X-ray) all-sky surveys. Over a five-year prime phase, the DSA-2000 will image the entire sky above declination -30 degrees every four months, detecting > 1 unique billion radio sources in a combined full-Stokes sky map with 500 nJy/beam rms noise. This all-sky survey will be complemented by intermediate and deep surveys, as well as spectral and polarization image cubes. The array will be a cornerstone for multi-messenger science, serving as the principal instrument for the US pulsar timing array community, and by searching for radio afterglows of compact object mergers detected by LIGO and Virgo. The array will simultaneously detect and localize ~10,000 fast radio bursts each year, realizing their ultimate use as a cosmological tool. The DSA-2000 will be proposed to the NSF Mid-Scale Research Infrastructure-2 program with a view to first light in 2026, Comment: Submitted as a Project White Paper for the Astro2020 Decadal Survey (12 pages, 7 figures, 2 table)

The Very Large Array Sky Survey (VLASS) is a synoptic, all-sky radio sky survey with a unique combination of high angular resolution ($\approx$2.5"), sensitivity (a 1$\sigma$ goal of 70 $\mu$Jy/beam in the coadded data), full linear Stokes polarimetry, time domain coverage, and wide bandwidth (2-4 GHz). The first observations began in September 2017, and observing for the survey will finish in 2024. VLASS will use approximately 5500 hours of time on the Karl G. Jansky Very Large Array (VLA) to cover the whole sky visible to the VLA (Declination $>-40^{\circ}$), a total of 33,885 deg$^2$. The data will be taken in three epochs to allow the discovery of variable and transient radio sources. The survey is designed to engage radio astronomy experts, multi-wavelength astronomers, and citizen scientists alike. By utilizing an "on the fly" interferometry mode, the observing overheads are much reduced compared to a conventional pointed survey. In this paper, we present the science case and observational strategy for the survey, and also results from early survey observations., Comment: 34 pages, accepted by PASP (modified from prior version to address referee's and coauthor comments). (v2) Minor fixes to author list

We present the results of a time-resolved photometric and time-series analysis of a RR Lyrae type star RR Gem. The main results are as follows: We found RR Gem's pulsation period, 0.39689 d, and its V and I mean magnitudes, 11.277 (V) and 11.063 (I) mag respectively. We confirm its variability type as RRab/BL because of its showing the Bla$\breve{z}$ko effect and it also shows asymmetric light curves (steep ascending branches), periods from 0.3 to 1.0 days, and amplitudes from 0.3 to 2 mag. in V. They are fundamental mode pulsators., Comment: 8 pages, 9 figures

We have carried out VLA imaging and a Fermi timing analysis of the 115 ms gamma-ray and radio pulsar PSR J0002+6216. We found that the pulsar lies at the apex of a narrowly collimated cometary-like 7 arcmin tail of non-thermal radio emission which we identify as a bow-shock pulsar wind nebula. The tail of the nebula points back toward the geometric center of the supernova remnant CTB 1 (G116.9+0.2) 28 arcmin away, at a position angle $\theta_\mu=113^\circ$. We measure a proper motion with 2.9$\sigma$ significance from a Fermi timing analysis giving $\mu$=115$\pm$33 mas yr$^{-1}$ and $\theta_\mu=121^\circ\pm{13}^\circ$, corresponding to a large transverse pulsar velocity of 1100 km s$^{-1}$ at a distance of 2 kpc. This proper motion is of the right magnitude and direction to support the claim that PSR J0002+6216 was born from the same supernova that produced CTB 1. We explore the implications for pulsar birth periods, asymmetric supernova explosions, and mechanisms for pulsar natal kick velocities., Comment: 14 pages, 4 figures, accepted for publication in ApJL

This paper presents the RF energy harvesting module (RECTENNA). The working range of this module includes multiple bands i.e. GSM, ISM, WLAN, and UWB band. To enhance the capturing RF power capability an array arrangement of coplanar monopole antenna has been proposed. Wilkinson power combiner has also been implemented to combine the powers of this antenna array. The RF DC converter circuit having seven stages has also been integrated with this structure. This module produces the DC voltage of 1.8V with respect to +40dB RF input. It is the unique module because it has no need of port connectors. The impedance matching of antenna and converter has been fulfilled by incorporating the passive component at the combiners branch. The value of this passive component is kept equal to the existing value of impedance at input port of converter circuit.

The effects of the antenna far-field power pattern limits the imaging performance of modern wide-bandwidth, high-sensitivity interferometric radio telescopes. Given a model for the aperture illumination pattern (AIP) of the antenna, referred to as the A-term, the wide-band (WB) A-Projection algorithm corrects for the effects of its time, frequency, and polarization structure. The level to which this correction is possible depends on how accurately the A-term, represents the true AIP. In this paper, we describe the A-Solver methodology that combines physical modeling with optimization to holographic measurements to build an accurate model for the AIP. Using a parametrized ray-tracing code as the predictor, we solve for the frequency dependence of the antenna optics and show that the resulting low-order model for the Karl G. Jansky Very Large Array (VLA) antenna captures the dominant frequency-dependent terms. The A-Solver methodology described here is generic and can be adapted for other types of antennas as well. The parameterization is based on the physical characteristics of the antenna structure and optics and is therefore arguably a compact representation (minimized degrees of freedom) of the frequency-dependent structure of the antenna A-term. In this paper, we also show that the parameters derived from A-Solver methodology are expected to improve sensitivity and imaging performance out to the first side-lobe of the antenna., Comment: 13pages, 10 Figures, Accepted for publication at the Astronomical Journal

This paper is concerned with algorithms for calibration of direction dependent effects (DDE) in aperture synthesis radio telescopes (ASRT). After correction of Direction Independent Effects (DIE) using self-calibration, imaging performance can be limited by the imprecise knowledge of the forward gain of the elements in the array. In general, the forward gain pattern is directionally dependent and varies with time due to a number of reasons. Some factors, such as rotation of the primary beam with Parallactic Angle for Azimuth-Elevation mount antennas are known a priori. Some, such as antenna pointing errors and structural deformation/projection effects for aperture-array elements cannot be measured {\em a priori}. Thus, in addition to algorithms to correct for DD effects known a priori, algorithms to solve for DD gains are required for high dynamic range imaging. Here, we discuss a mathematical framework for antenna-based DDE calibration algorithms and show that this framework leads to computationally efficient optimal algorithms which scale well in a parallel computing environment. As an example of an antenna-based DD calibration algorithm, we demonstrate the Pointing SelfCal algorithm to solve for the antenna pointing errors. Our analysis show that the sensitivity of modern ASRT is sufficient to solve for antenna pointing errors and other DD effects. We also discuss the use of the Pointing SelfCal algorithm in real-time calibration systems and extensions for antenna Shape SelfCal algorithm for real-time tracking and corrections for pointing offsets and changes in antenna shape., Comment: 11 pages, 5 figures, accepted for publication in The Astronomical Journal

Deep galaxy surveys have revealed that the global star formation rate (SFR) density in the Universe peaks at 1 < z < 2 and sharply declines towards z = 0. But a clear picture of the underlying processes, in particular the evolution of cold atomic (~100 K) and molecular gas phases, that drive such a strong evolution is yet to emerge. MALS is designed to use MeerKAT's L- and UHF-band receivers to carry out the most sensitive (N(HI)>10$^{19}$ cm$^{-2}$) dust-unbiased search of intervening HI 21-cm and OH 18-cm absorption lines at 0 < z < 2. This will provide reliable measurements of the evolution of cold atomic and molecular gas cross-sections of galaxies, and unravel the processes driving the steep evolution in the SFR density. The large sample of HI and OH absorbers obtained from the survey will (i) lead to tightest constraints on the fundamental constants of physics, and (ii) be ideally suited to probe the evolution of magnetic fields in disks of galaxies via Zeeman Splitting or Rotation Measure synthesis. The survey will also provide an unbiased census of HI and OH absorbers, i.e. cold gas associated with powerful AGNs (>10$^{24}$ W Hz$^{-1}$) at 0 < z < 2, and will simultaneously deliver a blind HI and OH emission line survey, and radio continuum survey. Here, we describe the MALS survey design, observing plan and the science issues to be addressed under various science themes., Comment: 16 pages, 3 figures Accepted for publication, Proceedings of Science, Workshop on "MeerKAT Science: On the Pathway to the SKA", held in Stellenbosch 25-27 May, 2016

Next generation radio telescope arrays are being designed and commissioned to accurately measure polarized intensity and rotation measures across the entire sky through deep, wide-field radio interferometric surveys. Radio interferometer dish antenna arrays are affected by direction-dependent (DD) gains due to both instrumental and atmospheric effects. In this paper we demonstrate the effect of DD errors for parabolic dish antenna array on the measured polarized intensities of radio sources in interferometric images. We characterize the extent of polarimetric image degradation due to the DD gains through wide-band VLA simulations of representative point source simulations of the radio sky at L-Band(1-2GHz). We show that at the 0.5 gain level of the primary beam (PB) there is significant flux leakage from Stokes $I$ to $Q$, $U$ amounting to 10\% of the total intensity. We further demonstrate that while the instrumental response averages down for observations over large parallactic angle intervals, full-polarization DD correction is required to remove the effects of DD leakage. We also explore the effect of the DD beam on the Rotation Measure(RM) signals and show that while the instrumental effect is primarily centered around 0 rad-m$^{-2}$, the effect is significant over a broad range of RM requiring full polarization DD correction to accurately reconstruct RM synthesis signal., Comment: version2, 13 pages, 12 figures, accepted for publication in The Astronomical Journal

We present $\simeq$0$.\!\!^{\prime\prime}4$-resolution extinction-independent distributions of star formation and dust in 11 star-forming galaxies (SFGs) at $z = 1.3-3.0$. These galaxies are selected from sensitive, blank-field surveys of the $2' \times 2'$ Hubble Ultra-Deep Field at $\lambda = 5$ cm and 1.3 mm using the Karl G. Jansky Very Large Array (VLA) and Atacama Large Millimeter/submillimeter Array (ALMA). They have star-formation rates (SFRs), stellar masses, and dust properties representative of massive main-sequence SFGs at $z \sim 2$. Morphological classification performed on spatially-resolved stellar mass maps indicates a mixture of disk and morphologically disturbed systems; half of the sample harbor X-ray active galactic nuclei (AGN), thereby representing a diversity of $z \sim 2$ SFGs undergoing vigorous mass assembly. We find that their intense star formation most frequently occurs at the location of stellar-mass concentration and extends over an area comparable to their stellar-mass distribution, with a median diameter of $4.2 \pm 1.8$ kpc. This provides direct evidence for galaxy-wide star formation in distant, blank-field-selected main-sequence SFGs. The typical galactic-average SFR surface density is 2.5 M$_{\odot}$yr$^{-1}$kpc$^{-2}$, sufficiently high to drive outflows. In X-ray-selected AGN where radio emission is enhanced over the level associated with star formation, the radio excess pinpoints the AGN, which are found to be co-spatial with star formation. The median extinction-independent size of main-sequence SFGs is two times larger than those of bright submillimeter galaxies whose SFRs are $3-8$ times larger, providing a constraint on the characteristic SFR ($\sim300$ M$_{\odot}$yr$^{-1}$) above which a significant population of more compact star-forming galaxies appears to emerge., Comment: 11 pages, 4 figures, accepted for publication in ApJ

Context. Most popular algorithms in use to remove the effects of a telescope's point spread function (PSF) in radio astronomy are variants of the CLEAN algorithm. Most of these algorithms model the sky brightness using the delta-function basis, which results in undesired artefacts when used on image extended emission. The adaptive scale pixel decomposition (Asp-Clean) algorithm models the sky brightness on a scale-sensitive basis and thus gives a significantly better imaging performance when imaging fields that contain both resolved and unresolved emission. Aims. However, the runtime cost of Asp-Clean is higher than that of scale-insensitive algorithms. In this paper, we identify the most expensive step in the original Asp-Clean algorithm and present an efficient implementation of it, which significantly reduces the computational cost while keeping the imaging performance comparable to the original algorithm. The PSF sidelobe levels of modern wide-band telescopes are significantly reduced, allowing us to make approximations to reduce the computing cost, which in turn allows for the deconvolution of larger images on reasonable timescales. Methods. As in the original algorithm, scales in the image are estimated through function fitting. Here we introduce an analytical method to model extended emission, and a modified method for estimating the initial values used for the fitting procedure, which ultimately leads to a lower computational cost. Results.The new implementation was tested with simulated EVLA data and the imaging performance compared well with the original Asp-Clean algorithm. Tests show that the current algorithm can recover features at different scales with lower computational cost., Comment: 6 pages; 4 figures

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy., Comment: 11 pages, 7 figures, 2 tables; accepted for publication in the Astrophysical Journal Letters; this version with small changes to affiliations

We present LOFAR Low Band observations of the Bootes and 3C295 fields. Our images made at 34, 46, and 62 MHz reach noise levels of 12, 8, and 5 mJy beam$^{-1}$, making them the deepest images ever obtained in this frequency range. In total, we detect between 300 and 400 sources in each of these images, covering an area of 17 to 52 deg$^{2}$. From the observations we derive Euclidean-normalized differential source counts. The 62 MHz source counts agree with previous GMRT 153 MHz and VLA 74 MHz differential source counts, scaling with a spectral index of $-0.7$. We find that a spectral index scaling of $-0.5$ is required to match up the LOFAR 34 MHz source counts. This result is also in agreement with source counts from the 38 MHz 8C survey, indicating that the average spectral index of radio sources flattens towards lower frequencies. We also find evidence for spectral flattening using the individual flux measurements of sources between 34 and 1400 MHz and by calculating the spectral index averaged over the source population. To select ultra-steep spectrum ($\alpha < -1.1$) radio sources, that could be associated with massive high redshift radio galaxies, we compute spectral indices between 62 MHz, 153 MHz and 1.4 GHz for sources in the Bo\"otes field. We cross-correlate these radio sources with optical and infrared catalogues and fit the spectral energy distribution to obtain photometric redshifts. We find that most of these ultra-steep spectrum sources are located in the $ 0.7 \lesssim z \lesssim 2.5$ range., Comment: 26 pages, 20 figures, ApJ in press

The new broad-band capabilities of large radio interferometers such as the GMRT and JVLA allow for long-integration mosaic imaging observations to create ultra-deep full-polarization images of the sky over wide frequency ranges. Achieving rms sensitivities of order 1 $\mu$Jy, these observations explore the radio source population at flux densities well below the regime dominated by classical radio galaxies and Active Galactic Nuclei. We present initial results from radio sources revealed with deep mosaicking observations with the GMRT and JVLA at respectively 0.6 and 5 GHz, and evidence that the $\mu$Jy sensitivity level marks the transition to detection of polarized emission from a population of sources dominated by emission from magnetic fields in the disks of starburst and normal galaxies.

The Sun has remained a difficult source to image for radio telescopes, especially at the low radio frequencies. Its morphologically complex emission features span a large range of angular scales, emission mechanisms involved and brightness temperatures. In addition, time and frequency synthesis, the key tool used by most radio interferometers to build up information about the source being imaged is not effective for solar imaging, because many of the features of interest are short lived and change dramatically over small fractional bandwidths. Building on the advances in radio frequency technology, digital signal processing and computing, the kind of instruments needed to simultaneously capture the evolution of solar emission in time, frequency, morphology and polarization over a large spectral span with the requisite imaging fidelity, and time and frequency resolution have only recently begun to appear. Of this class of instruments, the Murchison Widefield Array (MWA) is best suited for solar observations. The MWA has now entered a routine observing phase and here we present some early examples from MWA observations., Comment: 4 pages, 3 figures, submitted for the 31st URSI General Assembly and Scientific Symposium, to be held in Bejing, China from 16-23, August, 2014

Variations of the antenna primary beam (PB) pattern as a function of time, frequency and polarization form one of the dominant direction-dependent effects at most radio frequency bands. These gains may also vary from antenna to antenna. The A-Projection algorithm, published earlier, accounts for the effects of the narrow-band antenna PB in full polarization. In this paper we present the Wide-Band A-Projection algorithm (WB A-Projection) to include the effects of wide bandwidth in the A-term itself and show that the resulting algorithm simultaneously corrects for the time, frequency and polarization dependence of the PB. We discuss the combination of the WB A-Projection and the Multi-term Multi Frequency Synthesis (MT-MFS) algorithm for simultaneous mapping of the sky brightness distribution and the spectral index distribution across a wide field of view. We also discuss the use of the narrow-band A-Projection algorithm in hybrid imaging schemes that account for the frequency dependence of the PB in the image domain., Comment: 10 pages, 5 figures; accepted for publication in ApJ

The aimed high sensitivities and large fields of view of the new generation of interferometers impose to reach high dynamic range of order $\sim$1:$10^6$ to 1:$10^8$ in the case of the Square Kilometer Array. The main problem is the calibration and correction of the Direction Dependent Effects (DDE) that can affect the electro-magnetic field (antenna beams, ionosphere, Faraday rotation, etc.). As shown earlier the A-Projection is a fast and accurate algorithm that can potentially correct for any given DDE in the imaging step. With its very wide field of view, low operating frequency ($\sim30-250$ MHz), long baselines, and complex station-dependent beam patterns, the Low Frequency Array (LOFAR) is certainly the most complex SKA precursor. In this paper we present a few implementations of A-Projection applied to LOFAR that can deal with non-unitary station beams and non-diagonal Mueller matrices. The algorithm is designed to correct for all the DDE, including individual antenna, projection of the dipoles on the sky, beam forming and ionospheric effects. We describe a few important algorithmic optimizations related to LOFAR's architecture allowing us to build a fast imager. Based on simulated datasets we show that A-Projection can give dramatic dynamic range improvement for both phased array beams and ionospheric effects. We will use this algorithm for the construction of the deepest extragalactic surveys, comprising hundreds of days of integration.

The radio continuum emission from the Galaxy has a rich mix of thermal and non-thermal emission. This very richness makes their interpretation challenging since the low radio opacity means that a radio image represents the sum of all emission regions along the line-of-sight. These challenges make the existing narrow-band radio surveys of the Galactic plane difficult to interpret: e.g. a small region of emission might be a supernova remnant (SNR) or an HII region, or a complex combination of both. Instantaneous wide bandwidth radio observations in combination with the capability for high resolution spectral index mapping, can be directly used to disentangle these effects. Here we demonstrate simultaneous continuum and spectral index imaging capability at the full continuum sensitivity and resolution using newly developed wide-band wide-field imaging algorithms. Observations were done in the L- and C-Band with a total bandwidth of 1 and 2 GHz respectively. We present preliminary results in the form of a full-field continuum image covering the wide-band sensitivity pattern of the EVLA centered on a large but poorly studied SNR (G55.7+3.4) and relatively narrower field continuum and spectral index maps of three fields containing SNR and diffused thermal emission. We demonstrate that spatially resolved spectral index maps differentiates regions with emission of different physical origin (spectral index variation across composite SNRs and separation of thermal and non-thermal emission), superimposed along the line of sight. The wide-field image centered on the SNR G55.7+3.4 also demonstrates the excellent wide-field wide-band imaging capability of the EVLA., Comment: 6 pages; 4 figures; accepted for publication in ApJL, to appear in the ApJL EVLA special issue

Observations of the HI 21cm transition line promises to be an important probe into the cosmic dark ages and epoch of reionization. One of the challenges for the detection of this signal is the accuracy of the foreground source removal. This paper investigates the extragalactic point source contamination and how accurately the bright sources ($\gtrsim 1$ ~Jy) should be removed in order to reach the desired RMS noise and be able to detect the 21cm transition line. Here, we consider position and flux errors in the global sky-model for these bright sources as well as the frequency independent residual calibration errors. The synthesized beam is the only frequency dependent term included here. This work determines the level of accuracy for the calibration and source removal schemes and puts forward constraints for the design of the cosmic reionization data reduction scheme for the upcoming low frequency arrays like MWA,PAPER, etc. We show that in order to detect the reionization signal the bright sources need to be removed from the data-sets with a positional accuracy of $\sim 0.1$ arc-second. Our results also demonstrate that the efficient foreground source removal strategies can only tolerate a frequency independent antenna based mean residual calibration error of $\lesssim 0.2 %$ in amplitude or $\lesssim 0.2$ degree in phase, if they are constant over each days of observations (6 hours). In future papers we will extend this analysis to the power spectral domain and also include the frequency dependent calibration errors and direction dependent errors (ionosphere, primary beam, etc)., Comment: accepted by ApJ; 12 pages, 10 figures

Many scientific deliverables of the next generation low frequency radio telescopes require high dynamic range imaging. Next generation telescopes under construction indeed promise at least a ten-fold increase in the sensitivity compared with existing telescopes. The projected achievable RMS noise in the images from these telescopes is in the range of 1--10$\mu$Jy/beam corresponding to typical imaging dynamic ranges of $10^{6-7}$. High imaging dynamic range require removal of systematic errors to high accuracy and for long integration intervals. In general, many source of errors are directionally dependent and unless corrected for, will be a limiting factor for the imaging dynamic range of these next generation telescopes. This requires development of new algorithms and software for calibration and imaging which can correct for such direction and time dependent errors. In this paper, I discuss the resulting algorithmic and computing challenges and the recent progress made towards addressing these challenges., Comment: to appear in: The Low-Frequency Radio Universe ASP Conference Series, Vol. LFRU, 2009, Eds: D.J. Saikia, Dave Green, Y. Gupta and Tiziana Venturi

The nature of X-shaped sources is a matter of considerable debate: it has even been proposed that they provide evidence for black hole mergers$ / $spin reorientation, and therefore constrain the rate of strong gravitational wave events (Merritt & Ekers 2002). Based on morphological and spectral characteristics of these sources, currently a strong contender to explain the nature of these sources is the `alternative' model of Lal & Rao (2007), in which these sources consist of two pairs of jets, which are associated with two unresolved AGNs. Detailed morphological and spectral results on milliarcsecond-scales (mas) provide a crucial test of this model, and hence these sources are excellent candidates to study on mas; {\it i.e.}, to detect he presence/absence of double nuclei/AGNs, signs of helical/disrupted jets, thereby, to investigate spatially resolved/unresolved binary AGN systems and providing clues to understanding the physics of merging of AGNs on mas. We conducted a systematic study of a large sample of known X-shaped, comparison FR II radio galaxies, and newly discovered X-shaped candidate sources using Giant Metrewave Radio Telescope and Very Large Array at several radio frequencies. In our new observations of `comparison' FR II radio galaxies we find that almost all of our targets show standard spectral steepening as a function of distance from the hotspot. However, one source, 3C 321, has a low-surface-brightness extension that shows a flatter spectral index than the high-surface-brightness hotspots$ / $lobes, as found in `known' X-shaped sources., Comment: To appear in the proceedings of the Conference "The Central Kiloparsec: Active Galactic Nuclei and Their Hosts", Ierapetra 4-6 June 2008, Greece, Memorie della Societa Astronomica Italiana

This paper summarizes some of the major calibration and image reconstruction techniques used in radio interferometry and describes them in a common mathematical framework. The use of this framework has a number of benefits, ranging from clarification of the fundamentals, use of standard numerical optimization techniques, and generalization or specialization to new algorithms.

We consider a troublesome form of non-isoplanatism in synthesis radio telescopes: non-coplanar baselines. We present a novel interpretation of the non-coplanar baselines effect as being due to differential Fresnel diffraction in the neighborhood of the array antennas. We have developed a new algorithm to deal with this effect. Our new algorithm, which we call "W-projection", has markedly superior performance compared to existing algorithms. At roughly equivalent levels of accuracy, W-projection can be up to an order of magnitude faster than the corresponding facet-based algorithms. Furthermore, the precision of result is not tightly coupled to computing time. W-projection has important consequences for the design and operation of the new generation of radio telescopes operating at centimeter and longer wavelengths., Comment: Accepted for publication in "IEEE Journal of Selected Topics in Signal Processing"

Astronomical imaging using aperture synthesis telescopes requires deconvolution of the point spread function as well as calibration of instrumental and atmospheric effects. In general, such effects are time-variable and vary across the field of view as well, resulting in direction-dependent (DD), time-varying gains. Most existing imaging and calibration algorithms assume that the corruptions are direction independent, preventing even moderate dynamic range full-beam, full-Stokes imaging. We present a general framework for imaging algorithms which incorporate DD errors. We describe as well an iterative deconvolution algorithm that corrects known DD errors due to the antenna power patterns and pointing errors for high dynamic range full-beam polarimetric imaging. Using simulations we demonstrate that errors due to realistic primary beams as well as antenna pointing errors will limit the dynamic range of upcoming higher sensitivity instruments and that our new algorithm can be used to correct for such errors. We have applied this algorithm to VLA 1.4 GHz observations of a field that contains two ``4C'' sources and have obtained Stokes-I and -V images with systematic errors that are one order of magnitude lower than those obtained with conventional imaging tools. Our simulations show that on data with no other calibration errors, the algorithm corrects pointing errors as well as errors due to known asymmetries in the antenna pattern., Comment: submitted to A&A; some clarifications added in the text; accepted for publication in A&A

We propose an extension of a result by Repetowicz et al. about Wick's theorem and its applications: we first show that Wick's theorem can be extended to the uniform distribution on the sphere and then to the whole class of elliptical distributions. Then, as a special case, we detail this result for distributions that are Gaussian scale mixtures. Finally, we show that these results allow to recover easily a theorem by Folland about integration of polynomials over the sphere.

Radio and near-infared observations towards the steep spectrum Galactic plane radio source WKB 0314+57.8 are presented, in order to clarify the nature of this source. The radio observations include archival and survey data, together with new Giant Metrewave Radio Telescope observations at 617 MHz. The near-infrared observations are in the J and K bands, from the Gemini instrument on the Shane 3-m telescope. The radio observations show that WKB 0314+57.8 is extended, with an very steep spectrum (with flux density proportional to frequency to -2.5 power between 40 MHz and 1.5 GHz). The colour--magnitude diagram constructed from near-infrared observations of the field suggests the presence of a z approx 0.08 galaxy cluster behind the Galactic plane, reddened by about 6 magnitudes of visual extinction. Although the steep spectrum source has no obvious identification, two other radio sources in the field covered by the near-infrared observations have tentative identifications with galaxies. These observations indicate that WKB 0314+57.8 is a relic source in a cluster of galaxies, not a pulsar., Comment: 6 pages, to appear in MNRAS, typos corrected

Deconvolution of the telescope Point Spread Function (PSF) is necessary for even moderate dynamic range imaging with interferometric telescopes. The process of deconvolution can be treated as a search for a model image such that the residual image is consistent with the noise model. For any search algorithm, a parameterized function representing the model such that it fundamentally separates signal from noise will give optimal results. In this paper, the first in a series of forthcoming papers, we argue that in general, spatial correlation length (a measure of the scale of emission) is a stronger separator of the signal from the noise, compared to the strength of the signal alone. Consequently scale sensitive deconvolution algorithms result into more noise-like residuals. We present a scale-sensitive deconvolution algorithm for radio interferometric images, which models the image as a collection of Adaptive Scale Pixels (Asp). Some attempts at optimizing the runtime performance are also presented., Comment: 9 pages, 5 figures, accepted for publication in Astronomy & Astrophysics

We present the investigations of SN 1993J using low frequency observations with the Giant Meterwave Radio Telescope. We analyze the light curves of SN 1993J at 1420, 610, 325 and 243 MHz during $7.5-10$ years since explosion.The supernova has become optically thin early on in the 1420 MHz and 610 MHz bands while it has only recently entered the optically thin phase in the 325 MHz band. The radio light curve in the 235 MHz band is more or less flat. This indicates that the supernova is undergoing a transition from an optically thick to optically thin limit in this frequency band. In addition, we analyze the SN radio spectra at five epochs on day 3000, 3200, 3266, 3460 and 3730 since explosion. Day 3200 spectrum shows a synchrotron cooling break. SN 1993J is the only young supernova for which the magnetic field and the size of the radio emitting region are determined through unrelated methods. Thus the mechanism that controls the evolution of the radio spectra can be identified. We suggest that at all epochs, the synchrotron self absorption mechanism is primarily responsible for the turn-over in the spectra. Light curve models based on free free absorption in homogeneous or inhomogeneous media at high frequencies overpredict the flux densities at low frequencies. The discrepancy is increasingly larger at lower and lower frequencies. We suggest that an extra opacity, sensitively dependent on frequency, is likely to account for the difference at lower frequencies. The evolution of the magnetic field (determined from synchrotron self absorption turn-over) is roughly consistent with $B \propto t^{-1}$. Radio spectral index in the optically thin part evolves from $\alpha \sim 0.8-1.0$ at few tens of days to $\sim 0.6$ in about 10 years., Comment: 37 pages, 9 figures in LaTex; scheduled for ApJ 10 September 2004, v612 issue; send comments to: poonam@tifr.res.in

In this paper, we discuss the low frequency spectrum of SN 1993J with GMRT. We observed SN 1993J at several epochs in 20cm, 50cm, 90cm and 125cm wavelengths and achieved near simultaneous spectra. We fit synchrotron self absorption (SSA) and free-free models to the data. We compare the size of SN obtained using SSA fits to that of size extrapolated from VLBI measurements at various epochs using public data at earlier epochs. We find that the synchrotron self absorption process is insufficient to reproduce the observed size of the supernova under the assumption of equipartition between magnetic fields and relativistic electrons. We also derive the evolution of spectral index and magnetic field at several epochs., Comment: 5 pages, 12 figures. Uses svmult.cls. To appear in proceedings of IAU Colloquium 192 "Supernovae (10 years of SN 1993J)", April 2003, Valencia, Spain, eds. J. M. Marcaide and K. W. Weiler

Radio and X-ray studies of young supernovae probe the interaction between the supernova shock waves and the surrounding medium and give clues to the nature and past of the progenitor star. Here we discuss the early emission from type Ic SN 2002ap and argue that repeated Compton boosting of optical photons by hot electrons presents the most natural explanation of the prompt X-ray emission. We describe the radio spectrum of another type Ic SN 2003dh (GRB030329) obtained with combined GMRT and VLA data. We report on the low frequency radio monitoring of SN 1995N and our objectives of distinguishing between competing models of X-ray emission from this SN and the nature of its progenitor by X-ray spectroscopy. Radio studies on SN 2001gd, SN 2001ig and SN 2002hh are mentioned., Comment: 5 pages, 4 figures. Uses svmult.cls. To appear in proceedings of IAU Colloquium 192 "Supernovae (10 years of SN 1993J)", April 2003, Valencia, Spain, eds. J. M. Marcaide and K. W. Weiler

We report on the combined X-ray and radio observations of the type Ic SN 2002ap, using XMM-Newton ToO observation of M74 and the Giant Metrewave Radio Telescope (GMRT). We account for the presence of a nearby source in the pre-supernova Chandra field of view in our measurements of the X-ray flux (0.3 - 10 KeV) 5.2 days after the explosion. The X-ray spectrum is well fitted by a power law spectrum with photon index $\alpha= 2.6$. Our results suggest that the prompt X-ray emission originates from inverse Compton scattering of photospheric thermal emission by energetic electrons. Radio observations with the GMRT at 610 MHz (8 days after the explosion) and 1420 MHz (70 days after the explosion) are combined with the high frequency VLA observations of SN 2002ap reported earlier, and the early radiospheric properties of SN 2002ap are compared with similar data from two other supernovae. Finally, the GMRT radio map reveals four other X-ray sources in the field of view of M74 with radio counterparts., Comment: Replaced with the version accepted for publication in Astron.& Astrophys. (main journal). The text has been revised in places, but our conclusions remain unchanged