Flux-freezing breakdown in high-conductivity magnetohydrodynamic turbulence

Radiative Particle-in-Cell Simulations of Turbulent Comptonization in Magnetized Black-Hole Coronae

We report results from the first radiative particle-in-cell simulations of strong Alfvénic turbulence in plasmas of moderate optical depth. The simulations are performed in a local 3D periodic box and self-consistently follow the evolution of radiation as it interacts with a turbulent electron-positron plasma via Compton scattering. We focus on the conditions expected in magnetized coronae of accreting black holes and obtain an emission spectrum consistent with the observed hard state of Cyg X-1. Most of the turbulence power is transferred directly to the photons via bulk Comptonization, shaping the peak of the emission around 100 keV. The rest is released into nonthermal particles, which generate the MeV spectral tail. The method presented here shows promising potential for ab initio modeling of various astrophysical sources and opens a window into a new regime of kinetic plasma turbulence.

Extended Magnetic Reconnection in Kinetic Plasma Turbulence

Magnetic reconnection and plasma turbulence are ubiquitous processes important for laboratory, space, and astrophysical plasmas. Reconnection has been suggested to play an important role in the energetics and dynamics of turbulence by observations, simulations, and theory for two decades. The fundamental properties of reconnection at kinetic scales, essential to understanding the general problem of reconnection in magnetized turbulence, remain largely unknown at present. Here, we present an application of the magnetic flux transport method that can accurately identify reconnection in turbulence to a three-dimensional simulation. Contrary to ideas that reconnection in turbulence would be patchy and unpredictable, highly extended reconnection X lines, on the same order of magnitude as the system size, form at kinetic scales. Extended X lines develop through bidirectional reconnection spreading. They satisfy critical balance characteristic of turbulence, which predicts the X-line extent at a given scale. These results present a picture of fundamentally extended reconnection in kinetic-scale turbulence.

First-Principles Fermi Acceleration in Magnetized Turbulence

This Letter provides a concrete implementation of Fermi's model of particle acceleration in magnetohydrodynamic (MHD) turbulence, connecting the rate of energization to the gradients of the velocity of magnetic field lines, which it characterizes within a multifractal picture of turbulence intermittency. It then derives a transport equation in momentum space for the distribution function. This description is shown to be substantiated by a large-scale numerical simulation of strong MHD turbulence. The present general framework can be used to model particle acceleration in a variety of environments.

Energy transfer in reconnection and turbulence

Reconnection and turbulence are two of the most commonly observed dynamical processes in plasmas, but their relationship is still not fully understood. Using 2.5D kinetic particle-in-cell simulations of both strong turbulence and reconnection, we compare the cross-scale transfer of energy in the two systems by analyzing the generalization of the von Kármán Howarth equations for Hall magnetohydrodynamics, a formulation that subsumes the third-order law for steady energy transfer rates. Even though the large scale features are quite different, the finding is that the decomposition of the energy transfer is structurally very similar in the two cases. In the reconnection case, the time evolution of the energy transfer also exhibits a correlation with the reconnection rate. These results provide explicit evidence that reconnection dynamics fundamentally involves turbulence-like energy transfer.

Faster Form of Electron Magnetic Reconnection with a Finite Length X-Line

Observations in Earth's turbulent magnetosheath downstream of a quasiparallel bow shock reveal a prevalence of electron-scale current sheets favorable for electron-only reconnection where ions are not coupled to the reconnecting magnetic fields. In small-scale turbulence, magnetic structures associated with intense current sheets are limited in all dimensions. And since the coupling of ions are constrained by a minimum length scale, the dynamics of electron reconnection is likely to be 3D. Here, both 2D and 3D kinetic particle-in-cell simulations are used to investigate electron-only reconnection, focusing on the reconnection rate and associated electron flows. A new form of 3D electron-only reconnection spontaneously develops where the magnetic X-line is localized in the out-of-plane (z) direction. The consequence is an enhancement of the reconnection rate compared with two dimensions, which results from differential mass flux out of the diffusion region along z, enabling a faster inflow velocity and thus a larger reconnection rate. This outflow along z is due to the magnetic tension force in z just as the conventional exhaust tension force, allowing particles to leave the diffusion region efficiently along z unlike the 2D configuration.

Magnetic reconnection in partially ionized plasmas

Magnetic reconnection has been intensively studied in fully ionized plasmas. However, plasmas are often partially ionized in astrophysical environments. The interactions between the neutral particles and ionized plasmas might strongly affect the reconnection mechanisms. We review magnetic reconnection in partially ionized plasmas in different environments from theoretical, numerical, observational and experimental points of view. We focus on mechanisms which make magnetic reconnection fast enough to compare with observations, especially on the reconnection events in the low solar atmosphere. The heating mechanisms and the related observational evidence of the reconnection process in the partially ionized low solar atmosphere are also discussed. We describe magnetic reconnection in weakly ionized astrophysical environments, including the interstellar medium and protostellar discs. We present recent achievements about fast reconnection in laboratory experiments for partially ionized plasmas.

Statistical analysis of stochastic magnetic fields

Previous work has introduced scale-split energy density ψ_{l,L}(x,t)=1/2B_{l}·B_{L} for vector field B(x,t) coarse grained at scales l and L, in order to quantify the field stochasticity or spatial complexity. In this formalism, the L_{p} norms S_{p}(t)=1/2||1-B[over ̂]_{l}·B[over ̂]_{L}||_{p}, pth-order stochasticity level, and E_{p}(t)=1/2||B_{l}B_{L}||_{p}, pth order mean cross energy density, are used to analyze the evolution of the stochastic field B(x,t). Application to turbulent magnetic fields leads to the prediction that turbulence in general tends to tangle an initially smooth magnetic field increasing the magnetic stochasticity level, ∂_{t}S_{p}>0. An increasing magnetic stochasticity in turn leads to disalignments of the coarse-grained fields B_{d} at smaller scales, d≪L, thus they average to weaker fields B_{L} at larger scales upon coarse graining, i.e., ∂_{t}E_{p}<0. Magnetic field resists the tangling effect of the turbulence by means of magnetic tension force. This can lead at some point to a sudden slippage between the field and fluid, decreasing the stochasticity ∂_{t}S_{p}<0 and increasing the energy ∂_{t}E_{p}>0 by aligning small-scale fields B_{d}. Thus the maxima (minima) of magnetic stochasticity are expected to approximately coincide with the minima (maxima) of cross energy density, occurrence of which corresponds to slippage of the magnetic field through the fluid. In this formalism, magnetic reconnection and field-fluid slippage both correspond to T_{p}=∂_{t}S_{p}=0and∂_{t}T_{2}<0. Previous work has also linked field-fluid slippage to magnetic reconnection invoking totally different approaches. In this paper, (a) we test these theoretical predictions numerically using a homogeneous, incompressible magnetohydrodynamic (MHD) simulation. Apart from expected small-scale deviations, possibly due to, e.g., intermittency and strong field annihilation, the theoretically predicted global relationship between stochasticity and cross energy is observed in different subvolumes of the simulation box. This indicate ubiquitous local field-fluid slippage and reconnection events in MHD turbulence. In addition, (b) we show that the conditions T_{p}=∂_{t}S_{p}=0and∂_{t}T_{p}<0 lead to sudden increases in kinetic stochasticity level, i.e., τ_{p}=∂_{t}s_{p}(t)>0 with s_{p}(t)=1/2||1-u[over ̂]_{l}.u[over ̂]_{L}||_{p}, which may correspond to fluid jets spontaneously driven by sudden field-fluid slippage-magnetic reconnection. Otherwise, they may correspond only to field-fluid slippage without energy dissipation. This picture, therefore, suggests defining reconnection as field-fluid slippage (changes in S_{p}) accompanied with magnetic energy dissipation (changes in E_{p}). All in all, these provide a statistical approach to the reconnection in terms of the time evolution of magnetic and kinetic stochasticities, S_{p} and s_{p}, their time derivatives, T_{p}=∂_{t}S_{p}, τ_{p}=∂_{t}s_{p}, and corresponding cross energies, E_{p}, e_{p}(t)=1/2||u_{l}u_{L}||_{p}. Furthermore, (c) we introduce the scale-split magnetic helicity based on which we discuss the energy or stochasticity relaxation of turbulent magnetic fields-a generalized Taylor relaxation. Finally, (d) we construct and numerically test a toy model, which resembles a classical version of quantum mean field Ising model for magnetized fluids, in order to illustrate how turbulent energy can affect magnetic stochasticity in the weak field regime.

Magnetic stochasticity and diffusion

We develop a quantitative relationship between magnetic diffusion and the level of randomness, or stochasticity, of the diffusing magnetic field in a magnetized medium. A general mathematical formulation of magnetic stochasticity in turbulence has been developed in previous work in terms of the L_{p} norm S_{p}(t)=1/2∥1-B[over ̂]_{l}·B[over ̂]_{L}∥_{p}, pth-order magnetic stochasticity of the stochastic field B(x,t), based on the coarse-grained fields B_{l} and B_{L} at different scales l≠L. For laminar flows, the stochasticity level becomes the level of field self-entanglement or spatial complexity. In this paper, we establish a connection between magnetic stochasticity S_{p}(t) and magnetic diffusion in magnetohydrodynamic (MHD) turbulence and use a homogeneous, incompressible MHD simulation to test this prediction. Our results agree with the well-known fact that magnetic diffusion in turbulent media follows the superlinear Richardson dispersion scheme. This is intimately related to stochastic magnetic reconnection in which superlinear Richardson diffusion broadens the matter outflow width and accelerates the reconnection process.

Topology and stochasticity of turbulent magnetic fields

We present a mathematical formalism for the topology and stochasticity of vector fields based on renormalization group methodology. The concept of a scale-split energy density, ψ_{l,L}=B_{l}·B_{L}/2 for vector field B(x,t) renormalized at scales l and L, is introduced in order to quantify the notion of the field topological deformation, topology change, and stochasticity level. In particular, for magnetic fields, it is shown that the evolution of the field topology is directly related to the field-fluid slippage, which has already been linked to magnetic reconnection in previous work. The magnitude and direction of stochastic magnetic fields, shown to be governed, respectively, by the parallel and vertical components of the renormalized induction equation with respect to the magnetic field, can be studied separately by dividing ψ_{l,L} into two (3+1)-dimensional scalar fields. The velocity field can be approached in a similar way. Magnetic reconnection can then be defined in terms of the extrema of the L_{p} norms of these scalar fields. This formulation in fact clarifies different definitions of magnetic reconnection, which vaguely rely on the magnetic field topology, stochasticity, and energy conversion. Our results support the well-founded yet partly overlooked picture in which magnetic reconnection in turbulent fluids occurs on a wide range of scales as a result of nonlinearities at large scales (turbulence inertial range) and nonidealities at small scales (dissipative range). Lagrangian particle trajectories, as well as magnetic field lines, are stochastic in turbulent magnetized media in the limit of small resistivity and viscosity. The magnetic field tends to reduce its stochasticity induced by the turbulent flow by slipping through the fluid, which may accelerate fluid particles. This suggests that reconnection is a relaxation process by which the magnetic field lowers both its topological entanglements induced by turbulence and its energy level.

Local approach to the study of energy transfers in incompressible magnetohydrodynamic turbulence

We present a local approach to the study of scale-to-scale energy transfers in magnetohydrodynamic (MHD) turbulence. This approach is based on performing local averages of the physical fields, which amounts to filtering scales smaller than some parameter ℓ. A key step in this work is the derivation of a local Kármán-Howarth-Monin relation which provides a local form of Politano and Pouquet's 4/3 law, without any assumption of homogeneity or isotropy. Our approach is exact and nonrandom, and we show its connection to the usual statistical results of turbulence. Its implementation on data obtained via a three-dimensional direct numerical simulation of the forced incompressible MHD equations from the John Hopkins turbulence database constitutes the main part of our study. First, we show that the local Kármán-Howarth-Monin relation holds well. The space statistics of local cross-scale transfers are studied next, their means and standard deviations being maximum at inertial scales and their probability density functions (PDFs) displaying very wide tails. Events constituting the tails of the PDFs are shown to form structures of strong transfers, either positive or negative, which can be observed over the whole available range of scales. As ℓ is decreased, these structures become more and more localized in space while contributing to an increasing fraction of the mean energy cascade rate. Finally, we highlight their quasi-one-dimensional (filamentlike) or quasi-two-dimensional (sheetlike or ribbonlike) nature and show that they appear in areas of strong vorticity or electric current density.

From Primordial Seed Magnetic Fields to the Galactic Dynamo

The origin and maintenance of coherent magnetic fields in the Universe is reviewed with an emphasis on the possible challenges that arise in their theoretical understanding. We begin with the interesting possibility that magnetic fields originated at some level from the early universe. This could be during inflation, the electroweak, or the quark-hadron phase transitions. These mechanisms can give rise to fields which could be strong, but often with much smaller coherence scales than galactic scales. Their subsequent turbulent decay decreases their strength but increases their coherence. We then turn to astrophysical batteries which can generate seed magnetic fields. Here the coherence scale can be large, but the field strength is generally very small. These seed fields need to be further amplified and maintained by a dynamo to explain observed magnetic fields in galaxies. Basic ideas behind both small and large-scale turbulent dynamos are outlined. The small-scale dynamo may help to understand the first magnetization of young galaxies, while the large-scale dynamo is important for the generation of fields with scales larger than the stirring scale, as observed in nearby disk galaxies. The current theoretical challenges that turbulent dynamos encounter and their possible resolution are discussed.

Decoupled Cascades of Kinetic and Magnetic Energy in Magnetohydrodynamic Turbulence

Magnetic energy (ME) and kinetic energy (KE) in ideal incompressible magnetohydrodynamics are not global invariants and, therefore, it has been justified to discuss only the cascade of their sum total energy. We provide a physical argument based on scale locality, along with compelling evidence that ME and KE budgets statistically decouple beyond a transitional "conversion" range. This arises because magnetic field-line stretching is a large-scale process which vanishes on average at intermediate and small scales within the inertial-inductive range, thereby allowing each of the mean ME and KE to cascade conservatively and at an equal rate, yielding a turbulent magnetic Prandtl number of unity over these scales.

Effective Resistivity in Collisionless Magnetic Reconnection

An effective resistivity relevant to collisionless magnetic reconnection (MR) in plasma is presented. It is based on the argument that pitch angle scattering of electrons in the small electron diffusion region around the X line can lead to an effective, resistivity in collisionless plasma. The effective resistivity so obtained is in the form of a power law of the local plasma and magnetic field parameters. Its validity is confirmed by direct collisionless particle-in-cell (PIC) simulation. The result agrees very well with the resistivity (obtained from available data) of a large number of environments susceptible to MR: from the intergalactic and interstellar to solar and terrestrial to laboratory fusion plasmas. The scaling law can readily be incorporated into existing collisional magnetohydrodynamic simulation codes to investigate collisionless MR, as well as serve as a guide to ab initio theoretical investigations of the collisionless MR process.

Kinetic Turbulence in Relativistic Plasma: From Thermal Bath to Nonthermal Continuum

We present results from particle-in-cell simulations of driven turbulence in magnetized, collisionless, and relativistic pair plasmas. We find that the fluctuations are consistent with the classical k_{⊥}^{-5/3} magnetic energy spectrum at fluid scales and a steeper k_{⊥}^{-4} spectrum at sub-Larmor scales, where k_{⊥} is the wave vector perpendicular to the mean field. We demonstrate the development of a nonthermal, power-law particle energy distribution f(E)∼E^{-α}, with an index α that decreases with increasing magnetization and increases with an increasing system size (relative to the characteristic Larmor radius). Our simulations indicate that turbulence can be a viable source of energetic particles in high-energy astrophysical systems, such as pulsar wind nebulae, if scalings asymptotically become insensitive to the system size.

Perspectives on magnetic reconnection

Magnetic reconnection is a topological rearrangement of magnetic field that occurs on time scales much faster than the global magnetic diffusion time. Since the field lines break on microscopic scales but energy is stored and the field is driven on macroscopic scales, reconnection is an inherently multi-scale process that often involves both magnetohydrodynamic (MHD) and kinetic phenomena. In this article, we begin with the MHD point of view and then describe the dynamics and energetics of reconnection using a two-fluid formulation. We also focus on the respective roles of global and local processes and how they are coupled. We conclude that the triggers for reconnection are mostly global, that the key energy conversion and dissipation processes are either local or global, and that the presence of a continuum of scales coupled from microscopic to macroscopic may be the most likely path to fast reconnection.

Explosive Particle Dispersion in Plasma Turbulence

Particle dynamics are investigated in plasma turbulence, using self-consistent kinetic simulations, in two dimensions. In the steady state, the trajectories of single protons and proton pairs are studied, at different values of plasma β (ratio between kinetic and magnetic pressure). For single-particle displacements, results are consistent with fluids and magnetic field line dynamics, where particles undergo normal diffusion for very long times, with higher β's being more diffusive. In an intermediate time range, with separations lying in the inertial range, particles experience an explosive dispersion in time, consistent with the Richardson prediction. These results, obtained for the first time with a self-consistent kinetic model, are relevant for astrophysical and laboratory plasmas, where turbulence is crucial for heating, mixing, and acceleration processes.

The origin, evolution and signatures of primordial magnetic fields

The universe is magnetized on all scales probed so far. On the largest scales, galaxies and galaxy clusters host magnetic fields at the micro Gauss level coherent on scales up to ten kpc. Recent observational evidence suggests that even the intergalactic medium in voids could host a weak  ∼  10(-16) Gauss magnetic field, coherent on Mpc scales. An intriguing possibility is that these observed magnetic fields are a relic from the early universe, albeit one which has been subsequently amplified and maintained by a dynamo in collapsed objects. We review here the origin, evolution and signatures of primordial magnetic fields. After a brief summary of magnetohydrodynamics in the expanding universe, we turn to magnetic field generation during inflation and phase transitions. We trace the linear and nonlinear evolution of the generated primordial fields through the radiation era, including viscous effects. Sensitive observational signatures of primordial magnetic fields on the cosmic microwave background, including current constraints from Planck, are discussed. After recombination, primordial magnetic fields could strongly influence structure formation, especially on dwarf galaxy scales. The resulting signatures on reionization, the redshifted 21 cm line, weak lensing and the Lyman-α forest are outlined. Constraints from radio and γ-ray astronomy are summarized. Astrophysical batteries and the role of dynamos in reshaping the primordial field are briefly considered. The review ends with some final thoughts on primordial magnetic fields.

The heterogeneous energy landscape expression of KWW relaxation

Here we show a heterogeneous energy landscape approach to describing the Kohlrausch-Williams-Watts (KWW) relaxation function. For a homogeneous dynamic process, the distribution of free energy landscape is first proposed, revealing the significance of rugged fluctuations. In view of the heterogeneous relaxation given in two dynamic phases and the transmission coefficient in a rate process, we obtain a general characteristic relaxation time distribution equation for the KWW function in a closed, analytic form. Analyses of numerical computation show excellent accuracy, both in time and frequency domains, in the convergent performance of the heterogeneous energy landscape expression and shunning the catastrophic truncations reported in the previous work. The stretched exponential β, closely associated to temperature and apparent correlation with one dynamic phase, reveals a threshold value of 1/2 defining different behavior of the probability density functions. Our work may contribute, for example, to in-depth comprehension of the dynamic mechanism of glass transition, which cannot be provided by existing approaches.

Turbulent reconnection and its implications

Magnetic reconnection is a process of magnetic field topology change, which is one of the most fundamental processes happening in magnetized plasmas. In most astrophysical environments, the Reynolds numbers corresponding to plasma flows are large and therefore the transition to turbulence is inevitable. This turbulence, which can be pre-existing or driven by magnetic reconnection itself, must be taken into account for any theory of magnetic reconnection that attempts to describe the process in the aforementioned environments. This necessity is obvious as three-dimensional high-resolution numerical simulations show the transition to the turbulence state of initially laminar reconnecting magnetic fields. We discuss ideas of how turbulence can modify reconnection with the focus on the Lazarian & Vishniac (Lazarian & Vishniac 1999 Astrophys. J. 517, 700-718 (doi:10.1086/307233)) reconnection model. We present numerical evidence supporting the model and demonstrate that it is closely connected to the experimentally proven concept of Richardson dispersion/diffusion as well as to more recent advances in understanding of the Lagrangian dynamics of magnetized fluids. We point out that the generalized Ohm's law that accounts for turbulent motion predicts the subdominance of the microphysical plasma effects for reconnection for realistically turbulent media. We show that one of the most dramatic consequences of turbulence is the violation of the generally accepted notion of magnetic flux freezing. This notion is a cornerstone of most theories dealing with magnetized plasmas, and therefore its change induces fundamental shifts in accepted paradigms, for instance, turbulent reconnection entails reconnection diffusion process that is essential for understanding star formation. We argue that at sufficiently high Reynolds numbers the process of tearing reconnection should transfer to turbulent reconnection. We discuss flares that are predicted by turbulent reconnection and relate this process to solar flares and γ-ray bursts. With reference to experiments, we analyse solar observations in situ as measurements in the solar wind or heliospheric current sheet and show the correspondence of data with turbulent reconnection predictions. Finally, we discuss first-order Fermi acceleration of particles that is a natural consequence of the turbulent reconnection.

Asymptotic results for backwards two-particle dispersion in a turbulent flow

We derive an exact equation governing two-particle backwards mean-squared dispersion for both deterministic and stochastic tracer particles in turbulent flows. For the deterministic trajectories, we probe the consequences of our formula for short times and arrive at approximate expressions for the mean-squared dispersion which involve second order structure functions of the velocity and acceleration fields. For the stochastic trajectories, we analytically compute an exact t3 contribution to the squared separation of stochastic paths. We argue that this contribution appears also for deterministic paths at long times and present direct numerical simulation results for incompressible Navier-Stokes flows to support this claim. We also numerically compute the probability distribution of particle separations for the deterministic paths and the stochastic paths and show their strong self-similar nature.