Pith. Sign in

REVIEW

Diffuse Galactic continuum gamma rays. A model compatible with EGRET data and cosmic-ray measurements

Astrophysics · Submitted 9 Jun 2004 · 79 references

Pith's verdict: A solid GALPROP modeling paper whose 'optimized' fit is honestly labeled as a fit, though the proton normalization is less independent of the gamma rays than the abstract implies.

T0 · AI review

How this review was made: An AI referee (selfhosted/deepseek-v4.1-flash) read the paper on October 29, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 4 major objections · 8 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDastro-ph/0406254 (version 1)
Pith IDQKYY7R4R, Pith's permanent record of this paper and its review
CitationA. W. Strong (MPE, Garching), I. V. Moskalenko (NASA/GSFC), O. Reimer (Bochum). "Diffuse Galactic continuum gamma rays. A model compatible with EGRET data and cosmic-ray measurements." arXiv:astro-ph/0406254 (2004). https://pith.science/paper/QKYY7R4R
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

The Galaxy's diffuse gamma-ray glow can be reproduced from 30 MeV to 100 GeV by a cosmic-ray model whose electron and proton spectra sit only a modest factor above the locally measured ones.

The reading

For three decades the diffuse Galactic gamma-ray spectrum measured by EGRET has shown more emission between about 1 and 10 GeV than models built from the locally measured cosmic-ray spectra predict. This paper sets out to decide which of the proposed repairs is right, by testing each model against the whole sky rather than a single patch: the inner Galaxy, the outer plane, intermediate latitudes, and the poles, and by pushing the comparison out to 100 GeV. The authors show that the standard 'conventional' model fails everywhere, that adding hard sources in the inner Galaxy cannot rescue it because the excess also appears far from the plane, and that simply hardening the electron injection spectrum overshoots the data above 10 GeV while straining the local electron measurement. They then build a model in which the proton spectrum is pinned by secondary antiproton data and the electron spectrum is tuned so that inverse-Compton and bremsstrahlung emission fit the gamma rays themselves, and this single model reproduces the spectra, the longitude profile, and the latitude profile from the plane to the poles. If it holds, the diffuse glow stops being a discrepancy and becomes a measurement: it tells us the large-scale average cosmic-ray spectra, which differ from the local ones by factors of order two, and it makes secondary positrons and electrons a required ingredient of the Milky Way's gamma-ray output.

What carries the argument

The carrying mechanism is the GALPROP propagation and gamma-ray production code run in its two-dimensional (radially symmetric) mode, with gas traced by real H I and CO survey data in Galactocentric rings so that the pion-decay and bremsstrahlung skymaps retain Galactic structure while the cosmic-ray spectra are kpc-scale averages. The decisive move is the fitting order: the nucleon spectrum is not free, it is anchored by the secondary antiproton (and positron) fluxes that sample the same proton population over kiloparsecs, and only then is the electron spectrum adjusted until the sum of inverse Compton, bremsstrahlung, and pion decay matches the gamma-ray spectrum in all seven sky regions. Because the inverse-Compton term dominates at high latitudes and below 100 MeV, it is the piece that carries the fit wherever the gas-related emission cannot.

Extended reading notes

Core claim

The paper's claim is that one self-consistent propagation model can account for the diffuse Galactic continuum from 30 MeV to 100 GeV in every sky region tested. The construction has two steps. First, the Galactic-average proton and helium spectrum is fixed by the locally measured high-energy shape but renormalized upward by about 1.8 and flattened below 10 GeV, a choice constrained by secondary antiproton and positron data rather than by the gamma rays alone. Second, the electron injection spectrum is given a modest steepening (index 1.5 below a break at 20 GV, 2.42 above it) and normalized about a factor of four above the local electron flux, which supplies enough inverse-Compton curvature to match the observed GeV hump without the drastic hardening invoked in earlier work. The result is a full-sky $��^2$ of 462 over 796 data bins, against 1528 for the conventional model and 974 for the hard-electron model, with agreement extending to the 10–20, 20–50, and 50–120 GeV EGRET bands in the inner Galaxy where the hard-electron model is far too bright. The authors also report that secondary positrons and electrons, made in the same proton-gas collisions that produce the pions, contribute roughly half the leptons near 1 GeV, about 20% of the inverse-Compton emission at 1–10 MeV, and about 60% of the bremsstrahlung below 200 MeV, which is what brings the low-energy end of the spectrum into agreement.

Load-bearing premise

The inverse-Compton prediction depends on an interstellar radiation field the authors themselves say could be wrong by a factor of two, and since the electron spectrum is tuned to make that prediction match the gamma rays, an error in the radiation field and the fitted electron normalization cannot be separated.

Editorial extensions

If this is right

  • The diffuse Galactic gamma-ray spectrum becomes a usable measurement of the Galaxy-wide average electron and proton spectra, which come out a factor of a few above the locally measured values rather than wildly different in shape.
  • Secondary positrons and electrons must be included in any diffuse-emission calculation; they supply a large share of the bremsstrahlung below 200 MeV and a fifth of the inverse Compton at MeV energies.
  • The residual extragalactic gamma-ray background is reduced and its spectrum reshaped, because a larger inverse-Compton halo now accounts for part of what was previously attributed to unresolved extragalactic sources.
  • The remaining shortfall below about 30 MeV against COMPTEL points to unresolved compact sources rather than to a cosmic-ray origin, a prediction the authors connect to hard X-ray source counts.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the interstellar radiation field used to compute inverse Compton is really a factor of two higher than assumed, as the authors concede is possible, the required electron normalization drops toward the local value, which would make the model's central tension with local measurements disappear rather than deepen.
  • The flattening imposed on the proton spectrum below 10 GeV is doing double duty: it could equally be absorbing an error in low-energy pion production, so an independent measurement of that cross section would either confirm the spectral break or remove it.
  • A model in which the Solar neighbourhood sits in an interarm region while the gamma-ray-emitting average samples spiral-arm cosmic-ray densities makes a checkable prediction: the offset between local and large-scale spectra should correlate with arm/interarm position, which future spatially resolved gamma-ray data could test.
  • Because the fit leans on inverse Compton at high latitudes, a future direct measurement of the local interstellar electron spectrum above 10 GeV, free of solar modulation, would settle the normalization question without needing the gamma rays at all.
Original abstract

We present a study of the compatibility of some current models of the diffuse Galactic continuum gamma rays with EGRET data. A set of regions sampling the whole sky is chosen to provide a comprehensive range of tests. The range of EGRET data used is extended to 100 GeV. The models are computed with our GALPROP cosmic-ray propagation and gamma-ray production code. We confirm that the "conventional model" based on the locally observed electron and nucleon spectra is inadequate, for all sky regions. A conventional model plus hard sources in the inner Galaxy is also inadequate, since this cannot explain the GeV excess away from the Galactic plane. Models with a hard electron injection spectrum are inconsistent with the local spectrum even considering the expected fluctuations; they are also inconsistent with the EGRET data above 10 GeV. We present a new model which fits the spectrum in all sky regions adequately. Secondary antiproton data were used to fix the Galactic average proton spectrum, while the electron spectrum is adjusted using the spectrum of diffuse emission itself. The derived electron and proton spectra are compatible with those measured locally considering fluctuations due to energy losses, propagation, or possibly details of Galactic structure. This model requires a much less dramatic variation in the electron spectrum than models with a hard electron injection spectrum, and moreover it fits the gamma-ray spectrum better and to the highest EGRET energies. It gives a good representation of the latitude distribution of the gamma-ray emission from the plane to the poles, and of the longitude distribution. We show that secondary positrons and electrons make an essential contribution to Galactic diffuse gamma-ray emission.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-09-23: 1 request, 9,391 input tokens.

What the abstract itself states, as probabilities: claims first 0.04 · states limitation 0.34 · human participants 0.01 · controlled comparison 0.33 · randomized assignment 0.01 · states quantitative result 0.58 · compares with prior results 0.95 · mentions code or data availability 0.04.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

We present a study of the compatibility of some current models of the diffuse Galactic continuum gamma rays with EGRET data. A set of regions sampling the whole sky is chosen to provide a comprehensive range of tests. The range of EGRET data used is extended to 100 GeV. The models are computed with our GALPROP cosmic-ray propagation and gamma-ray production code. We confirm that the "conventional model" based on the locally observed electron and nucleon spectra is inadequate, for all sky regions. A conventional model plus hard sources in the inner Galaxy is also inadequate, since this cannot explain the GeV excess away from the Galactic plane. Models with a hard electron injection spectrum are inconsistent with the local spectrum even considering the expected fluctuations; they are also inconsistent with the EGRET data above 10 GeV. We present a new model which fits the spectrum in all sky regions adequately. Secondary antiproton data were used to fix the Galactic average proton spectrum, while the electron spectrum is adjusted using the spectrum of diffuse emission itself. The derived electron and proton spectra are compatible with those measured locally considering fluctuations due to energy losses, propagation, or possibly details of Galactic structure. This model requires a much less dramatic variation in the electron spectrum than models with a hard electron injection spectrum, and moreover it fits the gamma-ray spectrum better and to the highest EGRET energies. It gives a good representation of the latitude distribution of the gamma-ray emission from the plane to the poles, and of the longitude distribution. We show that secondary positrons and electrons make an essential contribution to Galactic diffuse gamma-ray emission.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 We present a study of the compatibility of some current models of the diffuse Galactic continuum gamma rays with EGRET data. contribution 0.99 0.85 0.04 unclear 0.52
s1 A set of regions sampling the whole sky is chosen to provide a comprehensive range of tests. method 0.98 0.90 0.02 asserted 0.94
s2 The range of EGRET data used is extended to 100 GeV. method 0.93 0.81 0.02 asserted 0.97
s3 The models are computed with our GALPROP cosmic-ray propagation and gamma-ray production code. method 1.00 0.78 0.02 asserted 0.99
s4 We confirm that the "conventional model" based on the locally observed electron and nucleon spectra is inadequate, for all sky regions. result 1.00 0.87 0.13 asserted 1.00
s5 A conventional model plus hard sources in the inner Galaxy is also inadequate, since this cannot explain the GeV excess away from the Galactic plane. result 1.00 0.80 0.06 asserted 1.00
s6 Models with a hard electron injection spectrum are inconsistent with the local spectrum even considering the expected fluctuations; they are also inconsistent with the EGRET data above 10 GeV. result 1.00 0.79 0.07 asserted 1.00
s7 We present a new model which fits the spectrum in all sky regions adequately. contribution 0.71 0.95 0.67 asserted 1.00
s8 Secondary antiproton data were used to fix the Galactic average proton spectrum, while the electron spectrum is adjusted using the spectrum of diffuse emission itself. method 0.99 0.90 0.07 asserted 1.00
s9 The derived electron and proton spectra are compatible with those measured locally considering fluctuations due to energy losses, propagation, or possibly details of Galactic structure. result 0.87 0.90 0.12 asserted 0.61
s10 This model requires a much less dramatic variation in the electron spectrum than models with a hard electron injection spectrum, and moreover it fits the gamma-ray spectrum better and to the highest EGRET energies. result 0.76 0.89 0.33 asserted 1.00
s11 It gives a good representation of the latitude distribution of the gamma-ray emission from the plane to the poles, and of the longitude distribution. result 1.00 0.90 0.08 asserted 1.00
s12 We show that secondary positrons and electrons make an essential contribution to Galactic diffuse gamma-ray emission. result 0.96 0.90 0.16 asserted 1.00

What kind of paper the classifier reports this is: simulation study 0.59 · evidence from analysis of existing data 0.75.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. — B/C ratio as calculated in reacceleration model. Lower curve – LIS, upper – modulated (Φ = 450 MV). Data below 200 MeV/nucleon: ACE (Davis et al. 2000), Ulysses (DuVernois, Simpson, & Thayer 1996), Voy￾ager (Lukasiak, McDonald, & Webber 1999); high energy data: HEAO-3 (Engelmann et al. 1990), for other references see Stephens & Streitmatter (1998). The parameters of the models are summarized in Ta￾ble 1. The models differ only in the injection spec￾tra of protons (and He) and electrons, while the injec￾tion spectra of heavier nuclei are assumed to have the same power-law in…
Figure 2. — Proton spectra as calculated in conventional (solid lines) and optimized (dots) models compared with the data (up￾per curve – LIS, lower – modulated to 650 MV). Thin dotted line shows the LIS spectrum best fitted to the data above 20 GeV (Moskalenko et al. 2002). Data: AMS (Alcaraz et al. 2000b), BESS 98 (Sanuki et al. 2000), CAPRICE 94 (Boezio et al. 1999), IMAX 92 (Menn et al. 2000), LEAP 87 (Seo et al. 1991). EGRET telescope above 10 GeV with an additional sys￾tematic error of ±5%. However, the actual number of photons >10 GeV is small: 1091, 362 and 53 events re￾spectively,…
Figure 4. — γ-ray spectrum of conventional model (44-500180) for the sky regions described in
Figure 5. — Antiproton flux as calculated in conventional and opti￾mized models compared with the data (upper curve – LIS, lower – modulated to 550 MV). The lines are coded as in
Figure 7. — γ-ray spectra of hard electron spectrum model (500181); regions and coding as for
Figure 8. The fits to the observed γ-ray spectra are bet￾ter than for the conventional and hard electron spectrum models, both in the 1–10 GeV region and above 20 GeV. The spectra in different regions are satisfactorily repro￾duced and there is no longer a significant GeV excess. Hence the spectrum can now be reproduced from 30 MeV to 100 GeV. The proposed scenario implies a substantial contribution from IC at all energies, but especially be￾5 Berezhko & V¨olk (2000) did not address the question of regions away from the inner Galaxy
Figure 9. — Longitude profiles (|b| < 5.5 ◦) for optimized model (500190), compared with EGRET data in 10 energy ranges, 30 MeV – 10 GeV. Lines are coded as in
Figure 10. — Latitude profiles for optimized model (500190), inner Galaxy (330◦ < l < 30◦), compared with EGRET data in 12 energy ranges 30 MeV – 50 GeV. Lines are coded as in
Figure 11. — Latitude profiles for optimized model (500190), except inner Galaxy (30◦ < l < 330◦), compared with EGRET data in 10 energy ranges 30 MeV – 10 GeV. Lines are coded as in
Figure 12. — Electron spectra for conventional (solid) and opti￾mized models (dots). Upper curve – LIS, lower – modulated to 600 MV. Secondary electrons are shown separately for the same models. tions, gas distribution in the Galaxy, ISRF, systematic errors in the CR measurements, heliospheric modulation, etc. Some possible errors and their effects have been dis￾cussed in Moskalenko et al. (2002, 2003). Here we qual￾itatively mention what we think may affect our conclu￾sions and what may not. Mori (1997, and references therein) has re-evaluated the π 0 -production in pp-collisions using…

Showing 10 of the paper's 12 figures. See all figures in the PDF ↗ Figures reproduced from arXiv:astro-ph/0406254 by the authors.

References
How later papers used it

110 later papers cite it

Evidence: substantial. 103 of the citing papers share no author with it, and 39 of those carry a confident label from Pith's citation classifier.

Pith read 25 of the 156 citing sentences, choosing only those whose citation markers point at this paper alone, and wrote the paraphrases from them. The summary, the counts, the evidence level, the links and the quotes come from code; the summary counts the papers listed below. A later paper is listed under a label when Pith's citation classifier gives one of its sentences that label at probability 0.8 or more (0.95 for support), Pith's reading finds that the citing paper itself did it, and a second check confirms the paraphrase against the sentence. This section lists uses only. It does not list criticism or disagreement, so their absence here says nothing about how the paper has been received. Papers that share an author with this one are listed apart, under "By the same authors". A quote is the citing authors' own sentence, shown when their paper's licence is CC BY, CC BY-SA or CC0. JSON

Uses its method or tool 3 of 9 papers with this label · 1 by the same authors

By the same authors

  • Calculated cosmic-ray propagation and the corresponding diffuse emission using this paper's GALPROP code.

    H I spin temperature in the Fermi-LAT era (2010) (same authors) · 1.00 · paraphrased: the licence does not allow quoting

Uses its data or values 5 of 26 papers with this label
How the literature uses it

88 arXiv papers cite this in 120 sentences

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

Papers citing it · 67 sentences without a confident label

Uses its method or tool 6 papers
  • In this paper, CR propagation and corresponding diffuse emission is calculated using the GALPROP code [ref] .

    H I spin temperature in the Fermi-LAT era · Introduction · 1.00

  • We run GALPROP in the 2-dimensional mode and use the “conventional" (ID:500180) model [ref] for the galactic set-up and cosmic ray source/injection parameters; however, we vary the injection spectral index $γ_0$ for the electron background (parameterized as a power-law $E^-γ_0$) in the range $2.2-2.9$.

    Dark matter and pulsar signals for Fermi LAT, PAMELA, ATIC, HESS and WMAP data · Introduction · 0.90

  • Similar to the model used by the Fermi Collaboration, we use the GALPROP code [ref] to calculate the $γ$-ray and neutrino spectra as well as their morphology.

    Galactic Neutrinos in the TeV to PeV Range · Galactic Diffuse Gamma-Ray Emission · 1.00

Uses its data or values 20 papers
Compares with it 9 papers
Cites it as background 8 papers

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.