An official website of the United States government
Abstract We study the propagation of mildly relativistic cosmic rays (CRs) in multiphase interstellar medium environments with conditions typical of nearby disk galaxies. We employ the techniques developed in Armillotta et al. to postprocess three high-resolution TIGRESS magnetohydrodynamic simulations modeling local patches of star-forming galactic disks. Together, the three simulations cover a wide range of gas surface density, gravitational potential, and star formation rate (SFR). Our prescription for CR propagation includes the effects of advection by the background gas, streaming along the magnetic field at the local ion Alfvén speed, and diffusion relative to the Alfvén waves, with the diffusion coefficient set by the balance between streaming-driven Alfvén wave excitation and damping mediated by local gas properties. We find that the combined transport processes are more effective in environments with higher SFR. These environments are characterized by higher-velocity hot outflows (created by clustered supernovae) that rapidly advect CRs away from the galactic plane. As a consequence, the ratio of midplane CR pressure to midplane gas pressures decreases with increasing SFR. We also use the postprocessed simulations to make predictions regarding the potential dynamical impacts of CRs. The relatively flat CR pressure profiles near the midplane argue that they would not provide significant support against gravity for most of the ISM mass. However, the CR pressure gradients are larger than the other pressure gradients in the extraplanar region (∣ z ∣ > 0.5 kpc), suggesting that CRs may affect the dynamics of galactic fountains and/or winds. The degree of this impact is expected to increase in environments with lower SFR.
more »
« less
This content will become publicly available on July 30, 2026
Cosmic-Ray Diffusion in the Turbulent Interstellar Medium: Effects of Mirror Diffusion and Pitch-angle Scattering
Abstract Cosmic rays (CRs) interact with turbulent magnetic fields in the interstellar medium (ISM), generating nonthermal emission. After many decades of studies, the theoretical understanding of their diffusion in the ISM continues to pose a challenge. This study numerically explores a recent prediction termed “mirror diffusion” and its synergy with the traditional diffusion mechanism based on gyroresonant scattering. Our study combines 3D MHD simulations of star-forming regions with test particle simulations to analyze CR diffusion. We demonstrate the significance of mirror diffusion in CR diffusion parallel to the magnetic field when the mirroring condition is satisfied. Our results support the theoretical expectation that the resulting particle propagation arising from mirror diffusion in combination with much faster diffusion induced by gyroresonant scattering resembles a Levy-flight-like propagation. Our study highlights the necessity to reevaluate the diffusion coefficients traditionally adopted in the ISM based on gyroresonant scattering alone. For instance, our simulations imply a diffusion coefficient ∼1027cm2s–1for particles with a few hundred TeV within regions spanning a few parsecs around the source. This estimate is in agreement with gamma-ray observations, which show the relevance of our results for the understanding of gamma-ray emission in star-forming regions.
more »
« less
- Award ID(s):
- 2307840
- PAR ID:
- 10632553
- Publisher / Repository:
- Astrophysical Journal
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 988
- Issue:
- 2
- ISSN:
- 0004-637X
- Page Range / eLocation ID:
- 269
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Geminga is the first pulsar around which a remarkable gamma-ray halo extending over a few degrees was discovered at TeV energies by MILAGRO and HAWC and later by H.E.S.S., and byFermi-LAT in the GeV band. More middle-aged pulsars have exhibited gamma-ray halos, and they are now recognised as an emerging class of Galactic gamma-ray sources. The emission appears in the late evolution stage of pulsars, and is most plausibly explained by inverse Compton scattering of CMB and interstellar photons by relativistic electrons and positrons escaping from the pulsar wind nebulae. These observations pose a number of theoretical challenges, particularly the origin of the inferred, significantly lower effective diffusion coefficients around the pulsar when compared to typical Galactic values. Tackling these questions requires constraining the ambient magnetic field properties, which can be achieved through X-ray observations. If the gamma-ray halos originate from a distribution of highly energetic electrons, synchrotron losses in the ambient magnetic fields of the same particles are expected to produce a diffuse X-ray emission with a similar spatial extension. We present the most comprehensive X-ray study of the Geminga pulsar halo to date, utilising archival data fromXMM-NewtonandNuSTAR. Our X-ray analysis covers a broad bandwidth (0.5 − 79 keV) and large field of view (θ ∼ 4°) for the first time. This was achieved by accurately measuring the background over the entire field of view, and taking into account both focused and stray-light X-ray photons from the pulsar halo withNuSTAR. We find no significant emission and set robust constraints on the X-ray halo flux. These are translated to stringent constraints on the ambient magnetic field strength and the diffusion coefficient by using a physical model considering particle injection, diffusion, and cooling over the pulsar’s lifetime, which is tuned by fitting multi-wavelength data. Our novel methodology for modelling and searching for synchrotron X-ray halos can be applied to other pulsar halo candidates.more » « less
-
Abstract We present the implementation and the first results of cosmic ray (CR) feedback in the Feedback In Realistic Environments (FIRE) simulations. We investigate CR feedback in non-cosmological simulations of dwarf, sub-L⋆ starburst, and L⋆ galaxies with different propagation models, including advection, isotropic and anisotropic diffusion, and streaming along field lines with different transport coefficients. We simulate CR diffusion and streaming simultaneously in galaxies with high resolution, using a two moment method. We forward-model and compare to observations of γ-ray emission from nearby and starburst galaxies. We reproduce the γ-ray observations of dwarf and L⋆ galaxies with constant isotropic diffusion coefficient κ ∼ 3 × 1029 cm2 s−1. Advection-only and streaming-only models produce order-of-magnitude too large γ-ray luminosities in dwarf and L⋆ galaxies. We show that in models that match the γ-ray observations, most CRs escape low-gas-density galaxies (e.g. dwarfs) before significant collisional losses, while starburst galaxies are CR proton calorimeters. While adiabatic losses can be significant, they occur only after CRs escape galaxies, so they are only of secondary importance for γ-ray emissivities. Models where CRs are “trapped” in the star-forming disk have lower star formation efficiency, but these models are ruled out by γ-ray observations. For models with constant κ that match the γ-ray observations, CRs form extended halos with scale heights of several kpc to several tens of kpc.more » « less
-
ABSTRACT The nature of cosmic ray (CR) transport in the Milky Way remains elusive. The predictions of current microphysical CR transport models in magnetohydrodynamic (MHD) turbulence are drastically different from what is observed. These models usually focus on MHD turbulence with a strong guide field and ignore the impact of turbulent intermittency on particle propagation. This motivates our studying the alternative regime of large-amplitude turbulence with δB/B0 ≫ 1, in which intermittent small-scale magnetic field reversals are ubiquitous. We study particle transport in such turbulence by integrating trajectories in stationary snapshots. To quantify spatial diffusion, we use a set-up with continuous particle injection and escape, which we term the turbulent leaky box. We find that particle transport is very different from the strong guide-field case. Low-energy particles are better confined than high-energy particles, despite less efficient pitch-angle isotropization at small energies. In the limit of weak guide field, energy-dependent confinement is driven by the energy-dependent (in)ability to follow reversing magnetic field lines exactly and by the scattering in regions of ‘resonant curvature’, where the field line bends on a scale that is of the order of the local particle gyro-radius. We derive a heuristic model of particle transport in magnetic folds that approximately reproduces the energy dependence of transport found numerically. We speculate that CR propagation in the Galaxy is regulated by the intermittent field reversals highlighted here and discuss the implications of our findings for CR transport in the Milky Way.more » « less
-
null (Ed.)ABSTRACT Blazars emit a highly variable non-thermal spectrum. It is usually assumed that the same non-thermal electrons are responsible for the IR-optical-UV emission (via synchrotron) and the gamma-ray emission (via inverse Compton). Hence, the light curves in the two bands should be correlated. Orphan gamma-ray flares (i.e. lacking a luminous low-frequency counterpart) challenge our theoretical understanding of blazars. By means of large-scale two-dimensional radiative particle-in-cell simulations, we show that orphan gamma-ray flares may be a self-consistent by-product of particle energization in turbulent magnetically dominated pair plasmas. The energized particles produce the gamma-ray flare by inverse Compton scattering an external radiation field, while the synchrotron luminosity is heavily suppressed since the particles are accelerated nearly along the direction of the local magnetic field. The ratio of inverse Compton to synchrotron luminosity is sensitive to the initial strength of turbulent fluctuations (a larger degree of turbulent fluctuations weakens the anisotropy of the energized particles, thus increasing the synchrotron luminosity). Our results show that the anisotropy of the non-thermal particle population is key to modelling the blazar emission.more » « less
