{"id":"c202db81-b78d-4475-b3da-c3fd0c4d41f3","arxiv_id":"2501.03912","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A fully relativistic (geodesic) pulsar timing model for the Galactic Center shows that 1PN-based timing formulas produce errors of up to seconds for tight orbits around Sgr A*.","lead":"Pulsars orbiting the supermassive black hole at the Galactic Center would be extremely precise probes of strong gravity, but this paper shows that current weak-field timing formulas fail badly in that regime. The authors propose a timing model that computes photon travel times by full geodesic integration, and find that 1PN approximations can miss by seconds, ruining phase connection.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The paper's 1PN-vs-geodesic comparison uses fixed parameters and a simplified delay formula, not a fitted timing code, so the 'current codes fail' claim is not yet established.","rationale":"Independent support: the numerical solver is validated against the Hackmann–Dhani exact solution to 10^-7 s, so the geodesic pipeline itself is credible. The load-bearing step is the extrapolation from 'Eq. (27) differs from exact for fixed parameters' to 'current codes fail in practice.' That extrapolation fails because (i) the delay model in real codes differs from Eq. (27) and (ii) parameter fitting can absorb model error. The paper even concedes that higher-order PN might alleviate the issue, which underscores that the 1PN version is not the only available approximation. The reader's weakest_assumption identifies the same gap. A mock TOA injection into TEMPO2/PINT is a cheap, decisive check. Until that is done, the conditional verdict is appropriate.","tokens_in":22217,"tokens_out":10290,"duration_ms":105121,"concrete_test":"Generate mock TOAs with the fully relativistic pipeline for a Table I orbit (e.g., Toy 2) with P=2 s and SKA-like 100 μs uncertainties, then fit them with TEMPO2 or PINT using the DD post-Keplerian model with the SMBH as the companion. If the best-fit residuals are white and within 100 μs, current timing codes can absorb the 1PN error; if residuals are phase-correlated or exceed P, the paper's claim is confirmed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section IV B compares the exact geodesic TOA against Eq. (27), a fixed-parameter sum of straight-line Rømer, Shapiro, and geometric delays. This is not what TEMPO2/PINT compute: the Damour–Deruelle model used by those codes includes post-Keplerian orbital dynamics, the Einstein delay γ sin u, aberration, and a Shapiro delay of different functional form. More importantly, a real timing fit adjusts all parameters to minimize residuals. The paper fixes the parameters at their input values; it never generates exact-geodesic TOAs and then fits them with a standard code. If the fitted mass, semi-major axis, or inclination shifts to absorb part of the ~2 s discrepancy, the post-fit residuals can remain well below the pulse period even though the delay formula is approximate. Figure 5's one-parameter-at-a-time sensitivity study is not a substitute for a full correlated least-squares fit. Therefore the headline inference that current timing codes lose phase connection is under-supported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a fully geodesic timing pipeline for future pulsars orbiting Sgr A*. For a given set of orbital and intrinsic parameters, it integrates the pulsar geodesic, solves the emitter-observer problem for each emitted photon, and inverts the TOA-versus-proper-time relation by spline interpolation to build a timing model. The authors compare this model with a fixed-parameter 1PN delay model composed of Rømer, Shapiro, and geometric delays (Eq. 27), reporting discrepancies of order seconds for a 100M circular orbit and much larger discrepancies for tighter and more eccentric orbits (Figure 6). They also compute residuals caused by one-at-a-time misestimations of parameters such as the pulsar period, spin-down rate, SMBH mass, semi-major axis, eccentricity, and inclination, and use these to argue for extreme future constraining power on the SMBH mass and orbital parameters.","tokens_in":22319,"tokens_out":6277,"duration_ms":67121,"significance":"If the comparison were made against the actual timing models used by TEMPO2/PINT and included a full parameter fit, this would be an important proof-of-concept for strong-field pulsar timing around Sgr A*. The numerical solver is a clear strength: it is validated against the independent Hackmann-Dhani exact solution for a circular-orbit case, and the framework is general enough to be extended beyond Schwarzschild to other spherically symmetric spacetimes. However, the paper's headline claim that 'current timing codes' fail is not yet established, because the 1PN comparison uses a simplified fixed-parameter delay sum rather than a fitted post-Keplerian model, and the sensitivity analysis is based on one-at-a-time parameter offsets rather than a correlated fit. These issues are load-bearing for the abstract and conclusions, so the paper needs a revision that either sharpens the claims or adds the missing comparison.","major_comments":[{"comment":"The comparison labeled 'current timing codes' tests a fixed-parameter sum ΔtPN = ΔtR + ΔtSh + Δtgeo, not the timing model implemented in TEMPO2/PINT. Standard codes use the Damour-Deruelle post-Keplerian model, which includes relativistic orbital dynamics, the Einstein delay γ sin u, aberration, and a different Shapiro delay parametrization, and they fit all parameters globally to minimize the residuals. A fixed-parameter comparison cannot determine whether a real timing code loses phase connection, because the best-fit mass, semi-major axis, inclination, and pulsar period can absorb part of the delay discrepancy. Please either reframe the claim as applying to the simplified 1PN model only, or generate geodesic TOAs and fit them with TEMPO2/PINT to test the 'current codes fail' conclusion directly.","section":"Section IV B, Eq. (27), Figure 6"},{"comment":"The claimed constraining power is derived from varying one parameter at a time while keeping all other parameters fixed, and no noise is included in the simulated residuals. In a real fit, parameter covariances and TOA uncertainties strongly affect the posterior, so the percentages reported in Table II should not be presented as projected measurement precisions. Please either run a full correlated fit (e.g., a Fisher-matrix or Markov-chain Monte Carlo forecast) or explicitly label these numbers as qualitative one-parameter sensitivities.","section":"Section IV A, Figure 5, Table II"},{"comment":"The text states that the exact Hackmann-Dhani solution is for a non-inclined circular orbit, but Figure 2 describes the toy model as having inclination i = 60°. Please clarify whether Eq. (20) applies to an inclined orbit; if it is restricted to zero inclination, the validation should be repeated for the inclined case or the caption corrected. This is the main external validation of the numerical pipeline, so the configuration must be unambiguous.","section":"Section II, Eq. (20), Figure 2"}],"minor_comments":[{"comment":"The sentence 'The results of this analysis are shown in Figure ??' contains an unresolved placeholder; the actual figure appears as Figure 6 and should be cited explicitly in the text.","section":"Section IV B"},{"comment":"There are several typographical errors, including 'Galctic' in the introduction to Section III A, 'Earth-bsaed' in the Conclusions, and 'Seciton' in the Conclusions; these should be corrected in a final proofreading pass.","section":"Throughout"},{"comment":"The inset axis label '5|Num-Ex| (10^-7 s)' is unclear; please label the quantity explicitly, for example '|Numerical - Exact| (10^-7 s)', and ensure the units are unambiguous.","section":"Figure 2"},{"comment":"The statement that Ndense roughly 50 times Ndata provides a robust reconstruction of TOA(τ) is not accompanied by a quantitative convergence test; the paper notes that increasing Ndense gives improvements below sensitivity, but showing the interpolation error as a function of Ndense would make the choice more defensible.","section":"Section III B"}],"recommendation":"major_revision","confidential_remarks":"The paper has a promising idea and a validated numerical core, but the abstract and conclusions currently overstate the case againt 'current timing codes' and the projected parameter constraints. I would be willing to reconsider after the authors either run a proper TEMPO2/PINT fit on geodesic TOAs or substantially soften the claims to refer specifically to simplified 1PN delay models. The Figure ?? placeholder and the inclined-orbit validation question should also be fixed promptly."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a useful proof-of-concept for building a fully geodesic pulsar timing model around Sgr A*, and the core numerical result—that simple 1PN delay formulas miss by seconds for tight orbits—is solid. The headline claim that current timing codes fail is not established by the comparison as written.\n\nWhat's genuinely new: they take their earlier emitter-observer solver [37] and turn it into a phase-connected timing pipeline, including interpolation of the TOA(τ) map and residual generation. The validation against the Hackmann–Dhani exact circular-orbit solution (Fig. 2, agreement at 1e-7 s) is the right kind of check and gives real weight to the 1PN-vs-geodesic discrepancy. The sensitivity maps in Fig. 5 and Table II are a reasonable qualitative tool for thinking about what such pulsars could measure; the order-of-magnitude claims about mass constraints are plausible if the pipeline is right.\n\nThe soft spots are real and concentrated in Section IV B. The 'current timing codes' comparison uses Eq. (27), a fixed-parameter sum of straight-line Rømer, Shapiro, and geometric delays. That is not what TEMPO2/PINT compute: they use Damour–Deruelle post-Keplerian orbital dynamics, Einstein delay, aberration, and a different Shapiro parametrization—and, critically, they fit the parameters to minimize residuals. The authors fix the parameters at input values and never generate geodesic TOAs and then fit them with a standard code. A real fit could absorb part of the ~2 s discrepancy by shifting mass or orbital elements, so the phase-connection loss is overclaimed. Also, the sensitivity analysis varies one parameter at a time, so the 'tremendous constraining power' phrasing is ahead of an actual covariance or Bayesian inference demonstration. No code or data is released, and the pipeline depends on the authors' own solver from [37], which is not public; that limits reproducibility but is not a fatal flaw given the external validation.\n\nThe paper knows it's a work in progress—the conclusion says so, which I take as honesty rather than an excuse.\n\nWho it's for: people planning GC pulsar searches or building timing codes for extreme-mass-ratio binaries. They will get a clear picture of the issues and a pipeline to build on. It deserves a serious referee, but the referee should ask for a realistic comparison to existing codes (or a fitted version of the 1PN model), a proper covariance/full-inference study, and public code/data before the strong claims are accepted.","headline":"Useful proof-of-concept for fully geodesic GC pulsar timing, but the claim that current codes fail is overreach: the 1PN comparison fixes parameters, never fits, and doesn't test real TEMPO2/PINT implementations.","tokens_in":22942,"tokens_out":2208,"would_cite":false,"duration_ms":21874,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C10","83C25","83C57"],"pacs":["04.20.-q","04.25.Nx","97.60.Gb"],"model":"deepseek-v4-flash","headline":"First-order post-Newtonian timing models fail by seconds for pulsars around supermassive black holes, so fully geodesic photon propagation is needed to keep phase connection.","keywords":["pulsar timing","Galactic Center","supermassive black hole","post-Newtonian approximation","photon propagation","arrival time residuals","Schwarzschild spacetime","gravitational lensing"],"falsifier":"Take the Toy 0 or S2-like orbit, generate noise-free TOAs from the geodesic pipeline, and fit them with the 1PN timing model of Eq. (27): a residual structure peaking near superior conjunction at the predicted 1-second-to-2-second level would confirm the paper's claim, whereas residuals below 100 microseconds across the full orbit would falsify it.","tokens_in":21862,"feed_emoji":"🕳️","tokens_out":5679,"duration_ms":56351,"temperature":0.7,"pith_summary":"The paper argues that standard pulsar timing models, which compute photon travel time at first post-Newtonian order and assume straight-line light paths, fail when the pulsar orbits a supermassive black hole like Sgr A* at the Galactic Center: the error in predicted pulse arrival time can reach about two seconds, far above the 100-microsecond accuracy goal of future telescopes. The authors propose replacing the approximate delay formulas with a fully relativistic computation that solves the photon geodesic numerically for each emitted pulse and inverts the resulting emission-time-to-arrival-time relation. They demonstrate the model on toy pulsar orbits of increasing compactness and show that misestimating parameters such as black hole mass, pulsar period, or orbital elements produces phase-dependent residuals detectable within months to a few years. If correct, future Galactic Center pulsar timing will require geodesic timing models rather than the post-Newtonian formulas now used in standard codes.","feed_headline":"Standard pulsar timing misses black-hole pulses by seconds","feed_subtitle":"For pulsars orbiting the Milky Way's central black hole, only full geodesic light paths keep arrival times in sync.","key_machinery":"The load-bearing object is the numerical solution of the emitter-observer problem for null geodesics in a static, spherically symmetric spacetime. For a chosen emission event, the impact parameter b is found by solving the angular integral that links the emitter and observer positions, and the coordinate travel time is the integral of dt/dr along the direct or indirect photon path; the paper composes this with a geodesically integrated pulsar orbit, samples emission times in proper time, computes the relativistic propagation time at each sample, and inverts the interpolated TOA(τ) function to define τ(TOA). This replaces the 1PN delay sum of Eq. (27), whose Rømer and Shapiro terms assume straight-line propagation and whose geometric delay is only a weak-lensing correction.","core_discovery":"The paper's central claim is that the 1PN timing model—straight-line Rømer delay, Shapiro delay, and weak-lensing geometric delay—cannot reproduce the relativistic photon propagation time for a pulsar orbiting a supermassive black hole. For a circular orbit at 100 gravitational radii, the discrepancy reaches about 2 seconds at superior conjunction, and for an S2-like orbit at the same semimajor axis as Sgr A* the discrepancy is about 0.1 seconds, three orders of magnitude above the nominal 100-microsecond sensitivity. The authors build a timing pipeline in which the pulsar orbit is integrated geodesically in Schwarzschild spacetime in harmonic coordinates, the emitter-observer problem is solved numerically for the primary photon image, and the monotonic TOA(τ) map is inverted by spline interpolation to generate phase-connected residuals. They then use this pipeline to show how misestimating intrinsic and orbital parameters imprints detectable, growing residuals, implying that a single pulsar could constrain the black hole mass far better than current stellar-orbit measurements.","pith_inferences":["Beyond the paper's claims: the same 1PN delay formulas would also corrupt timing of pulsars orbiting intermediate-mass black holes or dense cluster centers, so the geodesic pipeline likely has applications well outside the Galactic Center.","Beyond the paper's claims: the single-parameter sensitivity study hints at strong degeneracies between orbital parameters and the black hole mass; a full Bayesian fit on simulated TOAs would show how tightly these parameters can actually be separated rather than independently bounded.","Beyond the paper's claims: the claim that current codes fail depends on Eq. (27) faithfully representing them; testing the same orbit with a complete post-Keplerian timing model would either confirm the failure or localize it to the simplified straight-line delay model the paper adopts."],"forward_implications":["If the central claim is correct, 1PN timing residuals for a pulsar-SMBH system can exceed the pulsar period, so phase connection is lost and standard TOA fitting breaks down.","For an S2-like orbit around Sgr A*, the 1PN-versus-geodesic discrepancy is about 0.1 seconds, still three orders of magnitude above the 100-microsecond accuracy goal, so even mildly relativistic orbits require geodesic modeling.","Misestimating the black hole mass or the semimajor axis by about one part in 10^10 produces residuals above the detection threshold within months to a few years, depending on orbital compactness, which would sharply improve current Sgr A* mass measurements.","Extending the geodesic timing model to a rotating black hole spacetime is the stated next step, and in spherical symmetry the method already applies to any spacetime, not only Schwarzschild.","Indirect photon paths that dip toward the photon sphere are included in the propagation-time calculation, so the model captures strong lensing effects that 1PN formulas omit entirely."],"supporting_citations":[{"why":"Supplies the numerical emitter-observer methodology for computing fully relativistic photon propagation times in spherically symmetric spacetimes, which the paper extends into a complete timing model.","marker":"[37]"},{"why":"Supplies the exact analytical Schwarzschild propagation-time solution used to validate the numerical approach against a closed-form result for circular orbits.","marker":"[41]"},{"why":"Supplies the post-Newtonian orbital and delay formalism used to represent the weak-field timing model and to define the 1PN delay sum in Eq. (27).","marker":"[40]"},{"why":"Supplies the 1PN timing-model framework and delay parametrization that the paper identifies with current pulsar timing codes.","marker":"[43]"},{"why":"Supplies the geometric time-delay formula from weak gravitational lensing that appears as the third term in the 1PN delay model.","marker":"[44]"},{"why":"Supplies the 100-microsecond TOA accuracy goal for future Galactic Center pulsar observations, which sets the threshold that the 1PN discrepancies are compared against.","marker":"[33]"},{"why":"Supplies the standard treatment of relativistic Shapiro and lensing delays underlying the post-Newtonian comparison model.","marker":"[4]"}],"fun_headline_variants":["Pulsar timing fails near black holes by seconds","Black-hole pulsars demand full relativity for timing","Standard pulsar models wrong by seconds at galactic center","Geodesic timing solves pulsar errors near Sgr A*","Pulsar timing near black holes needs geodesic light paths"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper's conclusion that current timing codes fail at the Galactic Center rests on the assumption that the simplified 1PN delay formula of Eq. (27), with straight-line Rømer and Shapiro delays plus a weak-lensing geometric term, is what existing timing codes actually implement for a test particle around a point mass; if real codes use more complete post-Keplerian models, the claimed failure is not established.","fun_headline_variants_meta":{"raw":{"variants":["Pulsar timing fails near black holes by seconds","Black-hole pulsars demand full relativity for timing","Standard pulsar models wrong by seconds at galactic center","Geodesic timing solves pulsar errors near Sgr A*","Pulsar timing near black holes needs geodesic light paths"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000227,"raw_usage":{"total_tokens":1422,"prompt_tokens":844,"completion_tokens":578,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":460,"completion_tokens_details":{"reasoning_tokens":498}},"tokens_in":460,"tokens_out":578,"duration_ms":6304,"temperature":1.0,"reasoning_tokens":498,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:44:47.436821+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the Toy 0 or S2-like orbit, generate noise-free TOAs from the geodesic pipeline, and fit them with the 1PN timing model of Eq. (27): a residual structure peaking near superior conjunction at the predicted 1-second-to-2-second level would confirm the paper's claim, whereas residuals below 100 microseconds across the full orbit would falsify it.","supporting_citations":[{"cited_title":"Damour and N","cited_arxiv_id":null,"evidence_quote":"Supplies the post-Newtonian orbital and delay formalism used to represent the weak-field timing model and to define the 1PN delay sum in Eq. (27)."},{"cited_title":"The propagation delay in the timing of a pulsar orbiting a supermassive black hole","cited_arxiv_id":"1806.02547","evidence_quote":"Supplies the 1PN timing-model framework and delay parametrization that the paper identifies with current pulsar timing codes."}],"review_version":1}