Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userFirst-principles demonstration of diffusive-advective particle acceleration in kinetic simulations of relativistic plasma turbulence
119
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Nonthermal relativistic plasmas are ubiquitous in astrophysical systems like pulsar wind nebulae and active galactic nuclei, as inferred from their emission spectra. The underlying nonthermal particle acceleration (NTPA) processes have traditionally been modeled with a Fokker-Planck (FP) diffusion-advection equation in momentum space. In this paper, we directly test the FP framework in ab-initio kinetic simulations of driven magnetized turbulence in relativistic pair plasma. By statistically analyzing the motion of tracked particles, we demonstrate the diffusive nature of NTPA and measure the FP energy diffusion ($D$) and advection ($A$) coefficients as functions of particle energy $gamma m_e c^2$. We find that $D(gamma)$ scales as $gamma^2$ in the high-energy nonthermal tail, in line with 2nd-order Fermi acceleration theory, but has a much weaker scaling at lower energies. We also find that $A$ is not negligible and reduces NTPA by tending to pull particles towards the peak of the particle energy distribution. This study provides strong support for the FP picture of turbulent NTPA, thereby enhancing our understanding of space and astrophysical plasmas.
rate research
Read More
Numerical algorithms to load relativistic Maxwell distributions in particle-in-cell (PIC) and Monte-Carlo simulations are presented. For stationary relativistic Maxwellian, the inverse transform method and the Sobol algorithm are reviewed. To boost particles to obtain relativistic shifted-Maxwellian, two rejection methods are proposed in a physically transparent manner. Their acceptance efficiencies are ${approx}50%$ for generic cases and $100%$ for symmetric distributions. They can be combined with arbitrary base algorithms.
Magnetic reconnection, especially in the relativistic regime, provides an efficient mechanism for accelerating relativistic particles and thus offers an attractive physical explanation for nonthermal high-energy emission from various astrophysical sources. I present a simple analytical model that elucidates key physical processes responsible for reconnection-driven relativistic nonthermal particle acceleration (NTPA) in the large-system, plasmoid-dominated regime in two dimensions. The model aims to explain the numerically-observed dependencies of the power-law index $p$ and high-energy cutoff $gamma_c$ of the resulting nonthermal particle energy spectrum $f(gamma)$ on the ambient plasma magnetization $sigma$, and (for $gamma_c$) on the system size $L$. In this self-similar model, energetic particles are continuously accelerated by the out-of-plane reconnection electric field $E_{rm rec}$ until they become magnetized by the reconnected magnetic field and eventually trapped in plasmoids large enough to confine them. The model also includes diffusive Fermi acceleration by particle bouncing off rapidly moving plasmoids. I argue that the balance between electric acceleration and magnetization controls the power-law index, while trapping in plasmoids governs the cutoff, thus tying the particle energy spectrum to the plasmoid distribution.
Dispersive Alfven waves (DAWs) offer, an alternative to magnetic reconnection, opportunity to accelerate solar flare particles. We study the effect of DAW polarisation, L-, R-, circular and elliptical, in different regimes inertial and kinetic on the efficiency of particle acceleration. We use 2.5D PIC simulations to study how particles are accelerated when DAW, triggered by a solar flare, propagates in transversely inhomogeneous plasma that mimics solar coronal loop. (i) In inertial regime, fraction of accelerated electrons (along the magnetic field), in density gradient regions is ~20% by the time when DAW develops 3 wavelengths and is increasing to ~30% by the time DAW develops 13 wavelengths. In all considered cases ions are heated in transverse to the magnetic field direction and fraction of the heated particles is ~35%. (ii) The case of R-circular, L- and R- elliptical polarisation DAWs, with the electric field in the non-ignorable transverse direction exceeding several times that of in the ignorable direction, produce more pronounced parallel electron beams and transverse ion beams in the ignorable direction. In the inertial regime such polarisations yield the fraction of accelerated electrons ~20%. In the kinetic regime this increases to ~35%. (iii) The parallel electric field that is generated in the density inhomogeneity regions is independent of m_i/m_e and exceeds the Dreicer value by 8 orders of magnitude. (iv) Electron beam velocity has the phase velocity of the DAW. Thus electron acceleration is via Landau damping of DAWs. For the Alfven speeds of 0.3c the considered mechanism can accelerate electrons to energies circa 20 keV. (v) The increase of mass ratio from m_i/m_e=16 to 73.44 increases the fraction of accelerated electrons from 20% to 30-35% (depending on DAW polarisation). For the mass ratio m_i/m_e=1836 the fraction of accelerated electrons would be >35%.
A theoretical framework for low-frequency electromagnetic (drift-)kinetic turbulence in a collisionless, multi-species plasma is presented. The result generalises reduced magnetohydrodynamics (RMHD) and kinetic RMHD (Schekochihin et al. 2009) for pressure-anisotropic plasmas, allowing for species drifts---a situation routinely encountered in the solar wind and presumably ubiquitous in hot dilute astrophysical plasmas (e.g. intracluster medium). Two main objectives are achieved. First, in a non-Maxwellian plasma, the relationships between fluctuating fields (e.g., the Alfven ratio) are order-unity modified compared to the more commonly considered Maxwellian case, and so a quantitative theory is developed to support quantitative measurements now possible in the solar wind. The main physical feature of low-frequency plasma turbulence survives the generalisation to non-Maxwellian distributions: Alfvenic and compressive fluctuations are energetically decoupled, with the latter passively advected by the former; the Alfvenic cascade is fluid, satisfying RMHD equations (with the Alfven speed modified by pressure anisotropy and species drifts), whereas the compressive cascade is kinetic and subject to collisionless damping. Secondly, the organising principle of this turbulence is elucidated in the form of a generalised kinetic free-energy invariant. It is shown that non-Maxwellian features in the distribution function reduce the rate of phase mixing and the efficacy of magnetic stresses; these changes influence the partitioning of free energy amongst the various cascade channels. As the firehose or mirror instability thresholds are approached, the dynamics of the plasma are modified so as to reduce the energetic cost of bending magnetic-field lines or of compressing/rarefying them. Finally, it is shown that this theory can be derived as a long-wavelength limit of non-Maxwellian slab gyrokinetics.
Hot accretion flows contain collisionless plasmas that are believed to be capable of accelerating particles to very high energies, as a result of turbulence generated by the magnetorotational instability (MRI). We conduct unstratified shearing-box simulations of the MRI turbulence in ideal magnetohydrodynamic (MHD), and inject energetic (relativistic) test particles in simulation snapshots to conduct a detailed investigation on particle diffusion and stochastic acceleration. We consider different amount of net vertical magnetic flux to achieve different disk magnetizations levels at saturated states, with sufficiently high resolution to resolve the gyro-radii ($R_g$) of most particles. Particles with large $R_g$ ($gtrsim0.03$ disk scale height $H$) show spatial diffusion coefficients of $sim30$ and $sim5$ times Bohm values in the azimuthal and poloidal directions, respectively. We further measure particle momentum diffusion coefficient $D(p)$ by applying the Fokker-Planck equation to particle momentum evolution. For these particles, contribution from turbulent fluctuations scales as $D(p)propto p$, and shear acceleration takes over when $R_ggtrsim0.1H$, characterized by $D(p)propto p^3$. For particles with smaller $R_g$ ($lesssim0.03H$), their spatial diffusion coefficients roughly scale as $sim p^{-1}$, and show evidence of $D(p)propto p^2$ scaling in momentum diffusion but with large uncertainties. We find that multiple effects contribute to stochastic acceleration/deceleration, and the process is also likely affected by intermittency in the MRI turbulence. We also discuss the potential of accelerating PeV cosmic-rays in hot accretion flows around supermassive black holes.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
