Search

We investigate the role of the magnetohydrodynamic (MHD) turbulence measured by Voyager in the very local interstellar medium (VLISM) in modeling the Interstellar Boundary Explorer (IBEX) ribbon. We demonstrate that the mirroring by compressible modes of MHD turbulence dominates over that by the mean magnetic field. Based on the new mirror diffusion mechanism identified by Lazarian and Xu for particles with large pitch angles in MHD turbulence, we find that the mirror diffusion can both confine pickup ions and preserve their initial pitch angles, and thus accounts for the enhanced intensity of energetic neutral atoms that return to the heliosphere. The ribbon width is determined by both the range of pitch angles for effective turbulent mirroring and the field line wandering induced by Alfv\'{e}nic modes. It in turn provides a constraint on the amplitude of magnetic fluctuations of fast modes. The field line wandering also affects the coherence of the ribbon structure across the sky. By extrapolating the magnetic energy spectrum measured by Voyager, we find that the injection scale of the turbulence in the VLISM is less than $\sim 500$ au for the ribbon structure to be coherent., Comment: 8 pages, 3 figures, accepted for publication in ApJ

The theoretical prediction that magnetic reconnection spontaneously drives turbulence has been supported by magnetohydrodynamic (MHD) and kinetic simulations. While reconnection with externally driven turbulence is accepted as an effective mechanism for particle acceleration, the acceleration during the reconnection with self-driven turbulence is studied for the first time in this work. By using high-resolution 3D MHD simulations of reconnection with self-generated turbulence, we inject test particles into the reconnection layer to study their acceleration process. We find that the energy gain of the particles takes place when they bounce back and forth between converging turbulent magnetic fields. The particles can be efficiently accelerated in self-driven turbulent reconnection with the energy increase by about 3 orders of magnitude in the range of the box size. The acceleration proceeds when the particle gyroradii become larger than the thickness of the reconnection layer. We find that the acceleration in the direction perpendicular to the local magnetic field dominates over that in the parallel direction. The energy spectrum of accelerated particles is time-dependent with a slope that evolves toward -2.5. Our findings can have important implications for particle acceleration in high-energy astrophysical environments., Comment: 13 pages, 5 Figures. Accepted for publication in Journal of High Energy Astrophysics

In a weakly compressible high-$\beta$ medium, pitch-angle scattering and the associated scattering acceleration of cosmic rays (CRs) by anisotropic Alfv\'{e}n and slow modes of magnetohydrodynamic (MHD) turbulence is inefficient. To tap the energy from magnetic compressions for efficient particle acceleration, a diffusion mechanism that can effectively confine particles in space without causing their trapping or pitch-angle isotropization is needed. We find that the mirror diffusion in MHD turbulence recently identified in Lazarian and Xu (2021) satisfies all the above conditions and serves as a promising diffusion mechanism for efficient acceleration of CRs via their stochastic non-resonant interactions with magnetic compressions/expansions. The resulting mirror acceleration is dominated by the slow-mode eddies with their lifetime comparable to the mirror diffusion time of CRs. Consequently, we find that the acceleration time of mirror acceleration is independent of the spatial diffusion coefficient of CRs. The mirror acceleration brings new life for the particle acceleration in a weakly compressible/incompressible medium and has important implications for studying CR re-acceleration in the high-$\beta$ intracluster medium., Comment: 10 pages, 2 figures, submitted to ApJ

Propagation of cosmic rays (CRs) in turbulent and magnetized astrophysical media is a long-standing problem that requires both understanding of the properties of turbulent magnetic fields and their interaction with energetic particles. This review focuses on selected recent theoretical findings made based on the progress in understanding and simulating magnetohydrodynamic (MHD) turbulence. In particular, we address the problem of perpendicular and parallel propagation of CRs and identify the conditions when the perpendicular propagation is superdiffusive and diffusive. For the parallel diffusion, we discuss the problems of the traditionally used diffusion mechanism arising from pitch angle scattering and the possible solutions provided by the recently identified ``mirror diffusion" in the presence of turbulent magnetic mirrors., Comment: 21 pages, 4 figures, submitted to Frontiers in Astronomy and Space Sciences

The ubiquitous turbulence in astrophysical plasmas is important for both magnetic reconnection and reconnection acceleration. We study the particle acceleration during fast 3D turbulent reconnection with reconnection-driven turbulence. Particles bounce back and forth between the reconnection-driven inflows due to the mirror reflection and convergence of strong magnetic fields. Via successive head-on collisions, the kinetic energy of the inflows is converted into the accelerated particles. Turbulence not only regulates the inflow speed but also introduces various inflow obliquities with respect to the local turbulent magnetic fields. As both the energy gain and escape probability of particles depend on the inflow speed, the spectral index of particle energy spectrum is not universal. We find it in the range from $\approx 2.5$ to $4$, with the steepest spectrum expected at a strong guide field, i.e. a small angle between the total incoming magnetic field and the guide field. Without scattering diffusion needed for confining particles, the reconnection acceleration can be very efficient at a large inflow speed and a weak guide field., Comment: 16 pages, 8 figures, accepted for publication in ApJ

The cascade of magnetohydrodynamic (MHD) turbulence is subject to ion-neutral collisional damping and neutral viscous damping in the partially ionized interstellar medium. By examining the damping effects in the warm and partially ionized local interstellar medium, we find that the interstellar turbulence is damped by neutral viscosity at $\sim 261$ au and cannot account for the turbulent magnetic fluctuations detected by Voyager 1 and 2. The MHD turbulence measured by Voyager in the very local interstellar medium (VLISM) should be locally injected in the regime where ions are decoupled from neutrals for its cascade to survive the damping effects. With the imposed ion-neutral decoupling condition, and the strong turbulence condition for the observed Kolmogorov magnetic energy spectrum, we find that the turbulence in the VLISM is sub-Alfv\'{e}nic, and its largest possible injection scale is $\sim 194$ au., Comment: 9 pages, 2 figures, Accepted for publication in The Astrophysical Journal Letters

Delayed radio flares of optical tidal disruption events (TDEs) indicate the existence of non-relativistic outflows accompanying TDEs. The interaction of TDE outflows with the surrounding circumnuclear medium creates quasi-perpendicular shocks in the presence of toroidal magnetic fields. Because of the large shock obliquity and large outflow velocity, we find that the shock acceleration induced by TDE outflows generally leads to a steep particle energy spectrum, with the power-law index significantly larger than the "universal" index for a parallel shock. The measured synchrotron spectral indices of recently detected TDE radio flares are consistent with our theoretical expectation. It suggests that the particle acceleration at quasi-perpendicular shocks can be the general acceleration mechanism accounting for the delayed radio emission of TDEs., Comment: 8 pages, 5 figures; accepted for publication in The Astrophysical Journal

Alfv\'{e}nic component of MHD turbulence damps Alfv\'{e}nic waves. The consequences of this effect are important for many processes, from cosmic ray (CR) propagation to launching outflows and winds in galaxies and other magnetized systems. We discuss the differences in the damping of the streaming instability by turbulence and the damping of a plane parallel wave. The former takes place in the system of reference aligned with the local direction of magnetic field along which CRs stream. The latter is in the reference frame of the mean magnetic field and traditionally considered in plasma studies. We also compare the turbulent damping of streaming instability with ion-neutral collisional damping, which becomes the dominant damping effect at a sufficiently low ionization fraction. Numerical testing and astrophysical implications are also discussed., Comment: 19 pages, 3 figures, accepted for publication in Frontiers in Astronomy and Space Sciences

We study the streaming instability of GeV$-100~$GeV cosmic rays (CRs) and its damping in the turbulent interstellar medium (ISM). We find that the damping of streaming instability is dominated by ion-neutral collisional damping in weakly ionized molecular clouds, turbulent damping in the highly ionized warm medium, and nonlinear Landau damping in the Galactic halo. Only in the Galactic halo, is the streaming speed of CRs close to the Alfv\'{e}n speed. Alfv\'{e}nic turbulence plays an important role in both suppressing the streaming instability and regulating the diffusion of streaming CRs via magnetic field line tangling, with the effective mean free path of streaming CRs in the observer frame determined by the Alfv\'{e}nic scale in super-Alfv\'{e}nic turbulence. The resulting diffusion coefficient is sensitive to Alfv\'{e}n Mach number, which has a large range of values in the multi-phase ISM. Super-Alfv\'{e}nic turbulence contributes to additional confinement of streaming CRs, irrespective of the dominant damping mechanism., Comment: 11 pages, 4 figures, accepted for publication in ApJ

We investigate shock acceleration in a realistic astrophysical environment with density inhomogeneities. The turbulence induced by the interaction of the shock precursor with upstream density fluctuations amplifies both upstream and downstream magnetic fields via the turbulent dynamo. The dynamo-amplified turbulent magnetic fields (a) introduce variations of shock obliquities along the shock face, (b) enable energy gain through a combination of shock drift and diffusive processes, (c) give rise to various spectral indices of accelerated particles, (d) regulate the diffusion of particles both parallel and perpendicular to the magnetic field, and (e) increase the shock acceleration efficiency. Our results demonstrate that upstream density inhomogeneities and dynamo amplification of magnetic fields play an important role in shock acceleration, and thus shock acceleration depends on the condition of the ambient interstellar environment. The implications on understanding radio spectra of supernova remnants are also discussed., Comment: 17 pages, 12 figures, accepted for publication in ApJ

Small-scale turbulent dynamo is responsible for the amplification of magnetic fields on scales smaller than the driving scale of turbulence in diverse astrophysical media. Most earlier dynamo theories concern the kinematic regime and small-scale magnetic field amplification. Here we review our recent progress in developing the theories for the nonlinear dynamo and the dynamo regime in a partially ionized plasma. The importance of reconnection diffusion of magnetic fields is identified for both the nonlinear dynamo and magnetic field amplification during gravitational contraction. For the dynamo in a partially ionized plasma, the coupling state between neutrals and ions and the ion-neutral collisional damping can significantly affect the dynamo behavior and the resulting magnetic field structure. We present both our analytical predictions and numerical tests with a two-fluid dynamo simulation on the dynamo features in this regime. In addition, to illustrate the astrophysical implications, we discuss several examples for the applications of the dynamo theory to studying magnetic field evolution in both preshock and postshock regions of supernova remnants, in weakly magnetized molecular clouds, during the (primordial) star formation, and during the first galaxy formation., Comment: 28 pages, 12 figures, published in Reviews of Modern Plasma Physics

As the fundamental physical process with many astrophysical implications, the diffusion of cosmic rays (CRs) is determined by their interaction with magnetohydrodynamic (MHD) turbulence. We consider the magnetic mirroring effect arising from MHD turbulence on the diffusion of CRs. Due to the intrinsic superdiffusion of turbulent magnetic fields, CRs with large pitch angles that undergo mirror reflection, i.e., bouncing CRs, are not trapped between magnetic mirrors, but move diffusively along the turbulent magnetic field, leading to a new type of parallel diffusion, i.e., mirror diffusion. This mirror diffusion is in general slower than the diffusion of non-bouncing CRs with small pitch angles that undergo gyroresonant scattering. The critical pitch angle at the balance between magnetic mirroring and pitch-angle scattering is important for determining the diffusion coefficients of both bouncing and non-bouncing CRs and their scalings with the CR energy. We find non-universal energy scalings of diffusion coefficients, depending on the properties of MHD turbulence., Comment: 21 pages, 8 figures, accepted for publication in ApJ

Via amplification by turbulent dynamo, magnetic fields can be potentially important for the formation of the first stars. To examine the dynamo behavior during the gravitational collapse of primordial gas, we extend the theory of nonlinear turbulent dynamo to include the effect of gravitational compression. The relative importance between dynamo and compression varies during contraction, with the transition from dynamo- to compression-dominated amplification of magnetic fields with the increase of density. In the nonlinear stage of magnetic field amplification with the scale-by-scale energy equipartition between turbulence and magnetic fields, reconnection diffusion of magnetic fields in ideal magnetohydrodynamic (MHD) turbulence becomes important. It causes the violation of flux-freezing condition and accounts for (a) the small growth rate of nonlinear dynamo, (b) the weak dependence of magnetic energy on density during contraction, (c) the saturated magnetic energy, and (d) the large correlation length of magnetic fields. The resulting magnetic field structure and the scaling of magnetic field strength with density are radically different from the expectations of flux-freezing., Comment: 9 pages, 4 figures, accepted for publication in ApJ

Local magnetic reversals are an inseparable part of magnetohydrodynamic (MHD) turbulence whose collective outcome on an arbitrary scale in the inertial range may lead to a global stochastic reconnection event with a rate independent of small scale physics. We show that this picture is intimately related to the nanoflare theory of the solar corona. First, we argue that due to stochastic flux freezing, a generalized version of flux freezing in turbulence, the magnetic field follows the turbulent flow in a statistical sense. Bending and stretching an initially smooth field, therefore, the turbulence generally increases the magnetic spatial complexity. Strong magnetic shears associated with such a highly tangled field can trigger local reversals and field annihilations that convert magnetic energy into kinetic and thermal energy respectively. The former maintains the turbulence, which incidentally continues to entangle the field completing the cycle, while the latter enhances the heat generation in the dissipative range. We support this theoretical picture invoking recent analytical and numerical studies which suggest a correlation between magnetic complexity and magnetic energy dissipation. The amplification of multiple local, in-phase reversals by super-linear Richardson diffusion may initiate a global reconnection at larger scales, however, even in the absence of such a global stochastic reconnection, the small scale reversals will continue to interact with the turbulence. We employ conventional scaling laws of MHD turbulence to illustrate that these local events are indeed efficient in both enhancing the turbulence and generating heat. Finally, using an MHD numerical simulation, we show that the time evolution of the magnetic complexity is statistically correlated with the kinetic energy injection rate and/or magnetic-to-thermal energy conversion rate.

Astrophysical plasmas are turbulent and magnetized. The interaction between cosmic rays (CRs) and magnetohydrodynamic (MHD) turbulence is a fundamental astrophysical process. Based on the current understanding of MHD turbulence, we revisit the trapping of CRs by magnetic mirrors in the context of MHD turbulence. In compressible MHD turbulence, isotropic fast modes dominate both trapping and gyroresonant scattering of CRs. The presence of trapping significantly suppresses the pitch-angle scattering and the spatial diffusion of CRs along the magnetic field. The resulting parallel diffusion coefficient has a weaker dependence on CR energy at higher energies. In incompressible MHD turbulence, the trapping by pseudo-Alfv\'{e}n modes dominates over the gyroresonant scattering by anisotropic Alfv\'{e}n and pseudo-Alfv\'{e}n modes at all pitch angles and prevents CRs from diffusion., Comment: 8 pages, 8 figures, accepted for publication in ApJ

Magnetic reconnection, topological change in magnetic fields, is a fundamental process in magnetized plasmas. It is associated with energy release in regions of magnetic field annihilation, but this is only one facet of this process. Astrophysical flows normally have very large Reynolds numbers and are expected to be turbulent, in agreement with observations. In strong turbulence magnetic lines constantly reconnect everywhere at all scales, making magnetic reconnection an intrinsic part of turbulent cascade. We note that this is inconsistent with the usual practice of regarding magnetic lines as persistent dynamical elements. A number of theoretical, numerical, and observational studies, starting with Lazarian & Vishniac (1999), demonstrated that 3D turbulence makes magnetic reconnection fast and that these two processes are intrinsically connected. We discuss the dramatic violation of the textbook concept of magnetic flux-freezing in the presence of turbulence and demonstrate that in the presence of turbulence the plasma effects are subdominant to turbulence as far as the magnetic reconnection is concerned. This justifies an MHD-like treatment of magnetic reconnection at scales much larger than the relevant plasma scales. We discuss numerical and observational evidences supporting the turbulent reconnection model. In particular, we show that tearing reconnection is suppressed in 3D and, unlike the 2D case, the 3D reconnection induces turbulence that makes reconnection independent of resistivity. We show that turbulent reconnection dramatically affects the key astrophysical processes, e.g., star formation, turbulent dynamo, acceleration of cosmic rays. We provide criticism of the concept of "reconnection-mediated turbulence" and explain why turbulent reconnection is very different from enhanced turbulent resistivity and hyper-resistivity, and why the latter has fatal conceptual flaws., Comment: 94 pages, 39 figures

Particle acceleration is a fundamental process in many high-energy astrophysical environments and determines the spectral features of their synchrotron emission. We have studied the adiabatic stochastic acceleration (ASA) of electrons arising from the basic dynamics of magnetohydrodynamic (MHD) turbulence and found that the ASA acts to efficiently harden the injected electron energy spectrum. The dominance of the ASA at low energies and the dominance of synchrotron cooling at high energies result in a broken power-law shape of both electron energy spectrum and photon synchrotron spectrum. Furthermore, we have applied the ASA to studying the synchrotron spectra of the prompt emission of gamma-ray bursts (GRBs) and pulsar wind nebulae (PWNe). The good agreement between our theories and observations confirms that the stochastic particle acceleration is indispensable in explaining their synchrotron emission., Comment: 8 pages, 6 figures, conference proceedings paper for the 18th Annual International Astrophysics Conference (February 18-22, 2019, Pasadena, CA, USA)

The small-scale turbulent dynamo is an important process contributing to the cosmic magnetization. In partially ionized astrophysical plasmas, the dynamo growth of magnetic energy strongly depends on the coupling state between ions and neutrals and the ion-neutral collisional damping effect. A new damping stage of turbulent dynamo in a weakly ionized medium was theoretically predicted by Xu \& Lazarian (2016). By carrying out a 3D two-fluid dynamo simulation, here we for the first time numerically confirmed the physical conditions and the linear-in-time growth of magnetic field strength of the damping stage of dynamo. The dynamo-amplified magnetic field has a characteristic length as the damping scale, which increases with time and can reach the injection scale of turbulence after around eight largest eddy-turnover times given sufficiently low ionization fraction and weak initial magnetic field. Due to the weak coupling between ions and neutrals, most turbulent energy carried by neutrals cannot be converted to the magnetic energy, resulting in a relatively weak magnetic field at the end of dynamo. This result has important implications for the growth of magnetic fields in the partially ionized interstellar medium and shock acceleration of Galactic cosmic rays., Comment: Accepted for publication in ApJ; 11 pages, 5 figures

As shown by broad-band observations, pulsar wind nebulae (PWNe) are characterized by a broken power-law spectrum of synchrotron emission. Based on the modern magnetohydrodynamic (MHD) turbulence theories, we investigate the re-acceleration of electrons in the PWN through the adiabatic stochastic acceleration (ASA), which arises from fundamental dynamics of MHD turbulence. The ASA acts to flatten the injected energy spectrum of electrons at low energies, while synchrotron cooling results in a steep spectrum of electrons at high energies. Their dominance in different energy ranges leads to a flat radio spectrum ($F_\nu$) and a steep X-ray spectrum. Our analytical spectral shapes generally agree well with the observed synchrotron spectra of radio- and X-ray-bright PWNe. The spectral break corresponding to the balance between the ASA and synchrotron losses provides a constraint on the acceleration timescale of the ASA and the magnetic field strength in the PWN., Comment: Accepted for publication in ApJ, 11 pages, 5 figures

Astrophysical fluids are turbulent, magnetized and frequently partially ionized. As an example of astrophysical turbulence, the interstellar turbulence extends over a remarkably large range of spatial scales and participates in key astrophysical processes happening in different ranges of scales. A significant progress has been achieved in the understanding of the magnetohydrodynamic (MHD) turbulence since the turn of the century, and this enables us to better describe turbulence in magnetized and partially ionized plasmas. In fact, the modern revolutionized picture of the MHD turbulence physics facilitates the development of various theoretical domains, including the damping process for dissipating MHD turbulence and the dynamo process for generating MHD turbulence with many important astrophysical implications. In this paper, we review some important findings from our recent theoretical works to demonstrate the interconnection between the properties of MHD turbulence and those of turbulent dynamo in a partially ionized gas. We also briefly exemplify some new tentative studies on how the revised basic processes influence the associated outstanding astrophysical problems in, such as, magnetic reconnection, cosmic ray scattering, magnetic field amplification in both the early and the present-day universe., Comment: 25 pages, 9 figures, published in New Journal of Physics

As a fundamental astrophysical process, the scattering of particles by turbulent magnetic fields has its physical foundation laid by the magnetohydrodynamic (MHD) turbulence theory. In the framework of the modern theory of MHD turbulence, we derive a generalized broadened resonance function by taking into account both the magnetic fluctuations and nonlinear decorrelation of turbulent magnetic fields arising in MHD turbulence, and we specify the energy range of particles for the dominance of different broadening mechanisms. The broadened resonance allows for scattering of particles beyond the energy threshold of the linear resonance. By analytically determining the pitch-angle diffusion coefficients for transit time damping (TTD) with slow and fast modes, we demonstrate that the turbulence anisotropy of slow modes suppresses their scattering efficiency. Furthermore, we quantify the dependence of the relative importance between slow and fast modes in TTD scattering on (i) particle energy, (ii) plasma $\beta$ (the ratio of gas pressure to magnetic pressure), and (iii) damping of MHD turbulence, and we also provide the parameter space for the dominance of slow modes. To exemplify its applications, we find that among typical partially ionized interstellar phases, in the warm neutral medium slow and fast modes have comparable efficiencies in TTD scattering of cosmic rays. For low-energy particles, e.g., sub-Alfv\'{e}nic charged grains, we show that slow modes always dominate TTD scattering., Comment: 14 pages, 3 figures, accepted for publication in ApJ

We introduce a non-resonant acceleration mechanism arising from the second adiabatic invariant in magnetic turbulence and apply it to study the prompt emission spectra of gamma-ray bursts (GRBs). The mechanism contains both the first- and second-order Fermi acceleration, originating from the interacting turbulent reconnection and dynamo processes. It leads to a hard electron energy distribution up to a cutoff energy at the balance between the acceleration and synchrotron cooling. The sufficient acceleration rate ensures a rapid hardening of any initial energy distribution to a power-law distribution with the index $p \sim 1$, which naturally produces a low-energy photon index $\alpha \sim -1$ via the synchrotron radiation. For typical GRB parameters, the synchrotron emission can extend to a characteristic photon energy on the order of $\sim 100$ keV., Comment: 5 pages, 2 figures; accepted for publication in ApJL

By following the Kazantsev theory and taking into account both microscopic and turbulent diffusion of magnetic fields, we develop a unified treatment of the kinematic and nonlinear stages of turbulent dynamo, and study the dynamo process for a full range of magnetic Prandtl number Pm and ionization fractions. We find a striking similarity between the dependence of dynamo behavior on Pm in a conducting fluid and R (a function of ionization fraction) in partially ionized gas. In a weakly ionized medium, the kinematic stage is largely extended, including not only exponential growth but a new regime of dynamo characterized by linear-in-time growth of magnetic field strength, and the resulting magnetic energy is much higher than the kinetic energy carried by viscous-scale eddies. Unlike the kinematic stage, the subsequent nonlinear stage is unaffected by microscopic diffusion processes and has a universal linear-in-time growth of magnetic energy with the growth rate as a constant fraction $3/38$ of the turbulent energy transfer rate, showing good agreement with earlier numerical results. Applying the analysis to the first stars and galaxies, we find that the kinematic stage is able to generate a field strength only an order of magnitude smaller than the final saturation value. But the generation of large-scale magnetic fields can only be accounted for by the relatively inefficient nonlinear stage and requires longer time than the free-fall time. It suggests that magnetic fields may not have played a dynamically important role during the formation of the first stars., Comment: 21 pages, 2 figures, accepted for publication in ApJ

We study the damping from neutral-ion collisions of both incompressible and compressible magnetohydrodynamic (MHD) turbulence in partially ionized medium. We start from the linear analysis of MHD waves applying both single-fluid and two-fluid treatments. The damping rates derived from the linear analysis are then used in determining the damping scales of MHD turbulence. The physical connection between the damping scale of MHD turbulence and cutoff boundary of linear MHD waves is investigated. Our analytical results are shown to be applicable in a variety of partially ionized interstellar medium (ISM) phases and solar chromosphere. As a significant astrophysical utility, we introduce damping effects to propagation of cosmic rays in partially ionized ISM. The important role of turbulence damping in both transit-time damping and gyroresonance is identified., Comment: 29 pages, 16 figures