The Interaction of Cosmic Rays with Diffuse Clouds

We study the change in cosmic-ray pressure, the change in cosmic-ray density, and the level of cosmic-ray induced heating via Alfven-wave damping when cosmic rays move from a hot ionized plasma to a cool cloud embedded in that plasma. The general analysis method outlined here can apply to diffuse clouds in either the ionized interstellar medium or in galactic winds. We introduce a general-purpose model of cosmic-ray diffusion building upon the hydrodynamic approximation for cosmic rays (from McKenzie & Voelk and Breitschwerdt and collaborators). Our improved method self-consistently derives the cosmic-ray flux and diffusivity under the assumption that the streaming instability is the dominant mechanism for setting the cosmic-ray flux and diffusion. We find that, as expected, cosmic rays do not couple to gas within cool clouds (cosmic rays exert no forces inside of cool clouds), that the cosmic-ray density does not increase within clouds (it may slightly decrease in general, and decrease by an order of magnitude in some cases), and that cosmic-ray heating (via Alfven-wave damping and not collisional effects as for ~10 MeV cosmic rays) is only important under the conditions of relatively strong (10 micro-Gauss) magnetic fields or high cosmic-ray pressure (~10^{-11} ergs cm^{-3}).

💡 Research Summary

**
This paper investigates how ∼1 GeV cosmic rays (CRs) behave when they travel from a hot, ionized plasma (∼10⁶ K) into an embedded cool cloud (∼10²–10³ K). The authors develop a self‑consistent model of CR diffusion that is built on the cosmic‑ray hydrodynamic framework originally formulated by McKenzie & Völk and later extended by Breitschwerdt and collaborators. The key physical ingredient is the streaming instability: when CRs drift faster than the local Alfvén speed, they resonantly excite Alfvén waves whose growth is balanced by various damping processes (ion‑neutral friction, non‑linear wave–wave interactions, etc.). By solving the coupled equations for CR pressure, Alfvén‑wave pressure, and the effective diffusion coefficient κ_CR, the authors determine the CR flux and diffusivity as functions of the local plasma conditions.

In the hot ionized phase the streaming instability is active, limiting the bulk CR drift to roughly the Alfvén speed (v_A ≈ 20 km s⁻¹) and producing a short mean free path (λ_mfp ≈ 0.01–0.1 pc). When CRs encounter the cool cloud, the ionization fraction drops sharply, Alfvén waves are rapidly damped, and the instability shuts off. Consequently CRs become essentially free‑streaming, moving at speeds close to c. The CR pressure inside the cloud therefore remains comparable to, or slightly lower than, the upstream value; typical reductions are modest (∼20 % in most cases) but can reach an order‑of‑magnitude decrease for extreme parameter choices. Thus, contrary to earlier “exclusion” scenarios that predicted a strong depletion of CRs inside clouds, the present model finds only a mild reduction for GeV‑energy particles.

The authors also evaluate heating by Alfvén‑wave damping. The heating rate scales with the wave energy density, which in turn depends on the magnetic field strength and the upstream CR pressure. For typical interstellar values (B ≈ 3 µG, P_CR ≈ 10⁻¹² erg cm⁻³) the damping‑induced heating is negligible compared with radiative cooling. Only when the magnetic field is relatively strong (B ≳ 10 µG) and the CR pressure is high (P_CR ≳ 10⁻¹¹ erg cm⁻³) does the heating become comparable to cooling, potentially affecting the thermal balance at the cloud boundary.

The paper places its results in context with earlier work. Skilling & Strong (1976) and Cesarsky & Völk (1977) showed that low‑energy CRs (≲100 MeV) are efficiently excluded from clouds by ion‑neutral damping; the present study reproduces this behavior for low‑energy particles but demonstrates that GeV CRs are largely unaffected. Padoan & Scalo (2005) had suggested that CR density could increase as n_CR ∝ n^{1/2} inside dense clouds; the current analysis, which treats the streaming instability self‑consistently, finds the opposite for GeV CRs—no significant density enhancement, and sometimes a decrease.

Finally, the authors consider implications for multiphase galactic winds, where cool clouds are observed moving at hundreds of km s⁻¹ within a hot outflow. Since CRs do not exert a force inside the cool phase (they are not coupled to the gas), they cannot directly accelerate or disrupt the clouds. However, if the wind’s magnetic field is amplified near the cloud surface, Alfvén‑wave damping could provide localized heating that might alter cloud survival timescales. In typical wind conditions, this effect is expected to be minor.

In summary, the study concludes that for realistic interstellar and galactic‑wind environments, GeV cosmic rays pass through diffuse cool clouds with little change in pressure or density, exert no dynamical force inside the clouds, and contribute appreciable heating only under unusually strong magnetic fields or high CR pressures. This clarifies the role of CRs in multiphase media and provides a framework for future observational tests, such as gamma‑ray measurements of cloud‑embedded CR populations or diagnostics of CR‑driven heating at cloud boundaries.