Observations of the $gamma$-ray emission around star clusters, isolated supernova remnants, and pulsar wind nebulae indicate that the cosmic-ray (CR) diffusion coefficient near acceleration sites can be suppressed by a large factor compared to the Ga
laxy average. We explore the effects of such local suppression of CR diffusion on galaxy evolution using simulations of isolated disk galaxies with regular and high gas fractions. Our results show that while CR propagation with constant diffusivity can make gaseous disks more stable by increasing the midplane pressure, large-scale CR pressure gradients cannot prevent local fragmentation when the disk is unstable. In contrast, when CR diffusivity is suppressed in star-forming regions, the accumulation of CRs in these regions results in strong local pressure gradients that prevent the formation of massive gaseous clumps. As a result, the distribution of dense gas and star formation changes qualitatively: a globally unstable gaseous disk does not violently fragment into massive star-forming clumps but maintains a regular grand-design spiral structure. This effect regulates star formation and disk structure and is qualitatively different from and complementary to the global role of CRs in vertical hydrostatic support of the gaseous disk and in driving galactic winds.
We study diffusive shock acceleration (DSA) of electrons in non-relativistic quasi-perpendicular shocks using self-consistent one-dimensional particle-in-cell (PIC) simulations. By exploring the parameter space of sonic and Alfv{e}nic Mach numbers we
find that high Mach number quasi-perpendicular shocks can efficiently accelerate electrons to power-law downstream spectra with slopes consistent with DSA prediction. Electrons are reflected by magnetic mirroring at the shock and drive non-resonant waves in the upstream. Reflected electrons are trapped between the shock front and upstream waves and undergo multiple cycles of shock drift acceleration before the injection into DSA. Strong current-driven waves also temporarily change the shock obliquity and cause mild proton pre-acceleration even in quasi-perpendicular shocks, which otherwise do not accelerate protons. These results can be used to understand nonthermal emission in supernova remnants and intracluster medium in galaxy clusters.
We study the thermalization, injection, and acceleration of ions with different mass/charge ratios, $A/Z$, in non-relativistic collisionless shocks via hybrid (kinetic ions-fluid electrons) simulations. In general, ions thermalize to a post-shock tem
perature proportional to $A$. When diffusive shock acceleration is efficient, ions develop a non-thermal tail whose extent scales with $Z$ and whose normalization is enhanced as $(A/Z)^2$, so that incompletely-ionized heavy ions are preferentially accelerated. We discuss how these findings can explain observed heavy-ion enhancements in Galactic cosmic rays.
We formulate a magnetohydrodynamic-particle-in-cell (MHD-PIC) method for describing the interaction between collisionless cosmic ray (CR) particles and a thermal plasma. The thermal plasma is treated as a fluid, obeying equations of ideal MHD, while
CRs are treated as relativistic Lagrangian particles subject to the Lorentz force. Backreaction from CRs to the gas is included in the form of momentum and energy feedback. In addition, we include the electromagnetic feedback due to CR-induced Hall effect that becomes important when the electron-ion drift velocity of the background plasma induced by CRs approaches the Alfven velocity. Our method is applicable on scales much larger than the ion inertial length, bypassing the microscopic scales that must be resolved in conventional PIC methods, while retaining the full kinetic nature of the CRs. We have implemented and tested this method in the Athena MHD code, where the overall scheme is second-order accurate and fully conservative. As a first application, we describe a numerical experiment to study particle acceleration in non-relativistic shocks. Using a simplified prescription for ion injection, we reproduce the shock structure and the CR energy spectra obtained with more self-consistent hybrid-PIC simulations, but at substantially reduced computational cost. We also show that the CR-induced Hall effect reduces the growth rate of the Bell instability and affects the gas dynamics in the vicinity of the shock front. As a step forward, we are able to capture the transition of particle acceleration from non relativistic to relativistic regimes, with momentum spectrum $f(p)sim p^{-4}$ connecting smoothly through the transition, as expected from the theory of Fermi acceleration.
Particle acceleration at non-relativistic shocks can be very efficient, leading to the appearance of non-linear effects due to the dynamical reaction of the accelerated particles on the shock structure and to the non-linear amplification of the magne
tic field in the shock vicinity. The value of the maximum momentum $p_{max}$ in these circumstances cannot be estimated using the classical results obtained within the framework of test particle approaches. We provide here the first attempt at estimating $p_{max}$ in the cosmic ray modified regime, taking into account the non-linear effects mentioned above.
