Showing 1–50 of 51 results for author: Bhatnagar, S

Polystyrene based scintillation detector studies in Astro particle Physics and medical physics

Authors: Sonali Bhatnagar

Abstract: DEASA (Dayalbagh Educational Air Shower Array) consists of eight plastic scintillators each with an area of 1 square meter. The cosmic ray showers have been simulated in CORSIKA [1] for the different primary particles in the energy range of 1014- 1015 eV. The longitudinal and lateral profiles have been studied for Agra. The real-life applications of cosmic ray particles in space have been studied… ▽ More DEASA (Dayalbagh Educational Air Shower Array) consists of eight plastic scintillators each with an area of 1 square meter. The cosmic ray showers have been simulated in CORSIKA [1] for the different primary particles in the energy range of 1014- 1015 eV. The longitudinal and lateral profiles have been studied for Agra. The real-life applications of cosmic ray particles in space have been studied to protect the astronaut from the galactic cosmic rays [2]. A plastic scintillation detector is simulated in Geant4 to study applications in hadron and carbon ion therapy [3]. The proton and carbon beam are simulated through the tumour region to study the stopping power and depth dose distribution for different organs. The energy range for each study is optimized and the Bragg curve is then interpreted with Bragg peak position and range. △ Less

Submitted 26 January, 2024; originally announced January 2024.

Comments: 10 pages, 5 figures

MSC Class: 00A79 ACM Class: I.6

The MeerKAT Absorption Line Survey (MALS) data release I: Stokes I image catalogs at 1-1.4 GHz

Authors: P. P. Deka, N. Gupta, P. Jagannathan, S. Sekhar, E. Momjian, S. Bhatnagar, J. Wagenveld, H. -R. Klöckner, J. Jose, S. A. Balashev, F. Combes, M. Hilton, D. Borgaonkar, A. Chatterjee, K. L. Emig, A. N. Gaunekar, G. I. G. Józsa, D. Y. Klutse, K. Knowles, J-. K. Krogager, A. Mohapatra, K. Moodley, Sébastien Muller, P. Noterdaeme, P. Petitjean , et al. (2 additional authors not shown)

Abstract: The MeerKAT Absorption Line Survey (MALS) has observed 391 telescope pointings at L-band (900 - 1670 MHz) at $δ\lesssim$ $+20°$. 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 (… ▽ More The MeerKAT Absorption Line Survey (MALS) has observed 391 telescope pointings at L-band (900 - 1670 MHz) at $δ\lesssim$ $+20°$. 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 $μ$Jy beam$^{-1}$ (22 $μ$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 ($α$) 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 ($α<-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. △ Less

Submitted 23 August, 2023; originally announced August 2023.

Comments: 64 pages, 25 figures, accepted for publication in the ApJS (full version of the paper with complete tables is available at DR1 release notes)

Resolving the bow shock and tail of the cannonball pulsar PSR J0002+6216

Authors: P. Kumar, F. K. Schinzel, G. B. Taylor, M. Kerr, D. Castro, U. Rau, S. Bhatnagar

Abstract: 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 rad… ▽ More 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. △ Less

Submitted 18 February, 2023; v1 submitted 9 February, 2023; originally announced February 2023.

Comments: 15 pages, 7 figures Accepted for publication in ApJ

CASA, the Common Astronomy Software Applications for Radio Astronomy

Authors: THE CASA TEAM, Ben Bean, Sanjay Bhatnagar, Sandra Castro, Jennifer Donovan Meyer, Bjorn Emonts, Enrique Garcia, Robert Garwood, Kumar Golap, Justo Gonzalez Villalba, Pamela Harris, Yohei Hayashi, Josh Hoskins, Mingyu Hsieh, Preshanth Jagannathan, Wataru Kawasaki, Aard Keimpema, Mark Kettenis, Jorge Lopez, Joshua Marvil, Joseph Masters, Andrew McNichols, David Mehringer, Renaud Miel, George Moellenbrock , et al. (24 additional authors not shown)

Abstract: CASA, the Common Astronomy Software Applications, is the primary data processing software for the Atacama Large Millimeter/submillimeter Array (ALMA) and the Karl G. Jansky Very Large Array (VLA), and is frequently used also for other radio telescopes. The CASA software can handle data from single-dish, aperture-synthesis, and Very Long Baseline Interferometery (VLBI) telescopes. One of its core f… ▽ More CASA, the Common Astronomy Software Applications, is the primary data processing software for the Atacama Large Millimeter/submillimeter Array (ALMA) and the Karl G. Jansky Very Large Array (VLA), and is frequently used also for other radio telescopes. The CASA software can handle data from single-dish, aperture-synthesis, and Very Long Baseline Interferometery (VLBI) telescopes. One of its core functionalities is to support the calibration and imaging pipelines for ALMA, VLA, VLA Sky Survey (VLASS), and the Nobeyama 45m telescope. This paper presents a high-level overview of the basic structure of the CASA software, as well as procedures for calibrating and imaging astronomical radio data in CASA. CASA is being developed by an international consortium of scientists and software engineers based at the National Radio Astronomical Observatory (NRAO), the European Southern Observatory (ESO), the National Astronomical Observatory of Japan (NAOJ), and the Joint Institute for VLBI European Research Infrastructure Consortium (JIV-ERIC), under the guidance of NRAO. △ Less

Submitted 5 October, 2022; originally announced October 2022.

Comments: Accepted for publication in PASP (20 pages, 4 figures). Joint publication with CASA-VLBI paper

Direction Dependent Corrections in Polarimetric Radio Imaging III: A-to-Z Solver -- Modeling the full Jones antenna aperture illumination pattern

Authors: Srikrishna Sekhar, Preshanth Jagannathan, Brian Kirk, Sanjay Bhatnagar, Russ Taylor

Abstract: In this third paper of a series describing direction dependent corrections for polarimetric radio imaging, we present the the A-to-Z solver methodology to model the full Jones antenna aperture illumination pattern (AIP) with Zernike polynomials. In order to achieve thermal noise limited imaging with modern radio interferometers, it is necessary to correct for the instrumental effects of the antenn… ▽ More In this third paper of a series describing direction dependent corrections for polarimetric radio imaging, we present the the A-to-Z solver methodology to model the full Jones antenna aperture illumination pattern (AIP) with Zernike polynomials. In order to achieve thermal noise limited imaging with modern radio interferometers, it is necessary to correct for the instrumental effects of the antenna primary beam (PB) as a function of time, frequency, and polarization. The wideband AW projection algorithm enables those corrections provided an accurate model of the AIP is available. We present the A-to-Z solver as a more versatile algorithm for the modeling of the AIP. It employs the orthonormal circular Zernike polynomial basis to model the measured full Jones AIP. These full Jones models are then used to reconstruct the full Mueller AIP repsonse of an antenna, in principle accounting for all the off-axis leakage effects of the primary beam. The A-to-Z solver is general enough to accomodate any interferometer for which holographic measurements exist, we have successfully modelled the AIP of VLA, MeerKAT and ALMA as a demonstration of its versatility. We show that our models capture the PB morphology to high accuracy within the first 1-2 sidelobes, and show the viability of full Mueller gridding and deconvolution for any telescope given high quality holographic measurements. △ Less

Submitted 4 December, 2021; v1 submitted 21 July, 2021; originally announced July 2021.

Comments: 16 pages, 12 figures. Accepted for publication in AJ

Blind HI and OH absorption line search: first results with MALS and uGMRT processed using ARTIP

Authors: N. Gupta, P. Jagannathan, R. Srianand, S. Bhatnagar, P. Noterdaeme, F. Combes, P. Petitjean, J. Jose, S. Pandey, C. Kaski, A. J. Baker, S. A. Balashev, E. Boettcher, H. -W. Chen, C. Cress, R. Dutta, S. Goedhart, G. Heald, G. I. G. Józsa, E. Kamau, P. Kamphuis, J. Kerp, H. -R. Klöckner, K. Knowles, V. Krishnan , et al. (10 additional authors not shown)

Abstract: 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 col… ▽ More 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 $<ρ<$ 150 kpc) to be less than 0.022. Due to the small redshift path, $Δ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 $ρ>30$ kpc. MALS with the expected $Δ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$. △ Less

Submitted 17 November, 2020; v1 submitted 8 July, 2020; originally announced July 2020.

Comments: 30 pages, 15 figures, accepted in ApJ

The CASA software for radio astronomy: status update from ADASS 2019

Authors: B. Emonts, R. Raba, G. Moellenbrock, S. Castro, C. E. Garcia-Dabo, J. Donovan Meyer, P. Ford, R. Garwood, K. Golap, J. Gonzalez, W. Kawasaki, A. McNichols, D. Mehringer, R. Miel, F. Montesino Pouzols, T. Nakazato, S. Nishie, J. Ott, D. Petry, U. Rau, C. Reynolds, D. Schiebel, N. Schweighart, J. -W. Steeb, V. Suoranta , et al. (6 additional authors not shown)

Abstract: 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 sup… ▽ More 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. △ Less

Submitted 19 December, 2019; originally announced December 2019.

Comments: 4 pages, to appear in proceedings of ADASS XXIX, ASP Conf. Series

First Time-Resolved CCD Photometry and Time-series Analysis of NSV 13601

Authors: Virabhadrasinh A. Gohil, Sachchidanand Prakash Bhatnagar

Abstract: 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,… ▽ More 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. △ Less

Submitted 23 October, 2019; originally announced October 2019.

Journal ref: RAA-2019-0178.R2

The DSA-2000 -- A Radio Survey Camera

Authors: G. Hallinan, V. Ravi, S. Weinreb, J. Kocz, Y. Huang, D. P. Woody, J. Lamb, L. D'Addario, M. Catha, J. Shi, C. Law, S. R. Kulkarni, E. S. Phinney, M. W. Eastwood, K. L. Bouman, M. A. McLaughlin, S. M. Ransom, X. Siemens, J. M. Cordes, R. S. Lynch, D. L. Kaplan, S. Chatterjee, J. Lazio, A. Brazier, S. Bhatnagar , et al. (3 additional authors not shown)

Abstract: 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… ▽ More 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 △ Less

Submitted 17 July, 2019; originally announced July 2019.

Comments: Submitted as a Project White Paper for the Astro2020 Decadal Survey (12 pages, 7 figures, 2 table)

The Karl G. Jansky Very Large Array Sky Survey (VLASS). Science case and survey design

Authors: M. Lacy, S. A. Baum, C. J. Chandler, S. Chatterjee, T. E. Clarke, S. Deustua, J. English, J. Farnes, B. M. Gaensler, N. Gugliucci, G. Hallinan, B. R. Kent, A. Kimball, C. J. Law, T. J. W. Lazio, J. Marvil, S. A. Mao, D. Medlin, K. Mooley, E. J. Murphy, S. Myers, R. Osten, G. T. Richards, E. Rosolowsky, L. Rudnick , et al. (53 additional authors not shown)

Abstract: 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$σ$ goal of 70 $μ$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. VLAS… ▽ More 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$σ$ goal of 70 $μ$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. △ Less

Submitted 30 December, 2019; v1 submitted 3 July, 2019; originally announced July 2019.

Comments: 34 pages, accepted by PASP (modified from prior version to address referee's and coauthor comments). (v2) Minor fixes to author list

Time-resolved CCD photometry and time-series analysis of RR Lyrae type RR Gem star

Authors: V. A. Gohil, S. P. Bhatnagar

Abstract: 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 cur… ▽ More 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. △ Less

Submitted 25 April, 2019; originally announced April 2019.

Comments: 8 pages, 9 figures

Report number: MS 4367 (RAA-2019-0081)

Journal ref: RAA 2019 Vol. 19 No. 9, 125(8pp)

The Tail of PSR J0002+6216 and the Supernova Remnant CTB 1

Authors: F. K. Schinzel, M. Kerr, U. Rau, S. Bhatnagar, D. A. Frail

Abstract: 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… ▽ More 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 $θ_μ=113^\circ$. We measure a proper motion with 2.9$σ$ significance from a Fermi timing analysis giving $μ$=115$\pm$33 mas yr$^{-1}$ and $θ_μ=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. △ Less

Submitted 16 April, 2019; originally announced April 2019.

Comments: 14 pages, 4 figures, accepted for publication in ApJL

Direction Dependent Corrections in Polarimetric Radio Imaging II: A-Solver Methodology A low-order solver for the A-Term of the A-Projection algorithm

Authors: Preshanth Jagannathan, S. Bhatnagar, W. Brisken, A. R. Taylor

Abstract: 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 th… ▽ More 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. △ Less

Submitted 2 November, 2017; originally announced November 2017.

Comments: 13pages, 10 Figures, Accepted for publication at the Astronomical Journal

The Pointing Self Calibration algorithm for aperture synthesis radio telescopes

Authors: S. Bhatnagar, T. J. Cornwell

Abstract: 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… ▽ More 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. △ Less

Submitted 25 September, 2017; originally announced September 2017.

Comments: 11 pages, 5 figures, accepted for publication in The Astronomical Journal

The MeerKAT Absorption Line Survey (MALS)

Authors: N. Gupta, R. Srianand, W. Baan, A. Baker, R. Beswick, S. Bhatnagar, D. Bhattacharya, A. Bosma, C. Carilli, M. Cluver, F. Combes, C. Cress, R. Dutta, J. Fynbo, G. Heald, M. Hilton, T. Hussain, M. Jarvis, G. Jozsa, P. Kamphuis, A. Kembhavi, J. Kerp, H. -R. Klöckner, J. Krogager, V. Kulkarni , et al. (19 additional authors not shown)

Abstract: 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 ca… ▽ More 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. △ Less

Submitted 24 August, 2017; originally announced August 2017.

Comments: 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

Direction Dependent Corrections in Polarimetric Radio Imaging I : Characterizing the effects of the primary beam on full Stokes imaging

Authors: P. Jagannathan, S. Bhatnagar, U. Rau, A. R. Taylor

Abstract: 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… ▽ More 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. △ Less

Submitted 9 June, 2017; v1 submitted 5 June, 2017; originally announced June 2017.

Comments: version2, 13 pages, 12 figures, accepted for publication in The Astronomical Journal

Understanding Systematic Errors Through Modeling of ALMA Primary Beams

Authors: Kara Kundert, Urvashi Rau, Edwin Bergin, Sanjay Bhatnagar

Abstract: 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… ▽ More 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. △ Less

Submitted 3 February, 2017; v1 submitted 13 January, 2017; originally announced January 2017.

Comments: 13 pages, 5 figures

VLA and ALMA Imaging of Intense, Galaxy-Wide Star Formation in z ~ 2 Galaxies

Authors: W. Rujopakarn, J. S. Dunlop, G. H. Rieke, R. J. Ivison, A. Cibinel, K. Nyland, P. Jagannathan, J. D. Silverman, D. M. Alexander, A. D. Biggs, S. Bhatnagar, D. R. Ballantyne, M. Dickinson, D. Elbaz, J. E. Geach, C. C. Hayward, A. Kirkpatrick, R. J. McLure, M. J. Michalowski, N. A. Miller, D. Narayanan, F. N. Owen, M. Pannella, C. Papovich, A. Pope , et al. (8 additional authors not shown)

Abstract: 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 $λ= 5$ cm and 1.3 mm using the Karl G. Jansky Very Large Array (VLA) and Atacama Large Millimeter/submillimeter Arr… ▽ More 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 $λ= 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. △ Less

Submitted 27 October, 2016; v1 submitted 26 July, 2016; originally announced July 2016.

Comments: 11 pages, 4 figures, accepted for publication in ApJ

Efficient implementation of the adaptive scale pixel decomposition algorithm

Authors: L. Zhang, S. Bhatnagar, U. Rau, M. Zhang

Abstract: 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 brig… ▽ More 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. △ Less

Submitted 15 December, 2017; v1 submitted 25 June, 2016; originally announced June 2016.

Comments: 6 pages; 4 figures

Journal ref: A&A 592, A128 (2016)

Deep wideband single pointings and mosaics in radio interferometry - How accurately do we reconstruct intensities and spectral indices of faint sources?

Authors: Urvashi Rau, Sanjay Bhatnagar, Frazer N. Owen

Abstract: 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… ▽ More 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μ$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. △ Less

Submitted 24 May, 2016; originally announced May 2016.

Comments: 12 pages, 12 figures (accepted to AJ)

The Population of Compact Radio Sources in the Orion Nebula Cluster

Authors: Jan Forbrich, Victor M. Rivilla, Karl M. Menten, Mark J. Reid, Claire J. Chandler, Urvashi Rau, Sanjay Bhatnagar, Scott J. Wolk, Stefan Meingast

Abstract: 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 pr… ▽ More 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. △ Less

Submitted 17 March, 2016; originally announced March 2016.

Comments: ApJ, in press

An Overview of the 2014 ALMA Long Baseline Campaign

Authors: ALMA Partnership, E. B. Fomalont, C. Vlahakis, S. Corder, A. Remijan, D. Barkats, R. Lucas, T. R. Hunter, C. L. Brogan, Y. Asaki, S. Matsushita, W. R. F. Dent, R. E. Hills, N. Phillips, A. M. S. Richards, P. Cox, R. Amestica, D. Broguiere, W. Cotton, A. S. Hales, R. Hiriart, A. Hirota, J. A. Hodge, C. M. V. Impellizzeri, J. Kern , et al. (224 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 24 April, 2015; v1 submitted 19 April, 2015; originally announced April 2015.

Comments: 11 pages, 7 figures, 2 tables; accepted for publication in the Astrophysical Journal Letters; this version with small changes to affiliations

Journal ref: The Astrophysical Journal Letters, volume 808, issue 1, article id. L1, 11 pp. (2015)

Direction Dependent Effects In Wide-Field Wideband Full Stokes Radio Imaging

Authors: Preshanth Jagannathan, Sanjay Bhatnagar, Urvashi Rau, Russ Taylor

Abstract: 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 tha… ▽ More 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. △ Less

Submitted 17 November, 2014; v1 submitted 4 November, 2014; originally announced November 2014.

Comments: 4 pages, 2 Figures, Astronomical Data Analysis Software & Systems 2014, Calgary, Canada

LOFAR low-band antenna observations of the 3C295 and Bootes fields: source counts and ultra-steep spectrum sources

Authors: R. J. van Weeren, W. L. Williams, C. Tasse, H. J. A. Rottgering, D. A. Rafferty, S. van der Tol, G. Heald, G. J. White, A. Shulevski, P. Best, H. T. Intema, S. Bhatnagar, W. Reich, M. Steinmetz, S. van Velzen, T. A. Ensslin, I. Prandoni, F. de Gasperin, M. Jamrozy, G. Brunetti, M. J. Jarvis, J. P. McKean, M. W. Wise, C. Ferrari, J. Harwood , et al. (76 additional authors not shown)

Abstract: 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 dif… ▽ More 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 ($α< -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ötes 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. △ Less

Submitted 18 September, 2014; originally announced September 2014.

Comments: 26 pages, 20 figures, ApJ in press

Wideband VLA Observations of Abell 2256 I: Continuum, Rotation Measure and Spectral Imaging

Authors: Frazer Owen, Lawrence Rudnick, Jean Eilek, Urvashi Rau, Sanjay Bhatnagar, Leonid Kogan

Abstract: 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… ▽ More 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. △ Less

Submitted 25 August, 2014; originally announced August 2014.

Comments: 15 pages, 22 figures. Accepted from publication in the Astrophysical Journal. PDF of paper with full resolution figures available at http://www.aoc.nrao.edu/~fowen/papers/a2256/a2256paper1.pdf.gz

The Deep Full-Stokes Radio Sky

Authors: A. R. Taylor, S. Bhatnagar, J. Condon, D. A. Green, J. M. Stil, P. Jagannathan, N. Kantharia, R. Kothes, R. Perley, J. Wall, T. Willis

Abstract: 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 $μ$Jy, these observations explore the radio source population at flux densities well below the regime dominated by classical radio g… ▽ More 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 $μ$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 $μ$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. △ Less

Submitted 1 May, 2014; originally announced May 2014.

Observing the Sun with the Murchison Widefield Array

Authors: D. Oberoi, R. Sharma, S. Bhatnagar, C. J. Lonsdale, L. D. Matthews, I. H. Cairns, S. J. Tingay, L. Benkevitch, A. Donea, S. M. White, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, B. E. Corey, A. Deshpande, D. Emrich, B. M. Gaensler, R. Goeke, L. J. Greenhill, B. J. Hazelton, M. Johnston-Hollitt, D. L. Kaplan, J. C. Kasper, E. Kratzenberg , et al. (20 additional authors not shown)

Abstract: 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 i… ▽ More 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. △ Less

Submitted 25 March, 2014; originally announced March 2014.

Comments: 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

Wideband Mosaic Imaging with the VLA - quantifying faint source imaging accuracy

Authors: Urvashi Rau, Sanjay Bhatnagar, Frazer N. Owen

Abstract: 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 re… ▽ More 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. △ Less

Submitted 20 March, 2014; originally announced March 2014.

Comments: 4 pages 1 figure

A JVLA 10~degree^2 deep survey

Authors: Matt J. Jarvis, Sanjay Bhatnagar, Marcus Bruggen, Chiara Ferrari, Ian Heywood, Martin Hardcastle, Eric Murphy, Russ Taylor, Oleg Smirnov, Chris Simpson, Vernesa Smolcic, Jeroen Stil, Kurt van der Heyden

Abstract: (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 ga… ▽ More (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. △ Less

Submitted 16 January, 2014; originally announced January 2014.

Comments: 18 pages, 6 figures. White Paper submitted to the call for Very Large Array Sky Surveys

Wide-field wide-band interferometric imaging:The WB A-Projection and hybrid algorithms

Authors: S. Bhatnagar, U. Rau, K. Golap

Abstract: 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… ▽ More 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. △ Less

Submitted 17 April, 2013; originally announced April 2013.

Comments: 10 pages, 5 figures; accepted for publication in ApJ

Applying full polarization A-Projection to very wide field of view instruments: An imager for LOFAR

Authors: C. Tasse, B. van der Tol, J. van Zwieten, Ger van Diepen, S. Bhatnagar

Abstract: 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… ▽ More 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. △ Less

Submitted 26 December, 2012; originally announced December 2012.

Radio Continuum Surveys with Square Kilometre Array Pathfinders

Authors: Ray P. Norris, J. Afonso, D. Bacon, Rainer Beck, Martin Bell, R. J. Beswick, Philip Best, Sanjay Bhatnagar, Annalisa Bonafede, Gianfranco Brunetti, Tamas Budavari, Rossella Cassano, J. J. Condo, Catherine Cress, Arwa Dabbech, I. Feain, Rob Fender, Chiara Ferrari, B. M. Gaensler, G. Giovannini, Marijke Haverkorn, George Heald, Kurt van der Heyden, A. M. Hopkins, M. Jarvis , et al. (26 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 28 October, 2012; originally announced October 2012.

Comments: Accepted by PASA, 22 October 2012

EVLA Observations of Galactic Supernova Remnants: Wide-field Continuum and Spectral-index Imaging

Authors: S. Bhatnagar, U. Rau, D. A. Green, M. P. Rupen

Abstract: 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 regi… ▽ More 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. △ Less

Submitted 15 June, 2011; v1 submitted 14 June, 2011; originally announced June 2011.

Comments: 6 pages; 4 figures; accepted for publication in ApJL, to appear in the ApJL EVLA special issue

Detection of Signals from Cosmic Reionization using Radio Interferometric Signal Processing

Authors: A. Datta, S. Bhatnagar, C. L. Carilli

Abstract: 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 t… ▽ More 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). △ Less

Submitted 18 August, 2009; originally announced August 2009.

Comments: accepted by ApJ; 12 pages, 10 figures

Calibration and imaging challenges at low radio frequencies: An overview of the state of the art

Authors: S. Bhatnagar

Abstract: 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$μ$Jy/beam corresponding to typical imag… ▽ More 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$μ$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. △ Less

Submitted 2 June, 2009; originally announced June 2009.

Comments: 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

Nature of X-shaped sources

Authors: D. V. Lal, M. J. Hardcastle, R. P. Kraft, C. C. Cheung, A. P. Lobanov, J. A. Zensus, S. Bhatnagar, A. P. Rao

Abstract: 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… ▽ More 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. △ Less

Submitted 17 April, 2009; originally announced April 2009.

Comments: 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

Advances in Calibration and Imaging Techniques in Radio Interferometry

Authors: U. Rau, S. Bhatnagar, M. A. Voronkov, T. J. Cornwell

Abstract: 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. 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. △ Less

Submitted 26 November, 2009; v1 submitted 4 February, 2009; originally announced February 2009.

Journal ref: Proceedings of the IEEE, Vol.97, No.8, p-1472, August 2009

GMRT study of X-shaped radio sources

Authors: Dharam Vir Lal, A. Pramesh Rao, Martin J. Hardcastle, C. C. Cheung, Sanjay Bhatnagar, Ralph P. Kraft, Andrei P. Lobanov, Anton J. Zensus

Abstract: 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 di… ▽ More 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. △ Less

Submitted 6 October, 2008; originally announced October 2008.

Comments: To appear in the proceedings of "From Planets to Dark Energy: the Modern Radio Universe", eds. R. Beswick et al., published by PoS

The non-coplanar baselines effect in radio interferometry: The W-Projection algorithm

Authors: T. J. Cornwell, K. Golap, S. Bhatnagar

Abstract: 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 super… ▽ More 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. △ Less

Submitted 25 July, 2008; originally announced July 2008.

Comments: Accepted for publication in "IEEE Journal of Selected Topics in Signal Processing"

Correcting direction-dependent gains in the deconvolution of radio interferometric images

Authors: S. Bhatnagar, T. J. Cornwell, K. Golap, Juan M. Uson

Abstract: 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… ▽ More 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. △ Less

Submitted 3 June, 2008; v1 submitted 7 May, 2008; originally announced May 2008.

Comments: submitted to A&A; some clarifications added in the text; accepted for publication in A&A

A Faint, Steep Spectum Burst from the Radio Transient GCRT J1745-3009

Authors: Scott D. Hyman, Subhashis Roy, Sabyasachi Pal, T. Joseph W. Lazio, Paul S. Ray, Namir E. Kassim, Sanjay Bhatnagar

Abstract: 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 propert… ▽ More 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. △ Less

Submitted 4 January, 2007; originally announced January 2007.

Comments: 11 pages including 2 figures; submitted to ApJ Letters

Radio and near-infrared observations of the steep spectrum Galactic plane radio source WKB 0314+57.8

Authors: D. A. Green, M. Lacy, S. Bhatnagar, E. L. Gates, P. J. Warner

Abstract: 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-… ▽ More 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. △ Less

Submitted 26 August, 2004; v1 submitted 6 August, 2004; originally announced August 2004.

Comments: 6 pages, to appear in MNRAS, typos corrected

Journal ref: Mon.Not.Roy.Astron.Soc. 354 (2004) 1159-1164

Scale sensitive deconvolution of interferometric images I. Adaptive Scale Pixel (Asp) decomposition

Authors: S. Bhatnagar, T. J. Cornwell

Abstract: 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 sig… ▽ More 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. △ Less

Submitted 12 July, 2004; originally announced July 2004.

Comments: 9 pages, 5 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: Astron.Astrophys. 426 (2004) 747-754

The late time radio emission from SN 1993J at meter wavelengths

Authors: P. Chandra, A. Ray, S. Bhatnagar

Abstract: 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… ▽ More 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 $α\sim 0.8-1.0$ at few tens of days to $\sim 0.6$ in about 10 years. △ Less

Submitted 24 May, 2004; originally announced May 2004.

Comments: 37 pages, 9 figures in LaTex; scheduled for ApJ 10 September 2004, v612 issue; send comments to: poonam@tifr.res.in

Journal ref: Astrophys.J. 612 (2004) 974-987

Synchrotron aging and the radio spectrum of SN 1993J

Authors: Poonam Chandra, Alak Ray, Sanjay Bhatnagar

Abstract: 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, $Δα\sim 0.6$, in the optically thin part of the spe… ▽ More 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, $Δα\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. △ Less

Submitted 16 February, 2004; originally announced February 2004.

Comments: 12 pages, 2 figures. Accepted for publication in ApJL

Journal ref: Astrophys.J. 604 (2004) L97-L100

Low frequency radio and X-ray properties of core-collapse supernovae

Authors: A. Ray, P. Chandra, F. Sutaria, S. Bhatnagar

Abstract: 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 des… ▽ More 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. △ Less

Submitted 18 November, 2003; originally announced November 2003.

Comments: 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

Low frequency observations of SN 1993J with Giant Meterwave Radio Telescope

Authors: P. Chandra, A. Ray, S. Bhatnagar

Abstract: 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… ▽ More 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. △ Less

Submitted 18 November, 2003; originally announced November 2003.

Comments: 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

The nature of the prompt X-ray and radio emission from SN2002ap

Authors: F. K. Sutaria, P. Chandra, S. Bhatnagar, A. Ray

Abstract: 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 spectru… ▽ More 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 $α= 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. △ Less

Submitted 6 November, 2002; v1 submitted 5 July, 2002; originally announced July 2002.

Comments: Replaced with the version accepted for publication in Astron.& Astrophys. (main journal). The text has been revised in places, but our conclusions remain unchanged

Journal ref: Astron.Astrophys. 397 (2003) 1011-1018

Radio Studies of Young Core Collapse Supernovae

Authors: Poonam Chandra, Alak Ray, Sanjay Bhatnagar

Abstract: 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 SN199… ▽ More 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. △ Less

Submitted 16 May, 2002; originally announced May 2002.

Comments: 4 pages, 2 .ps figs, 2 tables: uses basi.sty (enclosed). Presented at the ASI meeting at Pune Feb 5-8, 2002; Bulletin of Astronomical Society of India (2002)

327-MHz GMRT observations of one candidate and three Galactic Supernova Remnants

Authors: Sanjay Bhatnagar

Abstract: 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… ▽ More 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. △ Less

Submitted 10 December, 2001; originally announced December 2001.

Comments: 7 pages, 6 Postscript figures, accepted for publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc. 332 (2002) 1