Search for: All records

Abstract Supernova energy drives interstellar medium (ISM) turbulence and can help launch galactic winds. What difference does it make if 10% of the energy is initially deposited into cosmic rays? To help answer this question and study cosmic-ray feedback, we perform galactic patch simulations of a stratified ISM in a low star formation rate, high magnetic field regime. We compare two magnetohydrodynamic and cosmic-ray simulations, which are identical except for how each supernova’s energy is injected. In one, 10% of the energy is injected as cosmic-ray energy. In the other case, energy injection is strictly thermal and kinetic. We find that cosmic-ray injections drive a faster, hotter, and more massive outflow long after the injections occur. Both simulations show the formation of cold clouds (with a total mass fraction > 50%) through the Parker instability and thermal instability. The Parker instability simultaneously produces high mass loading factorsη> 10³as it does not directly require star formation and supernovae. We also show how the Parker instability naturally leads to a decorrelation of cosmic-ray pressure and gas density. This decorrelation leads to a significant decrease in the calorimetric fraction for injected cosmic rays, but it depends on having a highly resolved magnetic field. 

Abstract We investigate the 1D plane-parallel front connecting the warm (10⁴K) and hot (10⁶K) phases of the circumgalactic medium (CGM), focusing on the influence of cosmic rays (CRs) in shaping these transition layers. We find that cosmic rays dictate the thermal balance while other fluxes (thermal conduction, radiative cooling, and gas flow) adjust to compensate. We compute column densities and ratios for the transition-temperature ions Siiv, Civ, Ovi, and Nv, and compare them with observational data. While most models struggle to simultaneously reproduce the observed Siiv/Civand Civ/Oviratios, a subset with intermediate magnetic field strength (e.g.,B= 20μG) shows overlap with the data, although we make no claims for their uniqueness. These discrepancies suggest that the models perform better at reproducing higher-temperature ions but underestimate the contribution from cooler, photoionized regions. Compared to models without CRs, CR-mediated fronts in sufficiently strong magnetic fields produce broader transition layers and higher ion ratios, indicating that CRs can significantly alter the thermal and ionization structure of the CGM. Our results suggest that CR heating may help explain some observed ion columns under specific conditions, though additional physics may be needed for full agreement with observations. 

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. 

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 kpc³scales. 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. 

Abstract Cosmic rays have been shown to be extremely important in the dynamics of diffuse gas in galaxies, helping to maintain hydrostatic equilibrium, and serving as a regulating force in star formation. In this paper, we address the influence of cosmic rays on galaxies by re-examining the theory of a cosmic ray Eddington limit, first proposed by Socrates et al. and elaborated upon by Crocker et al. and Huang & Davis. A cosmic ray Eddington limit represents a maximum cosmic ray energy density above which the interstellar gas cannot be in hydrostatic equilibrium, resulting in a wind. In this paper, we continue to explore the idea of a cosmic ray Eddington limit by introducing a general framework that accounts for the circumgalactic environment and applying it to five galaxies that we believe to be a good representative sample of the star-forming galaxy population, using different cosmic ray transport models to determine what gives each galaxy the best chance to reach this limit. We show that, while an Eddington limit for cosmic rays does exist, for our five galaxies, the limit either falls at star formation rates that are much larger or gas densities that are much lower than each galaxy’s measured values. This suggests that cosmic ray pressure is not the main factor limiting the luminosity of starburst galaxies. 

ABSTRACT Understanding the time-scales for diffusive processes and their degree of anisotropy is essential for modelling cosmic ray transport in turbulent magnetic fields. We show that the diffusion time-scales are isotropic over a large range of energy and turbulence levels, notwithstanding the high degree of anisotropy exhibited by the components of the diffusion tensor for cases with an ordered magnetic field component. The predictive power of the classical scattering relation as a description for the relation between the parallel and perpendicular diffusion coefficients is discussed and compared to numerical simulations. Very good agreement for a large parameter space is found, transforming classical scattering relation predictions into a computational prescription for the perpendicular component. We discuss and compare these findings, in particular, the time-scales to become diffusive with the time-scales that particles reside in astronomical environments, the so-called escape time-scales. The results show that, especially at high energies, the escape times obtained from diffusion coefficients may exceed the time-scales required for diffusion. In these cases, the escape time cannot be determined by the diffusion coefficients. 

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.