{"id":"f8951034-f29f-41fa-90d9-e00463b4e289","arxiv_id":"2608.09336","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Pulsar and supernova-remnant models yield TeV electron features with effective sharpness beta_eff of roughly 2 to 3.6, so a sharper, stable edge would favor exotic origins.","lead":"Nearby pulsars and supernova remnants can explain the high-energy electron features seen by DAMPE and H.E.S.S., but only with moderate sharpness. A truly sharp, stable edge in the spectrum would be hard to explain with ordinary astrophysics and would point toward dark matter or other new physics.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The diagnostic conclusion lacks a demonstrated upper bound on pulsar spectral sharpness: beta_eff in Table I is computed for one ISRF/B environment, and a higher-B/lower-ISRF path could restore a sharp cooling edge.","rationale":"The paper is a careful phenomenological study, and the qualitative broadening of the cooling edge from stochastic inverse-Compton losses is consistent with previous work and internally plausible. The reader's weakest assumption correctly identifies the environmental dependence. I go further: the missing parameter scan affects exactly the discriminating quantity, beta_eff, because the maximum sharpness achievable in ordinary pulsar propagation is left unbounded. The baseline Table I already reaches beta_eff=3.61, which is a super-exponential falloff; a plausible increase in B or decrease in ISRF could push this higher, weakening the paper's central diagnostic. However, this is a conditional concern rather than a fatal flaw: the stochastic-IC calculation itself is a genuine step, the authors explicitly limit their quantitative conclusions, and the qualitative statement that very sharp edges are hard for the specific Geminga-like environment remains valid. The verdict should remain CONDITIONAL, with the condition being a robustness scan over the local radiation and magnetic environment.","tokens_in":13746,"tokens_out":8873,"duration_ms":99706,"concrete_test":"Repeat the burst-like stochastic-IC Monte Carlo of Sec. III with B=10 uG (U_B=2.5 eV/cm^3) and the Popescu et al. ISRF scaled down by a factor of three, for T=100, 300, and 1000 kyr, and recompute beta_eff exactly as in Table I. If any beta_eff exceeds ~5, the claimed difficulty of producing sharp edges with ordinary pulsar propagation is not robust; also report the required W_e for those cases to check whether the energetics would already exclude them.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim of Sec. V is that a very sharp and stable edge-like feature would be difficult to accommodate with ordinary pulsar propagation. The quantitative support is Table I, where burst-like pulsars with stochastic inverse-Compton cooling give beta_eff from 1.29 to 3.61. These values are obtained for a single environmental setup: a uniform B=3 uG (U_B=0.225 eV/cm^3) and the path-averaged Geminga ISRF (rho_opt=0.435, rho_IR=0.504, rho_CMB=0.260 eV/cm^3; Sec. III). The stochastic broadening of the cooling edge is controlled by the fraction of energy lost in large discrete IC scatterings relative to continuous synchrotron losses. If a pulsar's propagation path has a stronger magnetic field and/or a weaker radiation field, synchrotron dominates and the cooling becomes more continuous, so the sharp cutoff at E_br is partially restored. The paper does not scan over B, ISRF normalization, injection index gamma, or Ecut, and therefore does not establish an upper bound on beta_eff for ordinary pulsar propagation. The highest baseline value, beta_eff=3.61 at T=1000 kyr, lies well above the continuous-injection value (~2) that the paper itself describes as 'sharp'; a moderate environmental variation could push it higher. Without this bound, the statement that a sharp edge would strengthen the case for exotic origins is not quantitatively supported.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript revisits the interpretation of TeV-scale features in the cosmic-ray all-electron spectrum. It constructs a diffuse electron/positron background from a retained ensemble of 64 GALPROP configurations, normalized at 100 GeV, and shows that the ensemble remains smooth without narrow structures. For local sources, it solves the diffusion-loss Green function for burst-like and continuous injection histories, and implements a Monte Carlo that treats synchrotron losses continuously while sampling individual inverse-Compton scatterings with the Klein-Nishina kernel on a path-averaged ISRF. The empirical post-peak sharpness parameter beta_eff is tabulated for the propagated templates. The paper then normalizes three illustrative source templates to the two highest-energy DAMPE points and concludes that smooth or cooling-broadened features are compatible with pulsars, whereas a very sharp and stable edge would be more difficult to accommodate with ordinary pulsar propagation and would strengthen the case for alternative origins.","tokens_in":14050,"tokens_out":7528,"duration_ms":84165,"significance":"The paper is valuable because it formulates the astrophysical-versus-exotic diagnostic in terms of a measurable shape parameter, uses a physically motivated GALPROP background rather than an ad hoc broken power law, and includes stochastic inverse-Compton cooling with a reasonably detailed treatment of the Klein-Nishina kernel and the local ISRF. The Monte Carlo description is transparent, convergence checks are reported, and the illustrative character of the DAMPE normalization is explicitly acknowledged. If the central claim survives a wider parameter-space exploration, the result would be a useful guide for interpreting upcoming DAMPE, HERD, and other high-energy lepton measurements. The main weakness is that the headline conclusion is currently supported by a single environmental setup, so the paper does not yet establish the claimed upper bound on spectral sharpness.","major_comments":[{"comment":"The central Sec. V claim that a very sharp and stable edge would be difficult to accommodate with ordinary pulsar propagation is supported only by beta_eff values computed for one environmental setup: B = 3 uG and the path-averaged Popescu et al. ISRF along the Geminga line of sight. Because the stochastic broadening of the cooling edge is controlled by the fraction of energy lost in discrete inverse-Compton scatterings relative to continuous synchrotron losses, a pulsar propagating through a region with stronger magnetic field and/or weaker radiation field would cool more continuously and could exhibit a sharper edge. The paper does not scan over B, ISRF normalization, injection index, or Ecut, so Table I does not provide an upper bound on beta_eff for ordinary pulsar propagation. I request a parameter scan (for example B = 1-10 uG and ISRF scaled by 0.1-3) with a discussion of how the conclusion shifts, or a careful restriction of the conclusion to the local Geminga-like environment.","section":"Sec. III and Table I"},{"comment":"The diagnostic statement 'a very sharp and stable edge-like feature' is not given a quantitative definition. The manuscript does not specify the beta_eff threshold above which a pulsar origin would be excluded, and beta_eff itself depends on the fitting window, the smoothing kernel, and the fixed p values used in Table I. The paper also does not quantify 'stable', although Table I shows beta_eff varying from 1.29 to 3.61 with source age. Without a target shape or a comparison of predicted beta_eff distributions against future data, the central claim is not falsifiable. The authors should either define an explicit quantitative criterion or soften the conclusion to a qualitative remark.","section":"Sec. V"},{"comment":"The comparison between the dashed continuous-loss curves and the solid stochastic-IC curves does not isolate the effect of stochasticity: the two calculations also differ in the energy dependence of the loss rate because the Monte Carlo uses a Klein-Nishina-reduced inverse-Compton rate while the dashed curves use b(E) = b0 E^2. The text acknowledges this in Sec. III, but Sec. V attributes the broadened decline to stochastic inverse-Compton cooling. To make the mechanism claim clean, the authors should compare the stochastic run against a Monte Carlo with the same Klein-Nishina loss rate treated continuously, or explicitly state throughout the conclusions that the broadening is a combined stochastic-plus-Klein-Nishina effect.","section":"Fig. 3"}],"minor_comments":[{"comment":"The phrase 'common propagation framework' in the abstract is stronger than what is implemented: the diffuse background uses GALPROP models, while the discrete-source calculation uses a fixed simple diffusion-loss setup with D0 = 4.3e28 cm^2/s, delta = 0.415, and b0 = 1e-16 GeV^-1 s^-1. Please clarify the relation between these two propagation descriptions.","section":"Abstract and Sec. II"},{"comment":"The label 'cutoff fit' in the legend could be confused with the dashed continuous-loss curves; consider renaming the empirical fits to 'empirical cutoff fit' or similar.","section":"Fig. 3 caption"},{"comment":"The heading 'T[kyr]' in the continuous-injection group is followed by values such as 0.1 and 0.3; adding units explicitly in the table header or a footnote would avoid ambiguity.","section":"Table I"},{"comment":"There is a missing space in 'Kobayashiet al.' in Sec. IV, and the title of Ref. [49] contains an unusual capitalization/style ('Dampe squib?'); please check against the published record.","section":"Sec. IV and Ref. [49]"}],"recommendation":"major_revision","confidential_remarks":"The paper's core calculation is sound and the authors are appropriately transparent about the illustrative DAMPE normalization and the fitting systematics. My main concern is that the headline conclusion goes beyond the demonstrated parameter space: the environmental scan is absent, and the sharpness threshold is not quantified. These are fixable issues rather than fundamental flaws, so I would be willing to accept after a parameter scan or a suitably restricted claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: nothing conceptually revolutionary here—the qualitative result that stochastic IC broadening erases sharp pulsar edges is already in John & Linden (2023), and the Green-function machinery is standard. What the paper adds is a systematic beta_eff map across ages and injection histories, plus energetics tied to the DAMPE residual. That is a genuinely useful reference for the DAMPE/HERD debate. The paper is honest and well-scoped. It uses a 64-model GALPROP ensemble for the diffuse background, derives burst and continuous Green-function solutions cleanly, and actually implements a Klein-Nishina Monte Carlo for stochastic IC rather than pasting in a fudge factor. The table of beta_eff values is the core deliverable, and the bootstrap uncertainties for the MC rows are appropriate. The text repeatedly flags that the DAMPE normalization is illustrative and that fitting systematics are not included. Credit where due. Soft spots, in order. First, the central diagnostic claim—that a sharp, stable edge would be difficult for pulsars—is supported only for one environmental setup: B=3 uG and the Geminga-path ISRF. The stress-test note is right: a stronger B or weaker ISRF makes synchrotron dominate, cooling becomes more continuous, and the edge sharpens. So Table I does not establish an upper bound on beta_eff for ordinary pulsar propagation. The paper's language ('more difficult', 'strengthen the case') is appropriately cautious, but the Sec. V sentence leans on a bound that is not demonstrated. Second, reproducibility: no code or GALPROP files are shipped, so the Monte Carlo results cannot be independently checked. Third, Table I reports uncertainties only for the stochastic IC rows; the continuous-injection and SNR rows have no errors, and the fixed-p fitting assumption is discussed but not varied. These are minor-to-moderate, not fatal. The central qualitative direction—stochastic IC broadens the edge, continuous injection inherits sharpness from the cutoff—holds. Who it is for: people doing phenomenology of TeV electron data, especially in the DAMPE/HERD context. It deserves a serious referee: it is a transparent, well-scoped calculation with a clear, stated limitation. I would accept it for review, and I would cite it in my own writing on electron spectral features.","headline":"Systematic beta_eff map of pulsar electron features; solid and honest, but the 'sharp edge' diagnostic lacks a demonstrated environmental bound.","tokens_in":14578,"tokens_out":1702,"would_cite":true,"duration_ms":19375,"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":"Nearby pulsars can make TeV electron features only as smooth, age-dependent cutoffs.","keywords":["cosmic-ray electrons","cosmic-ray positrons","pulsar wind nebulae","inverse-Compton cooling","Klein-Nishina scattering","diffusion-loss equation","TeV spectral features","dark matter interpretation"],"falsifier":"Take a high-resolution all-electron spectrum across a TeV feature and fit its post-peak decline with $A E^{-p}\\exp[-(E/E_{\\rm eff})^{\\beta_{\\rm eff}}]$; if $\\beta_{\\rm eff}$ comes out consistently above about 4 and independent of the fitted energy window or assumed source age, then the paper's claim that ordinary pulsar propagation cannot produce a sharp stable edge would be contradicted.","tokens_in":13538,"feed_emoji":"🔭","tokens_out":11088,"duration_ms":102813,"temperature":0.7,"pith_summary":"The paper asks whether the TeV-scale structures seen in the cosmic-ray electron and positron spectrum can be produced by ordinary nearby sources, pulsars and supernova remnants, within a single diffusion-loss propagation framework. It shows that once inverse-Compton cooling is treated stochastically with the Klein–Nishina kernel, the sharp cooling edge predicted by continuous-loss approximations becomes a broad, age-dependent decline with effective sharpness parameter $\\beta_{\\rm eff}$ between roughly 1.3 and 3.6. Continuous-injection pulsars inherit a moderately sharp falloff from an assumed super-exponential intrinsic cutoff, with $\\beta_{\\rm eff}$ near 2, while a supernova-remnant-like burst template is broader. The upshot is that only a very sharp and stable edge-like feature would be difficult to explain with ordinary pulsar propagation and would strengthen the case for dark-matter or other exotic origins.","feed_headline":"Sharp TeV electron edge would point beyond pulsars","feed_subtitle":"New modeling shows realistic pulsar and supernova-remnant signals are smooth; a sharp stable cutoff would favor exotic origins.","key_machinery":"The central object is the Green-function solution of the diffusion-loss equation for a point-like burst or continuous source, combined with a Monte Carlo treatment of inverse-Compton cooling. Instead of a deterministic loss law $b(E)=b_0 E^2$, each injected electron's trajectory samples individual inverse-Compton scatterings with the isotropic Klein–Nishina kernel against the tabulated interstellar radiation field, path-averaged optical, infrared, and CMB components along the line of sight of a nearby mature pulsar, while synchrotron losses are treated continuously with $B = 3\\,\\mu$G. Each trajectory carries a diffusion weight computed from the integrated diffusion coefficient along its cooled path. This machinery is what converts the sharp deterministic cooling boundary into the broadened, age-dependent declines characterized by $\\beta_{\\rm eff}$.","core_discovery":"On its own terms, this paper's central claim is that the apparent sharpness of any TeV-scale feature in the all-electron spectrum is a usable diagnostic: naive deterministic cooling predicts a narrow edge, but realistic stochastic inverse-Compton losses broaden it to effective indices $\\beta_{\\rm eff}\\sim 1.3$--$3.6$ that vary with source age, while continuous injection yields $\\beta_{\\rm eff}\\sim 2$ inherited from the source cutoff and a supernova-remnant-like burst source yields $\\beta_{\\rm eff}\\simeq 2.7$. Consequently, a very sharp and stable edge-like feature would be difficult to accommodate with ordinary pulsar propagation and would strengthen the case for alternative origins. The claim is established by constructing the diffuse background from a calibrated propagation-code ensemble, then adding discrete-source Green-function solutions with radiative losses, with stochastic inverse-Compton scattering sampled through a Klein–Nishina Monte Carlo using a path-averaged interstellar radiation field along the line of sight of a nearby mature pulsar.","pith_inferences":["Running the diagnostic in reverse: once future high-resolution measurements fit a feature's post-peak decline, $\\beta_{\\rm eff}$ becomes a cheap classification statistic, with values near 2–3 favoring local astrophysical sources and values above about 4 with no source-age dependence shifting the balance toward exotic origins.","The same stochastic inverse-Compton Monte Carlo could be applied to dark-matter-induced electron signals, since those signals also suffer radiative losses; this would put astrophysical and exotic interpretations on equal footing when comparing the shapes of their high-energy cutoffs.","The path-averaged interstellar radiation field is the main environmental input, so independent maps of the local radiation and magnetic environment from gamma-ray and radio observations could tighten the allowed $\\beta_{\\rm eff}$ range for pulsar interpretations."],"forward_implications":["If a TeV feature is really a cooling-broadened pulsar contribution, its post-peak sharpness should vary smoothly with source age, giving $\\beta_{\\rm eff}$ from about 1.3 for a young source to 3.6 for an old one rather than a stable universal edge.","The sharpest pulsar-like feature allowed by this treatment has $\\beta_{\\rm eff}\\simeq 2$, and that requires a super-exponential intrinsic cutoff in the injected pair spectrum; the sharpness is inherited, not produced by propagation.","A burst-like mature pulsar at roughly 500 pc needs about $5\\times10^{47}$ erg of injected pairs at 100 kyr and about $5\\times10^{48}$ erg at 300 kyr to reach 10% of the diffuse background at 1 TeV, so older examples are energetically demanding.","A young continuous-injection pulsar at roughly 100 pc can plausibly explain a TeV component for ages around 10 kyr, with a required power of about $1\\times10^{34}$ erg/s, but a 300-year-old source would require a power that only the most energetic pulsars could supply.","A supernova-remnant-like burst template gives the broadest contribution, $\\beta_{\\rm eff}\\simeq 2.7$, so it is the least able to mimic a particle-physics edge."],"supporting_citations":[{"why":"It supplies the Green-function solution for burst-like electron propagation with radiative losses that this paper builds on.","marker":"[7]"},{"why":"It establishes the framework of nearby discrete sources imprinting observable features on the high-energy electron spectrum.","marker":"[8]"},{"why":"It provides the supernova-remnant-like burst template and the argument that nearby remnants generate source-dependent 1–10 TeV electron structures.","marker":"[9]"},{"why":"It supplies the measured all-electron spectrum used to define the TeV break and to illustrate a possible residual local component.","marker":"[18]"},{"why":"It defines the pulsar pair-energy budget used to judge whether the required injected energies are plausible.","marker":"[27]"},{"why":"It supplies the numerical propagation code used to generate the diffuse-background ensemble.","marker":"[31]"},{"why":"It calibrates the propagation-model parameter set used for the diffuse background against gamma-ray observations.","marker":"[32]"},{"why":"It supplies the interstellar radiation field model and the line-of-sight-averaged energy densities used in the Monte Carlo.","marker":"[43]"},{"why":"It establishes that stochastic inverse-Compton cooling broadens pulsar spectral features, the effect this paper quantifies with a dedicated Monte Carlo.","marker":"[44]"}],"fun_headline_variants":["TeV electron edge too sharp for pulsars alone","Smooth pulsar glow vs sharp cosmic-ray edge","If TeV electron edge is real, pulsars can't explain it","Pulsar signals smooth; sharp TeV edge hints at exotic"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The quantitative sharpness values hinge on the assumed local radiation and magnetic environment—a path-averaged interstellar radiation field along one nearby pulsar's line of sight and a uniform 3 microgauss field—together with a fixed diffusion law, so if those environmental parameters differ, the fitted $\\beta_{\\rm eff}$ values shift, although the qualitative broadening is likely robust.","fun_headline_variants_meta":{"raw":{"variants":["TeV electron edge too sharp for pulsars alone","Smooth pulsar glow vs sharp cosmic-ray edge","If TeV electron edge is real, pulsars can't explain it","Pulsar signals smooth; sharp TeV edge hints at exotic"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000484,"raw_usage":{"total_tokens":2395,"prompt_tokens":954,"completion_tokens":1441,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":570,"completion_tokens_details":{"reasoning_tokens":1371}},"tokens_in":570,"tokens_out":1441,"duration_ms":9747,"temperature":1.0,"reasoning_tokens":1371,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T19:08:57.019958+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a high-resolution all-electron spectrum across a TeV feature and fit its post-peak decline with $A E^{-p}\\exp[-(E/E_{\\rm eff})^{\\beta_{\\rm eff}}]$; if $\\beta_{\\rm eff}$ comes out consistently above about 4 and independent of the fitted energy window or assumed source age, then the paper's claim that ordinary pulsar propagation cannot produce a sharp stable edge would be contradicted.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It supplies the Green-function solution for burst-like electron propagation with radiative losses that this paper builds on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It establishes the framework of nearby discrete sources imprinting observable features on the high-energy electron spectrum."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It provides the supernova-remnant-like burst template and the argument that nearby remnants generate source-dependent 1–10 TeV electron structures."},{"cited_title":"Ambrosiet al.(DAMPE Collaboration), Direct detec- tion of a break in the teraelectronvolt cosmic-ray spec- trum of electrons and positrons, Nature552, 63 (2017)","cited_arxiv_id":null,"evidence_quote":"It supplies the measured all-electron spectrum used to define the TeV break and to illustrate a possible residual local component."},{"cited_title":"Profumo, Dissecting cosmic-ray electron-positron data with occam’s razor: the role of known pulsars, Central European Journal of Physics10, 1 (2012)","cited_arxiv_id":null,"evidence_quote":"It supplies the numerical propagation code used to generate the diffuse-background ensemble."},{"cited_title":"Ackermann, M","cited_arxiv_id":null,"evidence_quote":"It calibrates the propagation-model parameter set used for the diffuse background against gamma-ray observations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It supplies the interstellar radiation field model and the line-of-sight-averaged energy densities used in the Monte Carlo."},{"cited_title":"The second case is a young nearby pulsar that continu- ously injects pairs over its finite lifetime","cited_arxiv_id":null,"evidence_quote":"It establishes that stochastic inverse-Compton cooling broadens pulsar spectral features, the effect this paper quantifies with a dedicated Monte Carlo."}],"review_version":1}