-
Rotation Measure study of FRB 20180916B with the uGMRT
Authors:
S. Bethapudi,
L. G. Spitler,
D. Z. Li,
V. R. Marthi,
M. Bause,
R. A. Main,
R. S. Wharton
Abstract:
Context. Fast Radio Burst 20180916B is a repeating FRB whose activity window has a 16.34 day periodicity that also shifts and varies in duration with the observing frequency. Recently, arxiv:2205.09221 reported the FRB has started to show secular Rotation Measure (RM) increasing trend after only showing stochastic variability around a constant value of $-114.6$ rad m$^{-2}$ since its discovery.…
▽ More
Context. Fast Radio Burst 20180916B is a repeating FRB whose activity window has a 16.34 day periodicity that also shifts and varies in duration with the observing frequency. Recently, arxiv:2205.09221 reported the FRB has started to show secular Rotation Measure (RM) increasing trend after only showing stochastic variability around a constant value of $-114.6$ rad m$^{-2}$ since its discovery.
Aims. We aim to further study the RM variability of FRB 20180916B. The data comes from the ongoing campaigns of FRB 20180916B using the upgraded Giant Metrewave Radio Telescope (uGMRT). The majority of the observations are in Band 4, which is centered at 650 MHz with 200 MHz bandwidth.
Methods. We apply a standard single pulse search pipeline to search for bursts. In total, we detect 116 bursts with $\sim$36 hours of on-source time spanning 1200 days, with two bursts detected during simultaneous frequency coverage observations. We develop and apply a polarization calibration strategy suited for our dataset. On the calibrated bursts, we use QU-fitting to measure RM. Lastly, we also measure various other properties such as rate, linear polarization fraction and fluence distribution.
Results. Of the 116 detected bursts, we could calibrate 79 of them. From which, we observed in our early observations the RM continued to follow linear trend as modeled by arxiv:2205.09221. However, our later observations suggest the source switch from the linear trend to stochastic variations around a constant value of $-58.75$ rad m$^{-2}$. We also study cumulative rate against fluence and note that rate at higher fluences (> 1.2 Jy ms) scales as $γ= -1.09(7)$ whereas that at lower fluences (between 0.2 and 1.2 Jy ms) only scales as $γ= -0.51(1)$, meaning rate at higher fluence regime is steeper than at lower fluence regime.
△ Less
Submitted 19 September, 2024;
originally announced September 2024.
-
Identifying pulsar candidates in interferometric radio images using scintillation
Authors:
Jitendra Salal,
Shriharsh P. Tendulkar,
Visweshwar Ram Marthi
Abstract:
Pulsars have been primarily detected by their narrow pulses or periodicity in time domain data. Interferometric surveys for pulsars are challenging due to the trade-off between beam sensitivity and beam size and the corresponding tradeoff between survey sensitivity (depth), sky coverage, and computational efforts. The detection sensitivity of time-domain searches for pulsars is affected by dispers…
▽ More
Pulsars have been primarily detected by their narrow pulses or periodicity in time domain data. Interferometric surveys for pulsars are challenging due to the trade-off between beam sensitivity and beam size and the corresponding tradeoff between survey sensitivity (depth), sky coverage, and computational efforts. The detection sensitivity of time-domain searches for pulsars is affected by dispersion smearing, scattering, and rapid orbital motion of pulsars in binaries. We have developed a new technique to select pulsar candidates in interferometric radio images by identifying scintillating sources and measuring their scintillation bandwidth and timescale. Identifying likely candidates allows sensitive, focused time-domain searches, saving computational effort. Pulsar scintillation is independent of its timing properties and hence offers a different selection of pulsars compared to time-domain searches. Candidates identified from this method could allow us to find hard-to-detect pulsars, such as sub-millisecond pulsars and pulsars in very compact, highly-accelerated binary orbits. We use uGMRT observations in the fields of PSR\,B1508+55, PSR\,J0437$-$4715, and PSR\,B0031$-$07 as test cases for our technique. We demonstrate that the technique correctly differentiates between the pulsar and other non-scintillating point sources and show that the extracted dynamic spectrum of the pulsar is equivalent to that extracted from the uGMRT phased array beam. We show the results from our analysis of known pulsar fields and discuss challenges in dealing with interference and instrumental effects.
△ Less
Submitted 25 July, 2024;
originally announced July 2024.
-
Simultaneous and panchromatic observations of the Fast Radio Burst FRB 20180916B
Authors:
M. Trudu,
M. Pilia,
L. Nicastro,
C. Guidorzi,
M. Orlandini,
L. Zampieri,
V. R. Marthi,
F. Ambrosino,
A. Possenti,
M. Burgay,
C. Casentini,
I. Mereminskiy,
V. Savchenko,
E. Palazzi,
F. Panessa,
A. Ridolfi,
F. Verrecchia,
M. Anedda,
G. Bernardi,
M. Bachetti,
R. Burenin,
A. Burtovoi,
P. Casella,
M. Fiori,
F. Frontera
, et al. (25 additional authors not shown)
Abstract:
Aims. Fast Radio Bursts are bright radio transients whose origin has not yet explained. The search for a multi-wavelength counterpart of those events can put a tight constrain on the emission mechanism and the progenitor source. Methods. We conducted a multi-wavelength observational campaign on FRB 20180916B between October 2020 and August 2021 during eight activity cycles of the source. Observati…
▽ More
Aims. Fast Radio Bursts are bright radio transients whose origin has not yet explained. The search for a multi-wavelength counterpart of those events can put a tight constrain on the emission mechanism and the progenitor source. Methods. We conducted a multi-wavelength observational campaign on FRB 20180916B between October 2020 and August 2021 during eight activity cycles of the source. Observations were led in the radio band by the SRT both at 336 MHz and 1547 MHz and the uGMRT at 400 MHz. Simultaneous observations have been conducted by the optical telescopes Asiago (Galileo and Copernico), CMO SAI MSU, CAHA 2.2m, RTT-150 and TNG, and X/Gamma-ray detectors on board the AGILE, Insight-HXMT, INTEGRAL and Swift satellites. Results. We present the detection of 14 new bursts detected with the SRT at 336 MHz and seven new bursts with the uGMRT from this source. We provide the deepest prompt upper limits in the optical band fro FRB 20180916B to date. In fact, the TNG/SiFAP2 observation simultaneous to a burst detection by uGMRT gives an upper limit E_optical / E_radio < 1.3 x 10^2. Another burst detected by the SRT at 336 MHz was also co-observed by Insight-HMXT. The non-detection in the X-rays yields an upper limit (1-30 keV band) of E_X-ray / E_radio in the range of (0.9-1.3) x 10^7, depending on which model is considered for the X-ray emission.
△ Less
Submitted 29 May, 2023;
originally announced May 2023.
-
Modelling Annual Scintillation Velocity Variations of FRB 20201124A
Authors:
R. A. Main,
S. Bethapudi,
V. R. Marthi,
M. L. Bause,
D. Z. Li,
H. -H. Lin,
L. G. Spitler,
R. S. Wharton
Abstract:
Compact radio sources exhibit scintillation, an interference pattern arising from propagation through inhomogeneous plasma, where scintillation patterns encode the relative distances and velocities of the source, scattering material, and Earth. In previous work, we showed that the scintillation velocity of the repeating fast radio burst FRB20201124A could be measured by correlating burst spectra p…
▽ More
Compact radio sources exhibit scintillation, an interference pattern arising from propagation through inhomogeneous plasma, where scintillation patterns encode the relative distances and velocities of the source, scattering material, and Earth. In previous work, we showed that the scintillation velocity of the repeating fast radio burst FRB20201124A could be measured by correlating burst spectra pairs, with low values of the scintillation velocity and scattering timescale suggesting scattering nearby the Earth at $\sim0.4\,$kpc. In this work, we have measured the scintillation velocity at 10 epochs spanning a year, observing an annual variation which strongly implies the screen is within the Milky Way. Modelling the annual variation with a 1D anisotropic or 2D isotropic screen results in a screen distance $d_{l} = 0.40\pm0.04\,$kpc or $d_{l} = 0.46\pm0.06\,$kpc from Earth respectively, possibly associated with material outside of the Local Bubble or the edge of the Orion-Eridanus Superbubble. Additional measurements particularly at times of low effective velocity will help probe changes in screen properties, and distinguish between screen models. Where scintillation of an FRB originates in its host galaxy or local environment, these techniques could be used to detect orbital motion, and probe the FRB's local ionised environment.
△ Less
Submitted 23 June, 2023; v1 submitted 9 December, 2022;
originally announced December 2022.
-
DM-power: an algorithm for high precision dispersion measure with application to fast radio bursts
Authors:
Hsiu-Hsien Lin,
Robert Main,
Ue-Li Pen,
Robert Wharton,
Marlon Luis Bause,
Suryarao Bethapudi,
Dongzi Li,
Fang Xi Lin,
Visweshwar Ram Marthi,
Laura G Spitler
Abstract:
We present DM-power, a new method for precisely determining the dispersion measure (DM) of radio bursts, and apply it to the Fast Radio Burst (FRB) source FRB~20180916B. Motivated by the complex structure on multiple time scales seen in FRBs, DM-power optimizes the DM by combining measurements at multiple Fourier frequencies in the power spectrum of the burst. By optimally weighting the measuremen…
▽ More
We present DM-power, a new method for precisely determining the dispersion measure (DM) of radio bursts, and apply it to the Fast Radio Burst (FRB) source FRB~20180916B. Motivated by the complex structure on multiple time scales seen in FRBs, DM-power optimizes the DM by combining measurements at multiple Fourier frequencies in the power spectrum of the burst. By optimally weighting the measurements at each Fourier frequency, DM-power finds a burst DM that effectively incorporates information on many different burst timescales. We validate this technique on simulated Gaussian pulse profiles with a precision down to $σ_{\rm DM} \sim 0.001~{\rm pc~cm}^{-3}$, and then apply it to bursts from pulsar B0329+54 and FRB~20180916B. The precision of these DM measurements are sufficient to measure a statistically significant variation in DM over a $\approx 2$ hr span. While this variation could be the result of electron density variations along the line of sight, it is more like that the observed variation is the result of intrinsic frequency-dependent burst structure that can mimic a dispersive delay.
△ Less
Submitted 21 June, 2023; v1 submitted 29 August, 2022;
originally announced August 2022.
-
Burst properties of the highly active FRB 20201124A using uGMRT
Authors:
V. R. Marthi,
S. Bethapudi,
R. A. Main,
H. -H. Lin,
L. G. Spitler,
R. S. Wharton,
D. Z. Li,
T. G. Gautam,
U. -L. Pen,
G. H. Hilmarsson
Abstract:
We report the observations of the highly active FRb20201124A with the upgraded Giant Metrewave Radio Telescope at 550-750~MHz. These observations in the incoherent array mode simultaneously provided an arcsecond localization of bursts from \rss, the discovery of persistent radio emission associated with the host galaxy, and the detection of 48 bursts. Using the brightest burst in the sample (…
▽ More
We report the observations of the highly active FRb20201124A with the upgraded Giant Metrewave Radio Telescope at 550-750~MHz. These observations in the incoherent array mode simultaneously provided an arcsecond localization of bursts from \rss, the discovery of persistent radio emission associated with the host galaxy, and the detection of 48 bursts. Using the brightest burst in the sample ($F= 108~{\rm Jy~ms}$) we find a structure-maximizing dispersion measure of $410.8 \pm 0.5~{\rm pc~cm}^{-3}$. We find that our observations are complete down to a fluence level of $10~{\rm Jy~ms}$, above which the cumulative burst rate scales as a power-law $R(>\!F) = 10~{\rm hr}^{-1} \left(F/10\mathrm{~Jy~ms}\right)^γ$ with $γ= -1.2 \pm 0.2$. We find that the bursts are on average wider than those reported for other repeating FRBs. We find that the waiting time between bursts is well approximated by an exponential distribution with a mean of $\sim 2.9$ min during our observations. We searched for periodicities using both a standard Fourier domain method and the Fast Folding Algorithm, but found no significant candidates. We measure bulk spectro-temporal drift rates between $-0.75$ and $-20~{\rm MHz~ms}^{-1}$. Finally, we use the brightest burst to set an upper limit to the scattering time of 11.1~ms at 550~MHz. The localization of FRB20201124A adds strength to the proof-of-concept method described in our earlier work and serves as a potential model for future localizations and follow-up of repeating FRBs with the uGMRT.
△ Less
Submitted 21 October, 2021; v1 submitted 2 August, 2021;
originally announced August 2021.
-
Scintillation timescale measurement of the highly active FRB20201124A
Authors:
R. A. Main,
G. H. Hilmarsson,
V. R. Marthi,
L. G. Spitler,
R. S. Wharton,
S. Bethapudi,
D. Z. Li,
H. -H. Lin
Abstract:
Scintillation of compact radio sources results from the interference between images caused by multipath propagation, and probes the intervening scattering plasma and the velocities of the emitting source and scattering screen. In FRB20201124A, a repeating fast radio burst (FRB) which entered a period of extreme activity, we obtained many burst detections in observations at the upgraded Giant Metre…
▽ More
Scintillation of compact radio sources results from the interference between images caused by multipath propagation, and probes the intervening scattering plasma and the velocities of the emitting source and scattering screen. In FRB20201124A, a repeating fast radio burst (FRB) which entered a period of extreme activity, we obtained many burst detections in observations at the upgraded Giant Metrewave Radio Telescope (uGMRT) and the Effelsberg 100\,m Radio Telescope. Bursts nearby in time show similar scintillation patterns, and we measure a scintillation timescale of $14.3\pm1.2$\,min and $7\pm2$\,min at Effelsberg (1370\,MHz) and uGMRT (650\,MHz), respectively, by correlating burst pair spectra. The scintillation bandwidth scaled to 1\,GHz is $0.5\pm0.1$\,MHz, and the inferred scintillation velocity at Effelsberg is $V_{\mathrm{ISS}}\approx (59\pm7) \sqrt{d_{l}/2\,\rm{kpc}}~{\rm km~s}^{-1}$, higher than Earth's velocity for any screen beyond a lens distance of $d_{l} \gtrsim 400\,$pc. From the measured scintillation bandwidth, FRB20201124A has comparatively lower scattering than nearby pulsars, and is underscattered by a factor of $\sim 30$ or $\sim 1200$ compared to the NE2001 and YMW16 model predictions respectively. This underscattering, together with the measured scintillation velocity are consistent with a scattering screen more nearby the Earth at $d_{l} \sim 400\,$pc, rather than at the 2\,kpc spiral arm which NE2001 predicts to be the dominant source of scattering. With future measurements, the distance, geometry, and velocity of the scattering screen could be obtained through modelling of the annual variation in $V_{\rm ISS}$, or through inter-station time delays or interferometric observations. Scintillation/scattering measurements of FRBs could help improve Galactic electron density models, particularly in the Galactic halo or at high Galactic latitudes.
△ Less
Submitted 22 November, 2021; v1 submitted 30 July, 2021;
originally announced August 2021.
-
Scintillation of PSR B1508+55 -- the view from a 10,000-km baseline
Authors:
V. R. Marthi,
D. Simard,
R. A. Main,
U. -L. Pen,
M. H. van Kerkwijk,
K. Vanderlinde,
Y. Gupta,
C. Roberts,
B. M. Quine
Abstract:
We report on the simultaneous Giant Metrewave Radio Telescope (GMRT) and Algonquin Radio Observatory (ARO) observations at 550-750 MHz of the scintillation of PSR B1508+55, resulting in a $\sim$10,000-km baseline. This regime of measurement lies between the shorter few 100-1000~km baselines of earlier multi-station observations and the much longer earth-space baselines. We measure a scintillation…
▽ More
We report on the simultaneous Giant Metrewave Radio Telescope (GMRT) and Algonquin Radio Observatory (ARO) observations at 550-750 MHz of the scintillation of PSR B1508+55, resulting in a $\sim$10,000-km baseline. This regime of measurement lies between the shorter few 100-1000~km baselines of earlier multi-station observations and the much longer earth-space baselines. We measure a scintillation cross-correlation coefficient of $0.22$, offset from zero time lag due to a $\sim 45$~s traversal time of the scintillation pattern. The scintillation time of 135~s is $3\times$ longer, ruling out isotropic as well as strictly 1D scattering. Hence, the low cross-correlation coefficient is indicative of highly anisotropic but 2D scattering. The common scintillation detected on the baseline is confined to low delays of $\lesssim 1 μ$s, suggesting that this correlation may not be associated with the parabolic scintillation arc detected at the GMRT. Detection of pulsed echoes and their direct imaging with the Low Frequency Array (LOFAR) by a different group enable them to measure a distance of 125~pc to the screen causing these echoes. These previous measurements, alongside our observations, lead us to propose that there are at least two scattering screens: the closer 125 pc screen causing the scintillation arc detected at GMRT, and a screen further beyond causing the scintillation detected on the GMRT-ARO baseline. We advance the hypothesis that the 125-pc screen partially resolves the speckle images on the screen beyond leading to loss of coherence in the scintillation dynamic spectrum, to explain the low cross-correlation coefficient.
△ Less
Submitted 8 July, 2021; v1 submitted 19 October, 2020;
originally announced October 2020.
-
Detection of 15 bursts from FRB 180916.J0158+65 with the uGMRT
Authors:
Visweshwar Ram Marthi,
Tasha Gautam,
Dongzi Li,
Hsiu-Hsien Lin,
Robert Main,
Arun Kumar Naidu,
Ue-Li Pen,
Robert Wharton
Abstract:
We report the findings of a uGMRT observing campaign on FRB 180916.J0158+65, discovered recently to show a 16.35-day periodicity of its active cycle. We observed the source at 550-750 MHz for $\sim 2$ hours each during three successive cycles at the peak of its expected active period. We find 0, 12, and 3 bursts respectively, implying a highly variable bursting rate even within the active phase. W…
▽ More
We report the findings of a uGMRT observing campaign on FRB 180916.J0158+65, discovered recently to show a 16.35-day periodicity of its active cycle. We observed the source at 550-750 MHz for $\sim 2$ hours each during three successive cycles at the peak of its expected active period. We find 0, 12, and 3 bursts respectively, implying a highly variable bursting rate even within the active phase. We consistently detect faint bursts with spectral energies only an order of magnitude higher than the Galactic burst source SGR~1935+2154. The times of arrival of the detected bursts rule out many possible aliased solutions, strengthening the findings of the 16.35-day periodicity. A short-timescale periodicity search returned no highly significant candidates. Two of the beamformer-detected bursts were bright enough to be clearly detected in the imaging data, achieving sub-arcsecond localization, and proving as a proof-of-concept for FRB imaging with the GMRT. We provide a $3σ$ upper limit of the persistent radio flux density at 650 MHz of $66~μ{\rm Jy}$ which, combined with the EVN and VLA limits at 1.6~GHz, further constrains any potential radio counterpart. These results demonstrate the power of uGMRT for targeted observations to detect and localize known repeating FRBs.
△ Less
Submitted 24 August, 2020; v1 submitted 28 July, 2020;
originally announced July 2020.
-
Simulated Predictions for HI at $z = 3.35$ with the Ooty Wide Field Array (OWFA) -- II: Foreground Avoidance
Authors:
Suman Chatterjee,
Somnath Bharadwaj,
Visweshwar Ram Marthi
Abstract:
Considering the upcoming OWFA, we use simulations of the foregrounds and the $z = 3.35$ HI 21-cm intensity mapping signal to identify the $(k_{\perp},k_{\parallel})$ modes where the expected 21-cm power spectrum $P(k_{\perp},k_{\parallel})$ is substantially larger than the predicted foreground contribution. Only these uncontaminated $k$-modes are used for measuring $P(k_{\perp},k_{\parallel})$ in…
▽ More
Considering the upcoming OWFA, we use simulations of the foregrounds and the $z = 3.35$ HI 21-cm intensity mapping signal to identify the $(k_{\perp},k_{\parallel})$ modes where the expected 21-cm power spectrum $P(k_{\perp},k_{\parallel})$ is substantially larger than the predicted foreground contribution. Only these uncontaminated $k$-modes are used for measuring $P(k_{\perp},k_{\parallel})$ in the "Foreground Avoidance" technique. Though the foregrounds are largely localised within a wedge. we find that the small leakage beyond the wedge surpasses the 21-cm signal across a significant part of the $(k_{\perp},k_{\parallel})$ plane. The extent of foreground leakage is extremely sensitive to the frequency window function used to estimate $P(k_{\perp},k_{\parallel})$. It is possible to reduce the leakage by making the window function narrower, however this comes at the expense of losing a larger fraction of the 21-cm signal. It is necessary to balance these competing effects to identify an optimal window function. Considering a broad class of cosine window functions, we identify a six term window function as optimal for 21-cm power spectrum estimation with OWFA. Considering only the $k$-modes where the expected 21-cm power spectrum exceeds the predicted foregrounds by a factor of $100$ or larger, a $5\,σ$ detection of the binned power spectrum is possible in the $k$ ranges $0.18 \leq k \leq 0.3 \, {\rm Mpc}^{-1}$ and $0.18 \le k \le 0.8 \, {\rm Mpc}^{-1}$ with $1,000-2,000$ hours and $10^4$ hours of observation respectively.
△ Less
Submitted 18 December, 2020; v1 submitted 13 November, 2019;
originally announced November 2019.
-
Disentangling interstellar plasma screens with pulsar VLBI: Combining auto- and cross-correlations
Authors:
Dana Simard,
Ue-Li Pen,
Visweshwar Ram Marthi,
Walter Brisken
Abstract:
Pulsar scintillation allows a glimpse into small-scale plasma structures in the interstellar medium, if we can infer their properties from the scintillation pattern. With Very Long Baseline Interferometry and working in delay-delay rate space, where the contributions of pairs of images to the interference pattern become localized, the scattering geometry and distribution of scattered images on the…
▽ More
Pulsar scintillation allows a glimpse into small-scale plasma structures in the interstellar medium, if we can infer their properties from the scintillation pattern. With Very Long Baseline Interferometry and working in delay-delay rate space, where the contributions of pairs of images to the interference pattern become localized, the scattering geometry and distribution of scattered images on the sky can be determined if a single, highly-anisotropic scattering screen is responsible for the scintillation. However, many pulsars are subject to much more complex scattering environments where this method cannot be used. We present a novel technique to reconstruct the scattered flux of the pulsar and solve for the scattering geometry in these cases by combining interferometric visibilities with cross-correlations of single-station intensities. This takes advantage of the fact that, considering a single image pair in delay-delay rate space, the visibilities are sensitive to the sum of the image angular displacements, while the cross-correlated intensities are sensitive to the difference, so that their combination can be used to localize both images of the pair. We show that this technique is able to reconstruct the published scattering geometry of PSR B0834+06, then apply it to simulations of more complicated scattering systems, where we find that it can distinguish features from different scattering screens even when the presence of multiple screens is not obvious in the Fourier transform of the dynamic spectrum. This technique will allow us to both better understand the distribution of scattering within the interstellar medium and to apply current scintillometry techniques, such as modelling scintillation and constraining the location of pulsar emission, to sources for which a current lack of understanding of the scattering environment precludes the use of these techniques. (abridged)
△ Less
Submitted 16 October, 2018;
originally announced October 2018.
-
An analytical method to simulate the HI 21-cm visibility signal for intensity mapping experiments
Authors:
Anjan Kumar Sarkar,
Somnath Bharadwaj,
Visweshwar Ram Marthi
Abstract:
Simulations play a vital role in testing and validating HI 21-cm power spectrum estimation techniques. Conventional methods use techniques like N-body simulations to simulate the sky signal which is then passed through a model of the instrument. This makes it necessary to simulate the HI distribution in a large cosmological volume, and incorporate both the light-cone effect and the telescope's chr…
▽ More
Simulations play a vital role in testing and validating HI 21-cm power spectrum estimation techniques. Conventional methods use techniques like N-body simulations to simulate the sky signal which is then passed through a model of the instrument. This makes it necessary to simulate the HI distribution in a large cosmological volume, and incorporate both the light-cone effect and the telescope's chromatic response. The computational requirements may be particularly large if one wishes to simulate many realizations of the signal. In this paper we present an analytical method to simulate the HI visibility signal. This is particularly efficient if one wishes to simulate a large number of realizations of the signal. Our method is based on theoretical predictions of the visibility correlation which incorporate both the light-cone effect and the telescope's chromatic response. We have demonstrated this method by applying it to simulate the HI visibility signal for the upcoming Ooty Wide Field Array Phase I.
△ Less
Submitted 12 September, 2017;
originally announced September 2017.
-
Simulated foreground predictions for HI at z = 3.35 with the Ooty Wide Field Array: I. Instrument and the foregrounds
Authors:
Visweshwar Ram Marthi,
Suman Chatterjee,
Jayaram Chengalur,
Somnath Bharadwaj
Abstract:
Foreground removal is the most important step in detecting the large-scale redshifted HI 21-cm signal. Modelling foreground spectra is challenging and is further complicated by the chromatic response of the telescope. We present a multi-frequency angular power spectrum (MAPS) estimator for use in a survey for redshifted HI 21-cm emission from z~3.35, and demonstrate its ability to accurately chara…
▽ More
Foreground removal is the most important step in detecting the large-scale redshifted HI 21-cm signal. Modelling foreground spectra is challenging and is further complicated by the chromatic response of the telescope. We present a multi-frequency angular power spectrum (MAPS) estimator for use in a survey for redshifted HI 21-cm emission from z~3.35, and demonstrate its ability to accurately characterize the foregrounds. This survey will be carried out with the two wide-field interferometer modes of the upgraded Ooty Radio Telescope, called the Ooty Wide Field Array (OWFA), at 326.5 MHz. We have tailored the two-visibility correlation for OWFA to estimate the MAPS and test it with simulated foregrounds. In the process, we describe a software model that encodes the geometry and the details of the telescope, and simulates a realistic model for the bright radio sky. This article presents simulations which include the full chromatic response of the telescope, in addition to the frequency dependence intrinsic to the foregrounds. We find that the visibility correlation MAPS estimator recovers the input angular power spectrum accurately, and that the instrument response to the foregrounds dominates the systematic errors in the recovered foreground power spectra.
△ Less
Submitted 17 July, 2017;
originally announced July 2017.
-
Prowess - a software model for the Ooty Wide Field Array
Authors:
Visweshwar Ram Marthi
Abstract:
One of the scientific objectives of the Ooty Wide Field Array(OWFA) is to observe the redshifted HI emission from z ~ 3.35. Although predictions spell out optimistic outcomes in reasonable integration times, these studies were based purely on analytical assumptions, without accounting for limiting systematics. A software model for OWFA has been developed with a view to understanding the instrument…
▽ More
One of the scientific objectives of the Ooty Wide Field Array(OWFA) is to observe the redshifted HI emission from z ~ 3.35. Although predictions spell out optimistic outcomes in reasonable integration times, these studies were based purely on analytical assumptions, without accounting for limiting systematics. A software model for OWFA has been developed with a view to understanding the instrument-induced systematics, by describing a complete software model for the instrument. This model has been implemented through a suite of programs, together called Prowess, which has been conceived with a dual role of an emulator as well as observatory data analysis software. The programming philosophy followed in building Prowess enables any user to define an own set of functions and add new functionality. This paper describes a co-ordinate system suitable for OWFA in which the baselines are defined. The foregrounds are simulated from their angular power spectra. The visibilities are then computed from the foregrounds. These visibilities are then used for further processing, such as calibration and power spectrum estimation. The package allows for rich visualisation features in multiple output formats in an interactive fashion, giving the user an intuitive feel for the data. Prowess has been extensively used for numerical predictions of the foregrounds for the OWFA HI experiment.
△ Less
Submitted 2 March, 2017;
originally announced March 2017.
-
Simulating the z = 3.35 HI 21-cm visibility signal for the Ooty Wide Field Array(OWFA)
Authors:
Suman Chatterjee,
S. Bharadwaj,
V. R. Marthi
Abstract:
The upcoming Ooty Wide Field Array (OWFA) will operate at $326.5 \, {\rm MHz}$ which corresponds to the redshifted 21-cm signal from neutral hydrogen (HI) at z = 3.35. We present two different prescriptions to simulate this signal and calculate the visibilities expected in radio-interferometric observations with OWFA. In the first method we use an input model for the expected 21-cm power spectrum…
▽ More
The upcoming Ooty Wide Field Array (OWFA) will operate at $326.5 \, {\rm MHz}$ which corresponds to the redshifted 21-cm signal from neutral hydrogen (HI) at z = 3.35. We present two different prescriptions to simulate this signal and calculate the visibilities expected in radio-interferometric observations with OWFA. In the first method we use an input model for the expected 21-cm power spectrum to directly simulate different random realizations of the brightness temperature fluctuations and calculate the visibilities. This method, which models the HI signal entirely as a diffuse radiation, is completely oblivious to the discrete nature of the astrophysical sources which host the HI. While each discrete source subtends an angle that is much smaller than the angular resolution of OWFA, the velocity structure of the HI inside the individual sources is well within reach of OWFA's frequency resolution and this is expected to have an impact on the observed HI signal. The second prescription is based on cosmological N-body simulations. Here we identify each simulation particle with a source that hosts the HI, and we have the freedom to implement any desired line profile for the HI emission from the individual sources. Implementing a simple model for the line profile, we have generated several random realizations of the complex visibilities. Correlations between the visibilities measured at different baselines and channels provides an unique method to quantify the statistical properties of the \HI signal. We have used this to quantify the results of our simulations, and explore the relation between the expected visibility correlations and the underlying HI power spectrum.
△ Less
Submitted 2 March, 2017;
originally announced March 2017.
-
Non-linear Redundancy Calibration
Authors:
Visweshwar Ram Marthi,
Jayaram Chengalur
Abstract:
For radio interferometric arrays with a sufficient number of redundant spacings the multiplicity of measurements of the same sky visibility can be used to determine both the antenna gains as well as the true visibilities. Many of the earlier approaches to this problem focused on linearized versions of the relation between the measured and the true visibilities. Here we propose to use a standard no…
▽ More
For radio interferometric arrays with a sufficient number of redundant spacings the multiplicity of measurements of the same sky visibility can be used to determine both the antenna gains as well as the true visibilities. Many of the earlier approaches to this problem focused on linearized versions of the relation between the measured and the true visibilities. Here we propose to use a standard non-linear minimization algorithm to solve for both the antenna gains as well as the true visibilities. We show through simulations done in the context of the ongoing upgrade to the Ooty Radio Telescope that the non-linear minimization algorithm is fast compared to the earlier approaches. Further, unlike the most straightforward linearized approach, which works with the logarithms of the visibilities and the gains, the non-linear minimization algorithm leads to unbiased solutions. Finally we present error estimates for the estimated gains and visibilities. Monte-Carlo simulations establish that the estimator is indeed statistically efficient, achieving the Cramer-Rao bound.
△ Less
Submitted 5 October, 2013;
originally announced October 2013.
-
The central point source in G76.9+1.0
Authors:
V. R. Marthi,
J. N. Chengalur,
Y. Gupta,
G. C. Dewangan,
D. Bhattacharya
Abstract:
We describe the serendipitous discovery of a radio point source in a 618 MHz image of the supernova remnant(SNR) G76.9+1.0. The SNR has a bipolar structure and the point source is located near a faint bridge of emission joining the two lobes of emission. The point source was also detected in follow-up higher frequency(1170 MHz) observations. The spectral index for the point source obtained from th…
▽ More
We describe the serendipitous discovery of a radio point source in a 618 MHz image of the supernova remnant(SNR) G76.9+1.0. The SNR has a bipolar structure and the point source is located near a faint bridge of emission joining the two lobes of emission. The point source was also detected in follow-up higher frequency(1170 MHz) observations. The spectral index for the point source obtained from the GMRT observations is alpha = -2.1. The steep spectrum, as well as the location of the point source near the centre of the SNR establish the fact that it is indeed the pulsar J2022+3842 associated with this SNR. Consistent with this, subsequent analysis of archival Chandra X-ray data shows a point source coincident with the radio point source, as well as diffuse extended X-ray emission surrounding the unresolved source. However, no pulsed emission was detected despite deep searches at both 610 MHz and 1160 MHz although pulsed emission has been seen at 2 GHz with the GBT. It appears that the most likely reason for not detecting the pulsed signal at the GMRT is temporal broadening: for the estimated DM towards this SNR, the pulse broadening time could be as large as tens of milliseconds. The diffuse X-ray emission is elongated along the same direction as the bipolar structure seen in the radio. We interpret the radio lobes as having been formed from an equatorial wind. Although direct detection of pulsed signal has not been possible, we show convincingly that sensitive, high-resolution, radio imaging at multiple frequencies is a useful method to search for pulsar candidates.
△ Less
Submitted 1 June, 2011;
originally announced June 2011.
-
A search for the 55 MHz OH line
Authors:
Visweshwar Ram Marthi,
Jayaram N. Chengalur
Abstract:
The OH molecule, found abundantly in the Milky Way, has four transitions at the ground state rotational level(J = 3/2) at cm wavelengths. These are E1 transitions between the F+ and F- hyperfine levels of the Lambda doublet of the J=3/2 state. There are also forbidden M1 transitions between the hyperfine levels within each of the doublet states occuring at frequencies 53.171 MHz and 55.128 MHz. Th…
▽ More
The OH molecule, found abundantly in the Milky Way, has four transitions at the ground state rotational level(J = 3/2) at cm wavelengths. These are E1 transitions between the F+ and F- hyperfine levels of the Lambda doublet of the J=3/2 state. There are also forbidden M1 transitions between the hyperfine levels within each of the doublet states occuring at frequencies 53.171 MHz and 55.128 MHz. These are extremely weak and hence difficult to detect. However there is a possibility that the level populations giving rise to these lines are inverted under special conditions, in which case it may be possible to detect them through their maser emission. We describe the observational diagnostics for determining when the hyperfine levels are inverted, and identify a region near W44 where these conditions are satisfied. A high-velocity-resolution search for these hyperfine OH lines using the low frequency feeds on four antennas of the GMRT and the new GMRT Software Backend(GSB) was performed on this target near W44. We place a 3-sigma upper limit of ~17.3 Jy (at 1 km/s velocity resolution) for the 55 MHz line from this region. This corresponds to an upper limit of 3 X 10^8 for the amplification of the Galactic synchrotron emission providing the background.
△ Less
Submitted 2 July, 2010; v1 submitted 26 April, 2010;
originally announced April 2010.
