Pith. sign in

REVIEW 3 major objections 5 minor 2 cited by

Constraining the Binarity of Massive Black Holes in the Galactic Center and Some Nearby Galaxies via Pulsar Timing Array Observations of Gravitational Waves

T0 review · 3 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash

Pith's one-line read Future pulsar timing array observations can independently reveal or exclude intermediate-mass black hole companions orbiting the massive black holes in the Galactic Center, M31, and M87, and, if pulsars are found near those nuclei, could…

desk verdict A careful, internally consistent sensitivity forecast for PTA detection of MBH companions in nearby nuclei; the headline mass-ratio thresholds are white-noise ideals that red spin noise will shift upward, but the qualitative rankings survive. read the letter →

arxiv 2411.14150 v1 pith:ZXW7Z3BM submitted 2024-11-21 astro-ph.HE astro-ph.GA

classification astro-ph.HEastro-ph.GA
keywords pulsartimingarraysgravitationalwavesmassiveblackholebinariesGalacticCenterintermediate-massholesSagittariusA*eccentricSKA-PTA
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper asks a concrete question: if the massive black hole at the Galactic Center has a lighter black hole companion on a months-to-years orbit, can pulsar timing arrays (PTAs) detect its gravitational waves, and what can a non-detection prove? It derives the expected signal-to-noise ratio for circular and eccentric binaries in the nanohertz band and maps the parameter space that a future Square Kilometer Array PTA (SKA-PTA) could probe in the Galactic Center and in the LMC, M31, M32, and M87. The central quantitative claims are that a 20-year SKA-PTA campaign with 1,000 pulsars at 10 ns timing precision could reveal or rule out a companion of roughly 500–5,000 solar masses (mass ratio $q\sim10^{-4}$–$10^{-3}$) around Sgr A* at separations of 20–3,000 AU, a companion with $q\gtrsim10^{-4}$ in M31, and one with $q\gtrsim10^{-5}$ in M87. The LMC and M32 would not be reachable with Earth-based pulsars for small mass ratios, but if a few millisecond pulsars are found within a parsec of the central black hole, a center-based PTA could push the detectable companion mass down to about 100 solar masses. The paper thus positions PTAs as an independent observational route to the binarity of nearby massive black holes, complementing stellar-orbital and dynamical constraints.

What carries the argument

The machinery that carries the argument is an effective characteristic strain $h_{c,\mathrm{eff}}$ defined at the peak harmonic frequency of an eccentric binary, $f_{\mathrm{pk}}=n_{\mathrm{pk}}f_{\mathrm{orb}}$ with $n_{\mathrm{pk}}=2(1+e)^{1.1954}/(1-e^2)^{1.5}$, together with two signal-to-noise estimators (matched-filtering and cross-correlation) that sum over all gravitational-wave harmonics inside the PTA frequency band. This definition converts any eccentric binary into a point on the same plot as a PTA sensitivity curve, letting the authors draw S/N contours in the mass-ratio versus semimajor-axis plane. A second key element is the near-field geometric factor $\chi$ that correctly handles the case where pulsars sit much closer to the gravitational-wave source than the Earth does; this factor is what gives center-based PTAs their large sensitivity boost over Earth-based arrays.

What would settle it

Observe Sgr A* with a 20-year, 1,000-pulsar SKA-PTA at 10 ns timing precision: a null search for a monochromatic or harmonic gravitational-wave signal would directly test the claimed exclusion region ($q\sim10^{-4}$–$10^{-3}$, $a\sim20$–$3{,}000$ AU), while a deep radio survey that finds far fewer than 5–10 pulsars within 1 pc of Sgr A* would falsify the GC-PTA sensitivity projections.

Watch

Extended reading notes

Core claim

The paper's central claim is that pulsar timing arrays, particularly the future SKA-PTA, can provide an independent way to reveal or exclude low-mass-ratio massive black hole binaries in the Galactic Center, M31, and M87. For the Galactic Center, a non-detection after 20 years of SKA-PTA observations would suggest the absence of an intermediate-mass black hole with mass $\gtrsim 500$–$5{,}000\,M_\odot$ (or $q\sim10^{-4}$–$10^{-3}$) on a semimajor axis of $\sim20$–$3{,}000$ AU; for M31 the reach is $q\gtrsim10^{-4}$ at $a\sim10^2$–$10^4$ AU, and for M87 it is $q\gtrsim10^{-5}$ at $a\sim10^3$–$2\times10^4$ AU. The LMC and M32 are expected to stay out of reach for the SKA-PTA when $q\ll1$, but hypothetical PTAs built from 5–10 millisecond pulsars within 0.1–1 pc of the central black hole would reveal companions with masses down to about $100\,M_\odot$, close to the stellar-mass regime. The paper also shows that high orbital eccentricity does not prevent detection, because the gravitational-wave power is distributed into high harmonics that can fall inside the PTA band.

Load-bearing premise

The most spectacular results assume future surveys find and stably time 5–10 millisecond pulsars within about 0.1–1 parsec of the central black hole in each galaxy, a population that has not yet been confirmed beyond a single magnetar near Sgr A*.

Editorial extensions

If this is right

  • A 20-year SKA-PTA non-detection toward the Galactic Center would independently rule out an intermediate-mass black hole of roughly 500–5,000 solar masses orbiting Sgr A* at 20–3,000 AU.
  • For M31 the same campaign would reveal or exclude companions with $q\gtrsim10^{-4}$ at $10^2$–$10^4$ AU, and for M87 companions with $q\gtrsim10^{-5}$ at $10^3$–$2\times10^4$ AU, within about 20 years.
  • If several millisecond pulsars are discovered within 0.1–1 pc of the central black hole in the GC, LMC, M31, or M32, a center-based PTA could detect companions down to roughly 100 solar masses, close to stellar masses.
  • Highly eccentric binaries remain detectable because their gravitational-wave power moves into high harmonics that can fall inside the PTA band, extending the reachable semimajor axis.
  • The effective characteristic strain formalism gives future searches a common yardstick for reporting detection or exclusion limits on the mass-ratio versus semimajor-axis plane.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If these sensitivity estimates hold, a null SKA-PTA result would become an independent, dynamics-free constraint on intermediate-mass black holes in the Galactic Center, complementing limits from stellar orbits and proper-motion measurements.
  • The same near-field enhancement that powers the GC-PTA idea implies that any future discovery of pulsars within a parsec of a galactic nucleus would make that nucleus a high-value gravitational-wave laboratory, even before a planet-scale array improves.
  • Because the paper deliberately omits red intrinsic spin noise, adding that noise to the model would raise the effective noise floor and shift the quoted mass-ratio boundaries upward; the 100-solar-mass reach should therefore be read as an optimistic ceiling.
  • The formalism is general enough to be turned on stellar-mass bodies: the same center-based method that reaches ~100 solar masses could eventually search for gravitational waves from stellar remnants or even S-stars orbiting Sgr A*, a direction the paper only notes in passing.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper develops a signal-to-noise ratio (S/N) formalism for pulsar timing array (PTA) searches for gravitational waves from an eccentric massive black hole binary (BBH) with a low-mass secondary, and applies it to hypothetical binaries in the Galactic center (GC), LMC, M31, M32, and M87. Using standard Peters-type harmonic decomposition, matched-filter and cross-correlation estimators, and assumed PTA configurations (IPTA, CPTA, SKA-PTA, and hypothetical center-based PTAs), the authors derive sensitivity contours in the mass-ratio versus semimajor-axis plane. Their central quantitative claims are that a future SKA-PTA with 20-year observations could reveal or constrain a GC companion with q ~ 1e-4 to 1e-3 (M_BH,2 ~ 500 to 5000 solar masses), a M31 companion with q ~ 1e-4, and a M87 companion with q ~ 1e-5, while center-based PTAs of a few MSPs within ~0.1-1 pc of the central black holes could push detectable companion masses down to about 100 solar masses.

Significance. If the assumptions hold, the paper provides a useful forecasting framework for an independent probe of low-mass-ratio massive black hole binaries in nearby galactic nuclei. The formalism is internally consistent: the eccentric-orbit harmonic decomposition follows standard Peters/Maggiore expressions, the matched-filter and cross-correlation S/N estimates are explicit, and the comparison with existing dynamical limits in Figure 5 gives context for the claimed constraints. The paper also credits and uses prior work for the near-field geometric factor and the GWB model, and it clearly labels the center-based PTA scenarios as conditional on future discoveries. The main value is the identification of parameter space where future PTAs could complement stellar-dynamics and S-star constraints on a hypothetical IMBH around Sgr A*, and similar bounds for M31 and M87.

major comments (3)
  1. [Sections 4.1 and 5.1, Eqs. (20)-(26)] The noise model is white shot noise plus the GWB only; red intrinsic spin noise is explicitly omitted as 'quite uncertain' (Section 4.1, Eqs. (20)-(22)). At the 1-30 nHz frequencies that set the quoted q ~ 1e-4 to 1e-3 thresholds, measured red spin noise in current PTA datasets is often comparable to or larger than a 10 ns white-timing floor. Because h_c,i is proportional to q for q << 1 and h_n appears in the denominator of Eqs. (23)-(26), including red noise raises the minimum detectable q by approximately h_n,red/h_n,white for the matched-filter method and by the square root of that ratio for the cross-correlation method. The headline statement in Section 5.1 that a 20-year SKA-PTA non-detection 'would suggest independently that there is no IMBH with mass greater than about 500 to 5,000 Msun' is therefore an idealized white-noise sensitivity bound, not a realistic non-detection forecast. This is not an internal inconsistency, but the quoted mass boundaries should be presented as optimistic limiting sensitivities and preferably supplemented with a representative red-noise case.
  2. [Table 1, Sections 5.2, 6.1-6.3] The center-based PTA scenarios (GC-PTA, LMCC-PTA, M31C-PTA, M32C-PTA) assume the existence of 5-10 timing-stable MSPs within 0.1-1 pc of each central MBH (Table 1, column r_p-BBH). This assumption is the decisive factor behind the headline result that companions down to about 100 solar masses could be revealed. For the GC only a single magnetar at roughly 0.1 pc is currently known; for M31, M32, and the LMC no such pulsars are confirmed. The paper's abstract does hedge with 'if a number of milli-second stable pulsars ... can be detected in future', but Sections 5.2 and 6.2-6.3 present the 100-solar-mass reach as a main result without a quantitative assessment of the probability or feasibility of discovering and timing such pulsars in the dense nuclear environment. These claims should be explicitly separated as a highly speculative projection, or accompanied by a discussion of scattering, dispersion measure variations, and timing stability requirements.
  3. [Section 4.4 and Figures 3-9] The interpretation of the S/N contours as exclusion or discovery boundaries uses a single deterministic threshold (rho_th = 1) with no accounting for detection probability, false-alarm rate, or noise realization. For example, the text after Figure 3 states that 'the parameter space below the brown curves cannot be ruled out' and that a non-detection 'would suggest' no IMBH above the quoted mass. A sensitivity curve with rho_th = 1 is a reasonable projection tool, but converting it into a non-detection statement requires a statistical framework (e.g., detection probability as a function of the source parameters, or a Bayesian upper limit). As written, the exclusion language overstates what a single S/N threshold actually delivers. I recommend rephrasing these conclusions as sensitivity boundaries or adding a brief statistical treatment of the non-detection case.
minor comments (5)
  1. [Page 10, text after Figure 5] There is a typo: 'can also be imited by some dynamical arguments' should read 'can also be limited by some dynamical arguments'.
  2. [Section 6.1, footnote 6] The footnote says 'the black hole in the LMC is in the mass rage of IMBHs'; 'rage' should be 'range'.
  3. [Section 6.2, first paragraph] The sentence 'We hypothesize that a BBH exists in the M31 center, with total mass M_BBH = 1.4 x 10^8 Msun and an mass ratio q' contains a grammatical error: 'an mass ratio' should be 'a mass ratio'.
  4. [Figures 3-13] The brown curves (GWB-removed cases) are described in the captions but not labeled directly in the figure panels; adding a legend entry or a label to the brown curves would improve readability.
  5. [Section 5.2, Eq. (33)] The remark that H_i in Eqs. (24) and (26) is 'the GW strain at Earth with r = 8 kpc' while Eq. (9) already includes this distance is a potential source of confusion; a short clarifying sentence on how the near-field chi modifies the S/N would help.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the quoted PTA thresholds follow from standard GW strain and PTA-sensitivity formulas with explicitly stated assumptions; the self-citations are to published derivations/inputs, not to the claimed result.

full rationale

The paper's derivation chain is self-contained: Eqs. (1)-(10) give the GW strain from standard quadrupole/Peters formulas; Eqs. (14)-(17) define the characteristic strain; Eqs. (18)-(22) define the noise from white timing residuals and the GWB; Eqs. (23)-(26) give the S/N for matched-filter and cross-correlation methods. The Section 5.1 thresholds (q ~ 1e-4 to 1e-3, M_BH,2 ~ 500-5000 Msun) are obtained by setting S/N=1 and inverting h_c proportional to q, not by fitting any parameter to the target claim. The PTA parameters in Table 1 are stated assumptions, and the paper explicitly flags the omission of red spin noise and the speculative nature of center pulsars; these are limitations on realism, not circularity. Citations to Guo et al. (2022) for the geometric factor chi and to Chen et al. (2020, 2023) for the GWB model are published derivations/measured inputs that do not incorporate the paper's conclusion, so they do not make the argument circular. There is no fitted input renamed as a prediction and no self-citation chain forcing the result.

Assumptions & free parameters 7 free parameters · 8 assumptions · 0 invented entities

The paper is a sensitivity forecast built from standard GW formulas and assumed instrument parameters. The free parameters are the chosen PTA timing precision, cadence, observation span, pulsar counts, pulsar distances, GWB noise shape, and detection threshold. No new physical entities are introduced: the hypothetical companions are existing astrophysical hypotheses, and the near-center pulsar PTAs are observational scenarios, not new physics.

free parameters (7)
  • PTA timing residual sigma_t = IPTA 100 ns; CPTA 20 ns; SKA-PTA 10 ns; hypothetical center PTAs 100 ns
    Chosen timing precision of pulsars (Table 1); directly sets the shot-noise floor and the sensitivity curves.
  • Number of MSPs N_pl = 49 (IPTA), 100 (CPTA), 1000 (SKA-PTA), 10 (GC-PTA), 5 (LMCC/M31C/M32C PTAs)
    Assumed pulsar counts; sensitivity scales roughly as sqrt(N_pl) in matched filtering and as N_pl^(1/4) to N_pl^(1/2) across methods.
  • Observation cadence Delta_t = 0.04 yr for IPTA/CPTA/SKA-PTA; 0.02 yr for hypothetical center PTAs
    Sets the high-frequency cutoff of the PTA band; chosen in Table 1.
  • Observation span T_obs = 20 yr for IPTA/CPTA/SKA-PTA/GC-PTA; 10 yr for LMCC/M31C/M32C PTAs
    Sets the low-frequency cutoff and the number of cycles; chosen in Table 1.
  • Pulsar distance to the source r_pl = 1 pc for GC/M31/M32; 0.1 pc for LMC
    Drives the near-field geometric factor chi, which scales as r/r_pl and produces the large S/N boost for center PTAs.
  • GWB amplitude and spectral shape parameters = A = 2.5e-15, f_bend = 1.15e-10 Hz, kappa_gw = 3.70, gamma_gw = 0.19
    Adopted from Chen et al. 2020 and 2023, from the same research group, to model GWB confusion noise; affects low-frequency sensitivity.
  • S/N threshold rho_th = 1, with contours also at 3, 10, and 100
    Detection threshold used to draw sensitivity curves; changing it shifts the quoted parameter-space boundaries.
assumptions (8)
  • standard math Standard gravitational-wave emission formulas for eccentric binaries, including Peters 1964 and Peters-Mathews harmonic decomposition with Bessel functions.
    Used throughout Sections 2 and 3 to compute strain amplitudes and harmonic power.
  • domain assumption Newtonian Keplerian orbits with slow GW-driven inspiral, so the orbital phase follows Keplerian motion and the binary is treated as stationary over the observation.
    Stated after Equation (4): the GW decay timescale is much longer than the orbital period for all systems considered.
  • domain assumption PTA noise consists only of white timing noise plus GWB confusion; red intrinsic spin noise is ignored.
    Section 4.1 states that red spin noise is uncertain and is not considered, which can make sensitivity forecasts optimistic.
  • domain assumption The geometric factor chi for the pulsar-Earth-source configuration follows Guo et al. 2022, including near-field enhancement when pulsars sit close to the GW source.
    Used in Sections 4.2 and 5.2, Equation (33); central to the large S/N boost for center-based PTAs.
  • domain assumption The GWB signal can be modeled and subtracted from PTA data in the optimistic scenarios.
    Section 4.1 and Figures 3 to 4: if the GWB can be well detected, it is treated as a known signal and removed; this has not yet been demonstrated at the required precision.
  • domain assumption Adopted central black hole masses and distances for the GC, LMC, M31, M32, and M87 from cited observations.
    Sections 5 and 6 use MBBH = 4.26e6, 2.4e4, 1.4e8, 2.4e6, and 6.5e9 solar masses respectively; these set the strain amplitudes.
  • ad hoc to paper Hypothetical future pulsar populations within about 0.1 to 1 pc of the central MBHs exist and can be discovered and timed.
    Defines GC-PTA, LMCC-PTA, M31C-PTA, and M32C-PTA in Table 1 and Sections 5.2 and 6.1 to 6.3; this is not yet established observationally.
  • domain assumption A fixed S/N threshold of 1 (or 3, 10, 100) marks a detection, without false-alarm statistics.
    Used to draw sensitivity curves and quoted constraints; practical searches typically require a higher threshold. Location: Section 4.4 and figures.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Constraining the Binarity of Massive Black Holes in the Galactic Center and Some Nearby Galaxies via Pulsar Timing Array Observations of Gravitational Waves." pith.science (2026). https://pith.science/paper/ZXW7Z3BM

@misc{pith2026241114150,
  author       = {Pith},
  title        = {Pith review of: Constraining the Binarity of Massive Black Holes in the Galactic Center and Some Nearby Galaxies via Pulsar Timing Array Observations of Gravitational Waves},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ZXW7Z3BM}},
  note         = {Machine review of arXiv:2411.14150}
}
abstract

Massive black holes (MBHs) exist in the Galactic center (GC) and other nearby galactic nuclei. As natural outcome of galaxy mergers, some MBHs may have a black hole (BH) companion. In this paper, assuming that the MBHs in the GC and some nearby galaxies are in binaries with orbital periods ranging from months to years (gravitational-wave frequency $\sim1-100$\,nHz), we investigate the detectability of gravitational-waves from these binary MBHs (BBHs) and constraints on the parameter space for the existence of BBHs in the GC, LMC, M31, M32, and M87, that may be obtained by current/future pulsar timing array (PTA) observations. We find that a BBH in the GC, if any, can be revealed by the Square Kilometer Array PTA (SKA-PTA) if its mass ratio $q\gtrsim10^{-4}-10^{-3}$ and semimajor axis $a\sim20-10^3$\,AU. The existence of a BH companion of the MBH can be revealed by SKA-PTA with $\sim20$-year observations in M31 if $q\gtrsim10^{-4}$ and $a\sim10^2-10^4$\,AU or in M87 if $q\gtrsim10^{-5}$ and $a\sim10^3-2\times10^4$\,AU, but not in LMC and M32 if $q\ll1$. If a number of milli-second stable pulsars with distances $\lesssim0.1-1$\,pc away from the central MBH in the GC, LMC, M32, or M31, can be detected in future and applied to PTAs, the BH companion with mass even down to $\sim100M_\odot$, close to stellar masses, can be revealed by such PTAs. Future PTAs are expected to provide an independent way to reveal BBHs and low-mass MBH companions in the GC and nearby galaxies, improving our understandings of the formation and evolution of MBHs and galaxies.

Figures

Figures reproduced from arXiv: 2411.14150 by the authors.

Figure 1
Figure 1. shows some examples for the GW waveform emitted from a hypothetical BBH in the GC, with q = 0.01, a = 100AU, e = 0 (top panel), 0.5 (middle panel), or 0.9 (bottom panel). As shown in this figure, with increasing eccentricity e, the amplitude of the GW strain emitted at the pericenter increases, but the time duration for the peak GW radiation amplitude becomes shorter. This suggests that the GW signal from a highly e… view at source ↗
Figure 2
Figure 2. Sensitivity curves (hsc) of different PTAs and the effective characteristic strains (hc,eff ) for GW signals radiated from a hypothetical GC BBH with different eccentricities. Left and right panels represent the results obtained by the matched-filtering (MF) and the cross-correlation (CC) method, respectively. The two vertical black dotted lines in each panel indicate the upper limits and the lower limits for the PT… view at source ↗
Figure 3
Figure 3. Contours of the expected S/N (ϱ) in the q/(1 + q) versus a parameter space of the hypothetical BBH in the GC, monitored by the assumed future SKA-PTA (see Tab. 1). The hypothetical BBH is assumed to have a mass ratio q, semimajor axis a, and eccentricity e. Top-left, top-right, bottom-left, and bottom-right panels show the cases with e = 0.2, 0.5, 0.7, and 0.9, respectively. For each panel, the left and the right ve… view at source ↗
Figures from the paper (10 more)
Figure 4
Figure 4. Figure 4: Expected S/N of a hypothetical BBH in the GC with mass ratio q, semimajor axis a, and eccentricity e, monitored by an assumed PTA composed of MSPs in the GC (see GC-PTA in Tab. 1). Legends are similar to those for [PITH_FULL_IMAGE:figures/full_fig_p009_4.png]
Figure 5
Figure 5. Figure 5: Possible constraints on the parameter space for the existence of a BBH in the GC expected from future PTA observations. The blue, red, and brown lines represent those BBHs (e = 0.2) with S/N ϱ = 1 estimated for the CPTA, SKA-PTA, and GC-PTA, respectively, with paramete…
Figure 6
Figure 6. Figure 6: Expected S/N of a hypothetical BBH in the LMC center with mass ratio q, semimajor axis a, and eccentricity e monitored by an assumed PTA composed of MSPs in the LMC center (see LMCC-PTA in Tab. 1). Legends are similar to those for [PITH_FULL_IMAGE:figures/full_fig_p01…
Figure 7
Figure 7. Figure 7: Expected S/N of a hypothetical BBH in the center of M31 with mass ratio q, semimajor axis a, and eccentricity e, monitored by the IPTA (see Tab. 1). Legends are similar to those for [PITH_FULL_IMAGE:figures/full_fig_p012_7.png]
Figure 8
Figure 8. Figure 8: Expected S/N of a hypothetical BBH in the center of M31 with mass ratio q, semimajor axis a, and eccentricity e, monitored by the CPTA (see Tab. 1). Legends are similar to those for [PITH_FULL_IMAGE:figures/full_fig_p013_8.png]
Figure 9
Figure 9. Figure 9: Expected S/N of a hypothetical BBH in the center of M31 with mass ratio q, semimajor axis a, and eccentricity e, monitored by the SKA-PTA (see Tab. 1). Legends are similar to those for [PITH_FULL_IMAGE:figures/full_fig_p013_9.png]
Figure 10
Figure 10. Figure 10: Expected S/N of a hypothetical BBH in the M31 center with mass ratio q, semimajor axis a, and eccentricity e, monitored by an assumed PTA composed of MSPs in the M31 center (see M31C–PTA in Tab. 1). Legends are similar to those for [PITH_FULL_IMAGE:figures/full_fig_p…
Figure 11
Figure 11. Figure 11: Expected S/N of a hypothetical equal-mass BBH in the center of M32 as a function of the orbital semimajor axis a, monitored by the SKA-PTA (see Tab. 1). Red, green, and cyan solid (or dashed) lines show the results obtained for a BBH with eccentricity e = 0, 0.5, and …
Figure 12
Figure 12. Figure 12: Expected S/N of a hypothetical BBH with mass ratio q, and semimajor axis a in the center of M32, monitored by a PTA with MSPs in the center of M32 (see M32C-PTA in Tab. 1). The legend is similar to that for [PITH_FULL_IMAGE:figures/full_fig_p016_12.png]
Figure 13
Figure 13. Figure 13: Expected S/N of a hypothetical BBH in the center of M87 with mass ratio q, semimajor axis a, and eccentricity e, monitored by the SKA-PTA. The legend is similar to that for [PITH_FULL_IMAGE:figures/full_fig_p016_13.png]

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Probing an Intermediate-Mass Black Hole Companion of Sagittarius A* with Pulsar Timing

    astro-ph.HE 2026-07 conditional novelty 6.0 of 10

    A 1PN numerical timing model shows that a pulsar orbiting Sgr A* produces large, distinctive post-fit residuals from an IMBH companion, enabling constraints that fill gaps left by existing S-star and proper-motion bounds.

  2. Prospects for Time-Domain and Multi-Messenger Science with eXTP

    astro-ph.HE 2025-06 conditional novelty 4.0 of 10

    A mission white paper forecasting the time-domain and multi-messenger science of the future X-ray observatory eXTP, with detection rates and model-discrimination estimates.

Reference graph

Works this paper leans on

88 extracted references · 61 canonical work pages · cited by 2 Pith papers

  1. [1]

    2005, Physics Reports, 419, 65

    Alexander, T. 2005, Physics Reports, 419, 65

  2. [2]

    2022, MNRAS, 510, 4873

    Antoniadis, J., Arzoumanian, Z., Babak, S., et al. 2022, MNRAS, 510, 4873

  3. [3]

    2023a, arXiv e-prints, arXiv:2306.16227 —

    Antoniadis, J., Arumugam, P., Arumugam, S., et al. 2023a, arXiv e-prints, arXiv:2306.16227 —. 2023b, arXiv e-prints, arXiv:2306.16226

  4. [4]

    T., Blumer, H., et al

    Arzoumanian, Z., Baker, P. T., Blumer, H., et al. 2020, ApJL, 905, L34

  5. [5]

    Balick, B., & Brown, R. L. 1974, ApJ, 194, 265

  6. [6]

    C., Blandford, R

    Begelman, M. C., Blandford, R. D., & Rees, M. J. 1980, Nature, 287, 307

  7. [7]

    2005, ApJ, 631, 280

    Bender, R., Kormendy, J., Bower, G., et al. 2005, ApJ, 631, 280

  8. [8]

    2015, Science China Physics, Mechanics, and Astronomy, 58, 5748

    Blair, D., Ju, L., Zhao, C., et al. 2015, Science China Physics, Mechanics, and Astronomy, 58, 5748

Show all 88 references
  1. [9]

    M., Sch¨odel, R., et al

    Boehle, A., Ghez, A. M., Sch¨odel, R., et al. 2016, ApJ, 830, 17

  2. [10]

    P., et al

    Boyce, H., L¨utzgendorf, N., van der Marel, R. P., et al. 2017, ApJ, 846, 14

  3. [11]

    2016, Monthly Notices of the Royal Astronomical Society, 458, 1267

    Brazier, A., Lassus, A., Petiteau, A., et al. 2016, Monthly Notices of the Royal Astronomical Society, 458, 1267

  4. [12]

    2018, Commun

    Cai, R.-G., Liu, T.-B., & Wang, S.-J. 2018, Commun. Theor. Phys., 70, 735

  5. [13]

    Chen, S., Sesana, A., & Conselice, C. J. 2019, MNRAS, 488, 401

  6. [14]

    2020, ApJ, 897, 86 —

    Chen, Y ., Yu, Q., & Lu, Y . 2020, ApJ, 897, 86 —. 2023, ApJ, 955, 132 —. 2024, ApJ, 974, 261

  7. [15]

    M., & Chatterjee, S

    Cordes, J. M., & Chatterjee, S. 2019, ARA&A, 57, 417

  8. [16]

    2011, Gravitational-Wave Physics and Astronomy: An Introduction to Theory, Experiment and Data Analysis

    Creighton, J., & Anderson, W. 2011, Gravitational-Wave Physics and Astronomy: An Introduction to Theory, Experiment and Data Analysis. (Wiley -VCH Verlag GmbH & Co. KGaA)

  9. [17]

    1979, ApJ, 234, 1100

    Detweiler, S. 1979, ApJ, 234, 1100

  10. [18]

    2019, Science, 365, 664

    Do, T., Hees, A., Ghez, A., et al. 2019, Science, 365, 664

  11. [19]

    W., & Antonini, F

    Erkal, D., Boubert, D., Gualandris, A., Evans, N. W., & Antonini, F. 2019, MNRAS, 483, 2007

  12. [20]

    A., Rasskazov, A., Remmelzwaal, A., et al

    Evans, F. A., Rasskazov, A., Remmelzwaal, A., et al. 2023, MNRAS, 525, 561 Event Horizon Telescope Collaboration, Akiyama, K., Alberdi, A., et al. 2019, ApJL, 875, L1

  13. [21]

    J., & Zimmerman, A

    Fouvry, J.-B., Bustamante-Rosell, M. J., & Zimmerman, A. 2023, MNRAS, 526, 1471

  14. [22]

    2011, ApJ, 729, 119

    Gebhardt, K., Adams, J., Richstone, D., et al. 2011, ApJ, 729, 119

  15. [23]

    2010, Rev

    Genzel, R., Eisenhauer, F., & Gillessen, S. 2010, Rev. Mod. Phys., 82, 3121

  16. [24]

    M., Duchˆene, G., Matthews, K., et al

    Ghez, A. M., Duchˆene, G., Matthews, K., et al. 2003, ApJL, 586, L127

  17. [25]

    M., Salim, S., Weinberg, N

    Ghez, A. M., Salim, S., Weinberg, N. N., et al. 2008, ApJ, 689, 1044

  18. [26]

    2009, ApJ, 692, 1075

    Gillessen, S., Eisenhauer, F., Trippe, S., et al. 2009, ApJ, 692, 1075

  19. [27]

    M., Eisenhauer, F., et al

    Gillessen, S., Plewa, P. M., Eisenhauer, F., et al. 2017, ApJ, 837, 30

  20. [28]

    2019, MNRAS, 482, 3669

    Girma, E., & Loeb, A. 2019, MNRAS, 482, 3669

  21. [29]

    M., Sesana, A., Holgado, A

    Goldstein, J. M., Sesana, A., Holgado, A. M., & Veitch, J. 2019, MNRAS, 485, 248 CONSTRAINING BINARITY OF NEARBY MBH S BY PTA O BSERVATIONS 19

  22. [30]

    2020, MNRAS, 497, 3264 GRA VITY Collaboration, Abuter, R., Amorim, A., et al

    Goncharov, B., Zhu, X.-J., & Thrane, E. 2020, MNRAS, 497, 3264 GRA VITY Collaboration, Abuter, R., Amorim, A., et al. 2020, A&A, 636, L5 Gravity Collaboration, Straub, O., Baub¨ock, M., et al. 2023, A&A, 672, A63

  23. [31]

    2010, MNRAS, 409, 1146

    Gualandris, A., Gillessen, S., & Merritt, D. 2010, MNRAS, 409, 1146

  24. [32]

    2009, ApJ, 705, 361

    Gualandris, A., & Merritt, D. 2009, ApJ, 705, 361

  25. [33]

    2007, MNRAS, 376, L29

    Gualandris, A., & Portegies Zwart, S. 2007, MNRAS, 376, L29

  26. [34]

    2022, ApJ, 939, 55

    Guo, X., Lu, Y ., & Yu, Q. 2022, ApJ, 939, 55

  27. [35]

    Hansen, B. M. S., & Milosavljevi´c, M. 2003, ApJL, 593, L77

  28. [36]

    A., McWilliams, S

    Huerta, E. A., McWilliams, S. T., Gair, J. R., & Taylor, S. R. 2015, PhRvD, 92, 063010

  29. [37]

    D., Frail, D

    Hyman, S. D., Frail, D. A., Deneva, J. S., et al. 2019, ApJ, 876, 20

  30. [38]

    2009, Analysis of gravitational-wave data, V ol

    Jaranowski, P., & Krolak, A. 2009, Analysis of gravitational-wave data, V ol. 29 (Cambridge: Cambridge Univ. Press), doi:10.1017/CBO9780511605482

  31. [39]

    C., Arumugasamy, P., Bagchi, M., et al

    Joshi, B. C., Arumugasamy, P., Bagchi, M., et al. 2018, Journal of Astrophysics and Astronomy, 39, 51

  32. [40]

    A., Burrows, D

    Kennea, J. A., Burrows, D. N., Kouveliotou, C., et al. 2013, ApJL, 770, L24

  33. [41]

    2012, ApJ, 752, 67

    Kocsis, B., Ray, A., & Portegies Zwart, S. 2012, ApJ, 752, 67

  34. [42]

    Kormendy, J., & Ho, L. C. 2013, ARA&A, 51, 511

  35. [43]

    Kramer, M., & Champion, D. J. 2013, Classical and Quantum Gravity, 30, 224009

  36. [44]

    Lazio, T. J. W. 2013, Classical and Quantum Gravity, 30, 224011

  37. [45]

    Lee, K. J. 2016, in Astronomical Society of the Pacific Conference

  38. [46]

    502, Frontiers in Radio Astronomy and FAST Early Sciences Symposium 2015, ed

    Series, V ol. 502, Frontiers in Radio Astronomy and FAST Early Sciences Symposium 2015, ed. L. Qain & D. Li, 19

  39. [47]

    J., Wex, N., Kramer, M., et al

    Lee, K. J., Wex, N., Kramer, M., et al. 2011, MNRAS, 414, 3251

  40. [48]

    M., Coles, W

    Lentati, L., Shannon, R. M., Coles, W. A., et al. 2016, MNRAS, 458, 2161

  41. [49]

    P., Kanekar, N., Frail, D

    Macquart, J. P., Kanekar, N., Frail, D. A., & Ransom, S. M. 2010, ApJ, 715, 939

  42. [50]

    2008, Gravitational waves: V olume 1: Theory and experiments (Oxford: Oxford University Press)

    Maggiore, M. 2008, Gravitational waves: V olume 1: Theory and experiments (Oxford: Oxford University Press)

  43. [51]

    Manchester, R. N. 2013, International Journal of Modern Physics D, 22, 1341007

  44. [52]

    N., Hobbs, G

    Manchester, R. N., Hobbs, G. B., Teoh, A., & Hobbs, M. 2005, AJ, 129, 1993

  45. [53]

    N., & IPTA

    Manchester, R. N., & IPTA. 2013, Classical and Quantum Gravity, 30, 224010

  46. [54]

    McLaughlin, M. A. 2013, Classical and Quantum Gravity, 30, 224008

  47. [55]

    2009, ApJL, 693, L35

    Merritt, D., Gualandris, A., & Mikkola, S. 2009, ApJL, 693, L35

  48. [56]

    T., Shannon, R

    Miles, M. T., Shannon, R. M., Bailes, M., et al. 2023, MNRAS, 519, 3976

  49. [57]

    Mingarelli, C. M. 2015, Gravitational wave astrophysics with pulsar timing arrays (Springer)

  50. [58]

    2011, International Journal of Modern Physics D, 20, 989

    Nan, R., Li, D., Jin, C., et al. 2011, International Journal of Modern Physics D, 20, 989

  51. [59]

    M., Ramirez-Ruiz, E., et al

    Naoz, S., Will, C. M., Ramirez-Ruiz, E., et al. 2020, ApJL, 888, L8

  52. [60]

    Perera, B. B. P., DeCesar, M. E., Demorest, P. B., et al. 2019, MNRAS, 490, 4666

  53. [61]

    Peters, P. C. 1964, Physical Review, 136, 1224

  54. [62]

    2004, ApJ, 615, 253 Pietrzy´nski, G., Graczyk, D., Gieren, W., et al

    Pfahl, E., & Loeb, A. 2004, ApJ, 615, 253 Pietrzy´nski, G., Graczyk, D., Gieren, W., et al. 2013, Nature, 495, 76 Portegies Zwart, S. F., Baumgardt, H., McMillan, S. L. W., et al. 2006, ApJ, 641, 319

  55. [63]

    2019, in Bulletin of the American Astronomical Society, V ol

    Ransom, S., Brazier, A., Chatterjee, S., et al. 2019, in Bulletin of the American Astronomical Society, V ol. 51, 195

  56. [64]

    J., Zic, A., Shannon, R

    Reardon, D. J., Zic, A., Shannon, R. M., et al. 2023, ApJL, 951, L6

  57. [65]

    J., & Brunthaler, A

    Reid, M. J., & Brunthaler, A. 2004, ApJ, 616, 872 —. 2020, ApJ, 892, 39

  58. [66]

    A., Sesana, A., & Gair, J

    Rosado, P. A., Sesana, A., & Gair, J. 2015, MNRAS, 451, 2417

  59. [67]

    2019, MNRAS, 488, L90

    Safarzadeh, M., Loeb, A., & Reid, M. 2019, MNRAS, 488, L90

  60. [68]

    J., & McWilliams, S

    Sampson, L., Cornish, N. J., & McWilliams, S. T. 2015, PhRvD, 91, 084055

  61. [69]

    Sazhin, M. V . 1978, Soviet Ast., 22, 36 Sch¨odel, R., Ott, T., Genzel, R., et al. 2003, ApJ, 596, 1015 —. 2002, Nature, 419, 694

  62. [70]

    2016, MNRAS, 459, 1737

    Schutz, K., & Ma, C.-P. 2016, MNRAS, 459, 1737

  63. [71]

    2013, Classical and Quantum Gravity, 30, 244009

    Sesana, A. 2013, Classical and Quantum Gravity, 30, 244009

  64. [72]

    2005, ApJ, 623, 23

    Sesana, A., Haardt, F., Madau, P., & V olonteri, M. 2005, ApJ, 623, 23

  65. [73]

    2010, Classical and Quantum Gravity, 27, 084016

    Sesana, A., & Vecchio, A. 2010, Classical and Quantum Gravity, 27, 084016

  66. [74]

    2009, MNRAS, 394, 2255

    Sesana, A., Vecchio, A., & V olonteri, M. 2009, MNRAS, 394, 2255

  67. [75]

    R., Kramer, M., et al

    Smits, R., Lorimer, D. R., Kramer, M., et al. 2009, A&A, 505, 919

  68. [76]

    2019, ApJL, 871, L1

    Takekawa, S., Oka, T., Iwata, Y ., Tsujimoto, S., & Nomura, M. 2019, ApJL, 871, L1

  69. [77]

    T., et al

    Taylor, S., Burke-Spolaor, S., Baker, P. T., et al. 2019, BAAS, 51, 336

  70. [78]

    Taylor, S. R. 2021, arXiv e-prints, arXiv:2105.13270

  71. [79]

    L., Dressler, A., Blakeslee, J

    Tonry, J. L., Dressler, A., Blakeslee, J. P., et al. 2001, ApJ, 546, 681

  72. [80]

    2017, ApJL, 850, L5 van Haasteren, R

    Tsuboi, M., Kitamura, Y ., Tsutsumi, T., et al. 2017, ApJL, 850, L5 van Haasteren, R. 2014, Gravitational Wave Detection and Data Analysis for Pulsar Timing Arrays (Springer)

  73. [81]

    Wang, Y ., & Mohanty, S. D. 2017, Physical Review Letters, 118, 151104

  74. [82]

    M., Naoz, S., Hees, A., et al

    Will, C. M., Naoz, S., Hees, A., et al. 2023, ApJ, 959, 58

  75. [83]

    2023, Research in Astronomy and Astrophysics, 23, 075024 20 G UO, YU, & L U

    Xu, H., Chen, S., Guo, Y ., et al. 2023, Research in Astronomy and Astrophysics, 23, 075024 20 G UO, YU, & L U

  76. [84]

    2002, MNRAS, 331, 935 —

    Yu, Q. 2002, MNRAS, 331, 935 —. 2020, The Innovation, 1, 100063

  77. [85]

    Yu, Q., Lu, Y ., & Lin, D. N. C. 2007, ApJ, 666, 919

  78. [86]

    2003, ApJ, 599, 1129

    Yu, Q., & Tremaine, S. 2003, ApJ, 599, 1129

  79. [87]

    2014, ApJ, 784, 106

    Zhang, F., Lu, Y ., & Yu, Q. 2014, ApJ, 784, 106

  80. [88]

    2014, Monthly Notices of the Royal Astronomical Society, 444, 3709

    Zhu, X.-J., Hobbs, G., Wen, L., et al. 2014, Monthly Notices of the Royal Astronomical Society, 444, 3709

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.