{"id":"50a797ec-2e29-4ecc-9bd3-54ae599c0677","arxiv_id":"2501.08957","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"A Crab-like synthetic population of young pulsar wind nebulae can contribute roughly 5% of the IceCube Galactic neutrino flux near 100 TeV in an optimistic hadronic-proton model.","lead":"The authors simulate a synthetic population of young pulsar wind nebulae in the Galaxy and estimate the neutrinos they could produce if protons are accelerated inside them. Reading it matters because it tests whether unresolved Galactic sources can explain part of the neutrino flux IceCube sees, beyond catalog-based estimates.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Crab-template assumption for proton fraction and filament gas density is the load-bearing step; the 5% neutrino contribution scales with this product and lacks independent calibration for non-Crab sources.","rationale":"The paper is a careful, honestly-hedged population estimate. The reader's weakest assumption — that all simulated PWNe share the Crab's proton fraction and filament amplification — is indeed the most load-bearing step: the 5% result scales linearly with the calibrated product of proton energy fraction and target gas density, and no independent constraint exists for non-Crab sources. My stress-test agrees with this identification. The alternative candidate, the optional flux cut that removes the three most contributing sources, is also noted by the authors as arbitrary, but it acts on only three outliers and is more a robust-estimation choice than a fundamental model input. In contrast, the Crab-template calibration sets the absolute scale for all ~4900 retained sources. The Vf inconsistency (0.15 for Crab, 0.5 for the population) strengthens the concern that the template is not applied consistently, though its numerical effect is not quantified in the paper. Because the claim is explicitly framed as an optimistic upper limit and the limitations are openly discussed, the conditional verdict remains appropriate. The proposed concrete test — a direct model comparison of leptonic-only versus leptonic+hadronic fits to the Crab SED — would settle whether the 1 PeV anchor is unique, and therefore whether the predicted flux is a genuine upper limit or merely one possible interpretation.","tokens_in":21848,"tokens_out":12422,"duration_ms":131506,"concrete_test":"Perform a model comparison on the Crab SED: fit the multi-wavelength spectrum including the LHAASO PeV points with a leptonic-only model (two electron populations and/or radially variable B, no pp contribution) and with the paper's leptonic+hadronic model; compute Δχ² or Bayesian evidence. If the leptonic-only model fits within 1σ of the hadronic model, then the η_p–fa calibration is not uniquely determined and the predicted neutrino flux should be viewed as unconstrained rather than as an upper limit.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central 5% figure at 100 TeV is obtained by calibrating the hadronic normalization on the Crab: η_p = 0.05 is fixed as the maximum allowed after requiring η_e ≥ 0.93 for the leptonic fit, and fa is then adjusted so that the Crab's pp gamma-ray flux matches the LHAASO 1 PeV measurement (Sections 3.1–3.2). This product (η_p times effective gas density) is then assumed identical for every simulated young PWN (Section 4.1). The neutrino flux is directly proportional to this product, so any error propagates linearly into the population sum. There is no independent measurement of the proton fraction in any other young PWN, and the LHAASO PeV excess itself could instead be leptonic (e.g., a second electron population or a spatially inhomogeneous magnetic field). If the hadronic attribution is wrong, the neutrino flux has no empirical anchor and could be arbitrarily smaller. Even if the Crab's hadronic component is real, the same fa is unlikely to hold across a population with different masses, ages, and filament geometries; the paper's own age-dependent scaling of filament saturation is a crude proxy. Additionally, the population calculation uses Vf = 0.5 while the Crab fit uses Vf = 0.15, so the 'Crab-like' template is not applied consistently, which further weakens the single-source calibration.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs a synthetic population of young, free-expanding pulsar wind nebulae in the Galaxy and computes their all-flavor high-energy neutrino emission using a time-dependent, spherically symmetric transport model. The Crab nebula is used as a template: its multi-wavelength spectrum is fitted with synchrotron and inverse-Compton emission from two electron populations plus a hadronic component, and the LHAASO 1 PeV gamma-ray point is used to fix the product of the proton energy fraction and the filament gas amplification factor. Applying the resulting Crab-like parameters to simulated young PWNe, the authors find that, in their optimistic case, the population contributes about 5% of the IceCube Galactic-plane flux around 100 TeV and that at PeV energies the contribution depends strongly on the proton injection index and on a few nearby bright sources.","tokens_in":22248,"tokens_out":5585,"duration_ms":65922,"significance":"If the estimate is taken as an upper bound, the paper provides a useful new estimate of an unresolved, population-level contribution to the Galactic high-energy neutrino flux that is not accessible to catalog stacking. The main strength is the forward-modeling approach: rather than relying on gamma-ray selected catalogs, the authors simulate the PWN population and propagate Crab-calibrated microphysics through a time-dependent transport equation. The paper also explicitly identifies several limitations, including the simplified filament treatment and the uncertain hadronic interpretation of the Crab's PeV excess. The central number, however, is not a measurement-based prediction: it is an optimistic upper limit under the assumption that the Crab's PeV gamma-ray excess is hadronic and that every young PWN shares the Crab's proton fraction and gas amplification. The manuscript would be strengthened by stating this status consistently and by quantifying the linear sensitivity of the neutrino flux to the calibrated parameters.","major_comments":[{"comment":"The synthetic population uses Vf = 0.5, while the Crab template fit that is supposed to calibrate the population uses Vf = 0.15. Since Vf enters Eq. (8) and therefore controls advective transport in Eq. (11), it directly sets the proton residence time and the pp neutrino yield. A source with Vf = 0.5 is not 'Crab-like' in this important respect, so the transfer of the calibrated ηp and fa from the Crab to the synthetic population is not internally consistent. Please either adopt the Crab value Vf = 0.15 for the population or show a sensitivity scan demonstrating that the 100 TeV and 1 PeV conclusions are insensitive to Vf.","section":"§4.1 and Table 3 vs. §3 and Table 1"},{"comment":"The headline '5% of the measured flux by IceCube' is not supported by the quantitative statement in §4.2, which says the post-flux-cut upper band is '∼ 5% relative to the lower limit of the KRAγ model'. The lower edge of a diffuse-model 1σ band is not the same as the measured IceCube flux, and the conclusion's wording ('about 5% of the best-fit results from IceCube') is again different. Moreover, this 5% number is the upper bound of the band after applying the optional flux cut F > 10^-5 erg cm^-2 s^-1 that removes the three dominant nearby sources. The abstract and conclusion should be rephrased to say: under the optimistic hadronic assumption and after excluding a few ultra-bright sources, the upper bound of the predicted young-PWN flux is about 5% of the lower edge of the KRAγ diffuse template.","section":"§4.2, Abstract, and §6"},{"comment":"The neutrino normalization is calibrated to the Crab under a specific hadronic interpretation: ηp = 0.05 is obtained as the maximum allowed after requiring ηe ≥ 0.93, and fa is then tuned so that the pp gamma-ray flux matches the LHAASO 1 PeV point. Because the neutrino and gamma-ray fluxes are produced by the same pp mechanism, this is a circular normalization in the sense that the neutrino flux inherits the assumed hadronic gamma-ray flux, and the result scales linearly with the product ηp × fa. If the PeV excess is leptonic, or if ηp is lower than its maximum in other PWNe, the neutrino flux decreases proportionally. The paper should explicitly label this as an upper-limit scenario and provide at least a one-parameter scaling (e.g., 'the contribution scales as (ηp/0.05) × (fa/fa,Crab)') or a short sensitivity table over ηp and fa.","section":"§3.1–3.2 and Eq. (27)"},{"comment":"The quoted 5% result is the upper envelope of 100 Monte Carlo realizations after the age cut and the optional flux cut, not a central value with propagated uncertainties. The text reports only min/max bands over realizations with fixed physics parameters (ηp, fa, Vf, Crab fitting parameters); the wide distributions in Table 3 for P0, Bs, Esn, and Mej are sampled, but the calibrated parameters themselves are not varied. Since the conclusion depends on three excluded bright sources and on the adopted Vf value, a proper estimate of the spread around the 5% number, including variation of the Crab-calibrated parameters, is needed before this can be presented as a quantitative contribution rather than an illustrative optimistic curve.","section":"§4.2 and Figures 4"}],"minor_comments":[{"comment":"The text contains the typo 'spacial variance' in the description of the diffusion coefficient; it should be 'spatial variance'.","section":"§2.2"},{"comment":"The word 'manetohydrodynamic' appears in the discussion of filament simulations; it should be 'magnetohydrodynamic'.","section":"§3.2"},{"comment":"The sentence describing the number of removed sources is easy to misread: 'Among the sample 71 sources are removed according to the 50-year criterion' should specify that this is the total across 100 realizations, i.e., about 0.7 sources per realization, before the optional flux criterion is applied.","section":"§4.2"},{"comment":"The definition of fa as a 'general effect of accumulation' differs from Atoyan & Aharonian (1996), but no explicit prescription is given for how fa multiplies n_m beyond Eq. (27). A one-sentence clarification that the effective density in the filament zone is neff = fa × nm would remove ambiguity.","section":"Eq. (27)"},{"comment":"The caption and text use 'Log10-normal' for Esn and Bs and 'normal' for P0; please state explicitly whether the parameters are drawn in log10 space and cite the source choices for the standard deviations, since the σ = 35 ms for the initial spin period is stated to be 'somehow arbitrary'.","section":"Table 3"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the forward-modeling approach is genuinely useful. The main issue is that the central 5% number is presented in the abstract and conclusion in a stronger form than the body of the paper supports: it is an upper bound under a hadronic assumption, defined relative to a model lower limit, and it depends on an optional flux cut and an inconsistent Vf choice. These are fixable with rephrasing, a sensitivity scan, and a clearer statement of the upper-limit nature, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The key thing to know is that this paper is a population-synthesis estimate, not a measurement. Its genuinely new piece is that it builds a synthetic population of young, free-expanding PWNe and computes their collective neutrino flux, instead of stacking observed gamma-ray catalog sources. That lets it include unresolved sources, which the catalog approaches miss. The machinery is standard – the Vorster/Moraal transport model with the Crab as the template – and the paper is careful about describing the input distributions and the selection criteria.\n\nWhat it does well: the Crab fit is thorough, with radio and wind electrons, inhomogeneous seed photon fields, and a physically motivated two-region gas density. The authors are transparent about the ad hoc elements: the fa amplification factor, the flux cut that excludes a handful of dominant synthetic sources, and the age cut. They present the result as an optimistic case, and they note that at 1 PeV the contribution depends strongly on the injection index and on the diffuse template.\n\nThe soft spots are real, though. The neutrino flux is proportional to the product of the proton energy fraction eta_p and the effective gas density (via fa), and both are calibrated on the Crab alone, with eta_p taken at its maximum allowed value and fa tuned to match the LHAASO 1 PeV flux. If that gamma-ray excess turns out to be leptonic rather than hadronic, the neutrino normalization loses its anchor and could drop by an order of magnitude or more. The paper acknowledges this, but it remains the load-bearing assumption. There is also an internal inconsistency: the Crab template uses Vf = 0.15, while the population uses Vf = 0.5 as the 'mean value'. That may be a reasonable choice for a population, but it is not justified from the Crab calibration, and it changes the advection/diffusion balance. The 5% figure also depends heavily on removing the 3 brightest sources and the 71 youngest ones; without those cuts, the flux overshoots the IceCube measurement. That doesn't invalidate the exercise, but it means the 5% should be read as 'upper end of a plausible range', not a central prediction.\n\nI was not able to verify the numerics since no code or data are shipped, but the equations are standard and the parameter choices are documented.\n\nWho is this for? Anyone working on Galactic neutrino sources or PWN physics. It is a useful reference for population-level estimates, and the discussion of the Crab fit and the filament treatment is worth reading. It deserves a serious referee: not because the prediction is firm, but because the methodology is clean enough to be tested and improved, and the uncertainty analysis (though partial) is honest.\n\nI would suggest accepting it for peer review with the expectation of major revisions, mainly to address the Vf inconsistency and to add a sensitivity study varying eta_p and fa away from the Crab-calibrated values.","headline":"A synthetic-population estimate of neutrino emission from young PWNe; the 5% IceCube contribution is an optimized upper bound, not a robust prediction, but the paper is honest about what it depends on.","tokens_in":22731,"tokens_out":3765,"would_cite":true,"duration_ms":34033,"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":"This paper claims that unresolved young pulsar wind nebulae, simulated as Crab-like systems, can contribute about 5 percent of the Galactic high-energy neutrino flux detected by IceCube at 100 TeV.","keywords":["pulsar wind nebulae","high-energy neutrinos","IceCube","Crab nebula","hadronic acceleration","Galactic diffuse emission","PeVatrons","population synthesis"],"falsifier":"A measurement showing that the Crab's 1 PeV gamma-ray excess is fully accounted for by inverse Compton emission from electrons, leaving no hadronic component, would remove the calibration of $\\eta_p$ and $f_a$ and collapse the template, falsifying the predicted 5% contribution. Alternatively, IceCube or KM3NeT measuring a Galactic-plane neutrino flux at 100 TeV that is statistically consistent with zero contribution from the three brightest young PWNe in the simulated sample would rule out the optimistic population prediction.","tokens_in":21623,"feed_emoji":"🌌","tokens_out":6605,"duration_ms":62356,"temperature":0.7,"pith_summary":"This paper asks whether young pulsar wind nebulae, the magnetized bubbles blown by young neutron stars while they are still in the free-expansion phase, contribute measurably to the high-energy neutrino sky detected by IceCube. Instead of stacking resolved gamma-ray sources, the authors simulate a synthetic population of Galactic young PWNe, taking the Crab nebula as the standard template for how much pulsar spin-down energy goes into protons and how much target gas is available for proton-proton collisions. In the optimistic case, the predicted all-flavor neutrino flux from this population is about 5% of IceCube's measured flux around 100 TeV. At about 1 PeV, the contribution depends strongly on the proton injection index; with a hard index of 1.5 it can rival the KRA-gamma diffuse template, while with the canonical index 2.0 it falls quickly. The result matters because unresolved discrete sources may be a hidden component of the Galactic neutrino signal, and because young PWNe would be sites of hadronic PeVatron acceleration.","feed_headline":"Young pulsar nebulae could yield 5% of galactic neutrino flux","feed_subtitle":"Simulating all young nebulae as Crab-like systems places unresolved sources in IceCube's 100 TeV signal.","key_machinery":"The carrying machinery is a spherically symmetric, time-dependent model of a free-expanding pulsar wind nebula: pulsar spin-down drives radial advection, a power-law magnetic field $B(r)=B_0(r/R_{\\rm ts})^{-0.5}$, and a transport equation for protons and electrons with diffusion, adiabatic and synchrotron losses, injection at the termination shock, and free escape at the outer boundary. The Crab nebula serves as the template fixing the model parameters, most importantly the proton energy fraction $\\eta_p=0.05$ and the filament gas amplification factor $f_a$ (about 15 for proton index $\\alpha_p=1.5$, about 60 for $\\alpha_p=2.0$), tuned so that the hadronic component reproduces the LHAASO 1 PeV gamma-ray flux. This template is then applied to a synthetic Galactic population of roughly 50 young PWNe per realization, laid down with spiral-arm structure and realistic pulsar birth properties, whose pp neutrino spectra are summed for comparison with the IceCube diffuse templates.","core_discovery":"The central claim is that unresolved young PWNe, simulated as Crab-like systems, produce a collective high-energy neutrino flux that is small but not negligible, reaching roughly 5% of the IceCube-measured Galactic flux at 100 TeV in the optimistic case. To get there, the paper time-evolves each PWN's radius, termination shock, magnetic field, and particle distribution through the free-expansion phase, injects protons at the termination shock, and computes pp neutrino emission from a two-zone gas density whose filament amplification factor is calibrated so that the Crab's hadronic gamma-ray flux matches the LHAASO 1 PeV point. For the synthetic population, only sources younger than their reverse-shock collision time and older than 50 years are kept, and a flux cut removes three extremely bright nearby objects. The authors conclude that while the bulk of the Galactic neutrino flux remains diffuse, a few luminous young PWNe can shape the observable sky at the highest energies, especially if the proton spectrum is hard.","pith_inferences":["If the Crab-like template is representative, future detectors such as KM3NeT, IceCube-Gen2, or Baikal-GVD may resolve individual young PWNe as neutrino point sources, and the same data would directly measure the proton fraction that this paper only calibrates from the Crab.","The model's gas-density prescription is testable at other wavelengths: mapping filament mass and clumpiness in young PWNe of different ages, for example with ALMA or JWST, would confirm or reject the age-dependent amplification factor assumed here.","The paper deliberately excludes reverberation-phase PWNe, but those nebulae can be compressed to higher densities; if protons survive the reverse shock, a follow-up population including that phase could raise the total Galactic PWN neutrino flux above the 5% estimate.","Because the proton fraction is the single most sensitive parameter, a future multi-wavelength fit of any young PWN other than the Crab, such as Vela-X, that includes a hadronic PeV component would provide an independent check on whether the 5% number is a floor or a ceiling."],"forward_implications":["If the optimistic prediction holds, unresolved young PWNe account for roughly one twentieth of IceCube's Galactic neutrino flux at 100 TeV, so diffuse-only templates slightly overestimate the truly diffuse component.","At about 1 PeV with a hard proton index of 1.5, the young-PWN population contributes about 45% of the KRA$^5_\\gamma$ template, meaning discrete hadronic sources can dominate the highest-energy Galactic neutrino sky even if they are minor at 100 TeV.","The three brightest simulated sources alone set the upper bound of the predicted flux; removing them cuts the upper bound by nearly an order of magnitude, so the Milky Way's PeV neutrino sky is likely clumpy rather than smooth.","For the canonical soft proton spectrum with index 2.0, the neutrino flux declines steeply toward PeV energies, making the injection spectral index a decisive observable for distinguishing hadronic acceleration models.","The population approach sidesteps the TeV sensitivity bias of catalog-based stacking, so the 5% estimate includes sources too faint to appear in current gamma-ray catalogs."],"supporting_citations":[{"why":"Supplies the spherically symmetric transport equation, velocity and magnetic-field profiles, and free-escape boundary condition used for the particle distributions.","marker":"(Vorster & Moraal 2013)"},{"why":"Provides the analytic approximations for the PWN outer boundary and reverse-shock radius that define the free-expansion phase for both the Crab and the synthetic population.","marker":"(Bandiera et al. 2023a)"},{"why":"Supplies the Crab's radio-electron distribution, dust and seed photon fields, and the multi-wavelength data used to calibrate the leptonic part of the template.","marker":"(Dirson & Horns 2023)"},{"why":"Provides the LHAASO 1 PeV gamma-ray flux from the Crab that fixes the hadronic normalization through the filament amplification factor.","marker":"(Cao et al. 2021)"},{"why":"Gives the IceCube best-fit Galactic diffuse neutrino templates (pi0, KRA5, KRA50) that the predicted young-PWN flux is compared against.","marker":"(Abbasi et al. 2023)"},{"why":"Supports the filament saturation level and Rayleigh-Taylor structure used to construct the two-zone gas density in the nebula.","marker":"(Porth et al. 2014b)"},{"why":"Supplies the total nebular mass and the fraction of ejecta in the filamentary outer region, which set the mean gas density in the hadronic calculation.","marker":"(Owen & Barlow 2015)"},{"why":"Provides the updated Galactic pulsar surface density and total pulsar count used to generate the synthetic young-pulsar population.","marker":"(Xie et al. 2024)"},{"why":"Supplies the aafragpy package for proton-proton interaction cross sections that converts proton distributions into gamma-ray and neutrino emissivities.","marker":"(Koldobskiy et al. 2021)"}],"fun_headline_variants":["Young pulsar nebulae: 5% of IceCube neutrinos at 100 TeV","Unresolved Crab-like nebulae might add 5% to IceCube signal","Crab-like nebulae could be 5% of the 100 TeV neutrino flux","Simulated pulsar nebulae add 5% to IceCube's high-energy signal"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central assumption is that every young pulsar wind nebula behaves like the Crab, converting the same 5% of spin-down power into protons and having the same filament-amplified target gas density, calibrated to the Crab's 1 PeV gamma-ray excess; if that excess is leptonic, or if other nebulae accelerate protons more weakly, the predicted neutrino flux drops sharply.","fun_headline_variants_meta":{"raw":{"variants":["Young pulsar nebulae: 5% of IceCube neutrinos at 100 TeV","Unresolved Crab-like nebulae might add 5% to IceCube signal","Crab-like nebulae could be 5% of the 100 TeV neutrino flux","Simulated pulsar nebulae add 5% to IceCube's high-energy signal"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001417,"raw_usage":{"total_tokens":5754,"prompt_tokens":1008,"completion_tokens":4746,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":624,"completion_tokens_details":{"reasoning_tokens":4653}},"tokens_in":624,"tokens_out":4746,"duration_ms":33207,"temperature":1.0,"reasoning_tokens":4653,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:13:30.033460+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A measurement showing that the Crab's 1 PeV gamma-ray excess is fully accounted for by inverse Compton emission from electrons, leaving no hadronic component, would remove the calibration of $\\eta_p$ and $f_a$ and collapse the template, falsifying the predicted 5% contribution. Alternatively, IceCube or KM3NeT measuring a Galactic-plane neutrino flux at 100 TeV that is statistically consistent with zero contribution from the three brightest young PWNe in the simulated sample would rule out the optimistic population prediction.","supporting_citations":[{"cited_title":"J., & Barlow, M","cited_arxiv_id":null,"evidence_quote":"Supplies the total nebular mass and the fraction of ejecta in the filamentary outer region, which set the mean gas density in the hadronic calculation."}],"review_version":1}