skip to main content


Title: Superdiffusion of cosmic rays in compressible magnetized turbulence
ABSTRACT

Owing to the complexity of turbulent magnetic fields, modelling the diffusion of cosmic rays is challenging. Based on the current understanding of anisotropic magnetohydrodynamic (MHD) turbulence, we use test particles to examine the cosmic rays’ superdiffusion in the direction perpendicular to the mean magnetic field. By changing Alfvén Mach number MA and sonic Mach number MS of compressible MHD simulations, our study covers a wide range of astrophysical conditions including subsonic warm gas phase and supersonic cold molecular gas. We show that freely streaming cosmic rays’ perpendicular displacement increases as 3/2 to the power of the time travelled along local magnetic field lines. This power-law index changes to 3/4 if the parallel propagation is diffusive. We find that the cosmic rays’ parallel mean free path decreases in a power-law relation of $M_\mathrm{ A}^{-2}$ in supersonic turbulence. We investigate the energy fraction of slow, fast, and Alfvénic modes and confirm the dominance of Alfvénic modes in the perpendicular superdiffusion. In particular, the energy fraction of fast mode, which is the main agent for pitch-angle scattering, increases with MA, but is insensitive to MS ≥ 2. Accordingly, our results suggest that the suppressed diffusion in supersonic molecular clouds arises primarily due to the variations of MA instead of MS.

 
more » « less
NSF-PAR ID:
10364327
Author(s) / Creator(s):
; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
512
Issue:
2
ISSN:
0035-8711
Page Range / eLocation ID:
p. 2111-2124
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We investigate the validity of Taylor’s hypothesis (TH) in the analysis of velocity and magnetic field fluctuations in Alfvénic solar wind streams measured by Parker Solar Probe (PSP) during the first four encounters. The analysis is based on a recent model of the spacetime correlation of magnetohydrodynamic (MHD) turbulence, which has been validated in high-resolution numerical simulations of strong reduced MHD turbulence. We use PSP velocity and magnetic field measurements from 24 h intervals selected from each of the first four encounters. The applicability of TH is investigated by measuring the parameter ϵ  =  δu 0 /√2 V ⊥ , which quantifies the ratio between the typical speed of large-scale fluctuations, δu 0 , and the local perpendicular PSP speed in the solar wind frame, V ⊥ . TH is expected to be applicable for ϵ ≲ 0.5 when PSP is moving nearly perpendicular to the local magnetic field in the plasma frame, irrespective of the Alfvén Mach number M A = V SW ∕ V A , where V SW and V A are the local solar wind and Alfvén speed, respectively. For the four selected solar wind intervals, we find that between 10 and 60% of the time, the parameter ϵ is below 0.2 and the sampling angle (between the spacecraft velocity in the plasma frame and the local magnetic field) is greater than 30°. For angles above 30°, the sampling direction is sufficiently oblique to allow one to reconstruct the reduced energy spectrum E ( k ⊥ ) of magnetic fluctuations from its measured frequency spectra. The spectral indices determined from power-law fits of the measured frequency spectrum accurately represent the spectral indices associated with the underlying spatial spectrum of turbulent fluctuations in the plasma frame. Aside from a frequency broadening due to large-scale sweeping that requires careful consideration, the spatial spectrum can be recovered to obtain the distribution of fluctuation’s energy across scales in the plasma frame. 
    more » « less
  2. We investigate three-dimensional (3-D) bow shocks in a highly collisional magnetized aluminium plasma, generated during the ablation phase of an exploding wire array on the MAGPIE facility (1.4 MA, 240 ns). Ablation of plasma from the wire array generates radially diverging, supersonic ( $M_S \sim 7$ ), super-Alfvénic ( $M_A > 1$ ) magnetized flows with frozen-in magnetic flux ( $R_M \gg 1$ ). These flows collide with an inductive probe placed in the flow, which serves both as the obstacle that generates the magnetized bow shock, and as a diagnostic of the advected magnetic field. Laser interferometry along two orthogonal lines of sight is used to measure the line-integrated electron density. A detached bow shock forms ahead of the probe, with a larger opening angle in the plane parallel to the magnetic field than in the plane normal to it. Since the resistive diffusion length of the plasma is comparable to the probe size, the magnetic field decouples from the ion fluid at the shock front and generates a hydrodynamic shock, whose structure is determined by the sonic Mach number, rather than the magnetosonic Mach number of the flow. The 3-D simulations performed using the resistive magnetohydrodynamic (MHD) code Gorgon confirm this picture, but under-predict the anisotropy observed in the shape of the experimental bow shock, suggesting that non-MHD mechanisms may be important for modifying the shock structure. 
    more » « less
  3. Abstract During its 10th orbit around the Sun, the Parker Solar Probe sampled two intervals where the local Alfvén speed exceeded the solar wind speed, lasting more than 10 hours in total. In this paper, we analyze the turbulence and wave properties during these periods. The turbulence is observed to be Alfvénic and unbalanced, dominated by outward-propagating modes. The power spectrum of the outward-propagating Elsässer z + mode steepens at high frequencies while that of the inward-propagating z − mode flattens. The observed Elsässer spectra can be explained by the nearly incompressible (NI) MHD turbulence model with both 2D and Alfvénic components. The modeling results show that the z + spectra are dominated by the NI/slab component, and the 2D component mainly affects the z − spectra at low frequencies. An MHD wave decomposition based on an isothermal closure suggests that outward-propagating Alfvén and fast magnetosonic wave modes are prevalent in the two sub-Alfvénic intervals, while the slow magnetosonic modes dominate the super-Alfvénic interval in between. The slow modes occur where the wavevector is nearly perpendicular to the local mean magnetic field, corresponding to nonpropagating pressure-balanced structures. The alternating forward and backward slow modes may also be features of magnetic reconnection in the near-Sun heliospheric current sheet. 
    more » « less
  4. Abstract

    Subsonic, compressive turbulence transfers energy to cosmic rays (CRs), a process known as nonresonant reacceleration. It is often invoked to explain the observed ratios of primary to secondary CRs at ∼GeV energies, assuming wholly diffusive CR transport. However, such estimates ignore the impact of CR self-confinement and streaming. We study these issues in stirring box magnetohydrodynamic (MHD) simulations using Athena++, with field-aligned diffusive and streaming CR transport. For diffusion only, we find CR reacceleration rates in good agreement with analytic predictions. When streaming is included, reacceleration rates depend on plasmaβ. Due to streaming-modified phase shifts between CR and gas variables, they are slower than canonical reacceleration rates in low-βenvironments like the interstellar medium but remain unchanged in high-βenvironments like the intracluster medium. We also quantify the streaming energy-loss rate in our simulations. For sub-Alfvénic turbulence, it is resolution dependent (hence unconverged in large-scale simulations) and heavily suppressed compared to the isotropic loss ratevA· ∇PCR/PCRvA/L0, due to misalignment between the mean field and isotropic CR gradients. Unlike acceleration efficiencies, CR losses are almost independent of magnetic field strength overβ∼ 1–100 and are, therefore, not the primary factor behind lower acceleration rates when streaming is included. While this paper is primarily concerned with how turbulence affects CRs, in a follow-up paper we consider how CRs affect turbulence by diverting energy from the MHD cascade, altering the pathway to gas heating and steepening the turbulent spectrum.

     
    more » « less
  5. Abstract

    We study the impact of compressibility on two-dimensional turbulent flows, such as those modeling astrophysical disks. We demonstrate that the direction of cascade undergoes continuous transition as the Mach numberMaincreases, from inverse atMa = 0, to direct atMa=. Thus, atMa1comparable amounts of energy flow from the pumping scale to large and small scales, in accord with previous data. For supersonic turbulence withMa1the cascade is direct, as in three dimensions, which results in multifractal density field. For that regime (Ma1) we derive a Kolmogorov-type law for potential forcing and obtain an explicit expression for the third order correlation tensor of the velocity. We further show that all third order structure functions are zero up to first order in the inertial range scales, which is in sharp contrast with incompressible turbulence where the third order structure function, that describes the energy flux associated with the energy cascade is non-zero. The properties of compressible turbulence have significant implications on the amplification of magnetic fields in conducting fluids. We thus demonstrate that imposing external magnetic field on compressible flows of conducting fluids allows to manipulate the flow producing possibly large changes even at small Mach numbers. Thus Zeldovich’s antidynamo theorem, by which atMa = 0 the magnetic field is zero in the steady state, must be used with caution. Real flows have finiteMaand, however small it is, for large enough values ofI, the magnetic flux through the disk, the magnetic field changes the flow appreciably, or rearranges it completely. This renders the limitMa → 0 singular for non-zero values ofI. Of particular interest is the effect of the density multifractality, atMa1which is relevant for astrophysical disks. We demonstrate that in that regime, in the presence of non-zeroIthe magnetic field energy is enhanced by a large factor as compared to its estimates based on the mean field. Finally, based on the insights described above, we propose a novel two-dimensional Burgers’ turbulence, whose three-dimensional counterpart is used for studies of the large-scale structure of the Universe, as a model for supersonic two-dimensional magnetohydrodynamic flows.

     
    more » « less