No Arabic abstract
Highly energetic, relativistic electrons are commonly present in many astrophysical systems, from solar flares to the intra-cluster medium, as indicated by observed electromagnetic radiation. However, open questions remain about the mechanisms responsible for their acceleration, and possible re-acceleration. Ubiquitous plasma turbulence is one of the possible universal mechanisms. We study the energization of transrelativistic electrons in turbulence using hybrid particle-in-cell, which provide a realistic model of Alfv{e}nic turbulence from MHD to sub-ion scales, and test particle simulations for electrons. We find that, depending on the electron initial energy and turbulence strength, electrons may undergo a fast and efficient phase of energization due to the magnetic curvature drift during the time they are trapped in dynamic magnetic structures. In addition, electrons are accelerated stochastically which is a slower process that yields lower maximum energies. The combined effect of these two processes determines the overall electron acceleration. With appropriate turbulence parameters, we find that superthermal electrons can be accelerated up to relativistic energies. For example, with heliospheric parameters and a relatively high turbulence level, rapid energization to MeV energies is possible.
The partition of irreversible heating between ions and electrons in compressively driven (but subsonic) collisionless turbulence is investigated by means of nonlinear hybrid gyrokinetic simulations. We derive a prescription for the ion-to-electron heating ratio $Q_rmi/Q_rme$ as a function of the compressive-to-Alfvenic driving power ratio $P_compr/P_AW$, of the ratio of ion thermal pressure to magnetic pressure $beta_rmi$, and of the ratio of ion-to-electron background temperatures $T_rmi/T_rme$. It is shown that $Q_rmi/Q_rme$ is an increasing function of $P_compr/P_AW$. When the compressive driving is sufficiently large, $Q_rmi/Q_rme$ approaches $simeq P_compr/P_AW$. This indicates that, in turbulence with large compressive fluctuations, the partition of heating is decided at the injection scales, rather than at kinetic scales. Analysis of phase-space spectra shows that the energy transfer from inertial-range compressive fluctuations to sub-Larmor-scale kinetic Alfven waves is absent for both low and high $beta_rmi$, meaning that the compressive driving is directly connected to the ion entropy fluctuations, which are converted into ion thermal energy. This result suggests that preferential electron heating is a very special case requiring low $beta_rmi$ and no, or weak, compressive driving. Our heating prescription has wide-ranging applications, including to the solar wind and to hot accretion disks such as M87 and Sgr A*.
A scaling theory of long-wavelength electrostatic turbulence in a magnetised, weakly collisional plasma (e.g., ITG turbulence) is proposed, with account taken both of the nonlinear advection of the perturbed particle distribution by fluctuating ExB flows and of its phase mixing, which is caused by the streaming of the particles along the mean magnetic field and, in a linear problem, would lead to Landau damping. It is found that it is possible to construct a consistent theory in which very little free energy leaks into high velocity moments of the distribution function, rendering the turbulent cascade in the energetically relevant part of the wave-number space essentially fluid-like. The velocity-space spectra of free energy expressed in terms of Hermite-moment orders are steep power laws and so the free-energy content of the phase space does not diverge at infinitesimal collisionality (while it does for a linear problem); collisional heating due to long-wavelength perturbations vanishes in this limit (also in contrast with the linear problem, in which it occurs at the finite rate equal to the Landau-damping rate). The ability of the free energy to stay in the low velocity moments of the distribution function is facilitated by the anti-phase-mixing effect, whose presence in the nonlinear system is due to the stochastic version of the plasma echo (the advecting velocity couples the phase-mixing and anti-phase-mixing perturbations). The partitioning of the wave-number space between the (energetically dominant) region where this is the case and the region where linear phase mixing wins its competition with nonlinear advection is governed by the critical balance between linear and nonlinear timescales (which for high Hermite moments splits into two thresholds, one demarcating the wave-number region where phase mixing predominates, the other where plasma echo does).
A Hamiltonian two-field gyrofluid model for kinetic Alfven waves (KAWs) in a magnetized electron-proton plasma, retaining ion finite-Larmor-radius corrections and parallel magnetic field fluctuations, is used to study the inverse cascades that develop when turbulence is randomly driven at sub-ion scales. In the directions perpendicular to the ambient field, the dynamics of the cascade turns out to be nonlocal and the ratio $chi_f$ of the wave period to the characteristic nonlinear time at the driving scale affect some of its properties. For example, at small values of $chi_f$, parametric decay instability of the modes driven by the forcing can develop, enhancing for a while inverse transfers. The balanced state, obtained at early time when the two counter-propagating waves are equally driven, also becomes unstable at small $chi_f$, leading to an inverse cascade. For $beta_e$ smaller than a few units, the cascade slows down when reaching the low-dispersion spectral range. For higher $beta_e$, the ratio of the KAW to the Alfven frequencies displays a local minimum. At the corresponding transverse wavenumber, a condensate is formed, and the cascade towards larger scales is then inhibited. Depending on the parameters, a parallel inverse cascade can develop, enhancing the elongation of the ion-scale magnetic vortices that generically form.
Natures most powerful high-energy sources are capable of accelerating particles to high energy and radiate it away on extremely short timescales, even shorter than the light crossing time of the system. It is yet unclear what physical processes can produce such an efficient acceleration, despite the copious radiative losses. By means of radiative particle-in-cell simulations, we show that magnetically dominated turbulence in pair plasmas subject to strong synchrotron cooling generates a nonthermal particle spectrum with a hard power-law range (slope $p sim 1$) within a few eddy turnover times. Low pitch-angle particles can significantly exceed the nominal radiation-reaction limit, before abruptly cooling down. The particle spectrum becomes even harder ($p < 1$) over time owing to particle cooling with an energy-dependent pitch-angle anisotropy. The resulting synchrotron spectrum is hard ($ u F_ u propto u^s$ with $s sim 1$). Our findings have important implications for understanding the nonthermal emission from high-energy astrophysical sources, most notably the prompt phase of gamma-ray bursts and gamma-ray flares from the Crab nebula.
The current understanding of MHD turbulence envisions turbulent eddies which are anisotropic in all three directions. In the plane perpendicular to the local mean magnetic field, this implies that such eddies become current-sheet-like structures at small scales. We analyze the role of magnetic reconnection in these structures and conclude that reconnection becomes important at a scale $lambdasim L S_L^{-4/7}$, where $S_L$ is the outer-scale ($L$) Lundquist number and $lambda$ is the smallest of the field-perpendicular eddy dimensions. This scale is larger than the scale set by the resistive diffusion of eddies, therefore implying a fundamentally different route to energy dissipation than that predicted by the Kolmogorov-like phenomenology. In particular, our analysis predicts the existence of the sub-inertial, reconnection interval of MHD turbulence, with the Fourier energy spectrum $E(k_perp)propto k_perp^{-5/2}$, where $k_perp$ is the wave number perpendicular to the local mean magnetic field. The same calculation is also performed for high (perpendicular) magnetic Prandtl number plasmas ($Pm$), where the reconnection scale is found to be $lambda/Lsim S_L^{-4/7}Pm^{-2/7}$.