- Physical Review Journals
All JournalsPhysics Magazine

Physical Review D

Reuse & Permissions

It is not necessary to obtain permission to reuse this article or its components as it is available under the terms of the Creative Commons Attribution 4.0 International license. This license permits unrestricted use, distribution, and reproduction in any medium, provided attribution to the author(s) and the published article's title, journal citation, and DOI are maintained. Please note that some figures may have been included with permission from other third parties. It is your responsibility to obtain the proper permission from the rights holder directly for these figures.

Export citation

Export citation

Choose format for download:

Download Citation
  • Open Access

Modeling non-Gaussianities in pulsar timing array data analysis using Gaussian mixture models

Mikel Falxa* and Alberto Sesana

  • *Contact author: mikel.falxa@unimib.it
PDF

Phys. Rev. D 113, 043047 – Published 23 February, 2026

DOI: https://doi.org/10.1103/9jdd-6ct8

Abstract

In pulsar timing array (PTA) data analysis, noise is typically assumed to be Gaussian, and the marginalized likelihood has a well-established analytical form derived within the framework of Gaussian processes. However, this Gaussianity assumption may break down for certain classes of astrophysical and cosmological signals, particularly for a gravitational wave background (GWB) generated by a population of supermassive black hole binaries (SMBHBs). In this work, we present a new method for testing the presence of non-Gaussian features in PTA data. We go beyond the Gaussian assumption by modeling the noise or signal statistics using a Gaussian mixture model (GMM). An advantage of this approach is that the marginalization of the likelihood remains fully analytical, expressed as a linear combination of Gaussian PTA likelihoods. This makes the method straightforward to implement within existing data analysis tools. Moreover, this method extends beyond the free spectrum analysis by producing posterior probability distributions of higher-order moments inferred from the data, which can be incorporated into spectral refitting techniques. We validate the model using simulations and demonstrate the sensitivity of PTAs to non-Gaussianity by computing the Bayes factor in favor of the GMM as a function of the injected excess moments. We apply the method to a more astrophysically motivated scenario where a single SMBHB is resolved on top of a Gaussian GWB and show that significant non-Gaussianities are introduced by the individual source. Finally, we test our model on a realistic GWB generated from a simulated population of SMBHBs.

View figure in article

Physics Subject Headings (PhySH)

Article Text

References (47)

  1. M. V. Sazhin, Opportunities for detecting ultralong gravitational waves, Sov. Astron. 22, 36 (1978), https://ui.adsabs.harvard.edu/abs/1978AZh....55...65S/abstract.
  2. S. Detweiler, Pulsar timing measurements and the search for gravitational waves, Astrophys. J. 234, 1100 (1979).
  3. R. S. Foster and D. C. Backer, Constructing a pulsar timing array, Astrophys. J. 361, 300 (1990).
  4. A. Sesana, A. Vecchio, and C. N. Colacino, The stochastic gravitational-wave background from massive black hole binary systems: Implications for observations with pulsar timing arrays, Mon. Not. R. Astron. Soc. 390, 192 (2008).
  5. A. Sesana, Insights into the astrophysics of supermassive black hole binaries from pulsar timing observations, Classical Quantum Gravity 30, 224014 (2013).
  6. C. Caprini and D. G. Figueroa, Cosmological backgrounds of gravitational waves, Classical Quantum Gravity 35, 163001 (2018).
  7. R. W. Hellings and G. S. Downs, Upper limits on the isotropic gravitational radiation background from pulsar timing analysis, Astrophys. J. Lett. 265, L39 (1983).
  8. B. Allen, Variance of the Hellings-Downs correlation, Phys. Rev. D 107, 043018 (2023).
  9. J. Antoniadis et al., The second data release from the European pulsar timing array: III. Search for gravitational wave signals, Astron. Astrophys. 678, A50 (2023).
  10. G. Agazie et al., The nanograv 15 yr data set: Evidence for a gravitational-wave background, Astrophys. J. Lett. 951, L8 (2023).
  11. D. J. Reardon et al., Search for an isotropic gravitational-wave background with the parkes pulsar timing array, Astrophys. J. Lett. 951, L6 (2023).
  12. M. T. Miles et al., The meerkat pulsar timing array: The first search for gravitational waves with the meerkat radio telescope, Mon. Not. R. Astron. Soc. 536, 1489 (2024).
  13. H. Xu et al., Searching for the nano-hertz stochastic gravitational wave background with the Chinese pulsar timing array data release I, Res. Astron. Astrophys. 23, 075024 (2023).
  14. R. van Haasteren and M. Vallisneri, New advances in the gaussian-process approach to pulsar-timing data analysis, Phys. Rev. D 90, 104012 (2014).
  15. N. Laal et al., Deep neural emulation of the supermassive black-hole binary population, Astrophys. J. 982, 55 (2025).
  16. L. Lentati, M. P. Hobson, and P. Alexander, Bayesian estimation of non-Gaussianity in pulsar timing analysis, Mon. Not. R. Astron. Soc. 444, 3863 (2014).
  17. U. von Toussaint, Bayesian inference in physics, Rev. Mod. Phys. 83, 943 (2011).
  18. M. Crisostomi et al., Beyond diagonal approximations: Improved covariance modeling for pulsar timing array data analysis, arXiv:2506.13866.
  19. G. Agazie et al., The nanograv 15 yr data set: Bayesian limits on gravitational waves from individual supermassive black hole binaries, Astrophys. J. Lett. 951, L50 (2023).
  20. N. Seto, Non-Gaussianity analysis of a gravitational wave background made by short-duration burst signals, Phys. Rev. D 80, 043003 (2009).
  21. B. Allen, D. Agarwal, J. D. Romano, and S. Valtolina, Source anisotropies and pulsar timing arrays, Phys. Rev. D 110, 123507 (2024).
  22. X. Xue, Z. Pan, and L. Dai, Non-Gaussian statistics of nanohertz stochastic gravitational waves, Phys. Rev. D 111, 043022 (2025).
  23. W. G. Lamb and S. R. Taylor, Spectral variance in a stochastic gravitational-wave background from a binary population, Astrophys. J. Lett. 971, L10 (2024).
  24. N. Laal, S. R. Taylor, R. van Haasteren, W. G. Lamb, and X. Siemens, Solving the PTA data analysis problem with a global Gibbs scheme, Phys. Rev. D 111, 063067 (2025).
  25. W. G. Lamb, S. R. Taylor, and R. van Haasteren, Rapid refitting techniques for Bayesian spectral characterization of the gravitational wave background using pulsar timing arrays, Phys. Rev. D 108, 103019 (2023).
  26. A. Mitridate et al., Ptarcade, arXiv:2306.16377.
  27. H. Quelquejay Leclere et al., Practical approaches to analyzing PTA data: Cosmic strings with six pulsars, Phys. Rev. D 108, 123527 (2023).
  28. J. A. Ellis, M. Vallisneri, S. R. Taylor, and P. T. Baker, Enterprise: Enhanced numerical toolbox enabling a robust pulsar inference suite, Zenodo, 10.5281/zenodo.4059815 (2020).
  29. S. Valtolina and R. van Haasteren, Regularizing the pulsar timing array likelihood: A path towards Fourier space, Phys. Rev. D 112, 043046 (2025).
  30. EPTA Collaboration and InPTA Collaboration, The second data release from the European pulsar timing array—IV. Implications for massive black holes, dark matter, and the early universe, Astron. Astrophys. 685, A94 (2024).
  31. B. P. Carlin and S. Chib, Bayesian model choice via Markov Chain Monte Carlo methods, J. R. Stat. Soc. Ser. B 57, 473 (1995).
  32. S. Hee et al., Bayesian model selection without evidences: Application to the dark energy equation-of-state, Mon. Not. R. Astron. Soc. 455, 2461 (2015).
  33. Z. Arzoumanian et al., The nanograv 11 year data set: Pulsar-timing constraints on the stochastic gravitational-wave background, Astrophys. J. 859, 47 (2018).
  34. J. S. Hazboun et al., The nanograv 11 yr data set: Evolution of gravitational-wave background statistics, Astrophys. J. 890, 108 (2020).
  35. É. Racine and C. Cutler, Gaussianity of LISA’s confusion backgrounds, Phys. Rev. D 76, 124033 (2007).
  36. A. Sesana, A. Vecchio, and M. Volonteri, Gravitational waves from resolvable massive black hole binary systems and observations with Pulsar Timing Arrays, Mon. Not. R. Astron. Soc. 394, 2255 (2009).
  37. S. Babak et al., European pulsar timing array limits on continuous gravitational waves from individual supermassive black hole binaries, Mon. Not. R. Astron. Soc. 455, 1665 (2015).
  38. J. Antoniadis et al., The second data release from the European pulsar timing array: V. Search for continuous gravitational wave signals, Astron. Astrophys. 690, A118 (2024).
  39. G. Agazie et al., The nanograv 15 yr data set: Bayesian limits on gravitational waves from individual supermassive black hole binaries, Astrophys. J. Lett. 951, L50 (2023).
  40. M. Falxa et al., Searching for continuous gravitational waves in the second data release of the international pulsar timing array, Mon. Not. R. Astron. Soc. 521, 5077 (2023).
  41. L. Boyle and U.-L. Pen, Pulsar timing arrays as imaging gravitational wave telescopes: Angular resolution and source (de)confusion, Phys. Rev. D 86, 124028 (2012).
  42. I. Ferranti, G. Shaifullah, A. Chalumeau, and A. Sesana, Separating deterministic and stochastic gravitational wave signals in realistic pulsar timing array datasets, Astron. Astrophys. 694, A194 (2025).
  43. J. S. Hazboun, J. D. Romano, and T. L. Smith, Realistic sensitivity curves for pulsar timing arrays, Phys. Rev. D 100, 104028 (2019).
  44. C. R. Contaldi, P. G. Ferreira, J. Magueijo, and K. M. Górski, A Bayesian estimate of the skewness of the cosmic microwave background, Astrophys. J. 534, 25 (2000).
  45. F.-S. Kitaura, Non-Gaussian gravitational clustering field statistics, Mon. Not. R. Astron. Soc. 420, 2737 (2012).
  46. S. Berkowitz and F. J. Garner, The calculation of multidimensional Hermite polynomials and gram-charlier coefficients, Math. Comput. 24, 537 (1970).
  47. S. Cook et al., Validation of software for Bayesian models using posterior quantiles, J. Comput. Graph. Stat. 15, 675 (2006).

Outline

Information

Phys. Rev. D 113, 043047– Published 23 February, 2026
  • Received 13 August 2025
  • Accepted 28 January 2026
Crossmark - Check for updates button
Creative Commons Logo

Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article’s title, journal citation, and DOI.

Published by the American Physical Society

Sign In to Your Journals Account

Filter

Filter

Article Lookup

Enter a citation