21 Search Results

Absmore » tract 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 $S^p\propto l^{{\zeta }_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 ₀ / B ₀ . We find that increasing values of σ and δ B ₀ / B ₀ 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.« less

Absmore » tract The recent discovery of astrophysical neutrinos from the Seyfert galaxy NGC 1068 suggests the presence of nonthermal protons within a compact “coronal” region close to the central black hole. The acceleration mechanism of these nonthermal protons remains elusive. We show that a large-scale magnetic reconnection layer, of the order of a few gravitational radii, may provide such a mechanism. In such a scenario, rough energy equipartition between magnetic fields, X-ray photons, and nonthermal protons is established in the reconnection region. Motivated by recent 3D particle-in-cell simulations of relativistic reconnection, we assume that the spectrum of accelerated protons is a broken power law, with the break energy being constrained by energy conservation (i.e., the energy density of accelerated protons is at most comparable to the magnetic energy density). The proton spectrum is ${\mathit{dn}}_p/{\mathit{dE}}_p\propto E_p^{-1}$ below the break and ${\mathit{dn}}_p/{\mathit{dE}}_p\propto E_p^{-s}$ above the break, with IceCube neutrino observations suggesting s ≃ 3. Protons above the break lose most of their energy within the reconnection layer via photohadronic collisions with the coronal X-rays, producing a neutrino signal in good agreement with the recent observations. Gamma rays injected in photohadronic collisions are cascaded to lower energies, sustaining the population of electron–positron pairs that makes the corona moderately Compton thick.« less

Turbulence is often invoked to explain the origin of nonthermal particles in space and astrophysical plasmas. By means of 3D fully kinetic particle-in-cell simulations, we demonstrate that turbulence in low-β plasmas (β is the ratio of plasma pressure to magnetic pressure) accelerates ions and electrons into a nonthermal energy distribution with a power-law energy range. The ion spectrum is harder than the electron one, and both distributions get steeper for higher β. We show that the energization of electrons is accompanied by a significant energy-dependent pitch-angle anisotropy, with most electrons moving parallel to the local magnetic field, while ions staymore » roughly isotropic. We demonstrate that particle injection from the thermal pool occurs in regions of high current density. Parallel electric fields associated with magnetic reconnection are responsible for the initial energy gain of electrons, whereas perpendicular electric fields control the overall energization of ions. Our findings have important implications for the origin of nonthermal particles in space and astrophysical plasmas.« less

Magnetohydrodynamic turbulence regulates the transfer of energy from large to small scales in many astrophysical systems, including the solar atmosphere. We perform three-dimensional magnetohydrodynamic simulations with unprecedentedly large magnetic Reynolds number to reveal how rapid reconnection of magnetic field lines changes the classical paradigm of the turbulent energy cascade. By breaking elongated current sheets into chains of small magnetic flux ropes (or plasmoids), magnetic reconnection leads to a previously undiscovered range of energy cascade, where the rate of energy transfer is controlled by the growth rate of the plasmoids. As a consequence, the turbulent energy spectra steepen and attain amore » spectral index of -2.2 that is accompanied by changes in the anisotropy of turbulence eddies. The omnipresence of plasmoids and their consequences on, for example, solar coronal heating, can be further explored with current and future spacecraft and telescopes.« less

Collisionless, magnetized turbulence offers a promising framework for the generation of nonthermal high-energy particles in various astrophysical sites. Yet, the detailed mechanism that governs particle acceleration has remained subject to debate. By means of 2D and 3D particle-in-cell, as well as 3D (incompressible) magnetohydrodynamic (MHD) simulations, we test here a recent model of nonresonant particle acceleration in strongly magnetized turbulence, which ascribes the energization of particles to their continuous interaction with the random velocity flow of the turbulence, in the spirit of the original Fermi model. To do so, we compare, for a large number of particles that were trackedmore » in the simulations, the predicted and the observed histories of particles momenta. The predicted history is that derived from the model, after extracting from the simulations, at each point along the particle trajectory, the three force terms that control acceleration: the acceleration of the field line velocity projected along the field line direction, its shear projected along the same direction, and its transverse compressive part. Overall, we find a clear correlation between the model predictions and the numerical experiments, indicating that this nonresonant model can successfully account for the bulk of particle energization through Fermi-type processes in strongly magnetized turbulence. Additionally we also observe that the parallel shear contribution tends to dominate the physics of energization in the particle-in-cell simulations, while in the magnetohydrodynamic incompressible simulation, both the parallel shear and the transverse compressive term provide about equal contributions.« less

Nature’s most powerful high-energy sources are capable of accelerating particles to high energy and radiating 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~1) within a few eddy turnover times. Low pitch-angle particles can significantly exceed the nominal radiation-reaction limit, before abruptly coolingmore » 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 (νF_ν ∝ ν^s with s~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.« less

Synchrotron emission from astrophysical nonthermal sources usually assumes that the emitting particles are isotropic. By means of large-scale two- and three-dimensional particle-in-cell simulations, we demonstrate that the dissipation of magnetically dominated ($${\sigma }_{0}\gg 1$$) turbulence in pair plasmas leads to strongly anisotropic particle distributions. At Lorentz factors $$\sim {\sigma }_{0}{\gamma }_{\mathrm{th}0}$$ (here, $${\gamma }_{\mathrm{th}0}$$ is the initial Lorentz factor), the particle velocity is preferentially aligned with the local magnetic field; instead, the highest energy particles are preferentially oriented in the plane perpendicular to the field. This energy-dependent anisotropy leads to a synchrotron spectral flux $$\nu {F}_{\nu }\propto {\nu }^{s}$$ thatmore » is much harder than for isotropic particles. Remarkably, for $${\sigma }_{0}\gg 1$$ we find that the angle-integrated spectral slope in the slow cooling regime is $$s\sim 0.5-0.7$$ for a wide range of turbulence fluctuations, $$0.25\lesssim \delta {B}_{\mathrm{rms}0}^{2}/{B}_{0}^{2}\lesssim 4$$, despite significant variations in the power-law energy spectrum of nonthermal particles. This is because weaker turbulence levels imprint a stronger degree of anisotropy, thereby counteracting the effect of the steeper particle spectrum. The synchrotron spectral slope may be even harder, $$s\gtrsim 0.7$$, if the observer is in the plane perpendicular to the mean magnetic field. Our results are independent of domain size and dimensionality. Further, our findings may help explain the origin of hard synchrotron spectra of astrophysical nonthermal sources, most notably the radio spectrum of pulsar wind nebulae.« less

We construct a turbulent model of the Crab Nebula's nonthermal emission. The present model resolves a number of long-standing problems of the Kennel–Coroniti model: (i) the sigma problem, (ii) the hard spectrum of radio electrons, (iii) the high peak energy of gamma-ray flares, (iv) and the spatial evolution of the infrared (IR) emission. The Nebula contains two populations of injected particles: Component-I, accelerated at the wind termination shock via the Fermi-I mechanism; and Component-II, accelerated in reconnecting turbulence in highly magnetized (σ >> 1) plasma in the central part of the Crab Nebula. The reconnecting turbulence in Component-II extends frommore » radio to gamma-rays: it accelerates radio electrons with a hard spectrum, destroys the large-scale magnetic flux (and thus resolves the sigma problem), and occasionally produces gamma-ray flares (from the largest-scale reconnection events). The model reproduces the broadband spectrum of the Crab Nebula, from low-frequency synchrotron emission in radio to inverse-Compton emission at TeV energies, as well as the spatially resolved evolution of the spectral indices in the IR and optical bands.« less

Magnetized turbulence and magnetic reconnection are often invoked to explain the nonthermal emission observed from a wide variety of astrophysical sources. By means of fully kinetic 2D and 3D particle-in-cell simulations, we investigate the interplay between turbulence and reconnection in generating nonthermal particles in magnetically dominated (or, equivalently, “relativistic”) pair plasmas. A generic by-product of the turbulence evolution is the generation of a nonthermal particle spectrum with a power-law energy range. The power-law slope p is harder for larger magnetizations and stronger turbulence fluctuations, and it can be as hard as p ≲ 2. The Larmor radius of particles atmore » the high-energy cutoff is comparable to the size l of the largest turbulent eddies. Plasmoid-mediated reconnection, which self-consistently occurs in the turbulent plasma, controls the physics of particle injection. Then, particles are further accelerated by stochastic scattering off turbulent fluctuations. The work done by parallel electric fields—naturally expected in reconnection layers—is responsible for most of the initial energy increase and is proportional to the magnetization σ of the system, while the subsequent energy gain, which dominates the overall energization of high-energy particles, is powered by the perpendicular electric fields of turbulent fluctuations. The two-stage acceleration process leaves an imprint in the particle pitch-angle distribution: lowenergy particles are aligned with the field, while the highest-energy particles move preferentially orthogonal to it. The energy diffusion coefficient of stochastic acceleration scales as Dγ ~ 0.1σ(c/l)γ², where γ is the particle Lorentz factor. This results in fast cceleration timescales t_acc ~ (3/σ)l/c. Our findings have important implications for understanding the generation of nonthermal particles in high-energy astrophysical sources.« less

Analytic scaling relations are derived for a phenomenological model of the plasmoid instability in an evolving current sheet, including the effects of reconnection outflow. Two scenarios are considered, where the plasmoid instability can be triggered either by an injected initial perturbation or by the natural noise of the system (here referred to as the system noise). The two scenarios lead to different scaling relations because the initial noise decays when the linear growth of the plasmoid instability is not sufficiently fast to overcome the advection loss caused by the reconnection outflow, whereas the system noise represents the lowest level ofmore » fluctuations in the system. The leading order approximation for the current sheet width at disruption takes the form of a power law multiplied by a logarithmic factor, and from that, the scaling relations for the wavenumber and the linear growth rate of the dominant mode are obtained. When the effects of the outflow are neglected, the scaling relations agree, up to the leading order approximation, with previously derived scaling relations based on a principle of least time. As a result, the analytical scaling relations are validated with numerical solutions of the model.« less