{"id":"09100f59-d994-42e7-add3-e709260c36a3","arxiv_id":"2411.09135","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"Interstellar object impacts onto neutron stars are a feasible source of some non-repeating fast radio bursts, with rates, durations, and energies consistent with observations.","lead":"Fast radio bursts could be caused by interstellar objects smashing into neutron stars. The paper finds the expected collision rate and burst properties roughly match telescope observations, making this a plausible source for some one-off FRBs.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The headline rate comparison ignores the energy threshold: Eq. (4) counts all 1I-sized ISO-NS collisions, while the observed FRB rate is for bursts above 2e35 erg; folding in the paper's own size distributions leaves only ~2–9% of collisions above threshold, so the predicted observable rate falls…","rationale":"The reader's verdict CONDITIONAL is appropriate. I examined the central rate claim and found that the weakest step is not (only) the uncertain ISO density, but the mismatch between the rate being computed and the rate being compared. Eq. (4) counts all 1I-sized collisions, but the observed FRB rate is threshold-limited. The paper's own Section 4 concedes this, yet the abstract's 'comparable' phrasing overstates the support. A back-of-the-envelope correction using the paper's Table 1 size distributions suggests the observable above-threshold rate is ~10^5–10^6 Gpc^-3 yr^-1, about two orders of magnitude below the Bochenek rate. This does not mean the model is impossible; stronger fields, flatter size distributions, or higher ISO densities could close the gap, but it does mean the headline claim is not yet established by the calculations presented. The reader identified n_ISO as the weakest assumption; I consider the threshold mismatch more load-bearing because it is a systematic bias rather than an order-of-magnitude uncertainty. The concrete test would settle it with modest additional computation. Therefore I leave the verdict at CONDITIONAL; the revision should either produce the threshold-matched rate or soften the claim.","tokens_in":11832,"tokens_out":9018,"duration_ms":194449,"concrete_test":"Recompute the observable rate with an energy threshold: adopt a power-law ISO size distribution dN/dR ∝ R^{-q} (q = 3.5, normalized to n_ISO at R > 0.1 km), use Eq. (6) to convert R to E_iso at B_surface = 10^12 G and the fiducial parameters, set E_th = 2×10^35 erg, and compute the fraction of collisions exceeding E_th. Multiply R_obs from Eq. (4) by this fraction. If the resulting rate is below 10^6 Gpc^-3 yr^-1, the abstract's 'comparable' claim is not supported and the paper should be revised to state a conditional feasibility with an explicit threshold-matched rate.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the ISO-NS collision rate is comparable to the observed FRB rate rests on Eq. (4), R_obs = f R_col ≈ 10^7 Gpc^-3 yr^-1, which is then compared to the Bochenek et al. (2020) rate of (7+9-6)×10^7 Gpc^-3 yr^-1. However, R_col in Eq. (3) is the rate of collisions with all ISOs of size ≳ 1I/'Oumuamua (R ≳ 0.1 km), whereas the Bochenek rate counts only bursts with isotropic-equivalent energy ≥ 2×10^35 erg. The paper's own Eq. (6) and Section 3 imply that reaching 10^35 erg requires R ≳ 0.5 km at typical pulsar fields, with E ∝ R^4. Using the size distributions in Table 1 (q = 2.5–3.5), the fraction of ISOs with R > 0.5 km relative to R > 0.1 km is (0.1/0.5)^(q-1) ≈ 0.02–0.09. Thus only a few percent of the collisions counted in Eq. (3) would be detectable in the Bochenek sample. The paper acknowledges this in Section 4 ('not directly comparable ... beyond the scope of this work'), yet the abstract and conclusion still assert comparability. This is a systematic overcount, not a mere uncertainty: with a rough threshold correction the predicted observable rate drops to ~10^5–10^6 Gpc^-3 yr^-1, roughly two orders of magnitude below the observed rate, which undermines the central feasibility claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes that collisions between interstellar objects (ISOs) and neutron stars (NS) could produce a subset of one-off fast radio bursts (FRBs). The authors calculate the ISO-NS collision rate using a Solar-neighborhood ISO number density, gravitational focusing, and counts of neutron stars and galaxies, obtaining a beaming-corrected rate of ~10^7 Gpc^-3 yr^-1, which they compare to the observed FRB rate. They then use the Dai et al. (2016) emission model with updated material strengths to show that observed FRB durations and energies correspond to planetesimal sizes of roughly 0.4–10 km, that the FRB energy distribution is consistent with a size distribution exponent q = 2.5–3.5, and that sub-bursts like those in FRB 200428 could be explained by binary ISOs. Finally, they argue that the ISO-NS collision rate should increase with cosmic time, matching the claimed redshift evolution of the FRB rate.","tokens_in":12137,"tokens_out":5496,"duration_ms":59355,"significance":"The paper is valuable as a concrete, quantitative exploration of a progenitor channel that is often mentioned but rarely calculated. It assembles independent inputs (ISO density, NS kick velocities, updated tensile/compressive strengths, power-law size distributions) rather than fitting to the FRB data, and it makes several testable predictions: lack of sub-burst structure from monolithic ISOs, a limit of two sub-bursts if binaries are the cause, a host-galaxy morphology dependence, and a cosmic evolution that decouples from the star formation rate. If the central rate comparison can be corrected for the energy threshold (a fix that is well within the paper's scope), the mechanism remains a plausible contributor to one-off FRBs. The paper is also transparent about the dominant uncertainty in the ISO number density, although it does not propagate uncertainties into the headline comparison.","major_comments":[{"comment":"The central rate comparison is not apples-to-apples. R_col in Eq. (3) counts collisions with all ISOs larger than 1I/'Oumuamua (R ≳ 0.1 km), while the Bochenek et al. (2020) rate used for comparison is for bursts with isotropic-equivalent energy ≥ 2×10^35 erg. Using the paper's own scalings (E ∝ L Δt ∝ R^4 from Eq. (6) and the text preceding it), a burst of 2×10^35 erg requires R ≳ 0.5 km at typical pulsar field strengths. For the size distributions in Table 1 (q = 2.5–3.5), the fraction of ISOs with R > 0.5 km relative to those with R > 0.1 km is (0.1/0.5)^(q-1) ≈ 0.02–0.09. Applying this to Eq. (4) reduces the predicted observable rate to roughly 2×10^5–9×10^5 Gpc^-3 yr^-1, about two orders of magnitude below the observed (7+9−6)×10^7 Gpc^-3 yr^-1. The paper explicitly acknowledges this mismatch in Section 4 ('not directly comparable... beyond the scope of this work'), yet the abstract and conclusion still state that the ISO-NS collision rate is 'comparable with' or 'consistent with' the observed FRB rate. This is a systematic overcount, not a random uncertainty. The authors should either perform a threshold-corrected rate calculation (even with a broad range of assumptions on the size distribution and magnetic field distribution) or carefully restrict the claim to the raw collision rate and state explicitly that the predicted observable FRB rate is substantially lower.","section":"§2 (Eq. 4) and §4"},{"comment":"The rate estimate in Eq. (3) stacks several factors (nISO, f, NNS, nGal, v∞) as point values, and the conclusion of 'within order-of-magnitude uncertainties' is never quantified. The dominant input, nISO ~ 10^15 pc^-3, is derived from a single detection and is acknowledged to carry an order-of-magnitude uncertainty, but no sensitivity analysis is presented. Since the headline claim is a rate comparison, the paper should show how R_obs varies when nISO is varied by 1–2 orders of magnitude, and when f is varied over a plausible range (e.g., 10^-3–10^-1). Without such a propagation, the reader cannot assess whether the central claim is robust even after the threshold correction in the first major comment is applied.","section":"§2 (Eq. 3)"}],"minor_comments":[{"comment":"The conversion from CHIME pulse width to planetesimal radius assumes s = 500 Pa and ρ = 3 g cm^-3 without displaying the sensitivity to these values; a brief demonstration that the inferred size range 0.4–10 km is robust to a factor of a few in these parameters would strengthen the figure.","section":"§3, Fig. 1"},{"comment":"The comparison of the predicted energy-distribution exponent γ = (3+q)/4 with the observed γ = 1.3+0.7−0.4 from Shin et al. (2023) is made using a sample that includes repeating FRBs and summed multi-pulse energies. The text acknowledges this but states without quantitative support that 'we do not expect this to change the value of γ by a significant amount.' Please either justify this expectation or cite an analysis that isolates non-repeating single pulses.","section":"§3, Table 1"},{"comment":"The argument that the FRB rate increases with cosmic time is based on a set of references, but the literature is mixed, with some analyses finding consistency with the star formation rate at low redshift. Please acknowledge the ongoing debate and clarify what is required for the ISO-NS scenario to match the observed evolution (e.g., a specific delay-time distribution or ISO production history).","section":"§4"},{"comment":"The duration expression is quoted from Colgate & Petschek (1981) with updated parameters, but it would be helpful to include a one-sentence derivation or a reference to the equation number in that paper, as the functional form (Δt ∝ R^4/3) is non-trivial and central to the duration-size mapping.","section":"§3, Eq. (5)"},{"comment":"The assumption that ISOs retain the same spatial distribution as their parent stars is cited to Hopkins et al. (2024, arXiv:2402.04904). It may be worth noting that this is a preprint at the time of writing, so that the reader can weigh the strength of the assumption accordingly.","section":"§2"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern about the energy-threshold mismatch is real and lands directly on the abstract's central claim. However, the paper's own Section 4 already identifies the mismatch, and the fix (a threshold-corrected rate estimate or a carefully qualified claim) is within the scope of the manuscript. The paper has several redeeming features—transparent order-of-magnitude reasoning, no fitting to FRB data, and concrete testable predictions—so I believe major revision is appropriate rather than rejection. The authors should also be encouraged to provide a simple uncertainty propagation for the rate comparison, since the current point estimate without error bars is difficult to evaluate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis paper is worth a look if you care about FRB progenitors or interstellar object demographics. The mechanism itself is not new — Geng & Huang and Dai et al. proposed planetesimal–neutron star collisions years ago, and Siraj & Loeb already applied it to ISOs. What this paper adds is a concrete quantitative mapping: updated tensile/compressive strengths (rubble-pile values rather than Colgate & Petschek's 10^9 Pa), a duration–size relation calibrated to CHIME Catalog 1, and a derived relation between the ISO size distribution exponent and the FRB energy distribution exponent. The comparison with observed FRB durations and energies is genuinely useful, and the suggestion that sub-bursts like FRB 200428 might come from binary ISOs is interesting, even if speculative. The paper is also honest about its biggest uncertainties — the ISO number density from a single detection and the lack of constraints on the ISO size distribution.\n\nThe soft spot is the headline rate claim. The abstract and conclusion say the ISO–NS collision rate is 'comparable' to the observed FRB rate, but the rate in Eq. (4) counts all collisions with 1I-sized ISOs, while the Bochenek rate applies only to bursts above ~2×10^35 erg. The paper's own Eq. (6) and size distributions imply that only a few percent of those collisions exceed threshold, so the predicted observable rate is ~10^5–10^6 Gpc^-3 yr^-1, about two orders of magnitude low. The paper acknowledges this in Section 4 and says a proper comparison is beyond scope, but then the abstract should not claim comparability. That is a presentation problem as much as a physical one — the mechanism might still contribute to some FRBs, but the central feasibility claim as stated is not supported.\n\nAlso, the rate estimate stacks n_ISO, f, N_NS, and galaxy density without propagated errors, and the cosmic-time argument (long-lived progenitors building up over time) is qualitative. None of this is fatal, but it needs to be tightened.\n\nVerdict: this deserves serious peer review, but the authors should either threshold-match the rates or reframe the headline as an upper limit. I would send it to a referee and expect a revision. I'd cite it for the updated strength mapping and the duration/energy relations.\n\nBest.","headline":"A useful feasibility study with an overstated headline rate: the energy-threshold mismatch cuts the predicted observable FRB rate by roughly two orders of magnitude, though the paper openly flags it in Section 4.","tokens_in":12813,"tokens_out":1447,"would_cite":true,"duration_ms":18304,"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":"Interstellar objects colliding with neutron stars are a feasible source of one-off fast radio bursts.","keywords":["fast radio bursts","interstellar objects","neutron stars","planetesimals","tidal disruption","collision rates","FRB progenitors","transient radio events"],"falsifier":"Measure the volumetric rate of single-pulse FRBs with a sample that cleanly excludes repeaters and compare it with $f\\,\\Gamma\\,N_{\\rm NS}\\,n_{\\rm Gal}$: if the observed rate exceeds the prediction even when every collision is assumed to produce a beamed burst ($f=1$), then ISO-neutron-star collisions cannot supply the FRB population. Alternatively, a precise measurement of the local ISO number density more than an order of magnitude below $10^{15}\\,\\mathrm{pc}^{-3}$ would remove the rate match, since $R_{\\rm obs}\\propto n_{\\rm ISO}$.","tokens_in":11523,"feed_emoji":"☄️","tokens_out":9158,"duration_ms":86740,"temperature":0.7,"pith_summary":"This paper argues that collisions between interstellar objects (ISOs), planetesimals ejected from their home planetary systems, and neutron stars occur often enough, and with the right burst properties, to be a feasible source of the observed one-off fast radio bursts (FRBs). The expected cosmological collision rate is $\\sim 10^7\\,\\mathrm{Gpc}^{-3}\\,\\mathrm{yr}^{-1}$ after beaming, within an order of magnitude of the observed rate of $(7^{+9}_{-6})\\times 10^7\\,\\mathrm{Gpc}^{-3}\\,\\mathrm{yr}^{-1}$. The paper also matches observed FRB durations to planetesimal radii of roughly $0.4$–$10\\,\\mathrm{km}$, consistent with the two known ISOs, and shows the FRB energy distribution is consistent with Solar System planetesimal size distributions. If right, FRBs would become a way to study interstellar objects far beyond the Solar System, and the FRB rate's cosmic evolution would trace accumulated neutron-star and ISO populations rather than star formation.","feed_headline":"Interstellar rocks crashing into neutron stars could explain FRBs","feed_subtitle":"Collision rates match observed bursts; durations and energies match known interstellar objects.","key_machinery":"The argument is carried by three linked pieces. (1) Gravitational focusing inflates the neutron-star collision cross-section to $\\sigma\\sim10^9\\,\\mathrm{km}^2$, giving a per-star rate $\\Gamma=n_{\\rm ISO}v_\\infty\\sigma\\sim10^{-7}\\,\\mathrm{yr}^{-1}$. (2) The Dai et al. (2016) emission model translates a planetesimal's radius $R$ into a burst: duration $\\Delta t\\propto R^{4/3}$ and isotropic energy $E_{\\rm iso}\\propto R^4$, with beaming $f\\sim10^{-2}$, so observed millisecond bursts correspond to $0.4$–$10\\,\\mathrm{km}$ bodies. (3) The assumption that ISOs are weak rubble-pile material (tensile strength $\\sim500\\,\\mathrm{Pa}$) means they fragment into a homogeneous stream during infall, yielding single pulses; sub-burst FRBs are then attributed to binary ISOs.","core_discovery":"The central claim is that interstellar objects are a viable reservoir of planetesimals for the proposed planetesimal–neutron-star FRB mechanism. Using a local ISO number density of $n_{\\rm ISO}\\sim10^{15}\\,\\mathrm{pc}^{-3}$, a relative speed of $v_\\infty\\sim100\\,\\mathrm{km\\,s^{-1}}$, and a gravitational-focusing cross-section of $\\sim10^9\\,\\mathrm{km}^2$, the per-neutron-star encounter rate is $\\Gamma\\sim10^{-7}\\,\\mathrm{yr}^{-1}$. Extrapolating to $\\sim10^9$ neutron stars per Milky-Way-type galaxy and $\\sim10^7$ galaxies per Gpc$^3$, and applying a beaming factor $f\\sim10^{-2}$, gives an observable rate $R_{\\rm obs}\\sim10^7\\,\\mathrm{Gpc}^{-3}\\,\\mathrm{yr}^{-1}$, comparable to the observed FRB rate. With updated tensile strength $s\\sim500\\,\\mathrm{Pa}$, the emission model maps the observed FRB duration distribution to planetesimal radii $\\sim0.4$–$10\\,\\mathrm{km}$, and the observed isotropic-energy power-law index $\\gamma\\approx1.3$ is consistent with known planetesimal size distributions via $\\gamma=(3+q)/4$ for $q\\approx2.5$–$3.5$.","pith_inferences":["Editorial inference: because the predicted observable rate scales linearly with $n_{\\rm ISO}$, a future, well-characterised measurement of the local interstellar-object number density from a wide-field survey would tighten or overturn the rate matching directly.","Editorial inference: if $E_{\\rm iso}\\propto R^4$, then the high-energy tail of the FRB energy distribution can be inverted to measure the size distribution of ejected planetesimals over a much larger size range than any direct Solar System census.","Editorial inference: the binary-ISO interpretation of two-subburst FRBs predicts that the subburst multiplicity distribution should fall sharply after two subbursts; counting subbursts in new CHIME catalogs offers a direct test, and the 29 ms gap in FRB 200428 can be used to estimate binary separation.","Editorial inference: the paper's cosmic-time argument implies FRBs from this channel should be hosted by older, more massive stellar populations than magnetar-driven FRBs, a distinction host-galaxy samples can test."],"forward_implications":["If the central claim holds, one-off FRBs can be produced by ordinary interstellar debris hitting neutron stars, so the FRB population need not be dominated by exotic magnetar-like engines.","The observed duration range of CHIME FRBs implies impacting planetesimals with radii between roughly 400 m and 10 km, sizes consistent with 'Oumuamua and Borisov.","The energy distribution of FRBs, with index $\\gamma\\approx1.3$, matches the size distributions of Solar System asteroid and trans-Neptunian populations through $\\gamma=(3+q)/4$.","ISO-neutron-star collisions cannot explain repeating FRBs, and the mechanism cannot be the only FRB source; the paper expects a one-off subset.","Because neutron stars and ISOs accumulate over cosmic time, the collision rate should rise with time rather than track the star-formation rate, matching the observed redshift evolution of the FRB rate."],"supporting_citations":[{"why":"Supplies the Solar-neighbourhood ISO number density $n_{\\rm ISO}\\sim10^{15}\\,\\mathrm{pc}^{-3}$ that anchors the collision-rate calculation.","marker":"(Do et al. 2018)"},{"why":"Provides the observed FRB volumetric rate $(7^{+9}_{-6})\\times10^7\\,\\mathrm{Gpc}^{-3}\\,\\mathrm{yr}^{-1}$ and the FRB 200428 energy baseline the paper compares against.","marker":"(Bochenek et al. 2020)"},{"why":"Proposes the planetesimal-neutron-star collision as an FRB mechanism whose emission model the paper adopts.","marker":"(Geng & Huang 2015)"},{"why":"Supplies the radiation mechanism linking planetesimal radius to FRB duration and luminosity.","marker":"(Dai et al. 2016)"},{"why":"Provides the tidal disruption cascade and the timescale formula $\\Delta t\\propto R^{4/3}$ used to convert burst durations to planetesimal sizes.","marker":"(Colgate & Petschek 1981)"},{"why":"Gives the beaming factor $f\\sim10^{-2}$ and the asteroidal interpretation of FRB 200428 that the paper revisits.","marker":"(Dai 2020)"},{"why":"Gives the observed FRB isotropic energy distribution exponent $\\gamma\\approx1.3$ compared with planetesimal size-distribution predictions.","marker":"(Shin et al. 2023)"},{"why":"Provides the FRB pulse-duration sample used to compare observed durations with planetesimal radii.","marker":"(CHIME/FRB Collaboration et al. 2021)"}],"fun_headline_variants":["Space rocks hitting neutron stars produce FRB-like bursts","Tiny interstellar objects may trigger fast radio bursts","Colliding space rocks and neutron stars match FRB rates","FRB mystery: interstellar objects collide with neutron stars","Neutron star hits by interstellar objects explain FRBs"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the interstellar-object density measured near the Sun, about $10^{15}\\,\\mathrm{pc}^{-3}$, also applies around neutron stars across the Universe; the predicted collision rate is directly proportional to this number.","fun_headline_variants_meta":{"raw":{"variants":["Space rocks hitting neutron stars produce FRB-like bursts","Tiny interstellar objects may trigger fast radio bursts","Colliding space rocks and neutron stars match FRB rates","FRB mystery: interstellar objects collide with neutron stars","Neutron star hits by interstellar objects explain FRBs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000239,"raw_usage":{"total_tokens":1551,"prompt_tokens":1021,"completion_tokens":530,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":637,"completion_tokens_details":{"reasoning_tokens":453}},"tokens_in":637,"tokens_out":530,"duration_ms":5767,"temperature":1.0,"reasoning_tokens":453,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T21:00:05.549584+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the volumetric rate of single-pulse FRBs with a sample that cleanly excludes repeaters and compare it with $f\\,\\Gamma\\,N_{\\rm NS}\\,n_{\\rm Gal}$: if the observed rate exceeds the prediction even when every collision is assumed to produce a beamed burst ($f=1$), then ISO-neutron-star collisions cannot supply the FRB population. Alternatively, a precise measurement of the local ISO number density more than an order of magnitude below $10^{15}\\,\\mathrm{pc}^{-3}$ would remove the rate match, since $R_{\\rm obs}\\propto n_{\\rm ISO}$.","supporting_citations":[{"cited_title":"A., & Petschek, A","cited_arxiv_id":null,"evidence_quote":"Provides the tidal disruption cascade and the timescale formula $\\Delta t\\propto R^{4/3}$ used to convert burst durations to planetesimal sizes."}],"review_version":1}