Search for: All records

Abstract The end of supernova remnant (SNR) evolution is characterized by a so-called “radiative” stage, in which efficient cooling of the hot bubble inside the forward shock slows expansion, leading to eventual shock breakup. Understanding SNR evolution at this stage is vital for predicting feedback in galaxies, since SNRs are expected to deposit their energy and momentum into the interstellar medium at the ends of their lives. A key prediction of SNR evolutionary models is the formation at the onset of the radiative stage of a cold, dense shell behind the forward shock. However, searches for these shells via their neutral hydrogen emission have had limited success. We instead introduce an independent observational signal of shell formation arising from the interaction between nonthermal particles accelerated by the SNR forward shock (cosmic rays) and the dense shell. Using a semi-analytic model of particle acceleration based on state-of-the-art simulations coupled with a high-resolution hydrodynamic model of SNR evolution, we predict the nonthermal emission that arises from this interaction. We demonstrate that the onset of the radiative stage leads to nonthermal signatures from radio to gamma rays, including radio and gamma-ray brightening by nearly 2 orders of magnitude. Such a signature may be detectable with current instruments, and will be resolvable with the next generation of gamma-ray telescopes (namely, the Cherenkov Telescope Array). 

Abstract The nonresonant streaming instability (Bell instability) plays a pivotal role in the acceleration and confinement of cosmic rays (CRs), yet the exact mechanism responsible for its saturation and the magnitude of the final amplified magnetic field have not been assessed from first principles. Using a survey of hybrid simulations (with kinetic ions and fluid electrons), we study the evolution of the Bell instability as a function of the parameters of the CR population. We find that at saturation, the magnetic pressure in the amplified field is comparable with the initial CR anisotropic pressure, rather than with the CR energy flux, as previously argued. These results provide a predictive prescription for the total magnetic field amplification expected in the many astrophysical environments where the Bell instability is important. 

Abstract Collisionless shocks tend to send charged particles into the upstream, driving electric currents through the plasma. Using kinetic particle-in-cell simulations, we investigate how the background thermal plasma neutralizes such currents in the upstream of quasi-parallel non-relativistic electron–proton shocks. We observe distinct processes in different regions: the far upstream, the shock precursor, and the shock foot. In the far upstream, the current is carried by nonthermal protons, which drive electrostatic modes and produce suprathermal electrons that move toward upstream infinity. Closer to the shock (in the precursor), both the current density and the momentum flux of the beam increase, which leads to electromagnetic streaming instabilities that contribute to the thermalization of suprathermal electrons. At the shock foot, these electrons are exposed to shock-reflected protons, resulting in a two-stream type instability. We analyze these processes and the resulting heating through particle tracking and controlled simulations. In particular, we show that the instability at the shock foot can make the effective thermal speed of electrons comparable to the drift speed of the reflected protons. These findings are important for understanding both the magnetic field amplification and the processes that may lead to the injection of suprathermal electrons into diffusive shock acceleration. 

Abstract We present a survey of 1D kinetic particle-in-cell simulations of quasi-parallel nonrelativistic shocks to identify the environments favorable for electron acceleration. We explore an unprecedented range of shock speedsv_sh≈ 0.067–0.267c, Alfvén Mach numbers ${\mathcal{M}}_{\mathrm{A}}=5–40$, sonic Mach numbers ${\mathcal{M}}_{\mathrm{s}}=5–160$, as well as the proton-to-electron mass ratiosm_i/m_e= 16–1836. We find that high Alfvén Mach number shocks can channel a large fraction of their kinetic energy into nonthermal particles, self-sustaining magnetic turbulence and acceleration to larger and larger energies. The fraction of injected particles is ≲0.5% for electrons and ≈1% for protons, and the corresponding energy efficiencies are ≲2% and ≈10%, respectively. The extent of the nonthermal tail is sensitive to the Alfvén Mach number; when ${\mathcal{M}}_{\mathrm{A}}\lesssim 10$, the nonthermal electron distribution exhibits minimal growth beyond the average momentum of the downstream thermal protons, independently of the proton-to-electron mass ratio. Acceleration is slow for shocks with low sonic Mach numbers, yet nonthermal electrons still achieve momenta exceeding the downstream thermal proton momentum when the shock Alfvén Mach number is large enough. We provide simulation-based parameterizations of the transition from thermal to nonthermal distribution in the downstream (found at a momentum around $p_{\mathrm{i},\mathrm{e}}/m_{\mathrm{i}}{v}_{\mathrm{sh}}\approx 3\sqrt{m_{\mathrm{i},\mathrm{e}}/m_{\mathrm{i}}}$), as well as the ratio of nonthermal electron to proton number density. The results are applicable to many different environments and are important for modeling shock-powered nonthermal radiation. 

Abstract GeV and TeV emission from the forward shocks of supernova remnants (SNRs) indicates that they are capable particle accelerators, making them promising sources of Galactic cosmic rays (CRs). However, it remains uncertain whether thisγ-ray emission arises primarily from the decay of neutral pions produced by very-high-energy hadrons, or from inverse-Compton and/or bremsstrahlung emission from relativistic leptons. By applying a semi-analytic approach to non-linear diffusive shock acceleration, and calculating the particle and photon spectra produced in different environments, we parameterize the relative strength of hadronic and leptonic emission. We show that even if CR acceleration is likely to occur in all SNRs, the observed photon spectra may primarily reflect the environment surrounding the SNR: the emission is expected to look hadronic unless the ambient density is particularly low (with proton number density ≲0.1 cm⁻³) or the photon background is enhanced with respect to average Galactic values (with radiation energy densityu_rad≳ 10 eV cm⁻³). We introduce a hadronicity parameter to characterize how hadronic or leptonic we expect a source to look based on its environment, which can be used to guide the interpretation of currentγ-ray observations and the detection of high-energy neutrinos from SNRs. 

Abstract In 2021 August, the Fermi Large Area Telescope, H.E.S.S., and MAGIC detected GeV and TeVγ-ray emission from an outburst of recurrent nova RS Ophiuchi. This detection represents the first very high-energyγ-rays observed from a nova, and it opens a new window to study particle acceleration. Both H.E.S.S. and MAGIC described the observedγ-rays as arising from a single, external shock. In this paper, we perform detailed, multi-zone modeling of RS Ophiuchi’s 2021 outburst, including a self-consistent prescription for particle acceleration and magnetic field amplification. We demonstrate that, contrary to previous work, a single shock cannot simultaneously explain RS Ophiuchi’s GeV and TeV emission, in particular the spectral shape and distinct light-curve peaks. Instead, we put forward a model involving multiple shocks that reproduces the observedγ-ray spectrum and temporal evolution. The simultaneous appearance of multiple distinct velocity components in the nova optical spectrum over the first several days of the outburst supports the presence of distinct shocks, which may arise either from the strong latitudinal dependence of the density of the external circumbinary medium (e.g., in the binary equatorial plane versus the poles) or due to internal collisions within the white dwarf ejecta (which power theγ-ray emission in classical novae). 

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