{"id":"6d8d4754-fead-4832-abe5-fa6078bd6a4b","arxiv_id":"2504.16293","paper_version":2,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":2,"one_line_summary":"CHAMPSS, a daily all-sky stacking search built on the CHIME/FRB datastream, discovered 11 new pulsars in a commissioning survey and is projected to reach about 30 microjansky at full scale.","lead":"This paper describes a new pulsar survey that reuses CHIME's fast radio burst data to watch the whole northern sky every day, stacking repeated observations to find faint and intermittent pulsars. It reports eleven newly discovered pulsars from a two-month commissioning run and projects that the full survey will reach about 30 microjansky sensitivity.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Flux calibration in Sec. 6.1 is the load-bearing uncertainty: a >2x systematic in Eq. 8 would invalidate the abstract's S600 > 0.1 mJy floor; the discovery claim itself is unaffected.","rationale":"The reader's weakest_assumption and my own analysis converge on the same point: Eq. 8's SEFD-based flux calibration is the least secure link in the quantitative version of the central claim. The eleven discoveries themselves are supported by folded profiles, timing residuals, known-source sifting, and CHIME/Pulsar follow-up, so I would not reject the paper. However, the abstract asserts a hard lower bound (S600 > 0.1 mJy) that could fail under a plausible 2–3x systematic calibration error, and the Table 1 S600 entries carry no uncertainties. The paper's own Section 6.2 check validates the sensitivity model against known pulsars but does not test Eq. 8's absolute scale. A cross-calibration against ATNF S600 values for known pulsars in the same footprint is a direct, inexpensive check. I therefore recommend conditional acceptance: publish the discoveries and system description, but require either the cross-check or a softened, qualified statement of the flux-density floor. A second stated limitation, the Section 3.2 normalization bug that adds red noise at f ≲ 1 Hz, reduces commissioning sensitivity but does not threaten the discovery claim because it is acknowledged and scheduled for a fix; similarly, the Section 3.8 injection-based completeness study is deferred, but the known-pulsar detection comparison in Section 6.2 partially fills that role.","tokens_in":29471,"tokens_out":11366,"duration_ms":116138,"concrete_test":"For every known pulsar in the commissioning footprint with an independent S600 estimate in ATNF and an archived CHIME/Pulsar fold-mode profile, recompute S600 exactly as in Eq. 8 using the same frequency-averaged SEFD calibration, and compare the ratios S600(Eq. 8)/S600(ATNF). Repeat the calculation using per-frequency SEFD values instead of a single band-averaged value. If the median ratio deviates from unity by more than a factor of two, or if J2302+4807 drops below 0.1 mJy under the corrected scale, the abstract and Table 1 need revised flux values or explicit error bars; if the median ratio is within a factor of two and the weakest source remains above threshold, the flux-calibration concern is resolved.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The discovery claim—eleven previously unknown pulsars confirmed by folding, timing, and CHIME/Pulsar follow-up—is well supported, and I find no soundness problem there. The quantitative layer of the central claim, namely 'each with S600 > 0.1 mJy' in the abstract and the S600 column of Table 1, rests on Eq. 8 in Section 6.1. That equation converts uncalibrated CHIME/Pulsar profiles to flux density using a frequency-averaged SEFD from calibrator monitoring, and it makes two unchecked assumptions: the off-pulse mean of the uncalibrated profile faithfully measures SEFD, and a single band-averaged SEFD is representative for steep-spectrum pulsars across 400–800 MHz. No uncertainties are quoted for any S600 value, and the faintest new pulsar (J2302+4807, 0.24 mJy) is only a factor of 2.4 above the abstract threshold. A correlated calibration error of order 2–3x would therefore push at least one source below 0.1 mJy and change the reported S600 range, though it would not affect the existence of the discoveries. The Section 6.2 predicted-versus-detected S/N comparison validates the survey sensitivity model but does not directly validate the absolute flux scale of Eq. 8 for the new sources.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents CHAMPSS, a pulsar periodicity search built on the CHIME/FRB intensity datastream. It describes the formation of quasi-tracking beams from adjacent static beams, position-dependent dedispersion and DM limits, power-spectrum generation, incoherent stacking over many days, candidate clustering, and follow-up via folding and timing. In a commissioning survey covering roughly 6% of the northern sky, the pipeline discovered 11 previously unknown isolated pulsars with spin periods 0.24-1.46 s and estimated S600 flux densities of 0.24-0.89 mJy, all confirmed by folded profiles and presented with timing residuals. The paper also validates the sensitivity model against known pulsars and reports three new pulsars with DMs in excess of NE2001/YMW16 predictions, including a scattered pulsar near the Cygnus region.","tokens_in":29713,"tokens_out":9741,"duration_ms":92247,"significance":"If the quantitative flux scale is confirmed, this is an important result: CHAMPSS demonstrates a new survey mode that revisits the full northern sky daily and reaches faint flux densities by incoherently stacking power spectra, with a real-time pipeline and a public codebase. The core discovery claim is credible because candidates were confirmed by folding and timing rather than by the search statistic alone, and the DM-excess pulsars provide falsifiable input for Galactic electron models. The paper's main weakness is the absolute calibration of the S600 values, which carries no quoted uncertainties and is load-bearing for the abstract's flux floor. The forecast sensitivity of the full survey is plausible but not yet demonstrated.","major_comments":[{"comment":"The S600 values are load-bearing for the abstract claim that all eleven new pulsars have S600 > 0.1 mJy and for the flux column of Table 1, but they are derived from uncalibrated CHIME/Pulsar profiles under two unchecked assumptions: that the off-pulse mean of the profile equals the SEFD and that a single frequency-averaged SEFD is representative across 400-800 MHz. No uncertainties are quoted for any S600 value, and the faintest source (J2302+4807, 0.24 mJy) is only a factor of 2.4 above the 0.1 mJy threshold. A correlated calibration error of order 2-3x would therefore invalidate the abstract's flux floor and change the reported range, even though it would not affect the existence of the discoveries. Please provide a systematic error budget for Eq. (8), validate the SEFD assumptions with data, or soften the abstract and table claims to explicitly reflect the calibration uncertainty.","section":"§6.1, Eq. (8), Table 1"},{"comment":"The text states that only 6 of the 11 pulsars have sufficient follow-up observations to derive timing solutions, yet Table 1 lists spin period and period-derivative values with uncertainties for all 11 pulsars and Figure 20 shows timing residuals for all 11. Please reconcile this inconsistency: either all 11 pulsars are phase-connected and timed, in which case the '6' is incorrect, or only 6 are, in which case the table and residual figure overstate the confirmation status. This matters because the paper's central claim is that all eleven sources are confirmed pulsars.","section":"§6.1, Table 1, Fig. 20"}],"minor_comments":[{"comment":"The pulsar labeled 'J2118+5001' in Figure 14 does not match 'J2118+5143' in Table 1; please correct the figure or the table.","section":"Figure 14 and Table 1"},{"comment":"The abstract says the commissioning data were searched 'over two months', while §6 describes a realtime survey running from October 2023 to June 2024 with each beamrow recorded for over a month; please clarify whether 'two months' refers to the stacking interval rather than the survey duration.","section":"Abstract and §6"},{"comment":"The time resolution is given as 0.98304 ms, 0.98306 ms, and 0.983 ms in different places; please use a single consistent value.","section":"§2.1 and §3.2"},{"comment":"Equation (4) is typeset incorrectly in the text ('f t 1≈'); it should read phi(t) = f0 t1 + (1/2) fdot t1^2, with t1 = t - tref.","section":"§4.2, Eq. (4)"},{"comment":"The 'good agreement' between predicted and detected S/N would be easier to assess if the plot included the number of known pulsars used and a quantitative measure of scatter or bias; as presented, the plot alone is not fully sufficient to validate the sensitivity model.","section":"§6.2, Fig. 15"},{"comment":"The 5-sigma and 6-sigma thresholds are quoted without a discussion of the number of independent trials; a sentence quantifying the effective trial factor and resulting false-positive rate would make the candidate filtering description more complete.","section":"§3.9"},{"comment":"In the candidate plot description, 'O Text Text field' appears to be a typo; it should read 'O Text field'.","section":"§8.1"}],"recommendation":"major_revision","confidential_remarks":"For the editor: this is a solid instrument and discovery paper, and I see no indication of circularity in the discovery claims — the candidates were confirmed by folding and timing. The main issue is the unsupported absolute flux scale in Eq. (8) and the inconsistency about the number of timing solutions; both are fixable within the scope of a revision. The discoveries themselves appear secure."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the discovery claim holds up. Eleven previously unknown pulsars, found in a two-month commissioning patch and confirmed by folding and phase-connected timing with small residuals, are real. The survey architecture is the actual novelty: it turns the CHIME/FRB static-beam datastream into a quasi-tracking, daily, all-sky periodicity search with incoherent power-spectrum stacking. Individually the pieces (FDMT dedispersion, harmonic summing, DBSCAN clustering, radiometer sensitivity estimates) are standard, but the combination on this datastream is new, and the code is on GitHub. The three DM-excess pulsars are a useful scientific payoff, especially J2108+5001 in the Cygnus region.\n\nSoft spots, in proportion. The one I would push on is the S600 column. Equation 8 converts uncalibrated CHIME/Pulsar profiles to flux density using a frequency-averaged SEFD from calibrator monitoring, with the off-pulse mean standing in for SEFD, and no uncertainties are quoted. The abstract's floor of S600 > 0.1 mJy and the faintest source at 0.24 mJy are only a factor of 2-3 above that floor, so a correlated calibration systematic of that size could drop at least one source below threshold and change the reported flux range. The existence of the pulsars does not depend on this calibration; folding and timing are calibration-independent. The Section 6.2 predicted-versus-detected S/N comparison validates the survey sensitivity model, not the absolute flux scale of Eq. 8. So this is a real caveat, not a fatal one. I would ask for per-source SEFD uncertainties and a propagation into S600 before publication.\n\nMinor things: the 30 uJy full-survey sensitivity is projected rather than measured; the normalization bug in Section 3.2 is disclosed honestly and reduced commissioning sensitivity but does not undermine the discoveries; there are small internal inconsistencies, like \"J2302+48\" versus \"J2302+4807\" in the text, that a referee should catch.\n\nThis paper deserves peer review. The discoveries are useful now, the survey will matter for pulsar and Galactic electron-model work, and the caveats are fixable. I would cite it and I would bring it to a pulsar-focused reading group.","headline":"A solid, honest commissioning paper: the 11 new pulsars are real, the survey architecture is genuinely new, and the only place I would push back is the unquoted S600 flux calibration before trusting the 0.1 mJy floor.","tokens_in":30413,"tokens_out":2384,"would_cite":true,"duration_ms":24734,"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":"Stacking daily sky scans yields 11 new pulsars","keywords":["pulsar surveys","radio astronomy","power spectrum stacking","CHIME","dispersion measure","neutron stars","fast Fourier transform","all-sky survey"],"falsifier":"Re-observe the eleven pulsars with an independently calibrated telescope at 600 MHz and compare the measured flux densities with the values in the paper's table; a systematic offset larger than a factor of two in the same direction would falsify the quantitative sensitivity claim, as would showing that the weakest pulsar, J2302+4807, falls below 0.1 mJy under a re-derived calibration.","tokens_in":1985,"feed_emoji":"📡","tokens_out":4139,"duration_ms":84840,"temperature":0.7,"pith_summary":"CHAMPSS is a new way of searching for pulsars with a radio telescope that already watches most of the northern sky every day. Instead of using the telescope's fixed beams as single snapshots, the survey stitches adjacent beams into quasi-tracking beams, dedisperses each pointing to a sky-position-dependent maximum dispersion measure, and converts each observation into a power spectrum. Repeated visits to the same point on the sky are added together incoherently, so a pulsar too faint to show up in any one day can accumulate signal over weeks or months, while intermittent sources can be caught on days when they happen to pulse. In a commissioning survey covering about 6% of the sky over roughly two months, the pipeline found eleven previously unknown isolated pulsars with spin periods between 0.24 and 1.46 seconds and estimated flux densities between 0.24 and 0.89 mJy at 600 MHz. If the survey scales as planned to the full northern sky with more than a year of stacking per sightline, it reaches $\\lesssim 30\\,\\mu$Jy sensitivity and becomes the deepest all-sky pulsar survey to date.","feed_headline":"Stacking daily sky scans yields 11 new pulsars","feed_subtitle":"A survey that combines days of power spectra promises a 30-microjansky all-sky pulsar search.","key_machinery":"The load-bearing mechanism is the power-spectrum stack. Each pointing's dedispersed time series is Fourier transformed, the resulting power spectra are barycentrically corrected, red-noise corrected, cleaned of radio-frequency interference using the zero-dispersion-measure spectrum as a reference, and then incoherently summed across days. Because the summed noise follows a $\\chi^2_{2m}$ distribution, the significance of a real periodic signal grows as more days are added while noise fluctuations average down; harmonic summing over 1, 2, 4, 8, 16, and 32 harmonics collects power from narrow pulse profiles. A second mechanism is the pointing map, which sets the maximum searched dispersion measure from the larger of two Galactic electron models for each sightline and uses that limit to decide how finely to channelize the data, keeping the daily data volume manageable at 52 TB per day.","core_discovery":"The central discovery is that a survey built on the daily, all-sky intensity stream of a transit radio telescope can discover pulsars by combining stationary beams into quasi-tracking beams, dedispersing to position-dependent dispersion-measure limits, and summing power spectra across many days. The eleven new pulsars demonstrate the concept: all were confirmed by folding, six of them have fitted timing solutions, and their periods, dispersion measures, and flux densities are tabulated. Three of the pulsars have dispersion measures above what standard Galactic electron models predict for their sightlines, which the paper interprets as evidence that the survey will reveal where those models fail. The intended endpoint is a full northern-sky survey that stacks more than a year of data per pointing and reaches $\\lesssim 30\\,\\mu$Jy, which the paper argues will make it deeper than any previous all-sky pulsar survey.","pith_inferences":["The same stacking design could be transferred to other large-format, transit-style radio arrays, turning fixed beams into effective tracking surveys at a fraction of the data rate of tied-array beamforming.","If the flux-density calibration holds up, the survey's per-pointing depth will keep growing with the square root of the number of stacked days, so the final sensitivity will depend on maintaining stable radio-frequency-interference statistics and calibration over years.","The three dispersion-measure-excess pulsars hint that the survey will double as a wide-field probe of Galactic ionized structures such as H II regions, potentially tracing their geometry through many new pulsar sightlines.","A testable extension would be to run the same power-spectrum stacking and clustering machinery on shorter time intervals, which could catch weakly periodic sources that are not stable enough to appear in month-long stacks."],"forward_implications":["The full survey, stacking more than one year per pointing, reaches $\\lesssim 30\\,\\mu$Jy for all sightlines above a declination of $10^\\circ$ and away from the Galactic plane, making it the deepest all-sky pulsar survey.","Daily repeated observations give the survey sensitivity to intermittent pulsars—nulling, eclipsing, scintillating, or precessing sources—that single-visit surveys can miss.","New pulsars on under-searched sightlines, especially those with dispersion measures in excess of model predictions, will help refine Galactic electron density models and clarify the boundary between high-dispersion-measure Galactic pulsars and low-dispersion-measure fast radio bursts.","The survey is complementary to targeted coherent follow-up observations: CHAMPSS finds faint candidates in the power-spectrum stack, and the follow-up timing pipeline converts them into pulsars with measured positions, spin periods, and period derivatives.","The eleven commissioning discoveries, confirmed by folding and partly by timing solutions, show that the full pipeline from data acquisition to candidate confirmation works end to end."],"supporting_citations":[{"why":"Describes the CHIME/FRB instrument and datastream, including the 1024 stationary beams, 0.983 ms sampling, and L1 RFI masking that CHAMPSS taps into.","marker":"CHIME/FRB Collaboration et al. 2018"},{"why":"Establishes the statistical basis for incoherently summing power spectra, specifically that the summed noise follows a $\\chi^2_{2m}$ distribution.","marker":"van der Klis 1989"},{"why":"Supplies the harmonic-summing methodology and significance evaluation that the CHAMPSS search uses to detect periodic signals in power spectra.","marker":"Ransom et al. 2002"},{"why":"Provides the fast dispersion-measure transform algorithm used to dedisperse each pointing to many trial dispersion measures.","marker":"Zackay & Ofek 2017"},{"why":"One of the two Galactic electron density models used to set the maximum searched dispersion measure per sightline.","marker":"Cordes & Lazio 2002"},{"why":"The other Galactic electron density model used alongside NE2001 to define the position-dependent dispersion-measure search limits.","marker":"Yao et al. 2017"},{"why":"Provides the CHIME gain value and receiver-temperature context used in the sensitivity estimates for known pulsars.","marker":"Good et al. 2021"},{"why":"Documents earlier pulsar discoveries from the CHIME/FRB system and supplies the known pulsar PSR J2208+4610 used as a test case for the multi-day phase-coherent search.","marker":"Dong et al. 2023"},{"why":"Describes the tracking-beam follow-up system used to confirm candidates and produce timing solutions.","marker":"CHIME/Pulsar Collaboration et al. 2021"}],"fun_headline_variants":["Daily all-sky stack finds 11 new pulsars","Revisiting the sky daily nets 11 pulsars","Sky scans stacked daily yield 11 new pulsars","All-sky stacking pings 11 new pulsars"],"cache_read_input_tokens":32384,"weakest_assumption_plain":"The reported flux-density range depends on converting uncalibrated pulse profiles into janskys using the telescope's noise temperature measured from calibrator radio sources; if that calibration carries a systematic error larger than roughly a factor of two, some of the reported values could drop below the 0.1 mJy lower bound stated in the abstract, although the existence of the eleven pulsars themselves does not depend on this calibration.","fun_headline_variants_meta":{"raw":{"variants":["Daily all-sky stack finds 11 new pulsars","Revisiting the sky daily nets 11 pulsars","Sky scans stacked daily yield 11 new pulsars","All-sky stacking pings 11 new pulsars"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000877,"raw_usage":{"total_tokens":3864,"prompt_tokens":1089,"completion_tokens":2775,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":705,"completion_tokens_details":{"reasoning_tokens":2709}},"tokens_in":705,"tokens_out":2775,"duration_ms":18567,"temperature":1.0,"reasoning_tokens":2709,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T11:08:45.920940+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-observe the eleven pulsars with an independently calibrated telescope at 600 MHz and compare the measured flux densities with the values in the paper's table; a systematic offset larger than a factor of two in the same direction would falsify the quantitative sensitivity claim, as would showing that the weakest pulsar, J2302+4807, falls below 0.1 mJy under a re-derived calibration.","supporting_citations":[],"review_version":1}