skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Anomalous Transport of Cosmic Rays in a Nonlinear Diffusion Model

Abstract

We investigate analytically and numerically the transport of cosmic rays following their escape from a shock or another localized acceleration site. Observed cosmic-ray distributions in the vicinity of heliospheric and astrophysical shocks imply that anomalous, superdiffusive transport plays a role in the evolution of the energetic particles. Several authors have quantitatively described the anomalous diffusion scalings, implied by the data, by solutions of a formal transport equation with fractional derivatives. Yet the physical basis of the fractional diffusion model remains uncertain. We explore an alternative model of the cosmic-ray transport: a nonlinear diffusion equation that follows from a self-consistent treatment of the resonantly interacting cosmic-ray particles and their self-generated turbulence. The nonlinear model naturally leads to superdiffusive scalings. In the presence of convection, the model yields a power-law dependence of the particle density on the distance upstream of the shock. Although the results do not refute the use of a fractional advection–diffusion equation, they indicate a viable alternative to explain the anomalous diffusion scalings of cosmic-ray particles.

Authors:
 [1]; ;  [2]
  1. Department of Mathematics, University of Waikato, P. B. 3105, Hamilton 3240 (New Zealand)
  2. Institut für Theoretische Physik IV, Ruhr-Universität Bochum, Universitätsstrasse 150, D-44780 Bochum (Germany)
Publication Date:
OSTI Identifier:
22663570
Resource Type:
Journal Article
Resource Relation:
Journal Name: Astrophysical Journal; Journal Volume: 841; Journal Issue: 1; Other Information: Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
79 ASTROPHYSICS, COSMOLOGY AND ASTRONOMY; ACCELERATION; ASTROPHYSICS; COMPUTERIZED SIMULATION; CONVECTION; COSMIC RADIATION; DENSITY; DIFFUSION; DIFFUSION EQUATIONS; DISTRIBUTION; EVOLUTION; NONLINEAR PROBLEMS; TRANSPORT THEORY; TURBULENCE

Citation Formats

Litvinenko, Yuri E., Fichtner, Horst, and Walter, Dominik. Anomalous Transport of Cosmic Rays in a Nonlinear Diffusion Model. United States: N. p., 2017. Web. doi:10.3847/1538-4357/AA71BA.
Litvinenko, Yuri E., Fichtner, Horst, & Walter, Dominik. Anomalous Transport of Cosmic Rays in a Nonlinear Diffusion Model. United States. doi:10.3847/1538-4357/AA71BA.
Litvinenko, Yuri E., Fichtner, Horst, and Walter, Dominik. Sat . "Anomalous Transport of Cosmic Rays in a Nonlinear Diffusion Model". United States. doi:10.3847/1538-4357/AA71BA.
@article{osti_22663570,
title = {Anomalous Transport of Cosmic Rays in a Nonlinear Diffusion Model},
author = {Litvinenko, Yuri E. and Fichtner, Horst and Walter, Dominik},
abstractNote = {We investigate analytically and numerically the transport of cosmic rays following their escape from a shock or another localized acceleration site. Observed cosmic-ray distributions in the vicinity of heliospheric and astrophysical shocks imply that anomalous, superdiffusive transport plays a role in the evolution of the energetic particles. Several authors have quantitatively described the anomalous diffusion scalings, implied by the data, by solutions of a formal transport equation with fractional derivatives. Yet the physical basis of the fractional diffusion model remains uncertain. We explore an alternative model of the cosmic-ray transport: a nonlinear diffusion equation that follows from a self-consistent treatment of the resonantly interacting cosmic-ray particles and their self-generated turbulence. The nonlinear model naturally leads to superdiffusive scalings. In the presence of convection, the model yields a power-law dependence of the particle density on the distance upstream of the shock. Although the results do not refute the use of a fractional advection–diffusion equation, they indicate a viable alternative to explain the anomalous diffusion scalings of cosmic-ray particles.},
doi = {10.3847/1538-4357/AA71BA},
journal = {Astrophysical Journal},
number = 1,
volume = 841,
place = {United States},
year = {Sat May 20 00:00:00 EDT 2017},
month = {Sat May 20 00:00:00 EDT 2017}
}
  • The original nonlinear guiding center (NLGC) theory by Matthaeus et al. was a breakthrough in establishing a theory that promised to reproduce for the first time the anomalous perpendicular diffusion results from test particle trajectory calculations in prescribed static magnetic field turbulence dominated by a two-dimensional component. The assumptions used in this approach guaranteed anomalous diffusion to be a classical process (the variance ({Delta}x {sup 2}) {proportional_to} ({Delta}t){sup {alpha}} with {alpha} = 1). However, Shalchi and Kourakis showed that similar calculations can be even better reproduced within the context of a generalized compound diffusion model for anomalous perpendicular diffusion wherebymore » anomalous diffusion is nonclassical (0 < {alpha} < 1 (subdiffusion) or 1 < {alpha} < 2 (superdiffusion)). In this paper, it is shown how NLGC theory can be generalized to model compound diffusion conditions consistent with the generalized compound diffusion model in terms of a nonuniform plasma medium. Such a medium is assumed to generate a distribution of cosmic-ray pitch-angle scattering times for cosmic rays interacting resonantly with a minor small-scale slab turbulence component. This suggests that the anomalous perpendicular diffusion from the test particle simulations might be explained best within the framework of anomalous diffusion in a nonuniform plasma medium. It is argued that, during intermediate times when the magnetic field appears to be static in the quiet solar wind near Earth, generalized NLGC theory predicts possibly subdiffusive anomalous perpendicular transport for intermediate-energy (E << 400-500 MeV) cosmic rays because they experience nonuniform scattering conditions. These conditions are suggested to be a product of stochastic wave growth of small-scale slab turbulence resulting intermittently in large patches of intense slab turbulence where scattering times are reduced and cosmic rays are trapped. Classical anomalous diffusion is expected to be restored at high cosmic-ray energies (E >> 400-500 MeV) in accordance with original NLGC theory because they encounter more uniform scattering conditions due to weaker stochastic wave growth. Magnetic turbulence observations near Earth during 2003 were used to show that the solar wind can be characterized as a nonuniform medium sustaining a power-law distribution of cosmic-ray pitch-angle scattering times. However, application of generalized NLGC theory to the 2003 observations produced unexpectedly a negative exponent for the variance of the anomalous diffusion ({alpha} < 0). This indicates that particle bunching occurs as cosmic rays are trapped in plasma regions associated with strong particle scattering. Such effects can plausibly be explained by the intermittent presence of strong compressive turbulence downstream of shocks, stream, and interaction regions.« less
  • We present the first results from a series of computer simulations of the solar modulation of galactic cosmic rays using a full three-dimensional model which incorporates all known important effects on particle transport, particle drifts, convection with the solar wind, energy loss, and anisotropic diffusion. The model is time-independent in the coordinate frame rotating with the Sun, so corotating effects can be studied. We consider modulation in an interplanetary magnetic field model in which the current sheet separating the northerm and southern solar hemispheres is warped and corotating with the Sun. The amplitude of the warp is varied to simulatemore » possible solar cycle variation of the magnetic field. We find substantial effects due to the warp of the current sheet. Comparison of the model results with various data is presented. Among other things, we show that the intensity may decrease away from the current sheet for both signs of the magnetic field, as suggested by recent observations, and in contrast with inferences from earlier, more approximate calculations.« less
  • We present a simple continuous-time random-walk model for the transport of energetic particles accelerated by a collection of supernova explosions in a galactic superbubble, developed to simulate and highlight signatures of anomalous transport on the particles' evolution and their spectra in a multi-shock context. We assume standard diffusive shock acceleration (DSA) theory for each shock encounter. The superbubble (an OB stars association) is idealized as a heterogeneous region of particle sources and sinks bounded by a random surface. The model is based on two coupled stochastic differential equations and is applied for protons and alpha particles. Using characteristic values formore » a typical bubble, our simulations suggest that acceleration and transport in the bubble may be sub-diffusive. In addition, a spectral break in the particles' evolution and spectra is evident located at Almost-Equal-To 10{sup 15} eV for protons and Almost-Equal-To 3 Multiplication-Sign 10{sup 15} eV for alphas. Our simulations are consistent with a bubble's mean magnetic field strength of Almost-Equal-To 1 {mu}G and a shock separation distance {approx}0.1 Multiplication-Sign the characteristic radius of the bubble. The simulations imply that the diffusion coefficient (for the elementary shock acceleration process) is {approx}< 10{sup 27} cm{sup 2} s{sup -1} at 1 GeV/c. While the sub-diffusive transport is readily attributed to the stochastic nature of the acceleration time according to DSA theory, the spectral break appears to be an artifact of transport in a finite medium. These simulations point to a new and intriguing phenomenon associated with the statistical nature of collective acceleration of high-energy cosmic rays in galactic superbubbles.« less
  • Using current modulation theory for the transport of low-energy cosmic ray particles and the model of Fisk (1976a) for their acceleration by transit time damping large-scale field variations in the outer solar system, we construct a set of parameters which reproduce the quiet time spectra of He, O, N, and Ne as observed at 1 AU during the time period 1973--1975. With an analytic approximation for the acceleration and a numerical solution for the steady state spherical symmetric transport equation, both the observed spectral shapes and the relative intensities for He, O, N, and Ne can be fitted simultaneously remarkablymore » well. For the fit, one set of parameters is used: the acceleration rate, the size of the heliosphere, and the total number of particles injected. To obtain this number, it is assumed that independent of the species, the same fraction of the particles, which are ionized inside the heliosphere (by UV radiation and charge exchange with solar wind particles), participate in the acceleration. For the calculation of the production rate of singly ionized ions the ionization rates provided by Siscoe and Mukherjee (1972) are used, and the solar system abundances of A. G. W. Cameron (unpublished manuscript, 1973) are assumed for the composition of the neutral wind.« less