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<z<1.5$.

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.

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.

We present the first results of a time-resolved CCD photometric and time-series analysis of NSV 13601, a variable star of Pegasus constellation. The 14" SCT of Maharaja Krishnakumarsinhji Bhavnagar University(MKBU) was used for observations. Analysis was performed both in Vstar and Period04 software and compared their results. The main results are as follows: From data analysis, it's period, $\sim$77.784 d(MKBU Data V)(VStar), $\sim$77.632 d(MKBU Data Ic)(VStar), $\sim$77.058 d(MKBU Data V)(Period04) and $\sim$49.560 d(MKBU Data Ic)(Period04). It's V and I mean magnitudes for MKBU Data are 12.203 (V) and 11.292 (I) mag respectively. We confirm it as a variable star. \keywords{Stars: individual: NSV 13601 --- stars: oscillations --- techniques: photometric.

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

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.

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.

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.

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.

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.

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.

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.

Many aspects of the Atacama Large Millimeter Array (ALMA) instrument are still unknown due to its young age. One such aspect is the true nature of the primary beam of each baseline, and how changes to the individual primary beams affect astronomical observations when said changes are ignored during imaging. This paper aims to create a more thorough understanding of the strengths and weaknesses of ALMA through realistic modeling of the primary beams and simulated observations, which in turn can inform the user of the necessity of implementing more computationally costly algorithms, such as A-Projection, and when simpler, quicker algorithms will suffice. We quantify our results by examining the dynamic range of each observation, along with the ability to reconstruct the Stokes I amplitude of the test sources. These tests conclude that for dynamic ranges of less than 1000, for point sources and sources much smaller than the main lobe of the primary beam, the accuracy of the primary beam model beyond the physical size of the aperture simply doesn't matter. In observations of large extended sources, deconvolution errors dominate the reconstructed images and the individual primary beam errors were indistinguishable from each other.

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.

Many deep wide-band wide-field radio interferometric surveys are being designed to accurately measure intensities, spectral indices and polarization properties of faint source populations. In this paper we compare various wideband imaging methods to evaluate the accuracy to which intensities and spectral indices of sources close to the confusion limit can be reconstructed. We simulated a wideband single-pointing (C-array, L-Band (1-2GHz)) and 46-pointing mosaic(D-array, C-Band (4-8GHz)) JVLA observation using realistic brightness distribution ranging from $1\mu$Jy to $100m$Jy and time-,frequency-, polarization- and direction-dependent instrumental effects. The main results from these comparisons are (a) errors in the reconstructed intensities and spectral indices are larger for weaker sources even in the absence of simulated noise, (b) errors are systematically lower for joint reconstruction methods (such as MT-MFS) along with A-Projection for accurate primary beam correction, and (c) use of MT-MFS for image reconstruction eliminates Clean-bias (which is present otherwise). Auxiliary tests include solutions for deficiencies of data partitioning methods (e.g. the use of masks to remove clean bias and hybrid methods to remove sidelobes from sources left undeconvolved), the effect of sources not at pixel centers and the consequences of various other numerical approximations within software implementations. This paper also demonstrates the level of detail at which such simulations must be done in order to reflect reality, enable one to systematically identify specific reasons for every trend that is observed and to estimate scientifically defensible imaging performance metrics and the associated computational complexity of the algorithms/analysis procedures.

We present a deep centimeter-wavelength catalog of the Orion Nebula Cluster (ONC), based on a 30h single-pointing observation with the Karl G. Jansky Very Large Array in its high-resolution A-configuration using two 1 GHz bands centered at 4.7 GHz and 7.3 GHz. A total of 556 compact sources were detected in a map with a nominal rms noise of 3 muJy/bm, limited by complex source structure and the primary beam response. Compared to previous catalogs, our detections increase the sample of known compact radio sources in the ONC by more than a factor of seven. The new data show complex emission on a wide range of spatial scales. Following a preliminary correction for the wideband primary-beam response, we determine radio spectral indices for 170 sources whose index uncertainties are less than +/-0.5. We compare the radio to the X-ray and near-infrared point-source populations, noting similarities and differences.

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.

Synthesis imaging in radio astronomy is affected by instrumental and atmospheric effects which introduce direction-dependent (DD) gains.The antenna power pattern varies both as a function of time and frequency. The broad band time varying nature of the antenna power pattern when not corrected leads to gross errors in full Stokes imaging and flux estimation. In this poster we explore the errors that arise in image deconvolution while not accounting for the time and frequency dependence of the antenna power pattern. Simulations were conducted with the wide-band full Stokes power pattern of the Karl G. Jansky Very Large Array (VLA) antennas to demonstrate the level of errors arising from direction-dependent gains and their non-neglegible impact on upcoming sky surveys such as the VLASS. DD corrections through hybrid projection algorithms are computationally expensive to perform. A highly parallel implementation through high performance computing architectures is the only feasible way of applying these corrections to the large data sizes of these upcoming surveys.

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.

We report new observations of Abell 2256 with the Karl G. Jansky Very Large Array (VLA) at frequencies between 1 and 8 GHz. These observations take advantage of the 2:1 bandwidths available for a single observation to study the spectral index, polarization and Rotation Measure as well as using the associated higher sensitivity to image total intensity features down to ~0.5" resolution. We find the Large Relic, which dominates the cluster, is made up of a complex of filaments which show correlated distributions in intensity, spectral index, and fractional polarization. The Rotation Measure varies across the face of the Large Relic but is not well correlated with the other properties of the source. The shape of individual filaments suggests that the Large Relic is at least 25 kpc thick. We detect a low surface brightness arc connecting the Large Relic to the Halo and other radio structures suggesting a physical connection between these features. The center of the F-complex is dominated by a very steep-spectrum, polarized, ring-like structure, F2, without an obvious optical identification, but the entire F-complex has interesting morphological similarities to the radio structure of NGC1265. Source C, the Long Tail, is unresolved in width near the galaxy core and is </~100pc in diameter there. This morphology suggests either that C is a one-sided jet or that the bending of the tails takes place very near the core, consistent with the parent galaxy having undergone extreme stripping. Overall it seems that many of the unusual phenomena can be understood in the context of Abell 2256 being near the pericenter of a slightly off-axis merger between a cluster and a smaller group. Given the lack of evidence for a strong shock associated with the Large Relic, other models should be considered, such as reconnection between two large-scale magnetic domains.

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.

A large number of deep and wide-field radio interferometric surveys are being designed to measure accurate statistics of faint source populations. Most require mosaic observations, and expect to benefit from the sensitivity provided by broad-band instruments. In this paper, we present preliminary results from a comparison of several wideband imaging methods in the context of how accurately they reconstruct the intensities and spectral indices of micro-Jy level sources.

(Abridged)One of the fundamental challenges for astrophysics in the 21st century is finding a way to untangle the physical processes that govern galaxy formation and evolution. Given the importance and scope of this problem, the multi-wavelength astronomical community has used the past decade to build up a wealth of information over specific extragalactic deep fields to address key questions in galaxy formation and evolution. These fields generally cover at least 10square degrees to facilitate the investigation of the rarest, typically most massive, galaxies and AGN. Furthermore, such areal coverage allows the environments to be fully accounted for, thereby linking the single halo to the two-halo terms in the halo occupation distribution. Surveys at radio wavelengths have begun to lag behind those at other wavelengths, especially in this medium-deep survey tier. However, the survey speed offered by the JVLA means that we can now reach a point where we can begin to obtain commensurate data at radio wavelengths to those which already exists from the X-ray through to the far-infrared over ~10 square degrees. We therefore present the case for a 10 square degree survey to 1.5uJy at L-band in A or B Array, requiring ~4000 hours to provide census of star-formation and AGN-accretion activity in the Universe. For example, the observations will allow galaxies forming stars at 10Msolar/yr to be detected out to z~1 and luminous infrared galaxies (1000Msolar/yr to be found out to z~6. Furthermore, the survey area ensures that we will have enough cosmic volume to find these rare sources at all epochs. The bandwidth will allow us to determine the polarisation properties galaxies in the high-redshift Universe as a function of stellar mass, morphology and redshift.

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.

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.

In the lead-up to the Square Kilometre Array (SKA) project, several next-generation radio telescopes and upgrades are already being built around the world. These include APERTIF (The Netherlands), ASKAP (Australia), eMERLIN (UK), VLA (USA), e-EVN (based in Europe), LOFAR (The Netherlands), Meerkat (South Africa), and the Murchison Widefield Array (MWA). Each of these new instruments has different strengths, and coordination of surveys between them can help maximise the science from each of them. A radio continuum survey is being planned on each of them with the primary science objective of understanding the formation and evolution of galaxies over cosmic time, and the cosmological parameters and large-scale structures which drive it. In pursuit of this objective, the different teams are developing a variety of new techniques, and refining existing ones. Here we describe these projects, their science goals, and the technical challenges which are being addressed to maximise the science return.

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.

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.

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).

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.

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.

Context. The nature of X-shaped sources is a matter of considerable debate in the literature: 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. Aim. To explore the nature of these X-shaped radio galaxies. Method. We conduct a systematic study of a large sample of known and newly discovered X-shaped sources along with a comparison sample. We used the Giant Metrewave Radio Telescope with resolution of $\sim6^{\prime\prime}$ to $\sim15^{\prime\prime}$ at 610 MHz and 240 MHz in the dual-frequency mode. Preliminary Result. Based on our careful analysis and estimation of the possible systematic errors, the known X-shaped sources divide into the following three categories: (i) the wings have flatter spectral indices than the active lobes, (ii) the wings and the active lobes have comparable spectral indices, and (iii) the wings have steeper spectral indices than the active lobes. In addition, based on our preliminary analysis, one sample source from our comparison sample shows a spectral index result belonging to category (i). Future. Milliarcsecond scale imaging will be conducted on some of these X-shaped sources to investigate if they are examples of binary AGN systems and thereby understand the nature of these sources.

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.

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.

GCRT J1745-3009 is a transient bursting radio source located in the direction of the Galactic center. It was discovered in a 330 MHz VLA observation from 2002 September 30--October 1 and subsequently rediscovered in a 330 MHz GMRT observation from 2003 September 28 by Hyman et al. Here we report a new radio detection of the source in 330 MHz GMRT data taken on 2004 March 20. The observed properties of the single burst detected differ significantly from those measured previously, suggesting that GCRT J1745-3009 was detected in a new physical state. The 2004 flux density was ~0.05 Jy, ~10x weaker than the single 2003 burst and ~30x weaker than the five bursts detected in 2002. We derive a very steep spectral index, alpha = -13.5 +/- 3.0, across the bandpass, a new result previously not detectable due to limitations in the analysis of the 2002 and 2003 observations. Also, the burst was detected for only ~2 min., in contrast to the 10 min. duration observed in the earlier bursts. Due to sparse sampling, only the single burst was detected in 2004, as in the 2003 epoch, and we cannot rule out additional undetected bursts that may have occurred with the same ~77 min. periodicity observed in 2002 or with a different periodicity. Considering our total time on source throughout both our archival and active monitoring campaigns, we estimate the source exhibits detectable bursting activity ~7% of the time.

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.

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.

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.

We combine the GMRT low frequency radio observations of SN 1993J with the VLA high frequency radio data to get a near simultaneous spectrum around day 3200 since explosion. The low frequency measurements of the supernova determine the turnover frequency and flux scale of the composite spectrum and help reveal a steepening in the spectral index, $\Delta \alpha \sim 0.6$, in the optically thin part of the spectrum. This is the first observational evidence of a break in the radio spectrum of a young supernova. We associate this break with the phenomenon of synchrotron aging of radiating electrons. From the break in the spectrum we calculate the magnetic field in the shocked region independent of the equipartition assumption between energy density of relativistic particles and magnetic energy density. We determine the ratio of these two energy densities and find that this ratio is in the range: $8\times 10^{-6}-5\times 10^{-4}$. We also predict the nature of the evolution of the synchrotron break frequency with time, with competing effects due to diffusive Fermi acceleration and adiabatic expansion of the radiative electron plasma.

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.

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.

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.

Multi-frequency observations of two radio bright SNe (SN1993J, SN1979C) carried out over a year with GMRT are presented. Their radio light curves trace the evolution of the mass loss in the stellar wind before the pre-supernova star exploded. Their spectra at low frequencies distinguish between the synchrotron self-absorption and free-free absorption based models. Two very similar type IIp SN1999gi and SN1999em were also observed and we comment on their progenitor mass. These observations of different subclasses of core-collapse supernovae will aid in defining their progenitors.

Results from the 327-MHz Giant Meterwave Radio Telescope (GMRT) observations of four fields - three fields containing Galactic supernova remnants (SNRs) (G001.4-0.1, G003.8+0.3, and G356.3-1.5) and one field containing a candidate SNR (G004.2+0.0), are reported in this paper. These fields were selected from the 843-MHz survey conducted by Gray (1994b) using the Molonglo Synthesis Telescope (MOST). All the three SNRs are detected in the GMRT images. Significant amount of thermal emission is seen at the location of the candidate SNR, and the GMRT image shows a discrete source of emission which is consistent with it being a flat spectrum thermal source. An incomplete arc of emission with a compact central source detected in one of the fields, is coincident with the extended OH(1720 MHz) emission reported earlier. Possible implications of this morphology and correlation with the OH(1720 MHz) emission are also discussed.

Mechanical and electronic imperfections can result into polarization leakage in individual antennas of a radio interferometer. Such leakages manifest themselves as closure errors even in co-polar visibility measurements of unpolarized sources. This paper describes and tests a method for the study of polarization leakage for radio interferometric telescopes using {\it only} the nominally co-polar visibilities for unpolarized calibrators. Interpretation of the resulting closure phases on the Poincar\'e sphere is presented. Since unpolarized sources are used, the actual solutions for leakage parameters is subject to a degeneracy which is discussed. This however, does not affect the correction of closure errors in our scheme.