Cosmic rays and molecular clouds

This paper deals with the cosmic-ray penetration into molecular clouds and with the related gamma–ray emission. High energy cosmic rays interact with the dense gas and produce neutral pions which in turn decay into two gamma rays. This makes molecular clouds potential sources of gamma rays, especially if they are located in the vicinity of a powerful accelerator that injects cosmic rays in the interstellar medium. The amplitude and duration in time of the cosmic–ray overdensity around a given source depend on how quickly cosmic rays diffuse in the turbulent galactic magnetic field. For these reasons, gamma-ray observations of molecular clouds can be used both to locate the sources of cosmic rays and to constrain the properties of cosmic-ray diffusion in the Galaxy.

💡 Research Summary

This paper investigates how high‑energy cosmic rays (CRs) penetrate dense molecular clouds (MCs) and produce observable gamma‑ray emission, and it shows how such observations can be used to locate CR accelerators and to constrain the diffusion properties of CRs in the Galaxy. The authors begin by reviewing the standard picture: Galactic CRs are thought to be accelerated mainly in supernova remnants and other powerful sources, after which they diffuse through the turbulent interstellar magnetic field. The diffusion coefficient D(E) is usually modeled as D₀(E/E₀)^δ, where the index δ reflects the underlying turbulence spectrum (Kolmogorov, Kraichnan, etc.). Because the diffusion process is poorly constrained on local scales, an independent probe is needed.

Molecular clouds provide exactly that probe. Their gas densities (10²–10⁴ cm⁻³) are orders of magnitude higher than the average interstellar medium, so CRs colliding with cloud material efficiently produce neutral pions (π⁰). The π⁰ decay almost instantly into two gamma photons, making the gamma‑ray luminosity of a cloud directly proportional to the product of the CR density and the cloud’s gas density. Consequently, a cloud illuminated by an excess of CRs appears as a bright gamma‑ray source, even if the accelerator itself is not directly visible.

The authors formulate a time‑dependent, spherically symmetric diffusion model. A point‑like accelerator injects a total CR energy W_CR with a power‑law spectrum N₀(E)∝E⁻ˢ at t = 0. The Green’s‑function solution of the diffusion equation yields the CR density at distance r and time t:

n_CR(r,t,E) = (W_CR / (4πDt)^{3/2}) E⁻ˢ exp