Pith. sign in

REVIEW 3 major objections 5 minor 3 cited by

Nanohertz gravitational waves can measure pulsar distances to sub-parsec precision from just a few sources.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-03 17:02 UTC pith:ZH57IP7M

load-bearing objection A useful 2D extension of pulsar-distance inference, but the headline precision numbers are best-of-84 partner-selection minima, so coverage needs to be shown before trusting them. the 3 major comments →

arxiv 2512.10729 v3 pith:ZH57IP7M submitted 2025-12-11 gr-qc astro-ph.COastro-ph.HEhep-ph

Two-Dimensional Pulsar Distance Inference from Nanohertz Gravitational Waves

Si-Ren Xiao , Ji-Yu Song , Yue Shao , Ling-Feng Wang , Jing-Fei Zhang , Xin Zhang This is my paper
classification gr-qc astro-ph.COastro-ph.HEhep-ph
keywords pulsar timing arrayscontinuous gravitational wavespulsar distancessub-parsec distancetwo-dimensional posteriorsupermassive black hole binariesSKA-era PTABayesian inference
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

This paper shows that the distance to a pulsar—long a limiting uncertainty for pulsar timing arrays—can be measured to sub-parsec precision by combining the pulsar-term phase information from a small number of continuous gravitational-wave sources. The key is to infer distances of two pulsars jointly, building a two-dimensional posterior that preserves the multi-peaked distance structure that one-dimensional marginalization destroys. Different sources imprint different periodic distance ladders, so combining them in two dimensions suppresses spurious peaks and isolates the true distance. In simulated SKA-era arrays, a ~1 kpc pulsar reaches sub-parsec distance uncertainty with only four or five sources at 50 ns timing noise, an order-of-magnitude improvement over existing one-dimensional combination approaches. If this holds on real data, precise pulsar distances will sharpen sky localization of nanohertz sources and help identify their host galaxies.

Core claim

The paper's central claim is that the distance–distance degeneracy shared by pairs of pulsars, which a common continuous gravitational-wave source induces, can be retained and exploited rather than averaged away. The authors construct a two-dimensional joint posterior over two pulsar distances, using the analytic mapping from each source's pulsar-term initial phase to a ladder of distance solutions. Each source yields a set of periodic, roughly parallel bands in the (L_p, L_q) plane; bands from different sources agree only near the true distance pair, so multiplying them narrows the posterior. The method then marginalizes over the partner pulsar to obtain a final one-dimensional distance con

What carries the argument

The load-bearing object is Eq. (4), the analytic pulsar-term initial phase–distance relation, which maps a sampled pulsar-term phase Φ_p to a discrete ladder of distance solutions L_p(Φ_p+2kπ, f0, M, Φ0). The spacing of the ladder is set by the GW frequency and the source–pulsar geometry. The two-dimensional posterior (Eq. 2) multiplies these ladders across all N sources and across a pair of pulsars, so that mismatched spacings destructively interfere while the true distance combination remains consistent. This preserves the distance–distance degeneracy structure that one-dimensional marginalization integrates out, and it is why a handful of sources can suppress the spurious peaks that degra

Load-bearing premise

The method assumes the sky positions of the continuous-wave sources are known exactly before the distance inference; if a source's true position differs from the assumed one, the phase–distance relation is systematically shifted and the inferred pulsar distances become biased.

What would settle it

Take the authors' own simulation and inject one source with a sky position offset by a few degrees (or sample source positions from a realistic PTA localization posterior instead of fixing them). If the sub-parsec fraction of the recovered pulsar distances collapses or the posteriors shift, the external-sky-position requirement is the limiting step. Alternatively, apply the method to a real CGW candidate with a known host galaxy and check whether the recovered pulsar distances agree with parallax distances to sub-parsec accuracy.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • Sub-parsec pulsar distances turn the pulsar term from a nuisance into a phase reference, which should sharply improve the sky localization of nanohertz continuous-gravitational-wave sources.
  • Improved localization makes host-galaxy identification for supermassive-black-hole binaries practical, opening the door to multimessenger follow-up and standard-siren cosmology.
  • The two-dimensional combination requires fewer continuous-wave sources than the one-dimensional method, so it remains usable with the sparse source populations expected in the SKA era.
  • Because the method selects the partner pulsar that gives the tightest constraint, it naturally identifies which pulsar pairs in an array carry the most distance information.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the sky-position assumption degrades in practice—say a source is misidentified or the host-galaxy association is ambiguous—the phase–distance relation shifts and the inferred pulsar distances will be biased; the method's gain is therefore contingent on externally secured source positions, which the paper fixes but does not itself produce.
  • The ladder-suppression mechanism suggests a straightforward generalization: joining three or more pulsars in a higher-dimensional posterior should suppress spurious peaks even further, at the cost of a steeper computational search; the paper only demonstrates pairs.
  • The simulations use white noise only; red noise, overlapping sources, and a stochastic background, which the authors list as necessary follow-ups, will likely weaken the quoted sub-parsec fractions, though the structural advantage of two-dimensional over one-dimensional inference should persist.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 5 minor

Summary. The paper proposes a two-dimensional method for inferring pulsar distances from pulsar timing array (PTA) data, using the pulsar-term phase information from multiple continuous gravitational-wave (CGW) sources. For each pair of pulsars, the method builds a joint posterior over the two distances by combining the phase-distance relations of all CGW sources, then marginalizes one distance to obtain a one-dimensional posterior for the other. The authors simulate an SKA-era PTA with white noise and known CGW sky locations, and report that with a few sources (N=4-5) pulsars at ~1 kpc can achieve sub-parsec distance precision, claimed to be a substantial improvement over the one-dimensional combination method of Yu & Pan. The central quantitative claim is that, e.g., four sources yield sub-parsec precision in about 90% of realizations for J0613-0200 at sigma_n=50 ns.

Significance. If the claims hold, the method would be a valuable contribution to PTA data analysis, potentially enabling high-precision pulsar distances without relying solely on parallax and thereby improving CGW host-galaxy identification and multimessenger studies. The paper provides an explicit simulation, including the source population model, noise level, and waveform treatment, and the core idea of retaining two-dimensional distance structure to avoid the loss of multimodality in one-dimensional marginalization is well motivated. The main reservation is that the headline results are computed after a data-dependent selection step and are not validated by coverage checks; without those, the quantitative improvement over prior work is not established.

major comments (3)
  1. [Simulations and distance inference (paragraph after Eq. (3))] The paper defines Delta L_p as the half-width of the 68% credible interval of the posterior PDF_p^(q) with the smallest half-width over all partners q. With ~84 partners, this is a post-hoc order statistic: the minimum of many correlated interval widths is systematically narrower than a typical pair posterior, even under a correct model. The reported distribution of Delta L_p in Fig. 2 and the text's 'four GW sources yield sub-parsec precision in about 90% of the realizations' is therefore a post-selection statistic, not a typical achievable precision. The paper never reports empirical coverage — the fraction of realizations in which the true injected L_p lies inside the selected 68% interval. Without coverage, 'sub-parsec precision' does not imply that the distance is correctly estimated. Moreover, the 1D comparison in Fig. 5 does not include an analogous partner-selection step, making
  2. [Simulations and distance inference / Supplementary Eq. (4)] The manuscript states that the pulsar-term initial phases {Phi_p} are sampled as free parameters, while simultaneously imposing Gaussian priors on the distances L_p. If posterior samples of Phi_p are mapped to L_p via Eq. (4) without a Jacobian or prior reweighting, the resulting L_p distribution is not the posterior with the stated Gaussian prior. The paper does not describe the required reweighting or demonstrate that the sampling is equivalent to the intended inference on L_p. This ambiguity affects the validity of all reported posteriors, including Fig. 1 and the Delta L_p values. Please specify the exact procedure (e.g., how the Gaussian prior enters the mapping, whether the Jacobian is included) or show that the result is independent of the prior choice.
  3. [Simulations and distance inference (third paragraph)] The analysis assumes that the sky locations of the CGW sources are known exactly. This is an acknowledged input, but the manuscript does not test the sensitivity of the inferred distances to errors in those positions. Since Eq. (4) depends on the dot product Omega-hat·p-hat, even a moderately inaccurate sky position can bias the phase-distance mapping and hence the inferred L_p. Given that the stated goal includes host-galaxy identification, the authors should quantify the robustness of Delta L_p and of the coverage to plausible positional uncertainties (e.g., by adding a small scatter to theta, phi in recovery). This is not a circularity issue, but it affects the practical applicability of the method.
minor comments (5)
  1. [Eq. (4)] The typesetting of Eq. (4) is difficult to parse; the numerator and denominator should be re-expressed or clarified, and the domain of validity (circular, GW-driven binary) should be stated explicitly.
  2. [Fig. 1] The axes are labeled L1 and L2 but the caption does not define which pulsar corresponds to which label. Please state L1 = J0030+0451 and L2 = J0613-0200 in the caption.
  3. [Simulations and distance inference] The priors on the source parameters are said to be the same as in Ref. [37], but they are not listed. For reproducibility, include the prior ranges or functional forms in the text or supplementary material.
  4. [Fig. 2 and Fig. 4] The violin plots would be more informative if they included the median and quartile markers, and if the fraction of realizations with Delta L_p < 1 pc were stated explicitly for each N and sigma_n.
  5. [Conclusions (last paragraph)] The paper mentions the simplified setting (white noise, known sky positions) but does not discuss the computational burden of the MCMC over the full parameter space (85 pulsars, N sources, plus pulsar-term phases). A brief note on expected runtime or scalability would be useful.

Circularity Check

0 steps flagged

No significant circularity: the 2D posterior construction is an independent statistical method; self-citations are contextual and the reported precision is an empirical simulation statistic, not an equation identity.

full rationale

The paper's derivation chain is a self-consistency injection study: simulated CGW signals are generated with a waveform model, and the same physical relation (Eq. 4) is used in the recovery to map pulsar-term phase posteriors to distances. This is standard model validation, not circularity, because the recovery does not use the true injected distances as an input; it infers them from noise-realized data. The two-dimensional posterior (Eq. 2) and its marginalization (Eq. 3) define a new analysis method, and the reported sub-parsec precision is an empirical distribution over simulated realizations, not a quantity forced by the equations. The analytic phase-distance relation is not fitted to the target distances; it is a physical waveform identity. Self-citations in the introduction (e.g., refs. [17], [19], [29]) are contextual references to prior standard-siren and dark-siren work by the same group; they are not load-bearing for the distance-inference claim. The choice to adopt the partner posterior with the smallest 68% credible interval is a post-selection statistic rather than a circular reduction: the smallest interval is not an input to Eq. (2), and the reported distribution still depends on the simulated data and noise. It could bias the quoted precision, but it is a soundness/statistics concern, not a derivation that reduces to its own inputs.

Axiom & Free-Parameter Ledger

6 free parameters · 5 axioms · 0 invented entities

All claimed precision numbers rest on the simulated catalog, the known-sky-position assumption, the waveform model, and the pair-selection rule. No new physical entities are introduced; the method is a statistical construction.

free parameters (6)
  • timing noise sigma_n = 50 ns and 100 ns
    Chosen from SKA forecasts; central precision values depend on this.
  • distance prior width sigma_p = O(10 pc) at 1 kpc
    Parallax prior from SKA timing-parallax forecast; nearly all reported sub-pc results are relative to this prior.
  • chirp mass M = 5e9 Msun
    Fixed in simulation though real CGW sources will have a mass distribution; affects phase evolution and distance resolution.
  • luminosity distance d_L = 1 Gpc
    Fiducial for loud detectable CGWs; amplitude and phase evolution depend on it.
  • observing campaign = 20 yr, 2-week cadence
    SKA-era assumed schedule; longer baseline would change precision.
  • number of CGW sources N = 2-5
    The method's performance is parameterized by N; central claim is 'few sources'.
axioms (5)
  • domain assumption Sky positions of CGW sources are known and fixed
    Section 'Simulations and distance inference': 'we assume that the sky locations of the CGW sources used for pulsar-distance inference are known and therefore treat their polar and azimuthal angles as fixed parameters.' If wrong, distances are biased.
  • domain assumption Circular, GW-driven SMBHB waveform with no eccentricity
    The residual model and Eq. (4) assume quasi-circular inspirals; eccentric or evolving binaries would break the phase-distance relation.
  • domain assumption White Gaussian timing noise only; no SGWB/overlapping CGWs/red noise
    'We include a CGW signal in white noise' (Conclusions); the impact of additional signals and noise is explicitly left to future work.
  • ad hoc to paper Pulsar-term initial phases {Phi_p} are sampled as free parameters
    Section 'Simulations': 'we sample them as free parameters to avoid the highly oscillatory likelihood structure'. This introduces unphysical degrees of freedom; the mapping Eq. (4) is used to recover the distance, but any mis-specification of the phase prior propagates in.
  • domain assumption Frequency-difference information between Earth and pulsar terms is negligible
    Section 'Simulations': 'We neglect the pulsar-distance information carried by the frequency differences...' without a quantitative demonstration in the main text; supplementary asserts it is negligible.

pith-pipeline@v1.3.0-alltime-deepseek · 10271 in / 15897 out tokens · 150876 ms · 2026-08-03T17:02:50.354561+00:00 · methodology

0 comments
read the original abstract

Pulsar timing arrays (PTAs) are limited in localizing nanohertz continuous gravitational waves (CGWs) by uncertainties in pulsar distances. We introduce a method to infer pulsar distances in two dimensions, using phase information from the pulsar terms of multiple CGW sources. Our approach can enhance distance precision and, in some cases, achieve order-of-magnitude improvements relative to existing one-dimensional distance-inference methods. Using simulations of an SKA-era PTA with realistic parallax-based distance priors, we demonstrate that pulsars at $\sim 1$ kpc can achieve sub-parsec distance precision with only a few CGW sources. Such improvements in pulsar-distance precision have important implications for CGW host-galaxy identification and multimessenger observational prospects.

Figures

Figures reproduced from arXiv: 2512.10729 by Jing-Fei Zhang, Ji-Yu Song, Ling-Feng Wang, Si-Ren Xiao, Xin Zhang, Yue Shao.

Figure 1
Figure 1. Figure 1: FIG. 1. Example of pulsar distance inference from the joint [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Distributions of the pulsar-distance uncertainty ∆ [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Posterior distributions of [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: shows that nearer pulsars can reach tighter distance constraints with fewer GW sources, whereas for more distant pulsars sub-parsec precision is achieved only in a smaller fraction of realizations even with an increased number of sources. 2 3 4 5 N 2.0 1.5 1.0 0.5 0.0 l o g 1 0 ( ¢ L p = p c ) (a) 2 3 4 5 N 1.0 0.5 0.0 0.5 1.0 1.5 (b) ¾n = 50 ns ¾n = 100 ns FIG. 4. Same as [PITH_FULL_IMAGE:figures/full_fi… view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Comparison between the 2D combination results and 1D combination results for pulsar-distance inference. Same as [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 3 Pith papers

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

  1. Archival Multiband Gravitational-Wave Signals from Massive Black Hole Binary Mergers

    astro-ph.HE 2026-04 unverdicted novelty 7.0

    Massive black hole binary mergers produce orphaned low-frequency signals in PTA pulsar terms that can be stacked for archival multiband gravitational-wave detection.

  2. Model-independent H0 from GWTC-4 standard sirens and TDCOSMO 2025 strong lensing time delays

    astro-ph.CO 2026-06 unverdicted novelty 5.0

    Combining GWTC-4 standard sirens with TDCOSMO2025 lensing data under the distance sum rule yields H0 = 83.78 +12.53/-10.23 km/s/Mpc (13.6% precision) in one configuration, consistent with both Planck and SH0ES.

  3. Synergy between CSST and future gravitational-wave detectors: Probing primordial black holes by cross-correlating dark sirens with galaxies

    astro-ph.CO 2026-06 unverdicted novelty 4.0

    Cross-correlating CSST galaxies with mock GW catalogs from ET2CE and BDET2CE networks can detect PBH merger fractions above ~40% and ~20% respectively via clustering bias differences.

Reference graph

Works this paper leans on

55 extracted references · 46 linked inside Pith · cited by 3 Pith papers

  1. [1]

    Agazie et al

    G. Agazie et al. (NANOGrav), The NANOGrav 15 yr Data Set: Evidence for a Gravitational-wave Background, Astrophys. J. Lett. 951, L8 (2023), arXiv:2306.16213 [astro-ph.HE]

    Pith/arXiv arXiv 2023

  2. [2]

    Antoniadis et al

    J. Antoniadis et al. (EPTA), The second data release from the European Pulsar Timing Array III. Search for gravitational wave signals, Astron. Astrophys. 678, A50 (2023), arXiv:2306.16214 [astro-ph.HE]

    Pith/arXiv arXiv 2023

  3. [3]

    D. J. Reardon et al. , Search for an Isotropic Gravitational-wave Background with the Parkes Pul- sar Timing Array, Astrophys. J. Lett. 951, L6 (2023), arXiv:2306.16215 [astro-ph.HE]

    Pith/arXiv arXiv 2023

  4. [4]

    Xu et al., Searching for the Nano-Hertz Stochastic Gravitational Wave Background with the Chinese Pulsar Timing Array Data Release I, Res

    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), arXiv:2306.16216 [astro-ph.HE]

    Pith/arXiv arXiv 2023

  5. [5]

    M. T. Miles et al., The MeerKAT Pulsar Timing Array: first data release, Mon. Not. Roy. Astron. Soc. 519, 3976 (2023), arXiv:2212.04648 [astro-ph.HE]

    Pith/arXiv arXiv 2023

  6. [6]

    Siemens, V

    X. Siemens, V. Mandic, and J. Creighton, Gravitational wave stochastic background from cosmic (super)strings, Phys. Rev. Lett. 98, 111101 (2007), arXiv:astro- ph/0610920

    arXiv 2007

  7. [7]

    Schwaller, Gravitational Waves from a Dark Phase Transition, Phys

    P. Schwaller, Gravitational Waves from a Dark Phase Transition, Phys. Rev. Lett. 115, 181101 (2015), arXiv:1504.07263 [hep-ph]

    Pith/arXiv arXiv 2015

  8. [8]

    Cai, Z.-K

    R.-G. Cai, Z.-K. Guo, J. Liu, L. Liu, and X.-Y. Yang, Pri- mordial black holes and gravitational waves from para- metric amplification of curvature perturbations, JCAP 06, 013, arXiv:1912.10437 [astro-ph.CO]

    Pith/arXiv arXiv 1912

  9. [9]

    Ellis and M

    J. Ellis and M. Lewicki, Cosmic String Interpretation of NANOGrav Pulsar Timing Data, Phys. Rev. Lett. 126, 041304 (2021), arXiv:2009.06555 [astro-ph.CO]

    Pith/arXiv arXiv 2021

  10. [10]

    Sesana, A

    A. Sesana, A. Vecchio, and M. Volonteri, Gravitational waves from resolvable massive black hole binary sys- tems and observations with Pulsar Timing Arrays, Mon. Not. Roy. Astron. Soc.394, 2255 (2009), arXiv:0809.3412 [astro-ph]

    Pith/arXiv arXiv 2009

  11. [11]

    Babak and A

    S. Babak and A. Sesana, Resolving multiple supermassive black hole binaries with pulsar timing arrays, Phys. Rev. D 85, 044034 (2012), arXiv:1112.1075 [astro-ph.CO]

    Pith/arXiv arXiv 2012

  12. [12]

    X.-J. Zhu, L. Wen, G. Hobbs, Y. Zhang, Y. Wang, D. R. Madison, R. N. Manchester, M. Kerr, P. A. Rosado, and J.-B. Wang, Detection and localization of single-source gravitational waves with pulsar timing arrays, Mon. Not. Roy. Astron. Soc. 449, 1650 (2015), arXiv:1502.06001 [astro-ph.IM]

    Pith/arXiv arXiv 2015

  13. [13]

    Wang and S

    Y. Wang and S. D. Mohanty, Pulsar Timing Array Based Search for Supermassive Black Hole Binaries in the Square Kilometer Array Era, Phys. Rev. Lett. 118, 151104 (2017), [Erratum: Phys.Rev.Lett. 124, 169901 (2020)], arXiv:1611.09440 [astro-ph.IM]

    Pith/arXiv arXiv 2017

  14. [14]

    C. M. F. Mingarelli, T. J. W. Lazio, A. Sesana, J. E. Greene, J. A. Ellis, C.-P. Ma, S. Croft, S. Burke-Spolaor, and S. R. Taylor, The Local Nanohertz Gravitational- Wave Landscape From Supermassive Black Hole Bina- 5 ries, Nature Astron. 1, 886 (2017), arXiv:1708.03491 [astro-ph.GA]

    Pith/arXiv arXiv 2017

  15. [15]

    D. L. Jow and U.-L. Pen, Measuring Cosmic Expansion with Diffractive Gravitational Scintillation of Nanohertz Gravitational Waves, Phys. Rev. Lett. 134, 131001 (2025), arXiv:2407.03214 [astro-ph.CO]

    Pith/arXiv arXiv 2025

  16. [16]

    C. Yan, W. Zhao, and Y. Lu, On Using Inspiraling Super- massive Binary Black Holes in the PTA Frequency Band as Standard Sirens to Constrain Dark Energy, Astrophys. J. 889, 10 (2020), arXiv:1912.04103 [astro-ph.CO]

    Pith/arXiv arXiv 2020

  17. [17]

    L.-F. Wang, Y. Shao, S.-R. Xiao, J.-F. Zhang, and X. Zhang, Ultra-low-frequency gravitational waves from individual supermassive black hole binaries as standard sirens, JCAP 05, 095, arXiv:2201.00607 [astro-ph.CO]

    Pith/arXiv arXiv

  18. [18]

    Jin, S.-S

    S.-J. Jin, S.-S. Xing, Y. Shao, J.-F. Zhang, and X. Zhang, Joint constraints on cosmological parame- ters using future multi-band gravitational wave standard siren observations*, Chin. Phys. C 47, 065104 (2023), arXiv:2301.06722 [astro-ph.CO]

    Pith/arXiv arXiv 2023

  19. [19]

    Jin, J.-Y

    S.-J. Jin, J.-Y. Song, T.-Y. Sun, S.-R. Xiao, H. Wang, L.- F. Wang, J.-F. Zhang, and X. Zhang, Gravitational wave standard sirens: A brief review of cosmological parameter estimation, Sci. China Phys. Mech. Astron. 69, 220401 (2026), arXiv:2507.12965 [astro-ph.CO]

    arXiv 2026

  20. [20]

    B. F. Schutz, Determining the Hubble Constant from Gravitational Wave Observations, Nature 323, 310 (1986)

    1986

  21. [21]

    Jin, Y.-Z

    S.-J. Jin, Y.-Z. Zhang, J.-Y. Song, J.-F. Zhang, and X. Zhang, Taiji-TianQin-LISA network: Precisely mea- suring the Hubble constant using both bright and dark sirens, Sci. China Phys. Mech. Astron.67, 220412 (2024), arXiv:2305.19714 [astro-ph.CO]

    Pith/arXiv arXiv 2024

  22. [22]

    Song, L.-F

    J.-Y. Song, L.-F. Wang, Y. Li, Z.-W. Zhao, J.-F. Zhang, W. Zhao, and X. Zhang, Synergy between CSST galaxy survey and gravitational-wave observation: Inferring the Hubble constant from dark standard sirens, Sci. China Phys. Mech. Astron.67, 230411 (2024), arXiv:2212.00531 [astro-ph.CO]

    Pith/arXiv arXiv 2024

  23. [23]

    Song, J.-Z

    J.-Y. Song, J.-Z. Qi, J.-F. Zhang, and X. Zhang, Model-independent H 0 within FLR W: Joint Constraints from GWTC-3 Standard Sirens and Strong Lensing Time Delays, Astrophys. J. Lett. 985, L44 (2025), arXiv:2503.10346 [astro-ph.CO]

    Pith/arXiv arXiv 2025

  24. [24]

    A. G. Abac et al. (LIGO Scientific, VIRGO, KAGRA), GWTC-4.0: Constraints on the Cosmic Expansion Rate and Modified Gravitational-wave Propagation, (2025), arXiv:2509.04348 [astro-ph.CO]

    Pith/arXiv arXiv 2025

  25. [25]

    Song, G.-H

    J.-Y. Song, G.-H. Du, T.-N. Li, L.-F. Wang, J.-Z. Qi, J.-F. Zhang, and X. Zhang, Gravitational wave stan- dard sirens from GWTC-3 combined with DESI DR2 and DESY5: A late-universe probe of the Hubble con- stant and dark energy, (2025), arXiv:2511.12017 [astro- ph.CO]

    arXiv 2025

  26. [26]

    S. R. Taylor, E. A. Huerta, J. R. Gair, and S. T. McWilliams, Detecting eccentric supermassive black hole binaries with pulsar timing arrays: Resolv- able source strategies, Astrophys. J. 817, 70 (2016), arXiv:1505.06208 [gr-qc]

    Pith/arXiv arXiv 2016

  27. [27]

    J. M. Goldstein, A. Sesana, A. M. Holgado, and J. Veitch, Associating host galaxy candidates to massive black hole binaries resolved by pulsar timing arrays, Mon. Not. Roy. Astron. Soc. 485, 248 (2019), arXiv:1812.02670 [astro- ph.IM]

    Pith/arXiv arXiv 2019

  28. [28]

    Petrov, S

    P. Petrov, S. R. Taylor, M. Charisi, and C.-P. Ma, Identi- fying the Host Galaxies of Supermassive Black Hole Bina- ries Found by Pulsar Timing Arrays, Astrophys. J. 976, 129 (2024), arXiv:2406.04409 [astro-ph.GA]

    Pith/arXiv arXiv 2024

  29. [29]

    S.-R. Xiao, Y. Shao, L.-F. Wang, J.-Y. Song, L. Feng, J.-F. Zhang, and X. Zhang, Nanohertz gravitational waves from a quasar-based supermassive black hole bi- nary population model as dark sirens, JCAP 04, 060, arXiv:2408.00609 [astro-ph.CO]

    Pith/arXiv arXiv

  30. [30]

    R. J. Truant, D. Izquierdo-Villalba, A. Sesana, G. M. Shaifullah, M. Bonetti, D. Spinoso, and S. Bonoli, Light- ing up the nano-hertz gravitational wave sky: oppor- tunities and challenges of multimessenger astronomy with PTA experiments, (2025), arXiv:2504.01074 [astro- ph.GA]

    Pith/arXiv arXiv 2025

  31. [31]

    Taylor, J

    S. Taylor, J. Ellis, and J. Gair, Accelerated Bayesian model-selection and parameter-estimation in continuous gravitational-wave searches with pulsar-timing arrays, Phys. Rev. D 90, 104028 (2014), arXiv:1406.5224 [gr-qc]

    Pith/arXiv arXiv 2014

  32. [32]

    Charisi, S

    M. Charisi, S. R. Taylor, C. A. Witt, and J. Runnoe, Ef- ficient Large-Scale, Targeted Gravitational-Wave Probes of Supermassive Black-Hole Binaries, Phys. Rev. Lett. 132, 061401 (2024), arXiv:2304.03786 [gr-qc]

    Pith/arXiv arXiv 2024

  33. [33]

    Grunthal, N

    K. Grunthal, N. Porayko, D. J. Champion, and M. Kramer, The role of distant pulsars in the de- tectability of continuous gravitational waves, (2025), arXiv:2512.04589 [astro-ph.HE]

    arXiv 2025

  34. [34]

    Boyle and U.-L

    L. Boyle and U.-L. Pen, Pulsar timing arrays as imag- ing gravitational wave telescopes: angular resolution and source (de)confusion, Phys. Rev. D 86, 124028 (2012), arXiv:1010.4337 [astro-ph.HE]

    Pith/arXiv arXiv 2012

  35. [35]

    Corbin and N

    V. Corbin and N. J. Cornish, Pulsar Timing Array Ob- servations of Massive Black Hole Binaries, arXiv e-prints (2010), arXiv:1008.1782 [astro-ph.HE]

    Pith/arXiv arXiv 2010

  36. [36]

    K. J. Lee et al., Gravitational wave astronomy of sin- gle sources with a pulsar timing array, Mon. Not. Roy. Astron. Soc. 414, 3251 (2011), arXiv:1103.0115 [astro- ph.HE]

    Pith/arXiv arXiv 2011

  37. [37]

    Kato and K

    R. Kato and K. Takahashi, Realistic assessment of a sin- gle gravitational wave source localization taking into ac- count precise pulsar distances with pulsar timing arrays, Phys. Rev. D 113, 022001 (2026), arXiv:2506.02819 [gr- qc]

    arXiv 2026

  38. [38]

    A. C. Tsai, D. L. Jow, and U.-L. Pen, Reach- ing diffraction-limited localization with coherent ptas, (2025), arXiv:2512.10795 [astro-ph.IM]

    Pith/arXiv arXiv 2025

  39. [39]

    A. T. Deller et al., Microarcsecond VLBI Pulsar Astrom- etry with PSR π II. Parallax Distances for 57 Pulsars, Astrophys. J. 875, 100 (2019), arXiv:1808.09046 [astro- ph.IM]

    Pith/arXiv arXiv 2019

  40. [40]

    Ding et al., The MSPSR π catalogue: VLBA astrom- etry of 18 millisecond pulsars, Mon

    H. Ding et al., The MSPSR π catalogue: VLBA astrom- etry of 18 millisecond pulsars, Mon. Not. Roy. Astron. Soc. 519, 4982 (2023), arXiv:2212.06351 [astro-ph.HE]

    Pith/arXiv arXiv 2023

  41. [41]

    D. J. Reardon et al., The Neutron Star Mass, Distance, and Inclination from Precision Timing of the Brilliant Millisecond Pulsar J0437-4715, Astrophys. J. Lett. 971, L18 (2024), arXiv:2407.07132 [astro-ph.HE]

    Pith/arXiv arXiv 2024

  42. [42]

    Smits, S

    R. Smits, S. J. Tingay, N. Wex, M. Kramer, and B. Stap- pers, Prospects for accurate distance measurements of pulsars with the SKA: enabling fundamental physics, Astron. Astrophys. 528, A108 (2011), arXiv:1101.5971 [astro-ph.IM]

    Pith/arXiv arXiv 2011

  43. [43]

    McGrath, D

    C. McGrath, D. J. D’Orazio, and J. Creighton, Measur- ing the Hubble constant with double gravitational wave sources in pulsar timing, Mon. Not. Roy. Astron. Soc. 6 517, 1242 (2022), arXiv:2208.06495 [astro-ph.CO]

    Pith/arXiv arXiv 2022

  44. [44]

    Yu and Z

    J. Yu and Z. Pan, Subparsec precision measurement of pulsar distances with nanohertz gravitational waves, Phys. Rev. D 112, 023012 (2025), arXiv:2503.23017 [astro-ph.HE]

    Pith/arXiv arXiv 2025

  45. [45]

    Agarwal et al., The NANOGrav 15 yr Data Set: Tar- geted Searches for Supermassive Black Hole Binaries, (2025), arXiv:2508.16534 [astro-ph.HE]

    N. Agarwal et al., The NANOGrav 15 yr Data Set: Tar- geted Searches for Supermassive Black Hole Binaries, (2025), arXiv:2508.16534 [astro-ph.HE]

    Pith/arXiv arXiv 2025

  46. [46]

    B´ ecsy, N

    B. B´ ecsy, N. J. Cornish, and M. C. Digman, Fast Bayesian analysis of individual binaries in pulsar tim- ing array data, Phys. Rev. D 105, 122003 (2022), arXiv:2204.07160 [gr-qc]

    Pith/arXiv arXiv 2022

  47. [47]

    Tian, Y.-C

    L.-W. Tian, Y.-C. Bi, Y.-M. Wu, and Q.-G. Huang, Tar- geted search for an individual SMBHB in NANOGrav 15- year and EPTA DR2 data sets, (2025), arXiv:2508.14742 [astro-ph.GA]

    arXiv 2025

  48. [48]

    Bardati, J

    J. Bardati, J. J. Ruan, D. Haggard, and M. Tremmel, Signatures of Massive Black Hole Merger Host Galax- ies from Cosmological Simulations. I. Unique Galaxy Morphologies in Imaging, Astrophys. J. 961, 34 (2024), arXiv:2308.03828 [astro-ph.GA]

    arXiv 2024

  49. [49]

    Horlaville, J

    P. Horlaville, J. J. Ruan, M. Eracleous, J. Bardati, J. C. Runnoe, and D. Haggard, Predicting Potential Host Galaxies of Supermassive Black Hole Binaries Based on Stellar Kinematics in Archival IFU Surveys, (2025), arXiv:2504.21145 [astro-ph.GA]

    Pith/arXiv arXiv 2025

  50. [50]

    Agazie et al

    G. Agazie et al. (NANOGrav), The NANOGrav 15 yr Data Set: Bayesian Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries, Astrophys. J. Lett. 951, L50 (2023), arXiv:2306.16222 [astro-ph.HE]

    Pith/arXiv arXiv 2023

  51. [51]

    Antoniadis et al

    J. Antoniadis et al. (EPTA), The second data release from the European Pulsar Timing Array - I. The dataset and timing analysis, Astron. Astrophys.678, A48 (2023), arXiv:2306.16224 [astro-ph.HE]

    Pith/arXiv arXiv 2023

  52. [52]

    R. N. Manchester, G. B. Hobbs, A. Teoh, and M. Hobbs, The Australia Telescope National Facility pul- sar catalogue, Astron. J. 129, 1993 (2005), arXiv:astro- ph/0412641

    arXiv 1993

  53. [53]

    Janssen et al., Gravitational wave astronomy with the SKA, PoS AASKA14, 037 (2015), arXiv:1501.00127 [astro-ph.IM]

    G. Janssen et al., Gravitational wave astronomy with the SKA, PoS AASKA14, 037 (2015), arXiv:1501.00127 [astro-ph.IM]

    Pith/arXiv arXiv 2015

  54. [54]

    Cella, S

    K. Cella, S. R. Taylor, and L. Z. Kelley, Host galaxy demographics of individually detectable supermassive black-hole binaries with pulsar timing arrays, Class. Quant. Grav. 42, 025021 (2025), arXiv:2407.01659 [astro- ph.GA]

    Pith/arXiv arXiv 2025

  55. [55]

    J. A. Ellis, A Bayesian analysis pipeline for continuous GW sources in the PTA band, Class. Quant. Grav. 30, 224004 (2013), arXiv:1305.0835 [astro-ph.IM]. 7 Supplementary Material I. ANAL YTIC RELA TION BETWEEN PULSAR DIST ANCE AND THE PULSAR-TERM INITIAL PHASE Here we provide the explicit expression for the analytic mapping between the pulsar distance Lp...

    Pith/arXiv arXiv 2013