We use kinetic hybrid simulations (kinetic ions-fluid electrons) to characterize the fraction of ions that are accelerated to non-thermal energies at non-relativistic collisionless shocks. We ...investigate the properties of the shock discontinuity and show that shocks propagating almost along the background magnetic field (quasi-parallel shocks) reform quasi-periodically on ion cyclotron scales. Ions that impinge on the shock when the discontinuity is the steepest are specularly reflected. This is a necessary condition for being injected, but it is not sufficient. Also, by following the trajectories of reflected ions, we calculate the minimum energy needed for injection into diffusive shock acceleration, as a function of the shock inclination. We construct a minimal model that accounts for the ion reflection from quasi-periodic shock barrier, for the fraction of injected ions, and for the ion spectrum throughout the transition from thermal to non-thermal energies. This model captures the physics relevant for ion injection at non-relativistic astrophysical shocks with arbitrary strengths and magnetic inclinations, and represents a crucial ingredient for understanding the diffusive shock acceleration of cosmic rays.
We propose that the inner engine of a type I binary-driven hypernova (BdHN) is composed of Kerr black hole (BH) in a non-stationary state, embedded in a uniform magnetic field B0 aligned with the BH ...rotation axis and surrounded by an ionized plasma of extremely low density of 10−14 g cm−3. Using GRB 130427A as a prototype, we show that this inner engine acts in a sequence of elementary impulses. Electrons accelerate to ultrarelativistic energy near the BH horizon, propagating along the polar axis, θ = 0, where they can reach energies of ∼1018 eV, partially contributing to ultrahigh-energy cosmic rays. When propagating with through the magnetic field B0, they produce GeV and TeV radiation through synchroton emission. The mass of BH, M = 2.31M , its spin, = 0.47, and the value of magnetic field B0 = 3.48 × 1010 G, are determined self consistently to fulfill the energetic and the transparency requirement. The repetition time of each elementary impulse of energy erg is ∼10−14 s at the beginning of the process, then slowly increases with time evolution. In principle, this "inner engine" can operate in a gamma-ray burst (GRB) for thousands of years. By scaling the BH mass and the magnetic field, the same inner engine can describe active galactic nuclei.
Abstract
We present an efficient algorithm to follow spectral evolution of cosmic rays (CR) coupled with an MHD system on Eulerian grids. The algorithm is designed for studies of CR energy spectrum ...evolution in MHD simulations of a galactic interstellar medium. The base algorithm for CR transport relies on the two-moment piece-wise power-law method, known also as coarse-grained momentum finite volume (CGMV), for solving the Fokker–Planck CR transport equation, with a low number of momentum bins extending over several decades of the momentum coordinate. We propose an extension of the CGMV with a novel feature that allows momentum boundaries to change in response to CR momentum gains or losses near the extremes of the population distribution. Our extension involves a special treatment of momentum bins containing spectral cutoff. Contrary to the regular bins of fixed width, those bins have variable width, and their outer edges coincide with spectral cutoffs. The cutoff positions are estimated from the particle number density and energy density in the outer bins for an assumed small value of an additional parameter representing the smallest physically significant level of CR spectral energy density. We performed a series of elementary tests to validate the algorithm and demonstrated, whenever possible, that results of the test simulations correspond, with a reasonable accuracy, to the results of analogous analytical solutions. In a more complex test of the galactic CR-driven wind problem, we obtained results consistent with expectations regarding the effects of advection, diffusion, adiabatic, and synchrotron cooling of a CR population.
Abstract
The past decade has brought impressive advances in the astrophysics of cosmic rays (CRs) and multiwavelength astronomy, thanks to the new instrumentation launched into space and built on the ...ground. Modern technologies employed by those instruments provide measurements with unmatched precision, enabling searches for subtle signatures of dark matter and new physics. Understanding the astrophysical backgrounds to better precision than the observed data is vital in moving to this new territory. A state-of-the-art CR propagation code, called
GalProp,
is designed to address exactly this challenge. Having 25 yr of development behind it, the
GalProp
framework has become a de facto standard in the astrophysics of CRs, diffuse photon emissions (radio to
γ
-rays), and searches for new physics.
GalProp
uses information from astronomy, particle physics, and nuclear physics to predict CRs and their associated emissions self-consistently, providing a unifying modeling framework. The range of its physical validity covers 18 orders of magnitude in energy, from sub-keV to PeV energies for particles and from
μ
eV to PeV energies for photons. The framework and the data sets are public and are extensively used by many experimental collaborations and by thousands of individual researchers worldwide for interpretation of their data and for making predictions. This paper details the latest release of the
GalProp
framework and updated cross sections, further developments of its initially auxiliary data sets for models of the interstellar medium that grew into independent studies of the Galactic structure—distributions of gas, dust, radiation, and magnetic fields—as well as the extension of its modeling capabilities. Example applications included with the distribution illustrating usage of the new features are also described.
Recent observations by PAMELA, Fermi-LAT, and AMS-02 have conclusively indicated a rise in the cosmic-ray positron fraction above 10 GeV, a feature which is impossible to mimic under the paradigm of ...secondary positron production with self-consistent Galactic cosmic-ray propagation models. A leading explanation for the positron fraction rise is an additional source of electron-positron pairs, for example one or more mature, energetic, and relatively nearby pulsars. We point out that any one of two well-known nearby pulsars, Geminga and Monogem, can satisfactorily provide enough positrons to reproduce AMS-02 observations. A smoking-gun signature of this scenario is an anisotropy in the arrival direction of the cosmic-ray electrons and positrons, which may be detectable by existing, or future, telescopes. The predicted anisotropy level is, at present, consistent with limits from Fermi-LAT and AMS-02. We argue that the large collecting area of atmospheric Cherenkov telescopes (ACTs) makes them optimal tools for detecting such an anisotropy. Specifically, we show that much of the proton and gamma -ray background which affects measurements of the cosmic-ray electron-position spectrum with ACTs may be controlled in the search for anisotiopies. We conclude that observations using archival ACT data could already constrain or substantiate the pulsar origin of the positron anomaly, while upcoming instruments (such as the Cherenkov Telescope Array) will provide strong constraints on the source of the rising position fraction.
Context. The AMS-02 collaboration recently released cosmic ray F/Si data with an unprecedented accuracy. Cosmic ray fluorine is predominantly produced by fragmentation of heavier progenitors, while ...silicon is mostly accelerated at source. This ratio is thus maximally sensitive to cosmic ray propagation. Aims. We study the compatibility of the transport parameters derived from the F/Si ratio with those obtained from the lighter Li/C, Be/C, and B/C ratios. We also inspect the cosmic ray source abundance of F, which is one of the few elements that has a high first ionisation potential but is only moderately volatile and is a potentially key element to study the acceleration mechanism of cosmic rays. Methods. We used the 1D diffusion model implemented in the USINE code and performed χ 2 analyses accounting for several systematic effects (energy correlations in data, nuclear cross sections, and solar modulation uncertainties). We also took advantage of the EXFOR nuclear database to update the F production cross sections for its most important progenitors (identified to be 56 Fe, 32 S, 28 Si, 27 Al, 24 Mg, 22 Ne, and 20 Ne). Results. The transport parameters obtained from AMS-02 F/Si data are compatible with those obtained from AMS-02 (Li,Be,B)/C data. The combined fit of all of these ratios leads to a χ min 2 /d.o.f. ≈ 1.1, with ≲10% adjustments of the B and F production cross sections (which are based on very few nuclear data points and would strongly benefit from new measurements). The F/Si ratio is compatible with a pure secondary origin of F, with a best-fit relative source abundance ( 19 F/ 28 Si) CRS ∼ 10 −3 and an upper limit of ∼5 × 10 −3 . Unfortunately, this limit is not sufficient to test global acceleration models of cosmic ray nuclei, for which values at the level of ∼10 −4 are required. Such levels could be attained with F/Si data with a few percent of accuracy at a few tens of TV, which is possibly within reach for the next generation of cosmic ray experiments.
Abstract We study sub-GeV dark matter (DM) particles that may annihilate or decay into Standard Model particles producing an exotic injection component in the Milky Way that leaves an imprint in both ...photon and cosmic-ray (CR) fluxes. Specifically, the DM particles may annihilate or decay into e + e − , μ + μ − , or π + π − and may radiate photons through their e ± products. The resulting e ± products can be directly observed in probes such as Voyager 1. Alternatively, the e ± products may produce bremsstrahlung radiation and upscatter the low-energy Galactic photon fields via the inverse Compton process, generating a broad emission from X-ray to γ -ray energies observable in experiments such as XMM-Newton. We find that we get a significant improvement in the DM annihilation and decay constraints from XMM-Newton (excluding thermally averaged cross sections of 10 −31 cm 3 s −1 ≲ 〈 σ v 〉 ≲ 10 −26 cm 3 s −1 and decay lifetimes of 10 26 s ≲ τ ≲ 10 28 s, respectively) by including best-fit CR propagation and diffusion parameters. This yields the strongest astrophysical constraints for this mass range of DM of 1 MeV to a few GeV and even surpasses cosmological bounds across a wide range of masses as well.
ABSTRACT We study the damping processes of both incompressible and compressible magnetohydrodynamic (MHD) turbulence in a partially ionized medium. We start from the linear analysis of MHD waves, ...applying both single-fluid and two-fluid treatments. The damping rates derived from the linear analysis are then used in determining the damping scales of MHD turbulence. The physical connection between the damping scale of MHD turbulence and the cutoff boundary of linear MHD waves is investigated. We find two branches of slow modes propagating in ions and neutrals, respectively, below the damping scale of slow MHD turbulence, and offer a thorough discussion of their propagation and dissipation behavior. Our analytical results are shown to be applicable in a variety of partially ionized interstellar medium (ISM) phases and the solar chromosphere. The importance of neutral viscosity in damping the Alfvenic turbulence in the interstellar warm neutral medium and the solar chromosphere is demonstrated. As a significant astrophysical utility, we introduce damping effects to the propagation of cosmic rays in partially ionized ISM. The important role of turbulence damping in both transit-time damping and gyroresonance is identified.