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)

   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. 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)

   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
3. Turbulent Reacceleration of Streaming Cosmic Rays

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

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

   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
4. Deep learning cosmic ray transport from density maps of simulated, turbulent gas

   https://doi.org/10.1088/2632-2153/ad262a  

   Bustard, Chad ; Wu, John ( February 2024 , Machine Learning: Science and Technology)

   The coarse-grained propagation of galactic cosmic rays (CRs) is traditionally constrained by phenomenological models of Milky Way CR propagation fit to a variety of direct and indirect observables; however, constraining the fine-grained transport of CRs along individual magnetic field lines—for instance, diffusive vs streaming transport models—is an unsolved challenge. Leveraging a recent training set of magnetohydrodynamic turbulent box simulations, with CRs spanning a range of transport parameters, we use convolutional neural networks (CNNs) trained solely on gas density maps to classify CR transport regimes. We find that even relatively simple CNNs can quite effectively classify density slices to corresponding CR transport parameters, distinguishing between streaming and diffusive transport, as well as magnitude of diffusivity, with class accuracies between 92% and 99%. As we show, the transport-dependent imprints that CRs leave on the gas are not all tied to the resulting density power spectra: classification accuracies are still high even when image spectra are flattened (85%–98% accuracy), highlighting CR transport-dependent changes to turbulent phase information. We interpret our results with saliency maps and image modifications, and we discuss physical insights and future applications.

    

   more » « less
5. JWST MIRI/Medium Resolution Spectrograph (MRS) Observations and Spectral Models of the Underluminous Type Ia Supernova 2022xkq

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

   DerKacy, J. M. ; Ashall, C. ; Hoeflich, P. ; Baron, E. ; Shahbandeh, M. ; Shappee, B. J. ; Andrews, J. ; Baade, D. ; Balangan, E. F. ; Bostroem, K. A. ; et al ( January 2024 , The Astrophysical Journal)

   We present a JWST mid-infrared (MIR) spectrum of the underluminous Type Ia Supernova (SN Ia) 2022xkq, obtained with the medium-resolution spectrometer on the Mid-Infrared Instrument (MIRI) ∼130 days post-explosion. We identify the first MIR lines beyond 14μm in SN Ia observations. We find features unique to underluminous SNe Ia, including the following: isolated emission of stable Ni, strong blends of [Tiii], and large ratios of singly ionized to doubly ionized species in both [Ar] and [Co]. Comparisons to normal-luminosity SNe Ia spectra at similar phases show a tentative trend between the width of the [Coiii] 11.888μm feature and the SN light-curve shape. Using non-LTE-multi-dimensional radiation hydro simulations and the observed electron capture elements, we constrain the mass of the exploding WD. The best-fitting model shows that SN 2022xkq is consistent with an off-center delayed-detonation explosion of a near-Chandrasekhar mass WD ($M_{\mathrm{WD}}$≈1.37M_⊙) of high central density (ρ_c≥ 2.0 × 10⁹g cm⁻³) seen equator-on, which producedM(⁵⁶Ni) =0.324M_⊙andM(⁵⁸Ni) ≥0.06M_⊙. The observed line widths are consistent with the overall abundance distribution; and the narrow stable Ni lines indicate little to no mixing in the central regions, favoring central ignition of subsonic carbon burning followed by an off-center deflagration-to-detonation transition beginning at a single point. Additional observations may further constrain the physics revealing the presence of additional species including Cr and Mn. Our work demonstrates the power of using the full coverage of MIRI in combination with detailed modeling to elucidate the physics of SNe Ia at a level not previously possible.

    

   more » « less