{"id":"d49f79e0-07cf-44c3-a460-76425e1fb8db","arxiv_id":"2602.19546","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"A 2PN-accurate pulsar timing model for a pulsar orbiting Sgr A* includes proper motion, aberration, and red noise, and forecasts that proper motion will not break the SMBH spin degeneracy for short-period orbits.","lead":"The authors build a detailed numerical model for timing radio pulses from a pulsar orbiting the Milky Way's supermassive black hole, including subtle light-bending and relativistic effects. It forecasts how precisely future telescopes could measure the black hole's mass, spin, and quadrupole moment, and shows proper motion alone won't lift a key degeneracy.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Neglected 3PN secular periastron advance is ~3 ms for fiducial P_b=0.5 yr (Eq. 13), above the 1 ms noise; <1% spin/quadrupole forecast is optimistic.","rationale":"The reader's weakest assumption correctly identified the 2PN truncation and order-of-magnitude justification as a key vulnerability. My concern is more specific: the paper's own Eq. (13) shows the omitted 3PN secular periastron advance produces timing residuals of ~3 ms for the fiducial P_b=0.5 yr system, exceeding the 1 ms precision used in the Fisher forecasts. This is not merely an environmental worry; it is a concrete known effect whose magnitude is above the noise floor. The paper's degeneracy argument does not neutralize the problem, because an unmodeled secular term of this size will bias the very parameters (χ, q) the forecast claims to measure. My recommendation remains CONDITIONAL (hence UNCHANGED): the model and code are serious contributions, but the central <1% claim depends on demonstrating that including 3PN (and similar) effects does not shift the recovered parameters beyond the statistical error. The proposed test provides that demonstration. I do not see a need to change the reader's verdict, only to sharpen the condition under which the forecast can be trusted.","tokens_in":35641,"tokens_out":15006,"duration_ms":125426,"concrete_test":"Generate mock TOAs for the fiducial system (and for P_b=0.3, 0.2 yr) using the same parameters but with a 3PN-accurate equations-of-motion (or exact Kerr geodesic) integrator, then fit with the paper's 2PN timing model assuming 1 ms white noise over 5 yr. If the best-fit χ or q differs from the injected value by more than the 1σ Fisher uncertainty, the 3PN secular effect is not negligible and the <1% forecast fails. Alternatively, add the 3PN periastron-advance term to the timing model and recompute the Fisher forecasts; a shift in predicted precision above 1% for P_b<0.5 yr would confirm the concern.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central forecast—<1% precision on SMBH mass, spin, and quadrupole for P_b<0.5 yr with 1 ms TOAs—rests on the 2PN truncation in Section 2.2. The paper's own estimate of the neglected 3PN secular periastron advance, Eq. (13), reads δt ~ 1 ms (T_obs/5 yr)(P_b/1 yr)^{-7/3} sin i. For the fiducial system (P_b=0.5 yr, i=π/5, T_obs=5 yr) this is ~3 ms, three times the assumed 1 ms timing precision. For P_b=0.3 and 0.2 yr it grows to ~10-100 ms. Table 1 lists the 'Higher PN' delay as 3e-3 s for the fiducial case, confirming the effect is not negligible. The paper argues this secular effect is 'largely degenerate with the secular effects caused by the environmental perturbations' and so omits it. That is not an adequate justification for the forecast: any unmodeled secular signal at the millisecond level will be partially absorbed by the fitted parameters (especially χ and q) and can bias the recovered values by more than the quoted statistical errors. The paper's own estimates thus undermine the central claim precisely in the orbital-period regime where <1% precision is asserted. A quantitative bias analysis is required.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":36054,"tokens_out":8059,"duration_ms":82451,"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":[{"comment":"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":"Section 2.2, Eq. (13); Table 1"},{"comment":"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.","section":"Section 2.2; Table 1; Section 3.2"}],"minor_comments":[{"comment":"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":"Section 2.3, Eq. (38)"},{"comment":"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":"Section 2.2, after Eq. (13)"},{"comment":"'Geodesics precession' should be 'Geodetic precession'.","section":"Section 2.2, Eq. (18)"},{"comment":"Minor typo: 'inverse the numerical relation' should be 'invert the numerical relation'.","section":"Appendix C or Section 2.4"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is technically solid and the convergence tests are a genuine strength. The recommendation is major revision rather than rejection because the 3PN bias issue is fixable within the paper's scope, either by including the term or by adding a mismodeling bias analysis and appropriately limiting the forecast. I did not find evidence of circularity or hidden free parameters; the main risk is that the headline precision claim is presented as unconditional when the model is only demonstrated for an isolated SMBH at 2PN order."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","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.\n\nThe 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.\n\nDoes 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.\n\nThe 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.\n\nWho 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.","headline":"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.","tokens_in":36547,"tokens_out":4104,"would_cite":true,"duration_ms":39446,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["pulsar timing","supermassive black holes","Sagittarius A*","post-Newtonian approximation","no-hair theorem","red noise","proper motion","parameter estimation"],"falsifier":"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.","tokens_in":35532,"feed_emoji":"🕳️","tokens_out":6399,"duration_ms":59197,"temperature":0.7,"pith_summary":"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.","feed_headline":"One pulsar around Sgr A* could measure mass, spin, quadrupole to <1%","feed_subtitle":"Tight-orbital timing at the millisecond level would give a no-hair test of the Galaxy's black hole.","key_machinery":"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.","core_discovery":"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","pith_inferences":["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."],"forward_implications":["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."],"fun_headline_variants":["Pulsar timing model predicts <1% measures of Sgr A* mass, spin, quadrupole","Tight pulsar orbits around Sgr A* could reveal black hole hair to <1%","Sub-percent Sgr A* parameters within reach via new pulsar timing model","SKA-era pulsar timing could test gravity with <1% on Sgr A* parameters","Model for pulsar-SMBH systems enables precise measurements of Sgr A*"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Pulsar timing model predicts <1% measures of Sgr A* mass, spin, quadrupole","Tight pulsar orbits around Sgr A* could reveal black hole hair to <1%","Sub-percent Sgr A* parameters within reach via new pulsar timing model","SKA-era pulsar timing could test gravity with <1% on Sgr A* parameters","Model for pulsar-SMBH systems enables precise measurements of Sgr A*"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000389,"raw_usage":{"total_tokens":1904,"prompt_tokens":776,"completion_tokens":1128,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":520,"completion_tokens_details":{"reasoning_tokens":1010}},"tokens_in":520,"tokens_out":1128,"duration_ms":10086,"temperature":1.0,"reasoning_tokens":1010,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T21:35:04.757674+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[],"review_version":1}