More Like this

1. Standard self-confinement and extrinsic turbulence models for cosmic ray transport are fundamentally incompatible with observations

   https://doi.org/10.1093/mnras/stac2909  

   Hopkins, Philip F.; Squire, Jonathan; Butsky, Iryna S.; Ji, Suoqing (October 2022, Monthly Notices of the Royal Astronomical Society)

   ABSTRACT Models for cosmic ray (CR) dynamics fundamentally depend on the rate of CR scattering from magnetic fluctuations. In the ISM, for CRs with energies ∼MeV-TeV, these fluctuations are usually attributed either to ‘extrinsic turbulence’ (ET) – a cascade from larger scales – or ‘self-confinement’ (SC) – self-generated fluctuations from CR streaming. Using simple analytic arguments and detailed ‘live’ numerical CR transport calculations in galaxy simulations, we show that both of these, in standard form, cannot explain even basic qualitative features of observed CR spectra. For ET, any spectrum that obeys critical balance or features realistic anisotropy, or any spectrum that accounts for finite damping below the dissipation scale, predicts qualitatively incorrect spectral shapes and scalings of B/C and other species. Even if somehow one ignored both anisotropy and damping, observationally required scattering rates disagree with ET predictions by orders of magnitude. For SC, the dependence of driving on CR energy density means that it is nearly impossible to recover observed CR spectral shapes and scalings, and again there is an orders-of-magnitude normalization problem. But more severely, SC solutions with super-Alfvénic streaming are unstable. In live simulations, they revert to either arbitrarily rapid CR escape with zero secondary production, or to bottleneck solutions with far-too-strong CR confinement and secondary production. Resolving these fundamental issues without discarding basic plasma processes requires invoking different drivers for scattering fluctuations. These must act on a broad range of scales with a power spectrum obeying several specific (but plausible) constraints. 

   more » « less
2. Spectrum and Location of Ongoing Extreme Particle Acceleration in Cassiopeia A

   https://doi.org/10.3847/1538-4357/ad98f1  

   Woo, Jooyun; Mori, Kaya; Hailey, Charles_J; Spira-Savett, Elizabeth; Bamba, Aya; Grefenstette, Brian_W; Humensky, Thomas_B; Mukherjee, Reshmi; Safi-Harb, Samar; Temim, Tea; et al (January 2025, The Astrophysical Journal)

   Abstract Young supernova remnants (SNRs) are believed to be the origin of energetic cosmic rays (CRs) below the “knee” of their spectrum at ∼3 PeV (10¹⁵eV). Nevertheless, the precise location, duration, and operation of CR acceleration in young SNRs are open questions. Here, we report on multiepoch X-ray observations of Cassiopeia A (Cas A), a 350 yr old SNR, in the 15–50 keV band that probes the most energetic CR electrons. The observed X-ray flux decrease (15% ± 1% over 10 yr), contrary to the expected >90% decrease based on previous radio, X-ray, and gamma-ray observations, provides unambiguous evidence for CR electron acceleration operating in Cas A. A temporal model for the radio and X-ray data accounting for electron cooling and continuous injection finds that the freshly injected electron spectrum is significantly harder (exponential cutoff power-law indexq= 2.15), and its cutoff energy is much higher (E_cut = 36 TeV), than the relic electron spectrum (q = 2.44 ± 0.03,E_cut = 4 ± 1 TeV). Both electron spectra are naturally explained by the recently developed modified nonlinear diffusive shock acceleration (mNLDSA) mechanism. The CR protons producing the observed gamma rays are likely accelerated at the same location by the same mechanism as the injected electrons. The Cas A observations and spectral modeling represent the first time radio, X-ray, gamma-ray, and CR spectra have been self-consistently tied to a specific acceleration mechanism—mNLDSA—in a young SNR. 

   more » « less
3. Particle Reacceleration by Turbulence and Radio Constraints on Multimessenger High-energy Emission from the Coma Cluster

   https://doi.org/10.3847/1538-4357/ac1cdb  

   Nishiwaki, Kosuke; Asano, Katsuaki; Murase, Kohta (November 2021, The Astrophysical Journal)

   Abstract Galaxy clusters are considered to be gigantic reservoirs of cosmic rays (CRs). Some of the clusters are found with extended radio emission, which provides evidence for the existence of magnetic fields and CR electrons in the intra-cluster medium. The mechanism of radio halo (RH) emission is still under debate, and it has been believed that turbulent reacceleration plays an important role. In this paper, we study the reacceleration of CR protons and electrons in detail by numerically solving the Fokker–Planck equation, and show how radio and gamma-ray observations can be used to constrain CR distributions and resulting high-energy emission for the Coma cluster. We take into account the radial diffusion of CRs and follow the time evolution of their one-dimensional distribution, by which we investigate the radial profile of the CR injection that is consistent with the observed RH surface brightness. We find that the required injection profile is nontrivial, depending on whether CR electrons have a primary or secondary origin. Although the secondary CR electron scenario predicts larger gamma-ray and neutrino fluxes, it is in tension with the observed RH spectrum for hard injection indexes, α < 2.45. This tension is relaxed if the turbulent diffusion of CRs is much less efficient than the fiducial model, or the reacceleration is more efficient for lower-energy CRs. In both the secondary and primary scenario, we find that galaxy clusters can make a sizable contribution to the all-sky neutrino intensity if the CR energy spectrum is nearly flat. 

   more » « less
4. Turbulent Reacceleration of Streaming Cosmic Rays

   https://doi.org/10.3847/1538-4357/aca021  

   Bustard, Chad; Oh, S. Peng (December 2022, The Astrophysical Journal)

   Abstract Subsonic, compressive turbulence transfers energy to cosmic rays (CRs), a process known as nonresonant reacceleration. It is often invoked to explain the observed ratios of primary to secondary CRs at ∼GeV energies, assuming wholly diffusive CR transport. However, such estimates ignore the impact of CR self-confinement and streaming. We study these issues in stirring box magnetohydrodynamic (MHD) simulations using Athena++, with field-aligned diffusive and streaming CR transport. For diffusion only, we find CR reacceleration rates in good agreement with analytic predictions. When streaming is included, reacceleration rates depend on plasmaβ. Due to streaming-modified phase shifts between CR and gas variables, they are slower than canonical reacceleration rates in low-βenvironments like the interstellar medium but remain unchanged in high-βenvironments like the intracluster medium. We also quantify the streaming energy-loss rate in our simulations. For sub-Alfvénic turbulence, it is resolution dependent (hence unconverged in large-scale simulations) and heavily suppressed compared to the isotropic loss ratev_A· ∇P_CR/P_CR∼v_A/L₀, due to misalignment between the mean field and isotropic CR gradients. Unlike acceleration efficiencies, CR losses are almost independent of magnetic field strength overβ∼ 1–100 and are, therefore, not the primary factor behind lower acceleration rates when streaming is included. While this paper is primarily concerned with how turbulence affects CRs, in a follow-up paper we consider how CRs affect turbulence by diverting energy from the MHD cascade, altering the pathway to gas heating and steepening the turbulent spectrum. 

   more » « less
5. Cosmic-Ray Drag and Damping of Compressive Turbulence

   https://doi.org/10.3847/1538-4357/aceef9  

   Bustard, Chad; Oh, S. Peng (September 2023, The Astrophysical Journal)

   Abstract While it is well known that cosmic rays (CRs) can gain energy from turbulence via second-order Fermi acceleration, how this energy transfer affects the turbulent cascade remains largely unexplored. Here, we show that damping and steepening of the compressive turbulent power spectrum are expected once the damping time $t_{\mathrm{damp}}\sim \rho v^2/{\dot{E}}_{\mathrm{CR}}\propto E_{\mathrm{CR}}^{-1}$becomes comparable to the turbulent cascade time. Magnetohydrodynamic simulations of stirred compressive turbulence in a gas-CR fluid with diffusive CR transport show clear imprints of CR-induced damping, saturating at ${\dot{E}}_{\mathrm{CR}}\sim \widetilde{ϵ}$, where $\widetilde{ϵ}$is the turbulent energy input rate. In that case, almost all of the energy in large-scale motions is absorbed by CRs and does not cascade down to grid scale. Through a Hodge–Helmholtz decomposition, we confirm that purely compressive forcing can generate significant solenoidal motions, and we find preferential CR damping of the compressive component in simulations with diffusion and streaming, rendering small-scale turbulence largely solenoidal, with implications for thermal instability and proposed resonant scattering ofE≳ 300 GeV CRs by fast modes. When CR transport is streaming dominated, CRs also damp large-scale motions, with kinetic energy reduced by up to 1 order of magnitude in realisticE_CR∼E_gscenarios, but turbulence (with a reduced amplitude) still cascades down to small scales with the same power spectrum. Such large-scale damping implies that turbulent velocities obtained from the observed velocity dispersion may significantly underestimate turbulent forcing rates, i.e., $\widetilde{ϵ}\gg \rho v^3/L$. 

   more » « less