A Realistic Pulsar -- Supermassive Black Hole Timing Model

Lijing Shao; Zexin Hu; Ziming Wang

REVIEW 2 major objections 4 minor 2 cited by

A realistic 2PN pulsar–SMBH timing model forecasts sub-percent measurements of Sgr A*'s mass, spin, and quadrupole for tight orbits, while showing proper motion will not break the spin degeneracy.

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-02 21:35 UTC pith:K53HFHO4

load-bearing objection A careful PN timing model for pulsars around Sgr A* with genuine new results on proper motion and aberration, but the <1% precision forecast ignores a ~3 ms 3PN effect and should be read as an upper bound. the 2 major comments →

arxiv 2602.19546 v2 pith:K53HFHO4 submitted 2026-02-23 astro-ph.HE gr-qc

A Realistic Pulsar -- Supermassive Black Hole Timing Model

Zexin Hu , Ziming Wang , Lijing Shao This is my paper

classification astro-ph.HE gr-qc

keywords pulsar timingsupermassive black holesSagittarius A*post-Newtonian approximationno-hair theoremred noiseproper motionparameter estimation

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.

The paper builds a realistic timing model for a pulsar orbiting the supermassive black hole at the Galactic Center, accurate to second post-Newtonian order. It includes next-to-leading-order Shapiro delay, frame dragging, aberration, and — for the first time — the proper motion of the black hole. The authors forecast that with 1 ms timing from future radio telescopes, a pulsar with an orbital period under about half a year would pin down the black hole's mass, spin, and quadrupole moment all to better than 1 percent, enabling a test of the no-hair theorem. They also find that the black hole's proper motion does not efficiently break the known degeneracy in spin measurement for these tight orbits, contrary to earlier hopes, and that red noise must be modeled jointly with the timing parameters to avoid biased results.

Core claim

The central claim is that Equations (21), (27), (33) and (34) together with the 2PN equations of motion (11) form a complete, self-consistent timing model for pulsar–SMBH systems ready for real data analysis. At 1 ms timing precision, the model predicts sub-percent measurements of Sgr A*'s mass, spin, and quadrupole moment for pulsars with orbital periods below about 0.5 yr. The paper's first-time inclusion of proper motion shows that the longitude of the ascending node remains poorly constrained for tight orbits, so the proper motion cannot efficiently break the leading-order spin degeneracy. The model also incorporates next-to-leading-order aberration delays, which for the first time allow

What carries the argument

The central object is the numerical timing model built from the 2PN equations of motion of the pulsar in the harmonic-coordinate spacetime of the SMBH, including leading spin-orbit and quadrupole terms, and the corresponding series of time delays: Rømer, 1PN and 2PN Shapiro, frame-dragging, Einstein, and aberration delays. The key mechanism is the inverse timing model, which integrates the equations of motion in arrival time rather than coordinate time, making the model computationally efficient for parameter estimation. The paper also uses a Fisher-matrix formalism with analytically-computed derivatives to forecast measurement precisions.

Load-bearing premise

The load-bearing premise is that environmental perturbations around Sgr A*—stars, dark matter, and the interstellar medium—produce timing effects smaller than the 1 ms precision over a 5-year span, so that the vacuum 2PN spacetime terms dominate; this is justified only by order-of-magnitude estimates and left unmodeled.

What would settle it

Fit five years of 1 ms TOAs from a real pulsar with orbital period about 0.5 yr around Sgr A* with this model, first allowing environmental terms (stellar encounters, dark matter) to be included; if these terms are required for a white-noise residual, the sub-percent forecasts and the degeneracy conclusions do not survive in that regime.

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

Share X Bluesky LinkedIn Reddit HN

If this is right

With a pulsar of orbital period below about 0.5 yr and 1 ms timing, the SMBH mass, spin, and quadrupole can each be measured to better than 1%, providing a quantitative basis for a no-hair theorem test.
The proper motion of Sgr A* will not efficiently break the spin degeneracy for tight orbits; real spin measurements will need multiple pulsars or complementary observations.
The timing model is ready to be applied to real data analysis of future Galactic Center pulsars, including all delays relevant at the 2PN level.
Red noise, if modeled jointly with timing parameters, leads only to a mild increase in parameter uncertainties for pulsars with orbital periods less than about a tenth of the observing span; ignoring it significantly underestimates uncertainties.
The next-to-leading-order aberration delay is potentially detectable, allowing the pulsar's spin axis direction to be constrained to about 1 rad for pulsars with orbital periods below about 0.5 yr.

Where Pith is reading between the lines

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

The model's numerical structure means additional physical effects—a stellar-mass perturber, a dark-matter spike, or modified gravity—can be inserted as extra terms in the equations of motion and the same inverse-timing machinery reused; the authors leave this to future work.
The weak constraints on the ascending node for tight orbits imply that a single close pulsar may not suffice for a full three-dimensional spin vector; combining two pulsars or adding astrometry of the SMBH may be required to lift the degeneracy.
Since red noise with a shallower spectrum (e.g., from unmodeled dispersion-measure variations) may leak more power into the orbital frequency band, the forecast sub-percent precision is likely optimistic for pulsars embedded in a strongly scattering Galactic Center medium.
The detectability of 3PN Shapiro and higher-order terms for extreme orbits (small orbital period, large eccentricity, near edge-on) suggests that the same model, extended by one order, could probe strong-field gravity deeper than GR's 2PN expansion.

Editorial analysis

A structured set of objections, weighed in public.

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

Desk Editor's Note

A careful PN timing model for pulsars around Sgr A* with genuine new results on proper motion and aberration, but the <1% precision forecast ignores a ~3 ms 3PN effect and should be read as an upper bound.

read the letter

Here's my take. The paper builds a 2PN numerical timing model for a pulsar around Sgr A* and does it more carefully than previous work. The new pieces are the inclusion of Sgr A* proper motion and the next-to-leading-order aberration delay. The proper-motion result is actually useful: timing-only measurements of proper motion will be far worse than VLBA, and the proper motion does not efficiently break the spin degeneracy for tight orbits. That is a concrete, reproducible negative result. The red-noise analysis is also sensible and shows that if you know the noise, the loss is mild; if you don't, you can be badly misled.

The main problem is the headline forecast. The paper claims <1% precision on M, spin, and quadrupole for P_b < 0.5 yr. But Table 1 lists the neglected 'Higher PN' term as 3e-3 s for the fiducial system—three times the 1 ms timing precision—and the estimate grows to tens of ms for tighter orbits. The paper's justification—that this is degenerate with environmental perturbations—doesn't make it safe. It means both are unmodeled. An unmodeled secular signal at the millisecond level will be absorbed partially by the spin and quadrupole parameters, and the bias can exceed the quoted statistical errors. The forecast is an optimistic upper bound, not a reliable prediction for real data.

Does this kill the paper? No. The model and the qualitative results stand. The proper-motion conclusion is about information content, and it doesn't depend on the 3PN term being included. The red-noise procedure is orthogonal. But the paper should either include the 3PN periastron term or do an explicit bias budget. As written, the abstract's '<1%' claim is too strong.

The numerical implementation isn't released, which is a minor issue but worth noting for reproducibility. The appendices show convergence at the microsecond level, so the codes clearly work.

Who should read it: anyone building pulsar-SMBH timing analyses for SKA-era GC observations. It deserves a serious referee, but with a request to address the 3PN bias. I'd cite it for the proper-motion and aberration work, not for the precision forecast.

Referee Report

2 major / 4 minor

Summary. The paper constructs a numerical pulsar timing model for a pulsar orbiting Sgr A*. The model combines post-Newtonian (2PN) equations of motion for the pulsar with R\"omer, 1PN/2PN Shapiro, frame-dragging, Einstein, aberration, and lensing delays, and it incorporates the measured proper motion of Sgr A*. An efficient inverse timing model is described and validated with convergence tests. Using Fisher-matrix forecasts for a fiducial 5-yr, weekly-cadence, 1-ms observation, the paper reports that for P_b < 0.5 yr the SMBH mass, spin, and quadrupole moment can be measured to better than 1%. It also concludes that the known proper motion will not efficiently break the spin degeneracy for tight orbits, and it analyzes the impact of timing red noise, recommending a Bayesian treatment rather than an effective white-noise assumption.

Significance. If the central forecast is correct, this is a valuable and timely tool for SKA-era searches for pulsars around Sgr A*. The manuscript is carefully built on standard PN expressions, includes explicit convergence tests (Appendices B and C), and introduces a variational method for Fisher derivatives that avoids finite-difference noise in poorly constrained angular parameters. The proper-motion analysis is a new quantitative result, and the red-noise section goes beyond earlier pulsar-SMBH timing studies. The paper is not circular: inputs are external measurements, standard equations, and clearly stated fiducial choices. However, the headline <1% precision forecast is conditional on the 2PN truncation and on neglecting environmental perturbations, and the manuscript's own estimates show that one omitted secular effect is at the millisecond level for the fiducial system. That issue needs to be addressed before the forecast can be accepted as a reliable guide for real data.

major comments (2)

[Section 2.2, Eq. (13); Table 1] The neglected 3PN secular periastron advance is a load-bearing omission. Eq. (13) gives ~1 ms (T_obs/5 yr)(P_b/1 yr)^{-7/3} sin i; for the fiducial system (P_b=0.5 yr, i=pi/5) this is ~3 ms, i.e., three times the assumed 1 ms TOA precision, and Table 1 lists the associated 'Higher PN' delay as 3e-3 s. The argument that this term is 'largely degenerate with environmental perturbations' is not sufficient: an unmodeled secular signal at the millisecond level will be partly absorbed by the fitted parameters, especially chi and q, and can bias the recovered values by more than the quoted statistical Fisher errors from Section 3.2. The <1% forecast for P_b<0.5 yr is therefore not yet supported. Please add a quantitative bias analysis, e.g., inject the 3PN/environmental secular terms into mock data and recover with the 2PN model, or explicitly restrict the forecast to the regime where Eq. (13)
[Section 2.2; Table 1; Section 3.2] The word 'realistic' in the title/abstract is stronger than what is demonstrated while environmental perturbations are unmodeled. The text acknowledges (citing Merritt et al. 2010; Hu et al. 2023) that stellar-mass perturbers and a dark-matter spike can affect the orbit, but Table 1 contains no entry for them and no mock-injection test is performed. Because such perturbations can produce secular orbital changes similar in character to the PN effects being measured, the Section 3.2 Fisher forecasts are conditional on an isolated SMBH. The paper explicitly defers environmental modeling, which is acceptable for a first model, but the 'realistic' claim and the data-analysis-readiness statement should be qualified accordingly.

minor comments (4)

[Section 2.3, Eq. (38)] The text states the apparent motion of Sgr A* is about -6.4 mas/yr along the Galactic plane and -0.22 mas/yr toward the North Galactic Pole, but Eq. (38) gives mu_alpha=-3.2 and mu_delta=-5.6 mas/yr. Please clarify the coordinate conversion between these two sets.
[Section 2.2, after Eq. (13)] The sentence 'This is different when considering secular effects, as they can have distinctive signatures and be separated' appears to say the opposite of the intended argument. Rephrase to avoid confusion.
[Section 2.2, Eq. (18)] 'Geodesics precession' should be 'Geodetic precession'.
[Appendix C or Section 2.4] Minor typo: 'inverse the numerical relation' should be 'invert the numerical relation'.

Circularity Check

0 steps flagged

No significant circularity: the timing model is assembled from standard PN equations plus external measurements; the precision forecasts are Fisher-matrix projections, not fitted predictions.

full rationale

The derivation chain is self-contained rather than circular. The equations of motion (Eq. 11) are the standard 1PN/2PN/SO/quadrupole accelerations cited to Barker & O'Connell (1975), Damour & Schaefer (1988), and Wex (1995); the propagation delays in Eqs. (2), (5)-(7), and (10) come from Shapiro (1964), Klioner & Zschocke (2010), and Wex & Kopeikin (1999). The SMBH proper motion is taken from the external VLBA measurement of Reid & Brunthaler (2020), not derived from the timing model. The new claims—proper motion of Sgr A* does not efficiently break the spin degeneracy, and aberration delays can constrain the pulsar spin orientation—are computed results of the model, not inputs. Self-citations to Hu et al. (2023) provide the numerical integration scheme and earlier forecast checks, but they do not force the new conclusions; the paper explicitly verifies consistency with previous independent frameworks (Liu et al. 2012; Psaltis et al. 2016; Zhang & Saha 2017). No parameter is fitted to a subset of data and then renamed a prediction; the forecasted precisions are Fisher-matrix covariances on noiseless simulated data. The paper's own limitation statement about the neglected 3PN secular periastron advance (Eq. 13; Table 1, 'Higher PN' row) is a genuine accuracy/bias concern for the <1% forecast, but it is an omitted-effect robustness issue, not a circular step: the effect is quantified, acknowledged, and not relabeled as part of the derivation. The manuscript explicitly leaves environmental modeling to future work, which further confirms that no input is being recycled as an output.

Axiom & Free-Parameter Ledger

7 free parameters · 7 axioms · 0 invented entities

The model is built from established PN equations and standard timing formulas. No parameters are fitted to real data; the central forecasts depend on hand-chosen fiducial parameters, the external VLBA proper-motion measurement, and assumed timing/noise properties. No new physical entities are introduced.

free parameters (7)

Fiducial SMBH spin magnitude χ = 0.6
Hand-chosen conservative value (Sec 2.2, Eq 17); Fisher forecasts of spin precision and degeneracy geometry depend on it.
Fiducial spin orientation (λ, η) = λ=π/6, η=5π/9
Hand-chosen to avoid special orientations that might render the estimation degenerate; spin-degeneracy conclusions may shift for other orientations.
Fiducial quadrupole q = -0.36
Set to -χ² for χ=0.6, consistent with GR expectation; treated as independent in the no-hair test.
Fiducial orbital parameters (P_b, e, i, ω, Ω, f0) = P_b=0.5 yr, e=0.8, i=π/5, ω=5π/7, Ω=0, f0=-3π/4
Hand-chosen fiducial system; the forecast curves for σ_M, σ_χ, σ_q and the cadence-related spikes depend on these choices.
Pulsar rotation parameters (ν, νdot, λp, ηp) = ν=1 Hz, νdot=1e-15 s^-2, λp=π/5, ηp=0
Normal-pulsar assumption; aberration amplitude and the claimed ~1 rad precision on pulsar spin direction scale with these.
Observation setup = weekly cadence, T_obs=5 yr, σ_TOA=1 ms
Assumed future SKA-like timing precision from Liu et al. (2012); all uncertainty forecasts scale linearly with σ_TOA, and weekly cadence creates resonance spikes.
Red noise model parameters (A, α, f_c) = A=1e-13 or 1e-15 yr^3, α=5, f_c=1/50 yr^-1
Chosen to represent strong and weak intrinsic red noise for normal pulsars; the recovery biases and uncertainty increases depend on these values.

axioms (7)

domain assumption Post-Newtonian equations of motion in Eqs (11)-(12) correctly describe 2PN test-particle dynamics in harmonic coordinates.
Invoked in Sec 2.2; the model neglects mass ratio, SSC ambiguity, 3PN terms, and radiation reaction based on order-of-magnitude estimates.
domain assumption Light-propagation time delays in Eqs (1)-(10) from Klioner & Zschocke (2010) and Wex & Kopeikin (1999) are correct to 2PN in the SMBH spacetime.
Used in Sec 2.1; ISM dispersion and scattering are ignored, and the Sgr A*-SSB distance is treated as effectively infinite.
domain assumption The SMBH spacetime is axisymmetric and the no-hair constraint q=-χ² is not imposed.
Stated in Sec 2.2; this permits an independent quadrupole measurement for testing the no-hair theorem.
domain assumption Environmental perturbations (stellar mass objects, dark matter spike, ISM) are negligible at the assumed timing precision.
Sec 2.2 and Sec 4 leave environmental modeling to future study; if these produce secular effects comparable to the 2PN spin/quadrupole signals, the forecasts would change.
domain assumption Proper motion of Sgr A* is known from VLBA and is applied as a secular rotation of the line-of-sight vector only.
Sec 2.3, Eq (36); PN-proper motion cross terms are argued to be below observational precision.
domain assumption Gaussian white timing noise and design-matrix linearization are valid for the Fisher and red-noise analyses.
Sec 3.1 and Sec 4.2; this is an approximation valid when noise and model inaccuracy are smaller than the pulsar rotation period.
domain assumption Red noise follows a power-law spectrum with cutoff f_c=1/50 yr^-1 and α=5 in the simulated examples.
Eqs (48)-(49); motivated by normal pulsar studies, but other red-noise shapes (e.g., DM variation) could affect shorter orbits.

pith-pipeline@v1.3.0-alltime-deepseek · 35309 in / 14295 out tokens · 117735 ms · 2026-08-02T21:35:04.757674+00:00 · methodology

0 comments

Cite this review

Pith. "Pith review of A Realistic Pulsar -- Supermassive Black Hole Timing Model." pith.science (2026). https://pith.science/paper/K53HFHO4

@misc{pith2026260219546,
  author       = {Pith},
  title        = {Pith review of: A Realistic Pulsar -- Supermassive Black Hole Timing Model},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/K53HFHO4}},
  note         = {Machine review of arXiv:2602.19546}
}

Share X Bluesky LinkedIn Reddit HN

read the original abstract

Timing observation of pulsars orbiting around a supermassive black hole (SMBH) can measure the spacetime around the SMBH to a high precision and thus be a novel probe of the gravity theory. Future high-frequency surveys of the Galactic Centre (GC) region to be performed by the next-generation radio telescopes, such as the SKA, may discover pulsars that orbit around Sagittarius A* (Sgr A*), the SMBH dwelling in our GC. In this paper, we present a realistic pulsar-SMBH timing model based on the post-Newtonian equations of motion of the pulsar. Considering the expected timing precision in the future, we take into account several next-to-leading order light propagation time delays in the timing model. For the first time, we include the effects of proper motion of Sgr A*, which were expected to break the spin measurement degeneracy. We forecast the measurement precision of various parameters of Sgr A*, and discuss the data analysis procedure in the presence of red noise, which can be strong if the pulsar is a normal pulsar. The realistic timing model constructed in this study will serve as a useful tool in future searching and timing of pulsar-SMBH systems in the GC.

Figures

Figures reproduced from arXiv: 2602.19546 by Lijing Shao, Zexin Hu, Ziming Wang.

**Figure 1.** Figure 1: Illustration of the pulsar orbit with system parameters given in Equation (17). The axes are in the unit of AU. For a pulsar in a pulsar-SMBH system with an orbital period Pb = 1 yr and orbital eccentricity e = 0.8, one has P˙ b ∼ −5 × 10−12. This leads to a cumulated time delay at the order of πP˙ b T 2 P 2 b a sin i c ∼ 30 µs Tobs 5 yr 2 Pb 1 yr−3 sin i . (16) Therefore, for a total observation ti… view at source ↗

**Figure 2.** Figure 2: Various time delays described in the text for the fiducial pulsar-SMBH system. A constant term is removed from the 1PN Shapiro delay with respect to Equation (2). every single pulse, we, in fact, do not know the pulsar rotation number N corresponding to the observed TOAs. Thus, the timing residual, which represents the difference between the model prediction and observation, is defined to be (Damour & De… view at source ↗

**Figure 3.** Figure 3: The fractional precision of SMBH parameters as functions of the pulsar’s (a) orbital period and (b) orbital eccentricity. 0.1 0.2 0.3 0.4 0.6 0.8 1.0 2.0 3.0 Pb [yr] 100 101 102 103 σ σΩ [rad] σµα [mas/yr] σµδ [mas/yr] 0.1 0.2 0.3 0.4 0.6 0.8 1.0 2.0 3.0 Pb [yr] 10−6 10−5 10−4 10−3 10−2 10−1 100 101 Fractional Precision σM/M w/o Ω σχ/χ w/o Ω σq/|q| w/o Ω σM/M w/ Ω σχ/χ w/ Ω σq/|q| w/ Ω (a) (b) [PITH_FULL_… view at source ↗

**Figure 4.** Figure 4: (a) The measurement precision of the proper motion parameters µα and µδ, and the longitude of the ascending node Ω. Due to the degeneracy caused by the rotation symmetry, for estimating the measurability of the proper motion parameters, we fix Ω, and vice versa, we fix the proper motion parameters when estimating the measurement precision of Ω. (b) Comparison of parameter estimation results with or without… view at source ↗

**Figure 5.** Figure 5: Expected measurement precision of λp and ηp as functions of the pulsar orbital period Pb. In Section 2.3, we discussed the aberration delays, including ∆A1, ∆A2 and ∆L, which come from the rotation origin of the pulsar’s “lighthouse” pulse. Though discussed in Damour & Deruelle (1986), the leadingorder aberration effect will be absorbed in a redefinition of various timing parameters, we have to include… view at source ↗

**Figure 6.** Figure 6: (a) Two noise realizations with different red noise amplitude A. The power-law index is chosen to be α = 5 for both panels. The white-noise components are shown by the blue dashed line and we have used σTOA = 1 ms. f2 denotes the timing residuals after a second-order polynomial fitting that subtracts the pulsar rotation’s contribution via N0, ν, and ˙ν. ffull denotes the timing residuals after a fitting of… view at source ↗

**Figure 7.** Figure 7: The evidence for four examples as functions of the Fourier bases number k. As studied by Lentati et al. (2013) and van Haasteren & Levin (2013), for estimating the red noise parameter Ξ, one can analytically marginalize over the Fourier amplitudes ϵ and the timing parameters δΘ when under the assumption of a flat prior. Integrating over ϵ gives P(δΘ,Ξ|δt) ∝ P(Ξ) exp h − 1 2 [PITH_FULL_IMAGE:figures/full_… view at source ↗

**Figure 8.** Figure 8: Posteriors of red-noise parameter for four cases with different pulsar orbital periods and red noise amplitudes. The number of Fourier bases is k = 6 for all these figures. The true values of the parameters are marked by the orange lines and squares. The dashed lines show the 1-σ intervals of the 1-D marginalized distributions, while the contours show the 68% and 90% credible regions in the 2-D parameter s… view at source ↗

**Figure 9.** Figure 9: (a) The absolute difference between the averaged values of the best-fit noise parameters and true parameters, as well as the standard deviation of the best-fit parameters for Nr = 10000 noise realizations. In this figure we have used A˜13 = 8×10−21 and ¯α = 5. (b) The fractional precision of SMBH parameters as functions of the pulsar orbital period in the presence of red noise. We have used A = 10−13 yr3 a… view at source ↗

**Figure 10.** Figure 10: The recovered posterior for the Pb = 0.5 yr and A = 10−13 yr3 case shown before. We only show the posteriors of the noise parameters, and M, χ, q of the SMBH. No clear correlation between noise parameters and timing parameters is found. we still suggest to perform full noise analysis for pulsarSMBH systems when possible. 5. CONCLUSIONS In this work, we construct a realistic pulsar-SMBH timing model based… view at source ↗

**Figure 11.** Figure 11: (a) Convergence test of t TOA as a function of t. We choose four different desired relative errors, namely ϵrel = 10−10 , 10−11, 10−12, and 10−13. One can see a clear convergence among the results. The system parameters used for this figure are the same as the parameters used for [PITH_FULL_IMAGE:figures/full_fig_p024_11.png] view at source ↗

**Figure 12.** Figure 12: (a) Comparison of the parameter estimation results obtained with two different methods. The dashed lines employ the finite difference method, while the solid lines are based on solving the differential equations. One clearly sees that, when the pulsar’s orbital periods are large, which leads to a worse measurement precision, the numerical errors in the results are largely suppressed with the method we pro… view at source ↗

discussion (0)

Forward citations

Cited by 2 Pith papers

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

Probing an Intermediate-Mass Black Hole Companion of Sagittarius A* with Pulsar Timing
astro-ph.HE 2026-07 conditional novelty 6.0

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.
Granular mass perturbations on the pulsar - supermassive black hole system
astro-ph.HE 2026-06 unverdicted novelty 6.0

Simulations find granular stellar-mass black hole cusp induces 10-100 s timing residuals in pulsar-SMBH orbits, with periastron-only analysis plus frame-dragging improving spin precision by ~10x.

Reference graph

Works this paper leans on

6 extracted references · 5 linked inside Pith · cited by 2 Pith papers

[1]

G., et al

Abac, A. G., et al. 2025a, Astrophys. J. Lett., 995, L18, doi: 10.3847/2041-8213/ae0c06 —. 2025b. https://arxiv.org/abs/2508.18082 Abbate, F., et al. 2025, Open J. Astrophys., 8, 54252, doi: 10.33232/001c.154252 Abbott, B. P., et al. 2016a, Phys. Rev. Lett., 116, 061102, doi: 10.1103/PhysRevLett.116.061102 —. 2016b, Phys. Rev. Lett., 116, 221101, doi: 10....

Pith/arXiv arXiv 2041
[5]

V., & Kopeikin, S

https://arxiv.org/abs/2502.01093 Doroshenko, O. V., & Kopeikin, S. M. 1995, Mon. Not. Roy. Astron. Soc., 274, 1029, doi: 10.1093/mnras/274.4.1029 Eatough, R. P., et al. 2013, Nature, 501, 391, doi: 10.1038/nature12499 Ferdman, R. D., et al. 2013, Astrophys. J., 767, 85, doi: 10.1088/0004-637X/767/1/85 Freire, P. C. C., & Wex, N. 2024, Living Rev. Rel., 27...

Pith/arXiv arXiv 1995
[6]

https://arxiv.org/abs/2508.09931 Shapiro, I. I. 1964, Phys. Rev. Lett., 13, 789, doi: 10.1103/PhysRevLett.13.789 Shklovskii, I. S. 1970, Soviet Ast., 13, 562 Stairs, I. H., Thorsett, S. E., & Arzoumanian, Z. 2004, Phys. Rev. Lett., 93, 141101, doi: 10.1103/PhysRevLett.93.141101 Taylor, J. H., Fowler, L. A., & McCulloch, P. M. 1979, Nature, 277, 437, doi: ...

arXiv 1964
[793]

N., et al

https://arxiv.org/abs/1810.06623 PULSAR-SMBH TIMING MODEL27 Caballero, R. N., et al. 2016, Mon. Not. Roy. Astron. Soc., 457, 4421, doi: 10.1093/mnras/stw179 Carter, B. 1971, Phys. Rev. Lett., 26, 331, doi: 10.1103/PhysRevLett.26.331 Chen, S., et al. 2025, Astron. Astrophys., 699, A165, doi: 10.1051/0004-6361/202452550 Cie´ slik, A., Hackmann, E., & Mach, ...

Pith/arXiv arXiv 2016
[2002]

1986, Ann

https://arxiv.org/abs/astro-ph/0207156 Damour, T., & Deruelle, N. 1986, Ann. Inst. Henri Poincar´ e Phys. Th´ eor., 44, 263 Damour, T., & Schaefer, G. 1988, Nuovo Cim. B, 101, 127, doi: 10.1007/BF02828697 Damour, T., & Taylor, J. H. 1992, Phys. Rev. D, 45, 1840, doi: 10.1103/PhysRevD.45.1840 Della Monica, R., & de Martino, I

Pith/arXiv arXiv 1986
[2025]

S., Cordes, J

https://arxiv.org/abs/2501.03912 Deneva, J. S., Cordes, J. M., & Lazio, T. J. W. 2009, Astrophys. J. Lett., 702, L177, doi: 10.1088/0004-637X/702/2/L177 Desvignes, G., Eatough, R. P., Men, Y., et al. 2026, Astron. Astrophys., 706, A113, doi: 10.1051/0004-6361/202556381 Dexter, J., & O’Leary, R. M. 2014, Astrophys. J. Lett., 783, L7, doi: 10.1088/2041-8205...

Pith/arXiv arXiv 2009

[1] [1]

G., et al

Abac, A. G., et al. 2025a, Astrophys. J. Lett., 995, L18, doi: 10.3847/2041-8213/ae0c06 —. 2025b. https://arxiv.org/abs/2508.18082 Abbate, F., et al. 2025, Open J. Astrophys., 8, 54252, doi: 10.33232/001c.154252 Abbott, B. P., et al. 2016a, Phys. Rev. Lett., 116, 061102, doi: 10.1103/PhysRevLett.116.061102 —. 2016b, Phys. Rev. Lett., 116, 221101, doi: 10....

Pith/arXiv arXiv 2041

[2] [5]

V., & Kopeikin, S

https://arxiv.org/abs/2502.01093 Doroshenko, O. V., & Kopeikin, S. M. 1995, Mon. Not. Roy. Astron. Soc., 274, 1029, doi: 10.1093/mnras/274.4.1029 Eatough, R. P., et al. 2013, Nature, 501, 391, doi: 10.1038/nature12499 Ferdman, R. D., et al. 2013, Astrophys. J., 767, 85, doi: 10.1088/0004-637X/767/1/85 Freire, P. C. C., & Wex, N. 2024, Living Rev. Rel., 27...

Pith/arXiv arXiv 1995

[3] [6]

https://arxiv.org/abs/2508.09931 Shapiro, I. I. 1964, Phys. Rev. Lett., 13, 789, doi: 10.1103/PhysRevLett.13.789 Shklovskii, I. S. 1970, Soviet Ast., 13, 562 Stairs, I. H., Thorsett, S. E., & Arzoumanian, Z. 2004, Phys. Rev. Lett., 93, 141101, doi: 10.1103/PhysRevLett.93.141101 Taylor, J. H., Fowler, L. A., & McCulloch, P. M. 1979, Nature, 277, 437, doi: ...

arXiv 1964

[4] [793]

N., et al

https://arxiv.org/abs/1810.06623 PULSAR-SMBH TIMING MODEL27 Caballero, R. N., et al. 2016, Mon. Not. Roy. Astron. Soc., 457, 4421, doi: 10.1093/mnras/stw179 Carter, B. 1971, Phys. Rev. Lett., 26, 331, doi: 10.1103/PhysRevLett.26.331 Chen, S., et al. 2025, Astron. Astrophys., 699, A165, doi: 10.1051/0004-6361/202452550 Cie´ slik, A., Hackmann, E., & Mach, ...

Pith/arXiv arXiv 2016

[5] [2002]

1986, Ann

https://arxiv.org/abs/astro-ph/0207156 Damour, T., & Deruelle, N. 1986, Ann. Inst. Henri Poincar´ e Phys. Th´ eor., 44, 263 Damour, T., & Schaefer, G. 1988, Nuovo Cim. B, 101, 127, doi: 10.1007/BF02828697 Damour, T., & Taylor, J. H. 1992, Phys. Rev. D, 45, 1840, doi: 10.1103/PhysRevD.45.1840 Della Monica, R., & de Martino, I

Pith/arXiv arXiv 1986

[6] [2025]

S., Cordes, J

https://arxiv.org/abs/2501.03912 Deneva, J. S., Cordes, J. M., & Lazio, T. J. W. 2009, Astrophys. J. Lett., 702, L177, doi: 10.1088/0004-637X/702/2/L177 Desvignes, G., Eatough, R. P., Men, Y., et al. 2026, Astron. Astrophys., 706, A113, doi: 10.1051/0004-6361/202556381 Dexter, J., & O’Leary, R. M. 2014, Astrophys. J. Lett., 783, L7, doi: 10.1088/2041-8205...

Pith/arXiv arXiv 2009

A Realistic Pulsar -- Supermassive Black Hole Timing Model

Core claim §

What carries the argument

Load-bearing premise §

What would settle it §

If this is right

Where Pith is reading between the lines

discussion (0)

Core claim

Load-bearing premise

What would settle it