Showing 1–50 of 83 results for author: Prabu, T

Low-frequency pulse-jitter measurement with the uGMRT I : PSR J0437$-$4715

Authors: Tomonosuke Kikunaga, Shinnosuke Hisano, Neelam Dhanda Batra, Shantanu Desai, Bhal Chandra Joshi, Manjari Bagchi, T. Prabu, Keitaro Takahashi, Swetha Arumugam, Adarsh Bathula, Subhajit Dandapat, Debabrata Deb, Churchil Dwivedi, Yashwant Gupta, Shebin Jose Jacob, Fazal Kareem, Nobleson K, Pragna Mamidipaka, Avinash Kumar Paladi, Arul Pandian B, Prerna Rana, Jaikhomba Singha, Aman Srivastava, Mayuresh Surnis, Pratik Tarafdar

Abstract: High-precision pulsar timing observations are limited in their accuracy by the jitter noise that appears in the arrival time of pulses. Therefore, it is important to systematically characterise the amplitude of the jitter noise and its variation with frequency. In this paper, we provide jitter measurements from low-frequency wideband observations of PSR J0437$-$4715 using data obtained as part of… ▽ More High-precision pulsar timing observations are limited in their accuracy by the jitter noise that appears in the arrival time of pulses. Therefore, it is important to systematically characterise the amplitude of the jitter noise and its variation with frequency. In this paper, we provide jitter measurements from low-frequency wideband observations of PSR J0437$-$4715 using data obtained as part of the Indian Pulsar Timing Array experiment. We were able to detect jitter in both the 300 - 500 MHz and 1260 - 1460 MHz observations of the upgraded Giant Metrewave Radio Telescope (uGMRT). The former is the first jitter measurement for this pulsar below 700 MHz, and the latter is in good agreement with results from previous studies. In addition, at 300 - 500 MHz, we investigated the frequency dependence of the jitter by calculating the jitter for each sub-banded arrival time of pulses. We found that the jitter amplitude increases with frequency. This trend is opposite as compared to previous studies, indicating that there is a turnover at intermediate frequencies. It will be possible to investigate this in more detail with uGMRT observations at 550 - 750 MHz and future high sensitive wideband observations from next generation telescopes, such as the Square Kilometre Array. We also explored the effect of jitter on the high precision dispersion measure (DM) measurements derived from short duration observations. We find that even though the DM precision will be better at lower frequencies due to the smaller amplitude of jitter noise, it will limit the DM precision for high signal-to-noise observations, which are of short durations. This limitation can be overcome by integrating for a long enough duration optimised for a given pulsar. △ Less

Submitted 18 April, 2024; v1 submitted 4 December, 2023; originally announced December 2023.

Comments: 13 pages, 12 figures, 3 tables, accepted for Publication of the Astronomical Society of Australia

Using low-frequency scatter-broadening measurements for precision estimates of dispersion measures

Authors: Jaikhomba Singha, Bhal Chandra Joshi, M. A. Krishnakumar, Fazal Kareem, Adarsh Bathula, Churchil Dwivedi, Shebin Jose Jacob, Shantanu Desai, Pratik Tarafdar, P. Arumugam, Swetha Arumugam, Manjari Bagchi, Neelam Dhanda Batra, Subhajit Dandapat, Debabrata Deb, Jyotijwal Debnath, A Gopakumar, Yashwant Gupta, Shinnosuke Hisano, Ryo Kato, Tomonosuke Kikunaga, Piyush Marmat, K. Nobleson, Avinash K. Paladi, Arul Pandian B. , et al. (6 additional authors not shown)

Abstract: A pulsar's pulse profile gets broadened at low frequencies due to dispersion along the line of sight or due to multi-path propagation. The dynamic nature of the interstellar medium makes both of these effects time-dependent and introduces slowly varying time delays in the measured times-of-arrival similar to those introduced by passing gravitational waves. In this article, we present a new method… ▽ More A pulsar's pulse profile gets broadened at low frequencies due to dispersion along the line of sight or due to multi-path propagation. The dynamic nature of the interstellar medium makes both of these effects time-dependent and introduces slowly varying time delays in the measured times-of-arrival similar to those introduced by passing gravitational waves. In this article, we present a new method to correct for such delays by obtaining unbiased dispersion measure (DM) measurements by using low-frequency estimates of the scattering parameters. We evaluate this method by comparing the obtained DM estimates with those, where scatter-broadening is ignored using simulated data. A bias is seen in the estimated DMs for simulated data with pulse-broadening with a larger variability for a data set with a variable frequency scaling index, $α$, as compared to that assuming a Kolmogorov turbulence. Application of the proposed method removes this bias robustly for data with band averaged signal-to-noise ratio larger than 100. We report, for the first time, the measurements of the scatter-broadening time and $α$ from analysis of PSR J1643$-$1224, observed with upgraded Giant Metrewave Radio Telescope as part of the Indian Pulsar Timing Array experiment. These scattering parameters were found to vary with epoch and $α$ was different from that expected for Kolmogorov turbulence. Finally, we present the DM time-series after application of the new technique to PSR J1643$-$1224. △ Less

Submitted 28 September, 2023; originally announced September 2023.

Comments: 9 pages, 8 figures, Submitted to MNRAS

Comparing recent PTA results on the nanohertz stochastic gravitational wave background

Authors: The International Pulsar Timing Array Collaboration, G. Agazie, J. Antoniadis, A. Anumarlapudi, A. M. Archibald, P. Arumugam, S. Arumugam, Z. Arzoumanian, J. Askew, S. Babak, M. Bagchi, M. Bailes, A. -S. Bak Nielsen, P. T. Baker, C. G. Bassa, A. Bathula, B. Bécsy, A. Berthereau, N. D. R. Bhat, L. Blecha, M. Bonetti, E. Bortolas, A. Brazier, P. R. Brook, M. Burgay , et al. (220 additional authors not shown)

Abstract: The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTA… ▽ More The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTAs that constitute the International Pulsar Timing Array (IPTA). We show that despite making different modeling choices, there is no significant difference in the GWB parameters that are measured by the different PTAs, agreeing within $1σ$. The pulsar noise parameters are also consistent between different PTAs for the majority of the pulsars included in these analyses. We bridge the differences in modeling choices by adopting a standardized noise model for all pulsars and PTAs, finding that under this model there is a reduction in the tension in the pulsar noise parameters. As part of this reanalysis, we "extended" each PTA's data set by adding extra pulsars that were not timed by that PTA. Under these extensions, we find better constraints on the GWB amplitude and a higher signal-to-noise ratio for the Hellings and Downs correlations. These extensions serve as a prelude to the benefits offered by a full combination of data across all pulsars in the IPTA, i.e., the IPTA's Data Release 3, which will involve not just adding in additional pulsars, but also including data from all three PTAs where any given pulsar is timed by more than as single PTA. △ Less

Submitted 1 September, 2023; originally announced September 2023.

Comments: 21 pages, 9 figures, submitted to ApJ

The second data release from the European Pulsar Timing Array: IV. Implications for massive black holes, dark matter and the early Universe

Authors: J. Antoniadis, P. Arumugam, S. Arumugam, P. Auclair, S. Babak, M. Bagchi, A. -S. Bak Nielsen, E. Barausse, C. G. Bassa, A. Bathula, A. Berthereau, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, C. Caprini, A. Chalumeau, D. J. Champion, S. Chanlaridis, S. Chen, I. Cognard, M. Crisostomi, S. Dandapat, D. Deb , et al. (89 additional authors not shown)

Abstract: The European Pulsar Timing Array (EPTA) and Indian Pulsar Timing Array (InPTA) collaborations have measured a low-frequency common signal in the combination of their second and first data releases respectively, with the correlation properties of a gravitational wave background (GWB). Such signal may have its origin in a number of physical processes including a cosmic population of inspiralling sup… ▽ More The European Pulsar Timing Array (EPTA) and Indian Pulsar Timing Array (InPTA) collaborations have measured a low-frequency common signal in the combination of their second and first data releases respectively, with the correlation properties of a gravitational wave background (GWB). Such signal may have its origin in a number of physical processes including a cosmic population of inspiralling supermassive black hole binaries (SMBHBs); inflation, phase transitions, cosmic strings and tensor mode generation by non-linear evolution of scalar perturbations in the early Universe; oscillations of the Galactic potential in the presence of ultra-light dark matter (ULDM). At the current stage of emerging evidence, it is impossible to discriminate among the different origins. Therefore, in this paper, we consider each process separately, and investigate the implications of the signal under the hypothesis that it is generated by that specific process. We find that the signal is consistent with a cosmic population of inspiralling SMBHBs, and its relatively high amplitude can be used to place constraints on binary merger timescales and the SMBH-host galaxy scaling relations. If this origin is confirmed, this is the first direct evidence that SMBHBs merge in nature, adding an important observational piece to the puzzle of structure formation and galaxy evolution. As for early Universe processes, the measurement would place tight constraints on the cosmic string tension and on the level of turbulence developed by first-order phase transitions. Other processes would require non-standard scenarios, such as a blue-tilted inflationary spectrum or an excess in the primordial spectrum of scalar perturbations at large wavenumbers. Finally, a ULDM origin of the detected signal is disfavoured, which leads to direct constraints on the abundance of ULDM in our Galaxy. △ Less

Submitted 15 May, 2024; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 30 pages, 23 figures, replaced to match the version published in Astronomy & Astrophysics, note the change in the numbering order in the series (now paper IV)

The second data release from the European Pulsar Timing Array V. Search for continuous gravitational wave signals

Authors: J. Antoniadis, P. Arumugam, S. Arumugam, S. Babak, M. Bagchi, A. S. Bak Nielsen, C. G. Bassa, A. Bathula, A. Berthereau, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, A. Chalumeau, D. J. Champion, S. Chanlaridis, S. Chen, I. Cognard, S. Dandapat, D. Deb, S. Desai, G. Desvignes, N. Dhanda-Batra, C. Dwivedi , et al. (75 additional authors not shown)

Abstract: We present the results of a search for continuous gravitational wave signals (CGWs) in the second data release (DR2) of the European Pulsar Timing Array (EPTA) collaboration. The most significant candidate event from this search has a gravitational wave frequency of 4-5 nHz. Such a signal could be generated by a supermassive black hole binary (SMBHB) in the local Universe. We present the results o… ▽ More We present the results of a search for continuous gravitational wave signals (CGWs) in the second data release (DR2) of the European Pulsar Timing Array (EPTA) collaboration. The most significant candidate event from this search has a gravitational wave frequency of 4-5 nHz. Such a signal could be generated by a supermassive black hole binary (SMBHB) in the local Universe. We present the results of a follow-up analysis of this candidate using both Bayesian and frequentist methods. The Bayesian analysis gives a Bayes factor of 4 in favor of the presence of the CGW over a common uncorrelated noise process, while the frequentist analysis estimates the p-value of the candidate to be 1%, also assuming the presence of common uncorrelated red noise. However, comparing a model that includes both a CGW and a gravitational wave background (GWB) to a GWB only, the Bayes factor in favour of the CGW model is only 0.7. Therefore, we cannot conclusively determine the origin of the observed feature, but we cannot rule it out as a CGW source. We present results of simulations that demonstrate that data containing a weak gravitational wave background can be misinterpreted as data including a CGW and vice versa, providing two plausible explanations of the EPTA DR2 data. Further investigations combining data from all PTA collaborations will be needed to reveal the true origin of this feature. △ Less

Submitted 25 June, 2024; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 13 figures, 15 pages, accepted

The second data release from the European Pulsar Timing Array II. Customised pulsar noise models for spatially correlated gravitational waves

Authors: J. Antoniadis, P. Arumugam, S. Arumugam, S. Babak, M. Bagchi, A. S. Bak Nielsen, C. G. Bassa, A. Bathula, A. Berthereau, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, A. Chalumeau, D. J. Champion, S. Chanlaridis, S. Chen, I. Cognard, S. Dandapat, D. Deb, S. Desai, G. Desvignes, N. Dhanda-Batra, C. Dwivedi , et al. (73 additional authors not shown)

Abstract: The nanohertz gravitational wave background (GWB) is expected to be an aggregate signal of an ensemble of gravitational waves emitted predominantly by a large population of coalescing supermassive black hole binaries in the centres of merging galaxies. Pulsar timing arrays, ensembles of extremely stable pulsars, are the most precise experiments capable of detecting this background. However, the su… ▽ More The nanohertz gravitational wave background (GWB) is expected to be an aggregate signal of an ensemble of gravitational waves emitted predominantly by a large population of coalescing supermassive black hole binaries in the centres of merging galaxies. Pulsar timing arrays, ensembles of extremely stable pulsars, are the most precise experiments capable of detecting this background. However, the subtle imprints that the GWB induces on pulsar timing data are obscured by many sources of noise. These must be carefully characterized to increase the sensitivity to the GWB. In this paper, we present a novel technique to estimate the optimal number of frequency coefficients for modelling achromatic and chromatic noise and perform model selection. We also incorporate a new model to fit for scattering variations in the pulsar timing package temponest and created realistic simulations of the European Pulsar Timing Array (EPTA) datasets that allowed us to test the efficacy of our noise modelling algorithms. We present an in-depth analysis of the noise properties of 25 millisecond pulsars (MSPs) that form the second data release (DR2) of the EPTA and investigate the effect of incorporating low-frequency data from the Indian PTA collaboration. We use enterprise and temponest packages to compare noise models with those reported with the EPTA DR1. We find that, while in some pulsars we can successfully disentangle chromatic from achromatic noise owing to the wider frequency coverage in DR2, in others the noise models evolve in a more complicated way. We also find evidence of long-term scattering variations in PSR J1600$-$3053. Through our simulations, we identify intrinsic biases in our current noise analysis techniques and discuss their effect on GWB searches. The results presented here directly help improve sensitivity to the GWB and are already being used as part of global PTA efforts. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 20 pages, 6 figures, 9 tables

Journal ref: A&A 678, A49 (2023)

The second data release from the European Pulsar Timing Array III. Search for gravitational wave signals

Authors: J. Antoniadis, P. Arumugam, S. Arumugam, S. Babak, M. Bagchi, A. -S. Bak Nielsen, C. G. Bassa, A. Bathula, A. Berthereau, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, A. Chalumeau, D. J. Champion, S. Chanlaridis, S. Chen, I. Cognard, S. Dandapat, D. Deb, S. Desai, G. Desvignes, N. Dhanda-Batra, C. Dwivedi , et al. (73 additional authors not shown)

Abstract: We present the results of the search for an isotropic stochastic gravitational wave background (GWB) at nanohertz frequencies using the second data release of the European Pulsar Timing Array (EPTA) for 25 millisecond pulsars and a combination with the first data release of the Indian Pulsar Timing Array (InPTA). We analysed (i) the full 24.7-year EPTA data set, (ii) its 10.3-year subset based on… ▽ More We present the results of the search for an isotropic stochastic gravitational wave background (GWB) at nanohertz frequencies using the second data release of the European Pulsar Timing Array (EPTA) for 25 millisecond pulsars and a combination with the first data release of the Indian Pulsar Timing Array (InPTA). We analysed (i) the full 24.7-year EPTA data set, (ii) its 10.3-year subset based on modern observing systems, (iii) the combination of the full data set with the first data release of the InPTA for ten commonly timed millisecond pulsars, and (iv) the combination of the 10.3-year subset with the InPTA data. These combinations allowed us to probe the contributions of instrumental noise and interstellar propagation effects. With the full data set, we find marginal evidence for a GWB, with a Bayes factor of four and a false alarm probability of $4\%$. With the 10.3-year subset, we report evidence for a GWB, with a Bayes factor of $60$ and a false alarm probability of about $0.1\%$ ($\gtrsim 3σ$ significance). The addition of the InPTA data yields results that are broadly consistent with the EPTA-only data sets, with the benefit of better noise modelling. Analyses were performed with different data processing pipelines to test the consistency of the results from independent software packages. The inferred spectrum from the latest EPTA data from new generation observing systems is rather uncertain and in mild tension with the common signal measured in the full data set. However, if the spectral index is fixed at 13/3, the two data sets give a similar amplitude of ($2.5\pm0.7)\times10^{-15}$ at a reference frequency of $1\,{\rm yr}^{-1}$. By continuing our detection efforts as part of the International Pulsar Timing Array (IPTA), we expect to be able to improve the measurement of spatial correlations and better characterise this signal in the coming years. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 21 pages, 14 figures, 4 appendix figures, accepted for publication in A&A

Journal ref: A&A 678, A50 (2023)

Antennas for low-frequency radio telescope of SKA

Authors: Agaram Raghunathan, Keerthipriya Satish, Arasi Sathyamurthy, T. Prabu, B. S. Girish, K. S. Srivani, Shiv K. Sethi

Abstract: The low-frequency radio telescope of the Square Kilometre Array (SKA) is being built by the international radio astronomical community to (i) have orders of magnitude higher sensitivity and (ii) be able to map the sky several hundred times faster, than any other existing facilities over the frequency range of 50 - 350 MHz. The sensitivity of a radio telescope array is in general, dependent upon th… ▽ More The low-frequency radio telescope of the Square Kilometre Array (SKA) is being built by the international radio astronomical community to (i) have orders of magnitude higher sensitivity and (ii) be able to map the sky several hundred times faster, than any other existing facilities over the frequency range of 50 - 350 MHz. The sensitivity of a radio telescope array is in general, dependent upon the number of electromagnetic sensors used to receive the sky signal. The total number of them is further constrained by the effects of mutual coupling between the sensor elements, allowable grating lobes in their radiation patterns, etc. The operating frequency band is governed by the desired spatial and spectral responses, acceptable sidelobe and backlobe levels, radiation efficiency, polarization purity and calibratability of sensors' response. This paper presents a brief review of several broadband antennas considered as potential candidates by various engineering groups across the globe, for the low-frequency radio telescope of SKA covering the frequency range of 50 - 350 MHz, on the basis of their suitability for conducting primary scientific objectives. △ Less

Submitted 26 May, 2023; v1 submitted 17 January, 2023; originally announced January 2023.

Comments: 14 pages, 33 figures, JoAA - Special issue on the SKA (2023) - Accepted for publication

Investigation of a Machine learning methodology for the SKA pulsar search pipeline

Authors: Shashank Sanjay Bhat, Thiagaraj Prabu, Ben Stappers, Atul Ghalame, Snehanshu Saha, T. S. B Sudarshan, Zafiirah Hosenie

Abstract: The SKA pulsar search pipeline will be used for real time detection of pulsars. Modern radio telescopes such as SKA will be generating petabytes of data in their full scale of operation. Hence experience-based and data-driven algorithms become indispensable for applications such as candidate detection. Here we describe our findings from testing a state of the art object detection algorithm called… ▽ More The SKA pulsar search pipeline will be used for real time detection of pulsars. Modern radio telescopes such as SKA will be generating petabytes of data in their full scale of operation. Hence experience-based and data-driven algorithms become indispensable for applications such as candidate detection. Here we describe our findings from testing a state of the art object detection algorithm called Mask R-CNN to detect candidate signatures in the SKA pulsar search pipeline. We have trained the Mask R-CNN model to detect candidate images. A custom annotation tool was developed to mark the regions of interest in large datasets efficiently. We have successfully demonstrated this algorithm by detecting candidate signatures on a simulation dataset. The paper presents details of this work with a highlight on the future prospects. △ Less

Submitted 17 January, 2023; v1 submitted 9 September, 2022; originally announced September 2022.

Journal ref: Journal of Astronomy and Astrophysics SKA special issue 2022-23 (Under review)

High-Performance Computing for SKA Transient Search: Use of FPGA based Accelerators -- a brief review

Authors: R. Aafreen, R. Abhishek, B. Ajithkumar, Arunkumar M. Vaidyanathan, Indrajit. V. Barve, Sahana Bhattramakki, Shashank Bhat, B. S. Girish, Atul Ghalame, Y. Gupta, Harshal G. Hayatnagarkar, P. A. Kamini, A. Karastergiou, L. Levin, S. Madhavi, M. Mekhala, M. Mickaliger, V. Mugundhan, Arun Naidu, J. Oppermann, B. Arul Pandian, N. Patra, A. Raghunathan, Jayanta Roy, Shiv Sethi , et al. (12 additional authors not shown)

Abstract: This paper presents the High-Performance computing efforts with FPGA for the accelerated pulsar/transient search for the SKA. Case studies are presented from within SKA and pathfinder telescopes highlighting future opportunities. It reviews the scenario that has shifted from offline processing of the radio telescope data to digitizing several hundreds/thousands of antenna outputs over huge bandwid… ▽ More This paper presents the High-Performance computing efforts with FPGA for the accelerated pulsar/transient search for the SKA. Case studies are presented from within SKA and pathfinder telescopes highlighting future opportunities. It reviews the scenario that has shifted from offline processing of the radio telescope data to digitizing several hundreds/thousands of antenna outputs over huge bandwidths, forming several 100s of beams, and processing the data in the SKA real-time pulsar search pipelines. A brief account of the different architectures of the accelerators, primarily the new generation Field Programmable Gate Array-based accelerators, showing their critical roles to achieve high-performance computing and in handling the enormous data volume problems of the SKA is presented here. It also presents the power-performance efficiency of this emerging technology and presents potential future scenarios. △ Less

Submitted 17 January, 2023; v1 submitted 14 July, 2022; originally announced July 2022.

Comments: Accepted for publication in the JoAA special issue on SKA (2023)

Nanohertz Gravitational Wave Astronomy during the SKA Era: An InPTA perspective

Authors: Bhal Chandra Joshi, Achamveedu Gopakumar, Arul Pandian, Thiagaraj Prabu, Lankeswar Dey, Manjari Bagchi, Shantanu Desai, Pratik Tarafdar, Prerna Rana, Yogesh Maan, Neelam Dhanda Batra, Raghav Girgaonkar, Nikita Agarwal, Paramasivan Arumugam, Sarmistha Banik, Avishek Basu, Adarsh Bathula, Subhajit Dandapat, Yashwant Gupta, Shinnosuke Hisano, Ryo Kato, Divyansh Kharbanda, Tomonosuke Kikunaga, Neel Kolhe, M. A. Krishnakumar , et al. (12 additional authors not shown)

Abstract: Decades long monitoring of millisecond pulsars, which exhibit highly stable rotational periods, in pulsar timing array experiments is on the threshold of discovering nanohertz stochastic gravitational wave background. This paper describes the Indian Pulsar timing array (InPTA) experiment, which employs the upgraded Giant Metrewave Radio Telescope (uGMRT) for timing an ensemble of millisecond pulsa… ▽ More Decades long monitoring of millisecond pulsars, which exhibit highly stable rotational periods, in pulsar timing array experiments is on the threshold of discovering nanohertz stochastic gravitational wave background. This paper describes the Indian Pulsar timing array (InPTA) experiment, which employs the upgraded Giant Metrewave Radio Telescope (uGMRT) for timing an ensemble of millisecond pulsars for this purpose. We highlight InPTA's observation strategies and analysis methods, which are relevant for a future PTA experiment with the more sensitive Square Kilometer Array (SKA) telescope. We show that the unique multi-sub-array multi-band wide-bandwidth frequency coverage of the InPTA provides Dispersion Measure estimates with unprecedented precision for PTA pulsars, e.g., ~ 2 x 10{-5} pc-cm{-3} for PSR J1909-3744. Configuring the SKA-low and SKA-mid as two and four sub-arrays respectively, it is shown that comparable precision is achievable, using observation strategies similar to those pursued by the InPTA, for a larger sample of 62 pulsars requiring about 26 and 7 hours per epoch for the SKA-mid and the SKA-low telescopes respectively. We also review the ongoing efforts to develop PTA-relevant general relativistic constructs that will be required to search for nanohertz gravitational waves from isolated super-massive black hole binary systems like blazar OJ 287. These efforts should be relevant to pursue persistent multi-messenger gravitational wave astronomy during the forthcoming era of the SKA telescope, the Thirty Meter Telescope, and the next-generation Event Horizon Telescope. △ Less

Submitted 13 July, 2022; originally announced July 2022.

Comments: Accepted for publication in Journal of Astronomy and Astrophysics for Special Issue on Indian Participation in the SKA (Editors : Abhirup Datta, Nirupam Roy, Preeti Kharb and Tirthankar Roy Choudhury)

The Indian Pulsar Timing Array: First data release

Authors: Pratik Tarafdar, Nobleson K., Prerna Rana, Jaikhomba Singha, M. A. Krishnakumar, Bhal Chandra Joshi, Avinash Kumar Paladi, Neel Kolhe, Neelam Dhanda Batra, Nikita Agarwal, Adarsh Bathula, Subhajit Dandapat, Shantanu Desai, Lankeswar Dey, Shinnosuke Hisano, Prathamesh Ingale, Ryo Kato, Divyansh Kharbanda, Tomonosuke Kikunaga, Piyush Marmat, B. Arul Pandian, T. Prabu, Aman Srivastava, Mayuresh Surnis, Sai Chaitanya Susarla , et al. (13 additional authors not shown)

Abstract: We present the pulse arrival times and high-precision dispersion measure estimates for 14 millisecond pulsars observed simultaneously in the 300-500 MHz and 1260-1460 MHz frequency bands using the upgraded Giant Metrewave Radio Telescope (uGMRT). The data spans over a baseline of 3.5 years (2018-2021), and is the first official data release made available by the Indian Pulsar Timing Array collabor… ▽ More We present the pulse arrival times and high-precision dispersion measure estimates for 14 millisecond pulsars observed simultaneously in the 300-500 MHz and 1260-1460 MHz frequency bands using the upgraded Giant Metrewave Radio Telescope (uGMRT). The data spans over a baseline of 3.5 years (2018-2021), and is the first official data release made available by the Indian Pulsar Timing Array collaboration. This data release presents a unique opportunity for investigating the interstellar medium effects at low radio frequencies and their impact on the timing precision of pulsar timing array experiments. In addition to the dispersion measure time series and pulse arrival times obtained using both narrowband and wideband timing techniques, we also present the dispersion measure structure function analysis for selected pulsars. Our ongoing investigations regarding the frequency dependence of dispersion measures have been discussed. Based on the preliminary analysis for five millisecond pulsars, we do not find any conclusive evidence of chromaticity in dispersion measures. Data from regular simultaneous two-frequency observations are presented for the first time in this work. This distinctive feature leads us to the highest precision dispersion measure estimates obtained so far for a subset of our sample. Simultaneous multi-band uGMRT observations in Band 3 and Band 5 are crucial for high-precision dispersion measure estimation and for the prospect of expanding the overall frequency coverage upon the combination of data from the various Pulsar Timing Array consortia in the near future. Parts of the data presented in this work are expected to be incorporated into the upcoming third data release of the International Pulsar Timing Array. △ Less

Submitted 25 October, 2022; v1 submitted 18 June, 2022; originally announced June 2022.

Comments: 23 pages, 21 figures, 3 tables. Published in PASA

Journal ref: Publications of the Astronomical Society of Australia, Volume 39, 2022, e053

Low-frequency wideband timing of InPTA pulsars observed with the uGMRT

Authors: K Nobleson, Nikita Agarwal, Raghav Girgaonkar, Arul Pandian, Bhal Chandra Joshi, M A Krishnakumar, Abhimanyu Susobhanan, Shantanu Desai, T Prabu, Adarsh Bathula, Timothy T Pennucci, Sarmistha Banik, Manjari Bagchi, Neelam Dhanda Batra, Arpita Choudhary, Subhajit Dandapat, Lankeswar Dey, Yashwant Gupta, Shinnosuke Hisano, Ryo Kato, Divyansh Kharbanda, Tomonosuke Kikunaga, Neel Kolhe, Yogesh Maan, Piyush Marmat , et al. (7 additional authors not shown)

Abstract: High-precision measurements of the pulsar dispersion measure (DM) are possible using telescopes with low-frequency wideband receivers. We present an initial study of the application of the wideband timing technique, which can simultaneously measure the pulsar times of arrival (ToAs) and DMs, for a set of five pulsars observed with the upgraded Giant Metrewave Radio Telescope (uGMRT) as part of the… ▽ More High-precision measurements of the pulsar dispersion measure (DM) are possible using telescopes with low-frequency wideband receivers. We present an initial study of the application of the wideband timing technique, which can simultaneously measure the pulsar times of arrival (ToAs) and DMs, for a set of five pulsars observed with the upgraded Giant Metrewave Radio Telescope (uGMRT) as part of the Indian Pulsar Timing Array (InPTA) campaign. We have used the observations with the 300-500 MHz band of the uGMRT for this purpose. We obtain high precision in DM measurements with precisions of the order 10^{-6}cm^{-3}pc. The ToAs obtained have sub-μs precision and the root-mean-square of the post-fit ToA residuals are in the sub-μs range. We find that the uncertainties in the DMs and ToAs obtained with this wideband technique, applied to low-frequency data, are consistent with the results obtained with traditional pulsar timing techniques and comparable to high-frequency results from other PTAs. This work opens up an interesting possibility of using low-frequency wideband observations for precision pulsar timing and gravitational wave detection with similar precision as high-frequency observations used conventionally. △ Less

Submitted 23 February, 2022; v1 submitted 13 December, 2021; originally announced December 2021.

Comments: Accepted for publication in MNRAS

Evidence for profile changes in PSR J1713+0747 using the uGMRT

Authors: Jaikhomba Singha, Mayuresh P Surnis, Bhal Chandra Joshi, Pratik Tarafdar, Prerna Rana, Abhimanyu Susobhanan, Raghav Girgaonkar, Neel Kolhe, Nikita Agarwal, Shantanu Desai, T Prabu, Adarsh Bathula, Subhajit Dandapat, Lankeswar Dey, Shinnosuke Hisano, Ryo Kato, Divyansh Kharbanda, Tomonosuke Kikunaga, Piyush Marmat, Sai Chaitanya Susarla, Manjari Bagchi, Neelam Dhanda Batra, Arpita Choudhury, A Gopakumar, Yashwant Gupta , et al. (7 additional authors not shown)

Abstract: PSR J1713+0747 is one of the most precisely timed pulsars in the international pulsar timing array experiment. This pulsar showed an abrupt profile shape change between April 16, 2021 (MJD 59320) and April 17, 2021 (MJD 59321). In this paper, we report the results from multi-frequency observations of this pulsar carried out with the upgraded Giant Metrewave Radio Telescope (uGMRT) before and after… ▽ More PSR J1713+0747 is one of the most precisely timed pulsars in the international pulsar timing array experiment. This pulsar showed an abrupt profile shape change between April 16, 2021 (MJD 59320) and April 17, 2021 (MJD 59321). In this paper, we report the results from multi-frequency observations of this pulsar carried out with the upgraded Giant Metrewave Radio Telescope (uGMRT) before and after the event. We demonstrate the profile change seen in Band 5 (1260 MHz - 1460 MHz) and Band 3 (300 MHz - 500 MHz). The timing analysis of this pulsar shows a disturbance accompanying this profile change followed by a recovery with a timescale of $\sim 159$ days. Our data suggest that a model with chromatic index as a free parameter is preferred over models with combinations of achromaticity with DM bump or scattering bump. We determine the frequency dependence to be $\simν^{+1.34}$. △ Less

Submitted 16 August, 2021; v1 submitted 9 July, 2021; originally announced July 2021.

Comments: Accepted for publication in MNRAS-Letters

High Precision Measurements of Interstellar Dispersion Measure with the upgraded GMRT

Authors: M. A. Krishnakumar, P. K. Manoharan, Bhal Chandra Joshi, Raghav Girgaonkar, Shantanu Desai, Manjari Bagchi, K. Nobleson, Lankeswar Dey, Abhimanyu Susobhanan, Sai Chaitanya Susarla, Mayuresh P. Surnis, Yogesh Maan, A. Gopakumar, Avishek Basu, Neelam Dhanda Batra, Arpita Choudhary, Kishalay De, Yashwant Gupta, Arun Kumar Naidu, Dhruv Pathak, Jaikhomba Singha, T. Prabu

Abstract: Pulsar radio emission undergoes dispersion due to the presence of free electrons in the interstellar medium (ISM). The dispersive delay in the arrival time of pulsar signal changes over time due to the varying ISM electron column density along the line of sight. Correcting for this delay accurately is crucial for the detection of nanohertz gravitational waves using Pulsar Timing Arrays. In this wo… ▽ More Pulsar radio emission undergoes dispersion due to the presence of free electrons in the interstellar medium (ISM). The dispersive delay in the arrival time of pulsar signal changes over time due to the varying ISM electron column density along the line of sight. Correcting for this delay accurately is crucial for the detection of nanohertz gravitational waves using Pulsar Timing Arrays. In this work, we present in-band and inter-band DM estimates of four pulsars observed with uGMRT over the timescale of a year using two different template alignment methods. The DMs obtained using both these methods show only subtle differences for PSR 1713+0747 and J1909$-$3744. A considerable offset is seen in the DM of PSR J1939+2134 and J2145$-$0750 between the two methods. This could be due to the presence of scattering in the former and profile evolution in the latter. We find that both methods are useful but could have a systematic offset between the DMs obtained. Irrespective of the template alignment methods followed, the precision on the DMs obtained is about $10^{-3}$ pc cm$^{-3}$ using only BAND3 and $10^{-4}$ pc cm$^{-3}$ after combining data from BAND3 and BAND5 of the uGMRT. In a particular result, we have detected a DM excess of about $5\times10^{-3}$ pc cm$^{-3}$ on 24 February 2019 for PSR J2145$-$0750. This excess appears to be due to the interaction region created by fast solar wind from a coronal hole and a coronal mass ejection (CME) observed from the Sun on that epoch. A detailed analysis of this interesting event is presented. △ Less

Submitted 1 July, 2021; v1 submitted 13 January, 2021; originally announced January 2021.

Comments: 11 pages, 6 figures, 2 tables. Accepted by A&A

Journal ref: A&A 651, A5 (2021)

pinta: The uGMRT Data Processing Pipeline for the Indian Pulsar Timing Array

Authors: Abhimanyu Susobhanan, Yogesh Maan, Bhal Chandra Joshi, T. Prabu, Shantanu Desai, K. Nobleson, Sai Chaitanya Susarla, Raghav Girgaonkar, Lankeswar Dey, Neelam Dhanda Batra, Yashwant Gupta, A. Gopakumar, Manjari Bagchi, Avishek Basu, Suryarao Bethapudi, Arpita Choudhary, Kishalay De, M. A. Krishnakumar, P. K. Manoharan, Arun Kumar Naidu, Dhruv Pathak, Jaikhomba Singha, Mayuresh P. Surnis

Abstract: We introduce pinta, a pipeline for reducing the upgraded Giant Metre-wave Radio Telescope (uGMRT) raw pulsar timing data, developed for the Indian Pulsar Timing Array experiment. We provide a detailed description of the workflow and usage of pinta, as well as its computational performance and RFI mitigation characteristics. We also discuss a novel and independent determination of the relative time… ▽ More We introduce pinta, a pipeline for reducing the upgraded Giant Metre-wave Radio Telescope (uGMRT) raw pulsar timing data, developed for the Indian Pulsar Timing Array experiment. We provide a detailed description of the workflow and usage of pinta, as well as its computational performance and RFI mitigation characteristics. We also discuss a novel and independent determination of the relative time offsets between the different back-end modes of uGMRT and the interpretation of the uGMRT observation frequency settings, and their agreement with results obtained from engineering tests. Further, we demonstrate the capability of pinta to generate data products which can produce high-precision TOAs using PSR J1909-3744 as an example. These results are crucial for performing precision pulsar timing with the uGMRT. △ Less

Submitted 16 April, 2021; v1 submitted 6 July, 2020; originally announced July 2020.

Journal ref: Publications of the Astronomical Society of Australia, 38, E017 (2021)

The Murchison Widefield Array Transients Survey (MWATS). A search for low frequency variability in a bright Southern hemisphere sample

Authors: M. E. Bell, Tara Murphy, P. J. Hancock, J. R. Callingham, S. Johnston, D. L. Kaplan, R. W. Hunstead, E. M. Sadler, S. Croft, S. V. White, N. Hurley-Walker, R. Chhetri, J. S. Morgan, P. G. Edwards, A. Rowlinson, A. R. Offringa, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, A. A. Deshpande, B. M. Gaensler, L. J. Greenhill, B. J. Hazelton, M. Johnston-Hollitt , et al. (16 additional authors not shown)

Abstract: We report on a search for low-frequency radio variability in 944 bright (> 4Jy at 154 MHz) unresolved, extragalactic radio sources monitored monthly for several years with the Murchison Widefield Array. In the majority of sources we find very low levels of variability with typical modulation indices < 5%. We detect 15 candidate low frequency variables that show significant long term variability (>… ▽ More We report on a search for low-frequency radio variability in 944 bright (> 4Jy at 154 MHz) unresolved, extragalactic radio sources monitored monthly for several years with the Murchison Widefield Array. In the majority of sources we find very low levels of variability with typical modulation indices < 5%. We detect 15 candidate low frequency variables that show significant long term variability (>2.8 years) with time-averaged modulation indices M = 3.1 - 7.1%. With 7/15 of these variable sources having peaked spectral energy distributions, and only 5.7% of the overall sample having peaked spectra, we find an increase in the prevalence of variability in this spectral class. We conclude that the variability seen in this survey is most probably a consequence of refractive interstellar scintillation and that these objects must have the majority of their flux density contained within angular diameters less than 50 milli-arcsec (which we support with multi-wavelength data). At 154 MHz we demonstrate that interstellar scintillation time-scales become long (~decades) and have low modulation indices, whilst synchrotron driven variability can only produce dynamic changes on time-scales of hundreds of years, with flux density changes less than one milli-jansky (without relativistic boosting). From this work we infer that the low frequency extra-galactic southern sky, as seen by SKA-Low, will be non-variable on time-scales shorter than one year. △ Less

Submitted 23 October, 2018; originally announced October 2018.

Comments: 19 pages, 11 figures

Pulsar Searches with the SKA

Authors: L. Levin, W. Armour, C. Baffa, E. Barr, S. Cooper, R. Eatough, A. Ensor, E. Giani, A. Karastergiou, R. Karuppusamy, M. Keith, M. Kramer, R. Lyon, M. Mackintosh, M. Mickaliger, R van Nieuwpoort, M. Pearson, T. Prabu, J. Roy, O. Sinnen, L. Spitler, H. Spreeuw, B. W. Stappers, W. van Straten, C. Williams , et al. (2 additional authors not shown)

Abstract: The Square Kilometre Array will be an amazing instrument for pulsar astronomy. While the full SKA will be sensitive enough to detect all pulsars in the Galaxy visible from Earth, already with SKA1, pulsar searches will discover enough pulsars to increase the currently known population by a factor of four, no doubt including a range of amazing unknown sources. Real time processing is needed to deal… ▽ More The Square Kilometre Array will be an amazing instrument for pulsar astronomy. While the full SKA will be sensitive enough to detect all pulsars in the Galaxy visible from Earth, already with SKA1, pulsar searches will discover enough pulsars to increase the currently known population by a factor of four, no doubt including a range of amazing unknown sources. Real time processing is needed to deal with the 60 PB of pulsar search data collected per day, using a signal processing pipeline required to perform more than 10 POps. Here we present the suggested design of the pulsar search engine for the SKA and discuss challenges and solutions to the pulsar search venture. △ Less

Submitted 4 December, 2017; originally announced December 2017.

Comments: 4 pages, 1 figure. To be published in Proceedings of IAU Symposium 337: Pulsar Astrophysics - The Next 50 Years

Spectral energy distribution and radio halo of NGC 253 at low radio frequencies

Authors: A. D. Kapinska, L. Staveley-Smith, R. Crocker, G. R. Meurer, S. Bhandari, N. Hurley-Walker, A. R. Offringa, D. J. Hanish, N. Seymour, R. D. Ekers, M. E. Bell, J. R. Callingham, K. S. Dwarakanath, B. -Q. For, B. M. Gaensler, P. J. Hancock, L. Hindson, M. Johnston-Hollitt, E. Lenc, B. McKinley, J. Morgan, P. Procopio, R. B. Wayth, C. Wu, Q. Zheng , et al. (45 additional authors not shown)

Abstract: We present new radio continuum observations of NGC253 from the Murchison Widefield Array at frequencies between 76 and 227 MHz. We model the broadband radio spectral energy distribution for the total flux density of NGC253 between 76 MHz and 11 GHz. The spectrum is best described as a sum of central starburst and extended emission. The central component, corresponding to the inner 500pc of the sta… ▽ More We present new radio continuum observations of NGC253 from the Murchison Widefield Array at frequencies between 76 and 227 MHz. We model the broadband radio spectral energy distribution for the total flux density of NGC253 between 76 MHz and 11 GHz. The spectrum is best described as a sum of central starburst and extended emission. The central component, corresponding to the inner 500pc of the starburst region of the galaxy, is best modelled as an internally free-free absorbed synchrotron plasma, with a turnover frequency around 230 MHz. The extended emission component of the NGC253 spectrum is best described as a synchrotron emission flattening at low radio frequencies. We find that 34% of the extended emission (outside the central starburst region) at 1 GHz becomes partially absorbed at low radio frequencies. Most of this flattening occurs in the western region of the SE halo, and may be indicative of synchrotron self-absorption of shock re-accelerated electrons or an intrinsic low-energy cut off of the electron distribution. Furthermore, we detect the large-scale synchrotron radio halo of NGC253 in our radio images. At 154 - 231 MHz the halo displays the well known X-shaped/horn-like structure, and extends out to ~8kpc in z-direction (from major axis). △ Less

Submitted 19 February, 2017; v1 submitted 8 February, 2017; originally announced February 2017.

Comments: Accepted for publication in Astrophysical Journal on 06 February 2017

A Matched Filter Technique For Slow Radio Transient Detection And First Demonstration With The Murchison Widefield Array

Authors: L. Feng, R. Vaulin, J. N. Hewitt, R. Remillard, D. L. Kaplan, Tara Murphy, N. Kudryavtseva, P. Hancock, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, A. A. Deshpande, B. M. Gaensler, L. J. Greenhill, B. J. Hazelton, M. Johnston-Hollitt, C. J. Lonsdale, S. R. McWhirter, D. A. Mitchell, M. F. Morales, E. Morgan, D. Oberoi, S. M. Ord, T. Prabu , et al. (8 additional authors not shown)

Abstract: Many astronomical sources produce transient phenomena at radio frequencies, but the transient sky at low frequencies (<300 MHz) remains relatively unexplored. Blind surveys with new widefield radio instruments are setting increasingly stringent limits on the transient surface density on various timescales. Although many of these instruments are limited by classical confusion noise from an ensemble… ▽ More Many astronomical sources produce transient phenomena at radio frequencies, but the transient sky at low frequencies (<300 MHz) remains relatively unexplored. Blind surveys with new widefield radio instruments are setting increasingly stringent limits on the transient surface density on various timescales. Although many of these instruments are limited by classical confusion noise from an ensemble of faint, unresolved sources, one can in principle detect transients below the classical confusion limit to the extent that the classical confusion noise is independent of time. We develop a technique for detecting radio transients that is based on temporal matched filters applied directly to time series of images rather than relying on source-finding algorithms applied to individual images. This technique has well-defined statistical properties and is applicable to variable and transient searches for both confusion-limited and non-confusion-limited instruments. Using the Murchison Widefield Array as an example, we demonstrate that the technique works well on real data despite the presence of classical confusion noise, sidelobe confusion noise, and other systematic errors. We searched for transients lasting between 2 minutes and 3 months. We found no transients and set improved upper limits on the transient surface density at 182 MHz for flux densities between ~20--200 mJy, providing the best limits to date for hour- and month-long transients. △ Less

Submitted 12 January, 2017; originally announced January 2017.

Comments: 16 pages, 14 figures, accepted to AJ

Wavelet-Based Characterization of Small-Scale Solar Emission Features at Low Radio Frequencies

Authors: Akshay Suresh, Rohit Sharma, Divya Oberoi, Srijan B. Das, Victor Pankratius, Brian Timar, Colin J. Lonsdale, Judd D. Bowman, Frank Briggs, Roger J. Cappallo, Brian E. Corey, Avinash A. Deshpande, David Emrich, Robert Goeke, Lincoln J. Greenhill, Bryna J. Hazelton, Melanie Johnston-Hollitt, David L. Kaplan, Justin C. Kasper, Eric Kratzenberg, Mervyn J. Lynch, S. Russell McWhirter, Daniel A. Mitchell, Miguel F. Morales, Edward Morgan , et al. (14 additional authors not shown)

Abstract: Low radio frequency solar observations using the Murchison Widefield Array have recently revealed the presence of numerous weak, short-lived and narrow-band emission features, even during moderately quiet solar conditions. These non-thermal features occur at rates of many thousands per hour in the 30.72 MHz observing bandwidth, and hence, necessarily require an automated approach for their detecti… ▽ More Low radio frequency solar observations using the Murchison Widefield Array have recently revealed the presence of numerous weak, short-lived and narrow-band emission features, even during moderately quiet solar conditions. These non-thermal features occur at rates of many thousands per hour in the 30.72 MHz observing bandwidth, and hence, necessarily require an automated approach for their detection and characterization. Here, we employ continuous wavelet transform using a mother Ricker wavelet for feature detection from the dynamic spectrum. We establish the efficacy of this approach and present the first statistically robust characterization of the properties of these features. In particular, we examine distributions of their peak flux densities, spectral spans, temporal spans and peak frequencies. We can reliably detect features weaker than 1 SFU, making them, to the best of our knowledge, the weakest bursts reported in literature. The distribution of their peak flux densities follows a power law with an index of -2.23 in the 12-155 SFU range, implying that they can provide an energetically significant contribution to coronal and chromospheric heating. These features typically last for 1-2 seconds and possess bandwidths of about 4-5 MHz. Their occurrence rate remains fairly flat in the 140-210 MHz frequency range. At the time resolution of the data, they appear as stationary bursts, exhibiting no perceptible frequency drift. These features also appear to ride on a broadband background continuum, hinting at the likelihood of them being weak type-I bursts. △ Less

Submitted 2 June, 2017; v1 submitted 3 December, 2016; originally announced December 2016.

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

GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey I: A low-frequency extragalactic catalogue

Authors: Natasha Hurley-Walker, Joseph R. Callingham, Paul J. Hancock, Thomas M. O. Franzen, Luke Hindson, Anna D. Kapinska, John Morgan, Andre R. Offringa, Randall B. Wayth, Chen Wu, Q. Zheng, Tara Murphy, Martin E. Bell, K. S. Dwarakanath, Bi-Qing For, Bryan M. Gaensler, Melanie Johnston-Hollitt, Emil Lenc, Pietro Procopio, Lister Staveley-Smith, Ron Ekers, Judd D. Bowman, Frank Briggs, R. J. Cappallo, Avinash A. Deshpande , et al. (18 additional authors not shown)

Abstract: Using the Murchison Widefield Array (MWA), the low-frequency Square Kilometre Array (SKA1 LOW) precursor located in Western Australia, we have completed the GaLactic and Extragalactic All-sky MWA (GLEAM) survey, and present the resulting extragalactic catalogue, utilising the first year of observations. The catalogue covers 24,831 square degrees, over declinations south of $+30^\circ$ and Galactic… ▽ More Using the Murchison Widefield Array (MWA), the low-frequency Square Kilometre Array (SKA1 LOW) precursor located in Western Australia, we have completed the GaLactic and Extragalactic All-sky MWA (GLEAM) survey, and present the resulting extragalactic catalogue, utilising the first year of observations. The catalogue covers 24,831 square degrees, over declinations south of $+30^\circ$ and Galactic latitudes outside $10^\circ$ of the Galactic plane, excluding some areas such as the Magellanic Clouds. It contains 307,455 radio sources with 20 separate flux density measurements across 72--231MHz, selected from a time- and frequency- integrated image centred at 200MHz, with a resolution of $\approx 2$'. Over the catalogued region, we estimate that the catalogue is 90% complete at 170mJy, and 50% complete at 55mJy, and large areas are complete at even lower flux density levels. Its reliability is 99.97% above the detection threshold of $5σ$, which itself is typically 50mJy. These observations constitute the widest fractional bandwidth and largest sky area survey at radio frequencies to date, and calibrate the low frequency flux density scale of the southern sky to better than 10%. This paper presents details of the flagging, imaging, mosaicking, and source extraction/characterisation, as well as estimates of the completeness and reliability. All source measurements and images are available online (http://www.mwatelescope.org/science/gleam-survey). This is the first in a series of publications describing the GLEAM survey results. △ Less

Submitted 26 October, 2016; originally announced October 2016.

Comments: 30 pages, 18 figures, 6 tables, published in Monthly Notices of the Royal Astronomical Society

Delay Spectrum with Phase-Tracking Arrays: Extracting the HI power spectrum from the Epoch of Reionization

Authors: Sourabh Paul, Shiv K. Sethi, Miguel F. Morales, K. S. Dwarkanath, N. Udaya Shankar, Ravi Subrahmanyan, N. Barry, A. P. Beardsley, Judd D. Bowman, F. Briggs, P. Carroll, A. de Oliveira-Costa, Joshua S. Dillon, A. Ewall-Wice, L. Feng, L. J. Greenhill, B. M. Gaensler, B. J. Hazelton, J. N. Hewitt, N. Hurley-Walker, D. J. Jacobs, Han-Seek Kim, P. Kittiwisit, E. Lenc, J. Line , et al. (29 additional authors not shown)

Abstract: The Detection of redshifted 21 cm emission from the epoch of reionization (EoR) is a challenging task owing to strong foregrounds that dominate the signal. In this paper, we propose a general method, based on the delay spectrum approach, to extract HI power spectra that is applicable to tracking observations using an imaging radio interferometer (Delay Spectrum with Imaging Arrays (DSIA)). Our met… ▽ More The Detection of redshifted 21 cm emission from the epoch of reionization (EoR) is a challenging task owing to strong foregrounds that dominate the signal. In this paper, we propose a general method, based on the delay spectrum approach, to extract HI power spectra that is applicable to tracking observations using an imaging radio interferometer (Delay Spectrum with Imaging Arrays (DSIA)). Our method is based on modelling the HI signal taking into account the impact of wide field effects such as the $w$-term which are then used as appropriate weights in cross-correlating the measured visibilities. Our method is applicable to any radio interferometer that tracks a phase center and could be utilized for arrays such as MWA, LOFAR, GMRT, PAPER and HERA. In the literature the delay spectrum approach has been implemented for near-redundant baselines using drift scan observations. In this paper we explore the scheme for non-redundant tracking arrays, and this is the first application of delay spectrum methodology to such data to extract the HI signal. We analyze 3 hours of MWA tracking data on the EoR1 field. We present both 2-dimensional ($k_\parallel,k_\perp$) and 1-dimensional (k) power spectra from the analysis. Our results are in agreement with the findings of other pipelines developed to analyse the MWA EoR data. △ Less

Submitted 22 October, 2016; originally announced October 2016.

Comments: 17 pages, 10 figures, accepted for publication in The Astrophysical Journal

Journal ref: ApJ, 2016, 833, 2, 213

First Season MWA EoR Power Spectrum Results at Redshift 7

Authors: A. P. Beardsley, B. J. Hazelton, I. S. Sullivan, P. Carroll, N. Barry, M. Rahimi, B. Pindor, C. M. Trott, J. Line, Daniel C. Jacobs, M. F. Morales, J. C. Pober, G. Bernardi, Judd D. Bowman, M. P. Busch, F. Briggs, R. J. Cappallo, B. E. Corey, A. de Oliveira-Costa, Joshua S. Dillon, D. Emrich, A. Ewall-Wice, L. Feng, B. M. Gaensler, R. Goeke , et al. (41 additional authors not shown)

Abstract: The Murchison Widefield Array (MWA) has collected hundreds of hours of Epoch of Reionization (EoR) data and now faces the challenge of overcoming foreground and systematic contamination to reduce the data to a cosmological measurement. We introduce several novel analysis techniques such as cable reflection calibration, hyper-resolution gridding kernels, diffuse foreground model subtraction, and qu… ▽ More The Murchison Widefield Array (MWA) has collected hundreds of hours of Epoch of Reionization (EoR) data and now faces the challenge of overcoming foreground and systematic contamination to reduce the data to a cosmological measurement. We introduce several novel analysis techniques such as cable reflection calibration, hyper-resolution gridding kernels, diffuse foreground model subtraction, and quality control methods. Each change to the analysis pipeline is tested against a two dimensional power spectrum figure of merit to demonstrate improvement. We incorporate the new techniques into a deep integration of 32 hours of MWA data. This data set is used to place a systematic-limited upper limit on the cosmological power spectrum of $Δ^2 \leq 2.7 \times 10^4$ mK$^2$ at $k=0.27$ h~Mpc$^{-1}$ and $z=7.1$, consistent with other published limits, and a modest improvement (factor of 1.4) over previous MWA results. From this deep analysis we have identified a list of improvements to be made to our EoR data analysis strategies. These improvements will be implemented in the future and detailed in upcoming publications. △ Less

Submitted 22 August, 2016; originally announced August 2016.

Comments: 21 pages, 16 figures, accepted to ApJ

Low frequency observations of linearly polarized structures in the interstellar medium near the south Galactic pole

Authors: Emil Lenc, B. M. Gaensler, X. H. Sun, E. M. Sadler, A. G. Willis, N. Barry, A. P. Beardsley, M. E. Bell, G. Bernardi, J. D. Bowman, F. Briggs, J. R. Callingham, R. J. Cappallo, P. Carroll, B. E. Corey, A. de Oliveira-Costa, A. A. Deshpande, J. S. Dillon, K. S. Dwarkanath, D. Emrich, A. Ewall-Wice, L. Feng, B. -Q. For, R. Goeke, L. J. Greenhill , et al. (54 additional authors not shown)

Abstract: We present deep polarimetric observations at 154 MHz with the Murchison Widefield Array (MWA), covering 625 deg^2 centered on RA=0 h, Dec=-27 deg. The sensitivity available in our deep observations allows an in-band, frequency-dependent analysis of polarized structure for the first time at long wavelengths. Our analysis suggests that the polarized structures are dominated by intrinsic emission but… ▽ More We present deep polarimetric observations at 154 MHz with the Murchison Widefield Array (MWA), covering 625 deg^2 centered on RA=0 h, Dec=-27 deg. The sensitivity available in our deep observations allows an in-band, frequency-dependent analysis of polarized structure for the first time at long wavelengths. Our analysis suggests that the polarized structures are dominated by intrinsic emission but may also have a foreground Faraday screen component. At these wavelengths, the compactness of the MWA baseline distribution provides excellent snapshot sensitivity to large-scale structure. The observations are sensitive to diffuse polarized emission at ~54' resolution with a sensitivity of 5.9 mJy beam^-1 and compact polarized sources at ~2.4' resolution with a sensitivity of 2.3 mJy beam^-1 for a subset (400 deg^2) of this field. The sensitivity allows the effect of ionospheric Faraday rotation to be spatially and temporally measured directly from the diffuse polarized background. Our observations reveal large-scale structures (~1 deg - 8 deg in extent) in linear polarization clearly detectable in ~2 minute snapshots, which would remain undetectable by interferometers with minimum baseline lengths >110 m at 154 MHz. The brightness temperature of these structures is on average 4 K in polarized intensity, peaking at 11 K. Rotation measure synthesis reveals that the structures have Faraday depths ranging from -2 rad m^-2 to 10 rad m^-2 with a large fraction peaking at ~+1 rad m^-2. We estimate a distance of 51+/-20 pc to the polarized emission based on measurements of the in-field pulsar J2330-2005. We detect four extragalactic linearly polarized point sources within the field in our compact source survey. Based on the known polarized source population at 1.4 GHz and non-detections at 154 MHz, we estimate an upper limit on the depolarization ratio of 0.08 from 1.4 GHz to 154 MHz. △ Less

Submitted 3 August, 2016; v1 submitted 19 July, 2016; originally announced July 2016.

Comments: 32 pages, 20 figures, accepted to ApJ

A High Reliability Survey of Discrete Epoch of Reionization Foreground Sources in the MWA EoR0 Field

Authors: P. A. Carroll, J. Line, M. F. Morales, N. Barry, A. P. Beardsley, B. J. Hazelton, D. C. Jacobs, J. C. Pober, I. S. Sullivan, R. L. Webster, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, B. E. Corey, A. de Oliveira-Costa, J. S. Dillon, D. Emrich, A. Ewall-Wice, L. Feng, B. M. Gaensler, R. Goeke, L. J. Greenhill, J. N. Hewitt, N. Hurley-Walker , et al. (39 additional authors not shown)

Abstract: Detection of the Epoch of Reionization HI signal requires a precise understanding of the intervening galaxies and AGN, both for instrumental calibration and foreground removal. We present a catalogue of 7394 extragalactic sources at 182 MHz detected in the RA=0 field of the Murchison Widefield Array Epoch of Reionization observation programme. Motivated by unprecedented requirements for precision… ▽ More Detection of the Epoch of Reionization HI signal requires a precise understanding of the intervening galaxies and AGN, both for instrumental calibration and foreground removal. We present a catalogue of 7394 extragalactic sources at 182 MHz detected in the RA=0 field of the Murchison Widefield Array Epoch of Reionization observation programme. Motivated by unprecedented requirements for precision and reliability we develop new methods for source finding and selection. We apply machine learning methods to self-consistently classify the relative reliability of 9490 source candidates. A subset of 7466 are selected based on reliability class and signal-to-noise ratio criteria. These are statistically cross-matched to four other radio surveys using both position and flux density information. We find 7369 sources to have confident matches, including 90 partially resolved sources that split into a total of 192 sub-components. An additional 25 unmatched sources are included as new radio detections. The catalogue sources have a median spectral index of -0.85. Spectral flattening is seen toward lower frequencies with a median of -0.71 predicted at 182 MHz. The astrometric error is 7 arcsec. compared to a 2.3 arcmin. beam FWHM. The resulting catalogue covers approximately 1400 sq. deg. and is complete to approximately 80 mJy within half beam power. This provides the most reliable discrete source sky model available to date in the MWA EoR0 field for precision foreground subtraction. △ Less

Submitted 13 July, 2016; originally announced July 2016.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal on June 30, 2016

The radio spectral energy distribution of infrared-faint radio sources

Authors: A. Herzog, R. P. Norris, E. Middelberg, N. Seymour, L. R. Spitler, B. H. C. Emonts, T. M. O. Franzen, R. Hunstead, H. T. Intema, J. Marvil, Q. A. Parker, S. K. Sirothia, N. Hurley-Walker, M. Bell, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, J. R. Callingham, A. A. Deshpande, K. S. Dwarakanath, B. -Q. For, L. J. Greenhill, P. Hancock, B. J. Hazelton , et al. (35 additional authors not shown)

Abstract: Infrared-faint radio sources (IFRS) are a class of radio-loud (RL) active galactic nuclei (AGN) at high redshifts (z > 1.7) that are characterised by their relative infrared faintness, resulting in enormous radio-to-infrared flux density ratios of up to several thousand. We aim to test the hypothesis that IFRS are young AGN, particularly GHz peaked-spectrum (GPS) and compact steep-spectrum (CSS) s… ▽ More Infrared-faint radio sources (IFRS) are a class of radio-loud (RL) active galactic nuclei (AGN) at high redshifts (z > 1.7) that are characterised by their relative infrared faintness, resulting in enormous radio-to-infrared flux density ratios of up to several thousand. We aim to test the hypothesis that IFRS are young AGN, particularly GHz peaked-spectrum (GPS) and compact steep-spectrum (CSS) sources that have a low frequency turnover. We use the rich radio data set available for the Australia Telescope Large Area Survey fields, covering the frequency range between 150 MHz and 34 GHz with up to 19 wavebands from different telescopes, and build radio spectral energy distributions (SEDs) for 34 IFRS. We then study the radio properties of this class of object with respect to turnover, spectral index, and behaviour towards higher frequencies. We also present the highest-frequency radio observations of an IFRS, observed with the Plateau de Bure Interferometer at 105 GHz, and model the multi-wavelength and radio-far-infrared SED of this source. We find IFRS usually follow single power laws down to observed frequencies of around 150 MHz. Mostly, the radio SEDs are steep, but we also find ultra-steep SEDs. In particular, IFRS show statistically significantly steeper radio SEDs than the broader RL AGN population. Our analysis reveals that the fractions of GPS and CSS sources in the population of IFRS are consistent with the fractions in the broader RL AGN population. We find that at least 18% of IFRS contain young AGN, although the fraction might be significantly higher as suggested by the steep SEDs and the compact morphology of IFRS. The detailed multi-wavelength SED modelling of one IFRS shows that it is different from ordinary AGN, although it is consistent with a composite starburst-AGN model with a star formation rate of 170 solar masses per year. △ Less

Submitted 10 July, 2016; originally announced July 2016.

Comments: 25 pages, 7 figures, 5 tables, accepted for publication in A&A

Journal ref: A&A 593, A130 (2016)

Time-domain and spectral properties of pulsars at 154 MHz

Authors: M. E. Bell, Tara Murphy, S. Johnston, D. L. Kaplan, S. Croft, P. Hancock, J. R. Callingham, A. Zic, D. Dobie, J. K. Swiggum, A. Rowlinson, N. Hurley-Walker, A. R. Offringa, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, A. A. Deshpande, B. M. Gaensler, L. J. Greenhill, B. J. Hazelton, M. Johnston-Hollitt, C. J. Lonsdale, S. R. McWhirter, D. A. Mitchell , et al. (13 additional authors not shown)

Abstract: We present 154 MHz Murchison Widefield Array imaging observations and variability information for a sample of pulsars. Over the declination range $-80^{\circ} < δ < 10^{\circ}$ we detect 17 known pulsars with mean flux density greater than 0.3 Jy. We explore the variability properties of this sample on timescales of minutes to years. For three of these pulsars, PSR J0953+0755, PSR J0437-4715 and P… ▽ More We present 154 MHz Murchison Widefield Array imaging observations and variability information for a sample of pulsars. Over the declination range $-80^{\circ} < δ < 10^{\circ}$ we detect 17 known pulsars with mean flux density greater than 0.3 Jy. We explore the variability properties of this sample on timescales of minutes to years. For three of these pulsars, PSR J0953+0755, PSR J0437-4715 and PSR J0630-2834 we observe interstellar scintillation and variability on timescales of greater than 2 minutes. One further pulsar, PSR J0034-0721, showed significant variability, the physical origins of which are difficult to determine. The dynamic spectra for PSR J0953+0755 and PSR J0437-4715 show discrete time and frequency structure consistent with diffractive interstellar scintillation and we present the scintillation bandwidth and timescales from these observations. The remaining pulsars within our sample were statistically non-variable. We also explore the spectral properties of this sample and find spectral curvature in pulsars PSR J0835-4510, PSR J1752-2806 and PSR J0437-4715. △ Less

Submitted 30 May, 2016; originally announced May 2016.

The Murchison Widefield Array 21 cm Power Spectrum Analysis Methodology

Authors: Daniel C. Jacobs, B. J. Hazelton, C. M. Trott, Joshua S. Dillon, B. Pindor, I. S. Sullivan, J. C. Pober, N. Barry, A. P. Beardsley, G. Bernardi, Judd D. Bowman, F. Briggs, R. J. Cappallo, P. Carroll, B. E. Corey, A. de Oliveira-Costa, D. Emrich, A. Ewall-Wice, L. Feng, B. M. Gaensler, R. Goeke, L. J. Greenhill, J. N. Hewitt, N. Hurley-Walker, M. Johnston-Hollitt , et al. (39 additional authors not shown)

Abstract: We present the 21 cm power spectrum analysis approach of the Murchison Widefield Array Epoch of Reionization project. In this paper, we compare the outputs of multiple pipelines for the purpose of validating statistical limits cosmological hydrogen at redshifts between 6 and 12. Multiple, independent, data calibration and reduction pipelines are used to make power spectrum limits on a fiducial nig… ▽ More We present the 21 cm power spectrum analysis approach of the Murchison Widefield Array Epoch of Reionization project. In this paper, we compare the outputs of multiple pipelines for the purpose of validating statistical limits cosmological hydrogen at redshifts between 6 and 12. Multiple, independent, data calibration and reduction pipelines are used to make power spectrum limits on a fiducial night of data. Comparing the outputs of imaging and power spectrum stages highlights differences in calibration, foreground subtraction and power spectrum calculation. The power spectra found using these different methods span a space defined by the various tradeoffs between speed, accuracy, and systematic control. Lessons learned from comparing the pipelines range from the algorithmic to the prosaically mundane; all demonstrate the many pitfalls of neglecting reproducibility. We briefly discuss the way these different methods attempt to handle the question of evaluating a significant detection in the presence of foregrounds. △ Less

Submitted 23 May, 2016; originally announced May 2016.

Comments: accepted to ApJ

A Large Scale, Low Frequency Murchison Widefield Array Survey of Galactic HII regions between 260< l <\340

Authors: L. Hindson, M. Johnston-Hollitt, N. Hurley-Walker, J. R. Callingham, H. Su, J. Morgan, M. Bell, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, A. A. Deshpande, K. S. Dwarakanath, B. -Q For, B. M. Gaensler, L. J. Greenhill, P. Hancock, B. J. Hazelton, A. D. Kapinska, D. L. Kaplan, E. Lenc, C. J. Lonsdale, B. Mckinley, S. R. McWhirter, D. A. Mitchell , et al. (18 additional authors not shown)

Abstract: We have compiled a catalogue of HII regions detected with the Murchison Widefield Array (MWA) between 72 and 231MHz. The multiple frequency bands provided by the MWA allow us identify the characteristic spectrum generated by the thermal Bremsstrahlung process in HII regions. We detect 302 HII regions between 260 < l < 340 and report on the positions, sizes, peak, integrated flux density, and spect… ▽ More We have compiled a catalogue of HII regions detected with the Murchison Widefield Array (MWA) between 72 and 231MHz. The multiple frequency bands provided by the MWA allow us identify the characteristic spectrum generated by the thermal Bremsstrahlung process in HII regions. We detect 302 HII regions between 260 < l < 340 and report on the positions, sizes, peak, integrated flux density, and spectral indices of these HII regions. By identifying the point at which HII regions transition from the optically thin to thick regime we derive the physical properties including the electron density, ionised gas mass and ionising photon flux, towards 61 HII regions. This catalogue of HII regions represents the most extensive and uniform low frequency survey of HII regions in the Galaxy to date. △ Less

Submitted 21 June, 2016; v1 submitted 10 May, 2016; originally announced May 2016.

Comments: 20 pages, 9 figures, 3 tables

First Limits on the 21 cm Power Spectrum during the Epoch of X-ray heating

Authors: A. Ewall-Wice, Joshua S. Dillon, J. N. Hewitt, A. Loeb, A. Mesinger, A. R. Neben, A. R. Offringa, M. Tegmark, N. Barry, A. P. Beardsley, G. Bernardi, Judd D. Bowman, F. Briggs, R. J. Cappallo, P. Carroll, B. E. Corey, A. de Oliveira-Costa, D. Emrich, L. Feng, B. M. Gaensler, R. Goeke, L. J. Greenhill, B. J. Hazelton, N. Hurley-Walker, M. Johnston-Hollitt , et al. (40 additional authors not shown)

Abstract: We present first results from radio observations with the Murchison Widefield Array seeking to constrain the power spectrum of 21 cm brightness temperature fluctuations between the redshifts of 11.6 and 17.9 (113 and 75 MHz). Three hours of observations were conducted over two nights with significantly different levels of ionospheric activity. We use these data to assess the impact of systematic e… ▽ More We present first results from radio observations with the Murchison Widefield Array seeking to constrain the power spectrum of 21 cm brightness temperature fluctuations between the redshifts of 11.6 and 17.9 (113 and 75 MHz). Three hours of observations were conducted over two nights with significantly different levels of ionospheric activity. We use these data to assess the impact of systematic errors at low frequency, including the ionosphere and radio-frequency interference, on a power spectrum measurement. We find that after the 1-3 hours of integration presented here, our measurements at the Murchison Radio Observatory are not limited by RFI, even within the FM band, and that the ionosphere does not appear to affect the level of power in the modes that we expect to be sensitive to cosmology. Power spectrum detections, inconsistent with noise, due to fine spectral structure imprinted on the foregrounds by reflections in the signal-chain, occupy the spatial Fourier modes where we would otherwise be most sensitive to the cosmological signal. We are able to reduce this contamination using calibration solutions derived from autocorrelations so that we achieve an sensitivity of $10^4$ mK on comoving scales $k\lesssim 0.5 h$Mpc$^{-1}$. This represents the first upper limits on the $21$ cm power spectrum fluctuations at redshifts $12\lesssim z \lesssim 18$ but is still limited by calibration systematics. While calibration improvements may allow us to further remove this contamination, our results emphasize that future experiments should consider carefully the existence of and their ability to calibrate out any spectral structure within the EoR window. △ Less

Submitted 16 May, 2016; v1 submitted 29 April, 2016; originally announced May 2016.

Comments: 31 pages, 23 figures, corrected spelling errors in metadata

The 154 MHz radio sky observed by the Murchison Widefield Array: noise, confusion and first source count analyses

Authors: T. M. O. Franzen, C. A. Jackson, A. R. Offringa, R. D. Ekers, R. B. Wayth, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, A. A. Deshpande, B. M. Gaensler, L. J. Greenhill, B. J. Hazelton, M. Johnston-Hollitt, D. L. Kaplan, C. J. Lonsdale, S. R. McWhirter, D. A. Mitchell, M. F. Morales, E. Morgan, J. Morgan, D. Oberoi, S. M. Ord, T. Prabu, N. Seymour , et al. (8 additional authors not shown)

Abstract: We analyse a 154 MHz image made from a 12 h observation with the Murchison Widefield Array (MWA) to determine the noise contribution and behaviour of the source counts down to 30 mJy. The MWA image has a bandwidth of 30.72 MHz, a field-of-view within the half-power contour of the primary beam of 570 deg^2, a resolution of 2.3 arcmin and contains 13,458 sources above 5 sigma. The rms noise in the c… ▽ More We analyse a 154 MHz image made from a 12 h observation with the Murchison Widefield Array (MWA) to determine the noise contribution and behaviour of the source counts down to 30 mJy. The MWA image has a bandwidth of 30.72 MHz, a field-of-view within the half-power contour of the primary beam of 570 deg^2, a resolution of 2.3 arcmin and contains 13,458 sources above 5 sigma. The rms noise in the centre of the image is 4-5 mJy/beam. The MWA counts are in excellent agreement with counts from other instruments and are the most precise ever derived in the flux density range 30-200 mJy due to the sky area covered. Using the deepest available source count data, we find that the MWA image is affected by sidelobe confusion noise at the ~3.5 mJy/beam level, due to incompletely-peeled and out-of-image sources, and classical confusion becomes apparent at ~1.7 mJy/beam. This work highlights that (i) further improvements in ionospheric calibration and deconvolution imaging techniques would be required to probe to the classical confusion limit and (ii) the shape of low-frequency source counts, including any flattening towards lower flux densities, must be determined from deeper ~150 MHz surveys as it cannot be directly inferred from higher frequency data. △ Less

Submitted 13 April, 2016; originally announced April 2016.

Comments: 13 pages, 13 figures, accepted for publication in MNRAS

High-energy sources at low radio frequency: the Murchison Widefield Array view of Fermi blazars

Authors: M. Giroletti, F. Massaro, R. D'Abrusco, R. Lico, D. Burlon, N. Hurley-Walker, M. Johnston-Hollitt, J. Morgan, V. Pavlidou, M. Bell, G. Bernardi, R. Bhat, J. D. Bowman, F. Briggs, R. J. Cappallo, B. E. Corey, A. A. Deshpande, A. Ewall-Rice, D. Emrich, B. M. Gaensler, R. Goeke, L. J. Greenhill, B. J. Hazelton, L. Hindson, D. L. Kaplan , et al. (32 additional authors not shown)

Abstract: Low-frequency radio arrays are opening a new window for the study of the sky, both to study new phenomena and to better characterize known source classes. Being flat-spectrum sources, blazars are so far poorly studied at low radio frequencies. We characterize the spectral properties of the blazar population at low radio frequency compare the radio and high-energy properties of the gamma-ray blaz… ▽ More Low-frequency radio arrays are opening a new window for the study of the sky, both to study new phenomena and to better characterize known source classes. Being flat-spectrum sources, blazars are so far poorly studied at low radio frequencies. We characterize the spectral properties of the blazar population at low radio frequency compare the radio and high-energy properties of the gamma-ray blazar population, and search for radio counterparts of unidentified gamma-ray sources. We cross-correlated the 6,100 deg^2 Murchison Widefield Array Commissioning Survey catalogue with the Roma blazar catalogue, the third catalogue of active galactic nuclei detected by Fermi-LAT, and the unidentified members of the entire third catalogue of gamma-ray sources detected by \fermilat. When available, we also added high-frequency radio data from the Australia Telescope 20 GHz catalogue. We find low-frequency counterparts for 186 out of 517 (36%) blazars, 79 out of 174 (45%) gamma-ray blazars, and 8 out of 73 (11%) gamma-ray blazar candidates. The mean low-frequency (120--180 MHz) blazar spectral index is $\langle α_\mathrm{low} \rangle=0.57\pm0.02$: blazar spectra are flatter than the rest of the population of low-frequency sources, but are steeper than at $\sim$GHz frequencies. Low-frequency radio flux density and gamma-ray energy flux display a mildly significant and broadly scattered correlation. Ten unidentified gamma-ray sources have a (probably fortuitous) positional match with low radio frequency sources. Low-frequency radio astronomy provides important information about sources with a flat radio spectrum and high energy. However, the relatively low sensitivity of the present surveys still misses a significant fraction of these objects. Upcoming deeper surveys, such as the GaLactic and Extragalactic All-Sky MWA (GLEAM) survey, will provide further insight into this population. △ Less

Submitted 29 February, 2016; originally announced February 2016.

Comments: accepted for publication in A&A

Limits on Fast Radio Bursts and other transient sources at 182 MHz using the Murchison Widefield Array

Authors: A. Rowlinson, M. E. Bell, T. Murphy, C. M. Trott, N. Hurley-Walker, S. Johnston, S. J. Tingay, D. L. Kaplan, D. Carbone, P. J. Hancock, L. Feng, A. R. Offringa, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, A. A. Deshpande, B. M. Gaensler, L. J. Greenhill, B. J. Hazelton, M. Johnston-Hollitt, C. J. Lonsdale, S. R. McWhirter, D. A. Mitchell, M. F. Morales , et al. (11 additional authors not shown)

Abstract: We present a survey for transient and variable sources, on timescales from 28 seconds to $\sim$1 year, using the Murchison Widefield Array (MWA) at 182 MHz. Down to a detection threshold of 0.285 Jy, no transient candidates were identified, making this the most constraining low-frequency survey to date and placing a limit on the surface density of transients of $<4.1 \times 10^{-7}$ deg$^{-2}$ for… ▽ More We present a survey for transient and variable sources, on timescales from 28 seconds to $\sim$1 year, using the Murchison Widefield Array (MWA) at 182 MHz. Down to a detection threshold of 0.285 Jy, no transient candidates were identified, making this the most constraining low-frequency survey to date and placing a limit on the surface density of transients of $<4.1 \times 10^{-7}$ deg$^{-2}$ for the shortest timescale considered. At these frequencies, emission from Fast Radio Bursts (FRBs) is expected to be detectable in the shortest timescale images without any corrections for interstellar or intergalactic dispersion. At an FRB limiting flux density of 7980 Jy, we find a rate of $<$82 FRBs per sky per day for dispersion measures $<$700 pc cm$^{-3}$. Assuming a cosmological population of standard candles, our rate limits are consistent with the FRB rates obtained by Thornton et al. (2013) if they have a flat spectral slope. Finally, we conduct an initial variability survey of sources in the field with flux densities $\gtrsim$0.5 Jy and identify no sources with significant variability in their lightcurves. However, we note that substantial further work is required to fully characterise both the short term and low level variability within this field. △ Less

Submitted 24 February, 2016; originally announced February 2016.

Comments: MNRAS Accepted, 17 pages, 10 figures

Beamforming Errors in Murchison Widefield Array Phased Array Antennas and their effects on Epoch of Reionization Science

Authors: A. R. Neben, J. N. Hewitt, R. F. Bradley, J. S. Dillon, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, B. E. Corey, A. A. Deshpande, R. Goeke, L. J. Greenhill, B. J. Hazelton, M. Johnston-Hollitt, D. L. Kaplan, C. J. Lonsdale, S. R. McWhirter, D. A. Mitchell, M. F. Morales, E. Morgan, D. Oberoi, S. M. Ord, T. Prabu, N. Udaya Shankar, K. S. Srivani , et al. (6 additional authors not shown)

Abstract: Accurate antenna beam models are critical for radio observations aiming to isolate the redshifted 21cm spectral line emission from the Dark Ages and the Epoch of Reionization and unlock the scientific potential of 21cm cosmology. Past work has focused on characterizing mean antenna beam models using either satellite signals or astronomical sources as calibrators, but antenna-to-antenna variation d… ▽ More Accurate antenna beam models are critical for radio observations aiming to isolate the redshifted 21cm spectral line emission from the Dark Ages and the Epoch of Reionization and unlock the scientific potential of 21cm cosmology. Past work has focused on characterizing mean antenna beam models using either satellite signals or astronomical sources as calibrators, but antenna-to-antenna variation due to imperfect instrumentation has remained unexplored. We characterize this variation for the Murchison Widefield Array (MWA) through laboratory measurements and simulations, finding typical deviations of order +/- 10-20% near the edges of the main lobe and in the sidelobes. We consider the ramifications of these results for image- and power spectrum-based science. In particular, we simulate visibilities measured by a 100m baseline and find that using an otherwise perfect foreground model, unmodeled beamforming errors severely limit foreground subtraction accuracy within the region of Fourier space contaminated by foreground emission (the "wedge"). This region likely contains much of the cosmological signal, and accessing it will require measurement of per-antenna beam patterns. However, unmodeled beamforming errors do not contaminate the Fourier space region expected to be free of foreground contamination (the "EOR window"), showing that foreground avoidance remains a viable strategy. △ Less

Submitted 16 February, 2016; originally announced February 2016.

Comments: Accepted to ApJ

The Importance of Wide-field Foreground Removal for 21 cm Cosmology: A Demonstration With Early MWA Epoch of Reionization Observations

Authors: J. C. Pober, B. J. Hazelton, A. P. Beardsley, N. A. Barry, Z. E. Martinot, I. S. Sullivan, M. F. Morales, M. E. Bell, G. Bernardi, N. D. R. Bhat, J. D. Bowman, F. Briggs, R. J. Cappallo, P. Carroll, B. E. Corey, A. de Oliveira-Costa, A. A. Deshpande, Joshua. S. Dillon, D. Emrich, A. M. Ewall-Wice, L. Feng, R. Goeke, L. J. Greenhill, J. N. Hewitt, L. Hindson , et al. (44 additional authors not shown)

Abstract: In this paper we present observations, simulations, and analysis demonstrating the direct connection between the location of foreground emission on the sky and its location in cosmological power spectra from interferometric redshifted 21 cm experiments. We begin with a heuristic formalism for understanding the mapping of sky coordinates into the cylindrically averaged power spectra measurements us… ▽ More In this paper we present observations, simulations, and analysis demonstrating the direct connection between the location of foreground emission on the sky and its location in cosmological power spectra from interferometric redshifted 21 cm experiments. We begin with a heuristic formalism for understanding the mapping of sky coordinates into the cylindrically averaged power spectra measurements used by 21 cm experiments, with a focus on the effects of the instrument beam response and the associated sidelobes. We then demonstrate this mapping by analyzing power spectra with both simulated and observed data from the Murchison Widefield Array. We find that removing a foreground model which includes sources in both the main field-of-view and the first sidelobes reduces the contamination in high k_parallel modes by several percent relative to a model which only includes sources in the main field-of-view, with the completeness of the foreground model setting the principal limitation on the amount of power removed. While small, a percent-level amount of foreground power is in itself more than enough to prevent recovery of any EoR signal from these modes. This result demonstrates that foreground subtraction for redshifted 21 cm experiments is truly a wide-field problem, and algorithms and simulations must extend beyond the main instrument field-of-view to potentially recover the full 21 cm power spectrum. △ Less

Submitted 22 January, 2016; originally announced January 2016.

Comments: 13 pages, 9 figures, matches version accepted to ApJ

CHIPS: The Cosmological HI Power Spectrum Estimator

Authors: Cathryn M. Trott, Bart Pindor, Pietro Procopio, Randall B. Wayth, Daniel A. Mitchell, Benjamin McKinley, Steven J. Tingay, N. Barry, A. P. Beardsley, G. Bernardi, Judd D. Bowman, F. Briggs, R. J. Cappallo, P. Carroll, A. de Oliveira-Costa, Joshua S. Dillon, A. Ewall-Wice, L. Feng, L. J. Greenhill, B. J. Hazelton, J. N. Hewitt, N. Hurley-Walker, M. Johnston-Hollitt, Daniel C. Jacobs, D. L. Kaplan , et al. (27 additional authors not shown)

Abstract: Detection of the cosmological neutral hydrogen signal from the Epoch of Reionization, and estimation of its basic physical parameters, is the principal scientific aim of many current low-frequency radio telescopes. Here we describe the Cosmological HI Power Spectrum Estimator (CHIPS), an algorithm developed and implemented with data from the Murchison Widefield Array (MWA), to compute the two-dime… ▽ More Detection of the cosmological neutral hydrogen signal from the Epoch of Reionization, and estimation of its basic physical parameters, is the principal scientific aim of many current low-frequency radio telescopes. Here we describe the Cosmological HI Power Spectrum Estimator (CHIPS), an algorithm developed and implemented with data from the Murchison Widefield Array (MWA), to compute the two-dimensional and spherically-averaged power spectrum of brightness temperature fluctuations. The principal motivations for CHIPS are the application of realistic instrumental and foreground models to form the optimal estimator, thereby maximising the likelihood of unbiased signal estimation, and allowing a full covariant understanding of the outputs. CHIPS employs an inverse-covariance weighting of the data through the maximum likelihood estimator, thereby allowing use of the full parameter space for signal estimation ("foreground suppression"). We describe the motivation for the algorithm, implementation, application to real and simulated data, and early outputs. Upon application to a set of 3 hours of data, we set a 2$σ$ upper limit on the EoR dimensionless power at $k=0.05$~h.Mpc$^{-1}$ of $Δ_k^2<7.6\times{10^4}$~mK$^2$ in the redshift range $z=[6.2-6.6]$, consistent with previous estimates. △ Less

Submitted 8 January, 2016; originally announced January 2016.

Comments: 34 pages, accepted for publication in ApJ

A search for Fast Radio Bursts at low frequencies with Murchison Widefield Array high time resolution imaging

Authors: S. J. Tingay, C. M. Trott, R. B. Wayth, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, A. A. Deshpande, L. Feng, B. M. Gaensler, L. J. Greenhill, P. J. Hancock, B. J. Hazelton, M. Johnston-Hollitt, D. L. Kaplan, C. J. Lonsdale, S. R. McWhirter, D. A. Mitchell, M. F. Morales, E. Morgan, T. Murphy, D. Oberoi, T. Prabu, N. Udaya Shankar, K. S. Srivani , et al. (4 additional authors not shown)

Abstract: We present the results of a pilot study search for Fast Radio Bursts (FRBs) using the Murchison Widefield Array (MWA) at low frequencies (139 - 170 MHz). We utilised MWA data obtained in a routine imaging mode from observations where the primary target was a field being studied for Epoch of Reionisation detection. We formed images with 2 second time resolution and 1.28~MHz frequency resolution for… ▽ More We present the results of a pilot study search for Fast Radio Bursts (FRBs) using the Murchison Widefield Array (MWA) at low frequencies (139 - 170 MHz). We utilised MWA data obtained in a routine imaging mode from observations where the primary target was a field being studied for Epoch of Reionisation detection. We formed images with 2 second time resolution and 1.28~MHz frequency resolution for 10.5 hours of observations, over 400 square degrees of the sky. We de-dispersed the dynamic spectrum in each of 372,100 resolution elements of 2$\times$2 arcmin$^{2}$, between dispersion measures of 170 and 675~pc~cm$^{-3}$. Based on the event rate calculations in Trott, Tingay & Wayth (2013), which assumes a standard candle luminosity of $8\times10^{37}$ Js$^{-1}$, we predict that with this choice of observational parameters, the MWA should detect ($\sim10$,$\sim2$,$\sim0$) FRBs with spectral indices corresponding to ($-$2, $-$1, 0), based on a 7$σ$ detection threshold. We find no FRB candidates above this threshold from our search, placing an event rate limit of $<700$ above 700 Jy.ms per day per sky and providing evidence against spectral indices $α<-1.2$ ($S\proptoν^α$). We compare our event rate and spectral index limits with others from the literature. We briefly discuss these limits in light of recent suggestions that supergiant pulses from young neutron stars could explain FRBs. We find that such supergiant pulses would have to have much flatter spectra between 150 and 1400 MHz than have been observed from Crab giant pulses to be consistent with the FRB spectral index limit we derive. △ Less

Submitted 9 November, 2015; originally announced November 2015.

Comments: 14 pages, 5 figures

Quantifying ionospheric effects on time-domain astrophysics with the Murchison Widefield Array

Authors: Shyeh Tjing Loi, Tara Murphy, Martin E. Bell, David L. Kaplan, Emil Lenc, Andre R. Offringa, Natasha Hurley-Walker, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, B. E. Corey, A. A. Deshpande, D. Emrich, B. M. Gaensler, R. Goeke, L. J. Greenhill, B. J. Hazelton, M. Johnston-Hollitt, J. C. Kasper, E. Kratzenberg, C. J. Lonsdale, M. J. Lynch, S. R. McWhirter, D. A. Mitchell , et al. (17 additional authors not shown)

Abstract: Refraction and diffraction of incoming radio waves by the ionosphere induce time variability in the angular positions, peak amplitudes and shapes of radio sources, potentially complicating the automated cross-matching and identification of transient and variable radio sources. In this work, we empirically assess the effects of the ionosphere on data taken by the Murchison Widefield Array (MWA) rad… ▽ More Refraction and diffraction of incoming radio waves by the ionosphere induce time variability in the angular positions, peak amplitudes and shapes of radio sources, potentially complicating the automated cross-matching and identification of transient and variable radio sources. In this work, we empirically assess the effects of the ionosphere on data taken by the Murchison Widefield Array (MWA) radio telescope. We directly examine 51 hours of data observed over 10 nights under quiet geomagnetic conditions (global storm index Kp < 2), analysing the behaviour of short-timescale angular position and peak flux density variations of around ten thousand unresolved sources. We find that while much of the variation in angular position can be attributed to ionospheric refraction, the characteristic displacements (10-20 arcsec) at 154 MHz are small enough that search radii of 1-2 arcmin should be sufficient for cross-matching under typical conditions. By examining bulk trends in amplitude variability, we place upper limits on the modulation index associated with ionospheric scintillation of 1-3% for the various nights. For sources fainter than ~1 Jy, this variation is below the image noise at typical MWA sensitivities. Our results demonstrate that the ionosphere is not a significant impediment to the goals of time-domain science with the MWA at 154 MHz. △ Less

Submitted 4 August, 2015; originally announced August 2015.

Comments: Accepted for publication in MNRAS

Murchison Widefield Array Observations of Anomalous Variability: A Serendipitous Night-time Detection of Interplanetary Scintillation

Authors: D. L. Kaplan, S. J. Tingay, P. K. Manoharan, J. -P. Macquart, P. Hancock, J. Morgan, D. A. Mitchell, R. D. Ekers, R. B. Wayth, C. Trott, T. Murphy, D. Oberoi, I. H. Cairns, L. Feng, N. Kudryavtseva, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, A. A. Deshpande, B. M. Gaensler, L. J. Greenhill, N. Hurley-Walker, B. J. Hazelton, M. Johnston-Hollitt , et al. (11 additional authors not shown)

Abstract: We present observations of high-amplitude rapid (2 s) variability toward two bright, compact extragalactic radio sources out of several hundred of the brightest radio sources in one of the 30x30 deg MWA Epoch of Reionization fields using the Murchison Widefield Array (MWA) at 155 MHz. After rejecting intrinsic, instrumental, and ionospheric origins we consider the most likely explanation for this… ▽ More We present observations of high-amplitude rapid (2 s) variability toward two bright, compact extragalactic radio sources out of several hundred of the brightest radio sources in one of the 30x30 deg MWA Epoch of Reionization fields using the Murchison Widefield Array (MWA) at 155 MHz. After rejecting intrinsic, instrumental, and ionospheric origins we consider the most likely explanation for this variability to be interplanetary scintillation (IPS), likely the result of a large coronal mass ejection propagating from the Sun. This is confirmed by roughly contemporaneous observations with the Ooty Radio Telescope. We see evidence for structure on spatial scales ranging from <1000 km to >1e6 km. The serendipitous night-time nature of these detections illustrates the new regime that the MWA has opened for IPS studies with sensitive night-time, wide-field, low-frequency observations. This regime complements traditional dedicated strategies for observing IPS and can be utilized in real-time to facilitate dedicated follow-up observations. At the same time, it allows large-scale surveys for compact (arcsec) structures in low-frequency radio sources despite the 2 arcmin resolution of the array. △ Less

Submitted 27 July, 2015; originally announced July 2015.

Comments: 7 pages, 4 figures. ApJ Letters, in press

Simultaneous Observations of Giant Pulses from the Crab Pulsar, with the Murchison Widefield Array and Parkes Radio Telescope: Implications for the Giant Pulse Emission Mechanism

Authors: S. I. Oronsaye, S. M. Ord, N. D. R. Bhat, S. E. Tremblay, S. J. McSweeney, S. J. Tingay, W. van Straten, A. Jameson, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, A. A. Deshpande, L. J. Greenhill, B. J. Hazelton, M. Johnston-Hollitt, D. L. Kaplan, C. J. Lonsdale, S. R. McWhirter, D. A. Mitchell, M. F. Morales, E. Morgan, D. Oberoi, T. Prabu, N. Udaya Shankar , et al. (6 additional authors not shown)

Abstract: We report on observations of giant pulses from the Crab pulsar performed simultaneously with the Parkes radio telescope and the incoherent combination of the Murchison Widefield Array (MWA) antenna tiles. The observations were performed over a duration of approximately one hour at a center frequency of 1382 MHz with 340 MHz bandwidth at Parkes, and at a center frequency of 193 MHz with 15 MHz band… ▽ More We report on observations of giant pulses from the Crab pulsar performed simultaneously with the Parkes radio telescope and the incoherent combination of the Murchison Widefield Array (MWA) antenna tiles. The observations were performed over a duration of approximately one hour at a center frequency of 1382 MHz with 340 MHz bandwidth at Parkes, and at a center frequency of 193 MHz with 15 MHz bandwidth at the MWA. Our analysis has led to the detection of 55 giant pulses at the MWA and 2075 at Parkes above a threshold of 3.5$σ$ and 6.5$σ$ respectively. We detected 51$\%$ of the MWA giant pulses at the Parkes radio telescope, with spectral indices in the range of $-3.6>α> -4.9$ ($S_{\rm ν} \propto ν^α$). We present a Monte Carlo analysis supporting the conjecture that the giant pulse emission in the Crab is intrinsically broadband, the less than $100\%$ correlation being due to the relative sensitivities of the two instruments and the width of the spectral index distribution. Our observations are consistent with the hypothesis that the spectral index of giant pulses is drawn from normal distribution of standard deviation 0.6, but with a mean that displays an evolution with frequency from -3.00 at 1382 MHz, to -2.85 at 192 MHz. △ Less

Submitted 17 July, 2015; originally announced July 2015.

Broadband Spectral Modeling of the Extreme Gigahertz-Peaked Spectrum Radio Source PKS B0008-421

Authors: J. R. Callingham, B. M. Gaensler, R. D. Ekers, S. J. Tingay, R. B. Wayth, J. Morgan, G. Bernardi, M. E. Bell, R. Bhat, J. D. Bowman, F. Briggs, R. J. Cappallo, A. A. Deshpande, A. Ewall-Wice, L. Feng, L. J. Greenhill, B. J. Hazelton, L. Hindson, N. Hurley-Walker, D. C. Jacobs, M. Johnston-Hollitt, D. L. Kaplan, N. Kudrayavtseva, E. Lenc, C. J. Lonsdale , et al. (18 additional authors not shown)

Abstract: We present broadband observations and spectral modeling of PKS B0008-421, and identify it as an extreme gigahertz-peaked spectrum (GPS) source. PKS B0008-421 is characterized by the steepest known spectral slope below the turnover, close to the theoretical limit of synchrotron self-absorption, and the smallest known spectral width of any GPS source. Spectral coverage of the source spans from 0.118… ▽ More We present broadband observations and spectral modeling of PKS B0008-421, and identify it as an extreme gigahertz-peaked spectrum (GPS) source. PKS B0008-421 is characterized by the steepest known spectral slope below the turnover, close to the theoretical limit of synchrotron self-absorption, and the smallest known spectral width of any GPS source. Spectral coverage of the source spans from 0.118 to 22 GHz, which includes data from the Murchison Widefield Array and the wide bandpass receivers on the Australia Telescope Compact Array. We have implemented a Bayesian inference model fitting routine to fit the data with various absorption models. We find that without the inclusion of a high-frequency exponential break the absorption models can not accurately fit the data, with significant deviations above and below the peak in the radio spectrum. The addition of a high-frequency break provides acceptable spectral fits for the inhomogeneous free-free absorption and double-component synchrotron self-absorption models, with the inhomogeneous free-free absorption model statistically favored. The requirement of a high-frequency spectral break implies that the source has ceased injecting fresh particles. Additional support for the inhomogeneous free-free absorption model as being responsible for the turnover in the spectrum is given by the consistency between the physical parameters derived from the model fit and the implications of the exponential spectral break, such as the necessity of the source being surrounded by a dense ambient medium to maintain the peak frequency near the gigahertz region. The discovery of PKS B0008-421 suggests that the next generation of low radio frequency surveys could reveal a large population of GPS sources that have ceased activity, and that a portion of the ultra-steep spectrum source population could be composed of these GPS sources in a relic phase. △ Less

Submitted 16 July, 2015; originally announced July 2015.

Comments: Accepted for publication in ApJ

Ionospheric modelling using GPS to calibrate the MWA. 1: Comparison of first order ionospheric effects between GPS models and MWA observations

Authors: B. S. Arora, J. Morgan, S. M. Ord, S. J. Tingay, N. Hurley-Walker, M. Bell, G. Bernardi, R. Bhat, F. Briggs, J. R. Callingham, A. A. Deshpande, K. S. Dwarakanath, A. Ewall-Wice, L. Feng, B. -Q. For, P. Hancock, B. J. Hazelton, L. Hindson, D. Jacobs, M. Johnston-Hollitt, A. D. Kapińska, N. Kudryavtseva, E. Lenc, B. McKinley, D. Mitchell , et al. (34 additional authors not shown)

Abstract: We compare first order (refractive) ionospheric effects seen by the Murchison Widefield Array (MWA) with the ionosphere as inferred from Global Positioning System (GPS) data. The first order ionosphere manifests itself as a bulk position shift of the observed sources across an MWA field of view. These effects can be computed from global ionosphere maps provided by GPS analysis centres, namely the… ▽ More We compare first order (refractive) ionospheric effects seen by the Murchison Widefield Array (MWA) with the ionosphere as inferred from Global Positioning System (GPS) data. The first order ionosphere manifests itself as a bulk position shift of the observed sources across an MWA field of view. These effects can be computed from global ionosphere maps provided by GPS analysis centres, namely the Center for Orbit Determination in Europe (CODE), using data from globally distributed GPS receivers. However, for the more accurate local ionosphere estimates required for precision radio astronomy applications, data from local GPS networks needs to be incorporated into ionospheric modelling. For GPS observations, the ionospheric parameters are biased by GPS receiver instrument delays, among other effects, also known as receiver Differential Code Biases (DCBs). The receiver DCBs need to be estimated for any non-CODE GPS station used for ionosphere modelling, a requirement for establishing dense GPS networks in arbitrary locations in the vicinity of the MWA. In this work, single GPS station-based ionospheric modelling is performed at a time resolution of 10 minutes. Also the receiver DCBs are estimated for selected Geoscience Australia (GA) GPS receivers, located at Murchison Radio Observatory (MRO1), Yarragadee (YAR3), Mount Magnet (MTMA) and Wiluna (WILU). The ionospheric gradients estimated from GPS are compared with the ionospheric gradients inferred from radio source position shifts observed with the MWA. The ionospheric gradients at all the GPS stations show a correlation with the gradients observed with the MWA. The ionosphere estimates obtained using GPS measurements show promise in terms of providing calibration information for the MWA. △ Less

Submitted 5 July, 2015; originally announced July 2015.

Comments: 24 pages

Confirmation of Wide-Field Signatures in Redshifted 21 cm Power Spectra

Authors: Nithyanandan Thyagarajan, Daniel C. Jacobs, Judd D. Bowman, N. Barry, A. P. Beardsley, G. Bernardi, F. Briggs, R. J. Cappallo, P. Carroll, A. A. Deshpande, A. de Oliveira-Costa, Joshua S. Dillon, A. Ewall-Wice, L. Feng, L. J. Greenhill, B. J. Hazelton, L. Hernquist, J. N. Hewitt, N. Hurley-Walker, M. Johnston-Hollitt, D. L. Kaplan, Han-Seek Kim, P. Kittiwisit, E. Lenc, J. Line , et al. (30 additional authors not shown)

Abstract: We confirm our recent prediction of the "pitchfork" foreground signature in power spectra of high-redshift 21 cm measurements where the interferometer is sensitive to large-scale structure on all baselines. This is due to the inherent response of a wide-field instrument and is characterized by enhanced power from foreground emission in Fourier modes adjacent to those considered to be the most sens… ▽ More We confirm our recent prediction of the "pitchfork" foreground signature in power spectra of high-redshift 21 cm measurements where the interferometer is sensitive to large-scale structure on all baselines. This is due to the inherent response of a wide-field instrument and is characterized by enhanced power from foreground emission in Fourier modes adjacent to those considered to be the most sensitive to the cosmological H I signal. In our recent paper, many signatures from the simulation that predicted this feature were validated against Murchison Widefield Array (MWA) data, but this key pitchfork signature was close to the noise level. In this paper, we improve the data sensitivity through the coherent averaging of 12 independent snapshots with identical instrument settings and provide the first confirmation of the prediction with a signal-to-noise ratio > 10. This wide-field effect can be mitigated by careful antenna designs that suppress sensitivity near the horizon. Simple models for antenna apertures that have been proposed for future instruments such as the Hydrogen Epoch of Reionization Array and the Square Kilometre Array indicate they should suppress foreground leakage from the pitchfork by ~40 dB relative to the MWA and significantly increase the likelihood of cosmological signal detection in these critical Fourier modes in the three-dimensional power spectrum. △ Less

Submitted 3 July, 2015; v1 submitted 19 June, 2015; originally announced June 2015.

Comments: Accepted for publication in ApJ Letters, 6 pages, 2 figures

Power spectrum analysis of ionospheric fluctuations with the Murchison Widefield Array

Authors: Shyeh Tjing Loi, Cathryn M. Trott, Tara Murphy, Iver H. Cairns, Martin Bell, Natasha Hurley-Walker, John Morgan, Emil Lenc, A. R. Offringa, L. Feng, P. J. Hancock, D. L. Kaplan, N. Kudryavtseva, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, B. E. Corey, A. A. Deshpande, D. Emrich, B. M. Gaensler, R. Goeke, L. J. Greenhill, B. J. Hazelton, M. Johnston-Hollitt , et al. (23 additional authors not shown)

Abstract: Low-frequency, wide field-of-view (FoV) radio telescopes such as the Murchison Widefield Array (MWA) enable the ionosphere to be sampled at high spatial completeness. We present the results of the first power spectrum analysis of ionospheric fluctuations in MWA data, where we examined the position offsets of radio sources appearing in two datasets. The refractive shifts in the positions of celesti… ▽ More Low-frequency, wide field-of-view (FoV) radio telescopes such as the Murchison Widefield Array (MWA) enable the ionosphere to be sampled at high spatial completeness. We present the results of the first power spectrum analysis of ionospheric fluctuations in MWA data, where we examined the position offsets of radio sources appearing in two datasets. The refractive shifts in the positions of celestial sources are proportional to spatial gradients in the electron column density transverse to the line of sight. These can be used to probe plasma structures and waves in the ionosphere. The regional (10-100 km) scales probed by the MWA, determined by the size of its FoV and the spatial density of radio sources (typically thousands in a single FoV), complement the global (100-1000 km) scales of GPS studies and local (0.01-1 km) scales of radar scattering measurements. Our data exhibit a range of complex structures and waves. Some fluctuations have the characteristics of travelling ionospheric disturbances (TIDs), while others take the form of narrow, slowly-drifting bands aligned along the Earth's magnetic field. △ Less

Submitted 5 June, 2015; originally announced June 2015.

Comments: Accepted for publication in Radio Science

Empirical Covariance Modeling for 21 cm Power Spectrum Estimation: A Method Demonstration and New Limits from Early Murchison Widefield Array 128-Tile Data

Authors: Joshua S. Dillon, Abraham R. Neben, Jacqueline N. Hewitt, Max Tegmark, N. Barry, A. P. Beardsley, J. D. Bowman, F. Briggs, P. Carroll, A. de Oliveira-Costa, A. Ewall-Wice, L. Feng, L. J. Greenhill, B. J. Hazelton, L. Hernquist, N. Hurley-Walker, D. C. Jacobs, H. S. Kim, P. Kittiwisit, E. Lenc, J. Line, A. Loeb, B. McKinley, D. A. Mitchell, M. F. Morales , et al. (30 additional authors not shown)

Abstract: The separation of the faint cosmological background signal from bright astrophysical foregrounds remains one of the most daunting challenges of mapping the high-redshift intergalactic medium with the redshifted 21 cm line of neutral hydrogen. Advances in mapping and modeling of diffuse and point source foregrounds have improved subtraction accuracy, but no subtraction scheme is perfect. Precisely… ▽ More The separation of the faint cosmological background signal from bright astrophysical foregrounds remains one of the most daunting challenges of mapping the high-redshift intergalactic medium with the redshifted 21 cm line of neutral hydrogen. Advances in mapping and modeling of diffuse and point source foregrounds have improved subtraction accuracy, but no subtraction scheme is perfect. Precisely quantifying the errors and error correlations due to missubtracted foregrounds allows for both the rigorous analysis of the 21 cm power spectrum and for the maximal isolation of the "EoR window" from foreground contamination. We present a method to infer the covariance of foreground residuals from the data itself in contrast to previous attempts at a priori modeling. We demonstrate our method by setting limits on the power spectrum using a 3 h integration from the 128-tile Murchison Widefield Array. Observing between 167 and 198 MHz, we find at 95% confidence a best limit of Delta^2(k) < 3.7 x 10^4 mK^2 at comoving scale k = 0.18 hMpc^-1 and at z = 6.8, consistent with existing limits. △ Less

Submitted 2 June, 2015; originally announced June 2015.

Comments: 18 pages, 9 figures. Accepted to Phys. Rev. D

Journal ref: Phys. Rev. D 91, 123011 (19 June 2015)

An analysis of the halo and relic radio emission from Abell 3376 from Murchison Widefield Array observations

Authors: Lijo T. George, K. S. Dwarakanath, M. Johnston-Hollitt, N. Hurley-Walker, L. Hindson, A. D. Kapińska, S. J. Tingay, M. Bell, J. R. Callingham, Bi-Qing For, P. J. Hancock, E. Lenc, B. McKinley, J. Morgan, A. Offringa, P. Procopio, L. Staveley-Smith, R. B. Wayth, Chen Wu, Q. Zheng, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, B. E. Corey , et al. (26 additional authors not shown)

Abstract: We have carried out multiwavelength observations of the near-by ($z=0.046$) rich, merging galaxy cluster Abell 3376 with the Murchison Widefield Array (MWA). As a part of the GaLactic and Extragalactic All-sky MWA survey (GLEAM), this cluster was observed at 88, 118, 154, 188 and 215 MHz. The known radio relics, towards the eastern and western peripheries of the cluster, were detected at all the f… ▽ More We have carried out multiwavelength observations of the near-by ($z=0.046$) rich, merging galaxy cluster Abell 3376 with the Murchison Widefield Array (MWA). As a part of the GaLactic and Extragalactic All-sky MWA survey (GLEAM), this cluster was observed at 88, 118, 154, 188 and 215 MHz. The known radio relics, towards the eastern and western peripheries of the cluster, were detected at all the frequencies. The relics, with a linear extent of $\sim$ 1 Mpc each, are separated by $\sim$ 2 Mpc. Combining the current observations with those in the literature, we have obtained the spectra of these relics over the frequency range 80 -- 1400 MHz. The spectra follow power laws, with $α$ = $-1.17\pm0.06$ and $-1.37\pm0.08$ for the west and east relics, respectively ($S \propto ν^α$). Assuming the break frequency to be near the lower end of the spectrum we estimate the age of the relics to be $\sim$ 0.4 Gyr. No diffuse radio emission from the central regions of the cluster (halo) was detected. The upper limit on the radio power of any possible halo that might be present in the cluster is a factor of 35 lower than that expected from the radio power and X-ray luminosity correlation for cluster halos. From this we conclude that the cluster halo is very extended ($>$ 500 kpc) and/or most of the radio emission from the halo has decayed. The current limit on the halo radio power is a factor of ten lower than the existing upper limits with possible implications for models of halo formation. △ Less

Submitted 1 June, 2015; originally announced June 2015.

Comments: 10 pages, 6 figures

Measuring Phased-Array Antenna Beampatterns with High Dynamic Range for the Murchison Widefield Array using 137 MHz ORBCOMM Satellites

Authors: A. R. Neben, R. F. Bradley, J. N. Hewitt, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, A. A. Deshpande, R. Goeke, L. J. Greenhill, B. J. Hazelton, M. Johnston-Hollitt, D. L. Kaplan, C. J. Lonsdale, S. R. McWhirter, D. A. Mitchell, M. F. Morales, E. Morgan, D. Oberoi, S. M. Ord, T. Prabu, N. Udaya Shankar, K. S. Srivani, R. Subrahmanyan, S. J. Tingay , et al. (4 additional authors not shown)

Abstract: Detection of the fluctuations in 21 cm line emission from neutral hydrogen during the Epoch of Reionization in thousand hour integrations poses stringent requirements on calibration and image quality, both of which necessitate accurate primary beam models. The Murchison Widefield Array (MWA) uses phased array antenna elements which maximize collecting area at the cost of complexity. To quantify th… ▽ More Detection of the fluctuations in 21 cm line emission from neutral hydrogen during the Epoch of Reionization in thousand hour integrations poses stringent requirements on calibration and image quality, both of which necessitate accurate primary beam models. The Murchison Widefield Array (MWA) uses phased array antenna elements which maximize collecting area at the cost of complexity. To quantify their performance, we have developed a novel beam measurement system using the 137 MHz ORBCOMM satellite constellation and a reference dipole antenna. Using power ratio measurements, we measure the {\it in situ} beampattern of the MWA antenna tile relative to that of the reference antenna, canceling the variation of satellite flux or polarization with time. We employ angular averaging to mitigate multipath effects (ground scattering), and assess environmental systematics with a null experiment in which the MWA tile is replaced with a second reference dipole. We achieve beam measurements over 30 dB dynamic range in beam sensitivity over a large field of view (65\% of the visible sky), far wider and deeper than drift scans through astronomical sources allow. We verify an analytic model of the MWA tile at this frequency within a few percent statistical scatter within the full width at half maximum. Towards the edges of the main lobe and in the sidelobes, we measure tens of percent systematic deviations. We compare these errors with those expected from known beamforming errors. △ Less

Submitted 26 May, 2015; originally announced May 2015.

Comments: 22 pages, 8 figures. Accepted for publication in Radio Science

GLEAM: The GaLactic and Extragalactic All-sky MWA survey

Authors: R. B. Wayth, E. Lenc, M. E. Bell, J. R. Callingham, K. S. Dwarakanath, T. M. O. Franzen, B. -Q. For, B. Gaensler, P. Hancock, L. Hindson, N. Hurley-Walker, C. A. Jackson, M. Johnston-Hollitt, A. D. Kapinska, B. McKinley, J. Morgan, A. R. Offringa, P. Procopio, L. Staveley-Smith, C. Wu, Q. Zheng, C. M. Trott, G. Bernardi, J. D. Bowman, F. Briggs , et al. (30 additional authors not shown)

Abstract: GLEAM, the GaLactic and Extragalactic All-sky MWA survey, is a survey of the entire radio sky south of declination +25 deg at frequencies between 72 and 231 MHz, made with the Murchison Widefield Array (MWA) using a drift scan method that makes efficient use of the MWA's very large field-of-view. We present the observation details, imaging strategies and theoretical sensitivity for GLEAM. The surv… ▽ More GLEAM, the GaLactic and Extragalactic All-sky MWA survey, is a survey of the entire radio sky south of declination +25 deg at frequencies between 72 and 231 MHz, made with the Murchison Widefield Array (MWA) using a drift scan method that makes efficient use of the MWA's very large field-of-view. We present the observation details, imaging strategies and theoretical sensitivity for GLEAM. The survey ran for two years, the first year using 40 kHz frequency resolution and 0.5 s time resolution; the second year using 10 kHz frequency resolution and 2 s time resolution. The resulting image resolution and sensitivity depends on observing frequency, sky pointing and image weighting scheme. At 154 MHz the image resolution is approximately 2.5 x 2.2/cos(DEC+26.7) arcmin with sensitivity to structures up to ~10 deg in angular size. We provide tables to calculate the expected thermal noise for GLEAM mosaics depending on pointing and frequency and discuss limitations to achieving theoretical noise in Stokes I images. We discuss challenges, and their solutions, that arise for GLEAM including ionospheric effects on source positions and linearly polarised emission, and the instrumental polarisation effects inherent to the MWA's primary beam. △ Less

Submitted 22 May, 2015; originally announced May 2015.

Comments: 17 pages, 7 Figures, 4 tables. Accepted for publication in the Publications of the Astronomical Society of Australia (PASA)

Real-time imaging of density ducts between the plasmasphere and ionosphere

Authors: Shyeh Tjing Loi, Tara Murphy, Iver H. Cairns, Frederick W. Menk, Colin L. Waters, Philip J. Erickson, Cathryn M. Trott, Natasha Hurley-Walker, John Morgan, Emil Lenc, Andre R. Offringa, Martin E. Bell, Ronald D. Ekers, B. M. Gaensler, Colin J. Lonsdale, Lu Feng, Paul J. Hancock, David L. Kaplan, G. Bernardi, J. D. Bowman, F. Briggs, R. J. Cappallo, A. A. Deshpande, L. J. Greenhill, B. J. Hazelton , et al. (16 additional authors not shown)

Abstract: Ionization of the Earth's atmosphere by sunlight forms a complex, multi-layered plasma environment within the Earth's magnetosphere, the innermost layers being the ionosphere and plasmasphere. The plasmasphere is believed to be embedded with cylindrical density structures (ducts) aligned along the Earth's magnetic field, but direct evidence for these remains scarce. Here we report the first direct… ▽ More Ionization of the Earth's atmosphere by sunlight forms a complex, multi-layered plasma environment within the Earth's magnetosphere, the innermost layers being the ionosphere and plasmasphere. The plasmasphere is believed to be embedded with cylindrical density structures (ducts) aligned along the Earth's magnetic field, but direct evidence for these remains scarce. Here we report the first direct wide-angle observation of an extensive array of field-aligned ducts bridging the upper ionosphere and inner plasmasphere, using a novel ground-based imaging technique. We establish their heights and motions by feature-tracking and parallax analysis. The structures are strikingly organized, appearing as regularly-spaced, alternating tubes of overdensities and underdensities strongly aligned with the Earth's magnetic field. These findings represent the first direct visual evidence for the existence of such structures. △ Less

Submitted 24 April, 2015; originally announced April 2015.

Comments: Accepted for publication in Geophysical Research Letters