Dynamics of charged particles and perpendicular diffusion in turbulent magnetic field

A test particle code is employed to explore the dynamics of charged particles and perpendicular diffusion in turbulent magnetic field, where a three-dimensional (3D) isotropic turbulence model is used in this paper. The obtained perpendicular diffusion at different particle energies is compared with that of the nonlinear guiding center (NLGC) theory. It is found that the NLGC theory is consistent with test particle simulations when the particle energies are small. However, the difference between the NLGC theory and test particle simulations tends to increase when the particle energy is sufficiently large, and the threshold is related to the turbulence bend-over length. In the NLGC theory, the gyrocenter of a charged particle is assumed to follow the magnetic field line. Therefore, when the particle has sufficiently large energy, its gyroradius will be larger than the turbulence bend-over length. Then the particle can cross the magnetic field lines, and the difference between the test particle simulations and NLGC theory occurs.     

Saturated magnetic field amplification at supernova shocks

Cosmic ray streaming instabilities at supernova shocks are discussed in the quasi-linear diffusion formalism which takes into account the feedback effect of wave growth on the cosmic ray streaming motion. In particular, the non-resonant instability that leads to magnetic field amplification in the short wavelength regime is considered. The linear growth rate is calculated using kinetic theory for a streaming distribution. We show that the non-resonant instability is actually driven by a compensating current in the background plasma. The non-resonant instability can develop into a non-linear regime generating turbulence. The saturation of the amplified magnetic fields due to particle diffusion in the turbulence is derived analytically. It is shown that the evolution of parallel and perpendicular cosmic ray pressures is predominantly determined by non-resonant diffusion. However, the saturation is determined by resonant diffusion which tends to reduce the streaming motion through pitch angle scattering. The saturated level can exceed the mean background magnetic field.     

Non Gaussian and Non Local Transport in the Earth's Distant Magnetotail

We study the ion dynamics in a magnetic field reversal with a constant electric field and with a model of three dimensional magnetic turbulence. By computing the mean square displacements in the plane of the current sheet we find superdiffusive and superballistic transport regimes. Since velocity increases with the length of the free path, we have accelerated Lévyflights. The possibility to generate power law velocity distribution functions is pointed out, as well as the long memory effects and non local properties of ion transport. This revised version was published online in July 2006 with corrections to the Cover Date.     

Compound diffusion of energetic particles: a Kappa model for the parallel distribution function

One can assume that energetic particles follow magnetic field lines while they propagate through a magnetized plasma. The latter scenario is usually described by the so-called field line random walk limit. This limit, however, is only valid if parallel diffusion is suppressed. As soon as the latter effect is taken into account, perpendicular transport becomes subdiffusive. This physical scenario is usually called compound diffusion or compound subdiffusion and can be described by a Chapman-Kolmogorov equation. In the latter equation the parallel distribution function is an essential ingredient. In the present paper we replace the standard Gaussian model by a Kappa distribution to compute distribution functions and mean square displacements across the field.     

A study of the magnetic helicity during eight HELIOS-observed solar proton events

We have studied the influence of the magnetic helicity on solar particle propagation using the IMF data observed by the HELIOS spacecraft in the range 0.31–0.95 AU, during eight solar proton events. For this, we have derived power and helicity spectra of the turbulence of the magnetic field during the time of the events. These are used to compute the particle pitch-angle scattering coefficients according to the quasi-linear theory (QLT) treatment of particle propagation in turbulent magnetic fields. The results show that in all the cases the helicity effects are negligible and the particle's mean free paths deduced from the pitch-angle diffusion coefficients are the same regardless of whether or not helicity effects are included in the calculations. The computed mean free paths are quite different in each case.Deceased 10 April, 1995.     

A graphics-card implementation of Monte-Carlo simulations for cosmic-ray transport

A graphics card implementation of a test-particle simulation code is presented that is based on the CUDA extension of the C/C++ programming language. The original CPU version has been developed for the calculation of cosmic-ray diffusion coefficients in artificial Kolmogorov-type turbulence. In the new implementation, the magnetic turbulence generation, which is the most time-consuming part, is separated from the particle transport and is performed on a graphics card. In this article, the modification of the basic approach of integrating test particle trajectories to employ the SIMD (single instruction, multiple data) model is presented and verified. The efficiency of the new code is tested and several language-specific accelerating factors are discussed. For the example of isotropic magnetostatic turbulence, sample results are shown and a comparison to the results of the CPU implementation is performed.     

Non-linear magnetic diffusivity in mean-field electrodynamics

We consider non-linear transport and drift processes caused by an inhomogeneous magnetic field in a turbulent fluid. The coefficients of magnetic diffusivity and drift velocity are calculated by making use of the second-order correlation approximation. Transport processes in the presence of a sufficiently strong magnetic field become anisotropic with larger diffusion rate and turbulent electrical resistivity across the field than along the field. Non-linear effects also lead to a drift of the magnetic field away from the regions with a higher magnetic energy.     

Analytic forms of the perpendicular cosmic ray diffusion coefficient for an arbitrary turbulence spectrum and applications on transport of Galactic protons and acceleration at interplanetary shocks

We investigate cosmic ray scattering in the direction perpendicular to a mean magnetic field. Unlike in previous articles we employ a general form of the turbulence wave spectrum with arbitrary behavior in the energy range. By employing an improved version of the nonlinear guiding center theory we compute analytically the perpendicular mean free path. As shown, the energy range spectral index, has a strong influence on the perpendicular diffusion coefficient. If this parameter is larger than one we find for some cases a perpendicular diffusion coefficient that is independent of the parallel mean free path and particle energy. Two applications are considered, namely transport of Galactic protons in the solar system and diffusive particle acceleration at highly perpendicular interplanetary shock waves.     

Diffusive propagation of fast particles in the presence of a moving shock wave

Based on an analytical model, we determined the temporal dynamics of the spectral shape and spatial distribution of the particles that were impulsively (in time) injected with a specified spectrum in the vicinity of a moving plane shock front. We obtained a condition to determine the influence of the shock front on the particle propagation, where the spatial diffusion coefficient of the particles plays a major role. Diffusive shock acceleration is shown to strongly affect low-energy particles (the intensity maximum coincides spatially with the shock front; hard and soft spectral regions are formed in the spectrum) and weakly affect high-energy particles (the time at which the intensity reaches its maximum is well ahead of the shock arrival time; the spectral shape does not change). In events accompanied by a significant increase in the turbulence level, the influence of the shock front on high-energy particles can change from weak to strong. This change shows up in the spatial distribution and spectral shape of the particles. The dynamics of the particle intensity, calculated with the diffusion coefficients that were determined in accordance with the quasi-linear theory for measured turbulence levels, qualitatively corresponds to the observed solar energetic-particle intensity.     

Stability of accretion discs threaded by a strong magnetic field

We study the stability of poloidal magnetic fields anchored in a thin accretion disc. The two-dimensional hydrodynamics in the disc plane is followed by a grid-based numerical simulation including the vertically integrated magnetic forces. The three-dimensional magnetic field outside the disc is calculated in a potential field approximation from the magnetic flux density distribution in the disc. For uniformly rotating discs we confirm numerically the existence of the interchange instability as predicted by Spruit, Stehle & Papaloizou . In agreement with predictions from the shearing sheet model, discs with Keplerian rotation are found to be stabilized by the shear, as long as the contribution of magnetic forces to support against gravity is small. When this support becomes significant, we find a global instability which transports angular momentum outwardly and allows mass to accrete inwardly. The instability takes the form of a m =1 rotating 'crescent', reminiscent of the purely hydrodynamic non-linear instability previously found in pressure-supported discs. A model where the initial surface mass density Σ( r ) and B _z ( r ) decrease with radius as power laws shows transient mass accretion during about six orbital periods, and settles into a state with surface density and field strength decreasing approximately exponentially with radius. We argue that this instability is likely to be the main angular momentum transport mechanism in discs with a poloidal magnetic field sufficiently strong to suppress magnetic turbulence. It may be especially relevant in jet-producing discs.     

Theoretical interpretation of cosmic magnetic fields

The paper discusses the possibilty of interpreting the magnetic fields of astronomical bodies in the framework of a unified field theory.Using one of the solutions of the generalized field theory, a direct relation between the polar magnetic field, the angular velocity and the gravitational potential of the body considered, is obtained. The geometric model used for applications has spherical symmetry and is of the type (FIGI).The predictions of the theoretical formula, obtained from the model, are compared with available observational data, and with the empirical formula of Blackett. The theoretical formula gives a possible interpretation of a seed magnetic field which will develop and produce the large-scale magnetic field observed for celestial objects. The formula shows that the field is generated as a result of rotation of a massive object.     

Magnetic field evolution in simulations with Euler potentials

Using two- and three-dimensional hydromagnetic simulations for a range of different flows, including laminar and turbulent ones, it is shown that solutions expressing the field in terms of Euler potentials (EP) are in general incorrect if the EP are evolved with an artificial diffusion term. In three dimensions, standard methods using the magnetic vector potential are found to permit dynamo action when the EP give decaying solutions. With an imposed field, the EP method yields excessive power at small scales. This effect is more exaggerated in the dynamic case, suggesting an unrealistically reduced feedback from the Lorentz force. The EP approach agrees with standard methods only at early times when magnetic diffusivity did not have time to act. It is demonstrated that the usage of EP with even a small artificial magnetic diffusivity does not converge to a proper solution of hydromagnetic turbulence. The source of this disagreement is not connected with magnetic helicity or the three-dimensionality of the magnetic field, but is simply due to the fact that the non-linear representation of the magnetic field in terms of EP that depend on the same coordinates is incompatible with the linear diffusion operator in the induction equation.     

The different transport regimes of pitch-angle scattering of energetic particles

Recently an advanced nonlinear diffusion theory for particle transport across the mean magnetic field has been developed. The method used in the derivation of the latter theory is based on the cosmic ray Fokker-Planck equation. In the present article we use the same approach to describe pitch-angle scattering and parallel spatial diffusion nonlinearly. Furthermore, we derive the quasilinear transport theory, the weakly nonlinear theory as well as the Bohm limit as special cases from our more general approach.     

Magnetic Field Line Transport in Anisotropic Magnetic Turbulence: Anomalous, Quasilinear, and Percolative Regimes Versus the Kubo Number

The transport of plasma and of energetic particles because of magnetic turbulence is relevant to many space plasmas, ranging from the planetary magnetospheres to the solar corona and to the heliosphere. Various transport regimes for magnetic field lines can be obtained depending on the Kubo number. Here we show, by means of a numerical simulation, that the Kubo number also determines the level of chaos of the field lines. Weak chaos, closed magnetic surfaces, and anomalous transport regimes are obtained for R≪ 1; widespread chaos, destroyed magnetic surfaces, and quasilinear scaling of the diffusion coefficient for R ≳ 0.3; and global stochasticity as well as percolation scaling of the diffusion coefficient for R≫ 1. This revised version was published online in July 2006 with corrections to the Cover Date.     

The generation of a large-scale Galactic magnetic field by electric currents of energetic particles

We consider the generation of a magnetic field in the Galaxy by the electric currents excited by cosmic-ray particles in the disk and halo. We assume that the sources of relativistic particles are distributed continuously and uniformly in the Galactic disk, their total power is equal to the observed value, and the particles themselves undergo anisotropic diffusion in a homogeneous medium. We take into account the differential rotation of the Galactic disk but disregard the turbulence gyrotropy (the α effect). The strength of the generated magnetic field in our model is shown to strongly depend on the symmetry of the relativistic proton and thermal electron diffusion tensors, as well as on the relations between the tensor components. In particular, if the diffusion is isotropic, then no magnetic field is generated. For the independent tensor components estimated from observed parameters of the Galactic medium and with a simultaneous allowance made for the turbulent field dissipation processes, the mechanism under consideration can provide an observable magnetic-field strength of the order of several microgauss. This mechanism does not require any seed magnetic field, which leads us to suggest that relativistic particles can give an appreciable and, possibly, determining contribution to the formation of the global Galactic magnetic field. However, a final answer can be obtained only from a nonlinear self-consistent treatment, in which the symmetry and magnitude of the particle diffusion tensor components should be determined together with the calculation of the magnetic field.     

Strong cosmic-ray scattering in an anisotropic random magnetic field

We calculate the kinetic coefficients and the transport mean free paths of high-energy particles parallel to the regular magnetic field in the approximation of a large-scale anisotropic random magnetic field by using a nonlinear collision integral, i.e., by taking into account the processes of strong random scattering. We consider the diffusion of solar and Galactic cosmic rays by two-dimensional turbulence. Strong random scattering by two-dimensional turbulence is shown to reduce the parallel transport mean free path several fold. The momentum dependence of the parallel mean free path does not change, Λ_∥ ∝ p^2?v. In the case of strong random scattering by turbulence formed by several modes, the parallel transport mean free path is Λ_∥ ∝ p. We show that two-dimensional turbulence can make a major contribution to the parallel transport mean free paths of cosmic rays in the heliosphere and the interstellar medium.     

Analytical forms of correlation functions and length scales of astrophysical turbulence

In this article correlation-functions, correlation-scales, and field line diffusion coefficients of magnetic turbulence are discussed. For the slab turbulence model these parameters are calculated analytically for different forms of the wave-spectrum of the magnetic fluctuations. It is also demonstrated in this article that the definition of the correlation-length used previously yields infinity in some cases. Furthermore, correlation-functions and correlation-scales are calculated analytically for the well-established slab/2D composite model. It is demonstrated that the 2D correlation function decays much slower than the slab correlation function. These results are essential for understanding previous theoretical results of superdiffusive field line random walk and very recent spacecraft measurements of magnetic correlations in the solar wind.