No Arabic abstract
The transport of charged energetic particles in the presence of strong intermittent heliospheric turbulence is computationally analyzed based on known properties of the interplanetary magnetic field and solar wind plasma at 1 Astronomical Unit (AU). The turbulence is assumed to be static, composite, and quasi-three-dimensional with a varying energy distribution between a one-dimensional Alfvenic (slab) and a structured two-dimensional component. The spatial fluctuations of the turbulent magnetic field are modeled either as homogeneous with a Gaussian probability distribution function (PDF), or as intermittent on large and small scales with a q-Gaussian PDF. Simulations showed that energetic particle diffusion coefficients both parallel and perpendicular to the background magnetic field are significantly affected by intermittency in the turbulence. This effect is especially strong for parallel transport where for large-scale intermittency results show an extended phase of subdiffusive parallel transport during which cross-field transport diffusion dominates. The effects of intermittency are found to depend on particle rigidity and the fraction of slab energy in the turbulence, yielding a perpendicular to parallel mean free path ratio close to 1 for large-scale intermittency. Investigation of higher order transport moments (kurtosis) indicates that non-Gaussian statistical properties of the intermittent turbulent magnetic field are present in the parallel transport, especially for low rigidity particles at all times.
Cosmic ray transport on galactic scales depends on the detailed properties of the magnetized, multiphase interstellar medium (ISM). In this work, we post-process a high-resolution TIGRESS magnetohydrodynamic simulation modeling a local galactic disk patch with a two-moment fluid algorithm for cosmic ray transport. We consider a variety of prescriptions for the cosmic rays, from a simple purely diffusive formalism with constant scattering coefficient, to a physically-motivated model in which the scattering coefficient is set by critical balance between streaming-driven Alfven wave excitation and damping mediated by local gas properties. We separately focus on cosmic rays with kinetic energies of $sim 1$ GeV (high-energy) and $sim 30$~MeV (low-energy), respectively important for ISM dynamics and chemistry. We find that simultaneously accounting for advection, streaming, and diffusion of cosmic rays is crucial for properly modeling their transport. Advection dominates in the high-velocity, low-density, hot phase, while diffusion and streaming are more important in higher density, cooler phases. Our physically-motivated model shows that there is no single diffusivity for cosmic-ray transport: the scattering coefficient varies by four or more orders of magnitude, maximal at density $n_mathrm{H} sim 0.01, mathrm{cm}^{-3}$. Ion-neutral damping of Alfven waves results in strong diffusion and nearly uniform cosmic ray pressure within most of the mass of the ISM. However, cosmic rays are trapped near the disk midplane by the higher scattering rate in the surrounding lower-density, higher-ionization gas. The transport of high-energy cosmic rays differs from that of low-energy cosmic rays, with less effective diffusion and greater energy losses for the latter.
A dramatic increase in the accuracy and statistics of space-borne cosmic ray (CR) measurements has yielded several breakthroughs over the last several years. The most puzzling is the rise in the positron fraction above ~10 GeV over the predictions of the propagation models assuming pure secondary production. The accuracy of the antiproton production cross section is critical for astrophysical applications and searches for new physics since antiprotons in CRs seem to hold the keys to many puzzles including the origin of those excess positrons. However, model calculations of antiproton production in CR interactions with interstellar gas are often employing parameterizations that are out of date or are using outdated physical concepts. That may lead to an incorrect interpretation of antiproton data which could have broad consequences for other areas of astrophysics. In this work, we calculate antiproton production in pp-, pA-, and AA-interactions using EPOS-LHC and QGSJET-II-04, two of the most advanced Monte Carlo (MC) generators tuned to numerous accelerator data including those from the Large Hadron Collider (LHC). We show that the antiproton yields obtained with these MC generators differ by up to an order of magnitude from yields of parameterizations commonly used in astrophysics.
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 in the downstream of a shock front, where both adiabatic and radiative losses are present. Here we calculate the spectrum of cosmic ray protons released during the evolution of supernovae of different types, accounting for the escape from upstream and for adiabatic losses of particles advected downstream of the shock and liberated at later times. The same calculation is carried out for electrons. The magnetic field in the post-shock region is calculated by using an analytic treatment of the magnetic field amplification due to non--resonant and resonant streaming instability and their saturation. We find that when the field is the result of the growth of the cosmic-ray--driven non--resonant instability alone, the spectrum of electrons and protons released by a supernova remnant are indeed different, but such a difference becomes appreciable only at energies $gtrsim 100-1000$ GeV, while observations of the electron spectrum require such a difference to be present at energies as low as $sim 10$ GeV. An effect at such low energies requires substantial magnetic field amplification in the late stages of the supernova remnant evolution (shock velocity $ll 1000$ km/s), perhaps not due to streaming instability but 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.
The status of the Greisen-Zatsepin-Kuzmin (GZK) cutoff and pair-production dip in Ultra High Energy Cosmic Rays (UHECR) is discussed.They are the features in the spectrum of protons propagating through CMB radiation in extragalactic space, and discovery of these features implies that primary particles are mostly extragalactic protons. The spectra measured by AGASA, Yakutsk, HiRes and Auger detectors are in good agreement with the pair-production dip, and HiRes data have strong evidences for the GZK cutoff. The Auger spectrum,as presented at the 30th ICRC 2007, agrees with the GZK cutoff, too. The AGASA data agree well with the beginning of the GZK cutoff at E leq 80 EeV, but show the excess of events at higher energies, the origin of which is not understood. The difference in the absolute fluxes measured by different detectors disappears after energy shift within the systematic errors of each experiment.
Cosmic-ray transport in astrophysical environments is often dominated by the diffusion of particles in a magnetic field composed of both a turbulent and a mean component. This process needs to be understood in order to properly model cosmic-ray signatures. One of the most important aspects in the modeling of cosmic-ray diffusion is that fully resonant scattering, the most effective such process, is only possible if the wave spectrum covers the entire range of propagation angles. By taking the wave spectrum boundaries into account, we quantify cosmic-ray diffusion parallel and perpendicular to the guide field direction at turbulence levels above 5% of the total magnetic field. We apply our results of the parallel and perpendicular diffusion coefficient to the Milky Way. We show that simple purely diffusive transport is in conflict with observations of the inner Galaxy, but that just by taking a Galactic wind into account, data can be matched in the central 5 kpc zone. Further comparison shows that the outer Galaxy at $>5$ kpc, on the other hand, should be dominated by perpendicular diffusion, likely changing to parallel diffusion at the outermost radii of the Milky Way.