Basic Properties of Compressible MHD Turbulence: Implications for Molecular Clouds

Recent advances in understanding of the basic properties of compressible Magnetohydrodynamic (MHD) turbulence call for revisions of some of the generally accepted concepts. First, the MHD turbulence is not so messy as it is usually believed. In fact, the notion of strong nonlinear coupling of compressible and incompressible motions 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 viscosity by neutrals does not suppress MHD turbulence in a partially ionized gas. Instead, MHD turbulence develops magnetic cascade at scales below the scale at which neutrals damp ordinary hydrodynamic motions. The implications of those changes of MHD turbulence paradigm for molecular clouds require further studies. Those studies can benefit from testing of theoretical predictions using new statistical techniques that utilize spectroscopic data. We briefly discuss advances in development of tools using which the statistics of turbulent velocity can be recovered from observations.     

Incompressible MHD Turbulence

The inertial range of incompressible MHD turbulence is most conveniently described in terms of counter propagating waves. Shear Alfvén waves control the cascade dynamics. Slow waves play a passive role and adopt the spectrum set by the shear Alfvén waves. Cascades composed entirely of shear Alfvén waves do not generate a significant measure of slow waves. MHD turbulence is anisotropic with energy cascading more rapidly along k _⊥ than along k _∥. Anisotropy increases with k _⊥ such that the excited modes are confined inside a cone bounded by k _∥∝ k _perp ^2/3. The opening angle of the cone, θ(k _⊥)∝ k _⊥ ^-1/3, defines the scale dependent anisotropy. MHD turbulence is generically strong in the sense that the waves which comprise it are critically damped. Nevertheless, deep inside the inertial range, turbulent fluctuations are small. Their energy density is less than that of the background field by a factor θ²(k _⊥)≪. MHD cascades are best understood geometrically. Wave packets suffer distortions as they move along magnetic field lines perturbed by counter propagating wave packets. Field lines perturbed by unidirectional waves map planes perpendicular to the local field into each other. Shear Alfvén waves are responsible for the mapping's shear and slow waves for its dilatation. The former exceeds the latter by θ^-1(k _⊥)≫ 1 which accounts for dominance of the shear Alfvén waves in controlling the cascade dynamics. This revised version was published online in July 2006 with corrections to the Cover Date.     

Acceleration of particles by electron plasma waves in a moderate magnetic field

A general scheme is established to examine any magnetohydrodynamic (MHD) configuration for its acceleration potential including the effects of various types of plasma waves. The analysis is restricted to plasma waves in a magnetic field with electron cyclotron frequency less than, but comparable to, the electron plasma frequency (moderate field). The general role of electron plasma waves is examined in this paper independent of a specific MHD configuration or generating mechanism in the weak turbulence limit. The evolution of arbitrary wave spectra in a non-relativistic plasma is examined, and it is shown that the nonlinear, process of induced scattering on the polarization clouds of ions leads to the collapse of the waves to an almost one-dimensional spectrum directed along the magnetic field. The subsequent acceleration of non-relativistic and relativistic particles is considered. It is shown for non-relativistic particles that when the wave distribution has a negative slope the acceleration is retarded for lower velocities and enhanced for higher velocities compared to acceleration by an isotropic distribution of electron plasma waves in a magnetic field. This change in behavior is expected to affect the development of wave spectra and the subsequent acceleration spectrum.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

HALL-MHD turbulence in the Solar atmosphere

When the velocity field of a magneto-fluid is accorded an intrinsically equal status with the magnetic field, standard magnetohydrodynamics (MHD) must be replaced by the dispersive or Hall MHD which retains some crucial two-fluid effects, in particular the physics on the ion skin depth scale. The larger system has three quadratic invariants: the generalized helicity (a sum of the cross and fluid helicities) is added to the ranks of the standard total energy and the magnetic helicity invariants. Based on this extended set, dimensional arguments à la Kolmogorov are invoked to derive the turbulent spectral distributions of the kinetic and the magnetic energy densities in the inertial range. By using the selective dissipation hypothesis, we construct the spectra on three different scales within the inertial range, and acknowledge the possible role of the dual cascade of the kinetic and the magnetic energy densities. The additional structure imparted to the spectral laws (by the inclusion of the generalized helicity) allows us to reproduce, remarkably well, the essentials as well as details of the observed spectra of the motions and of the magnetic fields of the solar atmosphere on the scales of a few thousand km.     

Spontaneous magnetic reconnection

The present review concerns the relevance of collisionless reconnection in the astrophysical context. Emphasis is put on recent developments in theory obtained from collisionless numerical simulations in two and three dimensions. It is stressed that magnetic reconnection is a universal process of particular importance under collisionless conditions, when both collisional and anomalous dissipation are irrelevant. While collisional (resistive) reconnection is a slow, diffusive process, collisionless reconnection is spontaneous. On any astrophysical time scale, it is explosive. It sets on when electric current widths become comparable to the leptonic inertial length in the so-called lepton (electron/positron) “diffusion region”, where leptons de-magnetise. Here, the magnetic field contacts its oppositely directed partner and annihilates. Spontaneous reconnection breaks the original magnetic symmetry, violently releases the stored free energy of the electric current, and causes plasma heating and particle acceleration. Ultimately, the released energy is provided by mechanical motion of either the two colliding magnetised plasmas that generate the current sheet or the internal turbulence cascading down to lepton-scale current filaments. Spontaneous reconnection in such extended current sheets that separate two colliding plasmas results in the generation of many reconnection sites (tearing modes) distributed over the current surface, each consisting of lepton exhausts and jets which are separated by plasmoids. Volume-filling factors of reconnection sites are estimated to be as large as \({<}10^{-5}\) per current sheet. Lepton currents inside exhausts may be strong enough to excite Buneman and, for large thermal pressure anisotropy, also Weibel instabilities. They bifurcate and break off into many small-scale current filaments and magnetic flux ropes exhibiting turbulent magnetic power spectra of very flat power-law shape \(W_b\propto k^{-\alpha }\) in wavenumber k with power becoming as low as \(\alpha \approx 2\). Spontaneous reconnection generates small-scale turbulence. Imposed external turbulence tends to temporarily increase the reconnection rate. Reconnecting ultra-relativistic current sheets decay into large numbers of magnetic flux ropes composed of chains of plasmoids and lepton exhausts. They form highly structured current surfaces, “current carpets”. By including synchrotron radiation losses, one favours tearing-mode reconnection over the drift-kink deformation of the current sheet. Lepton acceleration occurs in the reconnection-electric field in multiple encounters with the exhausts and plasmoids. This is a Fermi-like process. It results in power-law tails on the lepton energy distribution. This effect becomes pronounced in ultra-relativistic reconnection where it yields extremely hard lepton power-law energy spectra approaching \(F(\gamma )\propto \gamma ^{-1}\), with \(\gamma \) the lepton energy. The synchrotron radiation limit becomes substantially exceeded. Relativistic reconnection is a probable generator of current and magnetic turbulence, and a mechanism that produces high-energy radiation. It is also identified as the ultimate dissipation mechanism of the mechanical energy in collisionless magnetohydrodynamic turbulent cascades via lepton-inertial-scale turbulent current filaments. In this case, the volume-filling factor is large. Magnetic turbulence causes strong plasma heating of the entire turbulent volume and violent acceleration via spontaneous lepton-scale reconnection. This may lead to high-energy particle populations filling the whole volume. In this case, it causes non-thermal radiation spectra that span the entire interval from radio waves to gamma rays.     

Whistler wave cascades in solar wind plasma

Non-linear, three-dimensional, time-dependent fluid simulations of whistler wave turbulence are performed to investigate role of whistler waves in solar wind plasma turbulence in which characteristic turbulent fluctuations are characterized typically by the frequency and length-scales that are, respectively, bigger than ion gyrofrequency and smaller than ion gyroradius. The electron inertial length is an intrinsic length-scale in whistler wave turbulence that distinguishably divides the high-frequency solar wind turbulent spectra into scales smaller and bigger than the electron inertial length. Our simulations find that the dispersive whistler modes evolve entirely differently in the two regimes. While the dispersive whistler wave effects are stronger in the large-scale regime, they do not influence the spectral cascades which are describable by a Kolmogorov-like   k ^−7/3  spectrum. By contrast, the small-scale turbulent fluctuations exhibit a Navier–Stokes-like evolution where characteristic turbulent eddies exhibit a typical   k ^−5/3  hydrodynamic turbulent spectrum. By virtue of equipartition between the wave velocity and magnetic fields, we quantify the role of whistler waves in the solar wind plasma fluctuations.     

Diffusion and dispersion in magnetohydrodynamic turbulence: The influence of mean magnetic fields

The properties of magnetohydrodynamic (MHD) turbulence under the influence of a strong mean magnetic field are investigated from the Lagrangian viewpoint by tracking fluid particles in direct numerical simulations. The particle trajectories show characteristic bends near vortex sheets. A strong mean magnetic field leads to preferential diffusion parallel to the mean magnetic field. The two‐particle relative dispersion process shows a dependence on the orientation of the initial separation vector. The relative dispersion is slowed down for initial separation vectors aligned with the mean magnetic field. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Proton acceleration in flares and magnetohydrodynamic solar activity

The possibility of accelerated protons in solar flares having a sharp change in their spectral index is discussed. The analysis is based on the Tsytovich (1982, 1984, 1987a, b, c) acceleration model by MHD turbulence, which is shown to have different resonant conditions for non-relativistic and relativistic particles. The different resonant condition is shown to result in a sharp change in the accelerated proton spectral index, even in the absence of any peculiarity in the spectra of the MHD turbulence. Time scales for accelerated protons to relativistic energies are also derived, and shown to be consistent with observations. We also show that the threshold energy for electron acceleration by low frequency MHD turbulence is much greater than for proton acceleration. The turbulence therefore preferentially accelerates protons.     

Nonplanar magnetic reconnection associated with conical slow shocks

We study a nonplanar model of magnetic reconnection associated with conical slow shocks, assuming that the shock surfaces are two identical cones with circular cross sections symmetrical about the ±x-axis. In the inflow region upstream of the shocks, two oppositely directed magnetic fields are separated by a current sheet. The model treats the current sheet as a tangential discontinuity and treats shocks and tangential discontinuity as surfaces of zero thickness. The dynamical structure of the global magnetic field in the continuous regions is studied using compressible, non-resistive MHD equations. In the inflow region, nonplanar magnetic field lines first move toward the current sheet. Near the sheet, the middle sections of the field lines become highly flattened, almost parallel to the sheet. Eventually, then oppositely directed field lines merge across the tangential discontinuity between the two shocks, and the magnetic lines are reconnected at the intersection of the shock and the tangential discontinuity. Reconnected magnetic lines are carried away at high speeds by the MHD flow in the outflow region, downstream of the shocks.     

A coronal magnetic field model with horizontal volume and sheet currents

When globally mapping the observed photospheric magnetic field into the corona, the interaction of the solar wind and magnetic field has been treated either by imposing source surface boundary conditions that tacitly require volume currents outside the source surface (Schatten, Wilcox, and Ness, 1969) or by limiting the interaction to thin current sheets between oppositely directed field regions (Wolfson, 1985). Yet observations and numerical MHD calculations suggest the presence of non-force-free volume currents throughout the corona as well as thin current sheets in the neighborhoods of the interfaces between closed and open field lines or between oppositely directed open field lines surrounding coronal helmet-streamer structures. This work presents a model including both horizontal volume currents and streamer sheet currents. The present model builds on the magnetostatic equilibria developed by Bogdan and Low (1986) and the current-sheet modeling technique developed by Schatten (1971). The calculation uses synoptic charts of the line-of-sight component of the photospheric magnetic field measured at the Wilcox Solar Observatory. Comparison of an MHD model with the calculated model results for the case of a dipole field and comparison of eclipse observations with calculations for CR 1647 (near solar minimum) show that this horizontal current-current-sheet model reproduces polar plumes and axes of corona streamers better than the source-surface model and reproduces coronal helmet structures better than the current-sheet model.     

Physical conditions in radio jets and re-acceleration of relativistic electrons by MHD turbulence

Based on the observed features and physical conditions in the radio jets of several low-luminosity radio galaxies, I discuss the re-acceleration of relativistic electrons. On assuming a Fermi type acceleration, an acceleration coefficient of ～10^?15 s^?1 was obtained, which can well explain the radio brightness distribution in the jets and their spectrum. I further discuss the possibility of MHD turbulence providing the acceleration, and find that the turbulence energy spectral index must be restricted to the very narrow range 1.6–1.7.     

QPO as the Rosetta Stone for understanding black hole accretion

Quasi‐periodic oscillations (QPO) seen in the X‐ray fluxes of individual neutron stars and black hole sources are one of most intriguing phenomena in today's astrophysics. The QPO nature is visibly determined by super‐strong Einstein's gravity. I argue here that it also profoundly depends on the MRI turbulence in accretion flows. Understanding the QPO physics may therefore guide accretion theory out of its present state of confusion. Readers will find here an up‐to‐date, comprehensible account of what is known, and what is not, about the QPO physics. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Analytical theory of tachocline transport at mid latitudes

We provide a theory of magnetic diffusion, momentum transport, and mixing in the solar tachocline by considering magnetohydrodynamics (MHD) turbulence on a β plane subject to a large scale shear (provided by the latitudinal differential rotation). In the strong magnetic field regime, we find that the turbulent viscosity and diffusivity are reduced by magnetic fields only, similarly to the two-dimensional MHD case (without Rossby waves). In the weak magnetic field regime, we find a crossover scale (L_R) from a Alfvén dominated regime (on small scales) to a Rossby dominated regime (on large scales). For parameter values typical of the tachocline, L_R is larger than the solar radius so that Rossby waves are unlikely to play an important role in the transport of magnetic field and angular momentum. This is mainly due to the enhancement of magnetic back-reaction by shearing which efficiently generates small scales, thus strong currents. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Energetic particle acceleration in a three-dimensional magnetic field reconnection model: the role of magnetohydrodynamic turbulence

The role of magnetohydrodynamic (MHD) turbulence in the cosmic ray acceleration process in a volume with a reconnecting magnetic field is studied by means of Monte Carlo simulations. We performed modelling of proton acceleration, with the three-dimensional analytic model of stationary reconnection of Craig et al. providing the unperturbed background conditions. Perturbations of particle trajectories resulting from a turbulent magnetic field component were simulated using small-amplitude pitch-angle momentum scattering, enabling modelling of both small- and large-amplitude turbulence in a wide wavevector range. Within the approach, no second-order Fermi acceleration process is allowed. Comparison of the acceleration process in models involving particle trajectory perturbations with the unperturbed model reveals that the turbulence can substantially increase the acceleration efficiency, enabling much higher final particle energies and flat particle spectra.     

Narrowband dm-spikes in the 2 GHz frequency range and MHD cascading waves in reconnection outflows

Two examples of clouds of narrowband dm-spikes, observed by the Ondejov radiospectrograph in the 1–2 GHz frequency range, are analyzed. After transformation of the frequency scales to distances in the solar atmosphere, the power spectra analysis of size scales reveals a spectral index of –5/3, resembling that of Kolmogorov spectra of turbulent cascades. The narrowband dm-spikes are interpreted as radio emission from electrons accelerated in MHD cascading waves, probably generated in plasma outflows from magnetic field reconnection.     

The Hall effect as the producer of current sheets in astrophysical plasmas

Our solution of the MHD equations with the Hall effect shows that this effect can produce thin current sheets in stellar atmospheres at heights where the plasma is tenuous and the Hall effect can profoundly influence the magnetic field variations. The current in the sheets is directed oppositely to the local plasma density gradient. In partiuclar, such a phenomenon is possible on the Sun near the base of the corona.     

Can turbulent plane‐shear flows support large‐scale dynamos?

Turbulent plane‐shear flow is found to show same basic effects of mean‐fieldMHD as rotating turbulence. In particular, the mean electromotive force (EMF) includes highly anisotropic turbulent diffusion and alpha‐effect. Only magnetic diffusion remains for spatially‐uniform turbulence. The question is addressed whether in this case a self‐excitation of a magnetic field by so‐called sher‐current dynamo is possible and the quasilinear theory provides a negative answer. The streamaligned component of the EMF has the sign opposite to that required for dynamo. If, however, the turbulence is not uniform across the flow direction then a dynamo‐active α ‐effect emerges. The critical magnetic Reynolds number for the alpha‐shear dynamo is estimated to be slightly above ten. Possibilities for cross‐checking theoretical predictions with MHD experiments are discussed. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Compressible turbulence in Hall Magnetohydrodynamics

By direct numerical simulations we investigate the nonlinear dynamics of a compressible Hall Magnetohydrodynamic (MHD) plasma. At small scales, where the Hall effect dominates, we found an increase of the compressibility of the system and the breakdown of the strong link between velocity and magnetic fields, typical of usual MHD. Moreover, we find that small-scale fluctuations are characterized by an anti-correlation between density and magnetic field intensity. These features characterize the excitation of a quasi-perpendicular magnetosonic turbulence that can be interpreted as the small-scale signature of the break-down of the MHD nonlinear energy cascade due to Hall effect. Fluctuations with the same properties, based on measurements by Cluster spacecraft in space plasma turbulence during different magnetopause crossings, have been recently observed.     

Do mean‐field dynamos in nonrotating turbulent shear‐flows exist?

A plane‐shear flow in a fluid with forced turbulence is considered. If the fluid is electrically‐conducting then a mean electromotive force (EMF) results even without basic rotation and the magnetic diffusivity becomes a highly anisotropic tensor. It is checked whether in this case self‐excitation of a large‐scale magnetic field is possible (so‐called W̄ × J̄ ‐dynamo) and the answer is NO. The calculations reveal the cross‐stream components of the EMF perpendicular to the mean current having the wrong signs, at least for small magnetic Prandtl numbers. After our results numerical simulations with magnetic Prandtl number of about unity have only a restricted meaning as the Prandtl number dependence of the diffusivity tensor is rather strong. If, on the other hand, the turbulence field is strati.ed in the vertical direction then a dynamo‐active α ‐effect is produced. The critical magnetic Reynolds number for such a self‐excitation in a simple shear flow is slightly above 10 like for the other – but much more complicated – flow patterns used in existing dynamo experiments with liquid sodium or gallium. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)