An official website of the United States government
Abstract Cosmic-ray transport on galactic scales depends on the detailed properties of the magnetized, multiphase interstellar medium (ISM). In this work, we postprocess a high-resolution TIGRESS magnetohydrodynamic simulation modeling a local galactic disk patch with a two-moment fluid algorithm for cosmic-ray transport. We consider a variety of prescriptions for the cosmic rays, from a simple, purely diffusive formalism with constant scattering coefficient, to a physically motivated model in which the scattering coefficient is set by the critical balance between streaming-driven Alfvén wave excitation and damping mediated by local gas properties. We separately focus on cosmic rays with kinetic energies of ∼1 GeV (high-energy) and ∼30 MeV (low energy), respectively important for ISM dynamics and chemistry. We find that simultaneously accounting for advection, streaming, and diffusion of cosmic rays is crucial for properly modeling their transport. Advection dominates in the high-velocity, low-density hot phase, while diffusion and streaming are more important in higher-density, cooler phases. Our physically motivated model shows that there is no single diffusivity for cosmic-ray transport: the scattering coefficient varies by four or more orders of magnitude, maximal at density n H ∼ 0.01 cm −3 . The ion-neutral damping of Alfvén waves results in strong diffusion and nearly uniform cosmic-ray pressure within most of the mass of the ISM. However, cosmic rays are trapped near the disk midplane by the higher scattering rate in the surrounding lower-density, higher-ionization gas. The transport of high-energy cosmic rays differs from that of low-energy cosmic rays, with less effective diffusion and greater energy losses for the latter.
more »
« less
Cosmic-Ray Feedback on Bistable Interstellar Medium Turbulence
Abstract While cosmic rays (E≳ 1 GeV) are well coupled to a galaxy’s interstellar medium (ISM) at scales ofL> 100 pc, adjusting stratification and driving outflows, their impact on small scales is less clear. Based on calculations of the cosmic-ray diffusion coefficient from observations of the grammage in the Milky Way, cosmic rays have little time to dynamically impact the ISM on those small scales. Using numerical simulations, we explore how more complex cosmic-ray transport could allow cosmic rays to couple to the ISM on small scales. We create a two-zone model of cosmic-ray transport, with the cosmic-ray diffusion coefficient set at the estimated Milky Way value in cold gas but smaller in warm gas. We compare this model to simulations with a constant diffusion coefficient. Quicker diffusion through cold gas allows more cold gas to form compared to a simulation with a constant, small diffusion coefficient. However, slower diffusion in warm gas allows cosmic rays to take energy from the turbulent cascade anisotropically. This cosmic-ray energization comes at the expense of turbulent energy which would otherwise be lost during radiative cooling. Finally, we show our two-zone model is capable of matching observational estimates of the grammage for some transport paths through the simulation.
more »
« less
- Award ID(s):
- 2007323
- PAR ID:
- 10546129
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 974
- Issue:
- 1
- ISSN:
- 0004-637X
- Format(s):
- Medium: X Size: Article No. 17
- Size(s):
- Article No. 17
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract We use hydrodynamical simulations of two Milky Way–mass galaxies to demonstrate the impact of cosmic-ray pressure on the kinematics of cool and warm circumgalactic gas. Consistent with previous studies, we find that cosmic-ray pressure can dominate over thermal pressure in the inner 50 kpc of the circumgalactic medium (CGM), creating an overall cooler CGM than that of similar galaxy simulations run without cosmic rays. We generate synthetic sight lines of the simulated galaxies’ CGM and use Voigt profile-fitting methods to extract ion column densities, Doppler-bparameters, and velocity centroids of individual absorbers. We directly compare these synthetic spectral line fits with HST/COS CGM absorption-line data analyses, which tend to show that metallic species with a wide range of ionization potential energies are often kinematically aligned. Compared to the Milky Way simulation run without cosmic rays, the presence of cosmic-ray pressure in the inner CGM creates narrower Oviabsorption features and broader Siiiiabsorption features, a quality that is more consistent with observational data. Additionally, because the cool gas is buoyant due to nonthermal cosmic-ray pressure support, the velocity centroids of both cool and warm gas tend to align in the simulated Milky Way with feedback from cosmic rays. Our study demonstrates that detailed, direct comparisons between simulations and observations, focused on gas kinematics, have the potential to reveal the dominant physical mechanisms that shape the CGM.more » « less
-
We use hydrodynamical simulations of two Milky Way-mass galaxies to demonstrate the impact of cosmic-ray pressure on the kinematics of cool and warm circumgalactic gas. Consistent with previous studies, we find that cosmic-ray pressure can dominate over thermal pressure in the inner 50 kpc of the circumgalactic medium (CGM), creating an overall cooler CGM than that of similar galaxy simulations run without cosmic rays. We generate synthetic sightlines of the simulated galaxies' CGM and use Voigt profile fitting methods to extract ion column densities, Doppler-b parameters, and velocity centroids of individual absorbers. We directly compare these synthetic spectral line fits with HST/COS CGM absorption-line data analyses, which tend to show that metallic species with a wide range of ionization potential energies are often kinematically aligned. Compared to the Milky-Way simulation run without cosmic rays, the presence of cosmic-ray pressure in the inner CGM creates narrower OVI absorption features and broader SiIII absorption features, a quality which is more consistent with observational data. Additionally, because the cool gas is buoyant due to nonthermal cosmic-ray pressure support, the velocity centroids of both cool and warm gas tend to align in the simulated Milky Way with feedback from cosmic rays. Our study demonstrates that detailed, direct comparisons between simulations and observations, focused on gas kinematics, have the potential to reveal the dominant physical mechanisms that shape the CGM.more » « less
-
Abstract Cosmic-ray transport in astrophysical environments is often dominated by the diffusion of particles in a magnetic field composed of both a turbulent and a mean component. This process, which is two-fold turbulent mixing in that the particle motion is stochastic with respect to the field lines, needs to be understood in order to properly model cosmic-ray signatures. One of the most important aspects in the modeling of cosmic-ray diffusion is that fully resonant scattering, the most effective such process, is only possible if the wave spectrum covers the entire range of propagation angles. By taking the wave spectrum boundaries into account, we quantify cosmic-ray diffusion parallel and perpendicular to the guide field direction at turbulence levels above 5% of the total magnetic field. We apply our results of the parallel and perpendicular diffusion coefficient to the Milky Way. We show that simple purely diffusive transport is in conflict with observations of the inner Galaxy, but that just by taking a Galactic wind into account, data can be matched in the central 5 kpc zone. Further comparison shows that the outer Galaxy at $$>5$$ > 5 kpc, on the other hand, should be dominated by perpendicular diffusion, likely changing to parallel diffusion at the outermost radii of the Milky Way.more » « less
-
Abstract In a seminal paper, Parker showed the vertical stratification of the interstellar medium (ISM) is unstable if magnetic fields and cosmic rays provide too large a fraction of pressure support. Cosmic ray acceleration is linked to star formation, so Parker’s instability and its nonlinear outcomes are a type of star formation feedback. Numerical simulations have shown the instability can significantly restructure the ISM, thinning the thermal gas layer and thickening the magnetic field and cosmic ray layer. However, the timescale on which this occurs is rather long (∼0.4 Gyr). Furthermore, the conditions for instability depend on the model adopted for cosmic ray transport. In this work, we connect the instability and feedback problems by examining the effect of a single, spatially and temporally localized cosmic ray injection on the ISM over ∼1 kpc3scales. We perform cosmic ray magnetohydrodynamic simulations using theAthena++code, varying the background properties, dominant cosmic ray transport mechanism, and injection characteristics between our simulation runs. We find robust effects of buoyancy for all transport models, with disruption of the ISM on timescales as short as 100 Myr when the background equilibrium is dominated by cosmic ray pressure.more » « less
