Showing 1–50 of 280 results for author: Bhat, N

The neutron star mass, distance, and inclination from precision timing of the brilliant millisecond pulsar J0437$-$4715

Authors: Daniel J. Reardon, Matthew Bailes, Ryan M. Shannon, Chris Flynn, Jacob Askew, N. D. Ramesh Bhat, Zu-Cheng Chen, Małgorzata Curyło, Yi Feng, George B. Hobbs, Agastya Kapur, Matthew Kerr, Xiaojin Liu, Richard N. Manchester, Rami Mandow, Saurav Mishra, Christopher J. Russell, Mohsen Shamohammadi, Lei Zhang, Andrew Zic

Abstract: The observation of neutron stars enables the otherwise impossible study of fundamental physical processes. Timing of binary radio pulsars is particularly powerful, as it enables precise characterization of their (three-dimensional) positions and orbits. PSR J0437$-$4715 is an important millisecond pulsar for timing array experiments and is also a primary target for the Neutron Star Interior Compos… ▽ More The observation of neutron stars enables the otherwise impossible study of fundamental physical processes. Timing of binary radio pulsars is particularly powerful, as it enables precise characterization of their (three-dimensional) positions and orbits. PSR J0437$-$4715 is an important millisecond pulsar for timing array experiments and is also a primary target for the Neutron Star Interior Composition ExploreR (NICER). The main aim of the NICER mission is to constrain the neutron star equation of state by inferring the compactness ($M_p/R$) of the star. Direct measurements of the mass $M_p$ from pulsar timing therefore substantially improve constraints on the radius $R$, and the equation of state. Here we use observations spanning 26 years from Murriyang, the 64-m Parkes radio telescope, to improve the timing model for this pulsar. Among the new precise measurements are the pulsar mass $M_p=1.418\pm 0.044$ M$_{\odot}$, distance $D=156.96 \pm 0.11$ pc, and orbital inclination angle $i=137.506 \pm 0.016^\circ$, which can be used to inform the X-ray pulse profile models inferred from NICER observations. We demonstrate that these results are consistent between multiple data sets from the Parkes Pulsar Timing Array (PPTA), each modelled with different noise assumptions. Using the longest available PPTA data set, we measure an apparent second derivative of the pulsar spin frequency and discuss how this can be explained either by kinematic effects due to the proper motion and radial velocity of the pulsar, or excess low-frequency noise such as a gravitational-wave background. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: 13 pages, 3 figures, accepted for publication in Astrophysical Journal Letters

Increasing the detectability of long-period and nulling pulsars in next-generation pulsar surveys

Authors: G. Grover, N. D. R. Bhat, S. J. McSweeney

Abstract: Recent discoveries of multiple long-period pulsars (periods ${\sim}10\,$s or larger) are starting to challenge the conventional notion that coherent radio emission cannot be produced by objects that are below the many theorised death lines. Many of the past pulsar surveys and software have been prone to selection effects that restricted their sensitivities towards long-period and sporadically-emit… ▽ More Recent discoveries of multiple long-period pulsars (periods ${\sim}10\,$s or larger) are starting to challenge the conventional notion that coherent radio emission cannot be produced by objects that are below the many theorised death lines. Many of the past pulsar surveys and software have been prone to selection effects that restricted their sensitivities towards long-period and sporadically-emitting objects. Pulsar surveys using new-generation low-frequency facilities are starting to employ longer dwell times, which makes them significantly more sensitive in detecting long-period or nulling pulsars. There have also been software advancements to aid more sensitive searches towards long-period objects. Furthermore, recent discoveries suggest that nulling may be a key aspect of the long-period pulsar population. We simulate both long-period and nulling pulsar signals, using the Southern-sky MWA Rapid Two-meter (SMART) survey data as reference, and explore the detection efficacy of popular search methods such as the fast Fourier transform (FFT), fast-folding algorithm (FFA) and single pulse search (SPS). For FFT-based search and SPS, we make use of the PRESTO implementation, and for FFA we use RIPTIDE. We find RIPTIDE's FFA to be more sensitive; however, it is also the slowest algorithm. PRESTO's FFT, although faster than others, also shows some unexpected inaccuracies in detection properties. SPS is highly sensitive to long-period and nulling signals, but only for pulses with high intrinsic signal-to-noise ratios. We use these findings to inform current and future pulsar surveys that aim to uncover a large population of long-period or nulling objects and comment on how to make optimal use of these methods in unison. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 13 pages, 4 figures, 2 tables. Accepted for publication by PASA

Discovery and follow-up of a quasiperiodically nulling and sub-pulse drifting pulsar with the Murchison Widefield Array

Authors: G. Grover, N. D. R. Bhat, S. McSweeney, C. P. Lee, B. W. Meyers, C. M. Tan, S. S. Kudale

Abstract: The phenomenon of pulsar nulling, where pulsars temporarily and stochastically cease their radio emission, is thought to be indicative of a `dying' pulsar, where radio emission ceases entirely. Here we report the discovery of a long-period pulsar, PSR J0452-3418, from the ongoing Southern-sky MWA Rapid Two-meter (SMART) pulsar survey. The pulsar has a rotation period of ${\sim}$1.67\,s and a dispe… ▽ More The phenomenon of pulsar nulling, where pulsars temporarily and stochastically cease their radio emission, is thought to be indicative of a `dying' pulsar, where radio emission ceases entirely. Here we report the discovery of a long-period pulsar, PSR J0452-3418, from the ongoing Southern-sky MWA Rapid Two-meter (SMART) pulsar survey. The pulsar has a rotation period of ${\sim}$1.67\,s and a dispersion measure of 19.8\,\dmu, and it exhibits both quasi-periodic nulling and sub-pulse drifting. Periodic nulling is uncommon, only reported in $<1$\% of the pulsar population, with even a smaller fraction showing periodic nulling and sub-pulse drifting. We describe the discovery and follow-up of the pulsar, including a positional determination using high-resolution imaging with the upgraded Giant Metrewave Radio Telescope (uGMRT), initial timing analysis using the combination of MWA and uGMRT data, and detailed characterisation of the nulling and drifting properties in the MWA's frequency band (140-170\,MHz). Our analysis suggests a nulling fraction of 34$\pm6$\% and a nulling periodicity of 42$^{+1.5}_{-1.3}$ pulses. We measure the phase ($P_2$) and time modulation ($P_3$) caused by the sub-pulse drifting, with an average $P_2$ of 7.1$^{+26.3}_{-3.1}$ degrees and a $P_3$ of 4.8$^{+1.5}_{-0.9}$ pulses. We compare and contrast the observed properties with those of other pulsars that exhibit sub-pulse drifting and quasi-periodic nulling phenomena, and find that the majority of these objects tend to be in the `death valley' in the period-period derivative ($P$-$\dot{P}$) diagram. We also discuss some broader implications for pulsar emission physics and the detectability of similar objects using next-generation pulsar surveys. △ Less

Submitted 28 May, 2024; v1 submitted 26 May, 2024; originally announced May 2024.

Comments: 16 pages, 9 Figures, 4 Tables, Accepted for ApJ***

High-Time Resolution GPU Imager for FRB searches at low radio frequencies

Authors: M. Sokolowski, G. Aniruddha, C. Di Pietrantonio, C. Harris, D. C. Price, S. McSweeney, R. B. Wayth, N. D. R. Bhat

Abstract: Fast Radio Bursts (FRBs) are millisecond dispersed radio pulses of predominately extra-galactic origin. Although originally discovered at GHz frequencies, most FRBs have been detected between 400 to 800 MHz. Nevertheless, only a handful of FRBs were detected below 400 MHz. Searching for FRBs at low frequencies is computationally challenging due to increased dispersive delay that must be accounted… ▽ More Fast Radio Bursts (FRBs) are millisecond dispersed radio pulses of predominately extra-galactic origin. Although originally discovered at GHz frequencies, most FRBs have been detected between 400 to 800 MHz. Nevertheless, only a handful of FRBs were detected below 400 MHz. Searching for FRBs at low frequencies is computationally challenging due to increased dispersive delay that must be accounted for. However, the wide field of view (FoV) of low-frequency telescopes - such as the the Murchison Widefield Array (MWA), and prototype stations of the low-frequency Square Kilometre Array (SKA-Low) - makes them promising instruments to open a low-frequency window on FRB event rates, and constrain FRB emission models. The standard approach, inherited from high-frequencies, is to form multiple tied-array beams to tessellate the entire FoV and perform the search on the resulting time series. This approach, however, may not be optimal for low-frequency interferometers due to their large FoVs and high spatial resolutions leading to a large number of beams. Consequently, there are regions of parameter space in terms of number of antennas and resolution elements (pixels) where interferometric imaging is computationally more efficient. Here we present a new high-time resolution imager BLINK implemented on modern Graphical Processing Units (GPUs) and intended for radio astronomy. The main goal for this imager is to become part of a fully GPU-accelerated FRB search pipeline. We describe the imager and present its verification on real and simulated data processed to form all-sky and widefield images from the MWA and prototype SKA-Low stations. We also present and compare benchmarks of the GPU and CPU code executed on laptops, desktop computers, and Australian supercomputers. The code is publicly available at https://github.com/PaCER-BLINK-Project/imager and can be applied to data from any radio telescope. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: 19 pages, 9 figures, 6 tables. Accepted for publication in PASA

Constraining possible $γ$-ray burst emission from GW230529 using Swift-BAT and Fermi-GBM

Authors: Samuele Ronchini, Suman Bala, Joshua Wood, James Delaunay, Simone Dichiara, Jamie A. Kennea, Tyler Parsotan, Gayathri Raman, Aaron Tohuvavohu, Naresh Adhikari, Narayana P. Bhat, Sylvia Biscoveanu, Elisabetta Bissaldi, Eric Burns, Sergio Campana, Koustav Chandra, William H. Cleveland, Sarah Dalessi, Massimiliano De Pasquale, Juan García-Bellido, Claudio Gasbarra, Misty M. Giles, Ish Gupta, Dieter Hartmann, Boyan A. Hristov , et al. (13 additional authors not shown)

Abstract: GW230529 is the first compact binary coalescence detected by the LIGO-Virgo-KAGRA collaboration with at least one component mass confidently in the lower mass-gap, corresponding to the range 3-5$M_{\odot}$. If interpreted as a neutron star-black hole merger, this event has the most symmetric mass ratio detected so far and therefore has a relatively high probability of producing electromagnetic (EM… ▽ More GW230529 is the first compact binary coalescence detected by the LIGO-Virgo-KAGRA collaboration with at least one component mass confidently in the lower mass-gap, corresponding to the range 3-5$M_{\odot}$. If interpreted as a neutron star-black hole merger, this event has the most symmetric mass ratio detected so far and therefore has a relatively high probability of producing electromagnetic (EM) emission. However, no EM counterpart has been reported. At the merger time $t_0$, Swift-BAT and Fermi-GBM together covered 100$\%$ of the sky. Performing a targeted search in a time window $[t_0-20 \text{s},t_0+20 \text{s}]$, we report no detection by the Swift-BAT and the Fermi-GBM instruments. Combining the position-dependent $γ-$ray flux upper limits and the gravitational-wave posterior distribution of luminosity distance, sky localization and inclination angle of the binary, we derive constraints on the characteristic luminosity and structure of the jet possibly launched during the merger. Assuming a top-hat jet structure, we exclude at 90$\%$ credibility the presence of a jet which has at the same time an on-axis isotropic luminosity $\gtrsim 10^{48}$ erg s$^{-1}$, in the bolometric band 1 keV-10 MeV, and a jet opening angle $\gtrsim 15$ deg. Similar constraints are derived testing other assumptions about the jet structure profile. Excluding GRB 170817A, the luminosity upper limits derived here are below the luminosity of any GRB observed so far. △ Less

Submitted 17 May, 2024; originally announced May 2024.

Comments: 18 pages, 1 table, 11 figures

A commensal Fast Radio Burst search pipeline for the Murchison Widefield Array

Authors: M. Sokolowski, I. S. Morrison, D. Price, G. Sleap, B. Crosse, A. Williams, L. Williams, C. James, B. W. Meyers, S. McSweeney, N. D. R. Bhat, G. Anderson

Abstract: We present a demonstration version of a commensal pipeline for Fast Radio Burst (FRB) searches using a real-time incoherent beam from the Murchison Widefield Array (MWA). The main science target of the pipeline are bright nearby FRBs from the local Universe which are the best candidates to probe FRB progenitors and understand physical mechanisms powering these extremely energetic events. The new M… ▽ More We present a demonstration version of a commensal pipeline for Fast Radio Burst (FRB) searches using a real-time incoherent beam from the Murchison Widefield Array (MWA). The main science target of the pipeline are bright nearby FRBs from the local Universe which are the best candidates to probe FRB progenitors and understand physical mechanisms powering these extremely energetic events. The new MWA beamformer, known as the "MWAX multibeam beamformer", can form multiple incoherent and coherent beams commensally to any on-going MWA observations. One of the beams is currently used for FRB searches (tested in 10 kHz frequency resolution and time resolutions between 0.1 and 100 ms). A second beam is used for a Search for Extraterrestrial Intelligence (SETI). This paper focuses on the FRB search pipeline and its verification on selected known bright pulsars. The pipeline uses the FREDDA implementation of the Fast Dispersion Measure Transform algorithm (FDMT) for single pulse searches. Initially, it was tested during standard MWA observations, and more recently using dedicated observations of selected 11 bright pulsars. The pulsar PSR J0835-4510 (aka Vela) has been routinely used as the primary probe of the data quality because its folded profile was always detected in the frequency band 200 - 230 MHz with typical SNR >10. Similarly, the low DM pulsar PSR B0950+08 was always detected in folded profile in the frequency band 140 - 170 MHz, and so far has been the only object for which single pulses were detected. We present the estimated sensitivity of the search in the currently limited observing bandwidth of a single MWA coarse channel (1.28 MHz) and for the upgraded, future system with 12.8 MHz (10 channels) of bandwidth. Based on expected sensitivity and existing FRB rate measurements, we estimate an expected number of FRB detections to be between a few and a few tens per year. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: 22 pages, 12 figures, 4 tables. Accepted for publication in PASA

The High Time Resolution Universe Pulsar Survey -- XVIII. The reprocessing of the HTRU-S Low Lat survey around the Galactic centre using a Fast Folding Algorithm pipeline for accelerated pulsars

Authors: J. Wongphechauxsorn, D. J. Champion, M. Bailes, V. Balakrishnan, E. D. Barr, M. C. i Bernadich, N. D. R. Bhat, M. Burgay, A. D. Cameron, W. Chen, C. M. L. Flynn, A. Jameson, S. Johnston, M. J. Keith, M. Kramer, C. Ng, A. Possenti, R. Sengar, R. M. Shannon, B. Stappers, W. van Straten

Abstract: The HTRU-S Low Latitude survey data within 1$^{\circ}$of the Galactic Centre (GC) were searched for pulsars using the Fast Folding Algorithm (FFA). Unlike traditional Fast Fourier Transform (FFT) pipelines, the FFA optimally folds the data for all possible periods over a given range, which is particularly advantageous for pulsars with low-duty cycle. For the first time, a search over acceleration… ▽ More The HTRU-S Low Latitude survey data within 1$^{\circ}$of the Galactic Centre (GC) were searched for pulsars using the Fast Folding Algorithm (FFA). Unlike traditional Fast Fourier Transform (FFT) pipelines, the FFA optimally folds the data for all possible periods over a given range, which is particularly advantageous for pulsars with low-duty cycle. For the first time, a search over acceleration was included in the FFA to improve its sensitivity to binary pulsars. The steps in dispersion measure (DM) and acceleration were optimised, resulting in a reduction of the number of trials by 86 per cent. This was achieved over a search period range from 0.6-s to 432-s, i.e. 10 per cent of the observation time (4320s), with a maximum DM of 4000 pc cm$^{-3}$ and an acceleration range of $\pm 128$m s$^{-2}$. The search resulted in the re-detections of four known pulsars, including a pulsar which was missed in previous FFT processing of this survey. This result indicates that the FFA pipeline is more sensitive than the FFT pipeline used in the previous processing of the survey within our parameter range. Additionally, we discovered a 1.89-s pulsar, PSR J1746-2829, with a large DM, located~0.5 from the GC. Follow-up observations revealed that this pulsar has a relatively flat spectrum($α=-0.9\pm0.1$) and has a period derivative of $\sim1.3\times10^{-12}$ s s$^{-1}$, implying a surface magnetic field of $\sim5.2\times10^{13}$ G and a characteristic age of $\sim23000$ yr. While the period, spectral index, and surface magnetic field strength are similar to many radio magnetars, other characteristics such as high linear polarization are absent. △ Less

Submitted 21 October, 2023; originally announced October 2023.

Comments: 12 pages, 6 figures, 4 tables, Accepted for publication on Monthly Notices of the Royal Astronomical Society

Resolving moving heliospheric structures using interplanetary scintillation observations with the Murchison Widefield Array

Authors: A. Waszewski, J. S. Morgan, R. Chhetri, R. Ekers, M. C. M. Cheung, N. D. R Bhat, M. Johnston-Hollitt

Abstract: We have conducted a blind search in 49 consecutive days of interplanetary scintillation observations made by the Murchison Widefield Array from mid-2019, with overlapping daily observations approximately East and South-East of the Sun at an elongation of $\sim$30 degrees and a field of view of 30 degrees. These observations detect an unprecedented density of sources. In spite of these observations… ▽ More We have conducted a blind search in 49 consecutive days of interplanetary scintillation observations made by the Murchison Widefield Array from mid-2019, with overlapping daily observations approximately East and South-East of the Sun at an elongation of $\sim$30 degrees and a field of view of 30 degrees. These observations detect an unprecedented density of sources. In spite of these observations being taken at sunspot minimum, this search has revealed several interesting transitory features characterised by elevated scintillation levels. One solar wind enhancement is captured in two observations several hours apart, allowing its radial movement away from the Sun to be measured. We present here a methodology for measuring the plane-of-sky velocity for the moving heliospheric structure. The plane-of-sky velocity was inferred as $0.66\pm0.147\,^{\text{o}}\text{hr}^{-1}$, or $480\pm106\,\text{km}\,\text{s}^{-1}$ assuming a distance of 1AU. After cross-referencing our observed structure with multiple catalogues of heliospheric events, we propose that the likely source of our observed structure is a stream-interaction region originating from a low-latitude coronal hole. This work demonstrates the power of widefield interplanetary scintillation observations to capture detailed features in the heliosphere which are otherwise unresolvable and go undetected. △ Less

Submitted 19 September, 2023; originally announced September 2023.

Comments: 14 pages, 5 figures

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

PSR J0026-1955: A curious case of evolutionary subpulse drifting and nulling

Authors: Parul Janagal, Samuel J. McSweeney, Manoneeta Chakraborty, N. D. Ramesh Bhat

Abstract: PSR J0026-1955 was independently discovered by the Murchison Widefield Array (MWA) recently. The pulsar exhibits subpulse drifting, where the radio emission from a pulsar appears to drift in spin phase within the main pulse profile, and nulling, where the emission ceases briefly. The pulsar showcases a curious case of drift rate evolution as it exhibits rapid changes between the drift modes and a… ▽ More PSR J0026-1955 was independently discovered by the Murchison Widefield Array (MWA) recently. The pulsar exhibits subpulse drifting, where the radio emission from a pulsar appears to drift in spin phase within the main pulse profile, and nulling, where the emission ceases briefly. The pulsar showcases a curious case of drift rate evolution as it exhibits rapid changes between the drift modes and a gradual evolution in the drift rate within a mode. Here we report new analysis and results from observations of J0026-1955 made with the upgraded Giant Meterwave Radio Telescope (uGMRT) at 300-500 MHz. We identify two distinct subpulse drifting modes: A and B, with mode A sub-categorised into A0, A1, and A2, depending upon the drift rate evolutionary behaviour. Additionally, the pulsar exhibits short and long nulls, with an estimated overall nulling fraction of ~58%, which is lower than the previously reported value. Our results also provide evidence of subpulse memory across nulls and a consistent behaviour where mode A2 is often followed by a null. We investigate the drift rate modulations of J0026-1955 and put forward two different models to explain the observed drifting behaviour. We suggest that either a change in polar gap screening or a slow relaxation in the spark configuration could possibly drive the evolution in drift rates. J0026-1955 belongs to a rare subset of pulsars which exhibit subpulse drifting, nulling, mode changing, and drift rate evolution. It is, therefore, an ideal test bed for carousel models and to uncover the intricacies of pulsar emission physics. △ Less

Submitted 5 July, 2023; originally announced July 2023.

Comments: 14 pages, 10 figures, 1 table. Accepted for publication in MNRAS

The Parkes Pulsar Timing Array Third Data Release

Authors: Andrew Zic, Daniel J. Reardon, Agastya Kapur, George Hobbs, Rami Mandow, Małgorzata Curyło, Ryan M. Shannon, Jacob Askew, Matthew Bailes, N. D. Ramesh Bhat, Andrew Cameron, Zu-Cheng Chen, Shi Dai, Valentina Di Marco, Yi Feng, Matthew Kerr, Atharva Kulkarni, Marcus E. Lower, Rui Luo, Richard N. Manchester, Matthew T. Miles, Rowina S. Nathan, Stefan Osłowski, Axl F. Rogers, Christopher J. Russell , et al. (9 additional authors not shown)

Abstract: We present the third data release from the Parkes Pulsar Timing Array (PPTA) project. The release contains observations of 32 pulsars obtained using the 64-m Parkes "Murriyang" radio telescope. The data span is up to 18 years with a typical cadence of 3 weeks. This data release is formed by combining an updated version of our second data release with $\sim 3$ years of more recent data primarily ob… ▽ More We present the third data release from the Parkes Pulsar Timing Array (PPTA) project. The release contains observations of 32 pulsars obtained using the 64-m Parkes "Murriyang" radio telescope. The data span is up to 18 years with a typical cadence of 3 weeks. This data release is formed by combining an updated version of our second data release with $\sim 3$ years of more recent data primarily obtained using an ultra-wide-bandwidth receiver system that operates between 704 and 4032 MHz. We provide calibrated pulse profiles, flux-density dynamic spectra, pulse times of arrival, and initial pulsar timing models. We describe methods for processing such wide-bandwidth observations, and compare this data release with our previous release. △ Less

Submitted 17 October, 2023; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 15 pages, 6 figures. Accepted for publication in PASA

The gravitational-wave background null hypothesis: Characterizing noise in millisecond pulsar arrival times with the Parkes Pulsar Timing Array

Authors: Daniel J. Reardon, Andrew Zic, Ryan M. Shannon, Valentina Di Marco, George B. Hobbs, Agastya Kapur, Marcus E. Lower, Rami Mandow, Hannah Middleton, Matthew T. Miles, Axl F. Rogers, Jacob Askew, Matthew Bailes, N. D. Ramesh Bhat, Andrew Cameron, Matthew Kerr, Atharva Kulkarni, Richard N. Manchester, Rowina S. Nathan, Christopher J. Russell, Stefan Osłowski, Xing-Jiang Zhu

Abstract: The noise in millisecond pulsar (MSP) timing data can include contributions from observing instruments, the interstellar medium, the solar wind, solar system ephemeris errors, and the pulsars themselves. The noise environment must be accurately characterized in order to form the null hypothesis from which signal models can be compared, including the signature induced by nanohertz-frequency gravita… ▽ More The noise in millisecond pulsar (MSP) timing data can include contributions from observing instruments, the interstellar medium, the solar wind, solar system ephemeris errors, and the pulsars themselves. The noise environment must be accurately characterized in order to form the null hypothesis from which signal models can be compared, including the signature induced by nanohertz-frequency gravitational waves (GWs). Here we describe the noise models developed for each of the MSPs in the Parkes Pulsar Timing Array (PPTA) third data release, which have been used as the basis of a search for the isotropic stochastic GW background. We model pulsar spin noise, dispersion measure variations, scattering variations, events in the pulsar magnetospheres, solar wind variability, and instrumental effects. We also search for new timing model parameters and detected Shapiro delays in PSR~J0614$-$3329 and PSR~J1902$-$5105. The noise and timing models are validated by testing the normalized and whitened timing residuals for Gaussianity and residual correlations with time. We demonstrate that the choice of noise models significantly affects the inferred properties of a common-spectrum process. Using our detailed models, the recovered common-spectrum noise in the PPTA is consistent with a power law with a spectral index of $γ=13/3$, the value predicted for a stochastic GW background from a population of supermassive black hole binaries driven solely by GW emission. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 18 pages, 10 figures. Accepted for publication in ApJL

Search for an isotropic gravitational-wave background with the Parkes Pulsar Timing Array

Authors: Daniel J. Reardon, Andrew Zic, Ryan M. Shannon, George B. Hobbs, Matthew Bailes, Valentina Di Marco, Agastya Kapur, Axl F. Rogers, Eric Thrane, Jacob Askew, N. D. Ramesh Bhat, Andrew Cameron, Małgorzata Curyło, William A. Coles, Shi Dai, Boris Goncharov, Matthew Kerr, Atharva Kulkarni, Yuri Levin, Marcus E. Lower, Richard N. Manchester, Rami Mandow, Matthew T. Miles, Rowina S. Nathan, Stefan Osłowski , et al. (4 additional authors not shown)

Abstract: Pulsar timing arrays aim to detect nanohertz-frequency gravitational waves (GWs). A background of GWs modulates pulsar arrival times and manifests as a stochastic process, common to all pulsars, with a signature spatial correlation. Here we describe a search for an isotropic stochastic gravitational-wave background (GWB) using observations of 30 millisecond pulsars from the third data release of t… ▽ More Pulsar timing arrays aim to detect nanohertz-frequency gravitational waves (GWs). A background of GWs modulates pulsar arrival times and manifests as a stochastic process, common to all pulsars, with a signature spatial correlation. Here we describe a search for an isotropic stochastic gravitational-wave background (GWB) using observations of 30 millisecond pulsars from the third data release of the Parkes Pulsar Timing Array (PPTA), which spans 18 years. Using current Bayesian inference techniques we recover and characterize a common-spectrum noise process. Represented as a strain spectrum $h_c = A(f/1 {\rm yr}^{-1})^α$, we measure $A=3.1^{+1.3}_{-0.9} \times 10^{-15}$ and $α=-0.45 \pm 0.20$ respectively (median and 68% credible interval). For a spectral index of $α=-2/3$, corresponding to an isotropic background of GWs radiated by inspiraling supermassive black hole binaries, we recover an amplitude of $A=2.04^{+0.25}_{-0.22} \times 10^{-15}$. However, we demonstrate that the apparent signal strength is time-dependent, as the first half of our data set can be used to place an upper limit on $A$ that is in tension with the inferred common-spectrum amplitude using the complete data set. We search for spatial correlations in the observations by hierarchically analyzing individual pulsar pairs, which also allows for significance validation through randomizing pulsar positions on the sky. For a process with $α=-2/3$, we measure spatial correlations consistent with a GWB, with an estimated false-alarm probability of $p \lesssim 0.02$ (approx. $2σ$). The long timing baselines of the PPTA and the access to southern pulsars will continue to play an important role in the International Pulsar Timing Array. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 19 pages, 10 figures, Accepted for publication in ApJL

Single-pulse analysis and average emission characteristics of PSR J1820-0427 from observations made with the MWA and uGMRT

Authors: Parul Janagal, Manoneeta Chakraborty, N. D. Ramesh Bhat, Samuel J. McSweeney, Susmita Sett

Abstract: We have studied the pulse-to-pulse variability in PSR J1820--0427 and its frequency dependence using high-quality, wide-band observations made from the upgraded Giant Metrewave Radio Telescope (uGMRT; 300-750 MHz) and the Murchison Widefield Array ($\sim$170-200 MHz). The low-frequency data reveal a previously unreported feature in the average profile (at 185 MHz) after accounting for the effects… ▽ More We have studied the pulse-to-pulse variability in PSR J1820--0427 and its frequency dependence using high-quality, wide-band observations made from the upgraded Giant Metrewave Radio Telescope (uGMRT; 300-750 MHz) and the Murchison Widefield Array ($\sim$170-200 MHz). The low-frequency data reveal a previously unreported feature in the average profile (at 185 MHz) after accounting for the effects of temporal broadening arising from multi-path scattering due to the Interstellar Medium (ISM). We advance a new method for flux density calibration of beamformed data from the uGMRT and use it to measure the single pulse flux densities across the uGMRT band. Combined with previously published measurements, these flux densities are best fit with a power-law spectrum with a low-frequency turnover. We also use calibrated flux densities to explore the relationship between pulse-to-pulse variability and the spectral index of individual pulses. Our analysis reveals a large scatter in the single-pulse spectral indices and a general tendency for brighter pulses to show a steepening of the spectral index. We also examine the frequency-dependence of the pulse-fluence distribution and its relation to the Stochastic Growth Theory. △ Less

Submitted 13 June, 2023; originally announced June 2023.

Comments: 13 pages, 9 figures, 2 tables. Accepted for publication in MNRAS

Testing the circularity of PSR B0818-41's carousel

Authors: Samuel J. McSweeney, Lisa Smith, N. D. Ramesh Bhat, Geoff Wright

Abstract: The phenomenon of sub-pulse drifting is an important single-pulse phenomenon that can potentially provide important insights into the elusive radio emission mechanism in pulsars. We analyze the frequency behaviour of the single pulses of B0818-41, observed from 300 to 500 MHz (Band 3 of the uGMRT), and compare it to the evolution of the average profile to place constraints on the geometry of the p… ▽ More The phenomenon of sub-pulse drifting is an important single-pulse phenomenon that can potentially provide important insights into the elusive radio emission mechanism in pulsars. We analyze the frequency behaviour of the single pulses of B0818-41, observed from 300 to 500 MHz (Band 3 of the uGMRT), and compare it to the evolution of the average profile to place constraints on the geometry of the pulsar's emission beam. We show that a circular carousel of discrete beamlets, where each beamlet has radial symmetry, is not consistent with the observed behaviour, and describe an alternative, consistent range of possible elliptical carousel geometries. We also combine the uGMRT data with some archival MWA observations and several other published profiles to characterize the profile evolution across a frequency range spanning ~170 MHz to ~1.4 GHz. △ Less

Submitted 8 June, 2023; originally announced June 2023.

Extreme Variability in a Long Duration Gamma-ray Burst Associated with a Kilonova

Authors: P. Veres, P. N. Bhat, E. Burns, R. Hamburg, N. Fraija, D. Kocevski, R. Preece, S. Poolakkil, N. Christensen, M. A. Bizouard, T. Dal Canton, S. Bala, E. Bissaldi, M. S. Briggs, W. Cleveland, A. Goldstein, B. A. Hristov, C. M. Hui, S. Lesage, B. Mailyan, O. J. Roberts, C. A. Wilson-Hodge

Abstract: The recent discovery of a kilonova from the long duration gamma-ray burst, GRB 211211A, challenges classification schemes based on temporal information alone. Gamma-ray properties of GRB 211211A reveal an extreme event, which stands out among both short and long GRBs. We find very short variations (few ms) in the lightcurve of GRB 211211A and estimate ~1000 for the Lorentz factor of the outflow. W… ▽ More The recent discovery of a kilonova from the long duration gamma-ray burst, GRB 211211A, challenges classification schemes based on temporal information alone. Gamma-ray properties of GRB 211211A reveal an extreme event, which stands out among both short and long GRBs. We find very short variations (few ms) in the lightcurve of GRB 211211A and estimate ~1000 for the Lorentz factor of the outflow. We discuss the relevance of the short variations in identifying similar long GRBs resulting from compact mergers. Our findings indicate that in future gravitational wave follow-up campaigns, some long duration GRBs should be treated as possible strong gravitational wave counterparts. △ Less

Submitted 20 May, 2023; originally announced May 2023.

Comments: 10 pages 5 figures, submitted to AAS journals

Fermi-GBM Discovery of GRB 221009A: An Extraordinarily Bright GRB from Onset to Afterglow

Authors: S. Lesage, P. Veres, M. S. Briggs, A. Goldstein, D. Kocevski, E. Burns, C. A. Wilson-Hodge, P. N. Bhat, D. Huppenkothen, C. L. Fryer, R. Hamburg, J. Racusin, E. Bissaldi, W. H. Cleveland, S. Dalessi, C. Fletcher, M. M. Giles, B. A. Hristov, C. M. Hui, B. Mailyan, C. Malacaria, S. Poolakkil, O. J. Roberts, A. von Kienlin, J. Wood , et al. (115 additional authors not shown)

Abstract: We report the discovery of GRB 221009A, the highest flux gamma-ray burst ever observed by the Fermi Gamma-ray Burst Monitor (GBM). This GRB has continuous prompt emission lasting more than 600 seconds which smoothly transitions to afterglow visible in the GBM energy range (8 keV--40 MeV), and total energetics higher than any other burst in the GBM sample. By using a variety of new and existing ana… ▽ More We report the discovery of GRB 221009A, the highest flux gamma-ray burst ever observed by the Fermi Gamma-ray Burst Monitor (GBM). This GRB has continuous prompt emission lasting more than 600 seconds which smoothly transitions to afterglow visible in the GBM energy range (8 keV--40 MeV), and total energetics higher than any other burst in the GBM sample. By using a variety of new and existing analysis techniques we probe the spectral and temporal evolution of GRB 221009A. We find no emission prior to the GBM trigger time (t0; 2022 October 9 at 13:16:59.99 UTC), indicating that this is the time of prompt emission onset. The triggering pulse exhibits distinct spectral and temporal properties suggestive of the thermal, photospheric emission of shock-breakout, with significant emission up to $\sim$15 MeV. We characterize the onset of external shock at t0+600 s and find evidence of a plateau region in the early-afterglow phase which transitions to a slope consistent with Swift-XRT afterglow measurements. We place the total energetics of GRB 221009A in context with the rest of the GBM sample and find that this GRB has the highest total isotropic-equivalent energy ($\textrm{E}_{γ,\textrm{iso}}=1.0\times10^{55}$ erg) and second highest isotropic-equivalent luminosity ($\textrm{L}_{γ,\textrm{iso}}=9.9\times10^{53}$ erg/s) based on redshift of z = 0.151. These extreme energetics are what allowed us to observe the continuously emitting central engine of GBM from the beginning of the prompt emission phase through the onset of early afterglow. △ Less

Submitted 12 July, 2023; v1 submitted 24 March, 2023; originally announced March 2023.

Comments: 26 pages 7 figures - accepted for publication in ApJL

Searching for continuous Gravitational Waves in the second data release of the International Pulsar Timing Array

Authors: M. Falxa, S. Babak, P. T. Baker, B. Bécsy, A. Chalumeau, S. Chen, Z. Chen, N. J. Cornish, L. Guillemot, J. S. Hazboun, C. M. F. Mingarelli, A. Parthasarathy, A. Petiteau, N. S. Pol, A. Sesana, S. B. Spolaor, S. R. Taylor, G. Theureau, M. Vallisneri, S. J. Vigeland, C. A. Witt, X. Zhu, J. Antoniadis, Z. Arzoumanian, M. Bailes , et al. (102 additional authors not shown)

Abstract: The International Pulsar Timing Array 2nd data release is the combination of datasets from worldwide collaborations. In this study, we search for continuous waves: gravitational wave signals produced by individual supermassive black hole binaries in the local universe. We consider binaries on circular orbits and neglect the evolution of orbital frequency over the observational span. We find no evi… ▽ More The International Pulsar Timing Array 2nd data release is the combination of datasets from worldwide collaborations. In this study, we search for continuous waves: gravitational wave signals produced by individual supermassive black hole binaries in the local universe. We consider binaries on circular orbits and neglect the evolution of orbital frequency over the observational span. We find no evidence for such signals and set sky averaged 95% upper limits on their amplitude h 95 . The most sensitive frequency is 10nHz with h 95 = 9.1 10-15 . We achieved the best upper limit to date at low and high frequencies of the PTA band thanks to improved effective cadence of observations. In our analysis, we have taken into account the recently discovered common red noise process, which has an impact at low frequencies. We also find that the peculiar noise features present in some pulsars data must be taken into account to reduce the false alarm. We show that using custom noise models is essential in searching for continuous gravitational wave signals and setting the upper limit. △ Less

Submitted 19 March, 2023; originally announced March 2023.

The Southern-sky MWA Rapid Two-metre (SMART) pulsar survey -- II. Survey status, pulsar census, and first pulsar discoveries

Authors: N. D. R. Bhat, N. A. Swainston, S. J. McSweeney, M. Xue, B. W. Meyers, S. Kudale, S. Dai, S. E. Tremblay, W. van Straten, R. M. Shannon, K. R. Smith, M. Sokolowski, S. M. Ord, G. Sleap, A. Williams, P. J. Hancock, R. Lange, J. Tocknell, M. Johnston-Hollitt, D. L. Kaplan, S. J. Tingay, M. Walker

Abstract: In Paper I, we presented an overview of the Southern-sky MWA Rapid Two-metre (SMART) survey, including the survey design and search pipeline. While the combination of MWA's large field-of-view and the voltage capture system brings a survey speed of ~450 square degrees per hour, the survey progression relies on the availability of compact configuration of the Phase II array. Over the past few years… ▽ More In Paper I, we presented an overview of the Southern-sky MWA Rapid Two-metre (SMART) survey, including the survey design and search pipeline. While the combination of MWA's large field-of-view and the voltage capture system brings a survey speed of ~450 square degrees per hour, the survey progression relies on the availability of compact configuration of the Phase II array. Over the past few years, by taking advantage of multiple windows of opportunity when the compact configuration was available, we have advanced the survey to 75% completion. To date, about 10% of the data collected thus far have been processed for a first-pass search, where 10 minutes of observation is processed for dispersion measures out to 250 ${\rm pc\,cm^{-3}}$, to realise a shallow survey for long-period pulsars. The ongoing analysis has led to two new pulsar discoveries, as well as an independent discovery and a rediscovery of a previously incorrectly characterised pulsar, all from ~3% of the data for which candidate scrutiny is completed. Here we describe the strategies for further detailed follow-up including improved sky localisation and convergence to timing solution, and illustrate them using example pulsar discoveries. The processing has also led to re-detection of 120 pulsars in the SMART observing band, bringing the total number of pulsars detected to date with the MWA to 180, and these are used to assess the search sensitivity of current processing pipelines. The planned second-pass (deep survey) processing is expected to yield a three-fold increase in sensitivity for long-period pulsars, and a substantial improvement to millisecond pulsars by adopting optimal de-dispersion plans. The SMART survey will complement the highly successful Parkes High Time Resolution Universe survey at 1.2-1.5 GHz, and inform future large survey efforts such as those planned with the low-frequency Square Kilometre Array. △ Less

Submitted 23 February, 2023; originally announced February 2023.

Comments: 22 pages, 9 figures, 7 tables, Accepted for publication in PASA

The Southern-sky MWA Rapid Two-metre (SMART) pulsar survey -- I. Survey design and processing pipeline

Authors: N. D. R. Bhat, N. A. Swainston, S. J. McSweeney, M. Xue, B. W. Meyers, S. Kudale, S. Dai, S. E. Tremblay, W. van Straten, R. M. Shannon, K. R. Smith, M. Sokolowski, S. M. Ord, G. Sleap, A. Williams, P. J. Hancock, R. Lange, J. Tocknell, M. Johnston-Hollitt, D. L. Kaplan, S. J. Tingay, M. Walker

Abstract: We present an overview of the Southern-sky MWA Rapid Two-metre (SMART) pulsar survey that exploits the MWA's large field of view and voltage capture system to survey the sky south of 30 degree in declination for pulsars and fast transients in the 140-170 MHz band. The survey is enabled by the advent of the Phase II MWA's compact configuration, which offers an enormous efficiency in beam-forming an… ▽ More We present an overview of the Southern-sky MWA Rapid Two-metre (SMART) pulsar survey that exploits the MWA's large field of view and voltage capture system to survey the sky south of 30 degree in declination for pulsars and fast transients in the 140-170 MHz band. The survey is enabled by the advent of the Phase II MWA's compact configuration, which offers an enormous efficiency in beam-forming and processing costs, thereby making an all-sky survey of this magnitude tractable with the MWA. Even with the long dwell times of the survey (4800 s), data collection can be completed in < 100 hours of telescope time, while still retaining the ability to reach a limiting sensitivity of ~2-3 mJy. Each observation is processed to generate ~5000-8000 tied-array beams that tessellate the full ~610 square degree field of view, which are then processed to search for pulsars. The voltage-capture recording allows a multitude of post hoc processing options including the reprocessing of data for higher time resolution. Due to the substantial computational cost in pulsar searches at low frequencies, processing is undertaken in multiple passes: in the first pass, a shallow survey is performed, where 10 minutes of each observation is processed, reaching about one-third of the full search sensitivity. Here we present the system overview and initial results. Further details including first pulsar discoveries and a census of low-frequency detections are presented in a companion paper. Future plans include deeper searches to reach the full sensitivity and acceleration searches to target binary and millisecond pulsars. Simulation analysis forecasts ~300 new pulsars upon the completion of full processing. The SMART survey will also generate a complete digital record of the low-frequency sky, which will serve as a valuable reference for future pulsar searches planned with the low-frequency Square Kilometre Array. △ Less

Submitted 23 February, 2023; originally announced February 2023.

Comments: 22 pages, 12 figures, 2 tables, Accepted for publication in PASA

Image-based searches for pulsar candidates using MWA VCS data

Authors: S. Sett, N. D. R. Bhat, M. Sokolowski, E. Lenc

Abstract: Pulsars have proven instrumental in exploring a wide variety of physics. Pulsars at low radio frequencies is crucial to further our understanding of spectral properties and emission mechanisms.The Murchison Widefield Array Voltage Capture System (MWA-VCS) has been routinely used to study and discover pulsars at low frequencies, offering the unique opportunity of recording complex voltages ,which c… ▽ More Pulsars have proven instrumental in exploring a wide variety of physics. Pulsars at low radio frequencies is crucial to further our understanding of spectral properties and emission mechanisms.The Murchison Widefield Array Voltage Capture System (MWA-VCS) has been routinely used to study and discover pulsars at low frequencies, offering the unique opportunity of recording complex voltages ,which can be off-line beamformed or imaged at millisecond time resolution.Devising imaged-based methods for finding pulsar candidates, which can be verified in beamformed data, can accelerate the complete process and lead to more pulsar detections by reducing the number of tied-array beams required, increasing compute resource efficiency.Despite a factor of ~4 loss in sensitivity, searching for pulsar candidates in images from the MWA-VCS, we can explore a larger parameter space, potentially leading to discoveries of pulsars missed by high-frequency surveys such as pulsars obscured in high-time resolution timeseries data by propagation effects.Image-based searches are also essential to probing parts of parameter space inaccessible to traditional beamformed searches with the MWA.In this paper we describe the innovative approach and capability of dual-processing MWA VCS data, i.e. finding pulsar candidates in these images, and verifying by forming tied-array beam.We developed and tested image-based methods of finding pulsar candidates, based on pulsar properties such as spectral index, polarisation and variability.The efficiency of these methodologies has been verified on known pulsars, and the main limitations explained in terms of sensitivity and low-frequency spectral turnover of some pulsars.No candidates were confirmed to be a new pulsar.This new capability will now be applied to multiple observations to accelerate pulsar discoveries with MWA and speed up future searches with the SKA-Low. △ Less

Submitted 13 December, 2022; originally announced December 2022.

Comments: 12 pages, 9 figures, accepted for publication in Publications of the Astronomical Society of Australia (PASA)

A Cosmological Fireball with Sixteen-Percent Gamma-Ray Radiative Efficiency

Authors: Liang Li, Yu Wang, Felix Ryde, Asaf Pe'er, Bing Zhang, Sylvain Guiriec, Alberto J. Castro-Tirado, D. Alexander Kann, Magnus Axelsson, Kim Page, Peter Veres, P. N. Bhat

Abstract: Gamma-ray bursts (GRBs) are the most powerful explosions in the universe. How efficiently the jet converts its energy to radiation is a long-standing problem and it is poorly constrained. The standard model invokes a relativistic fireball with a bright photosphere emission component. A definitive diagnosis of GRB radiation components and measurement of GRB radiative efficiency require prompt emiss… ▽ More Gamma-ray bursts (GRBs) are the most powerful explosions in the universe. How efficiently the jet converts its energy to radiation is a long-standing problem and it is poorly constrained. The standard model invokes a relativistic fireball with a bright photosphere emission component. A definitive diagnosis of GRB radiation components and measurement of GRB radiative efficiency require prompt emission and afterglow data with high-resolution and wide-band coverage in time and energy. Here we report a comprehensive temporal and spectral analysis of the TeV-emitting bright GRB 190114C. Its fluence is one of the highest of all GRBs detected so far, which allows us to perform a high-resolution study of the prompt emission spectral properties and their temporal evolution down to a timescale of about 0.1 s. We observe that each of the initial pulses has a thermal component contributing $\sim20\%$ of the total energy, the corresponding temperature and the inferred Lorentz factor of the photosphere evolve following broken power-law shapes. From the observation of the non-thermal spectra and the light curve, the onset of afterglow corresponding to the deceleration of the fireball is considered at $\sim 6$~s. By incorporating the thermal and the non-thermal observations, as well as the photosphere and the synchrotron radiative mechanisms, we can directly derive the fireball energy budget with little dependence on hypothetical parameters and to measure a $\sim 16\%$ radiative efficiency for this GRB. With the fireball energy budget derived, the afterglow microphysics parameters can also be constrained directly from the data. △ Less

Submitted 12 February, 2023; v1 submitted 5 December, 2022; originally announced December 2022.

Comments: 27 pages, 8 figures (including 14 panels), 6 tables, accepted for publication in The Astrophysical Journal Letters

Wideband timing of the Parkes Pulsar Timing Array UWL data

Authors: Małgorzata Curyło, Timothy T. Pennucci, Matthew Bailes, N. D. Ramesh Bhat, Andrew D. Cameron, Shi Dai, George Hobbs, Agastya Kapur, Richard N. Manchester, Rami Mandow, Matthew T. Miles, Christopher J. Russell, Daniel J. Reardon, Ryan M. Shannon, Renée Spiewak, Andrew Zic, Xing-Jiang Zhu

Abstract: In 2018 an ultra-wide-bandwidth low-frequency (UWL) receiver was installed on the 64-m Parkes Radio Telescope enabling observations with an instantaneous frequency coverage from 704 to 4032 MHz. Here, we present the analysis of a three-year data set of 35 millisecond pulsars observed with the UWL by the Parkes Pulsar Timing Array (PPTA), using wideband timing methods. The two key differences compa… ▽ More In 2018 an ultra-wide-bandwidth low-frequency (UWL) receiver was installed on the 64-m Parkes Radio Telescope enabling observations with an instantaneous frequency coverage from 704 to 4032 MHz. Here, we present the analysis of a three-year data set of 35 millisecond pulsars observed with the UWL by the Parkes Pulsar Timing Array (PPTA), using wideband timing methods. The two key differences compared to typical narrow-band methods are, firstly, generation of two-dimensional templates accounting for pulse shape evolution with frequency and, secondly, simultaneous measurements of the pulse time-of-arrival (ToA) and dispersion measure (DM). This is the first time that wideband timing has been applied to a uniform data set collected with a single large-fractional bandwidth receiver, for which such techniques were originally developed. As a result of our study, we present a set of profile evolution models and new timing solutions including initial noise analysis. Precision of our ToA and DM measurements is in the range of 0.005 $-$ 2.08 $μ$s and (0.043$-$14.24)$\times10^{-4}$ cm$^{-3}$ pc, respectively, with 94% of the pulsars achieving a median ToA uncertainty of less than 1 $μ$s. △ Less

Submitted 23 November, 2022; originally announced November 2022.

Comments: 29 pages, 12 figures, Accepted for publication in ApJ

A targeted search for repeating fast radio bursts with the MWA

Authors: J. Tian, G. E. Anderson, P. J. Hancock, J. C. A. Miller-Jones, M. Sokolowski, C. W. James, N. D. R. Bhat, N. A. Swainston, D. Ung, B. W. Meyers

Abstract: We present a targeted search for low-frequency (144--215\,MHz) FRB emission from five repeating FRBs using 23.3\,hr of archival data taken with the Murchison Widefield Array (MWA) Voltage Capture System (VCS) between 2014 September and 2020 May. This is the first time that the MWA VCS has been used to search for FRB signals from known repeaters, which enables much more sensitive FRB searches than… ▽ More We present a targeted search for low-frequency (144--215\,MHz) FRB emission from five repeating FRBs using 23.3\,hr of archival data taken with the Murchison Widefield Array (MWA) Voltage Capture System (VCS) between 2014 September and 2020 May. This is the first time that the MWA VCS has been used to search for FRB signals from known repeaters, which enables much more sensitive FRB searches than previously performed with the standard MWA correlator mode. We performed a standard single pulse search with a temporal and spectral resolution of $400\,μ$s and 10\,kHz, respectively, over a $100\,\text{pc}\,\text{cm}^{-3}$ dispersion measure (DM) range centred at the known DM of each studied repeating FRB. No FRBs exceeding a $6σ$ threshold were detected. The fluence upper limits in the range of 32--1175\,Jy\,ms and 36--488\,Jy\,ms derived from 10 observations of FRB 20190711A and four observations of FRB 20201124A respectively, allow us to constrain the spectral indices of their bursts to $\gtrsim-1$ if these two repeaters were active during the MWA observations. If free-free absorption is responsible for our non-detection, we can constrain the size of the absorbing medium in terms of the electron temperature $T$ to $<1.00\times(T/10^4\text{K})^{-1.35}\,\text{pc}$, $<0.92\times(T/10^4\text{K})^{-1.35}\,\text{pc}$ and $<[0.22\text{--}2.50]\times(T/10^4\text{K})^{-1.35}\,\text{pc}$ for FRB 20190117A, 20190711A, and 20201124A, respectively. However, given that the activities of these repeaters are not well characterised, our non-detections could also suggest they were inactive during the MWA observations. △ Less

Submitted 21 November, 2022; originally announced November 2022.

Comments: Accepted for publication in MNRAS

Constraining ultralight vector dark matter with the Parkes Pulsar Timing Array second data release

Authors: Yu-Mei Wu, Zu-Cheng Chen, Qing-Guo Huang, Xingjiang Zhu, N. D. Ramesh Bhat, Yi Feng, George Hobbs, Richard N. Manchester, Christopher J. Russell, R. M. Shannon

Abstract: Composed of ultralight bosons, fuzzy dark matter provides an intriguing solution to challenges that the standard cold dark matter model encounters on sub-galactic scales. The ultralight dark matter with mass $m\sim10^{-23} \rm{eV}$ will induce a periodic oscillation in gravitational potentials with a frequency in the nanohertz band, leading to observable effects in the arrival times of radio pulse… ▽ More Composed of ultralight bosons, fuzzy dark matter provides an intriguing solution to challenges that the standard cold dark matter model encounters on sub-galactic scales. The ultralight dark matter with mass $m\sim10^{-23} \rm{eV}$ will induce a periodic oscillation in gravitational potentials with a frequency in the nanohertz band, leading to observable effects in the arrival times of radio pulses from pulsars. Unlike scalar dark matter, pulsar timing signals induced by the vector dark matter are dependent on the oscillation direction of the vector fields. In this work, we search for ultralight vector dark matter in the mass range of $[2\times 10^{-24}, 2\times 10^{-22}]{\rm{eV}}$ through its gravitational effect in the Parkes Pulsar Timing Array (PPTA) second data release. Since no statistically significant detection is made, we place $95\%$ upper limits on the local dark matter density as $ρ_{\rm{\tiny{VF}}} \lesssim 5{\rm{GeV/cm^{3}}}$ for $m\lesssim 10^{-23}{\rm{eV}}$. As no preferred direction is found for the vector dark matter, these constraints are comparable to those given by the scalar dark matter search with an earlier 12-year data set of PPTA. △ Less

Submitted 7 October, 2022; originally announced October 2022.

Comments: 9 pages, 3 figures, 2 tables; accepted for publication as a Letter in Phys. Rev. D

pulsar_spectra: A pulsar flux density catalogue and spectrum fitting repository

Authors: N. A. Swainston, C. P. Lee, S. J. McSweeney, N. D. R. Bhat

Abstract: We present the pulsar_spectra software repository, an open-source pulsar flux density catalogue and automated spectral fitting software that finds the best spectral model and produces publication-quality plots. The Python-based software includes features that enable users in the astronomical community to add newly published spectral measurements to the catalogue as they become available. The spect… ▽ More We present the pulsar_spectra software repository, an open-source pulsar flux density catalogue and automated spectral fitting software that finds the best spectral model and produces publication-quality plots. The Python-based software includes features that enable users in the astronomical community to add newly published spectral measurements to the catalogue as they become available. The spectral fitting software is an implementation of the method described in Jankowski et al. (2018) which uses robust statistical methods to decide on the best-fitting model for individual pulsar spectra. pulsar_spectra is motivated by the need for a centralised repository for pulsar flux density measurements to make published measurements more accessible to the astronomical community and provide a suite of tools for measuring spectra. △ Less

Submitted 27 September, 2022; originally announced September 2022.

Spectral analysis of 22 radio pulsars using SKA-Low precursor stations

Authors: C. P. Lee, N. D. R. Bhat, M. Sokolowski, N. A. Swainston, D. Ung, A. Magro, R. Chiello

Abstract: We present the first observational study of pulsars performed with the second-generation precursor stations to the low-frequency component of the Square Kilometre Array (SKA-Low): the Aperture Array Verification System 2 (AAVS2) and the Engineering Development Array 2 (EDA2). Using the SKA-Low stations, we have observed 100 southern-sky pulsars between 70-350 MHz, including follow-up observations… ▽ More We present the first observational study of pulsars performed with the second-generation precursor stations to the low-frequency component of the Square Kilometre Array (SKA-Low): the Aperture Array Verification System 2 (AAVS2) and the Engineering Development Array 2 (EDA2). Using the SKA-Low stations, we have observed 100 southern-sky pulsars between 70-350 MHz, including follow-up observations at multiple frequencies for a selected sample of bright pulsars. These observations have yielded detections of 22 pulsars, including the lowest-frequency detections ever published for 6 pulsars, despite the modest sensitivity of initial system where the recording bandwidth is limited to ~1 MHz. By comparing simultaneous flux density measurements obtained with the SKA-Low stations and performing rigorous electromagnetic simulations, we verify the accuracy of the SKA-Low sensitivity simulation code presented in Sokolowski et al. (2022). Furthermore, we perform model fits to the radio spectra of the detected pulsars using the method developed by Jankowski et al. (2018), including 9 pulsars which were not fitted in the original work. We robustly classify the spectra into 5 morphological classes and find that all but one pulsar exhibit deviations from simple power-law behaviour. These findings suggest that pulsars with well-determined spectra are more likely to show spectral flattening or turn-over than average. Our work demonstrates how SKA-Low stations can be meaningfully used for scientifically useful measurements and analysis of pulsar radio spectra, which are important inputs for informing pulsar surveys and science planned with the SKA-Low. △ Less

Submitted 15 August, 2022; originally announced August 2022.

Comments: 16 pages, 6 figures, 4 tables, accepted for publication in PASA

Evaluating the prevalence of spurious correlations in pulsar timing array datasets

Authors: Andrew Zic, George Hobbs, R. M. Shannon, Daniel Reardon, Boris Goncharov, N. D. Ramesh Bhat, Andrew Cameron, Shi Dai, J. R. Dawson, Matthew Kerr, R. N. Manchester, Rami Mandow, Tommy Marshman, Christopher J. Russell, Nithyanandan Thyagarajan, X. -J. Zhu

Abstract: Pulsar timing array collaborations have recently reported evidence for a noise process with a common spectrum among the millisecond pulsars in the arrays. The spectral properties of this common-noise process are consistent with expectations for an isotropic gravitational-wave background (GWB) from inspiralling supermassive black-hole binaries. However, recent simulation analyses based on Parkes Pu… ▽ More Pulsar timing array collaborations have recently reported evidence for a noise process with a common spectrum among the millisecond pulsars in the arrays. The spectral properties of this common-noise process are consistent with expectations for an isotropic gravitational-wave background (GWB) from inspiralling supermassive black-hole binaries. However, recent simulation analyses based on Parkes Pulsar Timing Array data indicate that such a detection may arise spuriously. In this paper, we use simulated pulsar timing array datasets to further test the robustness of the inference methods for spectral and spatial correlations from a GWB. Expanding on our previous results, we find strong support (Bayes factors exceeding $10^5$) for the presence of a common-spectrum noise process in datasets where no common process is present, under a wide range of timing noise prescriptions per pulsar. We show that these results are highly sensitive to the choice of Bayesian priors on timing noise parameters, with priors that more closely match the injected distributions of timing noise parameters resulting in diminished support for a common-spectrum noise process. These results emphasize shortcomings in current methods for inferring the presence of a common-spectrum process, and imply that the detection of a common process is not a reliable precursor to detection of the GWB. Future searches for the nanohertz GWB should remain focussed on detecting spatial correlations, and make use of more tailored specifications for a common-spectrum noise process. △ Less

Submitted 25 July, 2022; originally announced July 2022.

Comments: 12 pages, 12 figures. Accepted for publication in MNRAS

Consistency of the Parkes Pulsar Timing Array Signal with a Nanohertz Gravitational-wave Background

Authors: Boris Goncharov, Eric Thrane, Ryan M. Shannon, Jan Harms, N. D. Ramesh Bhat, George Hobbs, Matthew Kerr, Richard N. Manchester, Daniel J. Reardon, Christopher J. Russell, Xing-Jiang Zhu, Andrew Zic

Abstract: Pulsar timing array experiments have recently reported strong evidence for a common-spectrum stochastic process with a strain spectral index consistent with that expected of a nanohertz-frequency gravitational-wave background, but with negligible yet non-zero evidence for spatial correlations required for a definitive detection. However, it was pointed out by the Parkes Pulsar Timing Array (PPTA)… ▽ More Pulsar timing array experiments have recently reported strong evidence for a common-spectrum stochastic process with a strain spectral index consistent with that expected of a nanohertz-frequency gravitational-wave background, but with negligible yet non-zero evidence for spatial correlations required for a definitive detection. However, it was pointed out by the Parkes Pulsar Timing Array (PPTA) collaboration that the same models used in recent analyses resulted in strong evidence for a common-spectrum process in simulations where none is present. In this work, we introduce a methodology to distinguish pulsar power spectra with the same amplitude from noise power spectra of similar but distinct amplitudes. The former is the signature of a spatially uncorrelated pulsar term of a nanohertz gravitational-wave background, whereas the latter could represent ensemble pulsar noise properties. We test the methodology on simulated data sets. We find that the reported common process in PPTA pulsars is indeed consistent with the spectral feature of a pulsar term. We recommend this methodology as one of the validity tests that the real astrophysical and cosmological backgrounds should pass, as well as for inferences about the spatially uncorrelated component of the background. △ Less

Submitted 21 June, 2022; v1 submitted 8 June, 2022; originally announced June 2022.

Comments: 9 pages, 2 figures

Journal ref: The Astrophysical Journal Letters, Volume 932, Number 2 (2022)

Independent discovery of a nulling pulsar with unusual sub-pulse drifting properties with the Murchison Widefield Array

Authors: Samuel J. McSweeney, N. D. Ramesh Bhat, Nicholas A. Swainston, Keegan R. Smith, Sanjay Kudale, Paul Hancock, Willem van Straten, Shi Dai, Ryan M. Shannon, Steven J. Tingay, Melanie Johnston-Hollitt, David L. Kaplan, Mia Walker

Abstract: We report the independent discovery of PSR J0027-1956 with the Murchison Widefield Array (MWA) in the ongoing Southern-sky MWA Rapid Two-meter (SMART) pulsar survey. J0027-1956 has a period of ~1.306 s, a dispersion measure (DM) of ~20.869 pc cm^-3 , and a nulling fraction of ~77%. This pulsar highlights the advantages of the survey's long dwell times (~80 min), which, when fully searched, will be… ▽ More We report the independent discovery of PSR J0027-1956 with the Murchison Widefield Array (MWA) in the ongoing Southern-sky MWA Rapid Two-meter (SMART) pulsar survey. J0027-1956 has a period of ~1.306 s, a dispersion measure (DM) of ~20.869 pc cm^-3 , and a nulling fraction of ~77%. This pulsar highlights the advantages of the survey's long dwell times (~80 min), which, when fully searched, will be sensitive to the expected population of similarly bright, intermittent pulsars with long nulls. A single-pulse analysis in the MWA's 140-170 MHz band also reveals a complex sub-pulse drifting behavior, including both rapid changes of the drift rate characteristic of mode switching pulsars, as well as a slow, consistent evolution of the drift rate within modes. In some longer drift sequences, interruptions in the otherwise smooth drift rate evolution occur preferentially at a particular phase, typically lasting a few pulses. These properties make this pulsar an ideal test bed for prevailing models of drifting behavior such as the carousel model. △ Less

Submitted 1 June, 2022; originally announced June 2022.

Comments: Accepted for publication in ApJ

High time resolution search for prompt radio emission from the long GRB 210419A with the Murchison Widefield Array

Authors: J. Tian, G. E. Anderson, P. J. Hancock, J. C. A. Miller-Jones, M. Sokolowski, N. A. Swainston, A. Rowlinson, A. Williams, D. L. Kaplan, N. Hurley-Walker, J. Morgan, N. D. R. Bhat, D. Ung, S. Tingay, K. W. Bannister, M. E. Bell, B. W. Meyers, M. Walker

Abstract: We present a low-frequency (170\textendash200\,MHz) search for prompt radio emission associated with the long GRB 210419A using the rapid-response mode of the Murchison Widefield Array (MWA), triggering observations with the Voltage Capture System (VCS) for the first time. The MWA began observing GRB 210419A within 89\,s of its detection by \textit{Swift}, enabling us to capture any dispersion del… ▽ More We present a low-frequency (170\textendash200\,MHz) search for prompt radio emission associated with the long GRB 210419A using the rapid-response mode of the Murchison Widefield Array (MWA), triggering observations with the Voltage Capture System (VCS) for the first time. The MWA began observing GRB 210419A within 89\,s of its detection by \textit{Swift}, enabling us to capture any dispersion delayed signal emitted by this GRB for a typical range of redshifts. We conducted a standard single pulse search with a temporal and spectral resolution of $100\,\upmu$s and 10\,kHz over a broad range of dispersion measures from 1 to $5000\,\text{pc}\,\text{cm}^{-3}$, but none were detected. However, fluence upper limits of $77\text{--}224$\,Jy\,ms derived over a pulse width of $0.5\text{--}10$\,ms and a redshift of $0.6<z<4$ are some of the most stringent at low radio frequencies. We compared these fluence limits to the GRB jet-interstellar medium (ISM) interaction model, placing constraints on the fraction of magnetic energy ($ε_{\text{B}}\lesssim[0.05 \text{--} 0.1]$). We also searched for signals during the X-ray flaring activity of GRB 210419A on minute timescales in the image domain and found no emission, resulting in an intensity upper limit of $0.57\,\text{Jy}\,\text{beam}^{-1}$, corresponding to a constraint of $ε_{\text{B}}\lesssim10^{-3}$. Our non-detection could imply that GRB 210419A was at a high redshift, there was not enough magnetic energy for low-frequency emission, or that the radio waves did not escape from the GRB environment. △ Less

Submitted 26 May, 2022; originally announced May 2022.

Comments: Accepted for publication in MNRAS

Detection of frequency-dependent dispersion measure toward the millisecond pulsar PSR J2241-5236 from contemporaneous wide-band observations

Authors: Dilpreet Kaur, N. D. Ramesh Bhat, Shi Dai, Samuel J. McSweeney, Ryan M. Shannon, Sanjay Kudale, Willem van Straten

Abstract: Making precise measurements of pulsar dispersion measures (DMs) and applying suitable corrections for them is amongst the major challenges in high-precision timing programmes such as pulsar timing arrays (PTAs). While the advent of wide-band pulsar instrumentation can enable more precise DM measurements and thence improved timing precision, it also necessitates doing careful assessments of frequen… ▽ More Making precise measurements of pulsar dispersion measures (DMs) and applying suitable corrections for them is amongst the major challenges in high-precision timing programmes such as pulsar timing arrays (PTAs). While the advent of wide-band pulsar instrumentation can enable more precise DM measurements and thence improved timing precision, it also necessitates doing careful assessments of frequency-dependent (chromatic) DMs that was theorised by Cordes et al. (2016). Here we report the detection of such an effect in broadband observations of the millisecond pulsar PSR J2241-5236, a high-priority target for current and future PTAs. The observations were made contemporaneously using the wide-band receivers and capabilities now available at the Murchison Widefield Array (MWA), the upgraded Giant Metrewave Radio Telescope (uGMRT) and the Parkes telescopes, thus providing an unprecedentedly large frequency coverage from 80 MHz to 4 GHz. Our analysis shows the measurable changes in DM that scale with the observing frequency ($ν$) as $\rm δDM \sim ν^{2.5 \pm 0.1}$. We discuss the potential implications of such a frequency dependence in the measured DMs, and the likely impact on the timing noise budget, and comment on the usefulness of low-frequency observations in advancing PTA efforts. △ Less

Submitted 17 April, 2022; originally announced April 2022.

Comments: 14 pages, 5 figures, 3 tables, Accepted for publication in The ApJ Letters

The High Time Resolution Universe Pulsar Survey -- XVII. PSR J1325-6253, a low eccentricity double neutron star system from an ultra-stripped supernova

Authors: R. Sengar, V. Balakrishnan, S. Stevenson, M. Bailes, E. D. Barr, N. D. R. Bhat, M. Burgay, M. C. i Bernadich, A. D. Cameron, D. J. Champion, W. Chen, C. M. L. Flynn, A. Jameson, S. Johnston, M. J. Keith, M. Kramer, V. Morello, C. Ng, A. Possenti, B. Stappers, R. M. Shannon, W. van Straten, J. Wongphechauxsorn

Abstract: The observable population of double neutron star (DNS) systems in the Milky Way allow us to understand the nature of supernovae and binary stellar evolution. Until now, all DNS systems in wide orbits ($ P_{\textrm{orb}}>$ 1~day) have been found to have orbital eccentricities, $e > 0.1$. In this paper, we report the discovery of pulsar PSR J1325$-$6253: a DNS system in a 1.81 day orbit with a surpr… ▽ More The observable population of double neutron star (DNS) systems in the Milky Way allow us to understand the nature of supernovae and binary stellar evolution. Until now, all DNS systems in wide orbits ($ P_{\textrm{orb}}>$ 1~day) have been found to have orbital eccentricities, $e > 0.1$. In this paper, we report the discovery of pulsar PSR J1325$-$6253: a DNS system in a 1.81 day orbit with a surprisingly low eccentricity of just $e = 0.064$. Through 1.4 yr of dedicated timing with the Parkes radio telescope we have been able to measure its rate of advance of periastron, $\dotω=0.138 \pm 0.002$ $\rm deg$ $\rm yr^{-1}$. If this induced $\dotω$ is solely due to general relativity then the total mass of the system is, $M_{\rm sys} = 2.57 \pm 0.06$ M$_{\odot}$. Assuming an edge-on orbit the minimum companion mass is constrained to be $M_\mathrm{c,min}>0.98$ M$_{\odot}$ which implies the pulsar mass is $M_\mathrm{p,max}<1.59 $ M$_{\odot}$. Its location in the $P$-$\dot{P}$ diagram suggests that, like other DNS systems, PSR J1325$-$6253 is a recycled pulsar and if its mass is similar to the known examples ($>1.3$ M$_\odot$), then the companion neutron star is probably less than $\sim1.25$ M$_\odot$ and the system is inclined at about $50^{\circ}$-$60^{\circ}$. The low eccentricity along with the wide orbit of the system strongly favours a formation scenario involving an ultra-stripped supernova explosion. △ Less

Submitted 26 May, 2022; v1 submitted 14 April, 2022; originally announced April 2022.

Comments: 12 pages, 7 figures, 1 table, accepted for publication in MNRAS

The MeerTime Pulsar Timing Array -- A Census of Emission Properties and Timing Potential

Authors: R. Spiewak, M. Bailes, M. T. Miles, A. Parthasarathy, D. J. Reardon, M. Shamohammadi, R. M. Shannon, N. D. R. Bhat, S. Buchner, A. D. Cameron, F. Camilo, M. Geyer, S. Johnston, A. Karastergiou, M. Keith, M. Kramer, M. Serylak, W. van Straten, G. Theureau, V. Venkatraman Krishnan

Abstract: MeerTime is a five-year Large Survey Project to time pulsars with MeerKAT, the 64-dish South African precursor to the Square Kilometre Array. The science goals for the programme include timing millisecond pulsars (MSPs) to high precision (< 1 $μ$s) to study the Galactic MSP population and to contribute to global efforts to detect nanohertz gravitational waves with the International Pulsar Timing A… ▽ More MeerTime is a five-year Large Survey Project to time pulsars with MeerKAT, the 64-dish South African precursor to the Square Kilometre Array. The science goals for the programme include timing millisecond pulsars (MSPs) to high precision (< 1 $μ$s) to study the Galactic MSP population and to contribute to global efforts to detect nanohertz gravitational waves with the International Pulsar Timing Array (IPTA). In order to plan for the remainder of the programme and to use the allocated time most efficiently, we have conducted an initial census with the MeerKAT "L-band" receiver of 189 MSPs visible to MeerKAT and here present their dispersion measures, polarization profiles, polarization fractions, rotation measures, flux density measurements, spectral indices, and timing potential. As all of these observations are taken with the same instrument (which uses coherent dedispersion, interferometric polarization calibration techniques, and a uniform flux scale), they present an excellent resource for population studies. We used wideband pulse portraits as timing standards for each MSP and demonstrated that the MeerTime Pulsar Timing Array (MPTA) can already contribute significantly to the IPTA as it currently achieves better than 1 $μ$s timing accuracy on 89 MSPs (observed with fortnightly cadence). By the conclusion of the initial five-year MeerTime programme in July 2024, the MPTA will be extremely significant in global efforts to detect the gravitational wave background with a contribution to the detection statistic comparable to other long-standing timing programmes. △ Less

Submitted 8 April, 2022; originally announced April 2022.

Comments: Accepted to PASA. 27 figures. Data to be made available under the DOI 10.5281/zenodo.5347875 at the time of publication

MWA Tied-Array Processing IV: A Multi-Pixel Beamformer for Pulsar Surveys and Ionospheric Corrected Localisation

Authors: N. A. Swainston, N. D. R. Bhat, I. S. Morrison, S. J. McSweeney, S. M. Ord, S. E. Tremblay, M. Sokolowski

Abstract: The Murchison Widefield Array (MWA) is a low-frequency aperture array capable of high-time and frequency resolution astronomy applications such as pulsar studies. The large field-of-view of the MWA (hundreds of square degrees) can also be exploited to attain fast survey speeds for all-sky pulsar search applications, but to maximise sensitivity requires forming thousands of tied-array beams from ea… ▽ More The Murchison Widefield Array (MWA) is a low-frequency aperture array capable of high-time and frequency resolution astronomy applications such as pulsar studies. The large field-of-view of the MWA (hundreds of square degrees) can also be exploited to attain fast survey speeds for all-sky pulsar search applications, but to maximise sensitivity requires forming thousands of tied-array beams from each voltage-capture observation. The necessity of using calibration solutions that are separated from the target observation both temporally and spatially makes pulsar observations vulnerable to uncorrected, frequency-dependent positional offsets due to the ionosphere. These offsets may be large enough to move the source away from the centre of the tied-array beam, incurring sensitivity drops of $\sim$30-50\% in Phase II extended array configuration. We analyse these offsets in pulsar observations and develop a method for mitigating them, improving both the source position accuracy and the sensitivity. This analysis prompted the development of a multi-pixel beamforming functionality that can generate dozens of tied-array beams simultaneously, which runs a factor of ten times faster compared to the original single-pixel version. This enhancement makes it feasible to observe multiple pulsars within the vast field of view of the MWA and supports the ongoing large-scale pulsar survey efforts with the MWA. We explore the extent to which ionospheric offset correction will be necessary for the MWA Phase III and the low-frequency Square Kilometre Array (SKA-Low). △ Less

Submitted 29 March, 2022; v1 submitted 26 March, 2022; originally announced March 2022.

Comments: 10 pages, 5 figures

High-precision search for dark photon dark matter with the Parkes Pulsar Timing Array

Authors: Xiao Xue, Zi-Qing Xia, Xingjiang Zhu, Yue Zhao, Jing Shu, Qiang Yuan, N. D. Ramesh Bhat, Andrew D. Cameron, Shi Dai, Yi Feng, Boris Goncharov, George Hobbs, Eric Howard, Richard N. Manchester, Aditya Parthasarathy, Daniel J. Reardon, Christopher J. Russell, Ryan M. Shannon, Renée Spiewak, Nithyanandan Thyagarajan, Jingbo Wang, Lei Zhang, Songbo Zhang

Abstract: The nature of dark matter remains obscure in spite of decades of experimental efforts. The mass of dark matter candidates can span a wide range, and its coupling with the Standard Model sector remains uncertain. All these unknowns make the etection of dark matter extremely challenging. Ultralight dark matter, with $m \sim10^{-22}$ eV, is proposed to reconcile the disagreements between observations… ▽ More The nature of dark matter remains obscure in spite of decades of experimental efforts. The mass of dark matter candidates can span a wide range, and its coupling with the Standard Model sector remains uncertain. All these unknowns make the etection of dark matter extremely challenging. Ultralight dark matter, with $m \sim10^{-22}$ eV, is proposed to reconcile the disagreements between observations and predictions from simulations of small-scale structures in the cold dark matter paradigm, while remaining consistent with other observations. Because of its large de Broglie wavelength and large local occupation number within galaxies, ultralight dark matter behaves like a coherently oscillating background field with an oscillating frequency dependent on its mass. If the dark matter particle is a spin-1 dark photon, such as the $U(1)_B$ or $U(1)_{B-L}$ gauge boson, it can induce an external oscillating force and lead to displacements of test masses. Such an effect would be observable in the form of periodic variations in the arrival times of radio pulses from highly stable millisecond pulsars. In this study, we search for evidence of ultralight dark photon dark matter (DPDM) using 14-year high-precision observations of 26 pulsars collected with the Parkes Pulsar Timing Array. While no statistically significant signal is found, we place constraints on coupling constants for the $U(1)_B$ and $U(1)_{B-L}$ DPDM. Compared with other experiments, the limits on the dimensionless coupling constant $ε$ achieved in our study are improved by up to two orders of magnitude when the dark photon mass is smaller than $3\times10^{-22}$~eV ($10^{-22}$~eV) for the $U(1)_{B}$ ($U(1)_{B-L}$) scenario. △ Less

Submitted 26 December, 2021; v1 submitted 14 December, 2021; originally announced December 2021.

Comments: 12 pages, 4 figures

Journal ref: Phys. Rev. Research 4, L012022 (2022)

Revisiting the sub-pulse drifting phenomenon in PSR J1822-2256: Drift Modes, Sparks, and Emission Heights

Authors: Parul Janagal, Manoneeta Chakraborty, N. D. Ramesh Bhat, Bhaswati Bhattacharya, Samuel J. McSweeney

Abstract: Sub-pulse drifting in pulsar radio emission is considered to be one of the most promising phenomenon for uncovering the underlying physical processes. Here we present a detailed study of such a phenomenon in observations of PSR J1822$-$2256, made using the upgraded Giant Meterwave Radio Telescope (uGMRT). Observations were made simultaneously using the Band 3 (300-500 MHz) and Band 4 (550-750 MHz)… ▽ More Sub-pulse drifting in pulsar radio emission is considered to be one of the most promising phenomenon for uncovering the underlying physical processes. Here we present a detailed study of such a phenomenon in observations of PSR J1822$-$2256, made using the upgraded Giant Meterwave Radio Telescope (uGMRT). Observations were made simultaneously using the Band 3 (300-500 MHz) and Band 4 (550-750 MHz) receivers of the uGMRT. The pulsar is known to exhibit subpulse drifting, mode changing, and nulling. Our observations reveal four distinct sub-pulse drifting modes of emission (A, B, C, and D) for this pulsar, with the drift periodicities of 17.9 $P_1$, 5.8 $P_1$, 8 $P_1$, 14.1 $P_1$, respectively (where $P_1$ is the pulsar rotation period), two of which exhibit some new features that were not reported in the previous studies. We also investigate the possible spark configuration, characterised by the number of sparks ($n$) in the carousel patterns of these four drift modes, and our analysis suggests two representative solutions for the number of sparks for a carousel rotation period, $P_4$, which lies in the range of $13$ to $16$. The large frequency coverage of our data (300-750 MHz) is also leveraged to explore the frequency dependence of single-pulse characteristics of the pulsar emission, particularly the frequency-dependent subpulse behaviour and the emission heights for the observed drift modes. Our analysis suggests a clear modal dependence of inferred emission heights. We discuss the implications for the pulsar emission mechanism and its relation to the proposed spark configuration. △ Less

Submitted 13 November, 2021; originally announced November 2021.

Comments: 12 pages, 9 figures

Fermi-GBM Observations of GRB 210812A: Signatures of a Million Solar Mass Gravitational Lens

Authors: P. Veres, N. Bhat, N. Fraija, S. Lesage

Abstract: Observing gravitationally lensed objects in the time domain is difficult, and well-observed time-varying sources are rare. Lensed gamma-ray bursts (GRBs) offer improved timing precision to this class of objects complementing observations of quasars and supernovae. The rate of lensed GRBs is highly uncertain, approximately 1 in 1000. The Gamma-ray Burst Monitor (GBM) onboard the Fermi Gamma-ray Spa… ▽ More Observing gravitationally lensed objects in the time domain is difficult, and well-observed time-varying sources are rare. Lensed gamma-ray bursts (GRBs) offer improved timing precision to this class of objects complementing observations of quasars and supernovae. The rate of lensed GRBs is highly uncertain, approximately 1 in 1000. The Gamma-ray Burst Monitor (GBM) onboard the Fermi Gamma-ray Space Telescope has observed more than 3000 GRBs making it an ideal instrument to uncover lensed bursts. Here we present observations of GRB 210812A showing two emission episodes, separated by 33.3 s, and with flux ratio of about 4.5. An exhaustive temporal and spectral analysis shows that the two emission episodes have the same pulse and spectral shape, which poses challenges to GRB models. We report multiple lines of evidence for a gravitational lens origin. In particular, modeling the lightcurve using nested sampling we uncover strong evidence in favor of the lensing scenario. Assuming a point mass lens, the mass of the lensing object is about 1 million solar masses. High-resolution radio imaging is needed for future lens candidates to derive tighter constraints. △ Less

Submitted 12 October, 2021; originally announced October 2021.

Comments: 13 pages, 8 figures, AAS journals accepted (ApJL)

Constraining cosmological phase transitions with the Parkes Pulsar Timing Array

Authors: Xiao Xue, Ligong Bian, Jing Shu, Qiang Yuan, Xingjiang Zhu, N. D. Ramesh Bhat, Shi Dai, Yi Feng, Boris Goncharov, George Hobbs, Eric Howard, Richard N. Manchester, Christopher J. Russell, Daniel J. Reardon, Ryan M. Shannon, Renée Spiewak, Nithyanandan Thyagarajan, Jingbo Wang

Abstract: A cosmological first-order phase transition is expected to produce a stochastic gravitational wave background. If the phase transition temperature is on the MeV scale, the power spectrum of the induced stochastic gravitational waves peaks around nanohertz frequencies, and can thus be probed with high-precision pulsar timing observations. We search for such a stochastic gravitational wave backgroun… ▽ More A cosmological first-order phase transition is expected to produce a stochastic gravitational wave background. If the phase transition temperature is on the MeV scale, the power spectrum of the induced stochastic gravitational waves peaks around nanohertz frequencies, and can thus be probed with high-precision pulsar timing observations. We search for such a stochastic gravitational wave background with the latest data set of the Parkes Pulsar Timing Array. We find no evidence for a Hellings-Downs spatial correlation as expected for a stochastic gravitational wave background. Therefore, we present constraints on first-order phase transition model parameters. Our analysis shows that pulsar timing is particularly sensitive to the low-temperature ($T \sim 1 - 100$ MeV) phase transition with a duration $(β/H_*)^{-1}\sim 10^{-2}-10^{-1}$ and therefore can be used to constrain the dark and QCD phase transitions. △ Less

Submitted 25 November, 2022; v1 submitted 6 October, 2021; originally announced October 2021.

Comments: 9 pages, 4 figures, 2 tables

Journal ref: Phys. Rev. Lett. 127, 251303, Published 15 December 2021

The Parkes pulsar timing array second data release: Timing analysis

Authors: D. J. Reardon, R. M. Shannon, A. D. Cameron, B. Goncharov, G. B. Hobbs, H. Middleton, M. Shamohammadi, N. Thyagarajan, M. Bailes, N. D. R. Bhat, S. Dai, M. Kerr, R. N. Manchester, C. J. Russell, R. Spiewak, J. B. Wang, X. J. Zhu

Abstract: The main goal of pulsar timing array experiments is to detect correlated signals such as nanohertz-frequency gravitational waves. Pulsar timing data collected in dense monitoring campaigns can also be used to study the stars themselves, their binary companions, and the intervening ionised interstellar medium. Timing observations are extraordinarily sensitive to changes in path length between the p… ▽ More The main goal of pulsar timing array experiments is to detect correlated signals such as nanohertz-frequency gravitational waves. Pulsar timing data collected in dense monitoring campaigns can also be used to study the stars themselves, their binary companions, and the intervening ionised interstellar medium. Timing observations are extraordinarily sensitive to changes in path length between the pulsar and the Earth, enabling precise measurements of the pulsar positions, distances and velocities, and the shapes of their orbits. Here we present a timing analysis of 25 pulsars observed as part of the Parkes Pulsar Timing Array (PPTA) project over time spans of up to 24 yr. The data are from the second data release of the PPTA, which we have extended by including legacy data. We make the first detection of Shapiro delay in four Southern pulsars (PSRs J1017$-$7156, J1125$-$6014, J1545$-$4550, and J1732$-$5049), and of parallax in six pulsars. The prominent Shapiro delay of PSR J1125$-$6014 implies a neutron star mass of $M_p = 1.5 \pm 0.2 M_\odot$ (68% credibility interval). Measurements of both Shapiro delay and relativistic periastron advance in PSR J1600$-$3053 yield a large but uncertain pulsar mass of $M_p = 2.06^{+0.44}_{-0.41}$ M$_\odot$ (68% credibility interval). We measure the distance to PSR J1909$-$3744 to a precision of 10 lyr, indicating that for gravitational wave periods over a decade, the pulsar provides a coherent baseline for pulsar timing array experiments. △ Less

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

Comments: 44 pages, 32 figures, 7 tables. Accepted for publication in MNRAS; v2 includes corrected typos in parameter tables and reference added

The Location of Young Pulsar PSR J0837$-$2454: Galactic Halo or Local Supernova Remnant?

Authors: Nihan Pol, Sarah Burke-Spolaor, Natasha Hurley-Walker, Harsha Blumer, Simon Johnston, Michael Keith, Evan F. Keane, Marta Burgay, Andrea Possenti, Emily Petroff, N. D. Ramesh Bhat

Abstract: We present the discovery and timing of the young (age $\sim 28.6$ kyr) pulsar PSR J0837$-$2454. Based on its high latitude ($b = 9.8^{\circ}$) and dispersion measure (DM $ = 143$~pc~cm$^{-3}$), the pulsar appears to be at a $z$-height of $>$1 kpc above the Galactic plane, but near the edge of our Galaxy. This is many times the observed scale height of the canonical pulsar population, which suggest… ▽ More We present the discovery and timing of the young (age $\sim 28.6$ kyr) pulsar PSR J0837$-$2454. Based on its high latitude ($b = 9.8^{\circ}$) and dispersion measure (DM $ = 143$~pc~cm$^{-3}$), the pulsar appears to be at a $z$-height of $>$1 kpc above the Galactic plane, but near the edge of our Galaxy. This is many times the observed scale height of the canonical pulsar population, which suggests this pulsar may have been born far out of the plane. If accurate, the young age and high $z$-height imply that this is the first pulsar known to be born from a runaway O/B star. In follow-up imaging with the Australia Telescope Compact Array (ATCA), we detect the pulsar with a flux density $S_{1400} = 0.18 \pm 0.05$ mJy. We do not detect an obvious supernova remnant around the pulsar in our ATCA data, but we detect a co-located, low-surface-brightness region of $\sim$1.5$^\circ$ extent in archival Galactic and Extragalactic All-sky MWA Survey data. We also detect co-located H$α$ emission from the Southern H$α$ Sky Survey Atlas. Distance estimates based on these two detections come out to $\sim$0.9 kpc and $\sim$0.2 kpc respectively, both of which are much smaller than the distance predicted by the NE2001 model ($6.3$ kpc) and YMW model ($>25$ kpc) and place the pulsar much closer to the plane of the Galaxy. If the pulsar/remnant association holds, this result also highlights the inherent difficulty in the classification of transients as "Galactic" (pulsar) or "extragalactic" (fast radio burst) toward the Galactic anti-center based solely on the modeled Galactic electron contribution to a detection. △ Less

Submitted 23 April, 2021; originally announced April 2021.

Comments: Published in ApJ. 12 pages, 9 figures, 2 tables

Discovery of a steep-spectrum low-luminosity pulsar with the Murchison Widefield Array

Authors: N. A. Swainston, N. D. R. Bhat, M. Sokolowski, S. J. McSweeney, S. Kudale, S. Dai, K. R. Smith, I. S. Morrison, R. M. Shannon, W. van Straten, M. Xue, S. M. Ord, S. E. Tremblay, B. W. Meyers, A. Williams, G. Sleap, M. Johnston-Hollitt, D. L. Kaplan, S. J. Tingay, R. B. Wayth

Abstract: We report the discovery of the first new pulsar with the Murchison Widefield Array (MWA), PSR J0036$-$1033, a long-period (0.9 s) nonrecycled pulsar with a dispersion measure (DM) of 23.1 ${\rm pc\,cm^{-3}}$. It was found after processing only a small fraction ($\sim$1%) of data from an ongoing all-sky pulsar survey. Follow-up observations have been made with the MWA, the upgraded Giant Metrewave… ▽ More We report the discovery of the first new pulsar with the Murchison Widefield Array (MWA), PSR J0036$-$1033, a long-period (0.9 s) nonrecycled pulsar with a dispersion measure (DM) of 23.1 ${\rm pc\,cm^{-3}}$. It was found after processing only a small fraction ($\sim$1%) of data from an ongoing all-sky pulsar survey. Follow-up observations have been made with the MWA, the upgraded Giant Metrewave Radio Telescope (uGMRT), and the Parkes 64 m telescopes, spanning a frequency range from $\sim$150 MHz to 4 GHz. The pulsar is faint, with an estimated flux density ($S$) of $\sim$1 mJy at 400 MHz and a spectrum $S(ν)\,\propto\,ν^{-2.0 \pm 0.2}$, where $ν$ is frequency. The DM-derived distance implies that it is also a low-luminosity source ($\sim$ 0.1 ${\rm mJy\,kpc^2}$ at 1400 MHz). The analysis of archival MWA observations reveals that the pulsar's mean flux density varies by up to a factor of $\sim$5-6 on timescales of several weeks to months. By combining MWA and uGMRT data, the pulsar position was determined to arcsecond precision. We also report on polarization properties detected in the MWA and Parkes bands. The pulsar's nondetection in previous pulsar and continuum imaging surveys, the observed high variability, and its detection in a small fraction of the survey data searched to date, all hint at a larger population of pulsars that await discovery in the southern hemisphere, with the MWA and the future low-frequency Square Kilometre Array. △ Less

Submitted 20 April, 2021; originally announced April 2021.

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

A Southern-Hemisphere all-sky radio transient monitor for SKA-Low prototype stations

Authors: M. Sokolowski, R. B. Wayth, N. D. R. Bhat, D. Price, J. W. Broderick, G. Bernardi, P. Bolli, R. Chiello, G. Comoretto, B. Crosse, D. B. Davidson, G. Macario, A. Magro, A. Mattana, D. Minchin, A. McPhail, J. Monari, F. Perini, G. Pupillo, G. Sleap, S. Tingay, D. Ung, A. Williams

Abstract: We present the first southern-hemisphere all-sky imager and radio-transient monitoring system implemented on two prototype stations of the low-frequency component of the Square Kilometre Array. Since its deployment the system has been used for real-time monitoring of the recorded commissioning data. Additionally, a transient searching algorithm has been executed on the resulting all-sky images. It… ▽ More We present the first southern-hemisphere all-sky imager and radio-transient monitoring system implemented on two prototype stations of the low-frequency component of the Square Kilometre Array. Since its deployment the system has been used for real-time monitoring of the recorded commissioning data. Additionally, a transient searching algorithm has been executed on the resulting all-sky images. It uses a difference imaging technique, and has enabled identification of a wide variety of transient classes, ranging from human-made radio-frequency interference to genuine astrophysical events. Observations at the frequency 159.4 MHz and higher in a single coarse channel (0.926 MHz) were made with 2s time resolution, and multiple nights were analysed. Despite having modest sensitivity (~few Jy/beam), using a single coarse channel and 2-s imaging, the system detected bright transients from PSR B0950+08, proving that it can be used to detect bright transients of an astrophysical origin. The unusual, extreme activity of the pulsar PSR B0950+08 (up to ~155 Jy/beam) was initially detected in a "blind" search in the 2020-04-10/11 data and later assigned to this specific pulsar. The limitations of our data, however, prevent use from making firm conclusions of the effect being due to a combination of refractive and diffractive scintillation or intrinsic emission mechanisms. The system can routinely collect data over many days without interruptions; the large amount of recorded data at 159.4 and 229.7 MHz allowed us to determine a preliminary transient surface density upper limit of $1.32 \times 10^{-9} \text{deg}^{-2}$ for a timescale and limiting flux density of 2s and 42 Jy, respectively. We plan to extend the observing bandwidth to tens of MHz and improve time resolution to tens of milliseconds in order to increase the sensitivity and enable detections of Fast Radio Bursts below 300 MHz. △ Less

Submitted 7 April, 2021; originally announced April 2021.

Comments: 22 pages, 11 figures, 3 tables, Accepted for publication in PASA

The Fermi GBM Gamma-Ray Burst Spectral Catalog: 10 Years of Data

Authors: S. Poolakkil, R. Preece, C. Fletcher, A. Goldstein, P. N. Bhat, E. Bissaldi, M. S. Briggs, E. Burns, W. H. Cleveland, M. M. Giles, C. M. Hui, D. Kocevski, S. Lesage, B. Mailyan, C. Malacaria, W. S. Paciesas, O. J. Roberts, P. Veres, A. von Kienlin, C. A. Wilson-Hodge

Abstract: We present the systematic spectral analyses of gamma-ray bursts (GRBs) detected by the Fermi Gamma-Ray Burst Monitor (GBM) during its first ten years of operation. This catalog contains two types of spectra; time-integrated spectral fits and spectral fits at the brightest time bin, from 2297 GRBs, resulting in a compendium of over 18000 spectra. The four different spectral models used for fitting… ▽ More We present the systematic spectral analyses of gamma-ray bursts (GRBs) detected by the Fermi Gamma-Ray Burst Monitor (GBM) during its first ten years of operation. This catalog contains two types of spectra; time-integrated spectral fits and spectral fits at the brightest time bin, from 2297 GRBs, resulting in a compendium of over 18000 spectra. The four different spectral models used for fitting the spectra were selected based on their empirical importance to the shape of many GRBs. We describe in detail our procedure and criteria for the analyses, and present the bulk results in the form of parameter distributions both in the observer frame and in the GRB rest frame. 941 GRBs from the first four years have been re-fitted using the same methodology as that of the 1356 GRBs in years five through ten. The data files containing the complete results are available from the High-Energy Astrophysics Science Archive Research Center (HEASARC). △ Less

Submitted 24 March, 2021; originally announced March 2021.

Rapid Spectral Variability of a Giant Flare from a Magnetar in NGC 253

Authors: O. J. Roberts, P. Veres, M. G. Baring, M. S. Briggs, C. Kouveliotou, E. Bissaldi, G. Younes, S. I. Chastain, J. J. DeLaunay, D. Huppenkothen, A. Tohuvavohu, P. N. Bhat, E. Gogus, A. J. van der Horst, J. A. Kennea, D. Kocevski, J. D. Linford, S. Guiriec, R. Hamburg, C. A. Wilson-Hodge, E. Burns

Abstract: Magnetars are slowly-rotating neutron stars with extremely strong magnetic fields ($10^{13-15}$ G), episodically emitting $\sim100$ ms long X-ray bursts with energies of $\sim10^{40-41}$ erg. Rarely, they produce extremely bright, energetic giant flares that begin with a short ($\sim0.2$ s), intense flash, followed by fainter, longer lasting emission modulated by the magnetar spin period (typicall… ▽ More Magnetars are slowly-rotating neutron stars with extremely strong magnetic fields ($10^{13-15}$ G), episodically emitting $\sim100$ ms long X-ray bursts with energies of $\sim10^{40-41}$ erg. Rarely, they produce extremely bright, energetic giant flares that begin with a short ($\sim0.2$ s), intense flash, followed by fainter, longer lasting emission modulated by the magnetar spin period (typically 2-12 s), thus confirming their origin. Over the last 40 years, only three such flares have been observed in our local group; they all suffered from instrumental saturation due to their extreme intensity. It has been proposed that extra-galactic giant flares likely constitute a subset of short gamma-ray bursts, noting that the sensitivity of current instrumentation prevents us from detecting the pulsating tail, while the initial bright flash is readily observable out to distances $\sim 10-20$ Mpc. Here, we report X- and gamma-ray observations of GRB 200415A, which exhibits a rapid onset, very fast time variability, flat spectra and significant sub-millisecond spectral evolution. These attributes match well with those expected for a giant flare from an extra-galactic magnetar, noting that GRB 200415A is directionally associated with the galaxy NGC 253 ($\sim$3.5 Mpc away). The detection of $\sim3$ MeV photons provides definitive evidence for relativistic motion of the emitting plasma. The observed rapid spectral evolution can naturally be generated by radiation emanating from such rapidly-moving gas in a rotating magnetar. △ Less

Submitted 13 January, 2021; originally announced January 2021.

Identifying and mitigating noise sources in precision pulsar timing data sets

Authors: Boris Goncharov, D. J. Reardon, R. M. Shannon, Xing-Jiang Zhu, Eric Thrane, M. Bailes, N. D. R. Bhat, S. Dai, G. Hobbs, M. Kerr, R. N. Manchester, S. Osłowski, A. Parthasarathy, C. J. Russell, R. Spiewak, N. Thyagarajan, J. B. Wang

Abstract: Pulsar timing array projects measure the pulse arrival times of millisecond pulsars for the primary purpose of detecting nanohertz-frequency gravitational waves. The measurements include contributions from a number of astrophysical and instrumental processes, which can either be deterministic or stochastic. It is necessary to develop robust statistical and physical models for these noise processes… ▽ More Pulsar timing array projects measure the pulse arrival times of millisecond pulsars for the primary purpose of detecting nanohertz-frequency gravitational waves. The measurements include contributions from a number of astrophysical and instrumental processes, which can either be deterministic or stochastic. It is necessary to develop robust statistical and physical models for these noise processes because incorrect models diminish sensitivity and may cause a spurious gravitational wave detection. Here we characterise noise processes for the 26 pulsars in the second data release of the Parkes Pulsar Timing Array using Bayesian inference. In addition to well-studied noise sources found previously in pulsar timing array data sets such as achromatic timing noise and dispersion measure variations, we identify new noise sources including time-correlated chromatic noise that we attribute to variations in pulse scattering. We also identify "exponential dip" events in four pulsars, which we attribute to magnetospheric effects as evidenced by pulse profile shape changes observed for three of the pulsars. This includes an event in PSR J1713$+$0747, which had previously been attributed to interstellar propagation. We present noise models to be used in searches for gravitational waves. We outline a robust methodology to evaluate the performance of noise models and identify unknown signals in the data. The detection of variations in pulse profiles highlights the need to develop efficient profile domain timing methods. △ Less

Submitted 1 November, 2020; v1 submitted 12 October, 2020; originally announced October 2020.

Comments: 18 pages, 7 figures, 5 tables

Precision orbital dynamics from interstellar scintillation arcs for PSR J0437-4715

Authors: Daniel J. Reardon, William A. Coles, Matthew Bailes, N. D. Ramesh Bhat, Shi Dai, George B. Hobbs, Matthew Kerr, Richard N. Manchester, Stefan Oslowski, Aditya Parthasarathy, Christopher J. Russell, Ryan M. Shannon, Renee Spiewak, Lawrence Toomey, Artem V. Tuntsov, Willem van Straten, Mark A. Walker, Jingbo Wang, Lei Zhang, Xing-Jiang Zhu

Abstract: Intensity scintillations of radio pulsars are known to originate from interference between waves scattered by the electron density irregularities of interstellar plasma, often leading to parabolic arcs in the two-dimensional power spectrum of the recorded dynamic spectrum. The degree of arc curvature depends on the distance to the scattering plasma and its transverse velocity with respect to the l… ▽ More Intensity scintillations of radio pulsars are known to originate from interference between waves scattered by the electron density irregularities of interstellar plasma, often leading to parabolic arcs in the two-dimensional power spectrum of the recorded dynamic spectrum. The degree of arc curvature depends on the distance to the scattering plasma and its transverse velocity with respect to the line-of-sight. We report the observation of annual and orbital variations in the curvature of scintillation arcs over a period of 16 years for the bright millisecond pulsar, PSR J0437-4715. These variations are the signature of the relative transverse motions of the Earth, pulsar, and scattering medium, which we model to obtain precise measurements of parameters of the pulsar's binary orbit and the scattering medium itself. We observe two clear scintillation arcs in most of our $>$5000 observations and we show that they originate from scattering by thin screens located at distances $D_1 = 89.8 \pm 0.4$ pc and $D_2 = 124 \pm 3$ pc from Earth. The best-fit scattering model we derive for the brightest arc yields the pulsar's orbital inclination angle $i = 137.1 \pm 0.3^\circ$, and longitude of ascending node, $Ω=206.3\pm0.4^\circ$. Using scintillation arcs for precise astrometry and orbital dynamics can be superior to modelling variations in the diffractive scintillation timescale, because the arc curvature is independent of variations in the level of turbulence of interstellar plasma. This technique can be used in combination with pulsar timing to determine the full three-dimensional orbital geometries of binary pulsars, and provides parameters essential for testing theories of gravity and constraining neutron star masses. △ Less

Submitted 27 September, 2020; originally announced September 2020.

Comments: 19 pages, 10 figures. Accepted for publication in ApJ

MWA tied-array processing III: Microsecond time resolution via a polyphase synthesis filter

Authors: S. J. McSweeney, S. M. Ord, N. D. R. Bhat, B. W. Meyers, S. E. Tremblay, J. Jones, B. Crosse, K. R. Smith

Abstract: A new high time resolution observing mode for the Murchison Widefield Array (MWA) is described, enabling full polarimetric observations with up to 30.72 MHz of bandwidth and a time resolution of ~0.8 $μ$s. This mode makes use of a polyphase synthesis filter to "undo" the polyphase analysis filter stage of the standard MWA's Voltage Capture System (VCS) observing mode. Sources of potential error in… ▽ More A new high time resolution observing mode for the Murchison Widefield Array (MWA) is described, enabling full polarimetric observations with up to 30.72 MHz of bandwidth and a time resolution of ~0.8 $μ$s. This mode makes use of a polyphase synthesis filter to "undo" the polyphase analysis filter stage of the standard MWA's Voltage Capture System (VCS) observing mode. Sources of potential error in the reconstruction of the high time resolution data are identified and quantified, with the S/N loss induced by the back-to-back system not exceeding -0.65 dB for typical noise-dominated samples. The system is further verified by observing three pulsars with known structure on microsecond timescales. △ Less

Submitted 6 July, 2020; originally announced July 2020.

Comments: Accepted for publication in PASA

Observation of inverse Compton emission from a long $γ$-ray burst

Authors: V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, D. Baack, A. Babić, B. Banerjee, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Ž. Bošnjak, G. Busetto, R. Carosi, G. Ceribella, Y. Chai , et al. (279 additional authors not shown)

Abstract: Long-duration gamma-ray bursts (GRBs) originate from ultra-relativistic jets launched from the collapsing cores of dying massive stars. They are characterised by an initial phase of bright and highly variable radiation in the keV-MeV band that is likely produced within the jet and lasts from milliseconds to minutes, known as the prompt emission. Subsequently, the interaction of the jet with the ex… ▽ More Long-duration gamma-ray bursts (GRBs) originate from ultra-relativistic jets launched from the collapsing cores of dying massive stars. They are characterised by an initial phase of bright and highly variable radiation in the keV-MeV band that is likely produced within the jet and lasts from milliseconds to minutes, known as the prompt emission. Subsequently, the interaction of the jet with the external medium generates external shock waves, responsible for the afterglow emission, which lasts from days to months, and occurs over a broad energy range, from the radio to the GeV bands. The afterglow emission is generally well explained as synchrotron radiation by electrons accelerated at the external shock. Recently, an intense, long-lasting emission between 0.2 and 1 TeV was observed from the GRB 190114C. Here we present the results of our multi-frequency observational campaign of GRB~190114C, and study the evolution in time of the GRB emission across 17 orders of magnitude in energy, from $5\times10^{-6}$ up to $10^{12}$\,eV. We find that the broadband spectral energy distribution is double-peaked, with the TeV emission constituting a distinct spectral component that has power comparable to the synchrotron component. This component is associated with the afterglow, and is satisfactorily explained by inverse Compton upscattering of synchrotron photons by high-energy electrons. We find that the conditions required to account for the observed TeV component are not atypical, supporting the possibility that inverse Compton emission is commonly produced in GRBs. △ Less

Submitted 12 June, 2020; originally announced June 2020.

Journal ref: Nature 575 (2019) 459-463

The MeerKAT Telescope as a Pulsar Facility: System verification and early science results from MeerTime

Authors: M. Bailes, A. Jameson, F. Abbate, E. D. Barr, N. D. R. Bhat, L. Bondonneau, M. Burgay, S. J. Buchner, F. Camilo, D. J. Champion, I. Cognard, P. B. Demorest, P. C. C. Freire, T. Gautam, M. Geyer, J. M. Griessmeier, L. Guillemot, H. Hu, F. Jankowski, S. Johnston, A. Karastergiou, R. Karuppusamy, D. Kaur, M. J. Keith, M. Kramer , et al. (50 additional authors not shown)

Abstract: We describe system verification tests and early science results from the pulsar processor (PTUSE) developed for the newly-commissioned 64-dish SARAO MeerKAT radio telescope in South Africa. MeerKAT is a high-gain (~2.8 K/Jy) low-system temperature (~18 K at 20cm) radio array that currently operates from 580-1670 MHz and can produce tied-array beams suitable for pulsar observations. This paper pres… ▽ More We describe system verification tests and early science results from the pulsar processor (PTUSE) developed for the newly-commissioned 64-dish SARAO MeerKAT radio telescope in South Africa. MeerKAT is a high-gain (~2.8 K/Jy) low-system temperature (~18 K at 20cm) radio array that currently operates from 580-1670 MHz and can produce tied-array beams suitable for pulsar observations. This paper presents results from the MeerTime Large Survey Project and commissioning tests with PTUSE. Highlights include observations of the double pulsar J0737-3039A, pulse profiles from 34 millisecond pulsars from a single 2.5h observation of the Globular cluster Terzan 5, the rotation measure of Ter5O, a 420-sigma giant pulse from the Large Magellanic Cloud pulsar PSR J0540-6919, and nulling identified in the slow pulsar PSR J0633-2015. One of the key design specifications for MeerKAT was absolute timing errors of less than 5 ns using their novel precise time system. Our timing of two bright millisecond pulsars confirm that MeerKAT delivers exceptional timing. PSR J2241-5236 exhibits a jitter limit of <4 ns per hour whilst timing of PSR J1909-3744 over almost 11 months yields an rms residual of 66 ns with only 4 min integrations. Our results confirm that the MeerKAT is an exceptional pulsar telescope. The array can be split into four separate sub-arrays to time over 1000 pulsars per day and the future deployment of S-band (1750-3500 MHz) receivers will further enhance its capabilities. △ Less

Submitted 28 May, 2020; originally announced May 2020.

Comments: 20 pages, 16 Figures, 4 Tables, accepted for publication in PASA