{"id":"199a7a31-1a85-4414-8354-735c47481659","arxiv_id":"1908.03216","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"A sample of 27 Galactic pulsar wind nebulae shows a diffusion coefficient about two orders of magnitude below the Galactic average, implying large low-diffusion bubbles around these pulsars.","lead":"This paper uses gamma-ray maps from HESS and HAWC to measure how fast electrons and positrons drift away from 27 pulsars. It finds the particles move much slower than average for the Galaxy, implying each pulsar sits inside a low-diffusion bubble at least 80 parsecs wide.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central inference that low D0 is an ISM diffusion coefficient rests on the unproven claim that the TeV emission is from e± that have escaped the PWN; the paper's own two-zone fit shows a D0-rb degeneracy.","rationale":"The reader's weakest assumption is exactly the escape/ISM assumption, and I find no independent argument in the paper that neutralises it. The paper's one-zone model fit quality is real evidence that the ICS model can describe the data, but it is not discriminating evidence against confinement: a confined nebular profile could also be centrally peaked with similar angular scale, and the paper does not fit such a model. The age-split comparison (young vs old) is relevant but not decisive because confinement depends on SNR density and magnetic field, not only on age. The abstract's '15-80 pc' range also appears inconsistent with Table IV (e.g. HESS J1026-582 at 6±2 pc and HESS J1841-055 at 180±30 pc), but this is secondary. Because the central claim is explicitly conditional and the load-bearing condition is testable, the appropriate verdict remains CONDITIONAL. No change from the reader's verdict.","tokens_in":29306,"tokens_out":3980,"duration_ms":45422,"concrete_test":"Recompute, for each of the 27 sources in Table IV, the e± number density at the radius of the VHE emission using the paper's injection spectrum (Eq. 3 with the Table IV gamma_e and the same spin-down parameters) and the θICS/2 halo sizes, following the density comparison of Giacinti et al. (2019, Ref. [11]). If the majority of sources still have densities larger than the local ISM density, the VHE radiation is likely produced inside the PWN/SNR and the D0 values in Table IV are not ISM diffusion coefficients; if the densities become comparable to or below the ISM density with the corrected sizes, the escape assumption is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is the identification of D0. Equations (5), (9)-(11) and the surface-brightness fits in Sec. VI treat every source as a point-like injector of e± into a homogeneous ISM with one-zone diffusion; D0 is then read off from the profile. This is only the ISM diffusion coefficient if the radiating e± have actually left the PWN/SNR. The paper explicitly flags this as unresolved (Sec. I) and cites Ref. [11], which concludes that for most HGPS PWNe the e± density is above the ISM value, i.e. confinement. The rebuttal in Sec. VI argues that the halo sizes used by [11] are underestimated by about a factor of 2; but rescaling a density does not by itself prove escape, and the one-zone fits never include a confinement component. Moreover, the 'bubble size at least 80 pc' claim is derived by combining the fitted one-zone D0 with Fig. 4; the only direct two-zone fit (HESS J1825-137, Fig. 11) shows a flat chi2 over rb > 60 pc with D0 varying by an order of magnitude. So neither the ISM-escape assumption nor rb is directly constrained. If most electrons are nebular, the sample-wide low D0 is a property of the PWN environment, not evidence for interstellar low-diffusion bubbles.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper investigates whether the low-diffusion halos seen around Geminga and Monogem are a generic property of Galactic pulsar wind nebulae (PWNe). The authors model the inverse-Compton (IC) gamma-ray emission from electrons/positrons injected by PWNe, using a one-zone diffusion model with a two-zone variant containing an inner low-diffusion bubble of radius r_b, and study the 68% containment angle as a function of pulsar age, distance, D0 and r_b. They forecast the number of IC halos detectable by HAWC, HESS and CTA as a function of the pulsar-to-lepton efficiency η, and they construct a sample of 27 HGPS/HAWC sources associated with PWNe. Fitting the HGPS surface-brightness maps and the HAWC data for Geminga and Monogem source by source with D0 and η free, they obtain a sample-mean D0 ≈ 9×10^25 cm2/s at 1 GeV (D(1 TeV) ≈ 8×10^26 cm2/s), about two orders of magnitude below the average Galactic value, and infer low-diffusion bubbles of size at least about 80 pc around the sources. The paper also ranks pulsars by predicted IC flux and identifies the brightest expected halos.","tokens_in":29571,"tokens_out":11491,"duration_ms":114734,"significance":"If correct, the result would generalize the Geminga/Monogem low-diffusion halos to a substantial sample of Galactic PWNe, with implications for cosmic-ray transport in the inner Galaxy and for the pulsar contribution to the positron excess. The paper's strengths are its use of public HGPS and HAWC data, the source-by-source surface-brightness fits, the report of statistical errors, the explicit test of the correlation-radius choice, and the falsifiable predictions (source ranking and detection counts) that can be checked with future HAWC and CTA observations. The central inference, however, is conditional on the assumption that the TeV emission is produced by leptons that have escaped into the interstellar medium; the one-zone fits measure an effective diffusion scale in the emission region, and the two-zone fit shows a D0–r_b degeneracy. The sample-wide claim is therefore plausible but not yet demonstrated at the level claimed in the abstract.","major_comments":[{"comment":"The identification of the fitted D0 with the interstellar diffusion coefficient is the load-bearing step and is not established. The model assumes a point-like injector in a homogeneous ISM, so the profile fits in Sec. VI can only yield D0 in this interpretation if the radiating e± have actually left the PWN/SNR. The paper itself flags this as unresolved in Sec. I and cites Ref. [11], which concludes that for most HGPS PWNe the e± density exceeds the ISM value, i.e., confinement. The rebuttal in Sec. VI, that Ref. [11] underestimates halo sizes by about a factor of two, changes the density estimates but does not by itself demonstrate escape; a confined electron population radiating inside the PWN/SNR environment would also produce an extended profile that the one-zone diffusion fit could absorb into a low effective D0. This concern is concrete: the sample includes objects whose TeV emission is commonly attributed to the PWN itself, e.g., HESS J0835-455 (Vela X) in Tab. III. Please either soften the claim to an effective diffusion coefficient in the emission region or add a model component that distinguishes confined and escaped leptons.","section":"Secs. I and VI, Eqs. (5), (9), (11)"},{"comment":"The claim of a low-diffusion bubble size of at least 80 pc is not directly constrained by the data. The one-zone fits in Tab. IV fix r_b effectively infinite, and the only two-zone fit, HESS J1825-137 in Fig. 11, shows a strong D0–r_b degeneracy: the text reports r_b > 60 pc, with the best-fit D0 changing from about 16×10^26 cm2/s at r_b = 60 pc to about 2×10^26 cm2/s at r_b = 120 pc, and a chi2 that is flat for r_b > 70 pc. The 80 pc value therefore comes from combining the one-zone D0 with the model curves in Fig. 4, not from a fit. Please present the bubble-size result as a model-dependent lower limit under the one-zone assumption and report the degenerate D0–r_b range explicitly.","section":"Sec. VI, Fig. 11"},{"comment":"The surface-brightness fits do not include a background component. The paper states in Sec. V.B that, assuming the background is isotropic, it should act as a mere normalization without changing significantly the angular profile of the TeV surface brightness. This is not correct: a constant isotropic background adds a constant to the azimuthally averaged surface brightness, which flattens the profile at large angles and changes the fitted slope of the diffusion model unless the background is fitted simultaneously. Since the ROIs extend to 0.7–1.1 degrees and the HGPS maps contain interstellar emission and unresolved sources, the fitted D0 values may be biased. Please add a constant background parameter to the fits or demonstrate with off-source regions that the effect is negligible.","section":"Sec. V.B, Sec. VI"}],"minor_comments":[{"comment":"The abstract states that all sources are extended with gamma-ray emission of about 15–80 pc, but the sizes in Tab. IV range from about 6 pc (HESS J1026-582) to about 180 pc (HESS J1841-055); please clarify which statistic is quoted in the abstract and make the numbers consistent.","section":"Abstract, Tab. IV"},{"comment":"The same numerical value 8.2+20.9-5.9 appears both as D0 = 8.2×10^25 cm2/s for annuli of 0.1 degrees and as D(1 TeV) = 8.2×10^26 cm2/s, without stating the assumed spectral index delta used for the energy rescaling; please make the D0-to-D(E) relation explicit.","section":"Sec. VI"},{"comment":"Tab. III uses HESS J1837-069 / 2HWC J1837-065 while Tab. IV and parts of the text use HESS J1837-065; please unify the source naming.","section":"Tabs. III–IV"},{"comment":"The forecast uses a fixed HAWC sensitivity of 10^-14 (TeV cm2 s)^-1 and a fixed HESS sensitivity of 10^-12 (TeV cm2 s)^-1, with no declination or longitude dependence; the predicted detection numbers should be labeled as order-of-magnitude estimates.","section":"Sec. IV.A"},{"comment":"There are several typos and minor wording issues: 'PNW gamma-ray emission' in Sec. III.A, 'with an efficiency a slow as a few %' in Sec. IV.A, 'the constrain of the HAWC field of view' in Sec. IV.B, and 'an other one' before Fig. 4.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"To the editor: the manuscript makes good use of public HGPS and HAWC data and provides timely forecasts, but the central claim is stated more strongly than the fits justify given the known confinement/escape ambiguity. The authors should be asked to either reframe the D0 result as an effective diffusion coefficient in the emission region or add a confinement-aware model component. If the framing is corrected, the paper would be a solid contribution to the journal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things worth knowing. This is the first paper to derive D0 around a sample of 27 PWNe from the public HGPS/HAWC surface-brightness maps, and it finds a sample-wide suppression of about two orders of magnitude relative to the Galactic average. It also produces concrete detection forecasts for HAWC, HESS and CTA. The catch is one the authors themselves flag: the fitted D0 is only the interstellar diffusion coefficient if the TeV emission comes from e± that have escaped the PWN and are propagating in the ISM. If most of the leptons are still nebular, as Ref. [11] argued, the same fits measure the PWN environment, not ISM diffusion.\n\nWhat is genuinely good. The per-source fits are honest work. They use publicly released HGPS flux maps with Rc = 0.1 deg and cross-check with 0.2 deg, validate their surface-brightness extraction against the Gaussian sizes in the HGPS catalog, and report 1-sigma uncertainties. The old/young split and the discussion of selection effects are useful. The detection forecasts are falsifiable and could be checked against future HAWC/CTA data. That is real value.\n\nWhere it is soft. The load-bearing weakness is the escape assumption, and the rebuttal in Sec. VI does not fully close it. They show that using their larger halo sizes lowers the e± densities in [11], but a factor-two halo rescaling does not prove escape. The one-zone model, by construction, cannot separately constrain the bubble radius rb; the two-zone fit for HESS J1825-137 exhibits exactly that degeneracy, with D0 varying by an order of magnitude for rb > 60 pc. So the claim of low-diffusion bubbles with size at least 80 pc is inferred indirectly, not measured. There are smaller issues: no background component is modeled, the spectral index delta is never stated even though it enters the propagation, and the abstract's 15-80 pc spread does not match Table IV's 6-180 pc range.\n\nBottom line. This is a paper worth engaging with, not a result to take as settled. A serious referee should send it back with requests to confront the escape/confinement question head-on, quantify the two-zone degeneracy for more sources if possible, and fix the internal inconsistencies. It should not be desk-rejected.","headline":"First sample-wide D0 measurement around 27 PWNe, with real forecasts and fit work, but the ISM-escape interpretation is explicitly conditional and the two-zone degeneracy remains unresolved.","tokens_in":30133,"tokens_out":2694,"would_cite":true,"duration_ms":30226,"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":"Across a sample of 27 pulsar wind nebulae, gamma-ray halos imply diffusion coefficients two orders of magnitude below the Galactic average, making low-diffusion bubbles of at least 80 pc a general property of Galactic pulsars.","keywords":["pulsar wind nebulae","inverse Compton scattering","cosmic-ray diffusion","TeV gamma-ray halos","low-diffusion bubbles","gamma-ray astronomy","positron excess"],"falsifier":"Measure one bright halo, such as HESS J1825-137, out to radii well beyond its fitted bubble with CTA or HAWC: if the surface brightness breaks to the steeper, more extended profile predicted for Galactic diffusion beyond $r_b$, the two-zone picture is confirmed; if the steep low-diffusion profile continues to several hundred parsecs with no break, the bubble is much larger than 80 pc. If instead deep X-ray or radio mapping shows no escaping particles beyond the PWN boundary, the assumption that the TeV halo is made of interstellar-diffusing leptons collapses.","tokens_in":2060,"feed_emoji":"🌌","tokens_out":2277,"duration_ms":93399,"temperature":0.7,"pith_summary":"The paper asks whether the very slow cosmic-ray diffusion measured around the Geminga and Monogem pulsars is a quirk of those two objects or a general property of Galactic pulsar wind nebulae, and it argues for the second option. Using the extended gamma-ray halos of 27 pulsar wind nebulae from the HESS Galactic plane survey and the second HAWC catalog, the authors model the emission as inverse Compton scattering of electrons and positrons released into the interstellar medium. They find a mean diffusion coefficient of $D_0 = 9.1^{+17.4}_{-6.0}\\times10^{25}\\,\\mathrm{cm^2\\,s^{-1}}$ at 1 GeV, about two orders of magnitude below the average for the Galaxy, with halo sizes of tens of parsecs that require low-diffusion bubbles at least roughly 80 pc across. If right, the result turns the low-diffusion bubble from a local curiosity into a standard feature of pulsar environments, with consequences for the positron excess and for cosmic-ray transport modeling.","feed_headline":"Pulsar halos reveal diffusion two orders below the Galactic average","feed_subtitle":"Analysis of 27 pulsar wind nebulae finds the slow diffusion first seen near Geminga is a Galaxy-wide pattern.","key_machinery":"The load-bearing object is the two-zone diffusion model, in which the diffusion coefficient is $D(E,r)=D_0(E/1\\,\\mathrm{GeV})^\\delta$ inside a bubble of radius $r_b$ and jumps to the larger Galactic value outside. The electrons and positrons are injected continuously with a spin-down-powered spectrum $Q(E,t)\\propto L(t)E^{-\\gamma}e^{-E/E_c}$, propagated with energy losses and a diffusion length $\\lambda$, and converted into gamma rays by inverse Compton scattering off the cosmic microwave background, infrared, and starlight fields. Halo sizes are quantified by the 68% containment angle $\\theta_{68}$ and by the empirical $\\theta_{\\mathrm{ICS}}$ profile used for the catalog fits; the fitted surface-brightness profile is the mechanism that turns an observed angular extension into a diffusion coefficient.","core_discovery":"The central claim is that essentially every pulsar wind nebula in the sample sits inside a region where cosmic rays diffuse far more slowly than they do in the bulk of the Galaxy. The evidence is morphological: the TeV surface-brightness profiles of the 27 sources are well described by a leptonic inverse-Compton halo model with a low diffusion coefficient, whereas a simple Gaussian is often a worse fit. The mean fitted value is $D_0 = 9.1^{+17.4}_{-6.0}\\times10^{25}\\,\\mathrm{cm^2\\,s^{-1}}$ at 1 GeV, equivalent to $D(1\\,\\mathrm{TeV}) = 8.2^{+20.9}_{-5.9}\\times10^{26}\\,\\mathrm{cm^2\\,s^{-1}}$, with no significant difference between pulsars younger or older than 20 kyr and no trend with age. The inferred halo radii, typically about 35 pc and exceeding 100 pc for several sources, are then read as lower limits on the bubble radius $r_b$: for $D_0\\sim10^{26}\\,\\mathrm{cm^2\\,s^{-1}}$, bubbles smaller than about 80 pc would truncate the observed degree-scale halos. On this basis the paper concludes that low-diffusion bubbles of at least 80 pc are a general property of Galactic pulsar wind nebulae.","pith_inferences":["A decisive follow-up would measure one bright halo, such as HESS J1825-137, at GeV and TeV energies simultaneously with Fermi-LAT and CTA; the joint fit could break the $D_0$–$r_b$ degeneracy that the paper's own two-zone check exposes, turning the current lower limit on the bubble radius into a measured value.","If low-diffusion bubbles are ubiquitous, the standard practice of using a single average Galactic diffusion coefficient for cosmic-ray propagation may need to treat pulsar surroundings as a population of 'slow zones', which would alter predictions of the pulsar contribution to the AMS-02 positron excess.","The paper's finding of no age trend in $D_0$ contradicts the theoretical expectation of a strong increase in diffusion after roughly 100 kyr; if confirmed by a larger sample, transport models around old pulsars must be revised.","Because the TeV halo size is tied to the bubble radius, the same low-diffusion regions should also suppress the arrival of lower-energy cosmic rays from those directions, a signature that might be observable through the gamma-ray shadows of background sources seen behind the bubbles."],"forward_implications":["If the central claim holds, the low-diffusion bubble is a standard feature of pulsar wind nebulae, so models of cosmic-ray escape from pulsars must place every source in an $r_b \\gtrsim 80$ pc low-diffusion region rather than assigning only Geminga and Monogem a special environment.","With a few-percent conversion efficiency, HAWC and HESS should already have detected a few tens of inverse-Compton halos, and CTA should detect roughly four times more, about 100–130, making the population testable in the near future.","The ranking by predicted inverse-Compton flux is a strong selection tool: 21 of the 23 brightest predicted pulsars in the HAWC field are already in the 2HWC catalog, and PSR B1951+32 and PSR J1740+1000 are predicted as likely future detections.","The absence of an age dependence in $D_0$ argues against a strong evolution of the diffusion coefficient with pulsar age over the sampled 3–340 kyr range.","Reinterpreting halo sizes with $\\theta_{\\mathrm{ICS}}$ roughly doubles the halo volumes used in earlier confinement estimates; for at least one source, HESS J1825-137, the revised electron density drops below the interstellar value, weakening the case that all these halos are still nebular rather than interstellar."],"supporting_citations":[{"why":"Supplies the Fermi-LAT analysis of the Geminga halo and the electron propagation plus inverse-Compton formalism that this paper extends to a larger sample.","marker":"[3]"},{"why":"Provides the HAWC detection of Geminga and Monogem TeV halos and the benchmark low diffusion coefficient near $10^{26}$ cm$^2$/s that motivates the sample-wide study.","marker":"[1]"},{"why":"The HESS Galactic plane survey catalog from which the source sample and public flux maps used for the surface-brightness fits are taken.","marker":"[12]"},{"why":"The second HAWC catalog that supplies the extended source detections used for ranking predicted halos and comparing predicted bright sources with observed ones.","marker":"[16]"},{"why":"The competing interpretation that most electrons and positrons are still confined inside the PWNe; the paper revisits its halo-size assumption and shows the conclusion weakens with larger halo sizes.","marker":"[11]"},{"why":"Provides the AMS-02-based average Galactic diffusion coefficient to which the fitted low values are compared.","marker":"[9]"}],"fun_headline_variants":["Pulsar nebulae all sit in slow-diffusion bubbles","Low-diffusion bubbles aren't rare, they're the rule","Galactic pulsars are surrounded by slow-diffusion zones","Pulsar halos confirm slow diffusion is Galaxy-wide"],"cache_read_input_tokens":32256,"weakest_assumption_plain":"The load-bearing premise is that the observed TeV emission is inverse-Compton radiation from electrons and positrons that have escaped the pulsar wind nebula and are diffusing in the interstellar medium, so the fitted $D_0$ measures interstellar diffusion rather than confinement inside the nebula.","fun_headline_variants_meta":{"raw":{"variants":["Pulsar nebulae all sit in slow-diffusion bubbles","Low-diffusion bubbles aren't rare, they're the rule","Galactic pulsars are surrounded by slow-diffusion zones","Pulsar halos confirm slow diffusion is Galaxy-wide"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000345,"raw_usage":{"total_tokens":2030,"prompt_tokens":1219,"completion_tokens":811,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":835,"completion_tokens_details":{"reasoning_tokens":739}},"tokens_in":835,"tokens_out":811,"duration_ms":8459,"temperature":1.0,"reasoning_tokens":739,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:21:04.358403+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure one bright halo, such as HESS J1825-137, out to radii well beyond its fitted bubble with CTA or HAWC: if the surface brightness breaks to the steeper, more extended profile predicted for Galactic diffusion beyond $r_b$, the two-zone picture is confirmed; if the steep low-diffusion profile continues to several hundred parsecs with no break, the bubble is much larger than 80 pc. If instead deep X-ray or radio mapping shows no escaping particles beyond the PWN boundary, the assumption that the TeV halo is made of interstellar-diffusing leptons collapses.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the HAWC detection of Geminga and Monogem TeV halos and the benchmark low diffusion coefficient near $10^{26}$ cm$^2$/s that motivates the sample-wide study."}],"review_version":1}