Basic Physical Mechanism of Fast Reconnection Evolution in Space Plasmas

Large dissipative events, such as solar flares and geomagnetic substorms, may result from sudden onset of fast (explosive) magnetic reconnection. Hence, it is a long-standing problem to find the physical mechanism that makes magnetic reconnection explosive; in particular, how can the fast magnetic reconnection explosively evolve in space plasmas? In this respect, we have proposed the spontaneous fast reconnection model as a nonlinear instability that grows by the positive feedback between plasma microphysics (anomalous resistivity) and macrophysics (global reconnection flow). On the basis of MHD simulations, we demonstrate for a variety of physical situations that the fast reconnection mechanism involving slow shocks in fact evolves explosively as a nonlinear instability and is sustained quasi-steadily on the nonlinear saturation phase. Also, distinct plasma processes, such as large-scale plasmoid propagation, magnetic loop development and loop-top heating, and asymmetric fast reconnection evolution, directly result from the spontaneous fast reconnection model. Obviously, MHD simulations are very useful in understanding the basic physics of explosive fast reconnection evolution in space plasmas. However, they cannot treat the details of microphysics near an X neutral point, which should be precisely studied in the coming 21st century.     

Emerging Parameter Space Map of Magnetic Reconnection in Collisional and Kinetic Regimes

In large-scale systems of interest to solar physics, there is growing evidence that magnetic reconnection involves the formation of extended current sheets which are unstable to plasmoids (secondary magnetic islands). Recent results suggest that plasmoids may play a critical role in the evolution of reconnection, and have raised fundamental questions regarding the applicability of resistive MHD to various regimes. In collisional plasmas, where the thickness of all resistive layers remain larger than the ion gyroradius, simulations results indicate that plasmoids permit reconnection to proceed much faster than the slow Sweet-Parker scaling. However, it appears these rates are still a factor of ～10× slower than observed in kinetic regimes, where the diffusion region current sheet falls below the ion gyroradius and additional physics beyond MHD becomes crucially important. Over a broad range of interesting parameters, the formation of plasmoids may naturally induce a transition into these kinetic regimes. New insights into this scenario have emerged in recent years based on a combination of linear theory, fluid simulations and fully kinetic simulations which retain a Fokker-Planck collision operator to allow a rigorous treatment of Coulomb collisions as the reconnection electric field exceeds the runaway limit. Here, we present some new results from this approach for guide field reconnection. Based upon these results, a parameter space map is constructed that summarizes the present understanding of how reconnection proceeds in various regimes.     

Notes on Magnetohydrodynamics of Magnetic Reconnection in Turbulent Media

Astrophysical fluids have very large Reynolds numbers and therefore turbulence is their natural state. Magnetic reconnection is an important process in many astrophysical plasmas, which allows restructuring of magnetic fields and conversion of stored magnetic energy into heat and kinetic energy. Turbulence is known to dramatically change different transport processes and therefore it is not unexpected that turbulence can alter the dynamics of magnetic field lines within the reconnection process. We shall review the interaction between turbulence and reconnection at different scales, showing how a state of turbulent reconnection is natural in astrophysical plasmas, with implications for a range of phenomena across astrophysics. We consider the process of magnetic reconnection that is fast in magnetohydrodynamic (MHD) limit and discuss how turbulence—both externally driven and generated in the reconnecting system—can make reconnection independent on the microphysical properties of plasmas. We will also show how relaxation theory can be used to calculate the energy dissipated in turbulent reconnecting fields. As well as heating the plasma, the energy dissipated by turbulent reconnection may cause acceleration of non-thermal particles, which is briefly discussed here.     

Reconnection Diffusion in Turbulent Fluids and Its Implications for Star Formation

Astrophysical fluids are turbulent a fact which changes the dynamics of many key processes, including magnetic reconnection. Fast reconnection of magnetic field in turbulent fluids allows the field to change its topology and connections. As a result, the traditional concept of magnetic fields being frozen into the plasma is no longer applicable. Plasma associated with a given magnetic field line at one instant is distributed along a different set of magnetic field lines at the next instant. This diffusion of plasmas and magnetic field is enabled by reconnection and therefore is termed “reconnection diffusion”. The astrophysical implications of this concept include heat transfer in plasmas, advection of heavy elements in interstellar medium, magnetic field generation etc. However, the most dramatic implications of the concept are related to the star formation process. The reason is that magnetic fields are dynamically important for most of the stages of star formation. The existing theory of star formation has been developed ignoring the possibility of reconnection diffusion. Instead, it appeals to the decoupling of mass and magnetic field arising from neutrals drifting in respect to ions entrained on magnetic field lines, i.e. through the process that is termed “ambipolar diffusion”. The predictions of ambipolar diffusion and reconnection diffusion are very different. For instance, if the ionization of media is high, ambipolar diffusion predicts that the coupling of mass and magnetic field is nearly perfect. At the same time, reconnection diffusion is independent of the ionization but depends on the scale of the turbulent eddies and on the turbulent velocities. In the paper we explain the physics of reconnection diffusion both from macroscopic and microscopic points of view, i.e. appealing to the reconnection of flux tubes and to the diffusion of magnetic field lines. We make use of the Lazarian and Vishniac (Astrophys. J. 517:700, 1999) theory of magnetic reconnection and show that this theory is applicable to the partially ionized gas. We quantify the reconnection diffusion rate both for weak and strong MHD turbulence and address the problem of reconnection diffusion acting together with ambipolar diffusion. In addition, we provide a criterion for correctly representing the magnetic diffusivity in simulations of star formation. We discuss the intimate relation between the processes of reconnection diffusion, field wandering and turbulent mixing of a magnetized media and show that the role of the plasma effects is limited to “breaking up lines” on small scales and does not affect the rate of reconnection diffusion. We address the existing observational results and demonstrate how reconnection diffusion can explain the puzzles presented by observations, in particular, the observed higher magnetization of cloud cores in comparison with the magnetization of envelopes. We also outline a possible set of observational tests of the reconnection diffusion concept and discuss how the application of the new concept changes our understanding of star formation and its numerical modeling. Finally, we outline the differences of the process of reconnection diffusion and the process of accumulation of matter along magnetic field lines that is frequently invoked to explain the results of numerical simulations.     

Theory and Simulation of Reconnection

Reconnection is a major commonality of solar and magnetospheric physics. It was conjectured by Giovanelli in 1946 to explain particle acceleration in solar flares near magnetic neutral points. Since than it has been broadly applied in space physics including magnetospheric physics. In a special way this is due to Harry Petschek, who in 1994 published his ground breaking solution for a 2D magnetized plasma flow in regions containing singularities of vanishing magnetic field. Petschek’s reconnection theory was questioned in endless disputes and arguments, but his work stimulated the further investigation of this phenomenon like no other. However, there are questions left open. We consider two of them – “anomalous” resistivity in collisionless space plasma and the nature of reconnection in three dimensions. The CLUSTER and SOHO missions address these two aspects of reconnection in a complementary way -- the resistivity problem in situ in the magnetosphere and the 3D aspect by remote sensing of the Sun. We demonstrate that the search for answers to both questions leads beyond the applicability of analytical theories and that appropriate numerical approaches are necessary to investigate the essentially nonlinear and nonlocal processes involved. Necessary are both micro-physical, kinetic Vlasov-equation based methods of investigation as well as large scale (MHD) simulations to obtain the geometry and topology of the acting fields and flows.     

Turbulence, Magnetic Reconnection in Turbulent Fluids and Energetic Particle Acceleration

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én, 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 and Vishniac (Astrophys. J. 517:700–718, 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.     

Magnetosphere–Ionosphere Convection as a Compound System

Convection is the most fundamental process in understanding the structure of geospace and disturbances observed in the magnetosphere–ionosphere (M–I) system. In this paper, a self-consistent configuration of the global convection system is considered under the real topology as a compound system. Investigations are made based on the M–I coupling scheme by analyzing numerical results obtained from magnetohydrodynamic (MHD) simulations which guarantee the self-consistency in the whole system under the Bv (magnetic field and velocity) paradigm. It is emphasized in the M–I coupling scheme that convection and field-aligned current (FAC) are different aspects of same physical process characterizing the open magnetosphere. Special attention is given in this paper to the energy supplying (dynamo) process that drives the FAC system. In the convection system, the dynamo must be constructed from shear motion together with plasma population regimes to steadily drive the convection. Convection patterns observed in the ionosphere are also the manifestation of achievement in global self-consistency. A primary approach to apply these concepts to the study of geospace is to consider how the M–I system adjusts the relative motion between the compressible magnetosphere and the incompressible ionosphere when responding to given solar-wind conditions. The above principle is also applicable for the study of disturbance phenomena such as the substorm as well as for the study of apparently unique processes such as the flux transfer event (FTE), the sudden commencement (SC), and the theta aurora. Finally, an attempt is made to understand the substorm as the extension of enhanced convection under the southward interplanetary magnetic field (IMF) condition.     

MHD ANALYSIS OF PETSCHEK-TYPE RECONNECTION IN NON-UNIFORM FIELD AND FLOW GEOMETRIES

Analytical studies of reconnection have, for the most part, been confined to steady and uniform current sheet geometries. In contrast to these implifications, natural phenomena associated with the presence of current sheets indicate highly non-uniform structure and time-varying behaviour. Examples include the violent outbursts of energy on the Sun known as solar flares, and magnetospheric phenomena such as flux transfer events, plasmoids, and auroral activity. Unlike the theoretical models, reconnection therefore occurs in a highly dynamic and structured plasma environment. In this article we review the mathematical tools and techniques which are available to formulate models capable of describing the effects of reconnection in such situations. We confine attention to variants of the reconnection model first discussed by Petschek in the 1960s, in view of its successful application in predicting and interpreting phenomena in the terrestrial magnetosphere. The analysis of Petschek-type reconnection is based on the equations of ideal magnetohydrodynamics (MHD), which describe the large-scale behaviour of the magnetic field and plasma flow outside the diffusion region, which we determine as a localised part of the current sheet in which reconnection is initiated. The approach we adopt here is to transform the MHD equations into a Lagrangian or so-called 'frozen-in' coordinate system. In this coordinate system, the equation of motion transforms into a set of coupled nonlinear equations, in which the presence of inhomogeneous magnetic fields and/or plasma flows gives rise to a term similar to that which appears in the study of the ordinary string equation in a non-homogeneous medium. As demonstrated here, this approach not only clarifies and highlights the effects of such non-uniformities, it also simplifies the solution of the original set of MHD equations. In particular, this is true for those types of problem in which the total pressure can be considered as a known quantity from the outset. To illustrate the method, we solve several 2D problems involving magnetic field and flow non-uniformities: reconnection in a stagnation-point flow geometry with antiparallel magnetic fields; reconnection in a Y-type magnetic field geometry with and without velocity shear across the current sheet; and reconnection in a force-free magnetic field geometry with field lines of the form xy = const. These case examples, chosen for their tractability, each incorporate some aspects of the field and flow geomtries encountered in solar-terrestrial applications, and they provide a starting point for further analytical as well as numerical studies of reconnection.     

Electrodynamic coupling of the magnetosphere and ionosphere

The auroral zone ionosphere is coupled to the outer magnetosphere by means of field-aligned currents. Parallel electric fields associated with these currents are now widely accepted to be responsible for the acceleration of auroral particles. This paper will review the theoretical concepts and models describing this coupling. The dynamics of auroral zone particles will be described, beginning with the adiabatic motions of particles in the converging geomagnetic field in the presence of parallel potential drops and then considering the modifications to these adiabatic trajectories due to wave-particle interactions. The formation of parallel electric fields can be viewed both from microscopic and macroscopic viewpoints. The presence of a current carrying plasma can give rise to plasma instabilities which in a weakly turbulent situation can affect the particle motions, giving rise to an effective resistivity in the plasma. Recent satellite observations, however, indicate that the parallel electric field is organized into discrete potential jumps, known as double layers. From a macroscopic viewpoint, the response of the particles to a parallel potential drop leads to an approximately linear relationship between the current density and the potential drop.The currents flowing in the auroral circuit must close in the ionosphere. To a first approximation, the ionospheric conductivity can be considered to be constant, and in this case combining the ionospheric Ohm's Law with the linear current-voltage relation for parallel currents leads to an outer scale length, above which electric fields can map down to the ionosphere and below which parallel electric fields become important. The effects of particle precipitation make the picture more complex, leading to enhanced ionization in upward current regions and to the possibility of feedback interactions with the magnetosphere.Determining adiabatic particle orbits in steady-state electric and magnetic fields can be used to determine the self-consistent particle and field distributions on auroral field lines. However, it is difficult to pursue this approach when the fields are varying with time. Magnetohydrodynamic (MHD) models deal with these time-dependent situations by treating the particles as a fluid. This class of model, however, cannot treat kinetic effects in detail. Such effects can in some cases be modeled by effective transport coefficients inserted into the MHD equations. Intrinsically time-dependent processes such as the development of magnetic micropulsations and the response of the magnetosphere to ionospheric fluctuations can be readily treated in this framework.The response of the lower altitude auroral zone depends in part on how the system is driven. Currents are generated in the outer parts of the magnetosphere as a result of the plasma convection. The dynamics of this region is in turn affected by the coupling to the ionosphere. Since dissipation rates are very low in the outer magnetosphere, the convection may become turbulent, implying that nonlinear effects such as spectral transfer of energy to different scales become important. MHD turbulence theory, modified by the ionospheric coupling, can describe the dynamics of the boundary-layer region. Turbulent MHD fluids can give rise to the generation of field-aligned currents through the so-called -effect, which is utilized in the theory of the generation of the Earth's magnetic field. It is suggested that similar processes acting in the boundary-layer plasma may be ultimately responsible for the generation of auroral currents.     

Particle Acceleration in the Magnetotail and Aurora

This paper deals with acceleration processes in the magnetotail and the processes that enhance particle precipitation from the tail into the ionosphere through electric fields in the auroral acceleration region, generating or intensifying discrete auroral arcs. Particle acceleration in the magnetotail is closely related to substorms and the occurrence, and consequences, of magnetic reconnection. We discuss major advances in the understanding of relevant acceleration processes on the basis of simple analytical models, magnetohydrodynamic and test particle simulations, as well as full electromagnetic particle-in-cell simulations. The auroral acceleration mechanisms are not fully understood, although several, sometimes competing, theories and models received experimental support during the last decades. We review recent advances that emphasize the role of parallel electric fields produced by quasi-stationary or Alfvénic processes.     

Magnetic Reconnection in Extreme Astrophysical Environments

Magnetic reconnection is a fundamental plasma physics process in which ideal-MHD??s frozen-in constraints are broken and the magnetic field topology is dramatically re-arranged, which often leads to a violent release of the free magnetic energy. Most of the magnetic reconnection research done to date has been motivated by the applications to systems such as the solar corona, Earth??s magnetosphere, and magnetic confinement devices for thermonuclear fusion. These environments have relatively low energy densities and the plasma is adequately described as a mixture of equal numbers of electrons and ions and where the dissipated magnetic energy always stays with the plasma. In contrast, in this paper I would like to introduce a different, new direction of research??reconnection in high energy density radiative plasmas, in which photons play as important a role as electrons and ions; in particular, in which radiation pressure and radiative cooling become dominant factors in the pressure and energy balance. This research is motivated in part by rapid theoretical and experimental advances in High Energy Density Physics, and in part by several important problems in modern high-energy astrophysics. I first discuss some astrophysical examples of high-energy-density reconnection and then identify the key physical processes that distinguish them from traditional reconnection. Among the most important of these processes are: special-relativistic effects; radiative effects (radiative cooling, radiation pressure, and radiative resistivity); and, at the most extreme end??QED effects, including pair creation. The most notable among the astrophysical applications are situations involving magnetar-strength fields (10¹⁴?C10¹⁵ G, exceeding the quantum critical field B _??4×10¹³ G). The most important examples are giant flares in soft gamma repeaters (SGRs) and magnetic models of the central engines and relativistic jets of Gamma Ray Bursts (GRBs). The magnetic energy density in these environments is so high that, when it is suddenly released, the plasma is heated to ultra-relativistic temperatures. As a result, electron-positron pairs are created in copious quantities, dressing the reconnection layer in an optically thick pair coat, thereby trapping the photons. The plasma pressure inside the layer is then dominated by the combined radiation and pair pressure. At the same time, the timescale for radiation diffusion across the layer may, under some conditions, still be shorter than the global (along the layer) Alfvén transit time, and hence radiative cooling starts to dominate the thermodynamics of the problem. The reconnection problem then becomes essentially a radiative transfer problem. In addition, the high pair density makes the reconnection layer highly collisional, independent of the upstream plasma density, and hence radiative resistive MHD applies. The presence of all these processes calls for a substantial revision of our traditional physical picture of reconnection when applied to these environments and thus opens a new frontier in reconnection research.     

Magnetic Reconnection for Coronal Conditions: Reconnection Rates, Secondary Islands and Onset

Magnetic reconnection may play an important role in heating the corona through a release of magnetic energy. An understanding of how reconnection proceeds can contribute to explaining the observed behavior. Here, recent theoretical work on magnetic reconnection for coronal conditions is reviewed. Topics include the rate that collisionless (Hall) reconnection proceeds, the conditions under which Hall reconnection begins, and the effect of secondary islands (plasmoids) both on the scaling and properties of collisional (Sweet-Parker) reconnection and on the onset of Hall reconnection. Applications to magnetic energy storage and release in the corona are discussed.     

Current Status of Turbulent Dynamo Theory

Several recent advances in turbulent dynamo theory are reviewed. High resolution simulations of small-scale and large-scale dynamo action in periodic domains are compared with each other and contrasted with similar results at low magnetic Prandtl numbers. It is argued that all the different cases show similarities at intermediate length scales. On the other hand, in the presence of helicity of the turbulence, power develops on large scales, which is not present in non-helical small-scale turbulent dynamos. At small length scales, differences occur in connection with the dissipation cutoff scales associated with the respective value of the magnetic Prandtl number. These differences are found to be independent of whether or not there is large-scale dynamo action. However, large-scale dynamos in homogeneous systems are shown to suffer from resistive slow-down even at intermediate length scales. The results from simulations are connected to mean field theory and its applications. Recent work on magnetic helicity fluxes to alleviate large-scale dynamo quenching, shear dynamos, nonlocal effects and magnetic structures from strong density stratification are highlighted. Several insights which arise from analytic considerations of small-scale dynamos are discussed.     

Dayside Proton Aurora: Comparisons between Global MHD Simulations and IMAGE Observations

The IMAGE mission provides a unique opportunity to evaluate the accuracy of current global models of the solar wind interaction with the Earth's magnetosphere. In particular, images of proton auroras from the Far Ultraviolet Instrument (FUV) onboard the IMAGE spacecraft are well suited to support investigations of the response of the Earth's magnetosphere to interplanetary disturbances. Accordingly, we have modeled two events that occurred on June 8 and July 28, 2000, using plasma and magnetic field parameters measured upstream of the bow shock as input to three-dimensional magnetohydrodynamic (MHD) simulations. This paper begins with a discussion of images of proton auroras from the FUV SI-12 instrument in comparison with the simulation results. The comparison showed a very good agreement between intensifications in the auroral emissions measured by FUV SI-12 and the enhancement of plasma flows into the dayside ionosphere predicted by the global simulations. Subsequently, the IMAGE observations are analyzed in the context of the dayside magnetosphere's topological changes in magnetic field and plasma flows inferred from the simulation results. Finding include that the global dynamics of the auroral proton precipitation patterns observed by IMAGE are consistent with magnetic field reconnection occurring as a continuous process while the IMF changes in direction and the solar wind dynamic pressure varies. The global simulations also indicate that some of the transient patterns observed by IMAGE are consistent with sporadic reconnection processes. Global merging patterns found in the simulations agree with the antiparallel merging model, though locally component merging might broaden the merging region, especially in the region where shocked solar wind discontinuities first reach the magnetopause. Finally, the simulations predict the accretion of plasma near the bow shock in the regions threaded by newly open field lines on which plasma flows into the dayside ionosphere are enhanced. Overall the results of these initial comparisons between global MHD simulation results and IMAGE observations emphasize the interplay between reconnection and dynamic pressure processes at the dayside magnetopause, as well as the intricate connection between the bow shock and the auroral region.     

Low frequency turbulence and energy dissipation in the Solar Wind

The possibility to perform in-situ measurements of velocity, magnetic field, density and temperature fluctuations in the Solar Wind has greatly improved our knoweledge of MHD turbulence not only from the point of view of space physics but also from the more general point of view of plasma physics.These fluctuations on the one hand extend over a wide range of frequencies (about 5 decades), a fact which seems to be the signature of turbulent non-linear energy cascade, on the other hand display, mainly in the trailing edge of high speed streams, a number of striking features: (i) a high degree of correlation between magnetic and velocity field fluctuations, (ii) a very low level of fluctuations in mass density and magnetic field intensity, (iii) a considerable anisotropy revealed by minimum variance analysis of the magnetic field correlation tensor. More recently it has been stressed that MHD turbulence in the Solar Wind displays a clear intermittent character.The picture which emerges from the most recent analytical theories and numerical simulations is presented. In particular the observations which give us informations about the dissipation mechanism, which remains yet largely unknown, are discussed.     

Topics in Microphysics of Relativistic Plasmas

Astrophysical plasmas can have parameters vastly different from the more studied laboratory and space plasmas. In particular, the magnetic fields can be the dominant component of the plasma, with energy-density exceeding the particle rest-mass energy density. Magnetic fields then determine the plasma dynamical evolution, energy dissipation and acceleration of non-thermal particles. Recent data coming from astrophysical high energy missions, like magnetar bursts and Crab nebula flares, point to the importance of magnetic reconnection in these objects. In this review we outline a broad spectrum of problems related to the astrophysical relevant processes in magnetically dominated relativistic plasmas. We discuss the problems of large scale dynamics of relativistic plasmas, relativistic reconnection and particle acceleration at reconnecting layers, turbulent cascade in force-fee plasmas. A number of astrophysical applications are also discussed.     

Grid-Adaptive Computations of Magnetized Jets

We present grid-adaptive numerical simulations of magnetized plasma jets, modeled by means of the compressible magnetohydrodynamic equations. The Adaptive Mesh Refinement strategy makes it possible to investigate long-term jet dynamics where both large-scale and small-scale effects are at play. We extend recent findings for uniformly magnetized, periodic shear layers to planar and fully 3D extended jet segments. The jet lengths cover multiple, typically 10, axial wavelengths of the fastest growing Kelvin–Helmholtz (KH) like modes. The dominant linear MHD instabilities of the jet flows are quantified by means of MHD spectroscopic analysis. In cases characterized by sonic Mach numbers about unity and large plasma beta values, both single and double shear layers (planar jets) manifest self-organizing trends to large scales, e.g. by continuous pairing/merging between co-rotating vortices, simultaneously with the introduction of small-scale features by magnetic reconnection events. The vortices form as a result of KH unstable shear-flow layers, and their coalescence arises from the growth of subharmonic modes at multiple wavelengths of the fastest growing KH instability. In extended two-dimensional jet segments, we investigate how varying jet width alters this coalescence process occurring at both edges, e.g. by introducing Batchelor-like coupling between counter-rotating vortices formed at opposing weakly magnetized, close shear layers. Finally, periodic segments of supersonic magnetized jets are simulated in two- and three-dimensional cases, which are characterized by violent shock-dominated transients.