This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.

We review the current status of research in MHD turbulence theory and numerical experiments and their applications to astrophysics and solar science. We introduce general tools for studying turbulence, basic turbulence models, MHD equations and their wave modes. Subsequently, we cover the theories and numerics of Alfvenic turbulence, imbalanced turbulence, small-scale dynamos and models and numerics for supersonic MHD turbulence.

Magnetic reconnection is best known from observations of the Sun where it causes solar flares. Observations estimate the reconnection rate a small, but non-negligible fraction of the Alfv\'en speed, so-called fast reconnection. Until recently, the prevailing pictures of reconnection were referring to either resistivity or plasma microscopic effects, which was contradictory to the observed rates. The alternative picture was either reconnection due to the stochasticity of magnetic field lines in turbulence or the tearing instability of the thin current sheet. In this paper I simulated long-term three-dimensional nonlinear evolution of a thin, planar current sheet subject to fast oblique tearing instability using direct numerical simulations of resistive-viscous MHD. The late-time evolution resembles generic turbulence with -5/3 power spectrum and scale-dependent anisotropy, so I conclude that the tearing-driven reconnection becomes turbulent reconnection. The turbulence is local in scale, so microscopic diffusivity should not affect large-scale quantities. This is confirmed by convergence of the reconnection rate towards $\sim 0.015 v_A$ with increasing Lundquist number. In this spontaneous reconnection with mean field and without driving the dissipation rate per unit area also converge to $\sim 0.006 \rho v_A^3$, the dimensionless constants $0.015$ and $0.006$ are governed only by self-driven nonlinear dynamics of the sheared magnetic field. Remarkably, this also means that thin current sheet has a universal fluid resistance depending only on its length to width ratio and to $v_A/c$.

We demonstrate that particles are regularly accelerated while experiencing curvature drift in flows driven by magnetic tension. Some examples of such flows include spontaneous turbulent reconnection and decaying magnetohydrodynamic (MHD) turbulence, where magnetic field relaxes to a lower-energy configuration and transfers part of its energy to kinetic motions. The opposite process, such as dynamo, will actually result in the net cooling of particles by the curvature drift. Being very generic, this acceleration mechanism is likely to be responsible in production of non-thermal particle distribution in magnetically-dominant environments such as solar chromosphere, pulsar magnetosphere, jets from supermassive black holes, $\gamma$-ray bursts, etc.

Massive galaxy clusters are filled with a hot, turbulent and magnetized intra-cluster medium. Still forming under the action of gravitational instability, they grow in mass by accretion of supersonic flows. These flows partially dissipate into heat through a complex network of large-scale shocks [1], while residual transonic flows create giant turbulent eddies and cascades [2,3]. Turbulence heats the intra-cluster medium [4] and also amplifies magnetic energy by way of dynamo action [5-8]. However, the pattern regulating the transformation of gravitational energy into kinetic, thermal, turbulent and magnetic energies remains unknown. Here we report that the energy components of the intra-cluster medium are ordered according to a permanent hierarchy, in which the ratio of thermal to turbulent to magnetic energy densities remains virtually unaltered throughout the cluster's history, despite evolution of each individual component and the drive towards equipartition of the turbulent dynamo. This result revolves around the approximately constant efficiency of turbulence generation from the gravitational energy that is freed during mass accretion, revealed by our computational model of cosmological structure formation [3,9]. The permanent character of this hierarchy reflects yet another type of self-similarity in cosmology [10-13], while its structure, consistent with current data [14-18], encodes information about the efficiency of turbulent heating and dynamo action.

We compare DNS calculations of homogeneous isotropic turbulence with the statistical properties of intra-cluster turbulence from the Matryoshka Run (Miniati 2014) and find remarkable similarities between their inertial ranges. This allowed us to use the time dependent statistical properties of intra-cluster turbulence to evaluate dynamo action in the intra-cluster medium, based on earlier results from numerically resolved nonlinear magneto-hydrodynamic turbulent dynamo (Beresnyak 2012). We argue that this approach is necessary (a) to properly normalize dynamo action to the available intra-cluster turbulent energy and (b) to overcome the limitations of low Re affecting current numerical models of the intra-cluster medium. We find that while the properties of intra-cluster magnetic field are largely insensitive to the value and origin of the seed field, the resulting values for the Alfven speed and the outer scale of the magnetic field are consistent with current observational estimates, basically confirming the idea that magnetic field in today's galaxy clusters is a record of its past turbulent activity.

Anisotropy of MHD turbulence has been studied extensively for many years, most prominently by measurements in the solar wind and high resolution simulations. The spectrum parallel to the local magnetic field was observed to be steeper than perpendicular spectrum, typically $k^{-2}$, consistent with the widely accepted Goldreich-Sridhar (1995) model. In this Letter I looked deeper into the nature of the relation between parallel and perpendicular spectra and argue that this $k^{-2}$ scaling has the same origin as the $\omega^{-2}$ scaling of Lagrangian frequency spectrum in strong hydrodynamic turbulence. This follows from the fact that Alfv\'en waves propagate along magnetic field lines. It is now became clear that the observed anisotropy can be argued without invocation of the "critical balance" argument and is more robust that was previously thought. The relation between parallel (Lagrangian) and perpendicular (Eulerian) spectra is inevitable consequence of strong turbulence of Alfven waves, rather than a conjecture based on the uncertainty relation. I tested this using high-resolution simulations of MHD turbulence, in particular I verified that the cutoff of the parallel spectrum scales as Kolmogorov timescale, not lengthscale.

Cosmic strings are linear topological defects which are hypothesized to be produced during inflation. Most searches for strings have been relying on the string's lensing of background galaxies or CMB. In this paper I obtained the solution for the supersonic flow of the collisional gas past the cosmic string which has two planar shocks with shock compression ratio that depend on the angle defect of the string and its speed. The shocks result in compression and heating of the gas and, given favorable condition, particle acceleration. The gas heating and overdensity in an unusual wedge shape can be detected by observing HI line at high redshifts. The particle acceleration can occur in present-day Universe when the string crosses the hot gas contained in galaxy clusters and, since the consequences of such collision persist for cosmological timescales, could be located by looking at the unusual large-scale radio sources situated on a single spatial plane.

In a Comment by Perez et al (2014) it is claimed that recently published simulations of Beresnyak (2014) are grossly underresolved, compared to theirs, and that Beresnyak (2014) failed to estimate numerical error. Both claims are contrary to the fact. Firstly, as far as numerical resolution is concerned, Beresnyak (2014) was using $k_{\rm max}\eta_{41}$=1.06 resolution criterion, while Perez et al has been using $k_{\rm max}\eta_{41}$=0.8. Obviously, $1.06>0.8$. Secondly, Beresnyak (2014) have estimated numerical error and found it to be below $3\times 10^{-3}$, which is properly explained in the paper. On close inspection of the Comment I have not found a single numerical value or parameter pertaining to the criticized paper, Beresnyak (2014), and it is completely unclear how the Authors came to their conclusions.

Nonresonant current instability was identified by Bell (2004) as an important mechanism for magnetic field amplification in supernova remnants. In this paper we focus on studying the nonlinear stage of this instability using the incompressible MHD formulation. We demonstrate that the evolution of magnetic turbulence driven by the Bell instability resembles turbulence driven on large scales. More importantly, we demonstrate that the energy-containing scale for magnetic fields is proportional to the square root of magnetic energy density. Given the observational constraints of the possible field amplification, this new relation allows us to directly estimate the maximum energy of particles scattered by such fields and this estimate is normally below the average particle energy. This implies that, without taking into account the feedback to cosmic rays, the typical scales of Bell fields, in either linear or nonlinear regime, will be too small to affect high energy particle acceleration. We mention several scenarios of back-reaction to cosmic rays that could be important.

MHD Turbulence is common in many space physics and astrophysics environments. We first discuss the properties of incompressible MHD turbulence. A well-conductive fluid amplifies initial magnetic fields in a process called small-scale dynamo. Below equipartition scale for kinetic and magnetic energies the spectrum is steep (Kolmogorov -5/3) and is represented by critically balanced strong MHD turbulence. In this paper we report the basic reasoning behind universal nonlinear small-scale dynamo and the inertial range of MHD turbulence. We measured the efficiency of the small-scale dynamo $C_E=0.05$, Kolmogorov constant $C_K=4.2$ and anisotropy constant $C_A=0.63$ for MHD turbulence in high-resolution direct numerical simulations. We also discuss so-called imbalanced or cross-helical MHD turbulence which is relevant for in many objects, most prominently in the solar wind. We show that properties of incompressible MHD turbulence are similar to the properties of Alfv\'enic part of MHD cascade in compressible turbulence. The other parts of the cascade evolve according to their own dynamics. The slow modes are being cascaded by Alfv\'enic modes, while fast modes create an independent cascade. We show that different ways of decomposing compressible MHD turbulence into Alfv\'en, slow and fast modes provide consistent results and are useful in understanding not only turbulent cascade, but its interaction with fast particles.

Magnetohydrodynamic turbulence is a ubiquitous phenomenon in solar physics, plasma physics and astrophysics and governs many properties of the flows of well-conductive fluids. Recently, conflicting spectral slopes for the inertial range of MHD turbulence has been reported by different groups. Varying spectral shapes from earlier simulations hinted at a wider spectral locality of MHD, which necessitated higher resolution simulations and careful and rigorous numerical analysis. In this Letter we present two groups of simulations with resolution up to 4096^3 that are numerically well-resolved and has been analyzed with exact and well-tested method of scaling study. Our results from both simulation groups indicate that the power spectral slope for all energy-related quantities, such as total energy and residual energy are around -1.7, close to Kolmogorov's -5/3. This suggests that residual energy is the constant fraction, 0.15+-0.03 of the total energy and that in asymptotic regime magnetic and kinetic spectra have the same scaling. The -1.5 slope for energy and -2 slope for residual energy suggested by other groups seems to be completely inconsistent with numerics.

Stochasticity of magnetic field lines is important for particle transport properties. Magnetic field lines separate faster than diffusively in turbulent plasma, which is called superdiffusion. We discovered that this superdiffusion is pronouncedly asymmetric, so that the separation of field lines along the magnetic field direction is different from the separation in the opposite direction. While the symmetry of the flow is broken by the so-called imbalance or cross-helicity, the difference between forward and backward diffusion is not directly due to imbalance, but a non-trivial consequence of both imbalance and non-reversibility of turbulence. The asymmetric diffusion perpendicular to the mean magnetic field entails a variety of new physical phenomena, such as the production of parallel particle streaming in the presence of perpendicular particle gradients. Such streaming and associated instabilities could be significant for particle transport in laboratory, space and astrophysical plasmas.

Recently Perez et al [arXiv:1209.2011] wrote on the spectral slope of MHD turbulence claiming that it is consistent with -3/2. This work contains a number of errors, factual inaccuracies and puzzling methods in the interpretation of numerics. We argue that the numerical evidence that the authors presented is actually against -3/2 slope and the assumption of universal alignment, opposite to what they claim. Perhaps the most puzzling is the Fig.8 that claims that the authors measured the inertial range length and that it is consistent with Re^2/3 scaling. At a close inspection it appears that the datapoints are not the result of a measurement, but rather were calculated by a formula.

Cosmological MHD simulations of galaxy cluster formation show a significant amplification of seed magnetic fields. We developed a novel method to decompose cluster magnetized turbulence into modes and showed that the fraction of the fast mode is fairly large, around 1/4 in terms of energy. This is larger than that was estimated before, which implies that cluster turbulence interacts with cosmic rays rather efficiently. We propose a framework to deal with electron and proton reacceleration in galaxy clusters that includes feedback on turbulence. In particular, we establish a new upper limit on proton and electron fluxes based on turbulence intensity. These findings, along with detailed modeling of reacceleration, will help to reconcile the observed giant radio haloes and the unobserved diffuse gamma-ray emission from these clusters.

Turbulence is ubiquitous in astrophysics. It radically changes many astrophysical phenomena, in particular, the propagation and acceleration of cosmic rays. We present the modern understanding of compressible magnetohydrodynamic (MHD) turbulence, in particular its decomposition into Alfv\'en, slow and fast modes, discuss the density structure of turbulent subsonic and supersonic media, as well as other relevant regimes of astrophysical turbulence. All this information is essential for understanding the energetic particle acceleration that we discuss further in the review. For instance, we show how fast and slow modes accelerate energetic particles through the second order Fermi acceleration, while density fluctuations generate magnetic fields in pre-shock regions enabling the first order Fermi acceleration of high energy cosmic rays. Very importantly, however, the first order Fermi cosmic ray acceleration is also possible in sites of magnetic reconnection. In the presence of turbulence this reconnection gets fast and we present numerical evidence supporting the predictions of the Lazarian & Vishniac (1999) model of fast reconnection. The efficiency of this process suggests that magnetic reconnection can release substantial amounts of energy in short periods of time. As the particle tracing numerical simulations show that the particles can be efficiently accelerated during the reconnection, we argue that the process of magnetic reconnection may be much more important for particle acceleration than it is currently accepted. In particular, we discuss the acceleration arising from reconnection as a possible origin of the anomalous cosmic rays measured by Voyagers as well as the origin cosmic ray excess in the direction of Heliotail.

We revisit the issue of spectral slope of MHD turbulence in the inertial range and argue that numerics favor Goldreich-Sridhar -5/3 slope rather than -3/2 slope. We also did precision measurements of anisotropy of MHD turbulence and determined the anisotropy constant C_A=0.34 of Alfvenic turbulence. Together with previously measured Kolmogorov constant C_K=4.2, or 3.3 for purely Alfvenic case, it constitutes a full description of MHD cascade in terms of spectral quantities, which is of high practical value for astrophysics.

We consider astrophysically relevant nonlinear MHD dynamo at large Reynolds numbers (Re). We argue that it is universal in a sense that magnetic energy grows at a rate which is a constant fraction C_E of the total turbulent dissipation rate. On the basis of locality bounds we claim that this "efficiency of small-scale dynamo", C_E, is a true constant for large Re and is determined only by strongly nonlinear dynamics at the equipartition scale. We measured C_E in numerical simulations and observed a value around 0.05 in highest resolution simulations. We address the issue of C_E being small, unlike Kolmogorov constant which is of order unity.

The spectral slope of strong MHD turbulence has recently been a matter of controversy. While Goldreich-Sridhar model (1995) predicts Kolmogorov's -5/3 slope of turbulence, shallower slopes were often reported by numerical studies. We argue that earlier numerics was affected by driving due to a diffuse locality of energy transfer in MHD case. Our highest-resolution simulation (3072^2x1024) has been able to reach the asymptotic -5/3 regime of the energy slope. Additionally, we found that so-called dynamic alignment, proposed in the model with -3/2 slope, saturates and therefore can not affect asymptotic slope. The observation of the asymptotic regime allowed us to measure Kolmogorov constant C_KA=3.2+-0.2 for purely Alfv\'enic turbulence and C_K=4.1+-0.3 for full MHD turbulence. These values are much higher than the hydrodynamic value of 1.64. The larger value of Kolmogorov constant is an indication of a fairly inefficient energy transfer and, as we show in this Letter, is in theoretical agreement with our observation of diffuse locality. We also explain what has been missing in numerical studies that reported shallower slopes.

The search of ways to generalize the theory of strong MHD turbulence for the case of non-zero cross-helicity (or energy imbalance) has attracted considerable interest recently. In our earlier publications we performed three-dimensional numerical simulations and showed that some of existing models are inconsistent with numerics. In this paper we focused our attention on low-imbalance limit and performed new high-resolution simulations. The results strongly suggest that in the limit of small imbalances we smoothly transition to a standard Goldreich-Sridhar (1995) balanced model. We also claim that Perez-Boldyrev (2009) model that predicts the same nonlinear timescale for both components due to so-called "dynamic alignment" strongly contradicts numerical evidence.

We study diffusion of Cosmic Rays (CRs) in turbulent magnetic fields using test particle simulations. Electromagnetic fields are produced in direct numerical MHD simulations of turbulence and used as an input for particle tracing, particle feedback on turbulence being ignored. Statistical transport coefficients from the test particle runs are compared with earlier analytical predictions. We find qualitative correspondence between them in various aspects of CR diffusion. In the incompressible case, that we consider in this paper, the dominant scattering mechanism occurs to be the non-resonant mirror interactions with the slow-mode perturbations. Perpendicular transport roughly agrees with being produced by magnetic field wandering.

We propose a model for Diffusive Shock Acceleration (DSA) in which stochastic magnetic fields in the shock precursor are generated through purely fluid mechanisms of a so-called small-scale dynamo. This contrasts with previous DSA models that considered magnetic fields amplified through cosmic ray streaming instabilities; i.e., either by way of individual particles resonant scattering in the magnetic fields, or by macroscopic electric currents associated with large-scale cosmic ray streaming. Instead, in our picture, the solenoidal velocity perturbations that are required for the dynamo to work are produced through the interactions of the pressure gradient of the cosmic ray precursor and density perturbations in the inflowing fluid. Our estimates show that this mechanism provides fast growth of magnetic field and is very generic. We argue that for supernovae shocks the mechanism is capable of generating upstream magnetic fields that are sufficiently strong for accelerating cosmic rays up to around 10^16 eV. No action of any other mechanism is necessary.

In this paper we systematically study the spectrum and structure of incompressible MHD turbulence by means of high resolution direct numerical simulations. We considered both balanced and imbalanced (cross-helical) cases and simulated sub-Alfvenic as well as trans-Alfvenic turbulence. This paper extends numerics preliminarily reported in Beresnyak & Lazarian (2008). We confirm that driven imbalanced turbulence has a stationary state even for high degrees of imbalance. Our major finding is that the structure of the dominant and subdominant Alfvenic components are notably different. Using the most robust observed quantities, such as the energy ratio, we were able to reject several existing models of strong imbalanced turbulence.

We performed a series of high-resolution (up to 1024^3) direct numerical simulations of hydro and MHD strong turbulence. We found that for simulations with normal viscosity the slopes for spectra of MHD are similar, although slightly more shallower than for hydro simulations. However, for simulations with hyperviscosity the slopes were very different, for instance, the slopes for hydro simulations showed pronounced and well-defined bottleneck effect, while the MHD slopes were relatively much less affected. We believe that this is indicative of MHD strong turbulence being less local than Kolmogorov turbulence. This calls for revision of MHD strong turbulence models that assume local "as-in-hydro case" cascading. Nonlocality of MHD turbulence casts doubt on numerical determination of the slopes with currently available (512^3--1024^3) numerical resolutions, including simulations with normal viscosity. We also measure various so-called alignment effects and discuss their influence on the turbulent cascade.

We present numerical simulations of driven magnetohydrodynamic (MHD) turbulence with weak/moderate imposed magnetic fields. The main goal is to clarify dynamics of magnetic field growth. We also investigate the effects of the imposed magnetic fields on the MHD turbulence, including, as a limit, the case of zero external field. Our findings are as follows. First, when we start off simulations with weak mean magnetic field only (or with small scale random field with zero imposed field), we observe that there is a stage at which magnetic energy density grows linearly with time. Runs with different numerical resolutions and/or different simulation parameters show consistent results for the growth rate at the linear stage. Second, we find that, when the strength of the external field increases, the equilibrium kinetic energy density drops by roughly the product of the rms velocity and the strength of the external field. The equilibrium magnetic energy density rises by roughly the same amount. Third, when the external magnetic field is not very strong (say, less than ~0.2 times the rms velocity when measured in the units of Alfven speed), the turbulence at large scales remains statistically isotropic, i.e. there is no apparent global anisotropy of order B_0/v. We discuss implications of our results on astrophysical fluids.

Astrophysical fluids, including interstellar and interplanetary medium, are magnetized and turbulent. Their appearance, evolution, and overall properties are determined by the magnetic turbulence that stirs it. We argue that examining magnetic turbulence at a fundamental level is vital to understanding many processes. A point that frequently escapes the attention of researchers is that magnetic turbulence cannot be confidently understood only using "brute force" numerical approaches. In this review we illustrate this point on a number of examples, including intermittent heating of plasma by turbulence, interactions of turbulence with cosmic rays and effects of turbulence on the rate of magnetic reconnection. We show that the properties of magnetic turbulence may vary considerably in various environments, e.g. imbalanced turbulence in solar wind differs from balanced turbulence and both of these differ from turbulence in partially ionized gas. Appealing for the necessity of more observational data on magnetic fields, we discuss a possibility of studying interplanetary turbulence using alignment of Sodium atoms in the tail of comets.

We consider stationary, forced, imbalanced, or cross-helical MHD Alfvenic turbulence where the waves traveling in one direction have higher amplitudes than the opposite waves. This paper is dedicated to so-called strong turbulence, which cannot be treated perturbatively. Our main result is that the anisotropy of the weak waves is stronger than the anisotropy of a strong waves. We propose that critical balance, which was originally conceived as a causality argument, has to be amended by what we call a propagation argument. This revised formulation of critical balance is able to handle the imbalanced case and reduces to old formulation in the balanced case. We also provide phenomenological model of energy cascading and discuss possibility of self-similar solutions in a realistic setup of driven turbulence.

We present a model for nonlinear decay of the weak wave in three-dimensional incompressible magnetohydrodynamic (MHD) turbulence. We show that the decay rate is different for parallel and perpendicular waves. We provide a general formula for arbitrarily directed waves and discuss particular limiting cases known in the literature. We test our predictions with direct numerical simulations of wave decay in three-dimensional MHD turbulence, and discuss the influence of turbulent damping on the development of linear instabilities in the interstellar medium and on other important astrophysical processes.

Density is the turbulence statistics that is most readily available from observations. Different regimes of turbulence correspond to different density spectra. For instance, the viscosity-damped regime of MHD turbulence relevant, for instance, to partially ionized gas, can be characterized by shallow and very anisotropic spectrum of density. This spectrum can result in substantial variations of the column densities. Addressing MHD turbulence in the regime when viscosity is not important over the inertial range, we demonstrate with numerical simulations that it is possible to reproduce both the observed Kolmogorov spectrum of density fluctuations observed in ionized gas by measuring scintillations and more shallow spectra that are obtained through the emission measurements. We show that in supersonic turbulence the high density peaks dominate shallow isotropic spectrum, while the small-scale underlying turbulence that fills most of the volume has the Kolmogorov spectrum and demonstrates scale-dependent anisotropy. The limitations of the spectrum in studying turbulence induce searches of alternative statistics. We demonstrate that a measure called "bispectrum" may be a promising tool. Unlike spectrum, the bispectrum preserves the information about wave phases.

We perform three-dimensional (3D) compressible MHD simulations over many dynamical times for an extended range of sonic and Alfven Mach numbers and analyze the statistics of 3D density and 2D column density, which include probability distribution functions, spectra, skewness, kurtosis, She-Leveque exponents, and genus. In order to establish the relation between the statistics of the observables, i.e. column densities, and the underlying 3D statistics of density, we analyze the effects of cloud boundaries. We define the parameter space for 3D measures to be recovered from column densities. In addition, we show that for subsonic turbulence the spectra of density fluctuations are consistent with k^-7/3 in the case of a strong magnetic field and k^-5/3 in the case of a weak magnetic field. For supersonic turbulence we confirm the earlier findings of the shallow spectra of density and Kolmogorov spectra of the logarithm of density. We find that the intermittencies of the density and velocity are very different.

We study the cascading of fast MHD modes in magnetically dominated plasma by performing one-dimensional (1D) dynamical simulations. We find that the cascading becomes more efficient as an angle between wave vector and underlying magnetic field increases and fast mode becomes more compressive. We also study imbalanced turbulence, in which wave packets propagating in one-direction have more power than those in the opposite direction. Unlike imbalanced \Alfvenic turbulence, the imbalanced fast mode turbulence shows faster cascading as the degree of imbalance increases. We find that the spectral index of the velocity and magnetic field, which are carried by the fast-mode turbulence, quickly reaches stationary value of -2. Thus we conclude that the dissipation of fast mode, at least in 1D case, happens not due to weak or strong turbulent cascading, but mostly due to nonlinear steepening. The density fluctuations, which are carried by slow-mode perturbation in the larger scale and entropy-mode perturbation in the smaller scale, depend on initial driving spectrum and a ratio of specific heat.

It is most natural to assume that mysterious Small Ionized and Neutral Structures (SINS) indiffuse ISM arise from turbulence. There are two obvious problem with such an explanation, however. First of all, it is generally believed that at the small scales turbulence should be damped. Second, turbulence with Kolmogorov spectrum cannot be the responsible for the SINS. We consider, however, effects, that provide spectral index flatter than the Kolmogorov one and allow action at very small scales. These are the shocks that arise in high Mach number turbulence and transfer of energy to small scales by instabilities in cosmic rays. Our examples indicate that the origin of SINS may be discovered through systematic studies of astrophysical turbulence.

We study the scattering of low-energy Cosmic Rays (CRs) in a turbulent, compressive MHD fluid. We show that compressible MHD modes -- fast or slow waves with wave lengths smaller than CR mean free paths induce cyclotron instability in CRs. The instability feeds the new small-scale Alfvenic wave component with wave vectors mostly along magnetic field, which is not a part of the MHD turbulence cascade. This new component gives feedback on the instability through decreasing the CR mean free path. We show that the ambient turbulence fully suppresses the instability at large scales, while wave steepening constrains the amplitude of the waves at small scales. We provide the energy spectrum of the plane-parallel Alfvenic component and calculate mean free paths of CRs as a function of their energy. We find that for the typical parameters of turbulence in the interstellar medium and in the intercluster medium the new Alfvenic component provides the scattering of the low energy CRs that exceeds the direct resonance scattering by MHD modes. This solves the problem of insufficient scattering of low-energy CRs in the turbulent interstellar or intracluster medium that was reported in the literature.

Goldreich-Sridhar model of incompressible turbulence provides an elegant approach to describing strong MHD turbulence. It relies on the fact that interacting Alfvenic waves are independent and have random polarization. However, in case of strong interaction a spontaneous local assymetry can arise. We used direct numerical simulations to show that polarization alignment occurs and it grows larger at smaller scales. This effect would lead to a shallower spectrum and stronger anisotropy. Even small changes in these two properties will have important astrophysical consequencies, e.g. for the cosmic ray physics.

We summarise basic properties of MHD turbulence. First, MHD turbulence is not so messy as it is believed. In fact, the notion of strong non-linear coupling of compressible and incompressible motions along MHD cascade is not tenable. Alfven, slow and fast modes of MHD turbulence follow their own cascades and exhibit degrees of anisotropy consistent with theoretical expectations. Second, the fast decay of turbulence is not related to the compressibility of fluid. Rates of decay of compressible and incompressible motions are very similar. Third, the properties of Alfven and slow modes are similar to their counterparts in the incompressible MHD. The properties of fast modes are similar to accoustic turbulence, which does require more studies. Fourth, the density at low Mach numbers and logarithm of density at higher Mach numbers exhibit Kolmogorov-type spectrum.

We study the statistics of density for supersonic turbulence in a medium with magnetic pressure larger than the gaseous pressure. This study is motivated by molecular cloud research. Our simulations exhibit clumpy density structures, which contrast increases with the Mach number. At 10 Machs densities of some clumps are three orders of magnitude higher than the mean density. These clumps give rise to flat and approximately isotropic density spectrum corresponding to the random distribution of clumps in space. We claim that the clumps originate from our random, isotropic turbulence driving. When the contribution from those clumps is suppressed by studying logarithm of density, the density statistics exhibit scale-dependent anisotropy consistent with the models where density structures arise from shearing by Alfv\'en waves. It is noteworthy that originally such models were advocated for the case of low-Mach, nearly incompressible turbulence.

We present results of the simulation of the intensity distribution of radio pulses from the Moon due to interaction of EeV neutrinos with lunar regolith. The radiation mechanism is of coherent \^Cerenkov radiation of the negative charge excess in the shower, known as Askar'yan effect. Several realistic observational setups with ground radio telescopes are considered. Effective detector volume is calculated using maximum-knowledge Monte Carlo code, and the possibilities to set limits on the diffuse neutrino flux are discussed.

We consider the propagation and stability of the helical electromagnetic perturbations in a relativistic force-free Pointing flux dominated jet (magnetic helix). Perturbations are found to dissipate energy in the vicinity of the resonant surfaces inside the jet, where the resonance of the perturbation with the Alfven waves occurs. The energy of the perturbations is converted into the energy of accelerated particles, which emit synchrotron radiation. The intensity and polarization of this radiation is modulated with the perturbation. Perturbations with vanishing group velocity are the most important and could be associated with the bright knots in parsec scale jets.

We calculate the optically thin synchrotron emission of fast electrons and positrons in a spiral stationary magnetic field and a radial electric field of a rotating relativistic strongly magnetized force-free jet consisting of electron-positron pair plasma. The magnetic field has a helical structure with a uniform axial component and a toroidal component that is maximal inside the jet and decreasing to zero towards the boundary of the jet. Doppler boosting and swing of the polarization angle of synchrotron emission due to the relativistic motion of the emitting volume are calculated. The distribution of the plasma velocity in the jet is consistent with the electromagnetic field structure. Two spatial distributions of fast particles are considered: uniform, and concentrated in the vicinity of the Alfven resonance surface. The latter distribution corresponds to the regular acceleration by an electromagnetic wave in the vicinity of its Alfven resonance surface inside the jet. The polarization properties of the radiation have been obtained and compared with the existing VLBI polarization measurements of parsec-scale jets in BL Lac sources and quasars. Our results give a natural explanation of the observed bimodality in the alignment between the electric field vector of the polarized radiation and the projection of the jet axis on the plane of the sky. We interpret the motion of bright knots as a phase velocity of standing spiral eigenmodes of electromagnetic perturbations in a cylindrical jet. The degree of polarization and the velocity of the observed proper motion of bright knots depend upon the angular rotational velocity of the jet. The observed polarizations and velocities of knots indicate that the magnetic field lines are bent in the direction opposite to the direction of the jet rotation.

We develop a theory of particle acceleration inside relativistic rotating electron-positron force-free jets with spiral magnetic fields. We considered perturbation of the stationary magnetic field structure and found that acceleration takes place in the regions where the Alfven resonant condition with the eigenmodes in the jet is fulfilled, i.e. where the local Alfven speed is equal to the phase speed of an eigenmode. The acceleration mechanism is regular acceleration by the electric field of the electromagnetic wave, which is the eigenmode of the force-free cylindrical jet: particles drift out of the region of the large wave amplitude near the Alfven resonant surface and gain energy. Acceleration in the strong electric field near the Alfven resonance and synchrotron losses combine to form a power-law energy spectrum of ultra-relativistic electrons and positrons with index between 2 and 3 depending upon the initial energy of the injected particles. The power law distribution ranges from 10 MeV to 1000 MeV.

We developed theory of particle acceleration inside the relativistic rotating electron-positron force-free jet with spiral magnetic field. We considered perturbation of stationary magnetic field structure and found that acceleration takes place in the regions where the Alfven resonant condition with the eigenmodes in the jet is fulfilled. Acceleration process and synchrotron losses combined together form power law energy spectrum of ultrarelativistic electrons and positrons with index between 2 and 3 depending upon initial energy of injected particles. Synchrotron emission of these electrons and positrons in spiral magnetic field of rotating force-free jet has been calculated. Polarization properties of the radiation has been obtained and compared with existing VLBI polarization measurements of parsec-scale jets in BL Lac sources and quasars. Our results give explanation of observed bimodality in alignment between electric field vector and jet axis. Degree of polarization and velocity of observed proper motion of bright knots depend upon angular rotational velocity of the jet. Thus, comparing them to each other, we can estimate angular rotational velocity in jets. We obtain that the fact, that generally in BL Lac objects electric field vector is oriented parallel to the jet axis while in quasars perpendicular to the jet axis, may be due to intrinsically larger angular rotational velocity and large winding of magnetic field in BL Lac jets than in quasar jets.