{"id":"cf9f8ce9-37b8-4c4a-ae9c-f83bc5362b9c","arxiv_id":"2505.05021","paper_version":2,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"A new binary millisecond pulsar, M13I, was discovered in globular cluster M13 via a Fast Folding Algorithm search of FAST data, with a 6.37 ms spin period and an 18.23-day eccentric orbit.","lead":"Astronomers used a folding algorithm on FAST telescope data to search for pulsars in 30 globular clusters and found a new millisecond pulsar in M13. The new object, M13I, has an unusually wide and eccentric orbit, suggesting it formed differently from its cluster-mates.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"M13I's cluster membership rests on position and DM alone; a proper-motion fit to the existing 6-year timing data would settle whether the pulsar is truly in M13 and whether the evolutionary claims in Section 4.3 apply.","rationale":"Good-faith reading: the paper's main deliverable is a well-timed new 6.37 ms binary pulsar; the 68-ToA, 6-year coherent solution with 53 us RMS is strong evidence that the pulsar is real and that the orbital parameters are correctly measured. Independent support includes blind FFA re-detection of 38/39 isolated and 49/54 binary known pulsars, which validates the pipeline. The weakest link is association with M13: only angular proximity and DM consistency, not a dynamical measurement. This matches the reader's weakest_assumption, and I agree with that characterization. A proper-motion fit is the natural, low-cost check because the timing data already span 6 years; M13's PM is well measured by Gaia and should differ from the typical proper motion of field pulsars along this line of sight. If the test failed, the discovery of the pulsar itself would survive, but the claim of a new M13 cluster pulsar and the evolutionary discussion would need to be downgraded. Because the prior for chance alignment is small and such positional-plus-DM evidence is routinely used in globular cluster pulsar work, this concern does not change the ACCEPT verdict; it argues for reporting the PM fit or explicitly labeling membership as probable rather than certain. The FFA-versus-FFT sensitivity comparison is also not fully controlled, but that is secondary and does not affect the existence or timing of M13I.","tokens_in":12364,"tokens_out":11278,"duration_ms":130108,"concrete_test":"Re-fit the 68 ToAs used for Table 3 with TEMPO or PINT, adding proper motion (and optionally parallax) to the timing model, and compare the fitted proper motion with M13's Gaia DR3 proper motion (roughly -2.9, -4.9 mas/yr). If a proper-motion fit converges with values within 2-3 sigma of the cluster PM and residuals remain at the reported 53 us level, cluster membership is strongly supported; if the fit requires a PM differing by more than ~10 mas/yr, or fails to converge even with additional spin-frequency derivatives, the association should be treated as unconfirmed and the Section 4.3 conclusions qualified accordingly.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that M13I is a new pulsar in M13 (title, abstract, Section 3.2) depends on the assumption, made in Section 4.3, that M13I is physically associated with the cluster. The evidence is an angular offset of 11.7 arcsec from the cluster center and DM = 29.58 pc cm^-3 lying within the 27-34 pc cm^-3 range used for M13; there is no independent distance, parallax, or proper-motion measurement. If M13I were a foreground or background field binary pulsar, the object would still be a real pulsar, but the headline discovery 'in M13', the 'nearer hemisphere' inference from DM, and the comparison with the other M13 binaries (the eccentricity trend and the distinct evolutionary path) would lose their basis. The a priori chance of a field pulsar appearing within 11.7 arcsec of the cluster center is small, so this does not seriously threaten the reality of the pulsar, but it is the least secure link in the chain and the one on which the astrophysical interpretation rests.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports an FFA-based pulsar search of 16 globular clusters observed with FAST, using the RIPTIDE package. The pipeline re-detects 38/39 isolated and 49/54 binary pulsars among previously known sources, with non-detections attributed to orbital acceleration or low SNR. The main result is the discovery of M13I (PSR J1641+3627I), a 6.37 ms binary millisecond pulsar in M13, with a six-year phase-coherent timing solution yielding an 18.23-day orbit, eccentricity 0.064, and companion mass 0.45-1.37 Msun. The paper also compares FFA and FFT detection rates for M13I, reports a computational speed comparison, and discusses the pulsar's orbital properties as evidence for a distinct evolutionary path.","tokens_in":12568,"tokens_out":10072,"duration_ms":86179,"significance":"The discovery of M13I is a valuable addition to the known pulsar population in globular clusters, and the timing solution appears robust, with 68 ToAs and 53 us RMS residuals. The re-detection statistics provide a useful benchmark for FFA pipelines on FAST data. If the cluster membership is confirmed, the system's long orbital period and moderate eccentricity make it an interesting test case for binary evolution in dense cluster environments. The work also demonstrates the practical feasibility of FFA searches on FAST's high-rate data, which is relevant for future surveys. However, the significance of the evolutionary interpretation is conditional on the assumed cluster membership and on the validity of the FFA-vs-FFT comparison.","major_comments":[{"comment":"The physical association of M13I with the globular cluster M13 is asserted from the 11.7 arcsec angular offset and DM = 29.58 pc cm^-3 lying within the 27-34 pc cm^-3 range in Table 1, but no proper-motion measurement, parallax, or statistical chance-coincidence estimate is presented. This assumption is load-bearing: the title and abstract call M13I a pulsar 'in M13', and Section 4.3 uses the DM to place it in the nearer cluster hemisphere and uses the eccentricity to claim a contradiction with the Wang et al. (2020) trend. I request either a quantitative estimate of the probability that a field pulsar with this DM falls within 11.7 arcsec of the cluster center, or a proper-motion constraint from the 6-year timing data. At minimum, the wording in Section 4.3 should distinguish between the measured timing parameters and the assumed cluster membership.","section":"Section 3.2 and Section 4.3"},{"comment":"The paper claims that FFA detected M13I in 14 of 78 observations (17.9%) versus 8 of 78 (10.3%) for FFT, and Section 5 concludes that 'FFA demonstrates superior sensitivity to weakly accelerated signals.' However, the FFT search parameters are not specified: no harmonic-summing setting, acceleration-search zmax, period range, DM step, or SNR threshold is given for the comparison. Without a controlled setup, the 1.75x factor may reflect differences in search configuration rather than algorithmic sensitivity. Please state the exact FFT pipeline used (or clarify that the comparison is against archival PRESTO processing with a specific configuration) and temper the conclusion to a statement about this dataset. If a controlled comparison is not feasible, the sentence in the conclusion should be revised.","section":"Section 4.1 and Figure 4"},{"comment":"The explanation for the non-detection of M12B and M13H states that their period derivatives 'exceed 10^-11 s/s, far above the detectability threshold' derived from |\\dot{P}| ≤ 2(P/T)^2. For M13H, P = 11.21 ms and T = 1 hour give a threshold of approximately 1.9×10^-11 s/s, so a value just above 10^-11 is not 'far above.' Please report the actual measured \\dot{P} values (or at least their order of magnitude) for both pulsars, or rephrase the statement to be quantitatively accurate for each object.","section":"Section 4.1, detectability condition"}],"minor_comments":[{"comment":"The pulsar designation is inconsistent: the abstract gives PSR J1641+3627I, Section 3.2 gives J1641+3624I, and Table 3 coordinates (Dec +36:27:36.5) support the former. Please correct.","section":"Section 3.2"},{"comment":"Table 1 lists 9 known pulsars in M13 and 8 FFA detections, but the text in Section 3.1 says 93 known pulsars were analyzed, before the discovery of M13I. Please clarify whether the table includes the new pulsar and reconcile the total count.","section":"Section 3.1 and Table 1"},{"comment":"The statement that M13I was 'initially overlooked in earlier searches' is followed by the report of a single-day archival PRESTO detection; please clarify whether that detection was previously recognized as a pulsar or only in hindsight.","section":"Section 3.2, first paragraph"},{"comment":"The text first attributes their non-detection to weak signals (SNR<7 after coherent folding) and then to 'significant orbital acceleration effects'; please reconcile these explanations or state which is dominant.","section":"Section 4.1, M14D/M14E"},{"comment":"The speed comparison uses accelsearch with zmax=0, i.e., no acceleration search, while rffa is a full folding search; this is not a like-for-like benchmark and should be stated as a limitation.","section":"Section 4.2"},{"comment":"There are numerous typographical and grammatical errors, including 'puslars', 'de-dispered times series', 'redetectining', 'summerized', 'highghting', and 'privious'. A careful proofread is needed.","section":"General"},{"comment":"The two Losovsky & Dumsky 2014 entries appear to refer to the same work with different author-name formatting; merge or disambiguate them.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The timing solution for M13I is convincing and the re-detection statistics are useful, but the cluster-membership claim and the uncontrolled FFA-vs-FFT comparison are the main risks. Both can be addressed in revision, so I recommend major revision rather than rejection. Please also ensure the pulsar-name inconsistency is resolved before acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know this paper delivers a genuinely new binary millisecond pulsar, M13I, with a solid six-year phase-coherent timing solution (68 TOAs, 53 us RMS). It is a real addition to the globular cluster census. The FFA pipeline on FAST data re-detects most known pulsars in the searched clusters, which gives me confidence the search itself is sound. The discovery claim does not depend on the sensitivity comparison, and the timing parameters are well constrained. What is genuinely new: M13I is a 6.37 ms pulsar with an 18.23 day orbit and eccentricity 0.064, clearly different from the other M13 binaries. Applying RIPTIDE/FFA to FAST cluster data is a new combination, and the re-detection statistics are a useful benchmark for the method even if the algorithm is not new. The soft spots are minor but real. Cluster membership for M13I rests on an angular offset of 11.7 arcsec and a DM of 29.58 pc/cc within the cluster's range; there is no proper motion or parallax. If it is a field pulsar, the evolutionary discussion in Section 4.3 loses its basis. A proper-motion fit to the existing timing data would settle this cleanly. The FFA-versus-FFT sensitivity comparison is not a controlled experiment: the FFT search parameters and acceleration search settings are not specified, so the 17.9% versus 10.3% detection rate should be read as illustrative, not proof. There is also a minor inconsistency: the abstract gives J1641+3627I while Section 3.2 gives J1641+3624I, and the non-detections in Table 2 include two pulsars with SNR below threshold that are not high-acceleration cases despite the abstract's wording. None of this undermines the central result. The timing solution is reproducible, the re-detection statistics are strong, and the authors are appropriately honest about the RM uncertainty and the sensitivity limits. This is a solid, incremental contribution that deserves a serious referee. I would recommend acceptance after minor revisions that tighten the membership argument and present the sensitivity comparison with appropriate caveats.","headline":"A solid search paper that finds a genuinely new binary millisecond pulsar in M13 with a credible timing solution; the main caveat is that cluster membership rests on position and DM alone, not proper motion.","tokens_in":682,"tokens_out":1667,"would_cite":true,"duration_ms":32438,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A fast-folding search of 16 globular clusters has found a new 6.37-ms binary pulsar, M13I, and recovered 87 of 93 known pulsars.","keywords":["pulsars","globular clusters","Fast Folding Algorithm","millisecond pulsars","binary pulsars","FAST telescope","M13","timing solution"],"falsifier":"A proper-motion or parallax measurement of M13I that disagrees with M13's cluster motion would falsify the physical association; alternatively, an independent blind re-reduction of the same FAST observations with a public FFA implementation that fails to reproduce M13I's signal or the reported re-detection counts would falsify the survey result.","tokens_in":12178,"feed_emoji":"📡","tokens_out":6928,"duration_ms":64422,"temperature":0.7,"pith_summary":"The paper argues that applying the Fast Folding Algorithm (FFA) to L-band observations of globular clusters taken with the Five-hundred-meter Aperture Spherical radio Telescope recovers nearly every known pulsar and can uncover faint signals that Fourier-transform searches miss. On that basis it reports the discovery of M13I (PSR J1641+3624I), a 6.37-ms pulsar in an 18.23-day eccentric binary (e = 0.064) inside the cluster M13, with a phase-coherent timing solution spanning six years. The detection counts (38 of 39 isolated, 49 of 54 binary pulsars) and the new pulsar's orbital properties are the paper's central results. A sympathetic reader would care because the new system occupies a part of orbital-parameter space not previously seen in M13, suggesting an evolutionary path distinct from the cluster's other binary pulsars.","feed_headline":"New 6.37-ms binary pulsar found in cluster M13","feed_subtitle":"Fast-folding search of FAST data re-detects 87 of 93 known cluster pulsars and catches a wide-orbit pulsar FFT missed","key_machinery":"The central object is the Fast Folding Algorithm (FFA), a time-domain search that folds a dedispersed time series at each trial period and sums the folded profile over the full observation, implemented here in a software package whose search command operates on whitened time series. FFA's sensitivity advantage over FFT is largest for weak signals with long spin periods, and its computational cost scales as the inverse square of the minimum search period, which is why the search adopted cluster-specific minimum periods of 3-50 ms and a period ceiling of 100 s. The paper also uses the detectability condition $|\\dot P| \\le 2(P/T)^2$ to show that binaries with period derivatives above about $10^{-12}$ s/s accumulate more than one cycle of phase drift in a one-hour fold and therefore evade a fixed-period FFA search, quantitatively explaining the non-detected pulsars.","core_discovery":"The paper's central discovery is M13I, a millisecond pulsar with spin period 6.37 ms and dispersion measure 29.58 pc $cm^{-3}$, located 11.7 arcseconds from the center of globular cluster M13. Using 78 FAST observations over six years the authors obtained a phase-coherent timing solution showing an 18.233778-day orbit with eccentricity 0.064356, a projected semi-major axis of 20.71515 light-seconds, and a companion mass between 0.45 and 1.37 solar masses (median 0.54), consistent with a helium white dwarf. The same FFA pipeline re-detected 38 of 39 known isolated pulsars and 49 of 54 known binary pulsars across 16 clusters, with the six non-detections attributed to orbital acceleration, weak signal-to-noise, or use of a known ephemeris. The paper claims that FFA detected M13I in 14 of 78 observations (17.9%) versus 8 of 78 (10.3%) for FFT searches, and that M13I's wide, eccentric orbit breaks the eccentricity-versus-distance trend previously seen among M13 binaries.","pith_inferences":["If M13I is truly a cluster member, its low dispersion measure places it in the near hemisphere of M13; a future proper-motion or timing-acceleration measurement could test whether its wide eccentric orbit results from a recent dynamical encounter rather than from isolated binary evolution.","The inverse-square cost scaling implies that a rerun of the same pipeline with the minimum search period lowered to 1-2 ms, when computing budgets allow, could reveal sub-3-ms pulsars that the current search's harmonic-only sensitivity missed.","Extending the same FFA approach to FAST archival data for other globular clusters beyond the 30 studied here is a natural next step and could be done without new telescope time.","Adding an acceleration-search stage to the FFA pipeline would likely recover the five high-acceleration binaries missed here, closing most of the completeness gap."],"forward_implications":["For isolated cluster pulsars, FFA recovery is essentially complete: 38 of 39 were re-detected, so the one failure reflects signal weakness or a stacking-only detection, not a general method gap.","The six non-detections are quantitatively explained by the $|\\dot P| \\le 2(P/T)^2$ criterion, so future FFA searches can predict which binaries will be missed and can schedule shorter integrations or acceleration fits for them.","M13I's orbital parameters double the spread of M13 binary pulsars along both period and eccentricity, implying that M13's pulsar population is not a single coeval, tidally circularized family.","FFA's higher detection rate on M13I (14 of 78 versus 8 of 78 observations) suggests that reprocessing existing FAST globular-cluster data with FFA, rather than taking new observations, is a viable route to finding more faint pulsars.","No new pulsars emerged from 138,448 candidates in 14 clusters with no known pulsars, indicating that any pulsars there are fainter or shorter-period than this search could reach."],"supporting_citations":[{"why":"Introduces the Fast Folding Algorithm, the method the paper applies to cluster pulsar searches.","marker":"Staelin 1969"},{"why":"Provides the FFA implementation and the sensitivity and computational-scaling analysis that justify the pipeline's parameter choices.","marker":"Morello et al. 2020"},{"why":"Supplies the processing toolkit used for RFI masking, de-dispersion, folding, and initial orbital fitting.","marker":"Ransom et al. 2002"},{"why":"Reports the earlier M13 pulsar discoveries and the eccentricity-versus-distance trend that M13I contradicts.","marker":"Wang et al. 2020"},{"why":"Documents the FAST globular-cluster survey and the acceleration effects used to explain non-detections such as M14D and M14E.","marker":"Pan et al. 2021a"},{"why":"Provides the YMW16 Galactic electron-density model used to estimate dispersion measures for clusters without known pulsars.","marker":"Yao et al. 2017"},{"why":"Characterizes FAST's gain and system temperature, which enter the flux-density estimate for M13I.","marker":"Jiang et al. 2019"},{"why":"Supplies the dynamical-encounter framework for how cluster-core interactions can raise binary eccentricity, used to interpret M13I.","marker":"Heggie & Rasio 1996"}],"fun_headline_variants":["New 6.37-ms binary pulsar in M13 defies orbit trend","Fast folding finds pulsar in M13 that FFT missed","Wide-orbit millisecond pulsar discovered in globular M13","FFA re-detects 87 known pulsars, nabs new one in M13","M13I: wide-orbit pulsar with unexpected eccentricity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"That the new pulsar M13I is a member of globular cluster M13, inferred only from sky position (11.7 arcseconds from the center) and dispersion measure (29.58 pc $cm^{-3}$), with no independent distance or proper-motion measurement.","fun_headline_variants_meta":{"raw":{"variants":["New 6.37-ms binary pulsar in M13 defies orbit trend","Fast folding finds pulsar in M13 that FFT missed","Wide-orbit millisecond pulsar discovered in globular M13","FFA re-detects 87 known pulsars, nabs new one in M13","M13I: wide-orbit pulsar with unexpected eccentricity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000548,"raw_usage":{"total_tokens":2676,"prompt_tokens":1059,"completion_tokens":1617,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":675,"completion_tokens_details":{"reasoning_tokens":1517}},"tokens_in":675,"tokens_out":1617,"duration_ms":11271,"temperature":1.0,"reasoning_tokens":1517,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:14:29.195658+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A proper-motion or parallax measurement of M13I that disagrees with M13's cluster motion would falsify the physical association; alternatively, an independent blind re-reduction of the same FAST observations with a public FFA implementation that fails to reproduce M13I's signal or the reported re-detection counts would falsify the survey result.","supporting_citations":[],"review_version":1}