Showing 1–50 of 54 results for author: Mickaliger, M

Measuring glitch recoveries and braking indices with Bayesian model selection

Authors: Yang Liu, Michael J. Keith, Danai Antonopoulou, Patrick Weltevrede, Benjamin Shaw, Benjamin W. Stappers, Andrew G. Lyne, Mitchell B. Mickaliger, Avishek Basu

Abstract: For a selection of 35 pulsars with large spin-up glitches ($Δν/ν\geq10^{-6}$), which are monitored by the Jodrell Bank Observatory, we analyse 157 glitches and their recoveries. All parameters are measured consistently and we choose the best model to describe the post-glitch recovery based on Bayesian evidence. We present updated glitch epochs, sizes, changes of spin down rate, exponentially recov… ▽ More For a selection of 35 pulsars with large spin-up glitches ($Δν/ν\geq10^{-6}$), which are monitored by the Jodrell Bank Observatory, we analyse 157 glitches and their recoveries. All parameters are measured consistently and we choose the best model to describe the post-glitch recovery based on Bayesian evidence. We present updated glitch epochs, sizes, changes of spin down rate, exponentially recovering components (amplitude and corresponding timescale) when present, as well as pulsars' second frequency derivatives and their glitch associated changes if detected. We discuss the different observed styles of post-glitch recovery as well as some particularly interesting sources. Several correlations are revealed between glitch parameters and pulsar spin parameters, including a very strong correlation between a pulsar's interglitch $|\ddotν|$ and $\dotν$, as well as between the glitch-induced spin-down rate change $Δ\dotν_{\rm p}$ that does not relax exponentially and $\dotν$. We find that the ratio $\left|Δ\dotν_{\mathrm{p}}/\ddotν\right|$ can be used as an estimate of glitch recurrence times, especially for those pulsars for which there are indications of a characteristic glitch size and interglitch waiting time. We calculate the interglitch braking index $n$ and find that pulsars with large glitches typically have $n$ greater than $3$, suggesting that internal torques dominate the rotational evolution between glitches. The external torque, e.g. from electromagnetic dipole radiation, could dominate the observed $\ddotν$ for the youngest pulsars ($\lesssim10^{4}\;\mathrm{yr}$), which may be expected to display $n\sim3$. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: 24 pages, 14 figures, 8 tables. Accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 532, Issue 1, pp.859-882 July 2024

Improving pulsar timing precision through superior Time-of-Arrival creation

Authors: J. Wang, J. P. W. Verbiest, G. M. Shaifullah, I. Cognard, L. Guillemot, G. H. Janssen, M. B. Mickaliger, A. Possenti, G. Theureau

Abstract: The measurement of pulsar pulse times-of-arrival (ToAs) is a crucial step in detecting low-frequency gravitational waves. To determine ToAs, we can use template-matching to compare each observed pulse profile with a standard template. However, using different combinations of templates and template-matching methods (TMMs) without careful consideration may lead to inconsistent results. In pulsar tim… ▽ More The measurement of pulsar pulse times-of-arrival (ToAs) is a crucial step in detecting low-frequency gravitational waves. To determine ToAs, we can use template-matching to compare each observed pulse profile with a standard template. However, using different combinations of templates and template-matching methods (TMMs) without careful consideration may lead to inconsistent results. In pulsar timing array (PTA) experiments, distinct ToAs from the same observations can be obtained, due to the use of diverse templates and TMMs. In other words, employing diverse approaches can yield different timing results and would thus have a significant impact on subsequent gravitational wave searches. In this paper, we examine several commonly used combinations to analyze their effect on pulse ToAs. we evaluate the potential impact of template and TMM selection on thirteen typical millisecond pulsars within the European PTA. We employ pulsar timing methods, specifically the root mean square and reduced chi-square $χ_r^2$ of the residuals of the best timing solution to assess the outcomes. Additionally, we evaluate the system-limited noise floor (SLNF) for each pulsar at various telescopes operating around 1.4~GHz using frequency-resolved templates. Our findings suggest that utilizing data-derived and smoothed templates in conjunction with the Fourier-domain with Markov-chain Monte Carlo (FDM) TMM is generally the most effective approach, though there may be exceptions that require further attention. Furthermore, we determine that pulse phase jitter noise does not significantly limit the current precision of the European PTA's timing, as jitter levels derived from other studies are much smaller than the SLNF. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: Accepted by A&A

Radio Pulse Profile Evolution of Magnetar Swift J1818.0-1607

Authors: Rebecca Fisher, Elliot Butterworth, Kaustubh Rajwade, Ben Stappers, Gregory Desvignes, Ramesh Karuppusamy, Michael Kramer, Kuo Liu, Andrew Lyne, Mitchell Mickaliger, Benjamin Shaw, Patrick Weltevrede

Abstract: The shape and polarisation properties of the radio pulse profiles of radio-loud magnetars provide a unique opportunity to investigate their magnetospheric properties. Gaussian Process Regression analysis was used to investigate the variation in the total intensity shape of the radio pulse profiles of the magnetar Swift J1818.0-1607. The observed profile shape was found to evolve through three mode… ▽ More The shape and polarisation properties of the radio pulse profiles of radio-loud magnetars provide a unique opportunity to investigate their magnetospheric properties. Gaussian Process Regression analysis was used to investigate the variation in the total intensity shape of the radio pulse profiles of the magnetar Swift J1818.0-1607. The observed profile shape was found to evolve through three modes between MJDs 59104 and 59365. The times at which these transitions occurred coincided with changes in the amplitude of modulations in the spin-down rate. The amount of linear and circular polarisation was also found to vary significantly with time. Lomb-Scargle periodogram analysis of the spin-down rate revealed three possibly harmonically related frequencies. This could point to the magnetar experiencing seismic activity. However, no profile features exhibited significant periodicity, suggesting no simple correlations between the profile variability and fluctuations of the spin-down on shorter timescales within the modes. Overall, this implies the mode changes seen are a result of local magnetospheric changes, with other theories, such as precession, less able to explain these observations. △ Less

Submitted 23 January, 2024; originally announced January 2024.

Comments: 11 pages, 6 figures, accepted for publication in MNRAS

A Cacophony of Echoes from daily monitoring of the Crab Pulsar at Jodrell Bank

Authors: Thierry Serafin Nadeau, Marten H. van Kerkwijk, Cees G. Bassa, Ben W. Stappers, Mitchell B. Mickaliger, Andrew G. Lyne

Abstract: Using archival data from the 42 foot telescope at the Jodrell Bank Observatory, we produce daily stacks of aligned giant pulses for the Crab pulsar, to study changes to the daily profiles between April 2012 to December 2016. From these, we identify echoes, where intervening material away from the line of sight causes pulsed emission to be redirected towards the observer, with delay corresponding t… ▽ More Using archival data from the 42 foot telescope at the Jodrell Bank Observatory, we produce daily stacks of aligned giant pulses for the Crab pulsar, to study changes to the daily profiles between April 2012 to December 2016. From these, we identify echoes, where intervening material away from the line of sight causes pulsed emission to be redirected towards the observer, with delay corresponding to the increased distance of travel, resulting in additional profile components. These observations show that such echoes may be far more common than implied by the previous rate of detections. All the observed echoes are consistent with approaching zero-delay at their closest approach to the normal giant pulse emission. This indicates that the structures responsible for producing these events must be highly anisotropic, with typical lengths greater than $\sim 4\textrm{AU}$, typical widths on the sky of $\sim 0.1 \textrm{AU}$ and typical depths of $\sim 5\textrm{AU}$, given the previously observed electron densities of the nebular filaments, on the order of 1000 cm$^{-3}$. This suggests that these inhomogeneities are likely to be offshoot substructure from the larger nebular filaments of the Crab nebula. △ Less

Submitted 18 December, 2023; v1 submitted 10 October, 2023; originally announced October 2023.

Comments: 11 pages, 7 figures, accepted for publication in The Astrophysical Journal

Comparing recent PTA results on the nanohertz stochastic gravitational wave background

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

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

Submitted 1 September, 2023; originally announced September 2023.

Comments: 21 pages, 9 figures, submitted to ApJ

The second data release from the European Pulsar Timing Array: VI. Challenging the ultralight dark matter paradigm

Authors: Clemente Smarra, Boris Goncharov, Enrico Barausse, J. Antoniadis, S. Babak, A. -S. Bak Nielsen, C. G. Bassa, A. Berthereau, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, A. Chalumeau, D. J. Champion, S. Chanlaridis, S. Chen, I. Cognard, G. Desvignes, M. Falxa, R. D. Ferdman, A. Franchini, J. R. Gair, E. Graikou, J. -M. Grie , et al. (46 additional authors not shown)

Abstract: Pulsar Timing Array experiments probe the presence of possible scalar or pseudoscalar ultralight dark matter particles through decade-long timing of an ensemble of galactic millisecond radio pulsars. With the second data release of the European Pulsar Timing Array, we focus on the most robust scenario, in which dark matter interacts only gravitationally with ordinary baryonic matter. Our results s… ▽ More Pulsar Timing Array experiments probe the presence of possible scalar or pseudoscalar ultralight dark matter particles through decade-long timing of an ensemble of galactic millisecond radio pulsars. With the second data release of the European Pulsar Timing Array, we focus on the most robust scenario, in which dark matter interacts only gravitationally with ordinary baryonic matter. Our results show that ultralight particles with masses $10^{-24.0}~\text{eV} \lesssim m \lesssim 10^{-23.3}~\text{eV}$ cannot constitute $100\%$ of the measured local dark matter density, but can have at most local density $ρ\lesssim 0.3$ GeV/cm$^3$. △ Less

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

Comments: 5 pages + acknowledgements + refs, 2 figures. Text and figures match the version published in PRL

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

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

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

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

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

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

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

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

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

Comments: 13 figures, 15 pages, accepted

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

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

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

Submitted 28 June, 2023; originally announced June 2023.

Comments: 20 pages, 6 figures, 9 tables

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

The second data release from the European Pulsar Timing Array I. The dataset and timing analysis

Authors: J. Antoniadis, S. Babak, A. -S. Bak Nielsen, C. G. Bassa, A. Berthereau, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, A. Chalumeau, D. J. Champion, S. Chanlaridis, S. Chen, I. Cognard, G. Desvignes, M. Falxa, R. D. Ferdman, A. Franchini, J. R. Gair, B. Goncharov, E. Graikou, J. -M. Grießmeier, L. Guillemot, Y. J. Guo , et al. (44 additional authors not shown)

Abstract: Pulsar timing arrays offer a probe of the low-frequency gravitational wave spectrum (1 - 100 nanohertz), which is intimately connected to a number of markers that can uniquely trace the formation and evolution of the Universe. We present the dataset and the results of the timing analysis from the second data release of the European Pulsar Timing Array (EPTA). The dataset contains high-precision pu… ▽ More Pulsar timing arrays offer a probe of the low-frequency gravitational wave spectrum (1 - 100 nanohertz), which is intimately connected to a number of markers that can uniquely trace the formation and evolution of the Universe. We present the dataset and the results of the timing analysis from the second data release of the European Pulsar Timing Array (EPTA). The dataset contains high-precision pulsar timing data from 25 millisecond pulsars collected with the five largest radio telescopes in Europe, as well as the Large European Array for Pulsars. The dataset forms the foundation for the search for gravitational waves by the EPTA, presented in associated papers. We describe the dataset and present the results of the frequentist and Bayesian pulsar timing analysis for individual millisecond pulsars that have been observed over the last ~25 years. We discuss the improvements to the individual pulsar parameter estimates, as well as new measurements of the physical properties of these pulsars and their companions. This data release extends the dataset from EPTA Data Release 1 up to the beginning of 2021, with individual pulsar datasets with timespans ranging from 14 to 25 years. These lead to improved constraints on annual parallaxes, secular variation of the orbital period, and Shapiro delay for a number of sources. Based on these results, we derived astrophysical parameters that include distances, transverse velocities, binary pulsar masses, and annual orbital parallaxes. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 29 pages, 9 figures, 13 tables, Astronomy & Astrophysics in press

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

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

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

Submitted 28 June, 2023; originally announced June 2023.

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

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

Variable Scintillation Arcs of Millisecond Pulsars observed with the Large European Array for Pulsars

Authors: R. A. Main, J. Antoniadis, S. Chen, I. Cognard, H. Hu, J. Jang, R. Karuppusamy, M. Kramer, K. Liu, Y. Liu, G. Mall, J. W. McKee, M. B. Mickaliger, D. Perrodin, S. A. Sanidas, B. W. Stappers, T. Sprenger, O. Wucknitz, C. G. Bassa, M. Burgay, R. Concu, M. Gaikwad, G. H. Janssen, K. J. Lee, A. Melis , et al. (4 additional authors not shown)

Abstract: We present the first large sample of scintillation arcs in millisecond pulsars, analysing 12 sources observed with the Large European Array for Pulsars (LEAP), and the Effelsberg 100\,m telescope. We estimate the delays from multipath propagation, measuring significant correlated changes in scattering timescales over a 10-year timespan. Many sources show compact concentrations of power in the seco… ▽ More We present the first large sample of scintillation arcs in millisecond pulsars, analysing 12 sources observed with the Large European Array for Pulsars (LEAP), and the Effelsberg 100\,m telescope. We estimate the delays from multipath propagation, measuring significant correlated changes in scattering timescales over a 10-year timespan. Many sources show compact concentrations of power in the secondary spectrum, which in PSRs J0613$-$0200 and J1600$-$3053 can be tracked between observations, and are consistent with compact scattering at fixed angular positions. Other sources such as PSRs J1643$-$1224 and J0621+1002 show diffuse, asymmetric arcs which are likely related to phase-gradients across the scattering screen. PSR B1937+21 shows at least three distinct screens which dominate at different times and evidence of varying screen axes or multi-screen interactions. We model annual and orbital arc curvature variations in PSR J0613$-$0200, providing a measurement of the longitude of ascending node, resolving the sense of the orbital inclination, where our best fit model is of a screen with variable axis of anisotropy over time, corresponding to changes in the scattering of the source. Unmodeled variations of the screen's axis of anisotropy are likely to be a limiting factor in determining orbital parameters with scintillation, requiring careful consideration of variable screen properties, or independent VLBI measurements. Long-term scintillation studies such as this serve as a complementary tool to pulsar timing, to measure a source of correlated noise for pulsar timing arrays, solve pulsar orbits, and to understand the astrophysical origin of scattering screens. △ Less

Submitted 23 June, 2023; originally announced June 2023.

Comments: Submitted to MNRAS, first revision. 18 pages, 16 figures

Practical approaches to analyzing PTA data: Cosmic strings with six pulsars

Authors: Hippolyte Quelquejay Leclere, Pierre Auclair, Stanislav Babak, Aurélien Chalumeau, Danièle A. Steer, J. Antoniadis, A. -S. Bak Nielsen, C. G. Bassa, A. Berthereau, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, D. J. Champion, S. Chanlaridis, S. Chen, I. Cognard, G. Desvignes, M. Falxa, R. D. Ferdman, A. Franchini, J. R. Gair, B. Goncharov, E. Graikou , et al. (47 additional authors not shown)

Abstract: We search for a stochastic gravitational wave background (SGWB) generated by a network of cosmic strings using six millisecond pulsars from Data Release 2 (DR2) of the European Pulsar Timing Array (EPTA). We perform a Bayesian analysis considering two models for the network of cosmic string loops, and compare it to a simple power-law model which is expected from the population of supermassive blac… ▽ More We search for a stochastic gravitational wave background (SGWB) generated by a network of cosmic strings using six millisecond pulsars from Data Release 2 (DR2) of the European Pulsar Timing Array (EPTA). We perform a Bayesian analysis considering two models for the network of cosmic string loops, and compare it to a simple power-law model which is expected from the population of supermassive black hole binaries. Our main strong assumption is that the previously reported common red noise process is a SGWB. We find that the one-parameter cosmic string model is slightly favored over a power-law model thanks to its simplicity. If we assume a two-component stochastic signal in the data (supermassive black hole binary population and the signal from cosmic strings), we get a $95\%$ upper limit on the string tension of $\log_{10}(Gμ) < -9.9$ ($-10.5$) for the two cosmic string models we consider. In extended two-parameter string models, we were unable to constrain the number of kinks. We test two approximate and fast Bayesian data analysis methods against the most rigorous analysis and find consistent results. These two fast and efficient methods are applicable to all SGWBs, independent of their source, and will be crucial for analysis of extended data sets. △ Less

Submitted 3 May, 2024; v1 submitted 21 June, 2023; originally announced June 2023.

Comments: 14 pages, 6 figures; typo corrected in (5)

Journal ref: Phys. Rev. D 108 (2023), 123527

A renewed search for radio emission from the variable $γ$-ray pulsar PSR J2021$+$4026

Authors: B. Shaw, B. W. Stappers, P. Weltevrede, C. A. Jordan, M. B. Mickaliger, A. G. Lyne

Abstract: We undertake the first targeted search at 1.5 GHz for radio emission from the variable $γ$-ray pulsar PSR J2021$+$4026. This radio-quiet pulsar assumes one of two stable $γ$-ray emission states, between which it transitions on a timescale of years. These transitions, in both $γ$-ray flux and pulse profile shape, are accompanied by contemporaneous changes to the pulsar's spin-down rate. A number of… ▽ More We undertake the first targeted search at 1.5 GHz for radio emission from the variable $γ$-ray pulsar PSR J2021$+$4026. This radio-quiet pulsar assumes one of two stable $γ$-ray emission states, between which it transitions on a timescale of years. These transitions, in both $γ$-ray flux and pulse profile shape, are accompanied by contemporaneous changes to the pulsar's spin-down rate. A number of radio pulsars are known to exhibit similar correlated variability, which in some cases involves an emission state in which the radio emission ceases to be detectable. In this paper, we perform a search for radio emission from PSR J2021$+$4026, using archival radio observations recorded when the pulsar was in each of its emission/spin-down states. Using improved techniques, we search for periodic radio emission as well as single pulse phenomena such as giant radio pulses and RRAT-like emission. Our search reveals no evidence of radio emission from PSR J2021$+$4026. We estimate that the flux density for periodic emission from PSR J2021$+$4026 does not exceed 0.2 mJy at this frequency. We also estimate single-pulse flux limits for RRAT-like bursts and giant radio pulses to be 0.3 and 100 Jy respectively. We discuss the transitioning behaviour of PSR J2021$+$4026 in the context of pulsar glitches, intermittent pulsars and the increasingly common emission-rotation correlation observed in radio pulsars. △ Less

Submitted 11 May, 2023; v1 submitted 9 May, 2023; originally announced May 2023.

Comments: 10 pages, 5 figures, 3 tables

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

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

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

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

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

Detection of quasi-periodic micro-structure in three millisecond pulsars with the Large European Array for Pulsars

Authors: K. Liu, J. Antoniadis, C. G. Bassa, S. Chen, I. Cognard, M. Gaikwad, H. Hu, J. Jang, G. H. Janssen, R. Karuppusamy, M. Kramer, K. J. Lee, R. A. Main, G. Mall, J. W. McKee, M. B. Mickaliger, D. Perrodin, S. A. Sanidas, B. W. Stappers, L. Wang, W. W. Zhu, M. Burgay, R. Concu, A. Corongiu, A. Melis , et al. (2 additional authors not shown)

Abstract: We report on the detection of quasi-periodic micro-structure in three millisecond pulsars (MSPs), PSRs J1022+1001, J2145-0750 and J1744-1134, using high time resolution data acquired with the Large European Array for Pulsars at a radio frequency of 1.4 GHz. The occurrence rate of quasi-periodic micro-structure is consistent among pulses with different peak flux densities. Using an auto-correlation… ▽ More We report on the detection of quasi-periodic micro-structure in three millisecond pulsars (MSPs), PSRs J1022+1001, J2145-0750 and J1744-1134, using high time resolution data acquired with the Large European Array for Pulsars at a radio frequency of 1.4 GHz. The occurrence rate of quasi-periodic micro-structure is consistent among pulses with different peak flux densities. Using an auto-correlation analysis, we measure the periodicity and width of the micro-structure in these three pulsars. The detected micro-structure from PSRs J1022+1001 and J1744-1134 is often highly linearly polarised. In PSR J1022+1001, the linear polarisation position angles of micro-structure pulses are in general flat with a small degree of variation. Using these results, we further examine the frequency and rotational period dependency of micro-structure properties established in previous work, along with the angular beaming and temporal modulation models that explains the appearance of micro-structure. We also discuss a possible link of micro-structure to the properties of some of the recently discovered fast radio bursts which exhibit a very similar emission morphology. △ Less

Submitted 20 June, 2022; originally announced June 2022.

Comments: 8 pages, 6 figures, 1 table, published in MNRAS

Long-term scintillation studies of EPTA pulsars. I. Observations and basic results

Authors: Yulan Liu, Joris P. W. Verbiest, Robert A. Main, Ziwei Wu, Krishnakumar Moochickal Ambalappat, David J. Champion, Ismaël Cognard, Lucas Guillemot, Madhuri Gaikwad, Gemma H. Janssen, Michael Kramer, Michael J. Keith, Ramesh Karuppusamy, Lars Künkel, Kuo Liu, James W. McKee, Mitchell B. Mickaliger, Ben W. Stappers, Golam. M. Shaifullah, Gilles Theureau

Abstract: Interstellar scintillation analysis of pulsars allows us to probe the small-scale distribution and inhomogeneities of the ionized interstellar medium. Our priority is to present the data set and the basic measurements of scintillation parameters of pulsars employing long-term scintillation observations carried out from 2011 January to 2020 August by the European Pulsar Timing Array radio telesco… ▽ More Interstellar scintillation analysis of pulsars allows us to probe the small-scale distribution and inhomogeneities of the ionized interstellar medium. Our priority is to present the data set and the basic measurements of scintillation parameters of pulsars employing long-term scintillation observations carried out from 2011 January to 2020 August by the European Pulsar Timing Array radio telescopes in the 21-cm and 11-cm bands. Additionally, we aim to identify future possible lines of study using this long-term scintillation dataset. We present the long-term time series of $ν_{\rm d}$ and $τ_{\rm d}$ for 13 pulsars. Sanity-checks and comparisons indicate that the scintillation parameters of our work and previously published works are mostly consistent. For two pulsars, PSRs~J1857+0943 and J1939+2134, we were able to obtain measurements of the $ν_{\rm d}$ at both bands, which allows us to derive the time series of frequency scaling indices with a mean and a standard deviation of 2.82$\pm$1.95 and 3.18$\pm$0.60, respectively. We found some interesting features which will be studied in more detail in subsequent papers in this series: (i) in the time series of PSR~J1939+2134, where the scintillation bandwidth sharply increases or decreases associated with a sharp change of dispersion measure; (ii) PSR~J0613$-$0200 and PSR~J0636+5126 show a strong annual variation in the time series of the $τ_{\rm d}$; (iii) PSR~J1939+2134 shows a weak anti-correlation between scintillation timescale and dispersion in WSRT data. △ Less

Submitted 13 April, 2022; v1 submitted 31 March, 2022; originally announced March 2022.

Comments: 18 pages, 6 figures, 4 tables, accepted for publication in A&A

Journal ref: A&A 664, A116 (2022)

A search for planetary companions around 800 pulsars from the Jodrell Bank pulsar timing programme

Authors: Iuliana C. Niţu, Michael J. Keith, Ben W. Stappers, Andrew G. Lyne, Mitchell B. Mickaliger

Abstract: We have searched for planetary companions around 800 pulsars monitored at the Jodrell Bank Observatory, with both circular and eccentric orbits of periods between $20$ days and $17$ years and inclination-dependent planetary masses from $10^{-4}$ to $100\,\mathrm{M}_{\oplus}$. Using a Bayesian framework, we simultaneously model pulsar timing parameters and a stationary noise process with a power-la… ▽ More We have searched for planetary companions around 800 pulsars monitored at the Jodrell Bank Observatory, with both circular and eccentric orbits of periods between $20$ days and $17$ years and inclination-dependent planetary masses from $10^{-4}$ to $100\,\mathrm{M}_{\oplus}$. Using a Bayesian framework, we simultaneously model pulsar timing parameters and a stationary noise process with a power-law power spectral density. We put limits on the projected masses of any planetary companions, which reach as low as 1/100th of the mass of the Moon ($\sim 10^{-4}\,\mathrm{M}_{\oplus}$). We find that two-thirds of our pulsars are highly unlikely to host any companions above $2-8\,\mathrm{M}_{\oplus}$. Our results imply that fewer than $0.5\%$ of pulsars could host terrestrial planets as large as those known to orbit PSR B1257$+$12 ($\sim4\,\mathrm{M}_{\oplus}$); however, the smaller planet in this system ($\sim0.02\,\mathrm{M}_{\oplus}$) would be undetectable in $95\%$ of our sample, hidden by both instrumental and intrinsic noise processes, although it is not clear if such tiny planets could exist in isolation. We detect significant periodicities in 15 pulsars, however we find that intrinsic quasi-periodic magnetospheric effects can mimic the influence of a planet, and for the majority of these cases we believe this to be the origin of the detected periodicity. Notably, we find that the highly periodic oscillations in PSR B0144$+$59 are correlated with changes in the pulse profile and therefore can be attributed to magnetospheric effects. We believe the most plausible candidate for planetary companions in our sample is PSR J2007$+$3120. △ Less

Submitted 2 March, 2022; originally announced March 2022.

Comments: Accepted for publication in MNRAS

Long term radio and X-ray evolution of the magnetar Swift J1818.0-1607

Authors: K. M. Rajwade, B. W. Stappers, A. G. Lyne, B. Shaw, M. B. Mickaliger, K. Liu, M. Kramer, G. Desvignes, R. Karuppusamy, T. Enoto, T. Güver, Chin-Ping Hu, M. P. Surnis

Abstract: We report on the the long term monitoring campaign of the seemingly youngest magnetar Swift~J1818.0-1607 at radio and X-ray wavelengths over a span of one year. We obtained a coherent timing solution for the magnetar over the same time span. The frequency derivative of the magnetar shows systematic variation with the values oscillating about a mean value of $-$1.37$\times$10$^{-11}$ Hz s$^{-1}$. T… ▽ More We report on the the long term monitoring campaign of the seemingly youngest magnetar Swift~J1818.0-1607 at radio and X-ray wavelengths over a span of one year. We obtained a coherent timing solution for the magnetar over the same time span. The frequency derivative of the magnetar shows systematic variation with the values oscillating about a mean value of $-$1.37$\times$10$^{-11}$ Hz s$^{-1}$. The magnitude of the variation in the frequency derivative reduces with time before converging on the mean value. We were able to identify four states in the spin-frequency derivative that were quantified by the amount of modulation about the mean value and the transition between these states seem to be correlated with the change in the radio emission of the magnetar while no correlation is seen in the average radio profile variability on a shorter timescale (days). The 0.5--12 keV X-ray flux shows a monotonic decrease that can be attributed to thermal emission from a hot-spot on the surface of the neutron star that is reducing in size. Such decrease is consistent with what is seen in other magnetars. The potential correlation between the radio emission mode and the behaviour of the spin-down rate hints to a global change in the magnetopshere of the magnetar akin to the correlation seen in a subset of mode-changing radio pulsars and suggests a physical link between the two sub-populations. △ Less

Submitted 15 February, 2022; originally announced February 2022.

Comments: 9 pages, 6 figures, 2 tables, supplementary material. Accepted for publication in MNRAS

Discovery, Timing, and Multiwavelength Observations of the Black Widow Millisecond Pulsar PSR J1555-2908

Authors: Paul S. Ray, Lars Nieder, Colin J. Clark, Scott M. Ransom, H. Thankful Cromartie, Dale A. Frail, Kunal P. Mooley, Huib Intema, Preshanth Jagannathan, Paul Demorest, Kevin Stovall, Jules P. Halpern, Julia Deneva, Sebastien Guillot, Matthew Kerr, Samuel J. Swihart, Philippe Bruel, Ben W. Stappers, Andrew Lyne, Mitch Mickaliger, Fernando Camilo, Elizabeth C. Ferrara, Michael T. Wolff, P. F. Michelson

Abstract: We report the discovery of PSR J1555-2908, a 1.79 ms radio and gamma-ray pulsar in a 5.6 hr binary system with a minimum companion mass of 0.052 $M_\odot$. This fast and energetic ($\dot E = 3 \times 10^{35}$ erg/s) millisecond pulsar was first detected as a gamma-ray point source in Fermi LAT sky survey observations. Guided by a steep spectrum radio point source in the Fermi error region, we perf… ▽ More We report the discovery of PSR J1555-2908, a 1.79 ms radio and gamma-ray pulsar in a 5.6 hr binary system with a minimum companion mass of 0.052 $M_\odot$. This fast and energetic ($\dot E = 3 \times 10^{35}$ erg/s) millisecond pulsar was first detected as a gamma-ray point source in Fermi LAT sky survey observations. Guided by a steep spectrum radio point source in the Fermi error region, we performed a search at 820 MHz with the Green Bank Telescope that first discovered the pulsations. The initial radio pulse timing observations provided enough information to seed a search for gamma-ray pulsations in the LAT data, from which we derive a timing solution valid for the full Fermi mission. In addition to the radio and gamma-ray pulsation discovery and timing, we searched for X-ray pulsations using NICER but no significant pulsations were detected. We also obtained time-series r-band photometry that indicates strong heating of the companion star by the pulsar wind. Material blown off the heated companion eclipses the 820 MHz radio pulse during inferior conjunction of the companion for ~10% of the orbit, which is twice the angle subtended by its Roche lobe in an edge-on system. △ Less

Submitted 9 February, 2022; originally announced February 2022.

Comments: 15 pages, 6 figures, accepted by ApJ

Modelling annual scintillation arc variations in PSR J1643-1224 using the Large European Array for Pulsars

Authors: G. Mall, R. A. Main, J. Antoniadis, C. G. Bassa, M. Burgay, S. Chen, I. Cognard, R. Concu, A. Corongiu, M. Gaikwad, H. Hu, G. H. Janssen, R. Karuppusamy, K. J. Lee, K. Liu, J. W. McKee, A. Melis, M. B. Mickaliger, D. Perrodin, M. Pilia, A. Possenti, D. J. Reardon, S. A. Sanidas, T. Sprenger, B. W. Stappers , et al. (3 additional authors not shown)

Abstract: In this work we study variations in the parabolic scintillation arcs of the binary millisecond pulsar PSR J1643-1224 over five years using the Large European Array for Pulsars (LEAP). The 2D power spectrum of scintillation, called the secondary spectrum, often shows a parabolic distribution of power, where the arc curvature encodes the relative velocities and distances of the pulsar, ionised inter… ▽ More In this work we study variations in the parabolic scintillation arcs of the binary millisecond pulsar PSR J1643-1224 over five years using the Large European Array for Pulsars (LEAP). The 2D power spectrum of scintillation, called the secondary spectrum, often shows a parabolic distribution of power, where the arc curvature encodes the relative velocities and distances of the pulsar, ionised interstellar medium (IISM), and Earth. We observe a clear parabolic scintillation arc which varies in curvature throughout the year. The distribution of power in the secondary spectra are inconsistent with a single scattering screen which is fully 1D, or entirely isotropic. We fit the observed arc curvature variations with two models; an isotropic scattering screen, and a model with two independent 1D screens. We measure the distance to the scattering screen to be in the range 114-223 pc, depending on the model, consistent with the known distance of the foreground large-diameter HII region Sh 2-27 (112+/-17 pc), suggesting that it is the dominant source of scattering. We obtain only weak constraints on the pulsar's orbital inclination and angle of periastron, since the scintillation pattern is not very sensitive to the pulsar's motion, since the screen is much closer to the Earth than the pulsar. More measurements of this kind - where scattering screens can be associated with foreground objects - will help to inform the origins and distribution of scattering screens within our galaxy. △ Less

Submitted 11 January, 2022; originally announced January 2022.

Comments: 11 pages, 6 figures. Accepted by MNRAS

The International Pulsar Timing Array second data release: Search for an isotropic Gravitational Wave Background

Authors: J. Antoniadis, Z. Arzoumanian, S. Babak, M. Bailes, A. -S. Bak Nielsen, P. T. Baker, C. G. Bassa, B. Becsy, A. Berthereau, M. Bonetti, A. Brazier, P. R. Brook, M. Burgay, S. Burke-Spolaor, R. N. Caballero, J. A. Casey-Clyde, A. Chalumeau, D. J. Champion, M. Charisi, S. Chatterjee, S. Chen, I. Cognard, J. M. Cordes, N. J. Cornish, F. Crawford , et al. (101 additional authors not shown)

Abstract: We searched for an isotropic stochastic gravitational wave background in the second data release of the International Pulsar Timing Array, a global collaboration synthesizing decadal-length pulsar-timing campaigns in North America, Europe, and Australia. In our reference search for a power law strain spectrum of the form $h_c = A(f/1\,\mathrm{yr}^{-1})^α$, we found strong evidence for a spectrally… ▽ More We searched for an isotropic stochastic gravitational wave background in the second data release of the International Pulsar Timing Array, a global collaboration synthesizing decadal-length pulsar-timing campaigns in North America, Europe, and Australia. In our reference search for a power law strain spectrum of the form $h_c = A(f/1\,\mathrm{yr}^{-1})^α$, we found strong evidence for a spectrally-similar low-frequency stochastic process of amplitude $A = 3.8^{+6.3}_{-2.5}\times10^{-15}$ and spectral index $α= -0.5 \pm 0.5$, where the uncertainties represent 95\% credible regions, using information from the auto- and cross-correlation terms between the pulsars in the array. For a spectral index of $α= -2/3$, as expected from a population of inspiralling supermassive black hole binaries, the recovered amplitude is $A = 2.8^{+1.2}_{-0.8}\times10^{-15}$. Nonetheless, no significant evidence of the Hellings-Downs correlations that would indicate a gravitational-wave origin was found. We also analyzed the constituent data from the individual pulsar timing arrays in a consistent way, and clearly demonstrate that the combined international data set is more sensitive. Furthermore, we demonstrate that this combined data set produces comparable constraints to recent single-array data sets which have more data than the constituent parts of the combination. Future international data releases will deliver increased sensitivity to gravitational wave radiation, and significantly increase the detection probability. △ Less

Submitted 11 January, 2022; originally announced January 2022.

Comments: 17 pages, 12 figures, accepted in MNRAS

A MeerKAT, e-MERLIN, H.E.S.S. and Swift search for persistent and transient emission associated with three localised FRBs

Authors: James O. Chibueze, M. Caleb, L. Spitler, H. Ashkar, F. Schussler, B. W. Stappers, C. Venter, I. Heywood, A. M. S. Richards, D. R. A. Williams, M. Kramer, R. Beswick, M. C. Bezuidenhout, R. P. Breton, L. N. Driessen, F. Jankowski, E. F. Keane, M. Malenta, M. Mickaliger, V. Morello, H. Qiu, K. Rajwade, S. Sanidas, M. Surnis, T. W. Scragg , et al. (134 additional authors not shown)

Abstract: We report on a search for persistent radio emission from the one-off Fast Radio Burst (FRB) 20190714A, as well as from two repeating FRBs, 20190711A and 20171019A, using the MeerKAT radio telescope. For FRB 20171019A we also conducted simultaneous observations with the High Energy Stereoscopic System (H.E.S.S.) in very high energy gamma rays and searched for signals in the ultraviolet, optical, an… ▽ More We report on a search for persistent radio emission from the one-off Fast Radio Burst (FRB) 20190714A, as well as from two repeating FRBs, 20190711A and 20171019A, using the MeerKAT radio telescope. For FRB 20171019A we also conducted simultaneous observations with the High Energy Stereoscopic System (H.E.S.S.) in very high energy gamma rays and searched for signals in the ultraviolet, optical, and X-ray bands. For this FRB, we obtain a UV flux upper limit of 1.39x10^-16 erg/cm^-2/s/Amstrong, X-ray limit of ~ 6.6x10^-14 erg/cm^-2/s and a limit on the very-high-energy gamma-ray flux (Phi) (E > 120 GeV) < 1.7 x 10^-12 erg/cm^-2/s. We obtain a radio upper limit of ~15 microJy/beam for persistent emission at the locations of both FRBs 20190711A and 20171019A, but detect diffuse radio emission with a peak brightness of ~53 microJy/beam associated with FRB 20190714A at z = 0.2365. This represents the first detection of the radio continuum emission potentially associated with the host (galaxy) of FRB 20190714A, and is only the third known FRB to have such an association. Given the possible association of a faint persistent source, FRB 20190714A may potentially be a repeating FRB whose age lies between that of FRB 20121102A and FRB 20180916A. A parallel search for repeat bursts from these FRBs revealed no new detections down to a fluence of 0.08 Jy ms for a 1 ms duration burst. △ Less

Submitted 31 December, 2021; originally announced January 2022.

Comments: 15 pages, 9 figures

A comparative analysis of pulse time-of-arrival creation methods

Authors: J. Wang, G. M. Shaifullah, J. P. W. Verbiest, C. Tiburzi, D. J. Champion, I. Cognard, M. Gaikwad, E. Graikou, L. Guillemot, H. Hu, R. Karuppusamy, Michael J. Keith, Michael Kramer, Y. Liu, A. G. Lyne, M. B. Mickaliger, B. W. Stappers, G. Theureau

Abstract: Extracting precise pulse times of arrival (TOAs) and their uncertainties is the first and most fundamental step in high-precision pulsar timing. In the classical method, TOAs are derived from total intensity pulse profiles of pulsars via cross-correlation with an idealised `1D' template of that profile. While a number of results have been presented in the literature relying on the ever increasing… ▽ More Extracting precise pulse times of arrival (TOAs) and their uncertainties is the first and most fundamental step in high-precision pulsar timing. In the classical method, TOAs are derived from total intensity pulse profiles of pulsars via cross-correlation with an idealised `1D' template of that profile. While a number of results have been presented in the literature relying on the ever increasing sensitivity of such pulsar timing experiments, there is no consensus on the most reliable methods for TOA creation and, more importantly, the associated TOA uncertainties for each scheme. In this article, we present a comprehensive comparison of TOA determination practices, focusing on the creation of timing templates, TOA determination methods and the most useful TOA bandwidth. The aim is both to present a possible approach towards TOA optimisation as well as the (partial) identification of an optimal TOA-creation scheme and the demonstration of optimisation differences between pulsars and data sets. We compare the values of data-derived template profiles as compared to analytic profiles and evaluate the three most commonly used template-matching methods. Finally, we study the relation between timing precision and TOA bandwidth to identify any potential breaks in that relationship. As a practical demonstration, we apply our selected methods to European Pulsar Timing Array data on the three test pulsars PSRs\ J0218+4232, J1713+0747 and J2145$-$0750. △ Less

Submitted 26 November, 2021; originally announced November 2021.

Journal ref: A&A 658, A181 (2022)

The Jodrell Bank Glitch Catalogue: 106 new rotational glitches in 70 pulsars

Authors: Avishek Basu, Benjamin Shaw, Danai Antonopoulou, Michael J. Keith, Andrew G. Lyne, Mitchell B. Mickaliger, Benjamin W. Stappers, Patrick Weltevrede, Christine A. Jordan

Abstract: Pulsar glitches are rapid spin-up events that occur in the rotation of neutron stars, providing a valuable probe into the physics of the interiors of these objects. Long-term monitoring of a large number of pulsars facilitates the detection of glitches and the robust measurements of their parameters. The Jodrell Bank pulsar timing programme regularly monitors more than 800 radio pulsars and has ac… ▽ More Pulsar glitches are rapid spin-up events that occur in the rotation of neutron stars, providing a valuable probe into the physics of the interiors of these objects. Long-term monitoring of a large number of pulsars facilitates the detection of glitches and the robust measurements of their parameters. The Jodrell Bank pulsar timing programme regularly monitors more than 800 radio pulsars and has accrued, in some cases, over 50 years of timing history on individual objects. In this paper we present 106 new glitches in 70 radio pulsars as observed up to the end of 2018. For 70% of these pulsars, the event we report is its only known glitch. For each new glitch we provide measurements of its epoch, amplitude and any detected changes to the spin-down rate of the star. Combining these new glitches with those listed in the Jodrell Bank glitch catalogue we analyse a total sample of 543 glitches in 178 pulsars. We model the distribution of glitch amplitudes and spin-down rate changes using a mixture of two Gaussian components. We corroborate the known dependence of glitch rate and activity on pulsar spin-down rates and characteristic ages, and show that younger pulsars tend to exhibit larger glitches. Pulsars whose spin-down rates between $10^{-14}$ Hz s$^{-1}$ and $10^{-10.5}$ Hz s$^{-1}$ show a mean reversal of 1.8% of their spin-down as a consequence of glitches. Our results are qualitatively consistent with the superfluid vortex unpinning models of pulsar glitches. △ Less

Submitted 12 November, 2021; originally announced November 2021.

Comments: 14 pages, 17 figures, 3 tables. Accepted for publication in MNRAS

Noise analysis in the European Pulsar Timing Array data release 2 and its implications on the gravitational-wave background search

Authors: A. Chalumeau, S. Babak, A. Petiteau, S. Chen, A. Samajdar, R. N. Caballero, G. Theureau, L. Guillemot, G. Desvignes, A. Parthasarathy, K. Liu, G. Shaifullah, H. Hu, E. van der Wateren, J. Antoniadis, A. -S. Bak Nielsen, C. G. Bassa, A. Berthereau, M. Burgay, D. J. Champion, I. Cognard, M. Falxa, R. D. Ferdman, P. C. C. Freire, J. R. Gair , et al. (27 additional authors not shown)

Abstract: The European Pulsar Timing Array (EPTA) collaboration has recently released an extended data set for six pulsars (DR2) and reported evidence for a common red noise signal. Here we present a noise analysis for each of the six pulsars. We consider several types of noise: (i) radio frequency independent, "achromatic", and time-correlated red noise; (ii) variations of dispersion measure and scattering… ▽ More The European Pulsar Timing Array (EPTA) collaboration has recently released an extended data set for six pulsars (DR2) and reported evidence for a common red noise signal. Here we present a noise analysis for each of the six pulsars. We consider several types of noise: (i) radio frequency independent, "achromatic", and time-correlated red noise; (ii) variations of dispersion measure and scattering; (iii) system and band noise; and (iv) deterministic signals (other than gravitational waves) that could be present in the PTA data. We perform Bayesian model selection to find the optimal combination of noise components for each pulsar. Using these custom models we revisit the presence of the common uncorrelated red noise signal previously reported in the EPTA DR2 and show that the data still supports it with a high statistical significance. Next, we confirm that there is no preference for or against the Hellings-Downs spatial correlations expected for the stochastic gravitational-wave background. The main conclusion of the EPTA DR2 paper remains unchanged despite a very significant change in the noise model of each pulsar. However, modelling the noise is essential for the robust detection of gravitational waves and its impact could be significant when analysing the next EPTA data release, which will include a larger number of pulsars and more precise measurements. △ Less

Submitted 9 November, 2021; originally announced November 2021.

Comments: 21 pages, 11 figures, 7 tables, 1 appendix figure and 1 appendix table, accepted for publication to MNRAS

Common-red-signal analysis with 24-yr high-precision timing of the European Pulsar Timing Array: Inferences in the stochastic gravitational-wave background search

Authors: S. Chen, R. N. Caballero, Y. J. Guo, A. Chalumeau, K. Liu, G. Shaifullah, K. J. Lee, S. Babak, G. Desvignes, A. Parthasarathy, H. Hu, E. van der Wateren, J. Antoniadis, A. -S. Bak Nielsen, C. G. Bassa, A. Berthereau, M. Burgay, D. J. Champion, I. Cognard, M. Falxa, R. D. Ferdman, P. C. C. Freire, J. R. Gair, E. Graikou, L. Guillemot , et al. (27 additional authors not shown)

Abstract: We present results from the search for a stochastic gravitational-wave background (GWB) as predicted by the theory of General Relativity using six radio millisecond pulsars from the Data Release 2 (DR2) of the European Pulsar Timing Array (EPTA) covering a timespan up to 24 years. A GWB manifests itself as a long-term low-frequency stochastic signal common to all pulsars, a common red signal (CRS)… ▽ More We present results from the search for a stochastic gravitational-wave background (GWB) as predicted by the theory of General Relativity using six radio millisecond pulsars from the Data Release 2 (DR2) of the European Pulsar Timing Array (EPTA) covering a timespan up to 24 years. A GWB manifests itself as a long-term low-frequency stochastic signal common to all pulsars, a common red signal (CRS), with the characteristic Hellings-Downs (HD) spatial correlation. Our analysis is performed with two independent pipelines, \eprise{} and \tn{}+\ftwo{}, which produce consistent results. A search for a CRS with simultaneous estimation of its spatial correlations yields spectral properties compatible with theoretical GWB predictions, but does not result in the required measurement of the HD correlation, as required for GWB detection. Further Bayesian model comparison between different types of CRSs, including a GWB, finds the most favoured model to be the common uncorrelated red noise described by a power-law with $A = 5.13_{-2.73}^{+4.20} \times 10^{-15}$ and $γ= 3.78_{-0.59}^{+0.69}$ (95\% credible regions). Fixing the spectral index to $γ=13/3$ as expected from the GWB by circular, inspiralling supermassive black-hole binaries results in an amplitude of $A =2.95_{-0.72}^{+0.89} \times 10^{-15}$. We implement three different models, BAYESEPHEM, LINIMOSS and EPHEMGP, to address possible Solar-system ephemeris (SSE) systematics and conclude that our results may only marginally depend on these effects. This work builds on the methods and models from the studies on the EPTA DR1. We show that under the same analysis framework the results remain consistent after the data set extension. △ Less

Submitted 25 October, 2021; originally announced October 2021.

Comments: 25 pages, 10 figures, 7 tables, 2 appendix tables and 1 appendix figure

Timing observations of three Galactic millisecond pulsars

Authors: D. R. Lorimer, A. M. Kawash, P. C. C. Freire, D. A. Smith, M. Kerr, M. A. McLaughlin, M. B. Mickaliger, R. Spiewak, M. Bailes, E. Barr, M. Burgay, A. D. Cameron, F. Camilo, S. Johnston, F. Jankowski, E. F. Keane, M. Keith, M. Kramer, A. Possenti

Abstract: We report observed and derived timing parameters for three millisecond pulsars (MSPs) from observations collected with the Parkes 64-m telescope, Murriyang. The pulsars were found during re-processing of archival survey data by Mickaliger et al. One of the new pulsars (PSR J1546-5925) has a spin period $P=7.8$ ms and is isolated. The other two (PSR J0921-5202 with $P=9.7$ ms and PSR J1146-6610 wit… ▽ More We report observed and derived timing parameters for three millisecond pulsars (MSPs) from observations collected with the Parkes 64-m telescope, Murriyang. The pulsars were found during re-processing of archival survey data by Mickaliger et al. One of the new pulsars (PSR J1546-5925) has a spin period $P=7.8$ ms and is isolated. The other two (PSR J0921-5202 with $P=9.7$ ms and PSR J1146-6610 with $P=3.7$ ms) are in binary systems around low-mass ($>0.2 M_{\odot}$) companions. Their respective orbital periods are $38$.2 d and $62.8$ d. While PSR J0921-5202 has a low orbital eccentricity $e=1.3 \times 10^{-5}$, in keeping with many other Galactic MSPs, PSR J1146-6610 has a significantly larger eccentricity, $e = 7.4 \times 10^{-3}$. This makes it a likely member of a group of eccentric MSP-He white dwarf binary systems in the Galactic disk whose formation is poorly understood. Two of the pulsars are co-located with previously unidentified point sources discovered with the Fermi satellite's Large Area Telescope, but no $γ$-ray pulsations have been detected, likely due to their low spin-down powers. We also show that, particularly in terms of orbital diversity, the current sample of MSPs is far from complete and is subject to a number of selection biases. △ Less

Submitted 27 August, 2021; v1 submitted 9 August, 2021; originally announced August 2021.

Comments: 8 pages, 2 figures, accepted for publication in MNRAS

Study of 72 pulsars discovered in the PALFA survey: Timing analysis, glitch activity, emission variability, and a pulsar in an eccentric binary

Authors: E. Parent, H. Sewalls, P. C. C. Freire, T. Matheny, A. G. Lyne, B. B. P. Perera, F. Cardoso, M. A. McLaughlin, B. Allen, A. Brazier, F. Camilo, S. Chatterjee, J. M. Cordes, F. Crawford, J. S. Deneva, F. A. Dong, R. D. Ferdman, E. Fonseca, J. W. T. Hessels, V. M. Kaspi, B. Knispel, J. van Leeuwen, R. S. Lynch, B. M. Meyers, J. W. McKee , et al. (9 additional authors not shown)

Abstract: We present new discoveries and results from long-term timing of 72 pulsars discovered in the Arecibo PALFA survey, including precise determination of astrometric and spin parameters, and flux density and scatter broadening measurements at 1.4 GHz. Notable discoveries include two young pulsars (characteristic ages $\sim$30 kyr) with no apparent supernova remnant associations, three mode changing, 1… ▽ More We present new discoveries and results from long-term timing of 72 pulsars discovered in the Arecibo PALFA survey, including precise determination of astrometric and spin parameters, and flux density and scatter broadening measurements at 1.4 GHz. Notable discoveries include two young pulsars (characteristic ages $\sim$30 kyr) with no apparent supernova remnant associations, three mode changing, 12 nulling and two intermittent pulsars. We detected eight glitches in five pulsars. Among them is PSR J1939+2609, an apparently old pulsar (characteristic age $\sim$1 Gy), and PSR J1954+2529, which likely belongs to a newly-emerging class of binary pulsars. The latter is the only pulsar among the 72 that is clearly not isolated: a non-recycled neutron star with a 931-ms spin period in an eccentric ($e\,=\,0.114$) wide ($P_b\,=\,82.7\,$d) orbit with a companion of undetermined nature having a minimum mass of $\sim0.6\,M_{\odot}$. Since operations at Arecibo ceased in 2020 August, we give a final tally of PALFA sky coverage, and compare its 207 pulsar discoveries to the known population. On average, they are 50% more distant than other Galactic plane radio pulsars; PALFA millisecond pulsars (MSP) have twice the dispersion measure per unit spin period than the known population of MSP in the Plane. The four intermittent pulsars discovered by PALFA more than double the population of such objects, which should help to improve our understanding of pulsar magnetosphere physics. The statistics for these, RRATS, and nulling pulsars suggest that there are many more of these objects in the Galaxy than was previously thought. △ Less

Submitted 17 November, 2021; v1 submitted 4 August, 2021; originally announced August 2021.

Comments: 39 pages, 18 figures, 9 tables. Accepted for publication in ApJ

The slow rise and recovery of the 2019 Crab pulsar glitch

Authors: B. Shaw, M. J. Keith, A. G. Lyne, M. B. Mickaliger, B. W. Stappers, J. D. Turner, P. Weltevrede

Abstract: We present updated measurements of the Crab pulsar glitch of 2019 July 23 using a dataset of pulse arrival times spanning $\sim$5 months. On MJD 58687, the pulsar underwent its seventh largest glitch observed to date, characterised by an instantaneous spin-up of $\sim$1 $μ$Hz. Following the glitch the pulsar's rotation frequency relaxed exponentially towards pre-glitch values over a timescale of a… ▽ More We present updated measurements of the Crab pulsar glitch of 2019 July 23 using a dataset of pulse arrival times spanning $\sim$5 months. On MJD 58687, the pulsar underwent its seventh largest glitch observed to date, characterised by an instantaneous spin-up of $\sim$1 $μ$Hz. Following the glitch the pulsar's rotation frequency relaxed exponentially towards pre-glitch values over a timescale of approximately one week, resulting in a permanent frequency increment of $\sim$0.5 $μ$Hz. Due to our semi-continuous monitoring of the Crab pulsar, we were able to partially resolve a fraction of the total spin-up. This delayed spin-up occurred exponentially over a timescale of $\sim$18 hours. This is the sixth Crab pulsar glitch for which part of the initial rise was resolved in time and this phenomenon has not been observed in any other glitching pulsars, offering a unique opportunity to study the microphysical processes governing interactions between the neutron star interior and the crust. △ Less

Submitted 26 March, 2021; v1 submitted 24 March, 2021; originally announced March 2021.

Comments: 5 pages, 3 figures

Multi-frequency observations of SGR J1935+2154

Authors: M. Bailes, C. G. Bassa, G. Bernardi, S. Buchner, M. Burgay, M. Caleb, A. J. Cooper, G. Desvignes, P. J. Groot, I. Heywood, F. Jankowski, R. Karuppusamy, M. Kramer, M. Malenta, G. Naldi, M. Pilia, G. Pupillo, K. M. Rajwade, L. Spitler, M. Surnis, B. W. Stappers, A. Addis, S. Bloemen, M. C. Bezuidenhout, G. Bianchi , et al. (32 additional authors not shown)

Abstract: Magnetars are a promising candidate for the origin of Fast Radio Bursts (FRBs). The detection of an extremely luminous radio burst from the Galactic magnetar SGR J1935+2154 on 2020 April 28 added credence to this hypothesis. We report on simultaneous and non-simultaneous observing campaigns using the Arecibo, Effelsberg, LOFAR, MeerKAT, MK2 and Northern Cross radio telescopes and the MeerLICHT opt… ▽ More Magnetars are a promising candidate for the origin of Fast Radio Bursts (FRBs). The detection of an extremely luminous radio burst from the Galactic magnetar SGR J1935+2154 on 2020 April 28 added credence to this hypothesis. We report on simultaneous and non-simultaneous observing campaigns using the Arecibo, Effelsberg, LOFAR, MeerKAT, MK2 and Northern Cross radio telescopes and the MeerLICHT optical telescope in the days and months after the April 28 event. We did not detect any significant single radio pulses down to fluence limits between 25 mJy ms and 18 Jy ms. Some observing epochs overlapped with times when X-ray bursts were detected. Radio images made on four days using the MeerKAT telescope revealed no point-like persistent or transient emission at the location of the magnetar. No transient or persistent optical emission was detected over seven days. Using the multi-colour MeerLICHT images combined with relations between DM, NH and reddening we constrain the distance to SGR J1935+2154, to be between 1.5 and 6.5 kpc. The upper limit is consistent with some other distance indicators and suggests that the April 28 burst is closer to two orders of magnitude less energetic than the least energetic FRBs. The lack of single-pulse radio detections shows that the single pulses detected over a range of fluences are either rare, or highly clustered, or both. It may also indicate that the magnetar lies somewhere between being radio-quiet and radio-loud in terms of its ability to produce radio emission efficiently. △ Less

Submitted 10 March, 2021; originally announced March 2021.

Comments: Accepted for publication in MNRAS. Corresponding author B. W. Stappers

Measuring Interstellar Delays of PSR J0613-0200 over 7 years, using the Large European Array for Pulsars

Authors: R. A. Main, S. A. Sanidas, J. Antoniadis, C. Bassa, S. Chen, I. Cognard, M. Gaikwad, H. Hu, G. H. Janssen, R. Karuppusamy, M. Kramer, K. J. Lee, K. Liu, G. Mall, J. W. McKee, M. B. Mickaliger, D. Perrodin, B. W. Stappers, C. Tiburzi, O. Wucknitz, L. Wang, W. W. Zhu

Abstract: Using data from the Large European Array for Pulsars (LEAP), and the Effelsberg telescope, we study the scintillation parameters of the millisecond pulsar J0613-0200 over a 7 year timespan. The "secondary spectrum" -- the 2D power spectrum of scintillation -- presents the scattered power as a function of time delay, and contains the relative velocities of the pulsar, observer, and scattering mater… ▽ More Using data from the Large European Array for Pulsars (LEAP), and the Effelsberg telescope, we study the scintillation parameters of the millisecond pulsar J0613-0200 over a 7 year timespan. The "secondary spectrum" -- the 2D power spectrum of scintillation -- presents the scattered power as a function of time delay, and contains the relative velocities of the pulsar, observer, and scattering material. We detect a persistent parabolic scintillation arc, suggesting scattering is dominated by a thin, anisotropic region. The scattering is poorly described by a simple exponential tail, with excess power at high delays; we measure significant, detectable scattered power at times out to $\sim 5 μs$, and measure the bulk scattering delay to be between 50 to 200\,ns with particularly strong scattering throughout 2013. These delays are too small to detect a change of the pulse profile shape, yet they would change the times-of-arrival as measured through pulsar timing. The arc curvature varies annually, and is well fit by a one-dimensional scattering screen $\sim 40\%$ of the way towards the pulsar, with a changing orientation during the increased scattering in 2013. Effects of uncorrected scattering will introduce time delays correlated over time in individual pulsars, and may need to be considered in gravitational wave analyses. Pulsar timing programs would benefit from simultaneously recording in a way that scintillation can be resolved, in order to monitor the variable time delays caused by multipath propagation. △ Less

Submitted 3 November, 2020; v1 submitted 22 September, 2020; originally announced September 2020.

Comments: 12 pages, 10 figures. Published in MNRAS

Spectro-temporal analysis of a sample of bursts from FRB 121102

Authors: Kaustubh Rajwade, Mitchell Mickaliger, Benjamin Stappers, Manisha Caleb, Rene Breton, Aris Karastergiou, Evan Keane

Abstract: FRB~121102 was the first Fast Radio Burst (FRB) that was shown to repeat. Since its discovery in 2012, more than two hundred bursts have been detected from the source. These bursts exhibit a diverse range of spectral and temporal characteristics and many questions about their origin and form remain unanswered. Here, we present a sample of radio bursts from FRB 121102 detected using the Lovell tele… ▽ More FRB~121102 was the first Fast Radio Burst (FRB) that was shown to repeat. Since its discovery in 2012, more than two hundred bursts have been detected from the source. These bursts exhibit a diverse range of spectral and temporal characteristics and many questions about their origin and form remain unanswered. Here, we present a sample of radio bursts from FRB 121102 detected using the Lovell telescope at Jodrell Bank Observatory. We show four examples of bursts that show peculiar spectro-temporal characteristics and compare them with properties of bursts of FRB~121102 detected at other observatories. We report on a precursor burst that is separated by just 17~ms from the main burst, the shortest reported separation between two individual bursts to date. We also provide access to data for all the detections of FRB~121102 in this campaign. △ Less

Submitted 3 September, 2020; originally announced September 2020.

Comments: 3 pages, 1 Figure, to be published in RNAAS

High-cadence observations and variable spin behaviour of magnetar Swift J1818.0-1607 after its outburst

Authors: David Champion, Ismael Cognard, Marilyn Cruces, Gregory Desvignes, Fabian Jankowski, Ramesh Karuppusamy, Michael J. Keith, Chryssa Kouveliotou, Michael Kramer, Kuo Liu, Andrew G. Lyne, Mitchell B. Mickaliger, Brendan O'Connor, Aditya Parthasarathy, Nataliya Porayko, Kaustubh Rajwade, Ben W. Stappers, Pablo Torne, Alexander J. van der Horst, Patrick Weltevrede

Abstract: We report on multi-frequency radio observations of the new magnetar Swift J1818.0-1607, following it for more than one month with high cadence. The observations commenced less than 35 hours after its registered first outburst. We obtained timing, polarisation and spectral information. Swift J1818.0-1607 has an unusually steep spectrum for a radio emitting magnetar and also has a relatively narrow… ▽ More We report on multi-frequency radio observations of the new magnetar Swift J1818.0-1607, following it for more than one month with high cadence. The observations commenced less than 35 hours after its registered first outburst. We obtained timing, polarisation and spectral information. Swift J1818.0-1607 has an unusually steep spectrum for a radio emitting magnetar and also has a relatively narrow and simple pulse profile. The position angle swing of the polarisation is flat over the pulse profile, possibly suggesting that our line-of-sight grazes the edge of the emission beam. This may also explain the steep spectrum. The spin evolution shows large variation in the spin-down rate, associated with four distinct timing events over the course of our observations. Those events may be related to the appearance and disappearance of a second pulse component. The first timing event coincides with our actual observations, while we did not detect significant changes in the emission properties which could reveal further magnetospheric changes. Characteristic ages inferred from the timing measurements over the course of months vary by nearly an order of magnitude. A longer-term spin-down measurement over approximately 100 days suggests an characteristic age of about 500 years, larger than previously reported. Though Swift J1818.0-1607 could still be one of the youngest neutron stars (and magnetars) detected so far, we caution using the characteristic age as a true-age indicator given the caveats behind its calculation. △ Less

Submitted 8 September, 2020; originally announced September 2020.

Comments: 14 pages, 10 figures, accepted by MNRAS

Discovery of a Gamma-ray Black Widow Pulsar by GPU-accelerated Einstein@Home

Authors: L. Nieder, C. J. Clark, D. Kandel, R. W. Romani, C. G. Bassa, B. Allen, A. Ashok, I. Cognard, H. Fehrmann, P. Freire, R. Karuppusamy, M. Kramer, D. Li, B. Machenschalk, Z. Pan, M. A. Papa, S. M. Ransom, P. S. Ray, J. Roy, P. Wang, J. Wu, C. Aulbert, E. D. Barr, B. Beheshtipour, O. Behnke , et al. (17 additional authors not shown)

Abstract: We report the discovery of 1.97 ms period gamma-ray pulsations from the 75 minute orbital-period binary pulsar now named PSR J1653-0158. The associated Fermi Large Area Telescope gamma-ray source 4FGL J1653.6-0158 has long been expected to harbor a binary millisecond pulsar. Despite the pulsar-like gamma-ray spectrum and candidate optical/X-ray associations -- whose periodic brightness modulations… ▽ More We report the discovery of 1.97 ms period gamma-ray pulsations from the 75 minute orbital-period binary pulsar now named PSR J1653-0158. The associated Fermi Large Area Telescope gamma-ray source 4FGL J1653.6-0158 has long been expected to harbor a binary millisecond pulsar. Despite the pulsar-like gamma-ray spectrum and candidate optical/X-ray associations -- whose periodic brightness modulations suggested an orbit -- no radio pulsations had been found in many searches. The pulsar was discovered by directly searching the gamma-ray data using the GPU-accelerated Einstein@Home distributed volunteer computing system. The multi-dimensional parameter space was bounded by positional and orbital constraints obtained from the optical counterpart. More sensitive analyses of archival and new radio data using knowledge of the pulsar timing solution yield very stringent upper limits on radio emission. Any radio emission is thus either exceptionally weak, or eclipsed for a large fraction of the time. The pulsar has one of the three lowest inferred surface magnetic-field strengths of any known pulsar with $B_{\rm surf} \approx 4 \times 10^{7}\,$G. The resulting mass function, combined with models of the companion star's optical light curve and spectra, suggests a pulsar mass $\gtrsim 2\,M_{\odot}$. The companion is light-weight with mass $\sim 0.01\,M_{\odot}$, and the orbital period is the shortest known for any rotation-powered binary pulsar. This discovery demonstrates the Fermi Large Area Telescope's potential to discover extreme pulsars that would otherwise remain undetected. △ Less

Submitted 22 October, 2020; v1 submitted 3 September, 2020; originally announced September 2020.

Comments: 12 pages, 3 figures, published in ApJL

Possible periodic activity in the repeating FRB 121102

Authors: K. M. Rajwade, M. B. Mickaliger, B. W. Stappers, V. Morello, D. Agarwal, C. G. Bassa, R. P. Breton, M. Caleb, A. Karastergiou, E. F. Keane, D. R. Lorimer

Abstract: The discovery that at least some Fast Radio Bursts (FRBs) repeat has ruled out cataclysmic events as the progenitors of these particular bursts. FRB~121102 is the most well-studied repeating FRB but despite extensive monitoring of the source, no underlying pattern in the repetition has previously been identified. Here, we present the results from a radio monitoring campaign of FRB~121102 using the… ▽ More The discovery that at least some Fast Radio Bursts (FRBs) repeat has ruled out cataclysmic events as the progenitors of these particular bursts. FRB~121102 is the most well-studied repeating FRB but despite extensive monitoring of the source, no underlying pattern in the repetition has previously been identified. Here, we present the results from a radio monitoring campaign of FRB~121102 using the 76-m Lovell telescope. Using the pulses detected in the Lovell data along with pulses from the literature, we report a detection of periodic behaviour of the source over the span of five years of data. We predict that the source is currently `off' and that it should turn `on' for the approximate MJD range $59002-59089$ (2020-06-02 to 2020-08-28). This result, along with the recent detection of periodicity from another repeating FRB, highlights the need for long-term monitoring of repeating FRBs at a high cadence. Using simulations, we show that one needs at least 100 hours of telescope time to follow-up repeating FRBs at a cadence of 0.5--3 days to detect periodicities in the range of 10--150 days. If the period is real, it shows that repeating FRBs can have a large range in their activity periods that might be difficult to reconcile with neutron star precession models. △ Less

Submitted 29 April, 2020; v1 submitted 7 March, 2020; originally announced March 2020.

Comments: 9 pages, 6 Figures, 3 Tables. More detections added to the periodicity search. Updated period and activity period. Complete tables (Table 1 and Table 3) have been submitted as online supplementary materials

Limits on Absorption from a 332-MHz survey for Fast Radio Bursts

Authors: K. M. Rajwade, M. B. Mickaliger, B. W. Stappers, C. G. Bassa, R. P. Breton, A. Karastergiou, E. F. Keane

Abstract: Fast Radio Bursts (FRBs) are bright, extragalactic radio pulses whose origins are still unknown. Until recently, most FRBs have been detected at frequencies greater than 1 GHz with a few exceptions at 800 MHz. The recent discoveries of FRBs at 400 MHz from the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope has opened up possibilities for new insights about the progenitors while m… ▽ More Fast Radio Bursts (FRBs) are bright, extragalactic radio pulses whose origins are still unknown. Until recently, most FRBs have been detected at frequencies greater than 1 GHz with a few exceptions at 800 MHz. The recent discoveries of FRBs at 400 MHz from the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope has opened up possibilities for new insights about the progenitors while many other low frequency surveys in the past have failed to find any FRBs. Here, we present results from a FRB survey recently conducted at the Jodrell Bank Observatory at 332 MHz with the 76-m Lovell telescope for a total of 58 days. We did not detect any FRBs in the survey and report a 90$\%$ upper limit of 5500 FRBs per day per sky for a Euclidean Universe above a fluence threshold of 46 Jy ms. We discuss the possibility of absorption as the main cause of non-detections in low frequency (< 800 MHz) searches and invoke different absorption models to explain the same. We find that Induced Compton Scattering alone cannot account for absorption of radio emission and that our simulations favour a combination of Induced Compton Scattering and Free-Free Absorption to explain the non-detections. For a free-free absorption scenario, our constraints on the electron density are consistent with those expected in the post-shock region of the ionized ejecta in Super-Luminous SuperNovae (SLSNe). △ Less

Submitted 25 February, 2020; originally announced February 2020.

Comments: 12 pages, 9 Figures, 2 Tables, Second revision submitted to MNRAS

A fast radio burst in the direction of the Virgo cluster

Authors: Devansh Agarwal, Duncan R. Lorimer, Anastasia Fialkov, Keith W. Bannister, Ryan M. Shannon, Wael Farah, Shivani Bhandari, Jean-Pierre Macquart, Chris Flynn, Giuliano Pignata, Nicolas Tejos, Benjamin Gregg, Stefan Osłowski, Kaustubh Rajwade, Mitchell B. Mickaliger, Benjamin W. Stappers, Di Li, Weiwei Zhu, Lei Qian, Youling Yue, Pei Wang, Abraham Loeb

Abstract: The rate of fast radio bursts (FRBs) in the direction of nearby galaxy clusters is expected to be higher than the mean cosmological rate if intrinsically faint FRBs are numerous. In this paper, we describe a targeted search for faint FRBs near the core of the Virgo cluster using the Australian Square Kilometer Array Pathfinder telescope. During 300 hr of observations, we discovered one burst, FRB… ▽ More The rate of fast radio bursts (FRBs) in the direction of nearby galaxy clusters is expected to be higher than the mean cosmological rate if intrinsically faint FRBs are numerous. In this paper, we describe a targeted search for faint FRBs near the core of the Virgo cluster using the Australian Square Kilometer Array Pathfinder telescope. During 300 hr of observations, we discovered one burst, FRB 180417, with dispersion measure DM $=474.8$ cm$^{-3}$pc. The FRB was promptly followed up by several radio telescopes for 27 h, but no repeat bursts were detected. An optical follow-up of FRB 180417 using the PROMPT5 telescope revealed no new sources down to an $R$-band magnitude of 20.1. We argue that FRB 180417 is likely behind the Virgo cluster as the Galactic and intracluster DM contribution are small compared to the DM of the FRB, and there are no galaxies in the line of sight. The non-detection of FRBs from Virgo constrains the faint-end slope, $α<1.52$ (at 68\% confidence limit), and the minimum luminosity, $L_{\rm min}\gtrsim 2\times 10^{40}$ erg s$^{-1}$ (at 68\% confidence limit), of the FRB luminosity function assuming cosmic FRB rate of $10^4$ FRBs sky$^{-1}$ day$^{-1}$ with flux above 1 Jy located out to redshift of 1. Further FRB surveys of galaxy clusters with high-sensitivity instruments will tighten the constraints on the faint end of the luminosity function and, thus, are strongly encouraged. △ Less

Submitted 12 September, 2019; originally announced September 2019.

Comments: 9 pages, 6 figures, Accepted for publication in MNRAS

The dynamics of Galactic centre pulsars: constraining pulsar distances and intrinsic spin-down

Authors: B. B. P. Perera, E. D. Barr, M. B. Mickaliger, A. G. Lyne, D. R. Lorimer, B. W. Stappers, R. P. Eatough, M. Kramer, C. Ng, R. Spiewak, M. Bailes, D. J. Champion, V. Morello, A. Possenti

Abstract: Through high-precision radio timing observations, we show that five recycled pulsars in the direction of the Galactic Centre (GC) have anomalous spin period time derivative ($\dot P$) measurements -- PSRs J1748$-$3009, J1753$-$2819, J1757$-$2745, and J1804$-$2858 show negative values of $\dot P$ and PSR J1801$-$3210 is found to have an exceptionally small value of $\dot P$. We attribute these obse… ▽ More Through high-precision radio timing observations, we show that five recycled pulsars in the direction of the Galactic Centre (GC) have anomalous spin period time derivative ($\dot P$) measurements -- PSRs J1748$-$3009, J1753$-$2819, J1757$-$2745, and J1804$-$2858 show negative values of $\dot P$ and PSR J1801$-$3210 is found to have an exceptionally small value of $\dot P$. We attribute these observed $\dot P$ measurements to acceleration of these pulsars along their lines-of-sight (LOSs) due to the Galactic gravitational field. Using models of the Galactic mass distribution and pulsar velocities, we constrain the distances to these pulsars, placing them on the far-side of the Galaxy, providing the first accurate distance measurements to pulsars located in this region and allowing us to consider the electron density along these LOSs. We find the new electron density model YMW16 to be more consistent with these observations than the previous model NE2001. The LOS dynamics further constrain the model-dependent intrinsic $\dot P$ values for these pulsars and they are consistent with measurements for other known pulsars. In the future, the independent distance measurements to these and other pulsars near the GC would allow us to constrain the Galactic gravitational potential more accurately. △ Less

Submitted 25 April, 2019; originally announced April 2019.

Comments: 17 pages, 11 figures. Accepted 2019 April 25 for publication in MNRAS

Spin frequency evolution and pulse profile variations of the recently re-activated radio magnetar XTE J1810-197

Authors: L. Levin, A. G. Lyne, G. Desvignes, R. P. Eatough, R. Karuppusamy, M. Kramer, M. Mickaliger, B. W. Stappers, P. Weltevrede

Abstract: After spending almost a decade in a radio-quiet state, the Anomalous X-ray Pulsar XTE J1810-197 turned back on in early December 2018. We have observed this radio magnetar at 1.5 GHz with ~daily cadence since the first detection of radio re-activation on 8 December 2018. In this paper, we report on the current timing properties of XTE J1810-197 and find that the magnitude of the spin frequency der… ▽ More After spending almost a decade in a radio-quiet state, the Anomalous X-ray Pulsar XTE J1810-197 turned back on in early December 2018. We have observed this radio magnetar at 1.5 GHz with ~daily cadence since the first detection of radio re-activation on 8 December 2018. In this paper, we report on the current timing properties of XTE J1810-197 and find that the magnitude of the spin frequency derivative has increased by a factor of 2.6 over our 48-day data set. We compare our results with the spin-down evolution reported during its previous active phase in the radio band. We also present total intensity pulse profiles at five different observing frequencies between 1.5 and 8.4 GHz, collected with the Lovell and the Effelsberg telescopes. The profile evolution in our data set is less erratic than what was reported during the previous active phase, and can be seen varying smoothly between observations. Profiles observed immediately after the outburst show the presence of at least five cycles of a very stable ~50-ms periodicity in the main pulse component that lasts for at least tens of days. This remarkable structure is seen across the full range of observing frequencies. △ Less

Submitted 28 July, 2019; v1 submitted 6 March, 2019; originally announced March 2019.

Comments: 9 pages, 7 figures, updated with additional analysis of the 50-ms oscillation, accepted for publication in MNRAS

A Processing Pipeline for High Volume Pulsar Data Streams

Authors: R. J. Lyon, B. W. Stappers, L. Levin, M. B. Mickaliger, A. Scaife

Abstract: Pulsar data analysis pipelines have historically been comprised of bespoke software systems, supporting the off-line analysis of data. However modern data acquisition systems are making off-line analyses impractical. They often output multiple simultaneous high volume data streams, significantly increasing data capture rates. This leads to the accumulation of large data volumes, which are prohibit… ▽ More Pulsar data analysis pipelines have historically been comprised of bespoke software systems, supporting the off-line analysis of data. However modern data acquisition systems are making off-line analyses impractical. They often output multiple simultaneous high volume data streams, significantly increasing data capture rates. This leads to the accumulation of large data volumes, which are prohibitively expensive to retain. To maintain processing capabilities when off-line analysis becomes infeasible due to cost, requires a shift to on-line data processing. This paper makes four contributions facilitating this shift with respect to the search for radio pulsars: i) it characterises for the modern era, the key components of a pulsar search science (not signal processing) pipeline, ii) it examines the feasibility of implementing on-line pulsar search via existing tools, iii) problems preventing an easy transition to on-line search are identified and explained, and finally iv) it provides the design for a new prototype pipeline capable of overcoming such problems. Realised using Commercial off-the-shelf (COTS) software components, the deployable system is open source, simple, scalable, and cheap to produce. It has the potential to achieve pulsar search design requirements for the Square Kilometre Array (SKA), illustrated via testing under simulated SKA loads. △ Less

Submitted 14 October, 2018; originally announced October 2018.

Comments: pre-print submitted to Astronomy & Computing

A Study of Single Pulses in the Parkes Multibeam Pulsar Survey

Authors: Mitchell B. Mickaliger, Alex E. McEwen, Maura A. McLaughlin, Duncan R. Lorimer

Abstract: We reprocessed the Parkes Multibeam Pulsar Survey, searching for single pulses out to a DM of 5000 pc cm$^{-3}$ with widths of up to one second. We recorded single pulses from 264 known pulsars and 14 Rotating Radio Transients. We produced amplitude distributions for each pulsar which we fit with log-normal distributions, power-law tails, and a power-law function divided by an exponential function… ▽ More We reprocessed the Parkes Multibeam Pulsar Survey, searching for single pulses out to a DM of 5000 pc cm$^{-3}$ with widths of up to one second. We recorded single pulses from 264 known pulsars and 14 Rotating Radio Transients. We produced amplitude distributions for each pulsar which we fit with log-normal distributions, power-law tails, and a power-law function divided by an exponential function, finding that some pulsars show a deviation from a log-normal distribution in the form of an excess of high-energy pulses. We found that a function consisting of a power-law divided by an exponential fit the distributions of most pulsars better than either log-normal or power-law functions. For pulsars that were detected in a periodicity search, we computed the ratio of their single-pulse signal-to-noise ratios to their signal-to-noise ratios from a Fourier transform and looked for correlations between this ratio and physical parameters of the pulsars. The only correlation found is the expected relationship between this ratio and the spin period. Fitting log-normal distributions to the amplitudes of pulses from RRATs showed similar behaviour for most RRATs. Here, however, there seem to be two distinct distributions of pulses, with the lower-energy distribution being consistent with noise. Pulse-energy distributions for two of the RRATS processed were consistent with those found for normal pulsars, suggesting that pulsars and RRATs have a common emission mechanism, but other factors influence the specific emission properties of each source class. △ Less

Submitted 30 June, 2018; originally announced July 2018.

Comments: 11 pages, 6 figures, 3 tables, accepted for publication in MNRAS

The largest glitch observed in the Crab pulsar

Authors: B. Shaw, A. G. Lyne, B. W. Stappers, P. Weltevrede, C. G. Bassa, A. Y. Lien, M. B. Mickaliger, R. P. Breton, C. A. Jordan, M. J. Keith, H. A. Krimm

Abstract: We have observed a large glitch in the Crab pulsar (PSR B0531+21). The glitch occurred around MJD 58064 (2017 November 8) when the pulsar underwent an increase in the rotation rate of $Δν= 1.530 \times 10^{-5}$ Hz, corresponding to a fractional increase of $Δν/ ν= 0.516 \times 10^{-6}$ making this event the largest glitch ever observed in this source. Due to our high-cadence and long-dwell time ob… ▽ More We have observed a large glitch in the Crab pulsar (PSR B0531+21). The glitch occurred around MJD 58064 (2017 November 8) when the pulsar underwent an increase in the rotation rate of $Δν= 1.530 \times 10^{-5}$ Hz, corresponding to a fractional increase of $Δν/ ν= 0.516 \times 10^{-6}$ making this event the largest glitch ever observed in this source. Due to our high-cadence and long-dwell time observations of the Crab pulsar we are able to partially resolve a fraction of the total spin-up of the star. This delayed spin-up occurred over a timescale of $\sim$1.7 days and is similar to the behaviour seen in the 1989 and 1996 large Crab pulsar glitches. The spin-down rate also increased at the glitch epoch by $Δ\dotν / \dotν = 7 \times 10^{-3}$. In addition to being the largest such event observed in the Crab, the glitch occurred after the longest period of glitch inactivity since at least 1984 and we discuss a possible relationship between glitch size and waiting time. No changes to the shape of the pulse profile were observed near the glitch epoch at 610 MHz or 1520 MHz, nor did we identify any changes in the X-ray flux from the pulsar. The long-term recovery from the glitch continues to progress as $\dotν$ slowly rises towards pre-glitch values. In line with other large Crab glitches, we expect there to be a persistent change to $\dotν$. We continue to monitor the long-term recovery with frequent, high quality observations. △ Less

Submitted 14 May, 2018; originally announced May 2018.

Comments: 10 pages, 9 figures

LOFT-e: Localisation Of Fast Transients with e-MERLIN

Authors: C. R. H. Walker, R. P. Breton, P. A. Harrison, A. Holloway, M. J. Keith, M. Kramer, M. Malenta, M. B. Mickaliger, J. Roy, T. W. Scragg, B. W. Stappers

Abstract: The majority of fast radio bursts (FRBs) are poorly localised, hindering their potential scientific yield as galactic, intergalactic, and cosmological probes. LOFT-e, a digital backend for the U.K.'s e-MERLIN seven-telescope interferometer will provide commensal search and real-time detection of FRBs, taking full advantage of its field of view (FoV), sensitivity, and observation time. Upon burst d… ▽ More The majority of fast radio bursts (FRBs) are poorly localised, hindering their potential scientific yield as galactic, intergalactic, and cosmological probes. LOFT-e, a digital backend for the U.K.'s e-MERLIN seven-telescope interferometer will provide commensal search and real-time detection of FRBs, taking full advantage of its field of view (FoV), sensitivity, and observation time. Upon burst detection, LOFT-e will store raw data offline, enabling the sub-arcsecond localisation provided by e-MERLIN and expanding the pool of localised FRBs. The high-time resolution backend will additionally introduce pulsar observing capabilities to e-MERLIN. △ Less

Submitted 4 April, 2018; originally announced April 2018.

Comments: 2 pages, no figures, prepared for Proceedings of IAU Symposium 337: Pulsar Astrophysics - The Next 50 Years at Jodrell Bank Observatory, UK (September 4 - 8 2017)

Pulsar Searches with the SKA

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

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

Submitted 4 December, 2017; originally announced December 2017.

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

The SUrvey for Pulsars and Extragalactic Radio Bursts II: New FRB discoveries and their follow-up

Authors: S. Bhandari, E. F. Keane, E. D. Barr, A. Jameson, E. Petroff, S. Johnston, M. Bailes, N. D. R. Bhat, M. Burgay, S. Burke-Spolaor, M. Caleb, R. P. Eatough, C. Flynn, J. A. Green, F. Jankowski, M. Kramer, V. Venkatraman Krishnan, V. Morello, A. Possenti, B. Stappers, C. Tiburzi, W. van Straten, I. Andreoni, T. Butterley, P. Chandra , et al. (25 additional authors not shown)

Abstract: We report the discovery of four Fast Radio Bursts (FRBs) in the ongoing SUrvey for Pulsars and Extragalactic Radio Bursts (SUPERB) at the Parkes Radio Telescope: FRBs 150610, 151206, 151230 and 160102. Our real-time discoveries have enabled us to conduct extensive, rapid multi-messenger follow-up at 12 major facilities sensitive to radio, optical, X-ray, gamma-ray photons and neutrinos on time sca… ▽ More We report the discovery of four Fast Radio Bursts (FRBs) in the ongoing SUrvey for Pulsars and Extragalactic Radio Bursts (SUPERB) at the Parkes Radio Telescope: FRBs 150610, 151206, 151230 and 160102. Our real-time discoveries have enabled us to conduct extensive, rapid multi-messenger follow-up at 12 major facilities sensitive to radio, optical, X-ray, gamma-ray photons and neutrinos on time scales ranging from an hour to a few months post-burst. No counterparts to the FRBs were found and we provide upper limits on afterglow luminosities. None of the FRBs were seen to repeat. Formal fits to all FRBs show hints of scattering while their intrinsic widths are unresolved in time. FRB 151206 is at low Galactic latitude, FRB 151230 shows a sharp spectral cutoff, and FRB 160102 has the highest dispersion measure (DM = $2596.1\pm0.3$ pc cm$^{-3}$) detected to date. Three of the FRBs have high dispersion measures (DM >$1500$ pc cm$^{-3}$), favouring a scenario where the DM is dominated by contributions from the Intergalactic Medium. The slope of the Parkes FRB source counts distribution with fluences $>2$ Jyms is $α=-2.2^{+0.6}_{-1.2}$ and still consistent with a Euclidean distribution ($α=-3/2$). We also find that the all-sky rate is $1.7^{+1.5}_{-0.9}\times10^3$FRBs/($4π$ sr)/day above $\sim2$ Jyms and there is currently no strong evidence for a latitude-dependent FRB sky-rate. △ Less

Submitted 21 November, 2017; originally announced November 2017.

Comments: 21 pages, 8 figures and accepted for publication in MNRAS

A polarized fast radio burst at low Galactic latitude

Authors: E. Petroff, S. Burke-Spolaor, E. F. Keane, M. A. McLaughlin, R. Miller, I. Andreoni, M. Bailes, E. D. Barr, S. R. Bernard, S. Bhandari, N. D. R. Bhat, M. Burgay, M. Caleb, D. Champion, P. Chandra, J. Cooke, V. S. Dhillon, J. S. Farnes, L. K. Hardy, P. Jaroenjittichai, S. Johnston, M. Kasliwal, M. Kramer, S. P. Littlefair, J. P. Macquart , et al. (13 additional authors not shown)

Abstract: We report on the discovery of a new fast radio burst, FRB 150215, with the Parkes radio telescope on 2015 February 15. The burst was detected in real time with a dispersion measure (DM) of 1105.6$\pm$0.8 pc cm^{-3}, a pulse duration of 2.8^{+1.2}_{-0.5} ms, and a measured peak flux density assuming the burst was at beam center of 0.7^{+0.2}_{-0.1} Jy. The FRB originated at a Galactic longitude and… ▽ More We report on the discovery of a new fast radio burst, FRB 150215, with the Parkes radio telescope on 2015 February 15. The burst was detected in real time with a dispersion measure (DM) of 1105.6$\pm$0.8 pc cm^{-3}, a pulse duration of 2.8^{+1.2}_{-0.5} ms, and a measured peak flux density assuming the burst was at beam center of 0.7^{+0.2}_{-0.1} Jy. The FRB originated at a Galactic longitude and latitude of 24.66^{\circ}, 5.28^{\circ}, 25 degrees away from the Galactic Center. The burst was found to be 43$\pm$5% linearly polarized with a rotation measure (RM) in the range -9 < RM < 12 rad m^{-2} (95% confidence level), consistent with zero. The burst was followed-up with 11 telescopes to search for radio, optical, X-ray, gamma-ray and neutrino emission. Neither transient nor variable emission was found to be associated with the burst and no repeat pulses have been observed in 17.25 hours of observing. The sightline to the burst is close to the Galactic plane and the observed physical properties of FRB 150215 demonstrate the existence of sight lines of anomalously low RM for a given electron column density. The Galactic RM foreground may approach a null value due to magnetic field reversals along the line of sight, a decreased total electron column density from the Milky Way, or some combination of these effects. A lower Galactic DM contribution might explain why this burst was detectable whereas previous searches at low latitude have had lower detection rates than those out of the plane. △ Less

Submitted 8 May, 2017; originally announced May 2017.

Comments: 19 pages, 9 figures, accepted for publication in MNRAS

The repeating Fast Radio Burst FRB 121102: Multi-wavelength observations and additional bursts

Authors: P. Scholz, L. G. Spitler, J. W. T. Hessels, S. Chatterjee, J. M. Cordes, V. M. Kaspi, R. S. Wharton, C. G. Bassa, S. Bogdanov, F. Camilo, F. Crawford, J. Deneva, J. van Leeuwen, R. Lynch, E. C. Madsen, M. A. McLaughlin, M. Mickaliger, E. Parent, C. Patel, S. M. Ransom, A. Seymour, I. H. Stairs, B. W. Stappers, S. P. Tendulkar

Abstract: We report on radio and X-ray observations of the only known repeating Fast Radio Burst (FRB) source, FRB 121102. We have detected six additional radio bursts from this source: five with the Green Bank Telescope at 2 GHz, and one at 1.4 GHz at the Arecibo Observatory for a total of 17 bursts from this source. All have dispersion measures consistent with a single value ($\sim559$ pc cm$^{-3}$) that… ▽ More We report on radio and X-ray observations of the only known repeating Fast Radio Burst (FRB) source, FRB 121102. We have detected six additional radio bursts from this source: five with the Green Bank Telescope at 2 GHz, and one at 1.4 GHz at the Arecibo Observatory for a total of 17 bursts from this source. All have dispersion measures consistent with a single value ($\sim559$ pc cm$^{-3}$) that is three times the predicted maximum Galactic value. The 2-GHz bursts have highly variable spectra like those at 1.4 GHz, indicating that the frequency structure seen across the individual 1.4 and 2-GHz bandpasses is part of a wideband process. X-ray observations of the FRB 121102 field with the Swift and Chandra observatories show at least one possible counterpart; however, the probability of chance superposition is high. A radio imaging observation of the field with the Jansky Very Large Array at 1.6 GHz yields a 5$σ$ upper limit of 0.3 mJy on any point-source continuum emission. This upper limit, combined with archival WISE 22-$μ$m and IPHAS H$α$ surveys, rules out the presence of an intervening Galactic HII region. We update our estimate of the FRB detection rate in the PALFA survey to be 1.1$^{+3.7}_{-1.0} \times 10^4$ FRBs sky$^{-1}$ day$^{-1}$ (95% confidence) for peak flux density at 1.4 GHz above 300 mJy. We find that the intrinsic widths of the 12 FRB 121102 bursts from Arecibo are, on average, significantly longer than the intrinsic widths of the 13 single-component FRBs detected with the Parkes telescope. △ Less

Submitted 18 October, 2016; v1 submitted 29 March, 2016; originally announced March 2016.

Comments: 18 pages, 5 figures. Accepted for publication in ApJ

A Fast Radio Burst Host Galaxy

Authors: E. F. Keane, S. Johnston, S. Bhandari, E. Barr, N. D. R. Bhat, M. Burgay, M. Caleb, C. Flynn, A. Jameson, M. Kramer, E. Petroff, A. Possenti, W. van Straten, M. Bailes, S. Burke-Spolaor, R. P. Eatough, B. W. Stappers, T. Totani, M. Honma, H. Furusawa, T. Hattori, T. Morokuma, Y. Niino, H. Sugai, T. Terai , et al. (16 additional authors not shown)

Abstract: In recent years, millisecond duration radio signals originating from distant galaxies appear to have been discovered in the so-called Fast Radio Bursts. These signals are dispersed according to a precise physical law and this dispersion is a key observable quantity which, in tandem with a redshift measurement, can be used for fundamental physical investigations. While every fast radio burst has a… ▽ More In recent years, millisecond duration radio signals originating from distant galaxies appear to have been discovered in the so-called Fast Radio Bursts. These signals are dispersed according to a precise physical law and this dispersion is a key observable quantity which, in tandem with a redshift measurement, can be used for fundamental physical investigations. While every fast radio burst has a dispersion measurement, none before now have had a redshift measurement, due to the difficulty in pinpointing their celestial coordinates. Here we present the discovery of a fast radio burst and the identification of a fading radio transient lasting $\sim 6$ days after the event, which we use to identify the host galaxy; we measure the galaxy's redshift to be $z=0.492\pm0.008$. The dispersion measure and redshift, in combination, provide a direct measurement of the cosmic density of ionised baryons in the intergalactic medium of $Ω_{\mathrm{IGM}}=4.9 \pm 1.3\%$, in agreement with the expectation from WMAP, and including all of the so-called "missing baryons". The $\sim6$-day transient is largely consistent with a short gamma-ray burst radio afterglow, and its existence and timescale do not support progenitor models such as giant pulses from pulsars, and supernovae. This contrasts with the interpretation of another recently discovered fast radio burst, suggesting there are at least two classes of bursts. △ Less

Submitted 24 February, 2016; originally announced February 2016.

Comments: Published in Nature, 2016 Feb 25

Journal ref: Nature, 530, 453-456 (25 February 2016)

A Giant Sample of Giant Pulses from the Crab Pulsar

Authors: Mitchell B. Mickaliger, Maura A. McLaughlin, Duncan R. Lorimer, Glen I. Langston, Anna V. Bilous, Vlad I. Kondratiev, Maxim Lyutikov, Scott M. Ransom, Nipuni Palliyaguru

Abstract: We observed the Crab pulsar with the 43-m telescope in Green Bank, WV over a timespan of 15 months. In total we obtained 100 hours of data at 1.2 GHz and seven hours at 330 MHz, resulting in a sample of about 95000 giant pulses (GPs). This is the largest sample, to date, of GPs from the Crab pulsar taken with the same telescope and backend and analyzed as one data set. We calculated power-law fits… ▽ More We observed the Crab pulsar with the 43-m telescope in Green Bank, WV over a timespan of 15 months. In total we obtained 100 hours of data at 1.2 GHz and seven hours at 330 MHz, resulting in a sample of about 95000 giant pulses (GPs). This is the largest sample, to date, of GPs from the Crab pulsar taken with the same telescope and backend and analyzed as one data set. We calculated power-law fits to amplitude distributions for main pulse (MP) and interpulse (IP) GPs, resulting in indices in the range of 2.1-3.1 for MP GPs at 1.2 GHz and in the range of 2.5-3.0 and 2.4-3.1 for MP and IP GPs at 330 MHz. We also correlated the GPs at 1.2 GHz with GPs from the Robert C. Byrd Green Bank Telescope (GBT), which were obtained simultaneously at a higher frequency (8.9 GHz) over a span of 26 hours. In total, 7933 GPs from the 43-m telescope at 1.2 GHz and 39900 GPs from the GBT were recorded during these contemporaneous observations. At 1.2 GHz, 236 (3%) MP GPs and 23 (5%) IP GPs were detected at 8.9 GHz, both with zero chance probability. Another 15 (4%) low-frequency IP GPs were detected within one spin period of high-frequency IP GPs, with a chance probability of 9%. This indicates that the emission processes at high and low radio frequencies are related, despite significant pulse profile shape differences. The 43-m GPs were also correlated with Fermi gamma-ray photons to see if increased pair production in the magnetosphere is the mechanism responsible for GP emission. A total of 92022 GPs and 393 gamma-ray photons were used in this correlation analysis. No significant correlations were found between GPs and gamma-ray photons. This indicates that increased pair production in the magnetosphere is likely not the dominant cause of GPs. Possible methods of GP production may be increased coherence of synchrotron emission or changes in beaming direction. △ Less

Submitted 11 May, 2013; v1 submitted 1 October, 2012; originally announced October 2012.

Comments: 33 pages, 10 figures, 6 tables, accepted for publication in ApJ

Journal ref: The Astrophysical Journal, 760:64 (13pp), 2012