{"id":"1edca090-c0d5-4475-ac66-2d1f0d1898c0","arxiv_id":"2501.02720","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Using the first LHAASO source catalog, the authors find that unresolved sources can account for the outer-Galaxy diffuse gamma-ray excess above the cosmic-ray model, but not the inner-Galaxy excess.","lead":"Researchers estimated how many faint, unresolved gamma-ray sources hide in the LHAASO survey and how much they contribute to the diffuse emission measured along the Galactic plane. Their estimate can explain the outer Galaxy's excess but leaves the inner Galaxy requiring extra sources or new physics.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The outer-Galaxy consistency claim rests on an extrapolated KM2A source-count slope below the survey threshold; a 1σ change in β1 is not propagated into Figure 5 and could shift the unresolved SED by tens of percent.","rationale":"The paper does a transparent, standard calculation: use the bright catalog-resolved population to infer the unresolved contribution, then compare with the diffuse measurement. The method is publicly reproducible in principle, the source catalog is open, and the authors explicitly list catalog incompleteness and field-of-view issues in Section 6. I also credit the reported check that changing integration limits has a negligible effect. However, that convergence check does not address the more important uncertainty: the normalization and slope of dN/dF in the flux interval just below the detection threshold, which is exactly what Eq. 13 integrates. With β1≈1, the unresolved flux is not sharply divergent at the lowest fluxes, but it is assembled over roughly equal logarithmic bins below ~6×10^-15 ph cm^-2 s^-1; a ±0.2 change in β1 therefore changes the integral at the tens-of-percent level. Since the outer-Galaxy excess is only marginally covered by the combined model-plus-unresolved band, an unsupported slope change could flip the conclusion. The inner-Galaxy conclusion is more robust because the deficit is large, consistent with the reader's assessment. The conditional verdict is appropriate: I would not reject the paper, but the outer-Galaxy consistency claim should not be treated as settled until the Table 1 uncertainties are propagated through Eq. 13 and the faint-end slope sensitivity is quantified.","tokens_in":10310,"tokens_out":11780,"duration_ms":124648,"concrete_test":"Propagate the Table 1 fit covariance (β1, log10 A, Fb, and the WCDA equivalents) through Eq. 13 in a 1000-draw Monte Carlo and redraw Figure 5's outer-Galaxy combined band; check whether the Cao et al. (2023) red points remain inside the 1σ band in every energy bin. Then rerun the KM2A integral with β1 forced to 2.0 (a steep faint-end slope) while holding A and the break fixed; if the unresolved flux drops enough that the outer-Galaxy residual exceeds 1σ, the central claim is contingent on an untested extrapolation. As an independent cross-check, compare the predicted sub-threshold source counts with Fermi 4FGL or H.E.S.S. Galactic-plane counts where the flux ranges overlap.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equation 13 integrates dN/dF down to 10^-18 ph cm^-2 s^-1 for KM2A and 10^-16 ph cm^-2 s^-1 for WCDA, using the broken power law of Eq. 5 with Table 1 parameters. The catalog constrains N(>F) only down to the faintest detected source, so the sub-threshold slope is an assumption. KM2A efficiency reaches ~100% above ~6×10^-15 ph cm^-2 s^-1 (Fig. 4), and with β1≈1.0 each flux decade below that contributes roughly equally to Eq. 13; hence the unresolved flux is set by the extrapolated F≲6×10^-15 population. Table 1 gives β1=1.008±0.208 and log10 A=15.302±0.077, but the Figure 5 bands do not show propagation of these errors, and Section 6 itself concedes faint-source incompleteness as a systematic uncertainty. If the true population steepens below the observed range, the outer-Galaxy 'roughly consistent' conclusion could fail; if it flattens, the claim strengthens and the inner-Galaxy gap shrinks. This is the single most load-bearing link between the catalog extrapolation and the paper's main conclusion.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper estimates the contribution of unresolved gamma-ray sources to the LHAASO-KM2A diffuse Galactic emission by deriving the source count distribution and photon-index distribution from the first LHAASO catalog using the Lynden-Bell C− method. It fits a broken power law for dN/dF and a Gaussian for the photon-index distribution, constructs the detection efficiency, and integrates Eq. (13) down to very faint fluxes to obtain unresolved-source SEDs for the inner and outer Galaxy. Comparing the sum of the unresolved contribution and the CR propagation model of Zhang et al. (2023) with LHAASO-KM2A data, the paper concludes that the outer Galaxy (125°<l<235°) is roughly consistent within uncertainties, while the inner Galaxy requires additional components below about 30 TeV.","tokens_in":10585,"tokens_out":5086,"duration_ms":52717,"significance":"If robust, the result provides a concrete, catalog-based explanation of part of the TeV-PeV diffuse excess and sharpens the need for extra components in the inner Galaxy. The analysis is transparent, uses the public LHAASO catalog, and gives explicit formulas for each step; the authors also acknowledge catalog incompleteness and source-extension limitations in Section 6. However, the headline conclusion is currently supported mainly by central-value comparisons: the propagated uncertainties of the source-count fit are not shown in the SED bands, the detection efficiency is built from the same best-fit model used for the extrapolation, and the assumed source spatial distribution directly rescales the result. These issues are fixable but need to be addressed before the quantitative claims can be accepted.","major_comments":[{"comment":"The quoted SED bands for unresolved sources do not include the uncertainties of the fitted source-count parameters. In particular, the KM2A sub-break slope is β1=1.008±0.208, and because the detection efficiency saturates above F≈6×10^-15 ph cm^-2 s^-1 (Fig. 4), the unresolved flux integral is dominated by the extrapolated population below that flux. A 1σ change in β1 can shift the unresolved SED by tens of percent, which is enough to affect the 'roughly consistent' outer-Galaxy conclusion. Please propagate the Table 1 parameter uncertainties, including their covariance, into the Figure 5 bands and state the resulting range for the outer-Galaxy residual.","section":"Section 5, Figure 5, Eq. (13), Table 1"},{"comment":"The detection efficiency λ(F) is constructed as the ratio of the observed source count to the best-fit dN/dF, and the same best-fit broken power law is then used in Eq. (13) to extrapolate below the detection threshold. This makes the unresolved flux estimate dependent on the functional form assumed for the source count distribution in a way that is not captured by the statistical errors in Table 1. I ask for a robustness test with at least one alternative parameterization (e.g., a steeper or flatter sub-threshold slope, or a smoothly curved function) and a statement of how the outer-Galaxy conclusion changes under that alternative.","section":"Section 4.1, Eq. (13)"},{"comment":"The ROI photon fractions f_inner,ROI≈0.31 and f_outer,ROI≈0.22 are computed from an assumed spatial distribution of sources (Eq. 10) with fixed parameters (r⊙=8.5 kpc, zs=0.2 kpc) and no uncertainty. Since the unresolved SED in Eq. (13) is proportional to f_i,ROI, a change in the scale height or radial profile directly rescales the predicted contribution. The paper does not validate Eq. (10) against the longitude distribution of the LHAASO catalog sources. Please quantify the sensitivity of the unresolved flux to the assumed spatial distribution or calibrate it with the observed source positions.","section":"Section 4.2, Eqs. (10)–(13)"},{"comment":"The claim that the outer Galaxy is 'roughly consistent' is based on visual overlap of bands, but no quantitative measure (e.g., chi-square, residual significance, or pull) is given for the sum of the CR model and unresolved sources. Because the model and unresolved components are added linearly and the data uncertainties may be correlated across energy bins, a quantitative comparison would make the central claim testable and would clarify whether the agreement is actually within the propagated uncertainties.","section":"Section 5, Figure 5"}],"minor_comments":[{"comment":"There are small typographical errors: 'Lyden-Bell' should be 'Lynden-Bell' in the text and figure caption, and 'sorces' should be 'sources' in the Figure 3 caption.","section":"Section 3.2, Figure 3 caption"},{"comment":"Please define dN/dF dΓ explicitly as the joint differential source count and state the units of the integration limits; the current notation is ambiguous between a product and a joint distribution.","section":"Eq. (13)"},{"comment":"The detection efficiency is shown without error bars; adding Poisson uncertainties would help the reader assess at which flux the 100% completeness claim is reliable.","section":"Figure 4"},{"comment":"The sentence 'We have verified that varying the integration limits for the flux and photon index has a negligible effect' should be supported with a short quantitative demonstration, since the lower limits of 10^-18 and 10^-16 ph cm^-2 s^-1 are far below the catalog threshold and the verification is directly relevant to the extrapolation.","section":"Section 4.3"},{"comment":"The survival-product step of the Lynden-Bell method is terse; a short explanation of why n_j−1 and m_j−1 appear in the product, and how the product is ordered over sources, would improve readability for readers not familiar with the method.","section":"Equations (3) and (7)"},{"comment":"The lower subplots show fractional contributions without uncertainty bands; these fractions inherit the same unpropagated errors noted above and should be updated if the bands are updated.","section":"Figure 5, lower panels"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the journal's scope and the central idea is worth publishing after revision. The main risk is that the headline conclusion is presented with bands that do not include the dominant systematic uncertainties from the source-count extrapolation and the assumed spatial distribution; this is fixable with the requested tests and propagation. The shared authorship with the model paper (Zhang et al. 2023) is not in itself a concern, but the paper should state more explicitly that the model prediction is one class of CR propagation models and that other propagation assumptions could change the required unresolved contribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, this is the first paper I've seen that takes the LHAASO source catalog and uses a non-parametric Lynden-Bell method to estimate the unresolved-source contribution to the diffuse TeV–PeV gamma-ray emission. Second, the headline conclusion — that unresolved sources can explain the outer-Galaxy excess — is plausible but not settled. The inner-Galaxy result, which still leaves a large gap, is more secure.\n\nThe method is transparent and the equations check out. The authors are candid about the catalog's incompleteness, and they test the sensitivity to integration limits. Using the published catalog makes the analysis externally reproducible in principle, even though no code is released. The comparison with data in Figure 5 is clear.\n\nThe soft spots are in proportion. The biggest one is the extrapolated source-count slope below the detection threshold. The catalog constrains N(>F) only down to the faintest detected source, around 6e-15 ph cm^-2 s^-1 for KM2A, but the unresolved flux integral runs down to 1e-18. With beta1 about 1.0, each decade contributes about equally, so the result is set by a source population nobody has seen. A 1-sigma change in beta1 (0.208) could shift the unresolved SED by tens of percent, and that error is not propagated into the Figure 5 bands. The authors mention the systematic in the discussion but don't quantify it. The detection efficiency is reconstructed from the same best-fit model, which is a mild circularity, though the observed counts anchor it. The spatial distribution assumption (a smooth disk profile) is a reasonable first guess but adds unquantified uncertainty.\n\nThe inner-Galaxy conclusion doesn't depend on these issues: the gap is too large for unresolved sources to fill, given the 1.8 inner/outer ratio they derive. That part should survive scrutiny.\n\nOverall, this is a solid, honest paper, and the central outer-Galaxy claim is defensible as a 'roughly consistent' statement. It deserves a serious referee. I'd press the authors to propagate the beta1 uncertainty into the SED bands and run a sensitivity test where the faint-end slope is allowed to steepen or flatten below the observed flux range. That would make the result much more convincing.","headline":"Credible, transparent estimate of unresolved-source contribution to the LHAASO diffuse excess; outer-Galaxy claim is plausible but rests on an unquantified faint-end extrapolation, while the inner-Galaxy gap is robust.","tokens_in":11187,"tokens_out":2539,"would_cite":true,"duration_ms":24874,"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":"Unresolved sources can account for the outer Galaxy's diffuse gamma-ray excess, but not the inner Galaxy's.","keywords":["diffuse gamma-ray emission","unresolved sources","LHAASO","source count distribution","Lynden-Bell C-minus method","TeV–PeV gamma rays","Galactic plane","cosmic-ray propagation"],"falsifier":"A direct measurement of the source-count distribution at fluxes below the current LHAASO threshold, for example from a deeper survey with the same arrays, would settle the claim: if the cumulative N(>F) deviates from the broken power-law extrapolation with faint-end slope β1≈1.0, then the predicted unresolved flux, and with it the outer-Galaxy consistency, fails. A simpler proximate test is to propagate the quoted uncertainty on β1 (1.008 ± 0.208) into the reported SED bands and see whether the outer-Galaxy match survives.","tokens_in":10046,"feed_emoji":"🌌","tokens_out":5955,"duration_ms":53597,"temperature":0.7,"pith_summary":"LHAASO-KM2A has measured diffuse gamma-ray emission from the Galactic plane between 10 TeV and 1 PeV, and the flux exceeds predictions based on cosmic rays interacting with the interstellar medium. This paper asks how much of that excess comes from many individually undetected gamma-ray sources rather than from new physics or unusual cosmic-ray propagation. Using the first LHAASO source catalog, the authors reconstruct the intrinsic source-count distribution with a non-parametric method and extrapolate it below the detection threshold. They find that in the outer Galaxy the unresolved-source contribution, added to the standard cosmic-ray model, matches the measured diffuse flux within uncertainties; in the inner Galaxy it does not, leaving roughly half of the low-energy flux unexplained. The result matters because it separates a concrete, testable astrophysical population from the need for additional components.","feed_headline":"Faint sources explain outer Galaxy's diffuse gamma-ray excess","feed_subtitle":"A source-count extrapolation from the LHAASO catalog closes the outer-Galaxy gap but leaves the inner Galaxy needing extra components.","key_machinery":"The load-bearing tool is the Lynden-Bell C− method, a non-parametric estimator that recovers the intrinsic joint distribution of integrated flux and photon spectral index from a sample truncated by detection thresholds. The authors apply it to the 65 KM2A and 60 WCDA sources in the first LHAASO catalog with |b|<5°, fitting the cumulative flux distribution to a broken power law with faint-end slope β1≈1.0 and bright-end slope β2≈3.3 for KM2A, and a Gaussian photon-index distribution. They then form a flux-dependent detection efficiency λ(F) by comparing observed source counts with the best-fit distribution, model the spatial distribution of sources along the Galactic plane, and integrate the undetected fraction over flux and photon index to get the unresolved flux per solid angle in each ROI. That integral, Equation 13, is what converts an extrapolated catalog population into a prediction for the diffuse emission.","core_discovery":"The paper's central claim is that the diffuse gamma-ray excess seen by LHAASO-KM2A can be partially but not fully attributed to unresolved sources. For the outer Galaxy (125° < l < 235°, |b| < 5°), the flux from sources too faint to be detected, computed from the LHAASO catalog's source-count distribution extrapolated down to $10^{-18}$ ph $cm^{-2}$ $s^{-1}$ (KM2A) and $10^{-16}$ ph $cm^{-2}$ $s^{-1}$ (WCDA), plus the conventional cosmic-ray propagation model, is roughly consistent with the measured SED. The unresolved contribution falls from about 28% to 7% with energy. For the inner Galaxy (15° < l < 125°, |b| < 5°), the same calculation leaves about 50% of the measured flux unexplained below roughly 30 TeV, so the authors conclude that additional components, such as pulsar halos or photon leakage from pulsar wind nebulae, are required there. They also note the flux ratio of unresolved sources between inner and outer regions is about 1.8 from geometry alone.","pith_inferences":["If the outer-Galaxy match is real, the same extrapolation can be turned around: the diffuse emission itself becomes a measurement of the faint-end source-count slope, and future LHAASO data with a lower threshold will provide a direct cross-check.","The assumed spatial distribution of sources, anchored to supernova remnants and pulsars, is a free input; a different radial profile would change the inner-to-outer flux ratio of 1.8 and could weaken or strengthen the inner-Galaxy shortfall, so the conclusion is partly a statement about that distribution.","The WCDA sample shows a hard-spectrum subpopulation with Γ≲2 below 25 TeV; if those sources form a distinct class, the unresolved flux at 1–25 TeV could be dominated by a population not represented in the KM2A band, changing the predicted SED shape."],"forward_implications":["In the outer Galaxy, the diffuse excess between 10 TeV and 1 PeV can be explained without invoking new cosmic-ray physics; the required ingredient is a population of faint sources that follows the catalog's source-count slope.","The unresolved contribution is energy-dependent, decreasing from roughly 28% to 7% in the outer Galaxy and 17% to 5% in the inner Galaxy, so source confusion matters most at the lowest energies of the KM2A band.","For the inner Galaxy, some additional gamma-ray component below about 30 TeV is needed; candidate explanations include pulsar halos and leakage from pulsar wind nebulae.","The source-count distribution inferred from bright LHAASO sources has predictive power for future deeper surveys: a detection of the predicted faint population would confirm the interpretation, while a break in the slope would require revising it."],"supporting_citations":[{"why":"Supplies the measured diffuse gamma-ray SED and the inner and outer Galaxy regions of interest that the paper compares against.","marker":"Cao et al. 2023"},{"why":"Provides the first LHAASO source catalog, including the 65 KM2A and 60 WCDA sources with integrated flux and photon index used in the analysis.","marker":"Cao et al. 2024"},{"why":"Introduces the non-parametric C− method used to recover the intrinsic flux and photon-index distributions from the truncated catalog.","marker":"Lynden-Bell 1971"},{"why":"Gives the threshold-limit fitting procedure and the error treatment adopted when constructing the detection efficiency.","marker":"Singal et al. 2012"},{"why":"Supplies the conventional cosmic-ray propagation model prediction of diffuse emission that serves as the baseline in the comparison.","marker":"Zhang et al. 2023"},{"why":"Provides the radial and vertical source-distribution parameters used to weight the unresolved flux into the inner and outer Galaxy regions.","marker":"Trotta et al. 2011"}],"fun_headline_variants":["Faint sources close outer gamma excess, inner needs extras","Outer gamma excess matches faint sources; inner needs more","Unresolved sources fill outer gamma deficit, inner puzzle remains","LHAASO source counts resolve outer diffuse gamma-ray gap","Gamma excess: faint sources for outer, mystery for inner"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument stands on the assumption that the source-count distribution fitted to catalog sources above the detection threshold continues unchanged down to fluxes of $10^{-18}$ ph $cm^{-2}$ $s^{-1}$ (KM2A) and $10^{-16}$ ph $cm^{-2}$ $s^{-1}$ (WCDA), so the unresolved flux is set by a faint population that has not actually been observed.","fun_headline_variants_meta":{"raw":{"variants":["Faint sources close outer gamma excess, inner needs extras","Outer gamma excess matches faint sources; inner needs more","Unresolved sources fill outer gamma deficit, inner puzzle remains","LHAASO source counts resolve outer diffuse gamma-ray gap","Gamma excess: faint sources for outer, mystery for inner"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001225,"raw_usage":{"total_tokens":5042,"prompt_tokens":957,"completion_tokens":4085,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":573,"completion_tokens_details":{"reasoning_tokens":4002}},"tokens_in":573,"tokens_out":4085,"duration_ms":28211,"temperature":1.0,"reasoning_tokens":4002,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T22:06:10.365070+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct measurement of the source-count distribution at fluxes below the current LHAASO threshold, for example from a deeper survey with the same arrays, would settle the claim: if the cumulative N(>F) deviates from the broken power-law extrapolation with faint-end slope β1≈1.0, then the predicted unresolved flux, and with it the outer-Galaxy consistency, fails. A simpler proximate test is to propagate the quoted uncertainty on β1 (1.008 ± 0.208) into the reported SED bands and see whether the outer-Galaxy match survives.","supporting_citations":[{"cited_title":"2012, The Astrophysical Journal, 753, 45, doi: 10.1088/0004-637X/753/1/45","cited_arxiv_id":null,"evidence_quote":"Gives the threshold-limit fitting procedure and the error treatment adopted when constructing the detection efficiency."}],"review_version":1}