Turbulent Particle Acceleration Research Articles

Intermittency and Dissipative Structures Arising from Relativistic Magnetized Turbulence

Kinetic simulations of relativistic turbulence have significantly advanced our understanding of turbulent particle acceleration. Recent progress has highlighted the need for an updated acceleration theory that can account for particle acceleration within the plasma’s coherent structures. Here, we investigate how intermittency modeling connects statistical fluctuations in turbulence to regions of high-energy dissipation. This connection is established by employing a generalized She–Leveque model to characterize the exponents ζ p for the structure functions Sp∝lζp . The fitting of the scaling exponents provides us with a measure of the codimension of the dissipative structures, for which we subsequently determine the filling fraction. We perform our analysis for a range of magnetizations σ and relative fluctuation amplitudes δ B 0/B 0. We find that increasing values of σ and δ B 0/B 0 allow the turbulent cascade to break sheetlike structures into smaller regions of dissipation that resemble chains of flux ropes. However, as their dissipation measure increases, the dissipative regions become less volume filling. With this work, we aim to inform future turbulent acceleration theories that incorporate particle energization from interactions with coherent structures within relativistic turbulence.

System-size Convergence of Nonthermal Particle Acceleration in Relativistic Plasma Turbulence

Abstract We apply collisionless particle-in-cell simulations of relativistic pair plasmas to explore whether driven turbulence is a viable high-energy astrophysical particle accelerator. We characterize nonthermal particle distributions for varying system sizes up to L/2πρ e0 = 163, where L/2π is the driving scale and ρ e0 is the initial characteristic Larmor radius. We show that turbulent particle acceleration produces power-law energy distributions that, when compared at a fixed number of large-scale dynamical times, slowly steepen with increasing system size. We demonstrate, however, that convergence is obtained by comparing the distributions at different times that increase with system size (approximately logarithmically). We suggest that the system-size dependence arises from the time required for particles to reach the highest accessible energies via Fermi acceleration. The converged power-law index of the energy distribution, α ≈ 3.0 for magnetization σ = 3/8, makes turbulence a possible explanation for nonthermal spectra observed in systems such as the Crab Nebula.

Inertial-particle accelerations in turbulence: a Lagrangian closure

The distribution of particle accelerations in turbulence is intermittent, with non-Gaussian tails that are quite different for light and heavy particles. In this article we analyse a closure scheme for the acceleration fluctuations of light and heavy inertial particles in turbulence, formulated in terms of Lagrangian correlation functions of fluid tracers. We compute the variance and the flatness of inertial-particle accelerations and we discuss their dependency on the Stokes number. The closure incorporates effects induced by the Lagrangian correlations along the trajectories of fluid tracers, and its predictions agree well with results of direct numerical simulations of inertial particles in turbulence, provided that the effects induced by inertial preferential sampling of heavy/light particles outside/inside vortices are negligible. In particular, the scheme predicts the correct functional behaviour of the acceleration variance, as a function of $St$, as well as the presence of a minimum/maximum for the flatness of the acceleration of heavy/light particles, in good qualitative agreement with numerical data. We also show that the closure works well when applied to the Lagrangian evolution of particles using a stochastic surrogate for the underlying Eulerian velocity field. Our results support the conclusion that there exist important contributions to the statistics of the acceleration of inertial particles independent of the preferential sampling. For heavy particles we observe deviations between the predictions of the closure scheme and direct numerical simulations, at Stokes numbers of order unity. For light particles the deviation occurs for larger Stokes numbers.

Intrinsic rotation in tokamaks: theory

Self-consistent equations for intrinsic rotation in tokamaks with small poloidal magnetic field Bp compared to the total magnetic field B are derived. The model gives the momentum redistribution due to turbulence, collisional transport and energy injection. Intrinsic rotation is determined by the balance between the momentum redistribution and the turbulent diffusion and convection. Two different turbulence regimes are considered: turbulence with characteristic perpendicular lengths of the order of the ion gyroradius, ρi, and turbulence with characteristic lengths of the order of the poloidal gyroradius, (B/Bp)ρi. Intrinsic rotation driven by gyroradius scale turbulence is mainly due to the effect of neoclassical corrections and of finite orbit widths on turbulent momentum transport, whereas for the intrinsic rotation driven by poloidal gyroradius scale turbulence, the slow variation of turbulence characteristics in the radial and poloidal directions and the turbulent particle acceleration can be become as important as the neoclassical and finite orbit width effects. The magnetic drift is shown to be indispensable for the intrinsic rotation driven by the slow variation of turbulence characteristics and the turbulent particle acceleration. The equations are written in a form conducive to implementation in a flux tube code, and the effect of the radial variation of the turbulence is included in a novel way that does not require a global gyrokinetic formalism.

Particle Acceleration in Turbulence and Weakly Stochastic Reconnection

In this Letter we analyze the energy distribution evolution of test particles injected in three dimensional (3D) magnetohydrodynamic (MHD) simulations of different magnetic reconnection configurations. When considering a single Sweet-Parker topology, the particles accelerate predominantly through a first-order Fermi process, as predicted in and demonstrated numerically in . When turbulence is included within the current sheet, the acceleration rate is highly enhanced, because reconnection becomes fast and independent of resistivity and allows the formation of a thick volume filled with multiple simultaneously reconnecting magnetic fluxes. Charged particles trapped within this volume suffer several head-on scatterings with the contracting magnetic fluctuations, which significantly increase the acceleration rate and results in a first-order Fermi process. For comparison, we also tested acceleration in MHD turbulence, where particles suffer collisions with approaching and receding magnetic irregularities, resulting in a reduced acceleration rate. We argue that the dominant acceleration mechanism approaches a second order Fermi process in this case.

MULTIWAVELENGTH OBSERVATIONS OF THE GAMMA-RAY BLAZAR PKS 0528+134 IN QUIESCENCE

We present multiwavelength observations of the ultraluminous blazar-type radio loud quasar PKS 0528+134 in quiescence during the period July to December 2009. Significant flux variability on a time scale of several hours was found in the optical regime, accompanied by a weak trend of spectral softening with increasing flux. We suggest that this might be the signature of a contribution from the accretion disk at the blue end of the optical spectrum. The optical flux is weakly polarized with rapid variations of the degree and direction of polarization, while the polarization of the 43 GHz radio core remains steady. Optical spectropolarimetry suggests a trend of increasing degree of polarization with increasing wavelength, providing additional evidence for an accretion disc contribution towards the blue end of the optical spectrum. We constructed four SEDs indicating that even in the quiescent state, the bolometric luminosity of PKS 0528+134 is dominated by its gamma-ray emission. A leptonic single-zone jet model produced acceptable fits to the SEDs with contributions to the high-energy emission from synchrotron self-Compton radiation and Comptonization of direct accretion disk emission. Fit parameters close to equipartition were obtained. The moderate variability on long time scales implies the existence of on-going particle acceleration, while the observed optical polarization variability seems to point towards a turbulent acceleration process. Turbulent particle acceleration at stationary features along the jet therefore appears to be a viable possibility for the quiescent state of PKS 0528+134.

Multifractal analysis of fluid particle accelerations in turbulence

The probability density function (PDF) of accelerations in turbulence is derived analytically with the help of the multifractal analysis based on generalized entropy, i.e., the Tsallis or the Rényi entropy. It is shown that the derived PDF explains quite well the one obtained by Bodenschatz et al. in the measurement of fluid particle accelerations in the Lagrangian frame at R λ =690, and the one by Gotoh et al. in the DNS with the mesh size 1024 3 at R λ =380.

A role of the boundary shear layer in modeling of large scale jets

We discuss a role of the boundary shear layer of large scale jets for their observed multiwavelength emission. We consider a simple mechanism of a turbulent particle acceleration acting within such regions as an alternative mechanism with respect to the standard approaches involving internal shocks in the jet. We show that, under simple assumptions, the boundary layer acceleration can provide very high energy electrons possibly generating the X-ray emission observed by Chandra from some large scale jets in powerful radiogalaxies.