Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • Article
  • Published:

A radio transient with unusually slow periodic emission

Abstract

The high-frequency radio sky is bursting with synchrotron transients from massive stellar explosions and accretion events, but the low-frequency radio sky has, so far, been quiet beyond the Galactic pulsar population and the long-term scintillation of active galactic nuclei. The low-frequency band, however, is sensitive to exotic coherent and polarized radio-emission processes, such as electron-cyclotron maser emission from flaring M dwarfs1, stellar magnetospheric plasma interactions with exoplanets2 and a population of steep-spectrum pulsars3, making Galactic-plane searches a prospect for blind-transient discovery. Here we report an analysis of archival low-frequency radio data that reveals a periodic, low-frequency radio transient. We find that the source pulses every 18.18 min, an unusual periodicity that has, to our knowledge, not been observed previously. The emission is highly linearly polarized, bright, persists for 30–60 s on each occurrence and is visible across a broad frequency range. At times, the pulses comprise short-duration (<0.5 s) bursts; at others, a smoother profile is observed. These profiles evolve on timescales of hours. By measuring the dispersion of the radio pulses with respect to frequency, we have localized the source to within our own Galaxy and suggest that it could be an ultra-long-period magnetar.

This is a preview of subscription content, access via your institution

Access options

Rent or buy this article

Prices vary by article type

from$1.95

to$39.95

Prices may be subject to local taxes which are calculated during checkout

Fig. 1: Sixty-four of the 71 detected pulses of GLEAM-X J162759.5-523504.3 aligned to its measured period P and period derivative \(\dot{{\boldsymbol{P}}}\).
Fig. 2: Maximum brightness of each pulse and total pulse fluence as a function of time across the two observed intervals of activity.
Fig. 3: Dynamic spectrum of the observation recorded at 03:59 on 10 January 2018.
Fig. 4: Transient parameter space populated with common known radio transients.

Similar content being viewed by others

Data availability

Data that support this paper are available at the following public repository: https://github.com/nhurleywalker/GLEAM-X_Periodic_Transient. Further data products can be supplied by the authors on reasonable request.

Code availability

Code that supports this paper is available at the following public repository: https://github.com/nhurleywalker/GLEAM-X_Periodic_Transient. Figure 3 was generated using https://github.com/nhurleywalker/Transient_Phase_Space. Further code can be supplied by the authors on reasonable request.

References

  1. Lynch, C. R., Lenc, E., Kaplan, D. L., Murphy, T. & Anderson, G. E. 154 MHz detection of faint, polarized flares from UV Ceti. Astrophys. J. Lett. 836, L30 (2017).

    Article  ADS  Google Scholar 

  2. Vedantham, H. K. et al. Coherent radio emission from a quiescent red dwarf indicative of star–planet interaction. Nat. Astron. 4, 577–583 (2020).

    Article  ADS  Google Scholar 

  3. Swainston, N. A. et al. Discovery of a steep-spectrum low-luminosity pulsar with the Murchison Widefield Array. Astrophys. J. Lett. 911, L26 (2021).

    Article  ADS  Google Scholar 

  4. Yao, J. M., Manchester, R. N. & Wang, N. A new electron-density model for estimation of pulsar and FRB distances. Astrophys. J. 835, 29 (2017).

    Article  ADS  Google Scholar 

  5. Hutschenreuter, S. & Enßlin, T. A. The Galactic Faraday depth sky revisited. Astron. Astrophys. 633, A150 (2020).

    Article  ADS  CAS  Google Scholar 

  6. Vedantham, H. K. et al. Symmetric achromatic variability in active galaxies: a powerful new gravitational lensing probe? Astrophys. J. 845, 89 (2017).

    Article  ADS  Google Scholar 

  7. Bannister, K. W. et al. Real-time detection of an extreme scattering event: constraints on Galactic plasma lenses. Science 351, 354–356 (2016).

    Article  ADS  CAS  PubMed  Google Scholar 

  8. Benz, A. O. & Güdel, M. Physical processes in magnetically driven flares on the Sun, stars, and young stellar objects. Annu. Rev. Astron. Astrophys. 48, 241–287 (2010).

    Article  ADS  Google Scholar 

  9. Zarka, P., Treumann, R. A., Ryabov, B. P. & Ryabov, V. B. Magnetically-driven planetary radio emissions and application to extrasolar planets. Astrophys. Space Sci. 277, 293–300 (2001).

    Article  ADS  Google Scholar 

  10. Marsh, T. R. et al. A radio-pulsing white dwarf binary star. Nature 537, 374–377 (2016).

    Article  ADS  CAS  PubMed  Google Scholar 

  11. Hyman, S. D. et al. A powerful bursting radio source towards the Galactic Centre. Nature 434, 50–52 (2005).

    Article  ADS  CAS  PubMed  Google Scholar 

  12. Szary, A., Zhang, B., Melikidze, G. I., Gil, J. & Xu, R.-X. Radio efficiency of pulsars. Astrophys. J. 784, 59 (2014).

    Article  ADS  Google Scholar 

  13. Levin, L. et al. Radio emission evolution, polarimetry and multifrequency single pulse analysis of the radio magnetar PSR J1622–4950. Mon. Not. R. Astron. Soc. 422, 2489–2500 (2012).

    Article  ADS  Google Scholar 

  14. Rea, N., Pons, J. A., Torres, D. F. & Turolla, R. The fundamental plane for radio magnetars. Astrophys. J. Lett. 748, L12 (2012).

    Article  ADS  Google Scholar 

  15. Olausen, S. A. & Kaspi, V. M. The McGill magnetar catalog. Astrophys. J. Suppl. 212, 6 (2014).

    Article  ADS  Google Scholar 

  16. Carbone, D. et al. New methods to constrain the radio transient rate: results from a survey of four fields with LOFAR. Mon. Not. R. Astron. Soc. 459, 3161–3174 (2016).

    Article  ADS  CAS  PubMed  Google Scholar 

  17. Polisensky, E. et al. Exploring the transient radio sky with VLITE: early results. Astrophys. J. 832, 60 (2016).

    Article  ADS  Google Scholar 

  18. Stewart, A. J. et al. LOFAR MSSS: detection of a low-frequency radio transient in 400 h of monitoring of the North Celestial Pole. Mon. Not. R. Astron. Soc. 456, 2321–2342 (2016).

    Article  ADS  Google Scholar 

  19. Bell, M. E. et al. The Murchison Widefield Array Transients Survey (MWATS). A search for low-frequency variability in a bright Southern Hemisphere sample. Mon. Not. R. Astron. Soc. 482, 2484–2501 (2019).

    ADS  Google Scholar 

  20. Hajela, A., Mooley, K. P., Intema, H. T. & Frail, D. A. A GMRT 150 MHz search for variables and transients in Stripe 82. Mon. Not. R. Astron. Soc. 490, 4898–4906 (2019).

    Article  ADS  Google Scholar 

  21. Keith, M. J. et al. The High Time Resolution Universe Pulsar Survey – I. System configuration and initial discoveries. Mon. Not. R. Astron. Soc. 409, 619–627 (2010).

    Article  ADS  Google Scholar 

  22. Cameron, A. D., Barr, E. D., Champion, D. J., Kramer, M. & Zhu, W. W. An investigation of pulsar searching techniques with the fast folding algorithm. Mon. Not. R. Astron. Soc. 468, 1994–2010 (2017).

    Article  ADS  CAS  Google Scholar 

  23. Hurley-Walker, N. et al. Galactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey – I. A low-frequency extragalactic catalogue. Mon. Not. R. Astron. Soc. 464, 1146–1167 (2017).

    Article  ADS  Google Scholar 

  24. Pietka, M., Fender, R. P. & Keane, E. F. The variability time-scales and brightness temperatures of radio flares from stars to supermassive black holes. Mon. Not. R. Astron. Soc. 446, 3687–3696 (2015).

    Article  ADS  CAS  Google Scholar 

  25. Tingay, S. J. et al. The Murchison Widefield Array: the square kilometre array precursor at low radio frequencies. Publ. Astron. Soc. Aust. 30, e007 (2013).

    Article  ADS  Google Scholar 

  26. Wayth, R. B. et al. The Phase II Murchison Widefield Array: design overview. Publ. Astron. Soc. Aust. 35, e033 (2018).

    Article  ADS  Google Scholar 

  27. Wayth, R. B. et al. GLEAM: the Galactic and Extragalactic All-Sky MWA survey. Publ. Astron. Soc. Aust. 32, e025 (2015).

    Article  ADS  Google Scholar 

  28. Offringa, A. R. et al. Parametrizing Epoch of Reionization foregrounds: a deep survey of low-frequency point-source spectra with the Murchison Widefield Array. Mon. Not. R. Astron. Soc. 458, 1057–1070 (2016).

    Article  ADS  Google Scholar 

  29. Hurley-Walker, N. et al. Galactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey II: galactic plane 345° < l < 67°, 180° < l < 240°. Publ. Astron. Soc. Aust. 36, e047 (2019).

    Article  ADS  Google Scholar 

  30. Offringa, A. R. et al. WSCLEAN: an implementation of a fast, generic wide-field imager for radio astronomy. Mon. Not. R. Astron. Soc. 444, 606–619 (2014).

    Article  ADS  Google Scholar 

  31. Sutinjo, A. et al. Understanding instrumental Stokes leakage in Murchison Widefield Array polarimetry. Radio Sci. 50, 52–65 (2015).

    Article  ADS  Google Scholar 

  32. Lenc, E. et al. The challenges of low-frequency radio polarimetry: lessons from the Murchison Widefield Array. Publ. Astron. Soc. Aust. 34, e040 (2017).

    Article  ADS  Google Scholar 

  33. Riseley, C. J. et al. The Polarised GLEAM Survey (POGS) I: first results from a low-frequency radio linear polarisation survey of the southern sky. Publ. Astron. Soc. Aust. 35, e043 (2018).

    Article  ADS  Google Scholar 

  34. Riseley, C. J. et al. The Polarised GLEAM Survey (POGS) II: results from an all-sky rotation measure synthesis survey at long wavelengths. Publ. Astron. Soc. Aust. 37, e029 (2020).

    Article  ADS  Google Scholar 

  35. Brentjens, M. A. & De Bruyn, A. Faraday rotation measure synthesis. Astron. Astrophys. 441, 1217–1228 (2005).

    Article  ADS  Google Scholar 

  36. Lafler, J. & Kinman, T. D. An RR Lyrae star survey with Ihe Lick 20-INCH Astrograph II. The calculation of RR Lyrae periods by electronic computer. Astrophys. J. Suppl. 11, 216 (1965).

    Article  ADS  Google Scholar 

  37. Stellingwerf, R. F. Period determination using phase dispersion minimization. Astrophys. J. 224, 953–960 (1978).

    Article  ADS  Google Scholar 

  38. Schwarzenberg-Czerny, A. Fast and statistically optimal period search in uneven sampled observations. Astrophys. J. Lett. 460, L107 (1996).

    Article  ADS  Google Scholar 

  39. Clarke, D. String/Rope length methods using the Lafler-Kinman statistic. Astron. Astrophys. 386, 763–774 (2002).

    Article  ADS  Google Scholar 

  40. Manchester, R. N., Hobbs, G. B., Teoh, A. & Hobbs, M. The Australia Telescope National Facility Pulsar Catalogue. Astron. J 129, 1993–2006 (2005).

    Article  ADS  Google Scholar 

  41. Hurley-Walker, N. & Hancock, P. J. De-distorting ionospheric effects in the image plane. Astron. Comput. 25, 94–102 (2018).

    Article  ADS  Google Scholar 

  42. Gehrels, N. et al. The Swift gamma-ray burst mission. Astrophys. J. 611, 1005–1020 (2004).

    Article  ADS  CAS  Google Scholar 

  43. Burrows, D. N. et al. The Swift X-ray telescope. Space Sci. Rev. 120, 165–195 (2005).

    Article  ADS  Google Scholar 

  44. NASA High Energy Astrophysics Science Archive Research Center (HEASARC). HEAsoft: Unified Release of FTOOLS and XANADU (2014).

  45. Champion, D. et al. High-cadence observations and variable spin behaviour of magnetar Swift J1818.0–1607 after its outburst. Mon. Not. R. Astron. Soc. 498, 6044–6056 (2020).

    Article  ADS  CAS  Google Scholar 

  46. Lorimer, D. R. & Kramer, M. Handbook of Pulsar Astronomy (Cambridge Univ. Press, 2012).

  47. Zhang, B., Harding, A. K. & Muslimov, A. G. Radio pulsar death line revisited: is PSR J2144–3933 anomalous? Astrophys. J. Lett. 531, L135–L138 (2000).

    Article  ADS  CAS  Google Scholar 

Download references

Acknowledgements

N.H.-W. is the recipient of an Australian Research Council Future Fellowship (project number FT190100231) and G.E.A. is the recipient of an Australian Research Council Discovery Early Career Researcher Award (project number DE180100346) funded by the Australian Government. This scientific work makes use of the Murchison Radio-astronomy Observatory, operated by CSIRO. We acknowledge the Wajarri Yamatji people as the traditional owners of the observatory site. Support for the operation of the Murchison Widefield Array is provided by the Australian Government (NCRIS), under a contract to Curtin University administered by Astronomy Australia Limited. We acknowledge the Pawsey Supercomputing Centre, which is supported by the Western Australian Government and the Australian Government. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. This project was supported by resources and expertise provided by CSIRO IMT Scientific Computing. This work used resources of China SKA Regional Centre prototype funded by the National Key R&D Programme of China (2018YFA0404603) and Chinese Academy of Sciences (114231KYSB20170003).

Author information

Authors and Affiliations

Authors

Contributions

N.H.-W. calibrated and processed the data for the observations described herein, determined the position and flux density of the source and prepared the manuscript, with contributions from all co-authors. X.Z. processed all polarization data, including performing polarization calibration and analysis. A.B. and S.J.M. performed the analysis to derive the period and period derivative. A.B. performed the X-ray observations and analysis. S.J.M. calculated and applied the dispersion measure and barycentric corrections. T.N.O. developed the original detection methodology, performed the original archive search and made the initial discovery. P.J.H. helped to develop the detection methodology and provided supercomputing support. G.E.A. contributed to astrophysical calculations and interpretation of the data. T.J.G., G.H.H., J.S.M. and X.Z. determined polarization calibration methods. T.J.G. performed early refinement of the period estimate.

Corresponding author

Correspondence to N. Hurley-Walker.

Ethics declarations

Competing interests

The authors declare no competing interests.

Peer review information

Nature thanks the anonymous reviewers for their contribution to the peer review of this work. Peer reviewer reports are available.

Additional information

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Extended data figures and tables

Extended Data Fig. 1 Dynamic spectra of the observations used to calculate the dispersion measure.

From top to bottom, separated by horizontal black lines, we show the observations 1205008192 (72–103 MHz), 1205007112 (103–134 MHz), 1205011432 (139–170 MHz), 1205010352 (170–200 MHz) and 1205009272 (200–231 MHz). These observations were taken on 13 March 2018 between 20:12 and 21:24 (see Fig. 1). The left panel shows the data aligned using a period of 1091.1690 s, whereas the right panel shows the same, including a dispersion correction of 57 pc cm−3. Strong ionospheric scintillation is visible in the 72–103-MHz data, causing ripples in the brightness of the source over time.

Extended Data Fig. 2 The explored search space in P and \(\dot{{\boldsymbol{P}}}\) for the pulses recorded from GLEAM-X J162759.5-523504.3.

The contours show the peak flux density of the mean profile at 154 MHz recovered at each combination of P (y axis) and \(\dot{P}\) (x axis), in levels of 15, 14, 13, 12 and 11 Jy. The best-fit values of P = 1,091.1690 s and \(\dot{P}=6\times {10}^{-10}{{\rm{s}}{\rm{s}}}^{-1}\) are marked with a dark red ‘+’.

Extended Data Fig. 3 The flux density of GLEAM-X J162759.5-523504.3 as a function of frequency.

This is derived from the same observations shown in Extended Data Fig. 1. Points are determined via an average of the source profile in each frequency bin, weighting by the signal-to-noise ratio of the frequency-averaged profile. A power-law fit using the data spanning 95–195 MHz is shown in blue, with α = −1.16 ± 0.04.

Extended Data Fig. 4 X-ray luminosity and spectral properties of magnetars15 compared with the X-ray luminosity limits of GLEAM-X J162759.5-523504.3.

The two faintest known magnetars are labelled. The coloured contours represent expected Swift XRT count rates for putative X-ray luminosities and blackbody temperatures for a source at 1.3 kpc, assuming a hydrogen column density of NH ≈ 2 × 1021 cm−2. The grey dashed line represents the implied luminosity upper limit for GLEAM-X J162759.5-523504.3 based on the 3 − σ upper limit obtained from the Swift XRT observation. From the magnetar fundamental plane, we predict GLEAM-X J162759.5-523504.3 to have a quiescent luminosity 4.5 orders of magnitude lower than our limit from the Swift XRT.

Extended Data Fig. 5 Pulse profiles for the three detections on 3 January 2018 compared with a wide pulse detected on 14 March 2018.

Barycentric corrections and dedispersion have been applied. The data are all taken at the same frequency, 170–200 MHz. Vertical dashed lines encapsulating the profile found on 13 March 2018 are overplotted to guide the eye.

Extended Data Fig. 6 Images in full Stokes and polarized intensity (in Jy beam−1), and RM (in rad m−2), of the region 5° × 5° around GLEAM-X J162759.5-523504.3.

The images were made using observation 1200354592 on 18 January 2018 at 23:49, using only the interval where the source was producing emission. Faraday rotation over the imaged bandwidth of 30 MHz causes the Stokes Q, U and V emission to average to zero. The polarized intensity shows the maximum value of the RM spectrum. Where the polarized intensity is less than seven times the local noise, the corresponding RM value has been masked.

Extended Data Fig. 7 A scatter plot of period derivative \(\dot{{\boldsymbol{P}}}\) against period P.

GLEAM-X J162759.5-523504.3 is shown as a blue arrow with an upper limit on \(\dot{P}\) (see Methods) in context with the known pulsars40 (black dots), X-ray-detected magnetars15 (red dots and arrows) and magnetars known to emit in both X-ray and radio frequencies (red circles around black dots). The slowest (and radio-quiet) X-ray magnetar 1E 161348–5055 is also shown with an upper limit on \(\dot{P}\). The green dashed and dot-dashed lines correspond to the theoretical ‘death lines’ for pulsar radio emission for cases I and III calculated by Zhang et al.47.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hurley-Walker, N., Zhang, X., Bahramian, A. et al. A radio transient with unusually slow periodic emission. Nature 601, 526–530 (2022). https://doi.org/10.1038/s41586-021-04272-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/s41586-021-04272-x

This article is cited by

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Search

Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing