Search for Continuous Gravitational Wave Signals in Pulsar Timing Residuals: A New Scalable Approach with Diffusive Nested Sampling

Abstract

Detecting continuous nanohertz gravitational waves (GWs) generated by individual close binaries of supermassive black holes (CB-SMBHs) is one of the primary objectives of pulsar timing arrays (PTAs). The detection sensitivity is slated to increase significantly as the number of well-timed millisecond pulsars will increase by more than an order of magnitude with the advent of next-generation radio telescopes. Currently, the Bayesian analysis pipeline using parallel tempering Markov chain Monte Carlo has been applied in multiple studies for CB-SMBH searches, but it may be challenged by the high dimensionality of the parameter space for future large-scale PTAs. One solution is to reduce the dimensionality by maximizing or marginalizing over uninformative parameters semi-analytically, but it is not clear whether this approach can be extended to more complex signal models without making overly simplified assumptions. Recently, the method of diffusive nested (DNest) sampling shown the capability of coping with high dimensionality and multimodality effectively in Bayesian analysis. In this paper, we apply DNest to search for continuous GWs in simulated pulsar timing residuals and find that it performs well in terms of accuracy, robustness, and efficiency for a PTA including $\mathcal{O}(10^2)$ pulsars. DNest also allows a simultaneous search of multiple sources elegantly, which demonstrates its scalability and general applicability. Our results show that it is convenient and also high beneficial to include DNest in current toolboxes of PTA analysis.

arXiv comment20 pages, 10 figures, submitted to ApJ

References in corpus (45)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Observation of Gravitational Waves from a Binary Black Hole Merger
• GWTC-1: A Gravitational-Wave Transient Catalog of Compact Binary Mergers Observed by LIGO and Virgo during the First and Second Observing Runs
• GW170817: Measurements of Neutron Star Radii and Equation of State
• TianQin: a space-borne gravitational wave detector
• Using gravitational-wave standard sirens
• The NANOGrav 12.5-year Data Set: Search For An Isotropic Stochastic Gravitational-Wave Background
• Tempo2, a new pulsar timing package. II: The timing model and precision estimates
• A new route towards merging massive black holes
• High-precision timing of 42 millisecond pulsars with the European Pulsar Timing Array
• Fundamental Physics with the Square Kilometre Array
• Upper bounds on the low-frequency stochastic gravitational wave background from pulsar timing observations: current limits and future prospects
• A systematic search for close supermassive black hole binaries in the Catalina Real-Time Transient Survey
• Modeling reverberation mapping data II: dynamical modeling of the Lick AGN Monitoring Project 2008 dataset
• A Population of Short-Period Variable Quasars from PTF as Supermassive Black Hole Binary Candidates
• Markov Chain Monte Carlo Methods for Bayesian Data Analysis in Astronomy
• A Large Systematic Search for Recoiling and Close Supermassive Binary Black Holes
• European Pulsar Timing Array Limits on Continuous Gravitational Waves from Individual Supermassive Black Hole Binaries
• Pulsar searches and timing with the SKA
• PINT: A Modern Software Package for Pulsar Timing
• The NANOGrav 11-Year Data Set: Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries
• The NANOGrav 12.5 yr Data Set: Observations and Narrowband Timing of 47 Millisecond Pulsars
• An all-sky search for continuous gravitational waves in the Parkes Pulsar Timing Array data set
• A study of spatial correlations in pulsar timing array data
• The Parkes Pulsar Timing Array Project: Second data release
• Resolving multiple supermassive black hole binaries with pulsar timing arrays
• Optimal strategies for continuous gravitational wave detection in pulsar timing arrays
• Spectroscopic Indication for a Centi-parsec Supermassive Black Hole Binary in the Galactic Center of NGC 5548
• The First Spectroscopically Resolved Sub-parsec Orbit of a Supermassive Binary Black Hole
• Supermassive Black Holes with High Accretion Rates in Active Galactic Nuclei. VIII. Structure of the Broad-Line Region and Mass of the Central Black Hole in Mrk 142
• Pulsar Timing Array Based Search for Supermassive Black Hole Binaries in the Square Kilometer Array Era
• A Bayesian analysis pipeline for continuous GW sources in the PTA band
• A parallax distance to 3C 273 through spectroastrometry and reverberation mapping
• A Possible ~20 yr Periodicity in Long-term Optical Photometric and Spectral Variations of the Nearby Radio-Quiet Active Galactic Nucleus Ark 120
• Detection and localization of continuous gravitational waves with pulsar timing arrays: the role of pulsar terms
• Accelerated Bayesian model-selection and parameter-estimation in continuous gravitational-wave searches with pulsar-timing arrays
• Resolving multiple supermassive black hole binaries with pulsar timing arrays II: genetic algorithm implementation
• Joint search for isolated sources and an unresolved confusion background in pulsar timing array data
• Particle Swarm Optimization and gravitational wave data analysis: Performance on a binary inspiral testbed
• Kinematic signatures of reverberation mapping of close binaries of supermassive black holes in active galactic nuclei
• Multimessenger pulsar timing array constraints on supermassive black hole binaries traced by periodic light curves
• Extending the frequency reach of pulsar timing array based gravitational wave search without high cadence observations
• Differential Interferometric Signatures of Close Binaries of Supermassive Black Holes in Active Galactic Nuclei
• Continuous gravitational wave searches with pulsar timing arrays: Maximization versus marginalization over pulsar phase parameters
• Prospects for gravitational wave astronomy with next generation large-scale pulsar timing arrays

Cited by in corpus (1)

• Iterative time-domain method for resolving multiple gravitational wave sources in Pulsar Timing Array data