{"id":"36e277bb-52c7-49b8-b560-428dbcc1802e","arxiv_id":"2411.14895","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"low","formal_verification":"none","parameter_count":6,"one_line_summary":"Joint radio and gamma-ray timing of nine millisecond pulsars yields precise astrometric and spin parameters, including proper motions for four pulsars.","lead":"Astronomers combined radio and gamma-ray observations of nine newly discovered millisecond pulsars to measure their positions, spins, and orbits more precisely than either band alone. The joint analysis yielded proper motion measurements for four pulsars and a possible first detection of a relativistic signal in one binary system.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The reported proper-motion detections for two of the four claimed MSPs are not significant from the quoted 1σ errors; the headline overstates the joint-timing gain.","rationale":"I read the paper as claiming a methodological advance whose showcase result is precise proper motions for four MSPs. The reader's weakest assumption, concerning gamma-ray photon-weight mis-modelling, is a legitimate secondary concern, but it is anchored to J1709−0333, which is not one of the four proper-motion pulsars, and no evidence connects the weight problem to J1526−2744, J1623−6936, J1757−6032, or J1858−5422. By contrast, the significance mismatch is visible directly in Table 4 and affects the paper's own headline. I therefore disagree with the reader's identification of the most load-bearing weakness. This does not overturn the paper: the joint-timing method, and at least the J1623−6936 and J1858−5422 proper-motion detections, stand. However, the claim of four measured proper motions needs a formal detection-statistic check, so the conditional recommendation should remain, with that check made an explicit condition.","tokens_in":26746,"tokens_out":10402,"duration_ms":103757,"concrete_test":"Re-run the joint MCMC for each of the four pulsars with proper-motion components fixed to zero and compute the delta-log-likelihood / delta-chi² (or Bayes factor) between the null and full models using identical priors. Report the number of pulsars exceeding a pre-defined detection threshold (e.g., Δχ² ≥ 9 for two extra parameters, corresponding to ~3σ). Independently, split the gamma-ray data into two halves and fit proper motion separately; the two estimates should agree within their joint uncertainties. If fewer than four pulsars pass, revise the abstract and Sect. 7 claim to the actual number of secure detections and present the remaining proper motions as upper limits.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline claim that joint timing 'significantly measure[s] the proper motion of four MSPs' (Sect. 7) is not supported by the quoted 1σ errors in Table 4. For the four jointly timed binaries, the total proper motions are: J1526−2744: 9±7 mas/yr (~1.3σ), J1623−6936: 8.6±1.8 mas/yr (~4.8σ), J1757−6032: 8±4 mas/yr (~2σ), J1858−5422: 9±3 mas/yr (~3σ). Only two are secure at the conventional 3σ level; J1757 is consistent with zero, and the reported total of 8±4 is difficult to reconcile with its component values (−2.3±1.4, −3.0±2.3 mas/yr). J1526 has a 3σ PMRA component, but the 2D proper-motion vector is only ~2.5σ from zero because PMDEC carries a 14 mas/yr error, so it is not a precise total-proper-motion measurement. Since the paper's distinctive scientific payoff is precisely these four proper motions, the statistical evidence should be presented through a formal null-model comparison, with the number of claimed detections reduced if they do not pass a defined threshold.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper reports a multi-telescope radio and Fermi-LAT gamma-ray timing campaign on nine millisecond pulsars discovered in the TRAPUM L-band survey of unassociated Fermi sources. Gamma-ray pulsations are found for six targets; for five of them a new joint radio/gamma-ray timing analysis using MCMC is applied, and for the redback J1803−6707 a Gaussian-process model of orbital period variations is used to extend the timing solution across the Fermi mission. The paper presents spin, astrometric, and orbital parameters, proper-motion measurements for four pulsars, eclipse light curves for the two redbacks, phase-aligned radio/gamma-ray profiles, and polarimetry. The central claim is that joint timing yields precise short- and long-term parameters, in particular proper motions of four MSPs that a gamma-ray-only analysis would not have provided.","tokens_in":27073,"tokens_out":6364,"duration_ms":65910,"significance":"The joint timing code and its application to five pulsars is a useful methodological advance, and the treatment of J1803−6707's orbital phase variations with a Gaussian process is more sophisticated than the usual redback timing analysis. The paper is transparent about the marginal Shapiro delay and the photon re-weighting for J1709−0333, and it makes ephemerides publicly available. If the proper-motion detections survive a formal null-model test, the demonstration that 15 years of gamma-ray photons can be combined with radio ToAs to measure astrometric parameters would be valuable for the Fermi MSP population. The polarimetry and eclipse characterisation are useful additions. At present, however, the headline claim of four significant proper motions is not supported by the quoted uncertainties.","major_comments":[{"comment":"The summary states that joint timing 'significantly measure[s] the proper motion of four MSPs,' but the quoted 1σ errors in Table 4 do not support this. J1526−2744 has a total proper motion of 9±7 mas/yr (≈1.3σ), and J1757−6032 has 8±4 mas/yr (≈2σ); only J1623−6936 (8.6±1.8) and J1858−5422 (9±3) reach the 3σ level. Moreover, the tabulated total proper motions are internally inconsistent with the components: for J1757−6032, μ_α = −2.3±1.4 and μ_δ = −3.0±2.3 imply a magnitude ≈3.8 mas/yr, not 8; for J1526−2744 the components imply ≈12.6 mas/yr, not 9. The Shklovskii corrections in Table 4 are applied to J1526 and J1757 although their total proper motions are not significant. I request a formal null-model comparison (proper motion fixed to zero versus free) and a revision of the number of claimed detections and of the derived Pdot_int corrections.","section":"Section 7 / Table 4"},{"comment":"For PSR J1709−0333, the reported H = 76 is obtained after maximizing an energy-dependent re-weighting over the scale factor A and spectral index γ, and after an earlier tweak of fdot. No trials factor for this optimization is discussed, although the paper itself notes that 'much larger values are required to overcome the trial factors' arising from parameter searches. Please provide a trials-corrected significance or otherwise state how many independent re-weighting trials were performed; without this, the gamma-ray detection claim for this pulsar is not fully quantified.","section":"Section 3"},{"comment":"The text states that fitting both Shapiro-delay parameters yields no significant detection, yet immediately derives a strong inclination constraint (median 87.5° with a 68% interval of +1.4/−2.8 deg) from a χ² map in the M2–cos i plane. The dependence of this constraint on the assumed priors (pulsar mass range 1.17–3.2 M⊙, uniform cos i) and on the earlier assumption of i = 87.5° for the h3 measurement needs to be made explicit. As written, the paper risks overinterpreting a marginal signal as a measurement of the orbital inclination.","section":"Section 5.1"}],"minor_comments":[{"comment":"Section 4.1 says the full joint timing analysis was performed on four binary pulsars and one isolated pulsar, while Section 7 says the new joint technique was applied to the six sources with pulsations in both bands; please clarify how PSR J1803−6707 fits this description.","section":"Section 4.1 / Section 7"},{"comment":"The pulsar name is given as J1823−3543 in Table 4 but as J1823−3544 in the abstract, key words, and Section 5.3; this should be unified.","section":"Table 4 / Table 3"},{"comment":"The text refers to 'J1623−6939' in the first paragraph of Section 5.3; this should be J1623−6936.","section":"Section 5.3"},{"comment":"The pulsar is called J1803−6708 in the data availability section and in Section 7, whereas the rest of the paper uses J1803−6707.","section":"Section 7 / Data availability"},{"comment":"For PSR J1036−4353, the statement that W² = 134 means pulsations would not be expected on timescales shorter than about 3 years is unclear without a definition of the relationship between W² and the expected H-test power; please spell out the scaling or add a reference.","section":"Section 3"}],"recommendation":"major_revision","confidential_remarks":"The timing infrastructure and data products are sound, and the main result can be repaired with a statistical reanalysis. The proper-motion claim should be re-derived with a null-model comparison, and the number of claimed detections should be reduced if the comparison does not support them. The manuscript is within the scope of A&A, but I would not recommend acceptance before the proper-motion significance and the internal inconsistencies in Table 4 are addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a serious read. This is the first full joint radio/gamma-ray timing treatment of the TRAPUM Fermi-LAT discoveries, and it delivers coherent timing solutions for nine MSPs, proper motions for four, and a Gaussian-process model for the redback's orbital wander. The joint likelihood (PINT + emcee, combining radio chi-square with the unbinned gamma-ray likelihood) extends earlier work by Nieder et al. and Thongmeearkom et al., and the corner plots show well-behaved posteriors. The authors are transparent about the marginal Shapiro delay and about the photon re-weighting for J1709-0333.\n\nThe main soft spot is the proper-motion headline. Section 7 says four MSPs have significantly measured proper motions, but Table 4 does not support that for all four. J1526-2744's total is 9±7 mas/yr (about 1.3σ). J1757-6032's total is quoted as 8±4, while the components are -2.3±1.4 and -3.0±2.3—magnitude roughly 3.8, which is hard to reconcile with 8. Only J1623-6936 (8.6±1.8, ≈4.8σ) and J1858-5422 (9±3, ≈3σ) are secure. I'd ask for a formal null-model comparison or, at minimum, a correction of the J1757 numbers and a trimmed claim.\n\nOther issues are minor. The h3 Shapiro delay in J1757 is computed with inclination fixed at 87.5°, and the authors say so; a casual reader could still over-read the mass prospects. The J1709 re-weighting maximises the H-test over two parameters without a trials penalty; H going from 44 to 76 is probably fine, but it should be quantified. The Gaussian-process code for J1803-6707 is still 'in preparation' rather than released, which weakens reproducibility. And the broad worry that 4FGL weights are wrong for the other pulsars—there is no actual evidence of that in the corner plots or residuals, so I wouldn't hold it against the paper beyond the J1709 case.\n\nThis is a paper for pulsar timing specialists, Fermi MSP searchers, and anyone using MSP proper motions for population or gravitational-wave work. The core timing solutions look robust and the joint-timing implementation is a useful addition to the toolkit. I'd send it to a serious referee, with the J1757 proper-motion inconsistency flagged as a must-fix before acceptance.","headline":"Solid joint radio/gamma-ray timing paper whose headline proper-motion claims outrun the quoted errors; two of four are secure.","tokens_in":27689,"tokens_out":4954,"would_cite":true,"duration_ms":44131,"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":"Jointly fitting radio pulses and gamma-ray photons lets the authors measure proper motions of four millisecond pulsars that a gamma-ray-only analysis of 15 years of Fermi data could not detect.","keywords":["millisecond pulsars","gamma-ray timing","radio timing","proper motion","joint timing analysis","redback pulsars","Fermi-LAT","pulsar timing"],"falsifier":"Measure the proper motion of one of the four jointly timed pulsars (e.g., J1858−5422 or J1623−6936) with very long baseline interferometry. If the VLBI proper motion disagrees with the joint-timing value by more than the combined uncertainty, the gamma-ray weighting model is biased for that pulsar.","tokens_in":26576,"feed_emoji":"🛰️","tokens_out":8374,"duration_ms":72278,"temperature":0.7,"pith_summary":"This paper claims that a joint analysis of radio pulse arrival times and gamma-ray photons can extract timing parameters from millisecond pulsars that neither band alone provides. Using a newly developed Markov-chain Monte Carlo code that fits both data sets simultaneously, the authors measure the proper motion of four pulsars discovered in a radio survey of gamma-ray sources, even though a gamma-ray-only analysis of the full 15-year Fermi data span would not have been sensitive to it. The result matters because proper motions are needed to correct observed spin-down rates for kinematic effects and to obtain accurate pulsar energetics, and because precise long-baseline timing solutions enable multi-messenger follow-up such as continuous gravitational-wave searches.","feed_headline":"Joint radio and gamma-ray timing yields four pulsar proper motions","feed_subtitle":"Fifteen years of gamma-ray data plus radio pulses pin down motions that either band alone cannot measure.","key_machinery":"The central machinery is a joint log-likelihood function that adds the radio timing chi-squared (from pulse arrival-time residuals) to the unbinned gamma-ray template-based likelihood (from photon phases), plus priors, and optimizes all timing parameters together with a Markov-chain Monte Carlo sampler. A second piece of machinery is the Gaussian-process description of redback orbital-phase variations, with a Matérn covariance function, which lets the gamma-ray timing extend beyond the valid range of a deterministic radio ephemeris. The capacity that carries the argument is that radio data pin down short-term parameters and sky position with high signal-to-noise, while the long 15-year gamma-ray baseline drives the measurement of parameters that improve with time span, such as proper motion and spin-frequency derivative.","core_discovery":"On the paper's own terms, the key result is that joint radio and gamma-ray timing works: for six of nine millisecond pulsars discovered in the survey, gamma-ray pulsations were found using radio ephemerides, and a joint timing fit that combines the radio times of arrival with the sparse gamma-ray photons breaks the degeneracies that prevented either data set alone from converging. This yields significant proper-motion measurements for four pulsars, which the 15-year gamma-ray data alone could not provide, and allows the measured spin-down rates to be corrected for the Shklovskii effect. For the redback J1803−6707, a Gaussian-process treatment of its erratic orbital-period variations extends the gamma-ray timing solution across the full Fermi mission, and for J1757−6032 the timing reveals a nearly edge-on orbit with a hint of Shapiro delay, identifying it as a promising target for neutron-star mass measurement.","pith_inferences":["The joint timing technique should transfer to any Fermi-selected millisecond pulsar with detectable gamma-ray pulsations, potentially reducing the radio follow-up investment needed for proper-motion measurements across the whole population.","If the photon-weight bias seen for J1709−0333 is common in complex gamma-ray regions, some previously published gamma-ray-only timing positions and proper motions may need revision; re-processing with joint radio data could serve as a check.","The Gaussian-process treatment of orbital variations could be applied to other redbacks and spider pulsars where deterministic timing models fail, potentially recovering gamma-ray pulsations from systems currently thought to be gamma-ray quiet.","The shorter spin-period distribution of this sample, if it persists in larger targeted surveys, suggests that sensitive radio surveys of Fermi unassociated sources preferentially uncover fast millisecond pulsars, which are the best clocks for pulsar timing arrays."],"forward_implications":["The joint timing approach yields proper motions for millisecond pulsars that would otherwise require many more years of radio timing, substantially shortening the follow-up time needed to correct spin-down rates and derive pulsar energetics.","For redbacks with erratic orbital periods, the Gaussian-process timing method makes gamma-ray pulsations usable across the full Fermi data span, enabling long-term timing solutions that deterministic ephemerides cannot provide.","The high-inclination, high-companion-mass system J1757−6032 is now a strong candidate for a neutron-star mass measurement via Shapiro delay, which would test whether Case A Roche-lobe overflow produces massive neutron stars.","The joint timing solutions provide ephemerides that are valid over the Fermi mission, allowing multi-wavelength and multi-messenger searches, including continuous gravitational waves, to use these pulsars as targets."],"supporting_citations":[{"why":"Discovery paper for the nine MSPs and their initial radio ephemerides; this work extends those with long-term joint timing.","marker":"Clark et al. 2023"},{"why":"Supplies the 4FGL-DR3 gamma-ray source catalog and spectral models from which photon weights are computed.","marker":"Abdollahi et al. 2022"},{"why":"Provides the H-test statistic and weighted-photon formalism used to detect and evaluate gamma-ray pulsations.","marker":"Kerr 2011"},{"why":"Provides the pulsar timing infrastructure that the joint timing code uses to fold radio and gamma-ray data.","marker":"Luo et al. 2021"},{"why":"Supplies the affine-invariant MCMC sampler used to explore the joint parameter space.","marker":"Foreman-Mackey et al. 2013"},{"why":"Gives the kinematic proper-motion correction to the spin-down rate, which motivated the proper-motion measurements.","marker":"Shklovskii 1970"},{"why":"Introduces the Gaussian-process treatment of orbital phase variations used to time the redback J1803−6707 across the Fermi data.","marker":"Clark et al. 2021"},{"why":"Refines the Gaussian-process timing method with Gibbs sampling; this paper builds on that procedure.","marker":"Thongmeearkom et al. 2024"}],"fun_headline_variants":["Four pulsar proper motions from joint radio and gamma-ray timing","Joint timing breaks single-band deadlock for six millisecond pulsars","Radio+gamma-ray fit yields proper motions gamma-only data missed","New MCMC timing code unlocks gamma pulsations in six MSPs"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The gamma-ray photon weights, which say how likely each photon is to come from the pulsar rather than background, are assumed to be correct; if they are systematically wrong for any target, the joint timing fit returns biased positions and proper motions.","fun_headline_variants_meta":{"raw":{"variants":["Four pulsar proper motions from joint radio and gamma-ray timing","Joint timing breaks single-band deadlock for six millisecond pulsars","Radio+gamma-ray fit yields proper motions gamma-only data missed","New MCMC timing code unlocks gamma pulsations in six MSPs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000281,"raw_usage":{"total_tokens":1644,"prompt_tokens":906,"completion_tokens":738,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":522,"completion_tokens_details":{"reasoning_tokens":665}},"tokens_in":522,"tokens_out":738,"duration_ms":7926,"temperature":1.0,"reasoning_tokens":665,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T14:45:00.676063+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the proper motion of one of the four jointly timed pulsars (e.g., J1858−5422 or J1623−6936) with very long baseline interferometry. If the VLBI proper motion disagrees with the joint-timing value by more than the combined uncertainty, the gamma-ray weighting model is biased for that pulsar.","supporting_citations":[{"cited_title":"J., Breton, R","cited_arxiv_id":null,"evidence_quote":"Discovery paper for the nine MSPs and their initial radio ephemerides; this work extends those with long-term joint timing."},{"cited_title":"J., Nieder, L., Voisin, G., et al","cited_arxiv_id":null,"evidence_quote":"Introduces the Gaussian-process treatment of orbital phase variations used to time the redback J1803−6707 across the Fermi data."}],"review_version":1}