pith. sign in

arxiv: 2606.26945 · v1 · pith:IPS3TCFJnew · submitted 2026-06-25 · 🌀 gr-qc · astro-ph.HE

A geometric multimessenger consistency test of radiative and near-zone gravity with LISA and SKA

Bhooshan Gadre This is my paper

Pith reviewed 2026-06-26 03:44 UTC · model grok-4.3

classification 🌀 gr-qc astro-ph.HE
keywords multimessengergravitational wavespulsar timingorbital inclinationShapiro delayLISASKAconsistency test
0
0 comments X

The pith

A difference between radio-derived and gravitational-wave-derived orbital inclinations would indicate a mismatch between near-zone and radiative gravity.

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

The paper proposes using the same compact binary pulsar observed in radio timing and gravitational waves to compare its orbital inclination extracted from two independent messengers. Radio timing measures the inclination via the shape of the Shapiro delay, while the gravitational-wave signal measures it via the ratio of tensor polarization amplitudes. General relativity predicts these two values agree exactly, so the authors define a residual difference ε and estimate how precisely it can be constrained with LISA and SKA observations. For a double neutron star at 5 kpc observed for four years the precision reaches a few parts in a thousand; a pulsar-black-hole system could do better. The test is limited mainly by the gravitational-wave polarimetry, with radio timing supplying the near-zone reference.

Core claim

The central claim is that the common-epoch residual ε(t0) defined as the difference between the sine of the orbital inclination from the Shapiro-delay shape parameter and the sine of the orbital inclination from the gravitational-wave tensor polarization amplitude ratio must vanish in general relativity; any statistically significant nonzero value signals either an unmodeled systematic or an inconsistency between the near-zone and radiative descriptions of gravity. The paper quantifies the attainable uncertainty on this residual for representative LISA-SKA systems using a seven-parameter radio timing Fisher matrix and a sky-averaged LISA sensitivity that includes the Galactic foreground.

What carries the argument

The common-epoch residual ε(t0) = s_Shapiro(t0) - s_GW(t0) constructed from independent inclination measurements of the same binary.

If this is right

  • A null result on ε at the estimated precision would confirm that the near-zone geometry inferred from radio timing is consistent with the radiative geometry inferred from gravitational waves.
  • A nonzero ε at detectable levels would require either improved modeling of systematics or a revision of the assumption that the two descriptions of gravity must agree.
  • The test remains useful even if only a handful of double neutron star or pulsar-black-hole systems are detected, because radio timing supplies the reference measurement while gravitational-wave polarimetry sets the limiting uncertainty.
  • The reported uncertainties define quantitative targets that future joint Bayesian analyses of LISA and SKA data must reach.

Where Pith is reading between the lines

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

  • The same residual construction could be applied to other multimessenger binaries once additional detectors improve polarimetry sensitivity.
  • If a nonzero ε is found, distinguishing between an astrophysical systematic and a genuine mismatch in gravity descriptions would require cross-checks with independent strong-field observables.
  • The method supplies a direct geometric link between the two messengers that is independent of the specific waveform model details beyond the inclination extraction.

Load-bearing premise

The orbital inclination extracted from the Shapiro-delay shape parameter equals the inclination extracted from the tensor polarization amplitude ratio, with no additional modeling differences between the two messengers.

What would settle it

A statistically significant nonzero measurement of ε in any compact binary pulsar observed simultaneously in radio timing and gravitational waves with the stated precision would falsify the consistency prediction.

Figures

Figures reproduced from arXiv: 2606.26945 by Bhooshan Gadre.

Figure 1
Figure 1. Figure 1: FIG. 1. Single-channel inclination precision and cross-messenger disagreement sensitivity in the baseline model with [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. DNS sensitivity map in distance and inclination. The [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Illustrative eccentricity scan for the DNS and [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Stress tests for the main forecast assumptions. Left: multiplicative inflation of the LISA polarimetric inclination [PITH_FULL_IMAGE:figures/full_fig_p010_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. A simple time-averaged latitude proxy for Galactic [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
read the original abstract

Compact binary pulsars observed both through precision radio timing and low-frequency gravitational waves offer a direct way to compare the same binary geometry with two independent messengers. We propose a multimessenger consistency test based on the orbital inclination, measured from the Shapiro-delay shape parameter in radio timing and from the tensor polarization amplitude ratio in the gravitational-wave signal. Defining the common-epoch residual $\eps(t_0)=s_{\rm Shapiro}(t_0)-s_{\rm GW}(t_0)$, general relativity predicts $\eps=0$, while a nonzero value would indicate either an unmodeled systematic or a mismatch between the near-zone and radiative descriptions of gravity. We estimate the attainable precision on this quantity for representative LISA--SKA compact binary pulsars using a seven-parameter timing Fisher matrix and a sky-averaged LISA sensitivity curve including the Galactic foreground. We adopt a conservative radio baseline, $\sigma_{\rm TOA}=1\,\mu{\rm s}$ and $N_{\rm eff}=10^4$, intended to summarize radiometer noise, jitter, residual dispersion-measure and scattering effects, profile evolution, and cadence losses after wideband timing. For systems at $d=5\,{\rm kpc}$ observed for four years, we find $\sigma_\eps\simeq4\times10^{-3}$ for a favorable double neutron star and $\sigma_\eps\simeq9\times10^{-4}$ for a hypothetical pulsar--black-hole system. The former is the more robust astrophysical benchmark; the latter illustrates the reach if such a high-SNR chirping source is discovered. The useful cases remain limited mainly by gravitational-wave polarimetry, while radio timing supplies the near-zone reference measurement of the inclination. These results define a quantitative target for future joint Bayesian analyses of compact binary pulsars observed in both radio and gravitational waves.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript proposes a multimessenger consistency test for general relativity using orbital inclination from radio pulsar timing (Shapiro-delay shape parameter s_Shapiro) and gravitational-wave polarimetry (tensor amplitude ratio s_GW) for compact binaries observable by LISA and SKA. It defines the residual ε(t0) = s_Shapiro(t0) - s_GW(t0), states that GR predicts ε=0, and estimates attainable precision σ_ε using a seven-parameter radio timing Fisher matrix and sky-averaged LISA sensitivity curve with Galactic foreground, reporting σ_ε ≃ 4×10^{-3} for a double neutron star and ≃9×10^{-4} for a pulsar-black-hole system at 5 kpc over four years with conservative radio baseline σ_TOA=1 μs and N_eff=10^4.

Significance. If the equivalence of the two inclination extractions holds, the test would supply a direct geometric comparison between near-zone radio and radiative GW descriptions of gravity, with the provided Fisher estimates offering quantitative targets for joint analyses. The conservative radio baseline and explicit separation of messengers strengthen the proposal's robustness as a falsifiable consistency check.

major comments (2)
  1. [Abstract] Abstract / Definition of ε: The claim that GR predicts ε=0 (with nonzero ε indicating a near-zone/radiative mismatch) requires that s_Shapiro(t0) and s_GW(t0) extract identical values of the orbital inclination with no residual offsets from definition, coordinate choice, or higher-order corrections. The seven-parameter radio Fisher matrix and sky-averaged LISA polarimetry are constructed independently; no explicit cross-derivation or invariance check between the two extraction pipelines is provided, so any unaccounted mapping difference would produce spurious nonzero ε even in exact GR.
  2. [Abstract] Precision estimates paragraph: The reported σ_ε values rely on the seven-parameter timing Fisher matrix and sky-averaged LISA curve, but the abstract supplies no validation that the matrix marginalizes correctly over inclination covariances or that the foreground model does not introduce systematic offsets in the GW amplitude ratio; this directly affects the load-bearing claim that radio supplies the reference measurement.
minor comments (2)
  1. [Abstract] The abstract refers to 'the shape parameter s' without clarifying whether s is the raw Shapiro shape parameter or already mapped to sin i; explicit notation would aid comparison with the GW amplitude ratio.
  2. [Abstract] The conservative radio baseline is stated to summarize multiple effects, but a brief parenthetical list of the dominant contributors would improve reproducibility of the N_eff=10^4 choice.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thoughtful and constructive report. The comments correctly identify areas where the presentation of the ε definition and the supporting precision estimates can be strengthened for clarity. We address each major comment below and indicate the revisions that will be incorporated.

read point-by-point responses

  1. Referee: [Abstract] Abstract / Definition of ε: The claim that GR predicts ε=0 (with nonzero ε indicating a near-zone/radiative mismatch) requires that s_Shapiro(t0) and s_GW(t0) extract identical values of the orbital inclination with no residual offsets from definition, coordinate choice, or higher-order corrections. The seven-parameter radio Fisher matrix and sky-averaged LISA polarimetry are constructed independently; no explicit cross-derivation or invariance check between the two extraction pipelines is provided, so any unaccounted mapping difference would produce spurious nonzero ε even in exact GR.

    Authors: We agree that a clean interpretation of ε=0 in GR requires explicit justification that both parameters extract the same physical inclination without unaccounted offsets. Both s_Shapiro and s_GW are defined at leading order to measure sin i for the same orbital plane; the radio timing uses the standard Shapiro-delay shape parameter while the GW polarimetry uses the tensor amplitude ratio, both of which reduce to the identical geometric inclination in the GR limit. The independent construction of the two Fisher matrices is intentional, as each messenger supplies its own measurement, but the underlying GR prediction equates the extracted inclinations. We will revise the manuscript (including a brief addition near the definition of ε) to state this equivalence explicitly, note the coordinate conventions adopted, and confirm that higher-order post-Newtonian corrections lie below the reported σ_ε values. This addresses the concern without altering the core proposal. revision: yes

  2. Referee: [Abstract] Precision estimates paragraph: The reported σ_ε values rely on the seven-parameter timing Fisher matrix and sky-averaged LISA curve, but the abstract supplies no validation that the matrix marginalizes correctly over inclination covariances or that the foreground model does not introduce systematic offsets in the GW amplitude ratio; this directly affects the load-bearing claim that radio supplies the reference measurement.

    Authors: The seven-parameter radio timing Fisher matrix is constructed in the standard manner for pulsar timing arrays, with the inclination parameter (via s) included among the fitted parameters and the reported uncertainties obtained after marginalization over all covariances. The LISA analysis employs the sky-averaged sensitivity curve that already folds in the Galactic foreground as an additive noise component; no additional systematic offset is introduced into the amplitude ratio beyond the statistical uncertainty captured by the Fisher matrix. The abstract is necessarily concise, but the full manuscript details the matrix construction and noise model. We will revise the abstract to include a short parenthetical reference to the marginalization procedure and the foreground treatment, thereby making the validation more visible at the abstract level while leaving the quantitative results unchanged. revision: yes

Circularity Check

0 steps flagged

No circularity: residual defined between independent messengers; GR prediction is the testable hypothesis

full rationale

The paper defines ε(t0) explicitly as the difference s_Shapiro(t0) − s_GW(t0) between two separately extracted quantities (Shapiro-delay shape parameter from radio timing and tensor amplitude ratio from GW polarimetry). GR is asserted to predict ε = 0 precisely because the test checks whether the two messengers yield consistent geometry; this is the hypothesis under test rather than a derived result. Precision estimates rely on a seven-parameter radio Fisher matrix and a sky-averaged LISA curve, constructed independently with no parameter fitted to one dataset and then relabeled as a prediction on the other. No self-citations, uniqueness theorems, or ansatze appear in the text to support the central claim. The construction is therefore a direct consistency test between distinct observables and remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The proposal rests on standard radio timing and GW analysis assumptions plus the conservative baseline values σ_TOA=1 μs and N_eff=10^4; no new entities are introduced.

free parameters (2)
  • σ_TOA
    Conservative radio timing uncertainty of 1 μs chosen to summarize multiple noise sources; used to set the scale of the Fisher matrix.
  • N_eff
    Effective number of observations set to 10^4 to account for cadence losses and other effects.
axioms (2)
  • domain assumption Shapiro delay shape parameter directly yields orbital inclination without additional relativistic corrections beyond GR
    Invoked when equating the radio measurement to the geometric inclination used in the GW comparison.
  • domain assumption Tensor polarization amplitude ratio directly yields the same orbital inclination
    Central to defining the common geometric quantity s_GW.

pith-pipeline@v0.9.1-grok · 5871 in / 1334 out tokens · 40193 ms · 2026-06-26T03:44:16.042173+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

33 extracted references · 16 linked inside Pith

  1. [1]

    Damour and J

    T. Damour and J. H. Taylor, Strong-field tests of rela- tivistic gravity and binary pulsars, Phys. Rev. D45, 1840 (1992)

    1992

  2. [2]

    I. H. Stairs, Testing general relativity with pulsar timing, Living Reviews in Relativity6, 5 (2003)

    2003

  3. [3]

    P. C. C. Freire, N. Wex, G. Esposito-Far` ese, J. P. W. Ver- biest, M. Bailes, B. A. Jacoby, M. Kramer, I. H. Stairs, J. Antoniadis, and G. H. Janssen, The relativistic pulsar– white dwarf binary psr j1738+0333. ii. the most stringent test of scalar-tensor gravity, Mon. Not. R. Astron. Soc. 423, 3328 (2012), arXiv:1205.1450

    Pith/arXiv arXiv 2012

  4. [4]

    Kramer, I

    M. Kramer, I. H. Stairs, R. N. Manchester, N. Wex, A. T. Deller, W. A. Coles, M. Ali, M. Burgay, F. Camilo, I. Cognard,et al., Strong-field gravity tests with the double pulsar, Physical Review X11, 041050 (2021), arXiv:2112.06795

    arXiv 2021

  5. [5]

    C. M. Will, The confrontation between general relativity and experiment, Living Rev. Relativ.21, 3 (2018)

    2018

  6. [6]

    M. Y. M. Lau, I. Mandel, A. Vigna-G´ omez, C. J. Neijssel, S. Stevenson, and A. Sesana, Detecting double neutron stars with LISA, Mon. Not. R. Astron. Soc.492, 3061 (2020), arXiv:1910.12422

    arXiv 2020

  7. [7]

    Kyutoku, Y

    K. Kyutoku, Y. Nishino, and N. Seto, How to de- tect the shortest-period binary pulsars in the era of LISA, Mon. Not. Roy. Astron. Soc.483, 2615 (2019), arXiv:1812.02177 [astro-ph.HE]

    Pith/arXiv arXiv 2019

  8. [8]

    Amaro-Seoaneet al., Laser interferometer space antenna, arXiv e-prints , arXiv:1702.00786 (2017), arXiv:1702.00786

    P. Amaro-Seoaneet al., Laser interferometer space antenna, arXiv e-prints , arXiv:1702.00786 (2017), arXiv:1702.00786

    Pith/arXiv arXiv 2017

  9. [9]

    Cutler, Angular resolution of the lisa gravitational wave detector, Phys

    C. Cutler, Angular resolution of the lisa gravitational wave detector, Phys. Rev. D57, 7089 (1998), arXiv:gr- qc/9703068

    arXiv 1998

  10. [10]

    T. A. Prince, M. Tinto, S. L. Larson, and J. W. Arm- strong, The lisa optimal sensitivity, Phys. Rev. D66, 122002 (2002), arXiv:gr-qc/0209039

    Pith/arXiv arXiv 2002

  11. [11]

    N. J. Cornish and L. J. Rubbo, The lisa response function, Phys. Rev. D67, 022001 (2003), arXiv:gr- qc/0209011

    arXiv 2003

  12. [12]

    Kupfer, V

    T. Kupfer, V. Korol, S. Shah, G. Nelemans, T. R. Marsh, G. Ramsay, P. J. Groot, and D. T. H. Steeghs, Lisa ver- ification binaries with updated distances from gaia data release 2, Mon. Not. R. Astron. Soc.480, 302 (2018), arXiv:1805.00482

    Pith/arXiv arXiv 2018

  13. [13]

    Babak, M

    S. Babak, M. Hewitson, and A. Petiteau, Lisa sensitivity and signal extraction, arXiv e-prints , arXiv:2108.01167 (2021), arXiv:2108.01167

    arXiv 2021

  14. [14]

    Yunes and X

    N. Yunes and X. Siemens, Gravitational-wave tests of general relativity with ground-based detectors and pulsar-timing arrays, Living Reviews in Relativity16, 9 (2013), arXiv:1304.3473

    Pith/arXiv arXiv 2013

  15. [15]

    D. M. Eardley, D. L. Lee, and A. P. Lightman, Gravitational-wave observations as a tool for testing rel- ativistic gravity, Phys. Rev. D8, 3308 (1973)

    1973

  16. [16]

    Nishizawa, A

    A. Nishizawa, A. Taruya, K. Hayama, S. Kawamura, and M.-a. Sakagami, Probing non-tensorial polarizations of stochastic gravitational-wave backgrounds with ground- based laser interferometers, Phys. Rev. D79, 082002 (2009), arXiv:0903.0528

    Pith/arXiv arXiv 2009

  17. [17]

    Damour and G

    T. Damour and G. Esposito-Far` ese, Tensor-multi-scalar theories of gravitation, Class. Quantum Grav.9, 2093 (1992)

    2093

  18. [18]

    Damour and G

    T. Damour and G. Esposito-Far` ese, Tensor-scalar gravity and binary-pulsar experiments, Phys. Rev. D54, 1474 (1996), arXiv:gr-qc/9602056. 14

    Pith/arXiv arXiv 1996

  19. [19]

    Janssen, G

    G. Janssen, G. Hobbs, M. McLaughlin, C. Bassa, A. Deller, M. Kramer, K. Lee, C. M. F. Mingarelli, P. Rosado, S. Sanidas,et al., Gravitational wave as- tronomy with the SKA, inAdvancing Astrophysics with the Square Kilometre Array (AASKA14)(2015) p. 37, arXiv:1501.00127

    Pith/arXiv arXiv 2015

  20. [20]

    J. M. Cordes and R. M. Shannon, A measurement model for precision pulsar timing, arXiv e-prints , arXiv:1010.3785 (2010), arXiv:1010.3785

    Pith/arXiv arXiv 2010

  21. [21]

    X. P. You, G. B. Hobbs, W. A. Coles, R. N. Manch- ester, R. Edwards, M. Bailes, J. Sarkissian, J. P. W. Ver- biest, W. van Straten, A. Hotan,et al., Dispersion mea- sure variations and their effect on precision pulsar timing, Mon. Not. R. Astron. Soc.378, 493 (2007), arXiv:astro- ph/0702366

    arXiv 2007

  22. [22]

    Goncharov, D

    B. Goncharov, D. J. Reardon, R. M. Shannon, X.-J. Zhu, E. Thrane, M. Bailes, N. D. R. Bhat, S. Dai, G. Hobbs, M. Kerr,et al., Identifying and mitigating noise sources in precision pulsar timing data sets, Mon. Not. R. Astron. Soc.502, 478 (2021), arXiv:2010.06109

    arXiv 2021

  23. [23]

    R. T. Edwards, G. B. Hobbs, and R. N. Manchester, Tempo2, a new pulsar-timing package. ii. the timing model and precision estimates, Mon. Not. R. Astron. Soc. 372, 1549 (2006), arXiv:astro-ph/0607664

    Pith/arXiv arXiv 2006

  24. [24]

    G. B. Hobbs, R. T. Edwards, and R. N. Manchester, Tempo2, a new pulsar-timing package. i. an overview, Mon. Not. R. Astron. Soc.369, 655 (2006), arXiv:astro- ph/0603381

    arXiv 2006

  25. [25]

    van Straten and M

    W. van Straten and M. Bailes, Dspsr: Digital signal processing software for pulsar astronomy, Publications of the Astronomical Society of Australia28, 1 (2011), arXiv:1008.3973

    Pith/arXiv arXiv 2011

  26. [26]

    Lentati, R

    L. Lentati, R. M. Shannon, W. A. Coles, J. P. W. Ver- biest, R. van Haasteren, J. A. Ellis, R. N. Caballero, R. N. Manchester, Z. Arzoumanian, S. Babak,et al., From spin noise to systematics: stochastic processes in the first in- ternational pulsar timing array data release, Mon. Not. R. Astron. Soc.458, 2161 (2016), arXiv:1602.05570

    Pith/arXiv arXiv 2016

  27. [27]

    Robson, N

    T. Robson, N. J. Cornish, and C. Liu, The construction and use of lisa sensitivity curves, Class. Quantum Grav. 36, 105011 (2019)

    2019

  28. [28]

    Takahashi and N

    R. Takahashi and N. Seto, Parameter estimation for galactic binaries by the laser interferometer space an- tenna, Astrophys. J.575, 1030 (2002)

    2002

  29. [29]

    N. J. Cornish and S. L. Larson, Lisa data analysis: Source identification and subtraction, Phys. Rev. D67, 103001 (2003), arXiv:gr-qc/0301002

    Pith/arXiv arXiv 2003

  30. [30]

    M. R. Adams, N. J. Cornish, and T. B. Littenberg, Astro- physical model selection in gravitational-wave astronomy, Astrophys. J.791, 76 (2014), arXiv:1404.1546

    Pith/arXiv arXiv 2014

  31. [31]

    B. M. Barker and R. F. O’Connell, Gravitational two-body problem with arbitrary masses, spins, and quadrupole moments, Phys. Rev. D12, 329 (1975)

    1975

  32. [32]

    Wex and S

    N. Wex and S. M. Kopeikin, Frame dragging and other precessional effects in black hole-pulsar binaries, Astro- phys. J.514, 388 (1999)

    1999

  33. [33]

    Mirshekari, N

    S. Mirshekari, N. Yunes, and C. M. Will, Constrain- ing generic lorentz violation and the speed of the gravi- ton with gravitational waves, Phys. Rev. D85, 024041 (2012), arXiv:1110.2720

    Pith/arXiv arXiv 2012