Cosmic Rays in Bubbles: Why Heavy Nuclei Unexpectedly Rise with Energy ✛ ✎ ✦
Context
This work sits at the periphery of cosmic-ray origin studies. Nearby works cover galactic metal enrichment, neutrino tomography of the Earth, and black-hole winds, but there are almost no direct neighbors interpreting ATIC data. The closest arXiv papers (ATIC-2, CREAM) discuss nuclear spectra and secondary ratios but do not propose a closed-galaxy model. This work returns to the 1977 idea by Peters and Westergaard and applies it to new data.
Methods
The authors build a 'closed galaxy with bubbles' (CG+B) model: the galaxy is closed to diffusive escape, and all cosmic-ray sources reside in compact regions—bubbles created by supernovae. The Sun is located in a local bubble. The flux inside the bubble is a sum of a local component (standard leaky box) and a global equilibrium flux (bulk flux) described by a model with infinite escape length. The free parameter is the bulk flux fraction K. Calculations use path lengths from the book by Ginzburg and Syrovatskii, partial cross sections from a compilation, and the interstellar medium is taken as 90% protons and 10% helium. Flux ratios of nuclear groups (Ti+V+Cr)/Fe and H−/Fe (S–Cr) are compared with ATIC-2 data.
Results
The simple leaky-box model gives a monotonically decreasing ratio of secondary nuclei to iron and does not reproduce the rise above 100 GeV/nucleon. Adding primary fluxes (shaped like the oxygen spectrum) improves agreement below 100 GeV/nucleon, but the rise remains unexplained. The CG+B model with K=0.2 describes the data reasonably: the ratio first decreases, then rises due to competition between local and bulk fluxes. The model also predicts a universal hardening of primary nuclei spectra at 200–300 GeV/nucleon, consistent with ATIC and CREAM observations, but requires soft source spectra: an index of 2.55 for protons and 2.45 for oxygen. This may indicate a problem with the model. The estimated diffusion coefficient in the local bubble is 2–3 orders of magnitude lower than commonly accepted, assuming free escape of particles from the bubble boundary.
Implications
This work stands out by reviving the forgotten closed-galaxy model and applying it to new ATIC-2 data. Unlike studies focusing on individual element spectra or general propagation properties, it proposes a specific mechanism for the rise in heavy-to-iron ratios. The CG+B model explains the universal spectral hardening without invoking an additional soft source, but at the cost of soft source spectra themselves. This is the first attempt to quantitatively describe the rise within a local-bubble model.
Future development
The topic could develop toward clarifying the nature of local regions: whether they are spiral arms or supernova superbubbles. It is important to test the model against data from other experiments (AMS-02, CALET) and at higher energies. Alternative explanations for particle confinement in the bubble—reflection by walls or a combination with a low diffusion coefficient—also need investigation.
Impact
The results will affect models of cosmic-ray propagation in the Galaxy and the interpretation of secondary nuclear spectra. They are also important for understanding the role of local structures (bubbles) in shaping observed fluxes.
Next steps
The next step is to compare CG+B model predictions with data from new experiments on nuclear spectra and secondary ratios. It is also necessary to refine bubble parameters (size, density) and check whether agreement holds for other values of the bulk flux fraction K.
Key open problems
This work is connected to the unresolved problem of cosmic-ray origin and propagation in the Galaxy. The question of where particles are accelerated and how they leave their sources remains open. The closed-galaxy-with-bubbles model offers one possible answer but requires further testing.
The knowledge machine's view
This section is addressed to the author of the paper and their colleagues. Neighbouring works were found by meaning-based search across our archive, not by word overlap. We do not review: we show the field around the work.
What the work does
The paper revives the model of a closed Galaxy with local bubbles of cosmic-ray sources and applies it to ATIC-2 data. The authors compare the flux ratios of nuclear groups with the predictions of a simple 'leaky box' and the new CG+B model, showing that the latter describes the observed rise above 100 GeV/nucleon with a 20% contribution from the volume flux. The model also predicts a universal hardening of primary nuclear spectra at 200–300 GeV/nucleon, but requires soft source spectra with an index of ~2.5.
Where its strength is
A concrete mechanism for the rise in the ratios of heavy nuclei to iron through the competition between local and volume fluxes, which can explain the universal hardening of spectra without invoking an additional soft source.
Where it could go next
- It may be interesting to compare the proposed model with data on the spatial distribution of cosmic-ray sources, for example, with results on sky hot spots where ultra-high-energy particle emitters have been identified by anisotropy.based on: Sky Hotspots: The Senders of Ultra-Particles Identified arXiv:2510.06193
- Closely related is work on a bubble in a dark energy sea, which considers the dynamics of a local region with properties different from its surroundings; perhaps a similar approach to cosmic-ray bubble boundaries could refine the conditions for particle confinement.based on: Bubble in a Sea of Dark Energy arXiv:2607.18376v1
Worth trying
- Test how the CG+B model predictions change when accounting for the anisotropy of cosmic-ray arrival directions, using hot-spot analysis methods from work 2510.06193.
- Try an alternative description of the bubble boundary as a diffusion barrier, similar to the bubble boundary treatment in work 2607.18376v1, to assess the sensitivity of results to the assumption of free particle escape.
Why this may matter
If the closed Galaxy with bubbles model is confirmed, it will change our understanding of cosmic-ray propagation in the Galaxy and may explain the universal spectral hardening without additional sources, which is important for interpreting data from current and future experiments.
Nearby works in our archive
- Triple Thread: How Vector Dark Matter Gives Birth to Invisible Starssimilarity 0.65
- Bubble in a Sea of Dark Energysimilarity 0.65
- Dark Matter and the Dance of Stellar Swarmssimilarity 0.63
- Sky Hotspots: The Senders of Ultra-Particles Identifiedsimilarity 0.62
- Red dots and the magnetic cocoon of a black holesimilarity 0.62
- Gas Clouds Dance Around the Galactic Centersimilarity 0.62
The closest work in the archive sits at 0.65 by meaning — the usual level for our corpus.
The local bubble, which hosts the Sun, is a real region of tenuous hot gas about 100 parsecs across, created by supernova explosions. If the CG+B model is correct, we literally live inside a 'cauldron' where cosmic rays are brewed.
Key numbers
- rise energy: 100 GeV/nucleon
- predicted hardening: 200–300 GeV/nucleon
- bulk flux fraction K: 0.2
- proton spectral index: 2.55
- oxygen spectral index: 2.45