More Like this

1. Cosmic ray protons and electrons from supernova remnants

   https://doi.org/10.1051/0004-6361/202140448  

   Cristofari, P. ; Blasi, P. ; Caprioli, D. ( June 2021 , Astronomy & Astrophysics)

   null (Ed.)

   Context. The spectrum of cosmic ray protons and electrons released by supernova remnants throughout their evolution is poorly known because of the difficulty in accounting for particle escape and confinement downstream of a shock front, where both adiabatic and radiative losses are present. Since electrons lose energy mainly through synchrotron losses, it is natural to ask whether the spectrum released into the interstellar medium may be different from that of their hadronic counterpart. Independent studies of cosmic ray transport through the Galaxy require that the source spectrum of electrons and protons be very different. Therefore, the above question acquires a phenomenological relevance. Aims. Here we calculate the spectrum of cosmic ray protons released during the evolution of supernovae of different types, accounting for the escape from the upstream region and for adiabatic losses of particles advected downstream of the shock and liberated at later times. The same calculation is carried out for electrons, where in addition to adiabatic losses we take the radiative losses suffered behind the shock into account. These electrons are dominated by synchrotron losses in the magnetic field, which most likely is self-generated by cosmic rays accelerated at the shock. Methods. We use standard temporal evolution relations for supernova shocks expanding in different types of interstellar media together with an analytic description of particle acceleration and magnetic field amplification to determine the density and spectrum of cosmic ray particles. Their evolution in time is derived by numerically solving the equation describing advection with adiabatic and radiative losses for electrons and protons. The flux from particles continuously escaping the supernova remnants is also accounted for. Results. The magnetic field in the post-shock region is calculated by using an analytic treatment of the magnetic field amplification due to nonresonant and resonant streaming instability and their saturation. The resulting field is compared with the available set of observational results concerning the dependence of the magnetic field strength upon shock velocity. We find that when the field is the result of the growth of the cosmic-ray-driven nonresonant instability alone, the spectrum of electrons and protons released by a supernova remnant are indeed different; however, such a difference becomes appreciable only at energies ≳100−1000 GeV, while observations of the electron spectrum require such a difference to be present at energies as low as ∼10 GeV. An effect at such low energies requires substantial magnetic field amplification in the late stages of supernova remnant evolution (shock velocity ≪1000 km s −1 ); this may not be due to streaming instability but rather hydrodynamical processes. We comment on the feasibility of such conditions and speculate on the possibility that the difference in spectral shape between electrons and protons may reflect either some unknown acceleration effect or additional energy losses in cocoons around the sources. 

   more » « less
2. Cosmic-ray generated bubbles around their sources

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

   Schroer, B. ; Pezzi, O. ; Caprioli, D. ; Haggerty, C. C. ; Blasi, P. ( February 2022 , Monthly Notices of the Royal Astronomical Society)

   Cosmic rays (CRs) are thought to escape their sources streaming along the local magnetic field lines. We show that this phenomenon generally leads to the excitation of both resonant and non-resonant streaming instabilities. The self-generated magnetic fluctuations induce particle diffusion in extended regions around the source, so that CRs build up a large pressure gradient. By means of two-dimensional (2D) and three-dimensional (3D) hybrid particle-in-cell simulations, we show that such a pressure gradient excavates a cavity around the source and leads to the formation of a cosmic ray dominated bubble, inside which diffusivity is strongly suppressed. Based on the trends extracted from self-consistent simulations, we estimate that, in the absence of severe damping of the self-generated magnetic fields, the bubble should keep expanding until pressure balance with the surrounding medium is reached, corresponding to a radius of ∼10–50 pc. The implications of the formation of these regions of low diffusivity for sources of Galactic CRs are discussed. Special care is devoted to estimating the self-generated diffusion coefficient and the grammage that CRs might accumulate in the bubbles before moving into the interstellar medium. Based on the results of 3D simulations, general considerations on the morphology of the γ-ray and synchrotron emission from these extended regions also are outlined.

    

   more » « less
3. Numerical study of cosmic ray confinement through dust resonant drag instabilities

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

   Ji, Suoqing ; Squire, Jonathan ; Hopkins, Philip F ( April 2022 , Monthly Notices of the Royal Astronomical Society)

   ABSTRACT We investigate the possibility of cosmic ray (CR) confinement by charged dust grains through resonant drag instabilities (RDIs). We perform magnetohydrodynamic particle-in-cell simulations of magnetized gas mixed with charged dust and cosmic rays, with the gyro-radii of dust and GeV CRs on ∼au scales fully resolved. As a first study, we focus on one type of RDI wherein charged grains drift super-Alfvénically, with Lorentz forces strongly dominating over drag forces. Dust grains are unstable to the RDIs and form concentrated columns and sheets, whose scale grows until saturating at the simulation box size. Initially perfectly streaming CRs are strongly scattered by RDI-excited Alfvén waves, with the growth rate of the CR perpendicular velocity components equaling the growth rate of magnetic field perturbations. These rates are well-predicted by analytic linear theory. CRs finally become isotropized and drift at least at ∼vA by unidirectional Alfvén waves excited by the RDIs, with a uniform distribution of the pitch angle cosine μ and a flat profile of the CR pitch angle diffusion coefficient Dμμ around μ = 0, without the ‘90○ pitch angle problem.’ With CR feedback on the gas included, Dμμ decreases by a factor of a few, indicating a lower CR scattering rate, because the backreaction on the RDI from the CR pressure adds extra wave damping, leading to lower quasi-steady-state scattering rates. Our study demonstrates that the dust-induced CR confinement can be very important under certain conditions, e.g. the dusty circumgalactic medium around quasars or superluminous galaxies. 

   more » « less
4. Acceleration and Spectral Redistribution of Cosmic Rays in Radio-jet Shear Flows

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

   Webb, G. M. ; Xu, Y. ; Biermann, P. L. ; Al-Nussirat, S. ; Mostafavi, P. ; Li, G. ; Barghouty, A. F. ; Zank, G. P. ( November 2023 , The Astrophysical Journal)

   A steady-state, semi-analytical model of energetic particle acceleration in radio-jet shear flows due to cosmic-ray viscosity obtained by Webb et al. is generalized to take into account more general cosmic-ray boundary spectra. This involves solving a mixed Dirichlet–Von Neumann boundary value problem at the edge of the jet. The energetic particle distribution functionf₀(r,p) at cylindrical radiusrfrom the jet axis (assumed to lie along thez-axis) is given by convolving the particle momentum spectrum$f_0(\infty ,p\prime )$with the Green’s function$G(r,p;p\prime )$, which describes the monoenergetic spectrum solution in which$f_0\to \delta (p-p\prime )$asr→ ∞ . Previous work by Webb et al. studied only the Green’s function solution for$G(r,p;p\prime )$. In this paper, we explore for the first time, solutions for more general and realistic forms for$f_0(\infty ,p\prime )$. The flow velocityu=u(r)e_zis along the axis of the jet (thez-axis).uis independent ofz, andu(r) is a monotonic decreasing function ofr. The scattering time$\tau {(r,p)={\tau }_0(p/p_0)}^{\alpha }$in the shear flow region 0 <r<r₂, and$\tau {(r,p)={\tau }_0(p/p_0)}^{\alpha }{(r/r_2)}^s$, wheres> 0 in the regionr>r₂is outside the jet. Other original aspects of the analysis are (i) the use of cosmic ray flow lines in (r,p) space to clarify the particle spatial transport and momentum changes and (ii) the determination of the probability distribution${\psi }_p(r,p;p\prime )$that particles observed at (r,p) originated fromr→ ∞ with momentum$p\prime$. The acceleration of ultrahigh-energy cosmic rays in active galactic nuclei jet sources is discussed. Leaky box models for electron acceleration are described.

    

   more » « less
5. Core-collapse supernovae in dense environments – particle acceleration and non-thermal emission

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

   Brose, R. ; Sushch, I. ; Mackey, J. ( August 2022 , Monthly Notices of the Royal Astronomical Society)

   Supernova remnants (SNRs) are known to accelerate cosmic rays from the detection of non-thermal emission in radio waves, X-rays, and gamma-rays. However, the ability to accelerate cosmic rays up to PeV energies has yet to be demonstrated. The presence of cut-offs in the gamma-ray spectra of several young SNRs led to the idea that PeV energies might only be achieved during the first years of a remnant’s evolution. We use our time-dependent acceleration-code RATPaC to study the acceleration of cosmic rays in supernovae expanding into dense environments around massive stars. We performed spherically symmetric one-dimensional (1D) simulations in which we simultaneously solve the transport equations for cosmic rays, magnetic turbulence, and the hydrodynamical flow of the thermal plasma in the test-particle limit. We investigated typical circumstellar-medium (CSM) parameters expected around red supergiant (RSG) and luminous blue variable (LBV) stars for freely expanding winds and accounted for the strong γγ absorption in the first days after explosion. The maximum achievable particle energy is limited to below $600\,$TeV even for the largest considered values of the magnetic field and mass-loss rates. The maximum energy is not expected to surpass $\approx 200\,$ and $\approx 70\,$TeV for LBVs and RSGs that experience moderate mass-loss prior to the explosion. We find gamma-ray peak-luminosities consistent with current upper limits and evaluate that current-generation instruments are able to detect the gamma-rays from Type-IIP explosions at distances up to $\approx 60\,$ kpc and Type-IIn explosions up to $\approx 1.0\,$ Mpc. We also find a good agreement between the thermal X-ray and radio synchrotron emission predicted by our models with a range of observations.

    

   more » « less