Targeted search for eccentric supermassive binary black holes in OJ 287 and nearby galaxy clusters with PPTA DR3
Pith reviewed 2026-05-10 14:49 UTC · model grok-4.3
The pith
Pulsar timing array data sets upper limits on eccentric supermassive black hole binary masses in OJ 287 and nearby clusters, excluding equal-mass systems at 10 nHz.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We perform Bayesian targeted searches for continuous gravitational waves from eccentric supermassive binary black holes using the Parkes Pulsar Timing Array third data release. Six electromagnetically motivated sky directions are analyzed, including the blazar OJ 287 and five nearby galaxy clusters. No significant signals are found. For OJ 287, by explicitly incorporating orbital eccentricity up to e0 = 0.8 to robustly capture signal power spread across multiple harmonics, we constrain the total binary mass to M_tot < 5.25 × 10^10 M_⊙ at 95 percent credible level. We also place upper limits on the chirp mass of potential SMBBHs residing in galaxy clusters. By combining these limits with独立黑洞质
What carries the argument
Bayesian targeted search that includes orbital eccentricity in the waveform model up to e0 = 0.8 to account for gravitational-wave power distributed across multiple harmonics.
If this is right
- The total mass of any SMBBH in OJ 287 is limited to below 5.25 × 10^10 solar masses.
- Binaries with mass ratios q ≳ 10^{-2} are excluded at frequencies around 10 nHz in the sampled massive systems.
- Equal-mass black hole mergers are ruled out in the parameter space probed for these galaxy-cluster hosts.
- Pulsar timing arrays can place meaningful constraints on SMBBH populations even without a detection.
Where Pith is reading between the lines
- Extending the same targeted search method to additional electromagnetically identified SMBBH candidates could tighten population statistics.
- Future PTA data releases with improved sensitivity may either detect the excluded signals or push the mass-ratio limits lower.
- The mass-ratio exclusions could be cross-checked against predictions from galaxy-merger simulations for the same clusters.
Load-bearing premise
The chosen sky directions are accurate and the eccentricity model up to 0.8 fully captures signal power without missing harmonics or being biased by the PPTA DR3 noise model.
What would settle it
A future detection of a continuous gravitational-wave signal whose frequency evolution and sky location match an eccentric SMBBH in one of the six targeted directions would falsify the current non-detection and upper limits.
Figures
read the original abstract
We perform Bayesian targeted searches for continuous gravitational waves from eccentric supermassive binary black holes (SMBBHs) using the Parkes Pulsar Timing Array third data release (PPTA DR3). Six electromagnetically motivated sky directions are analyzed, including the blazar OJ~287 and five nearby galaxy clusters (Virgo, Fornax, Norma, Hercules, and Coma). No significant signals are found. For OJ 287, by explicitly incorporating orbital eccentricity (up to $e_0 = 0.8$) to robustly capture signal power spread across multiple harmonics, we constrain the total binary mass to $M_{\rm tot} < 5.25 \times 10^{10} M_{\odot}$ (95\% credible level). We also place upper limits on the chirp mass of potential SMBBHs residing in galaxy clusters. By combining these limits with independent black hole mass estimates, we place novel constraints on the allowed binary mass ratios for potential hosts such as M87 and NGC~4889. Specifically, our results exclude binaries with mass ratios $q \gtrsim 10^{-2}$ at around 10 nHz for these massive systems, effectively ruling out equal-mass black hole mergers in the sampled parameter space. These findings demonstrate the growing power of pulsar timing arrays to probe SMBBH populations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports Bayesian targeted searches for continuous gravitational waves from eccentric supermassive binary black holes using the PPTA DR3 timing residuals. Six electromagnetically motivated sky positions are analyzed (OJ 287 plus Virgo, Fornax, Norma, Hercules, and Coma clusters). No significant signals are detected. For OJ 287 the analysis incorporates orbital eccentricity up to e0 = 0.8 to spread signal power across harmonics, yielding a 95% credible upper limit M_tot < 5.25 × 10^10 M_⊙. Upper limits on chirp mass are placed for the clusters; when combined with independent black-hole mass estimates these exclude mass ratios q ≳ 10^{-2} near 10 nHz for systems such as M87 and NGC 4889, thereby ruling out equal-mass mergers in the sampled parameter space.
Significance. If the upper limits are robust, the work provides concrete, observationally grounded constraints on the allowed parameter space for SMBBHs in well-studied systems and demonstrates the utility of explicitly eccentric waveform models in PTA analyses. The combination of PTA non-detections with electromagnetic mass estimates to bound mass ratios is a useful addition to the literature on SMBBH populations.
major comments (2)
- [§3.2] §3.2 and Eq. (8): the statement that the eccentric waveform model with e0 ≤ 0.8 fully captures signal power requires explicit verification that higher-order harmonics (n > 5) remain negligible at the frequencies and eccentricities considered; otherwise the reported amplitude upper limits could be biased high.
- [§4.3] §4.3, Table 2: the 95% credible limits on M_tot for OJ 287 and on chirp mass for the clusters are derived under the assumption that the PPTA DR3 noise model (including red-noise parameters) is correctly specified at the target sky positions; any unmodeled common-mode or spatially correlated noise could systematically loosen these limits.
minor comments (3)
- [Abstract] The abstract and §5 should state the exact frequency range over which the q ≳ 10^{-2} exclusion applies, rather than the approximate phrase 'around 10 nHz'.
- [Figure 3] Figure 3 caption should clarify whether the plotted credible intervals are marginalized over all other parameters or conditioned on fixed sky position and eccentricity.
- [§2] A brief statement in §2 on the adopted priors for eccentricity and mass ratio would improve reproducibility.
Simulated Author's Rebuttal
We thank the referee for their constructive review and positive recommendation for minor revision. We address each major comment below and have incorporated revisions to strengthen the manuscript.
read point-by-point responses
-
Referee: [§3.2] §3.2 and Eq. (8): the statement that the eccentric waveform model with e0 ≤ 0.8 fully captures signal power requires explicit verification that higher-order harmonics (n > 5) remain negligible at the frequencies and eccentricities considered; otherwise the reported amplitude upper limits could be biased high.
Authors: We thank the referee for this observation. Our eccentric waveform implementation includes harmonics up to n=5, which we selected to capture the dominant power for e0 ≤ 0.8 at the nanohertz frequencies relevant to PPTA DR3. To address the request for explicit verification, we have performed additional calculations of the fractional power contained in n > 5 harmonics across the relevant frequency and eccentricity range. These show that the contribution is negligible (< 0.5% for e0 = 0.8 near 10 nHz). We will add a brief verification paragraph and a supporting plot to §3.2 in the revised manuscript. revision: yes
-
Referee: [§4.3] §4.3, Table 2: the 95% credible limits on M_tot for OJ 287 and on chirp mass for the clusters are derived under the assumption that the PPTA DR3 noise model (including red-noise parameters) is correctly specified at the target sky positions; any unmodeled common-mode or spatially correlated noise could systematically loosen these limits.
Authors: We agree that the reported upper limits are conditional on the PPTA DR3 noise model. The analysis uses the published red-noise parameters for each pulsar as provided in the data release. While unmodeled spatially correlated noise could in principle affect the results, the targeted searches at specific sky positions and the multi-pulsar combination provide some robustness against common-mode effects. We will revise the discussion in §4.3 to explicitly note this assumption and the conditional nature of the limits, including a short statement on potential impacts from additional noise components. revision: yes
Circularity Check
No significant circularity identified
full rationale
The paper reports a standard Bayesian targeted search for continuous-wave signals from eccentric SMBBHs in PPTA DR3 timing residuals at six fixed sky positions. No detections are claimed; upper limits (e.g., M_tot < 5.25e10 Msun at 95% for OJ287, and q ≳ 10^{-2} exclusions at ~10 nHz when combined with external BH mass estimates) follow directly from the posterior under the signal-plus-noise model. Eccentricity (e0 ≤ 0.8) is an explicit parameter of the waveform template that spreads power across harmonics; it is not fitted from the data and then used to 'predict' the same quantity. Mass-ratio bounds use independent galaxy mass catalogs, not internal fits. No self-definitional steps, no fitted inputs renamed as predictions, no load-bearing self-citations, and no ansatz or uniqueness claims that reduce the central result to its own inputs. The derivation chain is self-contained against the PTA data likelihood and external priors.
Axiom & Free-Parameter Ledger
axioms (2)
- standard math General relativity accurately describes the gravitational wave emission from eccentric SMBBHs
- domain assumption PPTA DR3 timing residuals are correctly characterized by the adopted noise model
Reference graph
Works this paper leans on
-
[1]
Detweiler S 1979 ApJ 234 1100
work page 1979
-
[2]
Hellings R W and Downs G S 1983 ApJ Lett. 265 L39
work page 1983
-
[3]
Manchester R N, Hobbs G, Bailes M, Coles W A, van Straten W, Keith M J, Shannon R M, Bhat N D R, Brown A, Burke-Spolaor S G, Champion D J, Chaudhary A, Ed- wards R T, Hampson G, Hotan A W, Jameson A, Jenet F A, Kesteven M J, Khoo J, Kocz J, Maciesiak K, Oslowski S, Ravi V, Reynolds J R, Sarkissian J M, Verbiest J P W, Wen Z L, Wilson W E, Yardley D, Yan W...
work page 2013
-
[4]
Agazie G, Anumarlapudi A, Archibald A M, Arzoumanian Z, Baker P T, Bécsy B, Blecha L, Brazier A, Brook P R, Burke-Spolaor S, Burnette R, Case R, Charisi M, Chatter- jee S, Chatziioannou K, Cheeseboro B D, Chen S, Cohen T, Cordes J M, Cornish N J, Crawford F, Cromartie H T, Crowter K, Cutler C J, Decesar M E, Degan D, Demorest P B, Deng H, Dolch T, Drachle...
work page 2023
-
[5]
Antoniadis J, Arumugam P, Arumugam S, Babak S, Bagchi M, Bak Nielsen A S, Bassa C G, Bathula A, Berthereau A, Bonetti M, Bortolas E, Brook P R, Bur- gay M, Caballero R N, Chalumeau A, Champion D J, Chanlaridis S, Chen S, Cognard I, Dandapat S, Deb D, Desai S, Desvignes G, Dhanda-Batra N, Dwivedi C, Falxa M, Ferdman R D, Franchini A, Gair J R, Goncharov B,...
work page 2023
-
[6]
Reardon D J, Zic A, Shannon R M, Hobbs G B, Bailes M, Di Marco V, Kapur A, Rogers A F, Thrane E, Askew J, Bhat N D R, Cameron A, Curylo M, Coles W A, Dai S, Goncharov B, Kerr M, Kulkarni A, Levin Y, Lower M E, Manchester R N, Mandow R, Miles M T, Nathan R S, Oslowski S, Russell C J, Spiewak R, Zhang S and Zhu X J 2023 ApJ Lett. 951 L6
work page 2023
-
[7]
Xu H, Chen S, Guo Y, Jiang J, Wang B, Xu J, Xue Z, Caballero R N, Yuan J, Xu Y, Wang J, Hao L, Luo J, Lee K, Han J, Jiang P, Shen Z, Wang M, Wang N, Xu R, Wu X, Manchester R, Qian L, Guan X, Huang M, Sun C, Zhu Y and CPTA Collaboration 2023 Res. Astron. Astrophys. 23 075024
work page 2023
-
[8]
Miles M T, Shannon R M, Reardon D J, Bailes M, Cham- pion D J, Geyer M, Gitika P, Grunthal K, Keith M J, Kramer M, Kulkarni A D, Nathan R S, Parthasarathy A, Singha J, Theureau G, Thrane E, Abbate F, Buchner S, Cameron A D, Camilo F, Moreschi B E, Shaifullah G, Shamohammadi M, Possenti A and Venkatraman Krish- nan V 2025 MNRAS 536 1489
work page 2025
-
[9]
Agazie G, Anumarlapudi A, Archibald A M, Arzoumanian Z, Baker P T, Bécsy B, Blecha L, Brazier A, Brook P R, Burke-Spolaor S, Case R, Casey-Clyde J A, Charisi M, Chatterjee S, Cohen T, Cordes J M, Cornish N J, Craw- ford F, Cromartie H T, Crowter K, DeCesar M, Demor- est P B, Digman M C, Dolch T, Drachler B, Ferrara E C, Fiore W, Fonseca E, Freedman G E, G...
work page 2026
-
[10]
Falxa M, Babak S, Baker P T, Bécsy B, Chalumeau A, Chen S, Chen Z, Cornish N J, Guillemot L, Hazboun J S, Mingarelli C M F, Parthasarathy A, Petiteau A, Pol N S, Sesana A, Spolaor S B, Taylor S R, Theureau G, Val- lisneri M, Vigeland S J, Witt C A, Zhu X, Antoniadis J, Arzoumanian Z, Bailes M, Bhat N D R, Blecha L, Brazier A, Brook P R, Caballero N, Camer...
work page 2023
-
[11]
Antoniadis J, Arumugam P, Arumugam S, Babak S, Bagchi M, Bak Nielsen A S, Bassa C G, Bathula A, Berthereau A, Bonetti M, Bortolas E, Brook P R, Burgay M, Caballero R N, Chalumeau A, Champion D J, Chan- laridis S, Chen S, Cognard I, Dandapat S, Deb D, Desai S, Desvignes G, Dhanda-Batra N, Dwivedi C, Falxa M, Fer- ranti I, Ferdman R D, Franchini A, Gair J R...
work page 2024
-
[12]
Zhao S Y, Chen Z C, Cardinal Tremblay J, Goncharov B, Zhu X J, Bhat N D R, Curylo M, Dai S, Di Marco V, Ding H, Hobbs G, Kapur A, Ling W, Liu T, Mandow R, Mishra S, Reardon D J, Russell C J, Shannon R M, Wang S, Zhang L and Zic A 2025 ApJ 992 181
work page 2025
-
[13]
Cardinal Tremblay J, Goncharov B, van Haasteren R, Bhat N D R, Chen Z C, Di Marco V, Iguchi S, Kapur A, Ling W, Mandow R, Mishra S, Reardon D J, Shannon R M, Sudou H, Wang J, Zhao S Y, Zhu X J and Zic A 2026 ApJ Lett. 998 L42
work page 2026
-
[14]
Sesana A 2013 Class. Quant. Grav. 30 224014
work page 2013
-
[15]
Blaes O, Lee M H and Socrates A 2002 ApJ 578 775
work page 2002
-
[16]
Zhu X J, Wen L, Hobbs G, Zhang Y, Wang Y, Madison D R, Manchester R N, Kerr M, Rosado P A and Wang J 2015 MNRAS 449 1650
work page 2015
-
[17]
Taylor S R, Huerta E A, Gair J R and McWilliams S T 2016 ApJ 817 70
work page 2016
-
[18]
Chen J W and Wang Y 2022 ApJ 929 168
work page 2022
-
[19]
Zic A, Reardon D J, Kapur A, Hobbs G, Mandow R, Curylo M, Shannon R M, Askew J, Bailes M, Bhat N D R, Cameron A, Chen Z C, Dai S, Di Marco V, Feng Y, Kerr M, Kulkarni A, Lower M E, Luo R, Manchester R N, Miles M T, Nathan R S, Oslowski S, Rogers A F, Russell C J, Sarkissian J M, Shamohammadi M, Spiewak R, Thyagarajan N, Toomey L, Wang S, Zhang L, Zhang S ...
work page 2023
-
[20]
Reardon D J, Zic A, Shannon R M, Di Marco V, Hobbs G B, Kapur A, Lower M E, Mandow R, Middleton H, Miles M T, Rogers A F, Askew J, Bailes M, Bhat N D R, Cameron A, Kerr M, Kulkarni A, Manchester R N, Nathan R S, Russell C J, Oslowski S and Zhu X J 2023 ApJ Lett. 951 L7
work page 2023
-
[21]
Valtonen M J, Zola S, Ciprini S, Gopakumar A, Mat- sumoto K, Sadakane K, Kidger M, Gazeas K, Nilsson K, Berdyugin A, Piirola V, Jermak H, Baliyan K S, Alicavus F, Boyd D, Campas Torrent M, Campos F, Carrillo Gómez J, Caton D B, Chavushyan V, Dalessio J, Debski B, Dim- itrov D, Drozdz M, Er H, Erdem A, Escartin Pérez A, Fallah Ramazani V, Filippenko A V, G...
work page 2016
-
[22]
Dey L, Valtonen M J, Gopakumar A, Zola S, Hudec R, Pihajoki P, Ciprini S, Matsumoto K, Sadakane K, Kidger M, Nilsson K, Mikkola S, Sillanpää A, Takalo L O, Lehto H J, Berdyugin A, Piirola V, Jermak H, Baliyan K S, Pursimo T, Caton D B, Alicavus F, Baransky A, Blay P, Boumis P, Boyd D, Campas Torrent M, Campos F, Carrillo Gómez J, Chandra S, Chavushyan V, ...
work page 2026
-
[23]
Susobhanan A, Gopakumar A and Taylor S R 2020 Phys. Rev. D 101 043022
work page 2020
-
[24]
Susobhanan A 2023 Class. Quant. Grav. 40 155014
work page 2023
-
[25]
Sillanpää A, Haarala S, Valtonen M J, Sundelius B and Byrd G G 1988 ApJ 325 628
work page 1988
-
[26]
Taylor S R, Baker P T, Hazboun J S, Simon J and Vige- land S J 2021 enterprise_extensions v2.4.3
work page 2021
-
[27]
Ellis J and van Haasteren R 2017 PTMCMCSampler ascl:1707.017
work page 2017
-
[28]
Vousden W D, Farr W M and Mandel I 2016 MNRAS 455 1919
work page 2016
-
[29]
Dickey J M 1971 Ann. Math. Stat. 42 204
work page 1971
-
[30]
Verdinelli I and Wasserman L 1995 J. Am. Stat. Assoc. 90 614
work page 1995
-
[31]
Wagenmakers E J, Lodewyckx T, Kuriyal H and Grasman R 2010 Cogn. Psychol. 60 158
work page 2010
-
[32]
Kushwaha P 2020 Galaxies 8 15
work page 2020
-
[33]
Mei S, Blakeslee J P, Côté P, Tonry J L, West M J, Fer- rarese L, Jordán A, Peng E W, Anthony A and Merritt D 2007 ApJ 655 144
work page 2007
-
[34]
Blakeslee J P, Jordán A, Mei S, Côté P, Ferrarese L, In- fante L, Peng E W, Tonry J L and West M J 2009 ApJ 694 556
work page 2009
-
[35]
Mutabazi T, Blyth S L, Woudt P A, Lucey J R, Jarrett T H, Bilicki M, Schröder A C and Moore S A W 2014 MNRAS 439 3666
work page 2014
-
[36]
Vaduvescu O, Petropoulou V, Reverte D and Pinter V 2018 A&A 616 A165
work page 2018
-
[37]
Agazie G, Arzoumanian Z, Baker P T, Bécsy B, Blecha L, Blumer H, Brazier A, Brook P R, Burke-Spolaor S, Casey-Clyde J A, Charisi M, Chatterjee S, Cheeseboro B D, Cohen T, Cordes J M, Cornish N J, Crawford F, Cro- martie H T, Decesar M E, Demorest P B, Dey L, Dolch T, Ellis J A, Ferdman R D, Ferrara E C, Fiore W, Fonseca E, Freedman G E, Garver-Daniels N, ...
work page 2024
-
[38]
Blanchet L 2014 Living Rev. Relativ. 17 2
work page 2014
-
[39]
Peters P C 1964 Phys. Rev. 136 B1224
work page 1964
-
[40]
Harms R J, Ford H C, Tsvetanov Z I, Hartig G F, Dressel L L, Kriss G A, Bohlin R and Davidsen A F 1994 ApJ Lett. 435 L35
work page 1994
-
[41]
Safarzadeh M, Loeb A and Reid M 2019 MNRAS 488 L90
work page 2019
-
[42]
Wagner S, Richtler T and Hopp U 1991 A&A 241 399
work page 1991
-
[43]
Sanders J S, Fabian A C, Sun M, Churazov E, Simionescu A, Walker S A and Werner N 2014 MNRAS 439 1182
work page 2014
-
[44]
Huang Z and Sarazin C L 1996 ApJ 461 622
work page 1996
-
[45]
Woudt P A, Kraan-Korteweg R C, Lucey J, Fairall A P and Moore S A W 2008 MNRAS 383 445
work page 2008
-
[46]
Agarwal N, Agazie G, Anumarlapudi A, Archibald A M, Arzoumanian Z, Baier J G, Baker P T, Bécsy B, Blecha L, Brazier A, Brook P R, Burke-Spolaor S, Burnette R, Case R, Casey-Clyde J A, Chang Y T, Charisi M, Chatterjee S, Cohen T, Coppi P, Cordes J M, Cornish N J, Crawford F, Cromartie H T, Crowter K, Decesar M E, Demorest P B, Deng H, Dey L, Dolch T, D’Ora...
work page 2026
-
[47]
Abdurro’uf, Accetta K, Aerts C, Silva Aguirre V, Ahu- mada R, Ajgaonkar N, Filiz Ak N, Alam S, Allende Prieto C, Almeida A, Anders F, Anderson S F, Andrews B H, An- guiano B, Aquino-Ortíz E, Aragón-Salamanca A, Argudo- Fernández M, Ata M, Aubert M, A vila-Reese V, Badenes C, Barbá R H, Barger K, Barrera-Ballesteros J K, Beaton R L, Beers T C, Belfiore F, ...
work page 2026
-
[48]
McConnell N J and Ma C P 2013 ApJ 764 184 xxxxxx-12 Chinese Physics Letters xx, xxxxxx (2026) New Submission Appendix A: Catalog of galactic central black hole mass estimates Table 3 compiles the catalog of supermassive black hole masses assumed for member galaxies across the five targeted clusters: Virgo, Fornax, Coma, Hercules, and Norma. These mass est...
work page 2013
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.